亚洲免费av电影一区二区三区,日韩爱爱视频,51精品视频一区二区三区,91视频爱爱,日韩欧美在线播放视频,中文字幕少妇AV,亚洲电影中文字幕,久久久久亚洲av成人网址,久久综合视频网站,国产在线不卡免费播放

        ?

        考慮重力及氣體非達西效應(yīng)影響的邊水氣藏邊水突破時間預(yù)測模型

        2016-07-21 07:52:08汪周華王子敦鄧丹郭平劉寧西南石油大學油氣藏地質(zhì)及開發(fā)國家重點實驗室中國石油華北油田公司
        石油鉆采工藝 2016年2期
        關(guān)鍵詞:邊水達西水氣

        汪周華王子敦鄧丹郭平劉寧.西南石油大學油氣藏地質(zhì)及開發(fā)國家重點實驗室;.中國石油華北油田公司

        考慮重力及氣體非達西效應(yīng)影響的邊水氣藏邊水突破時間預(yù)測模型

        汪周華1王子敦1鄧丹1郭平1劉寧2
        1.西南石油大學油氣藏地質(zhì)及開發(fā)國家重點實驗室;2.中國石油華北油田公司

        現(xiàn)有邊水突破時間預(yù)測模型都是基于平面,然而實際氣藏與水平面總存在一定的傾角,所以不能忽略重力作用的影響。基于氣水兩相滲流力學理論,綜合考慮地層傾角、氣體非達西流動效應(yīng)、氣水流度比和氣井距邊水的長度等因素,建立實際傾斜邊水氣藏新模型,并進行敏感性分析。實例分析表明,與已有預(yù)測模型相比,新模型預(yù)測的邊水突破時間更接近實際。敏感性分析表明,邊水突破時間與地層傾角呈二次拋物線關(guān)系,在地層傾角為90°時邊水突破時間達到最大;與氣井距邊水的長度呈冪函數(shù)關(guān)系,在氣井距邊水的長度達到500 m后,邊水突破時間的增加量每100 m增加7 d;與氣井產(chǎn)量及氣體非達西系數(shù)都呈反比關(guān)系,在開發(fā)前期邊水突破時間的降低程度接近90%,當氣井水淹后則基本對邊水突破時間無影響;與氣水流度比及儲層厚度都呈線性關(guān)系,邊水突破時間在氣水流度比值每增加1倍時增加30 d,在儲層厚度每增加1 m時則增大49 d。因此在氣井生產(chǎn)初期,準確確定這些參數(shù)顯得尤為重要。研究成果可對實際邊水氣藏的高效開發(fā)提供技術(shù)支撐。

        邊水氣藏;突破時間;預(yù)測模型;重力;非達西效應(yīng)

        對于邊底水氣藏,無論底水錐進還是邊水舌進,都會消耗舉升能量,減小氣井產(chǎn)量,因此,準確預(yù)測邊底水的突破時間是氣井高效開發(fā)的關(guān)鍵。許多學者運用不同的方法給出了不少預(yù)測邊底水突破時間的模型[1-3],其中對于底水突破時間的研究較多[4-6]。但對于邊水突破時間,一般采用數(shù)值模擬方法進行預(yù)測[7],史乃光等[8]運用系統(tǒng)辨識方法建立了預(yù)測公式;Sobocinski[9]、Kuo[10]、王會強[11]等運用滲流力學理論推導(dǎo)出了預(yù)測邊水突破時間的公式,楊芙蓉等[12]推導(dǎo)出適合高產(chǎn)氣井的邊水突破時間預(yù)測公式;修乃嶺[13]基于管流理論建立了一種適合于縫洞型儲層的新流動模型,得出裂縫寬度決定見水時間;Ahmadi[14]等人提出利用最小二乘支持向量機(LSSVM)預(yù)測水錐突破時間,計算精度高,可用于數(shù)值模擬軟件開發(fā)。

        上述邊水突破時間的預(yù)測模型都是基于平面的,但是實際氣藏必定與水平面存在一定的傾角,所以重力作用不能忽略。為此,以氣水兩相滲流力學理論為基礎(chǔ),建立了實際傾斜邊水氣藏新模型,并考慮地層傾角、氣體非達西流動效應(yīng)、氣水流度比和氣井距邊水的長度等因素,推導(dǎo)出傾斜邊水氣藏氣井邊水突破時間的計算公式,并進行了實例分析,研究結(jié)果對于實際邊水氣藏的高效開發(fā)具有指導(dǎo)作用。

        1 模型建立與公式推導(dǎo)Model building and formula derivation

        1.1物理模型

        Physical model

        有一邊水氣藏如圖1所示,儲層傾角為α。原始氣水邊界近似是一水平面(A0B0),與氣井相距L;若氣井開始生產(chǎn),則氣水分界面向井底運動,此時氣水界面變成曲面(A1B1),假設(shè)初始氣水界面處有一個水質(zhì)點A,在時間t后滲流到C點。為計算方便又可滿足物理變化過程,作以下假設(shè):(1)儲層等厚均質(zhì),水活塞式驅(qū)替氣體流向井底;(2)滲流過程中氣相、水相的密度和黏度都是定值;(3)水相服從達西滲流,氣相服從非達西滲流;(4)考慮重力而忽略毛管力作用。

        圖1 傾斜邊水氣藏生產(chǎn)示意圖Fig.1 Production in inclined gas reservoirs with edge water

        1.2公式推導(dǎo)

        Formula derivation

        根據(jù)氣水兩相滲流力學理論[15],氣水兩相滲流方程為

        不考慮毛管力作用,則在氣水界面處有

        根據(jù)水質(zhì)點在多孔介質(zhì)中運動時間與運動距離的關(guān)系,再考慮多孔介質(zhì)中存在的束縛水和殘余氣,并且初始時刻(t=0),x=0;氣井見水時(t=tbt),x=L,則氣井見水時間為

        其中

        式(6)即為傾斜均質(zhì)地層邊水氣藏氣井見水時間的預(yù)測公式。

        如果不考慮重力和氣體非達西效應(yīng)的影響,則式(6)變?yōu)槭剑?)[11]

        若不考慮重力作用,但考慮氣體非達西效應(yīng)的影響,則式(6)變?yōu)槭剑?)[12]

        2 實例分析Case study

        2.1氣井基礎(chǔ)數(shù)據(jù)

        Basic data of gas wells

        根據(jù)某海上油田2口氣井的基礎(chǔ)數(shù)據(jù)對本文的預(yù)測模型進行計算分析?;A(chǔ)數(shù)據(jù)見表1和表2。

        2.2邊水突破時間計算

        Determination of edge-water breakthrough time

        根據(jù)表1和表2的數(shù)據(jù),結(jié)合式(6)、(7)、(8),可以計算得到2口井的邊水突破時間。表3中模型1是在氣/水達西滲流條件下推導(dǎo)出的邊水突破時間預(yù)測公式,模型2是考慮氣體非達西效應(yīng)影響時推導(dǎo)出的邊水突破時間預(yù)測公式。

        表 1 氣井Ⅰ基礎(chǔ)數(shù)據(jù)Table 1 Basic data of Gas Well-Ⅰ

        表 2 氣井Ⅱ基礎(chǔ)數(shù)據(jù)Table 2 Basic data of Gas Well-Ⅱ

        表 3 邊水突破預(yù)測時間Table 3 Predicted edge-water breakthrough time

        由表3可以看出,由于氣體非達西效應(yīng)的影響,使氣體流速增大,所以邊水突破時間增快,因此模型2較模型1計算的結(jié)果低;但是重力作用使氣液流速都減小,所以邊水突破時間減緩,因此新模型較模型2計算的結(jié)果高。

        由圖2看出,在實際生產(chǎn)過程中,541 d左右Ⅰ井水氣比和氣井日產(chǎn)水量都開始增加,隨后氣井日產(chǎn)水量急劇增加,在800 d左右不得不采取措施(關(guān)閉出水層位的滑套)降低出水量;Ⅱ井在644 d左右水氣比和氣井日產(chǎn)水量都開始增加,但是沒有及時采取必要措施,在開發(fā)后期氣井產(chǎn)水量急劇上升,使氣井水淹,無法正常生產(chǎn)。所以,Ⅰ井實際邊水突破時間為541 d,Ⅱ井實際邊水突破時間為644 d,新模型計算結(jié)果與實際結(jié)果更為相近。

        圖 2 氣井實際產(chǎn)水曲線Fig.2 Actual in water productivity curve of the gas well

        3 敏感性分析Sensitivity analysis

        影響邊水突破時間的因素很多,如地層傾角、氣井距邊水長度、氣井產(chǎn)量、氣水流度比、氣體非達西系數(shù)和氣層厚度等等。以表1和表2中數(shù)據(jù)為基礎(chǔ),假定某一因素的變化值,利用新模型就關(guān)鍵的6項影響因素對邊水突破時間的影響程度進行討論。

        3.1地層傾角的影響

        Effects of dip angles

        由式(6)可知,C與α呈正比關(guān)系,而C與tbt呈反比,因此tbt與α呈反比,而α變大,則重力作用就變大,重力是氣液在多孔介質(zhì)中流動的阻力,所以圖3中顯示邊水突破時間隨地層傾角的增大而增大,在α<30°時,地層傾角每增加10°邊水突破時間增加量約為60 d;在α>30°時,隨地層傾角增加邊水突破時間的增加量逐漸減小,從80°到90°邊水突破時間僅增加3 d;在α=90°時邊水突破時間達到最大,此時邊水氣藏也變?yōu)榈姿畾獠亍?/p>

        3.2氣井距邊水長度的影響

        Effects of the distance between gas well and edge water

        圖3 邊水突破時間與地層傾角關(guān)系Fig.3Correlation between edge-water breakthrough time and dip angle

        氣井距離邊水的長度增長,則邊水到井底的流經(jīng)路程變長,在一定流度下邊水突破時間必然增大。由式(6)也可得知,tbt與L呈現(xiàn)正比關(guān)系,所以如圖4所示,邊水突破時間與氣井距離邊水的長度呈現(xiàn)冪函數(shù)關(guān)系,隨氣井距離邊水的長度增長而邊水突破時間增大,在500 m之前邊水突破時間緩慢增加,在500 m后邊水突破時間迅速增加,氣井距離邊水的長度每增加100 m時邊水突破時間的增加量增加7 d。

        圖4 邊水突破時間與氣井距離邊水長度關(guān)系Fig.4Correlation between edge-water breakthrough time and distance between gas well and edge water

        3.3氣井產(chǎn)量的影響

        Effects of gas well productivity

        由式(6)可知,氣井產(chǎn)量影響著A、B兩個參數(shù),并且與之呈反比關(guān)系。由實際生產(chǎn)過程也可知,在其他條件不變的情況下,氣井產(chǎn)量增大則氣體流速必然變快,地層中氣體減少速度變快,所以邊水的舌進速度就會增大,從而邊水突破時間變小。因此,如圖5所示,邊水突破時間隨氣井產(chǎn)量的增大而減小,氣井產(chǎn)量小于4×104m3/d時,隨氣井產(chǎn)量的增大而邊水突破時間迅速減小,降低程度達到88.64%;氣井產(chǎn)量大于4×104m3/d之后,邊水突破時間減小趨勢變緩;當氣井產(chǎn)量大于15×104m3/d時,由于邊水入侵速度很快,邊水突破時間基本不變,所以氣井產(chǎn)量繼續(xù)增大對邊水突破時間影響已經(jīng)不大。

        3.4氣水流度比的影響

        Effects of gas-water mobility ratio

        由式(6)可知,tbt與Mgw呈正比關(guān)系。氣水流度比增大,則氣水流度屬性差異增大,氣、水兩相在多孔介質(zhì)中的流動也就越困難?;谝陨蟽牲c,如圖6所示,邊水突破時間與氣水流度比呈現(xiàn)線性關(guān)系,隨氣水流度比的增加而邊水突破時間增大,氣水流度比值每增加1倍時邊水突破時間增加30 d。

        圖 5 邊水突破時間與氣井產(chǎn)量關(guān)系Fig.5Correlation between edge-water breakthrough time and productivity of gas well

        圖 6 邊水突破時間與氣水流度比關(guān)系Fig.6Correlation between edge-water breakthrough time and gas-water mobility ratio

        3.5氣體非達西系數(shù)的影響

        Effects of gas non-Darcy factor

        由式(6)可知,氣體非達西系數(shù)β主要影響參數(shù)B,二者呈現(xiàn)正比關(guān)系;并且β增大,則氣體流速必然增加,所以邊水舌進速度增加。因此,如圖7所示。

        圖 7 邊水突破時間與非達西系數(shù)關(guān)系Fig.7Correlation between edge-water breakthrough time and non-Darcy factor

        邊水突破時間隨氣體非達西系數(shù)的增大而減小,氣體非達西系數(shù)由0增大到1.31×1011m-1過程中,邊水突破時間從6 204 d驟降到617 d,降低程度達到90.05%;當氣體非達西系數(shù)大于5×1011m-1,邊水舌進速度很快,氣井開發(fā)即見水,致使邊水突破時間基本保持不變,這時氣體非達西系數(shù)繼續(xù)增大對邊水突破時間影響已經(jīng)不大了。

        3.6儲層厚度的影響

        Effects of reservoir thickness

        氣體流速vg與儲層厚度h為負相關(guān),且由式(6)可知,儲層厚度h影響參數(shù)A和B,儲層越厚,氣藏儲量就越大,在一定的產(chǎn)量下氣體流速就會越低,所以邊水舌進速度降低。因此,如圖8所示,邊水突破時間與儲層厚度呈現(xiàn)線性關(guān)系,邊水突破時間隨儲層厚度的增加而增大,儲層厚度每增加1 m則邊水突破時間增大49 d。

        圖 8 邊水突破時間與儲層厚度關(guān)系Fig.8Correlation between edge-water breakthrough time and reservoir thickness

        4 結(jié)論Conclusions

        (1)通過建立實際傾斜邊水氣藏模型,推導(dǎo)出綜合考慮地層傾角、氣體非達西流動效應(yīng)、氣水流度比和氣井距邊水的長度等影響因素的邊水突破時間預(yù)測新公式,較已有模型更接近邊水氣藏實際情況。

        (2)實例分析結(jié)果表明,Ⅰ井邊水突破時間的新模型計算值與實際值的相對誤差為3.14%,Ⅱ井邊水突破時間的新模型計算值與實際值的相對誤差僅為1.71%,說明新公式預(yù)測的邊水突破時間與氣井生產(chǎn)實際更接近,用新公式預(yù)測邊水突破時間更準確,對邊水氣藏的高效開發(fā)具有一定實際指導(dǎo)意義。

        (3)敏感性分析表明,影響邊水突破時間的因素主要有地層傾角、氣井距邊水長度、氣井產(chǎn)量、氣水流度比、氣體非達西系數(shù)和氣層厚度等。這些因素對邊水突破時間影響較大,因此在氣井生產(chǎn)初期,準確確定這些參數(shù)顯得尤為重要。

        符號說明:

        Nomenclature:

        Bg為氣體體積系數(shù);g為重力加速度,m/s2;h為儲層有效厚度,m;kg、kw分別為氣相滲透率和水相滲透率,mD;kgwi、kwgr分別為束縛水飽和度下氣相相對滲透率和殘余氣飽和度下水相相對滲透率;L為初始氣水邊界距離氣井距離,m;Mgw為氣水流度比;pg、pw分別為x處氣相、水相壓力,MPa;qsc為氣井產(chǎn)量,104m3/d;T為地層溫度,K;Swi和Sgr分別為束縛水飽和度和殘余氣飽和度,小數(shù);tbt為邊水突破時間,d;μg、μw分別為氣體黏度和水黏度,mPa·s;vg為氣相滲流速度,m/s;vw為氣相滲流速度,m/s;x為水質(zhì)點距離氣井距離,m;Z為偏差系數(shù),無量綱;α為儲層傾角,°;β為氣體非達西系數(shù),m-1;φ為儲層有效孔隙度,小數(shù);ρg、ρw分別為氣體密度和水密度,g/cm3。

        References:

        [1]LEE S H, TUNG W B. General coningCorrelations based on mechanistic studies[R]. SPE 20742, 1990.

        [2]趙新智,朱圣舉.低滲透帶隔板底水油藏油井見水時間預(yù)測[J].石油勘探與開發(fā),2012,39(4): 471-474.

        ZHAO Xinzhi, ZHU Shengju. Prediction of water breakthrough time for oil wells in low-permeability bottom water reservoirs with barrier [J]. Petroleum Exploration and Development, 2012, 39(4): 504-507.

        [3]ROOPA I V, BAKSH K. Predicting water breakthrough using pressure transient analysis in a gas condensate reservoir in the columbus basin, Trinidad[R]. SPE 158496, 2012.

        [4]黃炳光,劉蜀知.實用油藏工程與動態(tài)分析方法[M].北京:石油工業(yè)出版社,1997:199-202.

        HUANG Bingguang, LIU Shuzhi. Pragmatic methods of Petroleum reservoir engineering and dynamic analysis [M].Beijing: Petroleum Industry Press,1997: 199-202.

        [5]王會強,李曉平,楊琪,龔偉,曾玉強.底水氣藏見水時間預(yù)測方法[J].新疆石油地質(zhì),2007,28(1):92-93.

        WANG Huiqiang, LI Xiaoping, YANG Qi, GONG Wei,ZENG Yu-qiang. A prediction of water breakthrough time in gas reservoirs with bottom water [J].Xinjiang Petroleum Geology, 2007, 28(1): 92-93.

        [6]張烈輝,嚴謹,李允,李曉平,羅濤,胡勇,鐘兵,陳軍.考慮凝析油析出時的底水凝析氣藏見水時間預(yù)測新方法[J].天然氣工業(yè),2004,24(7):74-75.

        ZHANG Liehui, YAN Jin, LI Yun, LI Xiaoping, LUO Tao, HU Yong, ZHONG Bing, CHEN Jun. New method to predict condensing effect on water breakthrough in condensate reservoirs with bottom water [J].Natural Gas Industry, 2004, 24(7): 74-75.

        [7]石廣志,馮國慶,張烈輝.某邊水油藏開發(fā)數(shù)值模擬研究[J].天然氣勘探與開發(fā),2006,29(2):21-24.

        SHI Guangzhi, FENG Guoqing, ZHANG Liehui. Numerical simulation in an edge-water oil reservoir [J]. Natural Gas Exploration and Development, 2006, 29(2):21-24.

        [8]史乃光,楊維寧,楊正文.用系統(tǒng)辨識方法預(yù)測氣井見水時間[J].天然氣工業(yè),1992,12(1):25-31.

        SHI Naiguang, YANG Weining, YANG Zhengwen. Using system identification method to predict water breakthrough time in gas well [J]. Natural Gas Industry,1992, 12(1): 25-31.

        [9]SOBOCINSKI D P, CORNELIUS A J. ACorrelation for predicting water coning time [J]. Journal of Petroleum Technology, 1965, 17(5): 594-600.

        [10]KUO M C T. A simplified method for water coning predictions[R]. SPE 12067, 1983.

        [11]王會強,李曉平,吳鋒,何司平.邊水氣藏氣井見水時間預(yù)測方法[J].特種油氣藏,2008,15(4):73-74.

        WANG Huiqiang, LI Xiaoping, WU Feng, HE Siping. Prediction of water breakthrough time in gas reservoirs with edge water[J].Special Oil & Gas Reservoirs ,2008,15(4): 73-74.

        [12]楊芙蓉,樊平天,賀靜,吳克柳,李異才,李娜.邊水氣藏高產(chǎn)氣井見水時間預(yù)測方法[J].科學技術(shù)與工程,2013,13 (29):8745-8747,8754.

        YANG Furong, FAN Pingtian, HE jing,WU Keliu,LI Yicai, LI Na. An approach to predict water breakthrough time in high gas rate wells of gas reservoirs with edge water[J].Science Technology and Engineering,2013,13(29): 8745-8747, 8754.

        [13]修乃嶺,耿忠娟,熊偉,高樹生,薛惠.縫洞型碳酸鹽巖油藏開發(fā)特征和水動力學模擬[J].石油鉆采工藝,2008,30(2):72-74.

        XIU Nailing,GENG Zhongjuan, XIONG wei, GAO Shusheng, XUE Hui. Features and hydrogynamic simulation of fractured-vuggy reservoir development[J]. Oil Drilling & Production Technology, 2008, 30(2): 72-74.

        [14]AHMADI M A, EBADI M, HOSSEINI S M. Prediction breakthrough time of water coning in the fractured reservoirs by implementing low parameter support vector machine approach [J]. Fuel, 2014, 117: 579-589.

        [15]程林松.高等滲流力學[M].北京:石油工業(yè)出版社,2011: 54-60.

        CHEN Linsong. Higher fluid mechanics in porous medium [M]. Beijing: Petroleum Industry Press, 2011:54-60.

        (修改稿收到日期 2016-01-18)

        〔編輯 朱 偉〕

        Model for prediction of edge-water breakthrough time in reservoirs with edge water with consideration to effects of gravity and non-Darcy effect of gases

        WANG Zhouhua1, WANG Zidun1, DENG Dan1, GUO Ping1, LIU Ning2
        1. State Key Laboratory of Southwest Petroleum University, Chengdu, Sichuan 610500. China;2. PetroChina Huabei Oilfield Company, Renqiu, Hebei 062552, China

        All existing models for prediction of edge-water breakthrough time are based on horizontal plane. However, actual gas reservoirs always have certain dips from the horizontal plane. It is necessary to consider effects of the gravity. According to the theories for gas-water flow mechanics, and with consideration to dip angles, non-Darcy flow effects of gases, gas/water mobility ratios,distances between gas well and the edge water, and other factors, a new model for inclined gas reservoirs with edge water was built. Then, sensitivity analysis was performed for the new model. Compared with existing prediction models, the new model can predictbreakthrough time of edge water closer to actual time. Sensitivity analysis shows that edge-water breakthrough time is in quadratic parabola relationship with dip angles. The maximum edge water breakthrough time can be observed at the dip angle of 90°. The edgewater breakthrough time is in power-function relationship with the distance between the gas well and the edge water. With such distance of 500 m or more, the edge-water breakthrough time increases 7 d for every 100 m. Moreover, the edge-water breakthrough time is in inverse relationship with both gas well productivity and non-Darcy factors of relevant gases. In earlier stage of development, the breakthrough time of edge water reduces up to 90%. When the gas well is flooded, effects on breakthrough time of the edge water are ignorable. Furthermore, the edge-water breakthrough time is in liner relation with both gas/water flow rate and reservoir thickness. The edge-water breakthrough time increases 30 d when the gas/water flow rate is doubled, and increases 49 d for every additional meter in reservoir thickness. It is very important to determine these parameters accurately in early stage of gas well production. The research results can provide necessary technical supports for high-efficiency development of gas reservoirs with edge water.

        gas reservoirs with edge water; breakthrough time; prediction model; gravity; non-Darcy effect

        WANG Zhouhua, WANG Zidun, DENG Dan, GUO Ping, LIU Ning. Model for prediction of edge-water breakthrough time in reservoirs with edge water with consideration to effects of gravity and non-Darcy effect of gases[J]. Oil Drilling & Production Technology, 2016, 38(2): 210-215.

        TE349

        A

        1000 -7393( 2016 ) 02 -0210-06

        10.13639/j.odpt.2016.02.016

        國家科技重大專項“致密砂巖氣有效開發(fā)評價技術(shù)”(編號:2011ZX05013-002)。

        汪周華(1979-),現(xiàn)主要從事氣田開發(fā)理論與方法、注氣提高采收率、低滲致密氣藏滲流機理等方面研究工作,副研究員,博士。通訊地址:(610500)四川省成都市新都區(qū)西南石油大學油氣藏地質(zhì)及開發(fā)工程國家重點實驗室B407。。E-mail:wangzhouhua@126.com

        王子敦(1990-),在讀碩士研究生,主要研究方向為氣田開發(fā)理論與方法、低滲致密及頁巖氣藏滲流機理。通訊地址:(610500)四川省成都市新都區(qū)西南石油大學油氣藏地質(zhì)及開發(fā)工程國家重點實驗室B204。E-mail:1453708420@qq.com

        引用格式:汪周華,王子敦,鄧丹,郭平,劉寧.考慮重力及氣體非達西效應(yīng)影響的邊水氣藏邊水突破時間預(yù)測模型[J].石油鉆采工藝,2016,38(2):210-215.

        猜你喜歡
        邊水達西水氣
        遼中區(qū)患病草魚體內(nèi)嗜水氣單胞菌分離、鑒定與致病力測定
        海上邊水氣藏利用試井資料確定水侵狀況研究
        海洋石油(2021年3期)2021-11-05 07:42:26
        春光油田稀油小砂體控水技術(shù)應(yīng)用研究
        石油研究(2019年2期)2019-09-10 07:22:44
        某油田區(qū)塊油藏水平井開發(fā)問題的研究
        石油研究(2019年1期)2019-09-10 07:22:44
        傲慢與偏見
        GC-MS法分析藏藥坐珠達西中的化學成分
        中成藥(2016年4期)2016-05-17 06:07:46
        弱邊水扇形小斷塊油藏井網(wǎng)參數(shù)優(yōu)化
        堤壩Forchheimei型非達西滲流場特性分析
        氮氣抑水增油技術(shù)研究及現(xiàn)場應(yīng)用
        醫(yī)院感染嗜水氣單胞菌的臨床治療分析
        特黄做受又粗又长又大又硬| 一本色道久久88加勒比—综合| 中文人妻av久久人妻水蜜桃| 国产精品成人va在线观看| 亚洲男人的天堂在线aⅴ视频 | 华人在线视频精品在线| 一本色道久久hezyo无码| 色视频线观看在线网站| 久青草久青草视频在线观看| 国产精品高潮呻吟av久久无吗| 在线播放a欧美专区一区| 欧洲一卡2卡三卡4卡免费网站| 在线观看免费午夜大片| 99热在线精品播放| 色婷婷色99国产综合精品| 国内偷拍视频一区二区| 国产精品一区二区熟女不卡| 国产精品亚洲а∨无码播放| 精品一区二区三区免费播放 | 伊人色综合九久久天天蜜桃| 国产少妇露脸精品自拍网站| 日韩在线精品视频一区| 含紧一点h边做边走动免费视频 | 粉嫩少妇内射浓精videos| 久久精品国产亚洲av成人擦边| 国产精品一区av在线| 国产亚洲精品美女久久久m| 99精品国产一区二区| 亚洲综合一区无码精品| 国产午夜激情视频自拍| 久久本道久久综合一人| 久久亚洲av成人无码电影| 国产精品第一国产精品| 91亚洲国产成人aⅴ毛片大全| 日本一区二区在线播放观看| 五月开心六月开心婷婷网| 久久aaaa片一区二区| 久久无码专区国产精品s| 日韩中文字幕欧美亚洲第一区| 久久亚洲精彩无码天堂| 精品精品国产一区二区性色av|