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

    機械式旋沖工具對水平井送鉆摩阻的影響規(guī)律

    2023-05-30 10:48:04閆炎韓禮紅宋志龍劉永紅楊尚諭王建軍
    石油機械 2023年4期
    關(guān)鍵詞:摩阻水平井

    閆炎 韓禮紅 宋志龍 劉永紅 楊尚諭 王建軍

    摘要:托壓問題是制約水平井滑動鉆進機械鉆速的主要因素之一,使用振動減阻工具是水平井減阻最直接有效的方式。為分析振動減阻工具在水平井中的減阻效果,針對一種螺桿驅(qū)動的機械式旋轉(zhuǎn)沖擊鉆井工具,對工具在不同鉆進工況下的鉆壓波動、位移變化進行了模擬計算,解釋了低頻沖擊工具可以實現(xiàn)振動減阻的根本原因。計算結(jié)果表明:旋沖工具降低摩阻需要借助鉆壓;激勵頻率越大,鉆頭處的鉆壓波動越劇烈;存在最佳的激勵頻率使得鉆頭處鉆壓整體處于最高水平;旋沖工具工作時,鉆頭處的鉆壓波動不僅包含激勵本身的頻率,還會包含激勵頻率的2倍、3倍等更高倍數(shù)的頻率。研究結(jié)果可以為旋沖工具工作參數(shù)的優(yōu)化,以及工具在水平井中實現(xiàn)最佳提速效果提供理論指導(dǎo)。

    關(guān)鍵詞:水平井;機械式旋沖工具;摩阻;旋轉(zhuǎn)沖擊;沖擊頻率;鉆壓波動

    Weightonbit (WOB) congestion is one of the main factors restricting the penetration rate of sliding drilling in horizontal wells, and the use of vibration friction reduction tools is the most direct and effective way to reduce friction in horizontal wells.In order to analyze the friction reduction effect of vibration friction reduction tool in horizontal wells, the WOB fluctuation and displacement variation of a mechanical rotary percussion drilling tool driven by positive displacement motor (PDM) under different drilling conditions was calculated by simulation, and the reason that low frequency percussion tool can achieve vibration friction reduction was explained.The calculation results show that the rotary percussion tool needs to use WOB to reduce the friction;the greater the excitation frequency, the more severe the WOB fluctuation is;and there is an optimal excitation frequency to make the WOB at the highest level as a whole;moreover, when the rotary percussion tool works, the WOB fluctuation not only includes the frequency of the excitation itself, but also includes the frequency of 2 times and 3 times of the excitation frequency.The research results provide theoretical guidance for the optimization of working parameters of the rotary percussion tool and the realization of optimal penetration rate enhancement of the tool in horizontal wells.

    horizontal well;mechanical rotary percussion tool;friction;rotary percussion;percussion frequency;WOB fluctuation

    0 引 言

    鉆井作業(yè)中,托壓問題是制約定向井、水平井滑動鉆進機械鉆速的關(guān)鍵因素之一[1-4]。鉆進過程中,鉆柱與井壁間的摩擦力導(dǎo)致鉆壓無法完全施加到鉆頭上,這會使鉆頭的破巖效率降低。工程上通常通過降低鉆柱與井壁之間的摩擦因數(shù)[5-9]、優(yōu)化井眼軌跡[10-11]、使用降摩減阻工具[12-15]等方式來緩解托壓問題。

    在上述方法中,使用振動減阻工具是水平井減阻最直接有效的治理方式,不同石油公司與科研院所根據(jù)不同的工作原理研制了各種不同振動方式的振動減阻工具[16-19]。在這些振動減阻工具中,以螺桿驅(qū)動的減阻工具居多。為此,筆者以一種螺桿驅(qū)動的機械式旋轉(zhuǎn)沖擊鉆井工具[20]為例,進行水平井鉆進過程中送鉆摩阻的分析。機械式旋轉(zhuǎn)沖擊鉆井工具通過螺桿帶動下凸輪旋轉(zhuǎn),下凸輪帶動嚙合的上凸輪及整個鉆柱做軸向運動,完成對下凸輪及鉆頭的沖擊作用。在直井中使用該工具,可以為鉆頭提供沖擊力,提高鉆頭破巖效率。在定向井、水平井中使用該工具,可以有效解決托壓問題,并在現(xiàn)場應(yīng)用中取得了很好的提速效果[21]。但工具在水平井中的減阻效果及影響因素還需要進一步厘清。

    間隙元理論[22]、分段計算方法[23]、非線性有限單元法等方法都是研究水平井振動減阻技術(shù)的常用方法,這些方法大多忽略井底鉆頭,將鉆壓的施加、井壁對鉆柱的摩擦力、鉆柱底部的激勵等均簡化假設(shè)為一定的外部邊界條件;簡化后雖然可獲得鉆壓、摩擦因數(shù)、激勵頻率等參數(shù)對減阻的影響效果,但由于模型過于簡化,導(dǎo)致某些重要參數(shù)無法通過計算結(jié)果直觀地反應(yīng)出來,如鉆進過程中使用振動減阻工具前后鉆頭處的鉆壓增幅,鉆壓波動的幅值和頻率是否變化等?;诖?,筆者針對螺桿驅(qū)動的機械式旋轉(zhuǎn)沖擊鉆井工具工作原理,基于有限元軟件平臺ABAQUS建立了水平井軸向振動減阻模型,模型考慮了鉆柱與井壁的接觸,實際鉆井過程中上部鉆柱的鉆壓作用,并且把鉆頭上方旋沖工具對鉆柱的激勵直接設(shè)置為一個凸輪結(jié)構(gòu),以模擬出更加真實的激勵方式,以期為機械式旋轉(zhuǎn)沖擊鉆井工具在水平井中的減阻作用提供理論參考。

    1 數(shù)值模型

    1.1 幾何模型

    根據(jù)機械式旋轉(zhuǎn)沖擊鉆井工具舉升鉆柱的工作原理(見圖1)設(shè)計了仿真的幾何模型,如圖2所示。由圖2可以看出,模型包含井筒、鉆柱、鉆頭、凸輪、凸輪軸,鉆柱的軸線與井筒中心線重合。由于鉆柱在井筒內(nèi)沒有旋轉(zhuǎn)運動,且在仿真過程中只受重力、左端施加的鉆壓以及右端施加的位移激勵作用,所以在鉆柱與凸輪之間加了一個限位塊,用來約束鉆柱與鉆頭的橫向相對運動,且限位塊可將凸輪的旋轉(zhuǎn)運動轉(zhuǎn)化為對鉆柱的位移激勵。工具本體使用的是上下嚙合的凸輪,工具對鉆柱的激勵頻率由凸輪齒數(shù)與鉆頭轉(zhuǎn)速共同決定。為了研究方便,有限元模型使用圖2中所示的盤形凸輪,該凸輪繞凸輪軸旋轉(zhuǎn)1周即對鉆柱形成1次激勵,改變激勵頻率只需改變凸輪繞凸輪軸的旋轉(zhuǎn)速度即可。模型各部件的尺寸如下:鉆頭為215.9 mm(81/2 in)六刀翼PDC鉆頭,外徑215.9 mm;139.7 mm(51/2 in)鉆桿長1 000 m;井筒長1 001 m,內(nèi)徑216 mm;凸輪基圓半徑12 mm,升程為8 mm;凸輪軸直徑10 mm。

    施工作業(yè)過程中鉆柱與井壁、鉆柱右端與凸輪均存在接觸碰撞,所以在模擬過程中需要進行接觸設(shè)置。為了模擬鉆柱右端與凸輪的接觸碰撞行為,在鉆柱的右端與凸輪之間加了一個實體限位塊,且凸輪和限位塊均設(shè)置成剛體。因為模型中存在剛體與剛體的接觸行為,所以模擬過程中接觸算法使用罰函數(shù)法。

    在罰接觸算法中,彈簧單元作用在從結(jié)點上的法向接觸力為:

    式中:f為接觸點上的法向接觸力,N;l為從結(jié)點與主面之間的距離,m;k為彈簧剛度,N/m;n為接觸點處的法向單位矢量,1。

    罰函數(shù)法允許物體間發(fā)生一定的穿透。在求解過程中,在每一個時間步中,先檢查各從結(jié)點是否穿透主面,沒有穿透則不作任何處理,否則在該從結(jié)點與被穿透的主表面之間引入一個大小與穿透量及主片剛度成正比的接觸力。

    采用顯式算法求解運動方程時,在每個時間步中先不考慮物體間的接觸和碰撞,完全獨立地對各物體的運動方程進行求解,得到各物體運動狀態(tài)的試探解。如果試探解滿足物體之間的無接觸條件,則在本時間步內(nèi)物體之間沒有發(fā)生接觸,因此物體之間的運動完全獨立,試探解即為真實解;如果試探解不滿足物體之間的無接觸條件,則需要修正試探解,使其滿足接觸界面條件。

    1.3 材料參數(shù)與網(wǎng)格劃分

    由于本文只關(guān)注鉆柱在受到重力、鉆壓以及位移激勵情況下的變形和運動,除鉆柱外其他部件均設(shè)置為剛體。鉆柱的材料為鋼,密度為7 850 kg/m3,彈性模量為210 GPa,泊松比為0.3。由于鉆柱長1 000 m,且要與井筒發(fā)生碰撞,所以使用三維一階鐵木辛柯梁單元B31。仿真過程中,假設(shè)鉆頭、凸輪、凸輪軸、限位塊均不會發(fā)生變形;PDC鉆頭、凸輪和凸輪軸均使用4節(jié)點減縮積分殼單元S4R,并進行剛體約束;限位塊使用三維8節(jié)點減縮積分實體單元C3D8R,并進行剛體約束。模型各部件的網(wǎng)格劃分如圖3所示。

    1.4 分析步、邊界條件及接觸設(shè)置

    實際鉆井過程中,鉆柱一開始就因受重力和上部鉆柱所施加的鉆壓影響而貼于下井壁,且整個鉆柱中都存在軸向力,因此在仿真過程中,為了真實反映實際鉆井工況,仿真過程包含3個顯式動力學(xué)分析步:①對整個鉆柱施加重力,使其貼于下井壁,重力施加時間為65 s。②在鉆柱左端緩慢施加鉆壓。為了不使鉆柱產(chǎn)生過大振動,鉆壓施加時間設(shè)定為65 s。③使凸輪繞凸輪軸旋轉(zhuǎn),凸輪帶動限位塊做周期運動,對鉆柱施加位移激勵,激勵時間設(shè)置為20 s。

    模型中井筒的自由度全部被約束,對限位塊施加剛體約束,鉆柱右端、鉆頭與凸輪軸的自由度全部約束,約束凸輪與凸輪軸的相對自由度,使凸輪只能繞凸輪軸轉(zhuǎn)動。研究在沒有進尺的情況下,旋沖鉆井工具能產(chǎn)生的最大沖擊力與沖擊頻率,以及所施加鉆壓之間的關(guān)系。限位塊與凸輪外表面、限位塊與鉆頭外表面之間設(shè)置面面接觸,且無摩擦因數(shù);在模擬過程中鉆柱與井壁發(fā)生接觸碰撞,因此設(shè)置通用接觸,摩擦因數(shù)設(shè)置為0.3。

    2 旋沖工具對水平井送鉆摩阻的影響

    圖4為在鉆柱左端施加鉆壓100 kN、激勵頻率8 Hz的工況下(分析步②),鉆柱左端的位移與鉆頭處鉆壓隨時間變化的曲線。由圖4可以看出,施加鉆壓后鉆柱左端向右移動了0.086 6 m,鉆頭處的最終鉆壓只有23.3 kN,說明此時井筒內(nèi)出現(xiàn)了托壓。

    圖5為在鉆柱左端施加鉆壓100 kN、激勵頻率8 Hz工況下,旋沖工具工作過程中(分析步③)鉆柱左端位移與鉆頭處鉆壓隨時間變化的曲線。當(dāng)旋沖工具工作后,即凸輪給鉆柱右端一個頻率8 Hz、幅值8 mm的激勵位移,鉆柱左端的位移在5 s內(nèi)從0.086 6 m逐漸增大到0.136 1 m,向右移動了0.049 5 m,之后保持穩(wěn)定并伴隨有微小的波動;而鉆頭處的鉆壓也從23.3 kN逐漸增大,直到保持一個穩(wěn)定的波動狀態(tài),鉆壓谷值為45.5 kN,峰值為148.5 kN,平均值為94.3 kN。由圖4與圖5的參數(shù)對比可知,旋沖工具可以降低水平井中的送鉆摩阻,有效解決托壓問題。

    由圖5可知,鉆壓在135 s之后保持一個穩(wěn)定的波動狀態(tài),因此取135~150 s之間的鉆壓信號進行分析,對其進行快速傅里葉變換,提取鉆壓的頻譜,如圖6所示。由圖6可以看出,在旋沖工具穩(wěn)定工作的時候,鉆壓波動不僅包含凸輪本身對鉆柱的激勵頻率8 Hz,而且還包括16和24 Hz等2倍、3倍的倍頻。由此說明,旋沖工具單一頻率的激勵可以使鉆柱對鉆頭產(chǎn)生多種倍頻的沖擊。

    3 水平井送鉆摩阻影響因素分析

    3.1 激勵頻率對鉆頭處鉆壓幅值的影響

    圖7為鉆頭處鉆壓峰值、谷值和平均值隨激勵頻率的變化規(guī)律展示圖。由圖7a、圖7b可知,相同鉆壓下,隨著旋沖工具激勵頻率的增大,鉆頭處鉆壓峰值逐漸增大,鉆壓谷值逐漸減小,說明激勵頻率的升高加劇了鉆柱底端的振動;而鉆壓平均值隨著激勵頻率的升高先增后減,但均在激勵頻率為8 Hz時最大,分別為94.3和72.8 kN,說明此時鉆頭整體處于高壓狀態(tài)。

    由圖7c可看出,在鉆壓為60 kN工況下,無論激勵頻率為多少,鉆壓的谷值和平均值均接近于0,而鉆壓峰值約為30 kN;說明在鉆壓為60 kN工況下,由于送鉆摩阻的存在,鉆柱左端的鉆壓無法傳遞到鉆頭處,此時旋沖工具無法解決托壓問題,但由于凸輪對鉆柱的激勵作用,鉆頭處依然會產(chǎn)生30 kN左右的沖擊力。

    3.2 激勵頻率對鉆柱左端位移的影響

    圖8為不同激勵頻率時鉆柱左端的位移隨時間變化的規(guī)律。圖9為不同鉆壓下旋沖工具激勵引起的鉆柱左端最大位移與激勵頻率的關(guān)系。由圖8a、圖8b、圖9a、圖9b可知,在相同鉆壓下,鉆柱左端位移隨激勵頻率的升高先增后減,且均在8 Hz時最大,分別為0.049 5和0.049 3 m;說明旋沖工具在水平井的振動減阻作業(yè)中存在最佳工作頻率。

    由圖8c和圖9c可知,在鉆壓60 kN工況下,無論激勵頻率為多少,鉆柱左端位移均無變化,說明此時旋沖工具無法解決托壓問題。這也是圖7c中鉆頭處鉆壓平均值并沒有隨著激勵而增大的原因。

    3.3 激勵頻率對鉆頭處鉆壓波動頻率的影響

    圖10為不同鉆壓條件下不同激勵頻率時鉆頭處鉆壓的頻譜圖。由圖10a、圖10b可知,在100和80 kN鉆壓工況下,鉆頭處的鉆壓波動不僅包含旋沖工具本身的激勵頻率,而且還包含激勵頻率的2倍、3倍等高倍數(shù)頻率,說明低頻的旋沖工具也會帶來高頻沖擊,其幅值甚至不低于低頻沖擊的幅值;由圖10c可知,在60 kN鉆壓工況下,鉆頭處鉆壓波動各頻率的幅值均接近于0,結(jié)合圖7c可知,此時鉆頭處只存在幅值不大的沖擊力作用。

    4 結(jié) 論

    (1)旋沖工具在水平井中減摩阻作用受鉆壓的影響,當(dāng)鉆壓可以傳遞到鉆頭處時,旋沖工具可以有效解決托壓問題;當(dāng)鉆壓無法傳遞到鉆頭處時,旋沖工具無法解決托壓問題,但是其激勵作用會在鉆頭處產(chǎn)生一定的沖擊力。

    (2)鉆壓一定時隨著激勵頻率的增大,鉆頭處鉆壓峰值逐漸增大,鉆壓谷值逐漸減小,鉆壓波動會更加劇烈;鉆壓平均值隨著激勵頻率的增大先升高后降低,存在最佳的激勵頻率使得鉆頭處鉆壓整體處于最高水平。

    (3)在有旋沖工具的工況下,鉆頭處的鉆壓波動不僅包含旋沖工具本身的激勵頻率,還包含2倍、3倍甚至數(shù)倍的高頻,說明低頻的旋沖工具也會帶來高頻沖擊,其幅值甚至不低于低頻沖擊幅值。

    [1] NEWMAN K,BURNETT T,PURSELL J,et al.Modeling the affect of a downhole vibrator[C]∥SPE/ICoTA Coiled Tubing & Well Intervention Conference and Exhibition.The Woodlands,Texas:SPE,2009:SPE 121752-MS.

    [2] GOICOECHEA H E,BUEZAS F S,ROSALES M B.A novel tool to improve the rate of penetration by transferring drilling string vibration energy to hydraulic energy[J].International Journal of Mechanical Sciences,2019,157:98-110.

    [3] CHEN P J,GAO D L,WANG Z H,et al.Study on aggressively working casing string in extendedreach well[J].Journal of Petroleum Science and Engineering,2017,157:604-616.

    [4] 李勇政,陳濤,江川,等.四川盆地磨溪–高石梯區(qū)塊定向鉆井關(guān)鍵技術(shù)[J].石油鉆探技術(shù),2021,49(2):26-31.

    LI Y Z,CHEN T,JIANG C,et al.Key technologies of directional drilling in the MoxiGaoshiti Area of the Sichuan Basin[J].Petroleum Drilling Techniques,2021,49(2):26-31.

    [5] 陳樂亮.水平井鉆井液降摩阻問題[J].鉆井液與完井液,1992,9(4):1-6,75.

    CHEN L L.Problem of lowering borehole friction in horizontal wells[J].Drilling Fluid & Completion Fluid,1992,9(4):1-6,75.

    [6] 趙虎.鉆井液潤滑性評價方法對比研究[J].石油工業(yè)技術(shù)監(jiān)督,2021,37(4):24-27.

    ZHAO H.Comparison of evaluation methods for lubricity of drilling fluid[J].Technology Supervision in Petroleum Industry,2021,37(4):24-27.

    [7] LIU W L,NI H J,WANG P,et al.Investigation on the tribological performance of microdimples textured surface combined with longitudinal or transverse vibration under hydrodynamic lubrication[J].International Journal of Mechanical Sciences,2020,174:105474.

    [8] LIU W L,NI H J,WANG P,et al.Analytical investigation of the friction reduction performance of longitudinal vibration based on the modified elastoplastic contact model[J].Tribology International,2020,146:106237.

    [9] 白冬青.最小曲率法計算中的幾個問題[J].?dāng)鄩K油氣田,2007,14(5):67-69.

    BAI D Q.Several problems on application of minimum curvature method[J].FaultBlock Oil & Gas Field,2007,14(5):67-69.

    [10] 韋龍貴,王贊,陳立強,等.基于熵權(quán)法的大位移井井眼軌跡設(shè)計方案優(yōu)選[J].石油鉆采工藝,2021,43(3):295-301.

    WEI L G,WANG Z,CHEN L Q,et al.Selection of extended reach well trajectory design scheme based on entropy weight method[J].Oil Drilling & Production Technology,2021,43(3):295-301.

    [11] 趙靜.YT25-1S油田大位移井巖屑床清除器優(yōu)化研究[D].北京:中國石油大學(xué)(北京),2016.

    ZHAO J.Study on downhole cuttings bed cleaning tools for extendedreach drilling in YT25-1S oil field[D].Beijing:China University of Petroleum (Beijing),2016.

    [12] 王建龍,鄭鋒,劉學(xué)松,等.井眼清潔工具研究進展及展望[J].石油機械,2018,46(9):18-23.

    WANG J L,ZHENG F,LIU X S,et al.Advances and prospects of well cleaning tools[J].China Petroleum Machinery,2018,46(9):18-23.

    [13] 孔令镕,王瑜,鄒俊,等.水力振蕩減阻鉆進技術(shù)發(fā)展現(xiàn)狀與展望[J].石油鉆采工藝,2019,41(1):23-30.

    KONG L R,WANG Y,ZOU J,et al.Development status and prospect of hydrooscillation drag reduction drilling technology[J].Oil Drilling & Production Technology,2019,41(1):23-30.

    [14] 余志清.降摩阻短節(jié)在定向鉆井及水平鉆井中的應(yīng)用[J].鉆采工藝,1999,22(1):66-67.

    YU Z Q.Application of antifrictional short section in directional drilling and horizontal drilling[J].Drilling & Production Technology,1999,22(1):66-67.

    [15] CHEN J K,LIAO H L,ZHANG Y T,et al.A torsionalaxial vibration analysis of drill string endowed with kinematic coupling and stochastic approach[J].Journal of Petroleum Science and Engineering,2021,198:108157.

    [16] 張輝,吳仲華,蔡文軍.水力振蕩器的研制及現(xiàn)場試驗[J].石油機械,2014,42(6):12-15.

    ZHANG H,WU Z H,CAI W J.Development and field testing of hydraulic oscillator[J].China Petroleum Machinery,2014,42(6):12-15.

    [17] 張會增,管志川,劉永旺,等.基于旋轉(zhuǎn)激勵的鉆柱激振減阻工具的研制[J].石油機械,2015,43(5):9-12.

    ZHANG H Z,GUAN Z C,LIU Y W,et al.Development of drill string excitation drag reduction tool based on rotary excitation[J].China Petroleum Machinery,2015,43(5):9-12.

    [18] 王甲昌,滕春鳴,張海平,等.機械式旋轉(zhuǎn)沖擊鉆井工具研制及試驗[J].鉆采工藝,2020,43(6):68-71.

    WANG J C,TENG C M,ZHANG H P,et al.Development and experiment of rotary percussion positive displacement motor[J].Drilling & Production Technology,2020,43(6):68-71.

    [19] 席傳明,穆總結(jié),羅翼,等.GCY-I型沖擊螺桿鉆井提速技術(shù)研究與試驗[J].石油機械,2020,48(10):39-43,97.

    XI C M,MU Z J,LUO Y,et al.Research and field test of GCY-I screw impact drilling tools[J].China Petroleum Machinery,2020,48(10):39-43,97.

    [20] 張學(xué)鴻,董振剛,張雄,等.探井水平井整體鉆柱摩阻力分析的間隙元法[J].石油學(xué)報,2002,23(5):105-109.

    ZHANG X H,DONG Z G,ZHANG X,et al.Gap element method for frictional resistance analysis of whole drill string in prospecting horizontal well[J].Acta Petrolei Sinica,2002,23(5):105-109.

    [21] 秦永和,付勝利,高德利.大位移井摩阻扭矩力學(xué)分析新模型[J].天然氣工業(yè),2006,26(11):77-79.

    QIN Y H,F(xiàn)U S L,GAO D L.A new model for analyzing torque and drag in extended reach wells[J].Natural Gas Industry,2006,26(11):77-79.

    [22] 閆鐵,李慶明,王巖,等.水平井鉆柱摩阻扭矩分段計算模型[J].大慶石油學(xué)院學(xué)報,2011,35(5):69-72,83.

    YAN T,LI Q M,WANG Y,et al.Segmental calculation model for torque and drag of drillstring in horizontal wells[J].Journal of Daqing Petroleum Institute,2011,35(5):69-72,83.

    [23] 李文哲,代鋒,曾光,等.長寧超長水平井鉆井摩阻數(shù)值模擬分析[J].鉆采工藝,2021,44(5):39-44.

    LI W Z,DAI F,ZENG G,et al.Numerical simulation analysis of drilling friction in Changning ultralong horizontal well[J].Drilling & Production Technology,2021,44(5):39-44.

    第一閆炎,工程師,生于1993年,2021年畢業(yè)于中國石油大學(xué)(華東)油氣井工程專業(yè),獲博士學(xué)位,現(xiàn)從事油氣井巖石力學(xué)研究工作。地址:(030600)陜西省西安市。Email:yanyan3@cnpc.com.cn2022-10-16宋治國

    猜你喜歡
    摩阻水平井
    低滲透油田壓裂水平井生產(chǎn)動態(tài)分析
    云南化工(2020年11期)2021-01-14 00:50:42
    市政橋梁預(yù)應(yīng)力管道摩阻系數(shù)測試研究
    江西建材(2018年4期)2018-04-10 12:37:20
    提高胍膠壓裂液摩阻計算精度的方法
    液體管道瞬變流摩阻模型綜述
    基于水平井信息的單一河口壩內(nèi)部增生體識別
    一種計算水平井產(chǎn)能的新方法
    熱采水平井加熱半徑計算新模型
    致密油藏多級壓裂水平井?dāng)?shù)值模擬及應(yīng)用
    計算隱式摩阻系數(shù)方程數(shù)值解的簡便方法
    考慮扶正器影響的套管摩阻計算方法研究
    亚洲人成77777在线视频| 各种免费的搞黄视频| 日韩人妻精品一区2区三区| 亚洲第一av免费看| 国产亚洲av高清不卡| 国产99久久九九免费精品| 免费在线观看黄色视频的| 午夜福利网站1000一区二区三区| 国产在线免费精品| 国产成人午夜福利电影在线观看| 精品人妻一区二区三区麻豆| 日韩中文字幕欧美一区二区 | 18在线观看网站| 日韩,欧美,国产一区二区三区| 国产精品一二三区在线看| 亚洲欧美中文字幕日韩二区| 精品第一国产精品| 国产精品亚洲av一区麻豆 | 国产黄色视频一区二区在线观看| 啦啦啦视频在线资源免费观看| 国产一区二区在线观看av| 高清黄色对白视频在线免费看| 国产又色又爽无遮挡免| 日韩制服丝袜自拍偷拍| 成年人午夜在线观看视频| 激情五月婷婷亚洲| 99久久精品国产亚洲精品| 亚洲精品视频女| 午夜福利在线免费观看网站| 国产一卡二卡三卡精品 | h视频一区二区三区| 亚洲成人一二三区av| 午夜激情av网站| 久久97久久精品| 久久久精品94久久精品| 99热国产这里只有精品6| 亚洲成人av在线免费| 午夜影院在线不卡| 好男人视频免费观看在线| 欧美中文综合在线视频| 一区二区av电影网| 超色免费av| 51午夜福利影视在线观看| 99re6热这里在线精品视频| 国产日韩欧美亚洲二区| 妹子高潮喷水视频| 国产免费现黄频在线看| 伊人久久大香线蕉亚洲五| 亚洲精品在线美女| 欧美日韩av久久| 亚洲综合精品二区| 五月天丁香电影| 国精品久久久久久国模美| 国产麻豆69| 欧美激情极品国产一区二区三区| 80岁老熟妇乱子伦牲交| 国产精品香港三级国产av潘金莲 | 亚洲欧美日韩另类电影网站| 一本大道久久a久久精品| 国产免费一区二区三区四区乱码| 一区二区三区四区激情视频| 免费久久久久久久精品成人欧美视频| 少妇的丰满在线观看| 菩萨蛮人人尽说江南好唐韦庄| 欧美日韩视频精品一区| 香蕉国产在线看| 欧美精品一区二区免费开放| 只有这里有精品99| 欧美精品人与动牲交sv欧美| 制服丝袜香蕉在线| 欧美精品一区二区大全| 纵有疾风起免费观看全集完整版| 国产成人精品无人区| 夫妻性生交免费视频一级片| 精品一品国产午夜福利视频| 精品亚洲成国产av| 99热全是精品| 成年人午夜在线观看视频| 午夜福利影视在线免费观看| 中文字幕高清在线视频| 亚洲 欧美一区二区三区| 久久国产亚洲av麻豆专区| 久久精品久久精品一区二区三区| 亚洲国产精品一区三区| av一本久久久久| 亚洲欧洲日产国产| 国产男女内射视频| 国产视频首页在线观看| 亚洲精品国产一区二区精华液| 亚洲成国产人片在线观看| av片东京热男人的天堂| 人人妻人人澡人人看| 99香蕉大伊视频| 少妇被粗大猛烈的视频| 欧美日韩一级在线毛片| 91aial.com中文字幕在线观看| 高清黄色对白视频在线免费看| 最近最新中文字幕大全免费视频 | 亚洲,欧美精品.| 另类亚洲欧美激情| videos熟女内射| 女人被躁到高潮嗷嗷叫费观| 国产 精品1| 下体分泌物呈黄色| 亚洲一区中文字幕在线| 老司机深夜福利视频在线观看 | 精品国产乱码久久久久久男人| 欧美日韩亚洲综合一区二区三区_| 99精国产麻豆久久婷婷| 最近中文字幕2019免费版| 亚洲av电影在线进入| 亚洲中文av在线| 亚洲成人av在线免费| 高清黄色对白视频在线免费看| 纵有疾风起免费观看全集完整版| 国产黄色免费在线视频| 这个男人来自地球电影免费观看 | 国产精品嫩草影院av在线观看| 热99国产精品久久久久久7| 欧美精品亚洲一区二区| 下体分泌物呈黄色| 国产精品国产三级国产专区5o| 天堂俺去俺来也www色官网| 老司机靠b影院| 精品人妻一区二区三区麻豆| 人人妻,人人澡人人爽秒播 | 大片免费播放器 马上看| 亚洲综合精品二区| 欧美精品亚洲一区二区| 曰老女人黄片| 国产 精品1| 99热全是精品| 三上悠亚av全集在线观看| 婷婷成人精品国产| 黄网站色视频无遮挡免费观看| 中文乱码字字幕精品一区二区三区| 亚洲欧美一区二区三区国产| 麻豆精品久久久久久蜜桃| 亚洲成人一二三区av| 欧美精品一区二区大全| 国产高清不卡午夜福利| 18在线观看网站| 看非洲黑人一级黄片| 免费在线观看黄色视频的| 久久久久久久精品精品| 亚洲,一卡二卡三卡| 精品一区二区免费观看| 日韩制服骚丝袜av| 亚洲av电影在线进入| 天天影视国产精品| 久久久久国产精品人妻一区二区| 美国免费a级毛片| 欧美国产精品va在线观看不卡| 免费黄色在线免费观看| 99热全是精品| 精品少妇一区二区三区视频日本电影 | 最近最新中文字幕大全免费视频 | 精品国产乱码久久久久久男人| 波野结衣二区三区在线| 90打野战视频偷拍视频| 交换朋友夫妻互换小说| 午夜福利网站1000一区二区三区| 超色免费av| 制服诱惑二区| 欧美日韩视频精品一区| 视频区图区小说| 中国三级夫妇交换| 男女高潮啪啪啪动态图| 人人澡人人妻人| 啦啦啦在线免费观看视频4| 国产精品久久久久成人av| 老司机影院成人| 97人妻天天添夜夜摸| 亚洲精品在线美女| av在线app专区| 国产av一区二区精品久久| 午夜福利免费观看在线| 免费观看av网站的网址| 午夜日本视频在线| 久久精品久久久久久噜噜老黄| 七月丁香在线播放| 夫妻午夜视频| 精品卡一卡二卡四卡免费| 国产日韩欧美视频二区| 18禁动态无遮挡网站| 亚洲成人免费av在线播放| 在线观看免费视频网站a站| 尾随美女入室| 国产高清不卡午夜福利| 国产精品av久久久久免费| 电影成人av| 街头女战士在线观看网站| 免费高清在线观看视频在线观看| 精品一区二区三区四区五区乱码 | 人人妻人人澡人人爽人人夜夜| 亚洲av国产av综合av卡| 天天躁夜夜躁狠狠久久av| 亚洲人成电影观看| 国产成人91sexporn| 精品少妇久久久久久888优播| 2021少妇久久久久久久久久久| 不卡av一区二区三区| 在线观看www视频免费| 一区二区日韩欧美中文字幕| 精品人妻熟女毛片av久久网站| 美女脱内裤让男人舔精品视频| 香蕉国产在线看| 最近最新中文字幕大全免费视频 | 女人久久www免费人成看片| 久久精品人人爽人人爽视色| 男的添女的下面高潮视频| 国产精品一国产av| 国产一区亚洲一区在线观看| 国产日韩一区二区三区精品不卡| 午夜免费鲁丝| 久久精品亚洲熟妇少妇任你| 日韩中文字幕欧美一区二区 | 国产乱人偷精品视频| 免费看不卡的av| 亚洲成人av在线免费| 亚洲美女黄色视频免费看| 2018国产大陆天天弄谢| 黄片播放在线免费| 国精品久久久久久国模美| a 毛片基地| 99久久人妻综合| 精品人妻在线不人妻| 色网站视频免费| 国产熟女午夜一区二区三区| 国产精品三级大全| 2021少妇久久久久久久久久久| 纯流量卡能插随身wifi吗| 免费看不卡的av| 国产福利在线免费观看视频| 日韩精品有码人妻一区| 夫妻午夜视频| 国产精品久久久久成人av| 欧美日韩亚洲高清精品| 日韩精品免费视频一区二区三区| 国产在线免费精品| 男人爽女人下面视频在线观看| 国产精品国产三级国产专区5o| 涩涩av久久男人的天堂| 亚洲男人天堂网一区| 黄频高清免费视频| 极品人妻少妇av视频| 我要看黄色一级片免费的| 狠狠精品人妻久久久久久综合| 午夜影院在线不卡| 宅男免费午夜| 国产福利在线免费观看视频| 又大又黄又爽视频免费| 男女高潮啪啪啪动态图| 国产 精品1| 一个人免费看片子| 色94色欧美一区二区| 日韩伦理黄色片| 欧美另类一区| 不卡av一区二区三区| 亚洲图色成人| 日本黄色日本黄色录像| 一级,二级,三级黄色视频| 亚洲av福利一区| 日日爽夜夜爽网站| 黄网站色视频无遮挡免费观看| 国产日韩欧美在线精品| 国产成人精品在线电影| 看非洲黑人一级黄片| 午夜免费鲁丝| 在线免费观看不下载黄p国产| 国产免费又黄又爽又色| 国产男女超爽视频在线观看| 日韩一卡2卡3卡4卡2021年| 啦啦啦视频在线资源免费观看| 亚洲精品成人av观看孕妇| 日本欧美视频一区| 精品酒店卫生间| 欧美日韩视频精品一区| 国产高清国产精品国产三级| 欧美 亚洲 国产 日韩一| 乱人伦中国视频| 精品视频人人做人人爽| 免费高清在线观看日韩| 午夜av观看不卡| 精品第一国产精品| 亚洲国产精品成人久久小说| 丝袜喷水一区| 90打野战视频偷拍视频| 亚洲精品,欧美精品| av天堂久久9| 热99国产精品久久久久久7| 婷婷色av中文字幕| 久久婷婷青草| 久久人妻熟女aⅴ| 国产成人免费观看mmmm| 国产精品一区二区在线不卡| 色婷婷久久久亚洲欧美| av在线老鸭窝| 性高湖久久久久久久久免费观看| 国产极品天堂在线| 欧美日韩国产mv在线观看视频| 国产无遮挡羞羞视频在线观看| 精品免费久久久久久久清纯 | 97在线人人人人妻| 男女下面插进去视频免费观看| 在线观看国产h片| www.自偷自拍.com| 亚洲一区二区三区欧美精品| 久久久国产一区二区| 日韩中文字幕欧美一区二区 | 亚洲av日韩在线播放| 黄片播放在线免费| 热re99久久精品国产66热6| 黄片小视频在线播放| 日韩一卡2卡3卡4卡2021年| 天堂8中文在线网| 免费黄频网站在线观看国产| 91国产中文字幕| 成人18禁高潮啪啪吃奶动态图| 建设人人有责人人尽责人人享有的| 亚洲av电影在线进入| 午夜久久久在线观看| 亚洲国产欧美日韩在线播放| 国产成人精品无人区| 亚洲成国产人片在线观看| 久久久久国产精品人妻一区二区| 精品人妻熟女毛片av久久网站| 亚洲一区二区三区欧美精品| 欧美中文综合在线视频| 亚洲av综合色区一区| 巨乳人妻的诱惑在线观看| 日韩人妻精品一区2区三区| 日本av手机在线免费观看| 99久久人妻综合| 亚洲人成77777在线视频| 久久久久国产一级毛片高清牌| 不卡av一区二区三区| 国产又爽黄色视频| 国产熟女欧美一区二区| 久久精品久久久久久噜噜老黄| 日韩成人av中文字幕在线观看| 亚洲熟女毛片儿| a级片在线免费高清观看视频| 久久久亚洲精品成人影院| 自拍欧美九色日韩亚洲蝌蚪91| 日韩免费高清中文字幕av| 九九爱精品视频在线观看| 男女下面插进去视频免费观看| 日本猛色少妇xxxxx猛交久久| 9热在线视频观看99| 日韩精品免费视频一区二区三区| 国产福利在线免费观看视频| 中文字幕色久视频| 免费在线观看视频国产中文字幕亚洲 | 两个人免费观看高清视频| 男的添女的下面高潮视频| 亚洲国产精品一区三区| 美女视频免费永久观看网站| 成人亚洲精品一区在线观看| 亚洲国产中文字幕在线视频| 欧美日韩av久久| 久久婷婷青草| 免费高清在线观看视频在线观看| 久久久久久久久免费视频了| 在线观看免费视频网站a站| 欧美激情高清一区二区三区 | 日韩中文字幕欧美一区二区 | kizo精华| 日日撸夜夜添| 亚洲自偷自拍图片 自拍| 宅男免费午夜| 精品一区二区免费观看| 国产精品国产三级国产专区5o| 色网站视频免费| 人人澡人人妻人| 亚洲成人免费av在线播放| 婷婷色综合大香蕉| 国产亚洲av高清不卡| 天天添夜夜摸| www.熟女人妻精品国产| 不卡av一区二区三区| 免费少妇av软件| 高清黄色对白视频在线免费看| 亚洲国产中文字幕在线视频| 久久久久国产精品人妻一区二区| 亚洲欧美色中文字幕在线| 欧美老熟妇乱子伦牲交| 丝袜人妻中文字幕| 日本av手机在线免费观看| 午夜精品国产一区二区电影| 女人高潮潮喷娇喘18禁视频| 精品人妻熟女毛片av久久网站| 交换朋友夫妻互换小说| 男人操女人黄网站| 国产av国产精品国产| 国产欧美日韩一区二区三区在线| 国产欧美亚洲国产| 国产精品 欧美亚洲| 国产乱来视频区| 亚洲视频免费观看视频| 观看av在线不卡| 午夜日韩欧美国产| 无遮挡黄片免费观看| 纵有疾风起免费观看全集完整版| 亚洲欧洲日产国产| 一边摸一边做爽爽视频免费| 深夜精品福利| 国产精品人妻久久久影院| 如何舔出高潮| 1024视频免费在线观看| 热re99久久国产66热| 国产亚洲av片在线观看秒播厂| 亚洲精品自拍成人| 久久97久久精品| 新久久久久国产一级毛片| 中文字幕人妻丝袜制服| 大片电影免费在线观看免费| 国产精品国产三级国产专区5o| 性高湖久久久久久久久免费观看| 中国三级夫妇交换| 亚洲国产精品成人久久小说| 国产午夜精品一二区理论片| 日本午夜av视频| 国产在视频线精品| 成人午夜精彩视频在线观看| 在线观看免费日韩欧美大片| 亚洲国产欧美网| 欧美黑人精品巨大| 成人黄色视频免费在线看| 我要看黄色一级片免费的| 欧美激情 高清一区二区三区| 一区二区av电影网| 国产黄色视频一区二区在线观看| 九草在线视频观看| 美女扒开内裤让男人捅视频| 热re99久久精品国产66热6| 亚洲人成77777在线视频| 国产片内射在线| 18禁动态无遮挡网站| 女性生殖器流出的白浆| 男女边吃奶边做爰视频| 在线亚洲精品国产二区图片欧美| 精品国产一区二区三区久久久樱花| 久久久久久免费高清国产稀缺| 美女脱内裤让男人舔精品视频| 久久久久久久久免费视频了| 免费黄网站久久成人精品| 波野结衣二区三区在线| 国产人伦9x9x在线观看| videos熟女内射| 香蕉国产在线看| 侵犯人妻中文字幕一二三四区| 美女高潮到喷水免费观看| 色综合欧美亚洲国产小说| 蜜桃国产av成人99| 中文字幕人妻丝袜制服| 青青草视频在线视频观看| 人人妻,人人澡人人爽秒播 | 成年人午夜在线观看视频| 人人妻人人爽人人添夜夜欢视频| 国产麻豆69| 久热爱精品视频在线9| 欧美日本中文国产一区发布| 十分钟在线观看高清视频www| 最近2019中文字幕mv第一页| 国产欧美亚洲国产| 男女高潮啪啪啪动态图| 晚上一个人看的免费电影| 亚洲欧美一区二区三区黑人| 亚洲av电影在线观看一区二区三区| 中文字幕制服av| 男女免费视频国产| 黑人欧美特级aaaaaa片| 男女午夜视频在线观看| 国产精品 国内视频| av在线播放精品| 欧美97在线视频| 一区二区av电影网| 水蜜桃什么品种好| 国产精品.久久久| 精品一区在线观看国产| 熟妇人妻不卡中文字幕| 色网站视频免费| 久久ye,这里只有精品| 国产在线一区二区三区精| 亚洲激情五月婷婷啪啪| 午夜日韩欧美国产| 18禁国产床啪视频网站| 日本猛色少妇xxxxx猛交久久| 久久精品国产a三级三级三级| av免费观看日本| 国产精品99久久99久久久不卡 | 男女高潮啪啪啪动态图| www.自偷自拍.com| 日韩大码丰满熟妇| 亚洲综合精品二区| 一本一本久久a久久精品综合妖精| 日韩制服骚丝袜av| 亚洲精品成人av观看孕妇| 夫妻性生交免费视频一级片| av天堂久久9| tube8黄色片| 中文字幕人妻熟女乱码| 九色亚洲精品在线播放| 成人亚洲欧美一区二区av| 亚洲情色 制服丝袜| 欧美av亚洲av综合av国产av | 日韩 亚洲 欧美在线| 黑丝袜美女国产一区| 色视频在线一区二区三区| 亚洲精品久久成人aⅴ小说| 日本wwww免费看| 精品国产超薄肉色丝袜足j| 免费观看av网站的网址| 亚洲 欧美一区二区三区| 人妻 亚洲 视频| 少妇人妻久久综合中文| 久久精品久久久久久久性| 亚洲七黄色美女视频| 在线观看免费高清a一片| 1024香蕉在线观看| 中国三级夫妇交换| 午夜精品国产一区二区电影| 电影成人av| 亚洲精品国产av成人精品| 国产成人一区二区在线| 日韩av免费高清视频| av线在线观看网站| 两性夫妻黄色片| 亚洲欧美中文字幕日韩二区| 久久女婷五月综合色啪小说| 久久青草综合色| 久久天堂一区二区三区四区| xxx大片免费视频| 高清在线视频一区二区三区| 亚洲熟女精品中文字幕| 久热爱精品视频在线9| 欧美xxⅹ黑人| 99久久综合免费| 国产一区二区三区综合在线观看| 精品一区在线观看国产| 精品一品国产午夜福利视频| 少妇人妻久久综合中文| 在线亚洲精品国产二区图片欧美| 欧美亚洲 丝袜 人妻 在线| 曰老女人黄片| 天堂中文最新版在线下载| 亚洲婷婷狠狠爱综合网| 在线天堂中文资源库| 国产日韩一区二区三区精品不卡| 日韩一区二区视频免费看| 成年女人毛片免费观看观看9 | 美女高潮到喷水免费观看| 黄片播放在线免费| 国产午夜精品一二区理论片| 大片电影免费在线观看免费| 无遮挡黄片免费观看| 亚洲国产av影院在线观看| 国产欧美日韩一区二区三区在线| 国产精品偷伦视频观看了| 久久人人爽人人片av| 国产精品久久久久久人妻精品电影 | 美女大奶头黄色视频| 日本欧美国产在线视频| 国产精品av久久久久免费| 毛片一级片免费看久久久久| 亚洲男人天堂网一区| 精品一区在线观看国产| h视频一区二区三区| av又黄又爽大尺度在线免费看| 无遮挡黄片免费观看| 大码成人一级视频| 精品视频人人做人人爽| 中文字幕最新亚洲高清| kizo精华| 久久久久精品国产欧美久久久 | 岛国毛片在线播放| 9色porny在线观看| 国产探花极品一区二区| 日日爽夜夜爽网站| 国产免费现黄频在线看| av卡一久久| 国产成人精品久久久久久| 日韩 亚洲 欧美在线| 十八禁人妻一区二区| 在线免费观看不下载黄p国产| av.在线天堂| 亚洲欧美一区二区三区久久| 啦啦啦 在线观看视频| 精品第一国产精品| av视频免费观看在线观看| 操出白浆在线播放| 九色亚洲精品在线播放| 精品酒店卫生间| 国产男女内射视频| 中文字幕制服av| 纵有疾风起免费观看全集完整版| 亚洲欧洲日产国产| 女人被躁到高潮嗷嗷叫费观| 精品国产乱码久久久久久男人| 色视频在线一区二区三区| 日本欧美视频一区| 日韩不卡一区二区三区视频在线| 最近2019中文字幕mv第一页| 久久女婷五月综合色啪小说| 看免费av毛片| 精品人妻熟女毛片av久久网站| 国产精品一区二区在线不卡| videos熟女内射| 国产精品一国产av| 九色亚洲精品在线播放| 久久99一区二区三区| 亚洲自偷自拍图片 自拍| 啦啦啦 在线观看视频| av网站在线播放免费|