• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    氣井?dāng)y液用渦流工具結(jié)構(gòu)參數(shù)優(yōu)化

    2016-08-16 09:47:19陳德春付剛韓昊姚亞宋天驕謝雙喜
    石油鉆采工藝 2016年3期
    關(guān)鍵詞:導(dǎo)程攜液氣液

    陳德春付剛韓昊姚亞宋天驕謝雙喜

    1.中國石油大學(xué)(華東)石油工程學(xué)院;2.中海油能源發(fā)展股份有限公司工程技術(shù)分公司

    氣井?dāng)y液用渦流工具結(jié)構(gòu)參數(shù)優(yōu)化

    陳德春1付剛1韓昊1姚亞1宋天驕1謝雙喜2

    1.中國石油大學(xué)(華東)石油工程學(xué)院;2.中海油能源發(fā)展股份有限公司工程技術(shù)分公司

    為了提高氣井渦流排液采氣工具的工作效率,進(jìn)行了渦流工具結(jié)構(gòu)優(yōu)化研究。運(yùn)用AutoCAD和Fluent流體模擬軟件,建立了渦流工具氣液兩相流場模型,運(yùn)用正交試驗(yàn)確定了模擬計(jì)算方案,通過方差和極差分析得到了不同氣液條件下導(dǎo)程、導(dǎo)程數(shù)、槽深和槽寬等4個(gè)結(jié)構(gòu)參數(shù)對(duì)工具出口峰值速度的影響程度并確定了工具最優(yōu)結(jié)構(gòu)參數(shù)。結(jié)果表明,較小的槽深和槽寬、較大的導(dǎo)程和導(dǎo)程數(shù)可以起到很好的增速效果。槽深取5 mm,槽寬取40~45 mm,導(dǎo)程數(shù)取3倍,導(dǎo)程取172~219 mm時(shí),氣液混合物經(jīng)過渦流工具后氣相峰值速度最大,攜液效果最好。工具結(jié)構(gòu)參數(shù)對(duì)工具出口峰值速度的影響程度隨氣液條件而變化,液氣比較大時(shí),選擇較大的導(dǎo)程數(shù)具有更明顯的增速效果;入口速度較大時(shí),較小的槽寬具有更明顯的增速效果,而較大的導(dǎo)程增速效果并不明顯;較小的槽深始終具有明顯的增速效果。研究結(jié)果可以為現(xiàn)場渦流工具的優(yōu)選和使用效果的提高提供依據(jù)。

    氣井;排采技術(shù);渦流工具;結(jié)構(gòu)參數(shù)優(yōu)化;液氣比;數(shù)值模擬;正交試驗(yàn);敏感性分析

    井下渦流工具是一種新型的井下氣液分離裝置,不需要外界提供能量,無需增加管理人力?,F(xiàn)場試驗(yàn)證明[1-4]:使用渦流工具后氣井?dāng)y液量增大,井筒壓力損失降低,可為海上油氣田提供一種經(jīng)濟(jì)合理解決井底積液問題的方案[5]。張翠婷[6](2012)對(duì)渦流工具的分氣效果進(jìn)行數(shù)值模擬,并對(duì)其結(jié)構(gòu)尺寸進(jìn)行優(yōu)化;馮翠菊[7]等(2013)對(duì)影響渦流工具攜液效果的主要因素進(jìn)行了研究;吳丹[8]等(2014)研究了工具不同的結(jié)構(gòu)參數(shù)對(duì)氣井?dāng)y液效果的影響,并給出了工具參數(shù)的優(yōu)化結(jié)果。研究結(jié)果表明,影響渦流工具排水采氣效果的因素主要有2方面,一是渦流工具結(jié)構(gòu)參數(shù),包括導(dǎo)程、導(dǎo)程數(shù)(導(dǎo)程不變的情況下螺旋葉片的長度)、槽深和槽寬;二是井筒氣液條件,包括液氣比和氣液混合物速度。之前研究中使用單因素分析方法研究工具各結(jié)構(gòu)參數(shù)對(duì)攜液效果的影響,進(jìn)而確定各結(jié)構(gòu)參數(shù)的優(yōu)化結(jié)果。但由于工具各結(jié)構(gòu)參數(shù)之間存在相互影響,單因素分析中最優(yōu)的結(jié)構(gòu)參數(shù)的直接組合并不一定具有最好的攜液效果。此外,這些研究中的結(jié)構(gòu)參數(shù)優(yōu)化結(jié)果是在某一特定的氣液條件下得到的,具有局限性,在實(shí)際生產(chǎn)過程中,井的氣液條件不同,得到的最優(yōu)結(jié)構(gòu)參數(shù)也不相同。筆者通過建立渦流氣液兩相流場模型,采用數(shù)值模擬方法和正交試驗(yàn)方案,研究了不同氣液條件下渦流工具最優(yōu)結(jié)構(gòu)參數(shù)組合,分析渦流工具結(jié)構(gòu)參數(shù)對(duì)攜液效果影響的程度,為現(xiàn)場渦流工具結(jié)構(gòu)參數(shù)優(yōu)化提供依據(jù)。

    1 渦流氣液兩相流場模型

    Flow field model for swirl gas/liquid

    井下渦流工具主要由打撈頭、繞流器、導(dǎo)流筒、坐封器和接箍擋環(huán)組成,其中打撈頭用于配合打撈工具,繞流器用于改變流體的流態(tài),導(dǎo)流筒用于流體流通,坐封器和接箍擋環(huán)用于工具的固定[9]。利用AutoCAD建立繞流器和打撈頭結(jié)構(gòu)的幾何模型,如圖1所示。

    將建立的工具模型進(jìn)行網(wǎng)格劃分(圖2),并導(dǎo)入Fluent軟件中,設(shè)置模擬參數(shù)進(jìn)行計(jì)算。其中多相流模型選用歐拉模型,湍流模型選擇雷諾應(yīng)力模型[10]。

    圖1 螺旋變速體的幾何模型Fig. 1 Geometric model of helical gear

    圖2 網(wǎng)格模型圖Fig. 2 Network grid

    2 模擬計(jì)算方案

    Program for simulation and calculation

    使用渦流工具后,中心氣流阻力減小,氣體速度增大[11],臨界流速一定的情況下,氣相速度越大,攜液能力越強(qiáng),因此用出口氣相峰值速度作為工具攜液效果評(píng)判的指標(biāo)。

    將導(dǎo)程、導(dǎo)程數(shù)、槽深、槽寬4個(gè)參數(shù)作為研究因素??紤]到數(shù)值模擬的計(jì)算量較大,根據(jù)實(shí)際需要,每種因素取4 個(gè)水平,各因素各水平的取值見表1。由于工具結(jié)構(gòu)參數(shù)之間存在相互影響,理論上需要列舉所有的參數(shù)組合,得到256(44)種方案,進(jìn)行256次計(jì)算;為了減少工作量,根據(jù)正交表的設(shè)計(jì)方法,選擇L16(45)的正交表,只需計(jì)算16種方案即可得到精確的結(jié)果。

    表1 渦流工具結(jié)構(gòu)參數(shù)的因素水平表Table 1 Factor levels of structural parameters of swirl tools

    由于工具結(jié)構(gòu)參數(shù)的優(yōu)化是基于特定的井筒氣液條件,如果有多口井需要進(jìn)行工具結(jié)構(gòu)參數(shù)設(shè)計(jì),則每口井都需要做16組數(shù)值模擬計(jì)算,計(jì)算量很大。因此研究不同氣液條件下渦流工具的最優(yōu)結(jié)構(gòu)參數(shù)組合,應(yīng)用時(shí)只需根據(jù)井的實(shí)際氣液條件選擇對(duì)應(yīng)的優(yōu)化結(jié)果,方便快捷。

    3 結(jié)果分析

    Results

    3.1氣液條件一定時(shí)工具結(jié)構(gòu)參數(shù)優(yōu)化

    Optimization of structural parameters under designated gas/liquid conditions

    3.1.1極差分析 氣液條件取值:氣液混合物速度4 m/s,液氣比1∶50。A~D為渦流工具結(jié)構(gòu)參數(shù)的因素,由于在試驗(yàn)中留有空列可作為衡量試驗(yàn)因素和系統(tǒng)的誤差,并可以做方差分析,因此將E做為空列。將得到的16組出口氣相峰值速度進(jìn)行計(jì)算分析(見表2)。

    表2 極差分析數(shù)據(jù)(入口速度4 m/s,液氣比1∶50)Table 2 Data for range analysis (with inlet velocity of 4 m/s and liquid/gas ratio of 1∶50)

    通過表2可以看出影響出口氣相峰值速度的因素重要程度由大到小為:槽深>導(dǎo)程>導(dǎo)程數(shù)>槽寬。最優(yōu)的結(jié)構(gòu)參數(shù)組合為:導(dǎo)程218 mm,導(dǎo)程數(shù)3倍,槽深5 mm和槽寬40 mm。

    3.1.2方差分析 表3為方差分析結(jié)果,F(xiàn)值可以反映各因素水平改變時(shí)引起的差異顯著程度,越大說明因素對(duì)結(jié)果的影響越顯著。對(duì)于給定水平α=0.10,由P{F>λ}=0.10查F(3,12)分布表,得到λ=2.61??梢钥闯鰧?dǎo)程、導(dǎo)程數(shù)、槽深、槽寬的F>λ,即各工具結(jié)構(gòu)參數(shù)對(duì)出口速度峰值都有顯著影響。

    表3 方差分析數(shù)據(jù)(入口速度4 m/s,液氣比1∶50)Table 3 Data for variance analysis (with inlet velocity of 4 m/s and liquid/gas ratio of 1∶50)

    3.1.3最優(yōu)結(jié)構(gòu)參數(shù)下結(jié)果分析 建立最優(yōu)結(jié)構(gòu)參數(shù)下的工具模型,經(jīng)過Fluent模擬計(jì)算,出口峰值速度為4.421 m/s,比表2中所有方案的計(jì)算結(jié)果都要高,說明采用正交數(shù)值試驗(yàn)進(jìn)行工具參數(shù)優(yōu)化設(shè)計(jì)得到的結(jié)果是可靠的。得到的軸向方向水體分布見圖3。圖中中心紅色區(qū)域代表氣柱,邊緣綠色或者藍(lán)色區(qū)域代表液膜??梢钥闯黾尤霚u流工具之后,水的體積含量分布發(fā)生變化:管壁處的水的體積分?jǐn)?shù)增大,形成液膜;管道中心處的水的體積分?jǐn)?shù)減小,形成氣芯。

    圖3 軸向方向水體積分?jǐn)?shù)分布圖Fig. 3 Distribution of volume fractions of water along the axis

    圖4為氣液混合物經(jīng)過渦流工具后流動(dòng)軌跡曲線,可以看出:流經(jīng)渦流工具之后,井壁與螺旋葉片形成封閉空間,將流體轉(zhuǎn)變?yōu)檠芈菪€運(yùn)動(dòng),管中心附近會(huì)形成一個(gè)高速氣核區(qū),軸向速度增大。這是由于: (1)渦流工具減小流體的過流截面,氣體流量一定時(shí),橫截面積的減小必然造成流速的增大;在極高的速度下,液體獲得離心加速度被甩至管壁,氣體得以占據(jù)中心流道,形成氣核;(2)由于氣體具有較強(qiáng)的壓縮性,液體被甩至管壁造成中心區(qū)壓力降低,氣體膨脹,井筒過流截面不變,氣體流速增大。

    圖4 流體流動(dòng)軌跡示意圖Fig. 4 Flow trajectories of fluids

    3.2氣液條件變化時(shí)工具結(jié)構(gòu)參數(shù)優(yōu)化

    Optimization of structural parameters under changing gas/liquid conditions

    氣液條件取值:速度取2 m/s、4 m/s和6 m/s;液氣比取1∶50、1∶100、1∶200、1∶1 000和1∶2 000。進(jìn)行15組正交試驗(yàn),每組選擇L16(45)的正交表進(jìn)行16次模擬計(jì)算,并進(jìn)行極差分析和方差分析。通過極差分析可以得到設(shè)定氣液條件下工具的最優(yōu)結(jié)構(gòu)參數(shù),15種氣液條件下工具最優(yōu)結(jié)構(gòu)參數(shù)見表4。

    通過表4可以看出:槽寬范圍為40~80 mm,最優(yōu)結(jié)構(gòu)參數(shù)取40 mm和45 mm;槽深范圍為5~15 mm,15組氣液條件的最優(yōu)結(jié)構(gòu)參數(shù)均取5 mm;導(dǎo)程范圍為125~218 mm,最優(yōu)結(jié)構(gòu)參數(shù)取179 mm和218 mm;導(dǎo)程數(shù)范圍為1.5~3倍,15組氣液狀況的最優(yōu)結(jié)構(gòu)參數(shù)均取3倍。這說明較小的槽深和槽寬、較大的導(dǎo)程和導(dǎo)程數(shù)可以起到很好的增速效果。這是因?yàn)殡S著槽深和槽寬的減小,氣液混合物流經(jīng)截面積減小,速度增大;隨著導(dǎo)程和導(dǎo)程數(shù)的增大,螺旋葉片長度增加,流體沿螺旋葉片的加速作用增強(qiáng),流體速度增大,軸向分速度和徑向分速度增大,流體所受的離心力增大,氣液分離效果增強(qiáng),減小了氣液間的碰撞和摩擦,因此提高了中心氣柱的速度。

    通過方差分析判斷一定氣液條件下工具結(jié)構(gòu)參數(shù)對(duì)結(jié)果影響的顯著程度。對(duì)于給定水平α=0.10,由P{F>λ}=0.10查F(3,12)分布表,得到λ=2.61。為了更清楚地表現(xiàn)出工具結(jié)構(gòu)參數(shù)影響的顯著程度,取λ的一半即以λ=1.31為界,在“影響顯著”和“影響不顯著”之間加入“影響較顯著”。15種氣液條件得到結(jié)果匯總見表5。

    表5 不同氣液條件下方差分析結(jié)果Table 5 Optimal structural parameters for tools under different gas/liquid conditions

    通過表5可以看出:隨著氣液條件的改變,工具結(jié)構(gòu)參數(shù)對(duì)攜液效果影響的顯著程度隨之變化;液氣比相同時(shí),隨著入口速度的增大,導(dǎo)程對(duì)氣相出口峰值速度的影響逐漸不顯著,槽寬對(duì)氣相出口峰值速度的影響逐漸顯著;速度相同時(shí),隨著液氣比的增大,導(dǎo)程數(shù)對(duì)氣相出口峰值速度的影響逐漸顯著;槽深的變化對(duì)氣相出口峰值速度有顯著影響。

    4 結(jié)論

    Conclusions

    (1)應(yīng)用計(jì)算流體動(dòng)力學(xué)軟件Fluent和正交實(shí)驗(yàn)設(shè)計(jì)法,研究出一套采用正交數(shù)值試驗(yàn)對(duì)渦流工具結(jié)構(gòu)參數(shù)進(jìn)行優(yōu)化設(shè)計(jì)的方法。算例結(jié)果表明該方法得到的參數(shù)組合能達(dá)到最好攜液效果,結(jié)果合理可靠。

    (2)利用Fluent軟件模擬計(jì)算了15組氣液情況,通過極差分析得到不同氣液條件下工具最優(yōu)結(jié)構(gòu)參數(shù)組合,結(jié)果表明:槽深取5 mm,槽寬取40~45 mm,導(dǎo)程數(shù)取3倍,導(dǎo)程取172~219 mm時(shí),氣液混合物經(jīng)過渦流工具后氣相峰值速度最大,攜液效果最好。

    (3)通過方差分析研究了15組氣液情況下工具結(jié)構(gòu)參數(shù)對(duì)出口峰值速度影響顯著程度,結(jié)果表明:液氣比較大時(shí),較大的導(dǎo)程數(shù)增速效果更明顯;入口速度較大時(shí),較小的槽寬增速效果更明顯,而較大的導(dǎo)程增速效果并不明顯;較小的槽深始終具有明顯的增速效果。

    References:

    [1] 楊旭東,衛(wèi)亞明,肖述琴,于志剛,陳勇.井下渦流工具排水采氣在蘇里格氣田探索研究[J].鉆采工藝,2013,36(6):125-127. YANG Xudong, WEI Yaming, XIAO Shuqin, YU Zhigang, CHEN Yong. Downhole vortex tool drainage gas recovery research and exploration in Sulige gas field[J]. Drilling & Production Technology, 2013, 36(6): 125-127.

    [2] 楊銀山. 南八仙氣田排水采氣工藝先導(dǎo)性試驗(yàn)及效果研究[J]. 鉆采工藝,2013,36(4):44-47. YANY Yinshan. Draniage gas recovery technology test and effect analysis in Nanbaxian gas fieled[J]. Drilling & Production Technology, 2013, 36(4): 44-47.

    [3] 楊濤,余淑明,楊樺,李雋,李楠,曹光強(qiáng),王云.氣井渦流排水采氣新技術(shù)及其應(yīng)用[J].天然氣工業(yè),2012,32(8):63-66. YANY Tao, YU Shuming, YANG Hua, LI Jun, LI Nan,CAO Guangqing, WANG Yun. A new technology of vortex dewatering gas recovery in gas wells and its application[J], Natural Gas Industry, 2012, 32(8): 63-66.

    [4] 朱慶,張俊杰,謝飛,曾誠. 渦流排水采氣技術(shù)在四川氣田的應(yīng)用[J]. 天然氣技術(shù)與經(jīng)濟(jì),2013,7(1):37-39. ZHU Qing, ZHANG Junjie, XIE Fei, ZENG Cheng. Application of eddy drainage gas recovery to gasfields in Sichuan gas field[J]. Natural Gas Technology and Economy, 2013,7(1): 37-39.

    [5] 趙天沛.海上氣田氣井排液采氣技術(shù)應(yīng)用及發(fā)展淺談[J].海洋石油,2011,31(2):75-78. ZHAO Tianpei. Practice and direction of gas recovery by unloading in offshore gas fields[J]. Offshore Oil, 2011,31(2): 75-78.

    [6] 張翠婷.渦流工具在天然氣井井底排液中的應(yīng)用[D].大慶:東北石油大學(xué),2012. ZHANG Cuiting. The application of vortex tools in bottom’ s drainage of the natural gas wells[D]. Daqing: Northeast University of Petroleum, 2012.

    [7] 馮翠菊,王春生,張黌.天然氣井下渦流工具排液效果影響因素分析[J]. 石油機(jī)械,2013,41(1):78-81. FENG Cuiju, WANG Chunsheng, ZHANG Hong. Influencing factor analysis of the liquid discharge effect of downhole vortex tool in natural gas wells[J]. China Petroleum Machinery, 2013, 41(1): 78-81.

    [8] 吳丹,顏廷俊,謝雙喜,姜豐華.氣井渦流排液采氣工具參數(shù)仿真及結(jié)構(gòu)優(yōu)化[J].石油機(jī)械,2014,42(8):111-115. WU Dan, YAN Tingjun, XIE Shuangxi, JIANG Fenghua. Numerical simulation and structural optimization of vortex tool for liquid discharge in gas wells[J]. China Petroleum Machinery, 2014, 42(8): 111-115.

    [9] 張春,金大權(quán),王晉,張金波.蘇里格氣田井下渦流排水采氣工藝研究[J].天然氣技術(shù)與經(jīng)濟(jì),2012,6(5):45-48. ZHANG Chun, JING Daquan, WANG Jin, ZHANG Jinbo. Research of downhole eddy drainage gas recovery technology in Sulige gas field[J]. Natural Gas Technology and Economy, 2012, 6(5): 45-48.

    [10] 李雋,李楠,李佳宜,王云,曹光強(qiáng).渦流排水采氣技術(shù)數(shù)值模擬研究[J].石油鉆采工藝,2013,35(6):65-68. LI Jun, LI Nan, LI Jiayi, WANG Yun, CAO Guangqiang. Numerical simulation research on eddy current drainage gas recovery technology[J]. Oil Drilling & Production Technology, 2013, 35(6): 65-68.

    [11] 陳德春,姚亞,韓昊,付剛,宋天驕,謝雙喜.氣井渦流攜液機(jī)理和攜液效率數(shù)值模擬研究[J].石油機(jī)械,2015,43(9):91-94. CHEN Dechun, YAO Ya, HAN Hao, FU Gang, SONG Tianjiao, XIE Shuangxi. Numerical simulation on fluidcarrying mechanism and efficiency of vortex flow in gas wells[J]. China Petroleum Machinery, 2015, 43(9): 91-94.

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

    〔編輯 朱 偉〕

    Optimization of structural parameters for fluid-carrying swirl tool in gas wells

    CHEN Dechun1, FU Gang1, HAN Hao1, YAO Ya1, SONG Tianjiao1, XIE Shuangxi2
    1. College of Petroleum Engineering, China Uniνersity of Petroleum, Qingdao, Shandong 266580, China;2. CNOOC EnerTech-Drilling & Production Co., Tianjin, 300457, China

    A study was made on structural optimization of swirl tool to enhance operation efficiency of the swirl fluid-discharging and gas-producing tool in gas wells. AutoCAD and Fluent fluid simulation software were used in the study to build a flow field model for gas/liquid in swirl tool. In addition, orthogonal test was performed to clarify the program for simulation and calculation. Through variance and range analysis, impacts of four structural parameters (i.e. lead, number of leads, trough depth and trough width) on peak velocity at outlet of the tool were determined together with optimal structural parameters for such tool. The results show that minor trough depth and width, together with large lead and number of leads can effectively enhance velocities. With trough depth of 5 mm, trough width of 40-45 mm, lead number of 3 and leads of 172-219 mm, the maximum velocity of gas from gas/liquid mixture through swirl tool with the most desirable fluid-carrying capacity can be obtained. Impacts of structural parameters of such tool on the maximum outlet velocity may vary with conditions of gas/liquid. With relatively high liquid/gas ratio, higher number of leads may effectively enhance the velocity. With relatively high velocity at inlet, smaller trough width may obviously increase the velocity. Besides, larger leads may have

    gas well; discharging and production technique; swirl tool; optimization of structural parameter; liquid/gas ratio; numerical simulation; orthogonal test; sensitivity analysis

    付剛(1992-),在讀碩士研究生,研究方向:油氣田開發(fā)工程。通訊地址:(266580)山東省青島市經(jīng)濟(jì)技術(shù)開發(fā)區(qū)長江西路66號(hào)中國石油大學(xué)(華東)工科樓B座439。E-mail:1184404933@qq.cominsignificant impacts on increase in velocity. Smaller trough depth may constantly display obvious impacts on enhancement of velocity. The study results may provide necessary references for selecting optimal swirl tools and promoting performances of such tools on site.

    TE934

    A

    1000 - 7393( 2016 ) 03 - 0400- 05

    10.13639/j.odpt.2016.03.024

    CHEN Dechun, FU Gang, HAN Hao, YAO Ya , SONG Tianjiao, XIE Shuangxi. Optimization of structural parameters for fluid-carrying swirl tool in gas wells[J]. Oil Drilling & Production Technology, 2016, 38(3): 400-404.

    中海油能源發(fā)展股份有限公司科研項(xiàng)目“渦流工具攜液影響因素分析”(編號(hào):GC2014ZC2916)。

    陳德春(1969-),博士,教授,從事采油工程理論與技術(shù)的研究與教學(xué)。通訊地址:(266580)山東省青島市經(jīng)濟(jì)技術(shù)開發(fā)區(qū)長江西路66號(hào)中國石油大學(xué)(華東)工科樓B座412。E-mail:chendc@upc.edu.cn

    引用格式:陳德春,付剛,韓昊,姚亞,宋天驕,謝雙喜.氣井?dāng)y液用渦流工具結(jié)構(gòu)參數(shù)優(yōu)化[J].石油鉆采工藝,2016,38(3):400-404.

    猜你喜歡
    導(dǎo)程攜液氣液
    水平井組合管柱排水采氣研究及應(yīng)用:以鄂北D氣田為例
    微重力下兩相控溫型儲(chǔ)液器內(nèi)氣液界面仿真分析
    氣液分離罐液位計(jì)接管泄漏分析
    徐深氣田氣井臨界攜液影響因素研究
    特種油氣藏(2019年3期)2019-07-26 09:36:46
    論四頭雙導(dǎo)程蝸桿車削掛輪的選配
    論四頭雙導(dǎo)程蝸桿車削掛輪的選配
    普通車床車削多線螺紋的工藝探究
    CO2 驅(qū)低液量高氣液比井下氣錨模擬與優(yōu)化
    螺旋側(cè)板的導(dǎo)程對(duì)VIV影響的數(shù)值模擬
    氣井多液滴攜液模型實(shí)驗(yàn)研究
    斷塊油氣田(2014年6期)2014-03-11 15:33:59
    老熟妇乱子伦视频在线观看| 精品一区二区三区四区五区乱码| 18禁美女被吸乳视频| 亚洲av日韩精品久久久久久密| 国产主播在线观看一区二区| 啦啦啦免费观看视频1| ponron亚洲| 久久久久久免费高清国产稀缺| 国产又爽黄色视频| 中文字幕另类日韩欧美亚洲嫩草| 后天国语完整版免费观看| 中出人妻视频一区二区| 亚洲色图综合在线观看| 岛国在线观看网站| 波多野结衣高清无吗| 亚洲成人国产一区在线观看| 久久精品91无色码中文字幕| 美女扒开内裤让男人捅视频| 99国产精品一区二区三区| 99久久久亚洲精品蜜臀av| 午夜a级毛片| 亚洲人成77777在线视频| 嫁个100分男人电影在线观看| 亚洲一区高清亚洲精品| 精品国产国语对白av| 久久香蕉精品热| 一级a爱片免费观看的视频| 亚洲熟妇熟女久久| 亚洲五月婷婷丁香| 日本三级黄在线观看| 9191精品国产免费久久| 成人av一区二区三区在线看| 在线免费观看的www视频| 最近最新中文字幕大全免费视频| 狂野欧美激情性xxxx| 日日摸夜夜添夜夜添小说| 又紧又爽又黄一区二区| 一a级毛片在线观看| 国产精品爽爽va在线观看网站 | 久热爱精品视频在线9| 久久国产乱子伦精品免费另类| 国产一区二区激情短视频| 91国产中文字幕| 在线播放国产精品三级| 麻豆成人av在线观看| 国产不卡一卡二| 亚洲九九香蕉| 免费不卡黄色视频| 欧美乱色亚洲激情| 亚洲国产毛片av蜜桃av| 亚洲欧美一区二区三区黑人| 高清av免费在线| 老司机深夜福利视频在线观看| 国产深夜福利视频在线观看| 久久中文字幕人妻熟女| 最新在线观看一区二区三区| 韩国av一区二区三区四区| 国产深夜福利视频在线观看| 免费观看人在逋| 久久人人精品亚洲av| www.自偷自拍.com| 深夜精品福利| 日韩有码中文字幕| 丁香六月欧美| 91成人精品电影| 成人亚洲精品一区在线观看| 欧美人与性动交α欧美精品济南到| 中出人妻视频一区二区| 日本撒尿小便嘘嘘汇集6| 久久人人97超碰香蕉20202| 色在线成人网| 国产av一区二区精品久久| 亚洲久久久国产精品| 少妇裸体淫交视频免费看高清 | 99热只有精品国产| 男女之事视频高清在线观看| 黄片大片在线免费观看| 精品一区二区三区四区五区乱码| 91麻豆av在线| 夜夜躁狠狠躁天天躁| 国产成人精品久久二区二区免费| 水蜜桃什么品种好| 少妇粗大呻吟视频| 国产亚洲精品第一综合不卡| 亚洲片人在线观看| 精品乱码久久久久久99久播| 久久久国产精品麻豆| 亚洲精品中文字幕一二三四区| 久久性视频一级片| 国产麻豆69| 免费久久久久久久精品成人欧美视频| 9色porny在线观看| 美国免费a级毛片| 在线观看免费视频日本深夜| 自线自在国产av| 婷婷六月久久综合丁香| 亚洲久久久国产精品| av视频免费观看在线观看| 日韩中文字幕欧美一区二区| 亚洲三区欧美一区| 久久精品国产清高在天天线| 欧美色视频一区免费| 久久国产精品人妻蜜桃| 天天躁夜夜躁狠狠躁躁| 天天添夜夜摸| 中出人妻视频一区二区| 亚洲午夜精品一区,二区,三区| 国产不卡一卡二| 老鸭窝网址在线观看| 久久伊人香网站| 亚洲专区国产一区二区| 亚洲熟妇中文字幕五十中出 | 天天躁夜夜躁狠狠躁躁| 亚洲人成77777在线视频| 国产麻豆69| 一本大道久久a久久精品| 亚洲欧美日韩无卡精品| av超薄肉色丝袜交足视频| 天堂影院成人在线观看| 99久久久亚洲精品蜜臀av| 亚洲av日韩精品久久久久久密| 97超级碰碰碰精品色视频在线观看| 一本综合久久免费| 一区二区日韩欧美中文字幕| 欧美日韩视频精品一区| 99久久99久久久精品蜜桃| 在线视频色国产色| 国产精品久久视频播放| 国产精品野战在线观看 | 狂野欧美激情性xxxx| 美女国产高潮福利片在线看| 国产成人系列免费观看| 亚洲中文av在线| 精品免费久久久久久久清纯| 久久久精品国产亚洲av高清涩受| 久久久国产一区二区| 精品国产乱子伦一区二区三区| 国产亚洲精品一区二区www| 欧美精品亚洲一区二区| 国产亚洲精品第一综合不卡| 亚洲成av片中文字幕在线观看| 少妇裸体淫交视频免费看高清 | 中国美女看黄片| 黄色片一级片一级黄色片| 精品久久久久久成人av| 成人永久免费在线观看视频| 身体一侧抽搐| 两个人免费观看高清视频| 久久久久久久午夜电影 | 亚洲专区国产一区二区| 男女午夜视频在线观看| 日本精品一区二区三区蜜桃| 91av网站免费观看| 搡老乐熟女国产| 妹子高潮喷水视频| 中文字幕人妻丝袜制服| 国产精品一区二区在线不卡| 亚洲精品在线观看二区| 在线看a的网站| 91大片在线观看| 国产一区二区三区视频了| videosex国产| 少妇的丰满在线观看| 老司机在亚洲福利影院| 一区福利在线观看| 女警被强在线播放| 国产乱人伦免费视频| 日本黄色日本黄色录像| 亚洲国产中文字幕在线视频| 天天躁夜夜躁狠狠躁躁| 一边摸一边抽搐一进一小说| 9191精品国产免费久久| 在线免费观看的www视频| 精品一区二区三区四区五区乱码| 亚洲欧美激情在线| 91成年电影在线观看| 中出人妻视频一区二区| 纯流量卡能插随身wifi吗| 欧美激情高清一区二区三区| 高清在线国产一区| 日本五十路高清| a级毛片在线看网站| 91老司机精品| av中文乱码字幕在线| 99国产精品一区二区蜜桃av| 久久久久久久久中文| 美女午夜性视频免费| 亚洲五月婷婷丁香| 高清欧美精品videossex| 他把我摸到了高潮在线观看| 国产麻豆69| 欧美精品亚洲一区二区| 99国产综合亚洲精品| 日本三级黄在线观看| 午夜精品国产一区二区电影| a级毛片在线看网站| 黄色毛片三级朝国网站| 亚洲第一av免费看| 免费在线观看黄色视频的| 国产精品乱码一区二三区的特点 | 国产黄色免费在线视频| 99国产精品一区二区三区| 日本精品一区二区三区蜜桃| 午夜成年电影在线免费观看| 日本五十路高清| 这个男人来自地球电影免费观看| 午夜福利影视在线免费观看| 久久精品影院6| 欧美一级毛片孕妇| 夜夜看夜夜爽夜夜摸 | 欧美乱妇无乱码| 免费av中文字幕在线| 一边摸一边抽搐一进一出视频| 精品人妻1区二区| 欧美激情极品国产一区二区三区| 99久久国产精品久久久| 啦啦啦在线免费观看视频4| 中文亚洲av片在线观看爽| 久久狼人影院| 中亚洲国语对白在线视频| 老司机午夜福利在线观看视频| 麻豆一二三区av精品| 妹子高潮喷水视频| 日本精品一区二区三区蜜桃| 国产av又大| 国产精品一区二区在线不卡| 国产亚洲精品久久久久久毛片| 日韩 欧美 亚洲 中文字幕| 丝袜美腿诱惑在线| 欧美日韩亚洲综合一区二区三区_| 成人亚洲精品一区在线观看| 国产成人欧美在线观看| 热99国产精品久久久久久7| 久久国产乱子伦精品免费另类| 黄色视频,在线免费观看| 国产成人一区二区三区免费视频网站| 在线观看www视频免费| 极品教师在线免费播放| 看免费av毛片| 国产不卡一卡二| 黄色片一级片一级黄色片| 侵犯人妻中文字幕一二三四区| 淫秽高清视频在线观看| 国产高清国产精品国产三级| 精品人妻1区二区| 久久精品亚洲av国产电影网| av免费在线观看网站| 久久精品人人爽人人爽视色| 搡老乐熟女国产| 天天躁狠狠躁夜夜躁狠狠躁| 视频区图区小说| 在线观看免费视频日本深夜| 欧美日韩乱码在线| 淫妇啪啪啪对白视频| 精品国产一区二区三区四区第35| 亚洲欧美激情综合另类| 国产亚洲精品久久久久久毛片| 国产亚洲精品久久久久5区| 精品久久久久久久久久免费视频 | 久久久精品国产亚洲av高清涩受| 999久久久精品免费观看国产| 午夜亚洲福利在线播放| 热re99久久国产66热| 黄片小视频在线播放| 女人被狂操c到高潮| 成人av一区二区三区在线看| 大香蕉久久成人网| 麻豆成人av在线观看| 中文字幕人妻丝袜制服| 久久精品国产综合久久久| 亚洲久久久国产精品| 国产成人啪精品午夜网站| 日韩精品青青久久久久久| 日韩高清综合在线| 日韩免费高清中文字幕av| 99在线人妻在线中文字幕| 免费久久久久久久精品成人欧美视频| 女人被躁到高潮嗷嗷叫费观| 国产精品免费一区二区三区在线| 一边摸一边做爽爽视频免费| 日韩 欧美 亚洲 中文字幕| 一区二区三区精品91| 日日爽夜夜爽网站| 黄色 视频免费看| 国产黄a三级三级三级人| 如日韩欧美国产精品一区二区三区| 午夜免费鲁丝| 老熟妇乱子伦视频在线观看| 久久精品国产亚洲av香蕉五月| 国产黄色免费在线视频| 中文字幕最新亚洲高清| 国产无遮挡羞羞视频在线观看| 黄色视频不卡| av福利片在线| 麻豆一二三区av精品| 精品久久久久久,| 亚洲七黄色美女视频| av网站在线播放免费| 国产精品一区二区精品视频观看| 岛国视频午夜一区免费看| 国产极品粉嫩免费观看在线| 高清欧美精品videossex| 午夜福利一区二区在线看| 99在线视频只有这里精品首页| 美国免费a级毛片| 人人妻人人爽人人添夜夜欢视频| 色婷婷av一区二区三区视频| 一二三四在线观看免费中文在| 国产三级在线视频| 老司机靠b影院| 天天添夜夜摸| 91老司机精品| 国产成人系列免费观看| av福利片在线| 一边摸一边抽搐一进一出视频| 亚洲av美国av| 如日韩欧美国产精品一区二区三区| 成人永久免费在线观看视频| avwww免费| 在线观看日韩欧美| 色播在线永久视频| 亚洲av成人不卡在线观看播放网| 国产伦人伦偷精品视频| 可以免费在线观看a视频的电影网站| 丰满饥渴人妻一区二区三| 久久久久久久精品吃奶| 亚洲一区中文字幕在线| 亚洲av第一区精品v没综合| 欧美日韩国产mv在线观看视频| 老汉色av国产亚洲站长工具| 国产99久久九九免费精品| 国产欧美日韩一区二区三| 亚洲欧美精品综合久久99| 久久人人精品亚洲av| 午夜老司机福利片| 亚洲全国av大片| 欧美人与性动交α欧美精品济南到| 欧美国产精品va在线观看不卡| 国产精品久久久人人做人人爽| 国产成人欧美在线观看| 国产成人系列免费观看| 成熟少妇高潮喷水视频| 国产真人三级小视频在线观看| 欧美中文日本在线观看视频| 一级a爱片免费观看的视频| 亚洲一区高清亚洲精品| 美女大奶头视频| 99久久国产精品久久久| 99国产极品粉嫩在线观看| 无人区码免费观看不卡| 精品人妻1区二区| 免费在线观看亚洲国产| 一进一出抽搐gif免费好疼 | 女同久久另类99精品国产91| 午夜免费激情av| 欧美精品亚洲一区二区| www.www免费av| 老司机在亚洲福利影院| 视频在线观看一区二区三区| 天堂俺去俺来也www色官网| 嫁个100分男人电影在线观看| 国产亚洲av高清不卡| 国产黄色免费在线视频| 亚洲精品在线美女| av福利片在线| 亚洲精品在线观看二区| 亚洲成国产人片在线观看| av中文乱码字幕在线| 亚洲成a人片在线一区二区| 精品人妻1区二区| 日韩一卡2卡3卡4卡2021年| 人人妻,人人澡人人爽秒播| 很黄的视频免费| 人人妻,人人澡人人爽秒播| 日韩一卡2卡3卡4卡2021年| av天堂在线播放| 免费不卡黄色视频| 亚洲第一av免费看| 国产精品 欧美亚洲| 丰满的人妻完整版| av天堂久久9| 精品乱码久久久久久99久播| 免费观看精品视频网站| 国产一区二区激情短视频| 国产精品野战在线观看 | 亚洲一区中文字幕在线| 国产精品久久视频播放| 激情视频va一区二区三区| 午夜免费激情av| 精品卡一卡二卡四卡免费| 欧美老熟妇乱子伦牲交| 国产日韩一区二区三区精品不卡| 欧美黑人欧美精品刺激| 脱女人内裤的视频| 一个人免费在线观看的高清视频| 国产熟女xx| 99精国产麻豆久久婷婷| 无限看片的www在线观看| 亚洲av第一区精品v没综合| 欧美日韩瑟瑟在线播放| av超薄肉色丝袜交足视频| 女同久久另类99精品国产91| 后天国语完整版免费观看| 久久久久久人人人人人| 亚洲成a人片在线一区二区| 久久国产乱子伦精品免费另类| 久久草成人影院| 搡老岳熟女国产| 老司机午夜福利在线观看视频| 99re在线观看精品视频| 中文字幕最新亚洲高清| 日韩三级视频一区二区三区| 欧美日韩国产mv在线观看视频| 91av网站免费观看| 97人妻天天添夜夜摸| 色精品久久人妻99蜜桃| 日本黄色视频三级网站网址| 国产极品粉嫩免费观看在线| 亚洲自拍偷在线| 人妻久久中文字幕网| 日韩精品免费视频一区二区三区| 正在播放国产对白刺激| 在线天堂中文资源库| 色综合婷婷激情| 中文字幕人妻丝袜制服| 一二三四在线观看免费中文在| 大香蕉久久成人网| 好男人电影高清在线观看| 亚洲成a人片在线一区二区| 亚洲第一av免费看| av天堂久久9| 亚洲精品美女久久av网站| 看黄色毛片网站| 欧美日韩乱码在线| 一级黄色大片毛片| 99精品久久久久人妻精品| 精品人妻1区二区| 欧美成人免费av一区二区三区| 男女下面进入的视频免费午夜 | 一级,二级,三级黄色视频| 天堂俺去俺来也www色官网| 别揉我奶头~嗯~啊~动态视频| 18禁观看日本| 十分钟在线观看高清视频www| 99re在线观看精品视频| xxx96com| 精品久久久久久久毛片微露脸| 超碰成人久久| 国产精品久久视频播放| 成人av一区二区三区在线看| 欧美丝袜亚洲另类 | 久久久精品欧美日韩精品| 欧美日韩视频精品一区| 久久久精品欧美日韩精品| 脱女人内裤的视频| 中文字幕人妻丝袜一区二区| 精品第一国产精品| 99热只有精品国产| 两个人看的免费小视频| 另类亚洲欧美激情| 国产精品 欧美亚洲| 精品国产美女av久久久久小说| 美女扒开内裤让男人捅视频| 在线视频色国产色| 欧美黄色淫秽网站| 97碰自拍视频| 亚洲人成电影免费在线| 亚洲 欧美 日韩 在线 免费| 免费人成视频x8x8入口观看| 中文字幕高清在线视频| 天堂√8在线中文| 亚洲在线自拍视频| 99国产精品免费福利视频| 人人澡人人妻人| 另类亚洲欧美激情| 亚洲欧美精品综合久久99| 日本五十路高清| 美女 人体艺术 gogo| 悠悠久久av| 桃红色精品国产亚洲av| 无限看片的www在线观看| 午夜福利免费观看在线| 午夜免费鲁丝| 精品乱码久久久久久99久播| www.精华液| 精品一区二区三卡| 亚洲成a人片在线一区二区| 自线自在国产av| 欧美人与性动交α欧美精品济南到| 激情视频va一区二区三区| 日韩大码丰满熟妇| 老司机亚洲免费影院| 18禁美女被吸乳视频| 一本综合久久免费| 亚洲人成77777在线视频| 在线观看免费高清a一片| 亚洲精品美女久久久久99蜜臀| 成人18禁高潮啪啪吃奶动态图| 国产av一区二区精品久久| 国产成人啪精品午夜网站| 国产一区在线观看成人免费| 国产亚洲欧美98| 久久欧美精品欧美久久欧美| 三级毛片av免费| 91字幕亚洲| 成人黄色视频免费在线看| 亚洲成人久久性| 久99久视频精品免费| 999久久久国产精品视频| 久9热在线精品视频| 99久久久亚洲精品蜜臀av| 国产单亲对白刺激| 欧美老熟妇乱子伦牲交| 国产熟女xx| 日韩高清综合在线| 免费看a级黄色片| 国产高清videossex| 90打野战视频偷拍视频| 亚洲欧美精品综合一区二区三区| 中文字幕人妻熟女乱码| 日韩成人在线观看一区二区三区| 高清在线国产一区| 久久久久久大精品| 亚洲午夜精品一区,二区,三区| 亚洲色图 男人天堂 中文字幕| 久久久久亚洲av毛片大全| 欧美日韩亚洲国产一区二区在线观看| 亚洲精品美女久久久久99蜜臀| 成人永久免费在线观看视频| 超色免费av| 国产深夜福利视频在线观看| 色尼玛亚洲综合影院| 欧美日本中文国产一区发布| 色婷婷久久久亚洲欧美| 日韩高清综合在线| 亚洲精品一区av在线观看| 欧洲精品卡2卡3卡4卡5卡区| 99在线人妻在线中文字幕| 日韩国内少妇激情av| 国产成人影院久久av| 精品人妻1区二区| 亚洲精品国产色婷婷电影| 90打野战视频偷拍视频| 一级a爱视频在线免费观看| 国产国语露脸激情在线看| 国产真人三级小视频在线观看| 夜夜躁狠狠躁天天躁| 黄频高清免费视频| 国产三级在线视频| 亚洲情色 制服丝袜| 黄频高清免费视频| 99久久久亚洲精品蜜臀av| 欧美另类亚洲清纯唯美| 成人免费观看视频高清| 亚洲色图av天堂| 日韩有码中文字幕| av视频免费观看在线观看| 亚洲av片天天在线观看| 久久久久九九精品影院| 最好的美女福利视频网| 久久影院123| 久久人妻熟女aⅴ| 免费在线观看影片大全网站| 精品午夜福利视频在线观看一区| 19禁男女啪啪无遮挡网站| 欧美精品亚洲一区二区| 日韩一卡2卡3卡4卡2021年| 欧美激情 高清一区二区三区| 又黄又粗又硬又大视频| 久久香蕉激情| 婷婷六月久久综合丁香| 91九色精品人成在线观看| 国产精品偷伦视频观看了| 在线观看一区二区三区激情| 一级毛片精品| 99国产精品一区二区蜜桃av| 97碰自拍视频| xxx96com| 69av精品久久久久久| ponron亚洲| 国产高清videossex| 中文字幕另类日韩欧美亚洲嫩草| 老司机靠b影院| 多毛熟女@视频| 极品教师在线免费播放| 欧美亚洲日本最大视频资源| 国产成人啪精品午夜网站| 日韩视频一区二区在线观看| 久久久久久免费高清国产稀缺| 一级片'在线观看视频| videosex国产| 夜夜看夜夜爽夜夜摸 | 亚洲av五月六月丁香网| 精品久久久久久电影网| 亚洲成人国产一区在线观看| 不卡av一区二区三区| bbb黄色大片| 免费av毛片视频| 免费在线观看视频国产中文字幕亚洲| 欧美日韩视频精品一区| 日韩欧美一区视频在线观看| 美国免费a级毛片| 午夜精品久久久久久毛片777| 男女下面进入的视频免费午夜 | 性少妇av在线| 日本精品一区二区三区蜜桃| 亚洲av五月六月丁香网| 两人在一起打扑克的视频| 国产一区二区三区综合在线观看| 亚洲专区国产一区二区| 国产成人欧美在线观看| 国产亚洲欧美在线一区二区| 黑人巨大精品欧美一区二区mp4| 亚洲第一欧美日韩一区二区三区| 视频区欧美日本亚洲|