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

    連續(xù)管作業(yè)氣液兩相流多模型組合計算方法

    2024-12-31 00:00:00朱碩劉閱祝葉黃裕航鄭磊栗文銀
    石油機(jī)械 2024年11期
    關(guān)鍵詞:優(yōu)化模型

    目前連續(xù)管作業(yè)井筒氣液兩相流壓力計算以單一模型為主,該模型難以有效模擬整個井筒的實(shí)際流動狀況,亟需解決單一氣液兩相流動模型僅適用于一種流動工況,與實(shí)際井筒同時存在多種流動工況兩者難以匹配的問題。為此,通過試驗數(shù)據(jù)建立不同流動工況下的最優(yōu)流動模型,在此基礎(chǔ)上形成多模型組合優(yōu)化的氣液兩相流計算方法。該方法能夠針對不同的流動工況將整個井筒(管內(nèi)/環(huán)空)合理劃分為多個井段,每個井段采用對應(yīng)流動工況下的最佳氣液兩相模型,最終完成整個井筒的壓耗計算。通過將環(huán)空壓力與泵壓的相對誤差分別與其他5種單一氣液兩相模型相比可知,優(yōu)化模型綜合平均相對誤差最低為8.23%。該方法可為連續(xù)管作業(yè)施工方案的合理設(shè)計提供理論參考。

    連續(xù)管作業(yè);井筒氣液兩相流;組合優(yōu)化模型;流動工況;環(huán)空壓力;泵壓

    中圖分類號:TE319

    文獻(xiàn)標(biāo)識碼:A

    DOI:10.16082/i.cnki.issn.1001-4578.2024.11.010

    基金項目:中國石油天然氣集團(tuán)有限公司科學(xué)研究與技術(shù)開發(fā)項目“200 ℃/105 MPa連續(xù)管鉆磨射流作業(yè)工具研制”(2020b-4018);中國石油工程技術(shù)研究院有限公司產(chǎn)業(yè)基金項目“SmartCT系列連續(xù)管作業(yè)軟件”(CPETCY202411);中石油江漢機(jī)械研究所有限公司研究課題“連續(xù)管技術(shù)力學(xué)及巖屑運(yùn)移理論研究”(JJY202305)。

    Gas-Liquid Two-Phase Flow Combined Model

    for Coiled Tubing Operation

    Zhu Shuo1" Liu Yue1" Zhu Ye1" Huang Yuhang2" Zheng Lei1" Li Wenyin1

    (1.Jianghan Machinery Research Institute Limited Company of CNPC;2.Coiled Tubing Operation Company of CCDC Downhole Service Company)

    A single model is usually used for calculating the gas-liquid two-phase flow pressure in wellbore with coiled tubing (CT) operation.However,this model is only applicable to one flow condition in the wellbore,but cannot effectively simulate the actual multiple flow conditions in the whole wellbore.In this paper,the optimal flow model for each flow condition was built using the test data,and accordingly an optimized gas-liquid two-phase flow combined model was formed.This combined model can divide the whole wellbore (in string/annulus) into multiple sections depending on flow conditions.The optimal gas-liquid two-phase model for corresponding flow condition is used for each section,and ultimately the pressure loss calculation is completed for the whole wellbore.The combined model was compared with each of the five optimal models in terms of relative error between annulus pressure and pump pressure.It is found that the average relative error of the combined models is 8.23%.This proposed combined model provides a theoretical reference for the rational design of CT operation plan.

    coiled tubing operation;gas-liquid two-phase flow in wellbore;combined model;flow condition;annulus pressure;pump pressure

    0" 引" 言

    氣液兩相流模型廣泛應(yīng)用于連續(xù)管氣舉排液、速度管柱、井筒清理、小井眼側(cè)鉆等鉆修井作業(yè)場景中[1-3],在作業(yè)前、作業(yè)過程中需準(zhǔn)確分析作業(yè)流體在井筒中產(chǎn)生的壓耗,計算全井筒壓力分布,避免地層漏失、井筒溢流等復(fù)雜情況發(fā)生。隨著大位移井、長水平井等復(fù)雜結(jié)構(gòu)井型增多,連續(xù)管作業(yè)環(huán)境更加惡劣,氣液兩相壓降計算更加困難。因此,亟需針對復(fù)雜井型中氣液兩相流壓力分布的準(zhǔn)確計算進(jìn)行研究。

    朱碩,等:連續(xù)管作業(yè)氣液兩相流多模型組合計算方法

    目前,氣液兩相流模型自身已經(jīng)發(fā)展得較為成熟,常見的氣液兩相流模型有經(jīng)驗?zāi)P团c機(jī)理模型2大類。經(jīng)驗?zāi)P屯ㄟ^室內(nèi)試驗數(shù)據(jù)擬合相關(guān)經(jīng)驗關(guān)系式,具有效率高、穩(wěn)定性強(qiáng)、精度尚可的特點(diǎn)。機(jī)理模型需要多次迭代求解封閉微分方程組,具有效率低、穩(wěn)定性差,但精度高的特點(diǎn)。在石油工程應(yīng)用中,為同時保證計算效率和精度,通常采用經(jīng)驗?zāi)P陀嬎闳矇毫Ψ植?。對于常用的幾種經(jīng)驗?zāi)P?,根?jù)是否考慮氣液兩相滑移、是否考慮流型可分為3類。以Poettmann amp; Carpenter[4]為代表的是不考慮兩相滑移、不考慮流型的氣液兩相模型,該類模型基于注入氣液兩相的氣液比計算混合流體的密度,即認(rèn)為氣液兩相以相同的速度在井筒內(nèi)運(yùn)移;以Hagedorn amp; Brown[5]為代表的考慮兩相滑移、未考慮流型的氣液兩相模型,該類模型認(rèn)為氣相和液相以不同速度在井筒內(nèi)運(yùn)移,通過試驗數(shù)據(jù)擬合得到持液率計算公式,在所有流動類型下均采用相同的持液率和摩阻系數(shù)計算公式;以Duns amp; Ros[6]、Mukherjee amp; Brill等[7-9]為代表的是考慮兩相滑移、考慮流型的氣液兩相模型,該類模型通過室內(nèi)試驗研究得到不同流型下的持液率和摩阻系數(shù)計算公式,分流型進(jìn)行計算。在機(jī)理模型方面,ZHANG H.Q.等[10-11]基于段塞流流動機(jī)理建立的統(tǒng)一的氣液兩相流動模型,適用于不同傾斜角、不同流動方向的氣液兩相流動計算;H.V.NICKENS[12]直接從多相流偏微分方程組出發(fā),建立鉆井過程中考慮氣體侵入的氣液兩相流計算方法。

    目前,在計算整個井筒氣液兩相流動壓力分布時,均采用其中一類模型進(jìn)行整個井筒的壓力模擬計算。但工作流體在整個井筒中會同時存在向上流、向下流、水平流動等多種工況,而每個模型均在一定的試驗條件內(nèi)所建立,都有各自的適用范圍,單一的模型僅能滿足其中一種流動工況,因此,常常導(dǎo)致計算結(jié)果與現(xiàn)場實(shí)測數(shù)據(jù)存在一定的誤差。

    本文建立的氣液兩相流多模型組合計算方法,首先通過試驗數(shù)據(jù)分析各種流動工況(向上流動、向下流動、水平流動)下的最佳氣液兩相流動模型,準(zhǔn)確計算相應(yīng)工況下的流動壓耗。其次以井斜角、流動方向為劃分依據(jù),建立多模型組合優(yōu)化的氣液兩相流計算方法,該方法能夠針對復(fù)雜井型同時存在的多種流動工況進(jìn)行差異化計算,即將整個井筒依據(jù)不同流動工況合理劃分為若干個井段,每個井段采用最佳的氣液兩相模型,最終完成整個井筒的壓耗計算。

    1" 氣液兩相流模型建立與評估

    選用5種經(jīng)典氣液兩相模型作為組合優(yōu)化模型的待選模型,如表1所示。首先梳理各個模型的適用范圍、關(guān)鍵參數(shù)的求解方法,編程求解5種氣液兩相模型,最后利用一口試驗井?dāng)?shù)據(jù)及相關(guān)商業(yè)軟件計算結(jié)果驗證5種氣液兩相模型的計算精度與可靠性。

    1.1" 氣液兩相模型建立

    氣液兩相流模型的壓降梯度計算公式為:

    dpdL=dpdLf+dpdLe+dpdLa(1)

    式中:dpdL為總壓力梯度,Pa/m;dpdLf為摩擦壓力梯度,Pa/m;dpdLe為靜液柱壓力梯度,Pa/m;dpdLa為加速度壓力梯度,Pa/m。

    式(1)可進(jìn)一步變化為:

    dpdL=-τπdA-ρgsin θ-ρvdvdL(2)

    式中:τ為剪切應(yīng)力,N/m2;d為管道直徑,m;A為管道截面積,m2;ρ為流體密度,kg/m3;g為重力加速度,m/s2;θ為管柱與水平方向的夾角,rad;v為流體流速,m/s。

    表1中的不同氣液兩相模型均遵循式(2),僅僅是在計算持液率、摩擦因數(shù)等關(guān)鍵參數(shù)時各模型有各自的計算方法,且靜液柱壓力梯度的準(zhǔn)確計算決定了整個井筒壓降計算的誤差。

    靜液柱壓力梯度的計算核心在于持液率/持氣率的準(zhǔn)確計算,不同的模型有不同計算準(zhǔn)則與方法:

    dpdLe=ρmg(3)

    ρm=ρLHL+ρGHG(4)

    式中:ρm為氣液兩相混合密度,kg/m3;ρL為液體的密度,kg/m3;ρG為氣體的密度,kg/m3;HL為持液率,%;HG為持氣率,%。

    摩擦壓力梯度采用如下公式計算:

    dpdLf=fρv22d(5)

    式中:f為摩擦因數(shù),無量綱;d為管道直徑,m。

    關(guān)于摩擦因數(shù)f、流體流速v,不同的模型有不同的計算方法。

    1.2" 氣液兩相模型評估

    為驗證1.1中氣液兩相模型的可靠性與穩(wěn)定性,利用大尺寸氣液兩相流試驗數(shù)據(jù)[15]以及商業(yè)軟件Cerburs計算結(jié)果進(jìn)行對比驗證,為后續(xù)組合優(yōu)化模型的建立提供可靠的模型基礎(chǔ)。

    1.2.1" 試驗數(shù)據(jù)

    試驗井井身如圖1所示。該試驗井為一口垂直井,井深1 278 m,套管內(nèi)徑159.4 mm,油管外徑88.9 mm、內(nèi)徑70.2 mm。通過油管注入氮?dú)夂退?,分離器處的壓力為0.515 MPa,詳細(xì)試驗參數(shù)如表2所示。該試驗井更貼近實(shí)際油氣井的尺寸,具有可靠的對比驗證價值。

    1.2.2" 模型對比評估

    利用1.1節(jié)所建立的5種模型以及商業(yè)軟件Cerberus中相應(yīng)算法模擬試驗井中的10種工況。由于已知井口分離器處的壓力為0.515MPa,所以應(yīng)以已知井口壓力為起始計算點(diǎn),分段迭代計算得到井底壓力。

    5種氣液兩相模型及Cerberus商業(yè)軟件詳細(xì)計算結(jié)果如表3所示。表3中DR為Duns amp; Ros模型,HB為Hagedorn amp; Brown模型,BB為Beggs amp; Brill模型,MB為Mukherjee amp; Brill模型。

    各個模型、Cerberus模擬結(jié)果與試驗數(shù)據(jù)的對比如圖2~圖6所示。

    從圖中2可以看出:相較Cerberus計算結(jié)果,DR計算結(jié)果變化趨勢及數(shù)值大小基本符合;與試驗數(shù)據(jù)相比,其與試驗數(shù)據(jù)數(shù)值大小均基本一致。相較試驗值,DR計算結(jié)果平均相對誤差為5.1%,Cerberus計算結(jié)果平均相對誤差為4.2%。基于上述分析,可充分證明DR算法的可靠性。

    從圖3可以看出:相較Cerberus計算結(jié)果,HB計算結(jié)果變化趨勢基本一致,部分?jǐn)?shù)值存在較小差距;與試驗數(shù)據(jù)相比,2種方法計算結(jié)果均小于試驗值?;谏鲜龇治?,可保證所建立的HB算法基本可靠,但在試驗井條件下,該算法計算精度稍差。

    從圖4可以看出,BB與Cerberus計算結(jié)果在變化趨勢上與試驗數(shù)據(jù)基本一致,但數(shù)值上均差異明顯。BB模型計算結(jié)果較試驗值偏小,而Cerberus計算結(jié)果相較試驗值均偏大(該現(xiàn)象與Payne的結(jié)論一致,即原始的BB模型計算的持液率偏大,因此導(dǎo)致試驗井工況下井底壓力計算偏大)。本文的BB模型采用了Payne給出的修正系數(shù),因此相同條件下,計算結(jié)果偏小。相較試驗值,本文的BB計算結(jié)果平均相對誤差為12.1%,而Cerberus計算結(jié)果平均相對誤差為20.2%?;谏鲜龇治觯疚牡腂B算法基本可靠。

    從圖5可以看出:Gray與Cerberus計算結(jié)果相比較,其變化趨勢基本一致,部分?jǐn)?shù)值存在較小差距;與試驗數(shù)據(jù)相比,2種方法計算結(jié)果均小于試驗值?;谏鲜龇治?,可保證本文所建立的Gray算法基本可靠,但在試驗井條件下,該算法計算精度較差。

    從圖6可知,MB計算結(jié)果與試驗數(shù)據(jù)數(shù)值變化趨勢基本一致。相較試驗值,平均相對誤差為13.5%,可保證本文所建立的MB算法基本可靠。

    綜合上述分析,所建立的5種氣液兩相模型基本能夠滿足工程計算需求。

    2" 氣液兩相流模型適應(yīng)性分析

    在連續(xù)管作業(yè)過程中,流體在井筒內(nèi)的流動同時存在不同流動工況,包括向下流動、向上流動、水平流動。因此,應(yīng)明確在不同流動工況下采用何種模型計算最為準(zhǔn)確。從表1可知,Duns amp; Ros(DR)、Hagedorn amp; Brown(HB)、Gray均為在向上流動工況下所建立的模型,而Beggs amp; Brill(BB)、Mukherjee amp; Brill(MB)為在全工況下所建立的模型,而MB是對BB模型的進(jìn)一步優(yōu)化與完善。為減少模型對比優(yōu)選工作量,選取在向上流動工況中表現(xiàn)較好的DR模型與全流動工況的MB模型進(jìn)行對比分析。

    2.1" 向上流動

    向上流動工況下的模型對比仍采用1.2.2節(jié)中試驗井試驗數(shù)據(jù),該井在環(huán)空部分為向上流動,并沿用其模型對比結(jié)果,分別如圖7、圖8所示。

    從圖7可見,5種模型計算結(jié)果變化趨勢均與試驗值基本一致,其中,DR、BB、MB模型計算結(jié)果與試驗值符合較好。

    從圖8可見,DR模型的計算結(jié)果與試驗值的相對誤差均遠(yuǎn)小于其他模型的相對誤差,且平均相對誤差接近5%。因此,通過與試驗數(shù)據(jù)的對比可確定,對于向上流動的工況,應(yīng)該選擇DR模型計算。

    2.2" 向下流動

    對于向下流動與水平流動工況,選擇了Stanford[16]氣液兩相流數(shù)據(jù)集中的部分?jǐn)?shù)據(jù)。該試驗數(shù)據(jù)由室內(nèi)試驗獲得,傾斜角度為-90°~90°,相較Lage試驗數(shù)據(jù),管道試驗段尺寸小、長度短,試驗結(jié)果包括持液率、壓降、流動類型。這里分別對比向下流動工況下DR、MB的計算結(jié)果,對比結(jié)果如圖9、圖10所示。

    從圖9可知,在壓降計算方面,雖然在個別數(shù)據(jù)點(diǎn)上,DR與MB模型計算結(jié)果與試驗值存在較大誤差(誤差可能來自于試驗數(shù)據(jù)記錄不準(zhǔn),或試驗測量誤差,或模型不能準(zhǔn)確計算所有參數(shù)組合下的壓降),大部分?jǐn)?shù)據(jù)點(diǎn)均與試驗值符合較好。其中,DR與試驗值的平均相對誤差為3.43%,MB的平均相對誤差為1.21%。

    從圖10可知,在持液率計算方面,整體上2個模型計算值均與試驗值符合較好,且2個模型的平均相對誤差基本一致,DR平均相對誤差為0.43%,MB平均相對誤差為0.47%。

    綜上所述,在向下流動工況條件下,MB模型表現(xiàn)優(yōu)于DR模型,因此,在向下流動工況下應(yīng)盡量選取MB模型進(jìn)行模擬計算。

    2.3" 水平流動

    水平流動下DR與MB持液率計算結(jié)果與試驗值對比如圖11、表4所示。從圖11和表4可知,雖然在某些數(shù)據(jù)點(diǎn)上,DR模型的計算值更接近試驗值,但整體上MB模型持液率計算值更穩(wěn)定,也更貼近試驗值,且該模型本身是在-90°~90°室內(nèi)試驗條件下所建立的。MB能夠準(zhǔn)確預(yù)測更多的流動型態(tài)。因此,在水平流動工況下應(yīng)采用MB模型。

    3" 氣液兩相組合優(yōu)化模型建立與案例分析

    基于對不同流動工況下的氣液兩相流模型適應(yīng)性分析,可建立氣液兩相組合優(yōu)化模型,如表5、圖12所示。

    利用1口連續(xù)管欠平衡鉆井的井底壓力隨鉆測量數(shù)據(jù)對模型進(jìn)行驗證。該井為1口連續(xù)管欠平衡鉆井的水平井,井深1 597 m,油管外徑88.9 mm、內(nèi)徑70.2 mm,油管總長3 246 m,油管下深1 300 m。一開套管下深1 121 m,內(nèi)徑125.73 mm,裸眼段內(nèi)徑120.65 mm。 油管注入氮?dú)夂望}水,井口溫度15 ℃,井底溫度39 ℃。部分隨鉆測量井底壓力數(shù)據(jù)如表6所示。

    分別利用5種氣液兩相模型和優(yōu)化模型計算相同工況下的環(huán)空壓力與泵壓,計算結(jié)果與測量值的對比如表7、表8所示,平均相對誤差如表9所示。優(yōu)化模型的環(huán)空壓力與流速如圖13、圖14所示。

    在環(huán)空壓力計算方面,從表9可見,計算誤差較小的為DR模型與優(yōu)化模型,且這2種方法計算得到的環(huán)空壓力值(見表7)在數(shù)值上一致;從圖13可知,在環(huán)空0~1 050 m區(qū)間內(nèi),井斜角小于80°,垂直向上流動,2種模型均采用DR計算方法,在該區(qū)間內(nèi),因采用同一公式計算,所以壓降計算結(jié)果一致;在1 050~1 300 m區(qū)間內(nèi),該段為水平段,靜液柱壓降為0,優(yōu)化模型采用MB,計算的流型為彈狀流,采用DR模型計算的流型為氣泡流,雖然預(yù)測的流型不一致,但2種方法在相應(yīng)流型下,采用的摩擦壓降公式基本一致,因此在1 050~1 300 m內(nèi),2種模型計算的壓降基本一致。綜上所述,在該流動工況下,2種模型計算得出的環(huán)空壓力基本一致。

    在泵壓計算方面,從表9可見,計算誤差較小的有BB、MB與優(yōu)化模型,雖然優(yōu)化模型并不是表現(xiàn)最好的,但其誤差在10%左右,基本能夠滿足工程上的應(yīng)用需求。

    此外,綜合環(huán)空壓力與泵壓2個關(guān)鍵壓力點(diǎn)的預(yù)測誤差,優(yōu)化模型的平均相對誤差最小,為8.23%,表現(xiàn)最好。由此說明依據(jù)不同流動工況,采用不同氣液兩相模型的優(yōu)化思路能夠提高壓降的計算精度。

    4" 結(jié)" 論

    針對單一氣液兩相流動模型僅適用于一種流動工況,與實(shí)際井筒同時存在多種流動工況兩者難以匹配的問題,提出一種井筒氣液兩相流多模型組合優(yōu)化計算方法,并得出以下結(jié)論。

    (1)前人已建立了多種不同工況下的氣液兩相流經(jīng)驗?zāi)P?,在工程?yīng)用時不應(yīng)盲目采用某種單一模型,應(yīng)具體評估不同流動工況下各種氣液兩相模型的表現(xiàn)。

    (2)以井斜角、流動方向為劃分依據(jù),根據(jù)不同流動工況將整個井筒合理劃分為多個井段,每個井段采用最佳氣液兩相模型,建立多模型組合優(yōu)化的氣液兩相流動壓力梯度計算方法,實(shí)現(xiàn)井筒壓力的準(zhǔn)確計算。

    (3)將該方法用于一口連續(xù)管欠平衡鉆井案例分析中,結(jié)果表明:通過不同流動工況選擇不同氣液兩相流動模型計算流動壓力梯度,本文模型綜合平均相對誤差最低為8.23%,相較采用單一模型計算,組合優(yōu)化模型能夠更好地模擬井筒復(fù)雜流動工況下的流動壓力梯度,有效規(guī)避欠平衡鉆井溢漏風(fēng)險,保證連續(xù)管側(cè)鉆安全、高效進(jìn)行。

    [1]" AYDIN H,MEREY S.Design of electrical submersible pump system in geothermal wells:a case study from West Anatolia,Turkey[J].Energy,2021,230:120891.

    [2]" ESKIN D.An engineering model of a gas bubble-liquid turbulent flow in coiled tubing wound on a reel[J].Chemical Engineering and Processing:Process Intensification,2010,49(8):843-851.

    [3]" 邱小雪,鐘光海,李賢勝,等.不同井斜頁巖氣水平井流動特征的CFD模擬研究[J].油氣藏評價與開發(fā),2023,13(3):340-347.

    QIU X X,ZHONG G H,LI X S,et al.CFD simulation of flow characteristics of shale gas horizontal wells with different inclination[J].Petroleum Reservoir Evaluation and Development,2023,13(3):340-347.

    [4]" POETTMAN F H,CARPENTER P G.The multiphase flow of gas,oil,and water through vertical flow strings with application to the design of gas-lift installations[Z].Drilling and Production Practice,1952:257-317.

    [5]" HAGEDORN A R,BROWN K E.Experimental study of pressure gradients occurring during continuous two-phase flow in small-diameter vertical conduits[J].Journal of Petroleum Technology,1965,17(4):475-484.

    [6]" DUNS H,ROS N C J.Vertical flow of gas and liquid mixtures in wells[C]∥The 6th World Petroleum Congress.Frankfurt am Main,Germany:WPC,1963:WPC-10132.

    [7]" MUKHERJEE H,BRILL J P.Empirical equations to predict flow patterns in two-phase inclined flow[J].International Journal of Multiphase Flow,1985,11(3):299-315.

    [8]" MUKHERJEE H,BRILL J P.Liquid holdup correlations for inclined two-phase flow[J].Journal of Petroleum Technology,1983,35(5):1003-1008.

    [9]" MUKHERJEE H,BRILL J P.Pressure drop correlations for inclined two-phase flow[J].Journal of Energy Resources Technology,1985,107(4):549-554.

    [10]" ZHANG H Q,WANG Q,SARICA C,et al.Unified model for gas-liquid pipe flow via slug dynamics-part 1:model development[J].Journal of Energy Resources Technology,2003,125(4):266-273.

    [11]" ZHANG H Q,WANG Q,SARICA C,et al.Unified model for gas-liquid pipe flow via slug dynamics-part 2:model validation[J].Journal of Energy Resources Technology,2003,125(4):274-283.

    [12]" NICKENS H V.A dynamic computer model of a kicking well[J].SPE Drilling Engineering,1987,2(2):159-173.

    [13]" OYEWOLE A.Extension of the gray correlation to inclination angles[C]∥SPE Annual Technical Conference and Exhibition.Houston:SPE,2015:SPE 178727-STU.

    [14]" BEGGS D H,BRILL J P.A Study of two-phase flow in inclined pipes[J].Journal of Petroleum Technology,1973,25(5):607-617.

    [15]" LAGE A C V M,TIME R W.Mechanistic model for upward two-phase flow in annuli[C]∥SPE Annual Technical Conference and Exhibition.Dallas,Texas:SPE,2000:SPE 63127-MS.

    [16]" PETALAS N,AZIZ K.Stanford multiphase flow database-users manual,version 0.2[Z].Petroleum Engineering Department,Stanford University:Stanford,CA,USA,1995.

    第一朱碩,工程師,生于1996年,2022年畢業(yè)于中國石油大學(xué)(北京)油氣井工程專業(yè),現(xiàn)從事連續(xù)管技術(shù)理論與連續(xù)管智能化研究等工作。地址:(430024)湖北省武漢市江漢區(qū)。email:1309656323@qq.com。2024-08-14楊曉峰

    猜你喜歡
    優(yōu)化模型
    一半模型
    超限高層建筑結(jié)構(gòu)設(shè)計與優(yōu)化思考
    民用建筑防煙排煙設(shè)計優(yōu)化探討
    關(guān)于優(yōu)化消防安全告知承諾的一些思考
    一道優(yōu)化題的幾何解法
    由“形”啟“數(shù)”優(yōu)化運(yùn)算——以2021年解析幾何高考題為例
    重要模型『一線三等角』
    重尾非線性自回歸模型自加權(quán)M-估計的漸近分布
    3D打印中的模型分割與打包
    FLUKA幾何模型到CAD幾何模型轉(zhuǎn)換方法初步研究
    亚洲第一青青草原| 国产一区二区三区av在线| 一本—道久久a久久精品蜜桃钙片| 日本av手机在线免费观看| 亚洲欧美日韩另类电影网站| 免费少妇av软件| 90打野战视频偷拍视频| 一区在线观看完整版| 久久久久人妻精品一区果冻| 一级a爱视频在线免费观看| av在线老鸭窝| 满18在线观看网站| 日韩三级伦理在线观看| 日韩免费高清中文字幕av| 亚洲美女搞黄在线观看| 国产在线免费精品| 欧美日韩精品成人综合77777| 国产亚洲最大av| 亚洲成人av在线免费| 制服诱惑二区| 天美传媒精品一区二区| 蜜桃国产av成人99| 在线精品无人区一区二区三| 久久精品国产自在天天线| 综合色丁香网| 午夜精品国产一区二区电影| 看免费成人av毛片| 亚洲av中文av极速乱| 大片电影免费在线观看免费| 亚洲国产色片| 午夜免费男女啪啪视频观看| 国产精品偷伦视频观看了| 亚洲欧洲精品一区二区精品久久久 | 中文欧美无线码| 男女国产视频网站| 秋霞在线观看毛片| 午夜福利一区二区在线看| 午夜福利网站1000一区二区三区| 欧美另类一区| 一级a爱视频在线免费观看| 一区二区日韩欧美中文字幕| 国产极品天堂在线| 青春草亚洲视频在线观看| 少妇人妻精品综合一区二区| 伦精品一区二区三区| 一区二区三区精品91| 搡老乐熟女国产| 亚洲欧美日韩另类电影网站| 国产探花极品一区二区| 69精品国产乱码久久久| 亚洲色图综合在线观看| 纯流量卡能插随身wifi吗| 老汉色av国产亚洲站长工具| 一区二区av电影网| 王馨瑶露胸无遮挡在线观看| 国产日韩欧美在线精品| 国产成人精品无人区| 在线观看三级黄色| 99久国产av精品国产电影| 天天影视国产精品| 亚洲国产精品一区二区三区在线| 欧美精品av麻豆av| 天堂俺去俺来也www色官网| 1024香蕉在线观看| 成人手机av| 美女xxoo啪啪120秒动态图| 亚洲精品成人av观看孕妇| 国产午夜精品一二区理论片| 免费高清在线观看视频在线观看| 免费观看性生交大片5| 如何舔出高潮| 丝袜美腿诱惑在线| 精品午夜福利在线看| www日本在线高清视频| 国产精品久久久av美女十八| 亚洲av中文av极速乱| 亚洲av电影在线观看一区二区三区| 免费黄色在线免费观看| 伦精品一区二区三区| 免费观看a级毛片全部| 精品一区在线观看国产| 热re99久久国产66热| 国产欧美日韩一区二区三区在线| 午夜免费男女啪啪视频观看| 国产精品偷伦视频观看了| 欧美日韩视频高清一区二区三区二| 免费看av在线观看网站| 国产成人精品福利久久| 亚洲综合色惰| 日韩 亚洲 欧美在线| 少妇被粗大的猛进出69影院| 欧美日韩视频精品一区| 少妇猛男粗大的猛烈进出视频| 亚洲美女搞黄在线观看| 成年美女黄网站色视频大全免费| 久久ye,这里只有精品| 丰满饥渴人妻一区二区三| 18+在线观看网站| 国产精品偷伦视频观看了| 看免费av毛片| 久久久久久久亚洲中文字幕| 成人毛片a级毛片在线播放| 亚洲欧洲日产国产| 久久久精品区二区三区| av有码第一页| av网站在线播放免费| 97精品久久久久久久久久精品| a级片在线免费高清观看视频| 两性夫妻黄色片| 高清不卡的av网站| 在线观看人妻少妇| 成人毛片a级毛片在线播放| 国产视频首页在线观看| 亚洲国产欧美日韩在线播放| 永久网站在线| 青青草视频在线视频观看| 久久久久久久大尺度免费视频| videosex国产| 夫妻午夜视频| 成年av动漫网址| 在线天堂最新版资源| 亚洲综合色网址| 97人妻天天添夜夜摸| 宅男免费午夜| 亚洲国产av影院在线观看| 超碰成人久久| 日本黄色日本黄色录像| 亚洲综合色惰| 亚洲三区欧美一区| 国产免费现黄频在线看| 高清av免费在线| 少妇熟女欧美另类| 女人久久www免费人成看片| 国产成人免费观看mmmm| 黄色视频在线播放观看不卡| 成年人免费黄色播放视频| 一区二区av电影网| 曰老女人黄片| 国产精品二区激情视频| 亚洲av成人精品一二三区| 2021少妇久久久久久久久久久| 成人午夜精彩视频在线观看| 亚洲欧洲国产日韩| 美女国产高潮福利片在线看| 亚洲精品美女久久久久99蜜臀 | av在线观看视频网站免费| av国产精品久久久久影院| 亚洲精品一区蜜桃| 一级片免费观看大全| 亚洲情色 制服丝袜| 精品99又大又爽又粗少妇毛片| 国产有黄有色有爽视频| 国产片内射在线| 丝袜在线中文字幕| 免费观看在线日韩| 久久综合国产亚洲精品| 少妇熟女欧美另类| 女性生殖器流出的白浆| 啦啦啦在线观看免费高清www| 青春草国产在线视频| 日本黄色日本黄色录像| 中文字幕色久视频| 天堂俺去俺来也www色官网| 国产免费现黄频在线看| 男女边吃奶边做爰视频| 久久久久久久国产电影| 人成视频在线观看免费观看| 色视频在线一区二区三区| 97在线人人人人妻| 日韩人妻精品一区2区三区| 丰满乱子伦码专区| 韩国精品一区二区三区| 国产精品 国内视频| 久久青草综合色| 91成人精品电影| 国产精品 国内视频| 国产av国产精品国产| 亚洲国产毛片av蜜桃av| 国产精品成人在线| 日韩欧美精品免费久久| 午夜福利一区二区在线看| 看非洲黑人一级黄片| 在现免费观看毛片| 亚洲国产日韩一区二区| 777米奇影视久久| 亚洲精品在线美女| 日本-黄色视频高清免费观看| 成人18禁高潮啪啪吃奶动态图| 国产爽快片一区二区三区| 男女免费视频国产| 免费观看无遮挡的男女| 成年动漫av网址| 日韩中字成人| 少妇人妻精品综合一区二区| 一级毛片黄色毛片免费观看视频| 一区二区三区四区激情视频| 国产av一区二区精品久久| av福利片在线| 在线免费观看不下载黄p国产| 国产一区二区激情短视频 | 中文字幕人妻丝袜一区二区 | 中文字幕精品免费在线观看视频| 看免费成人av毛片| √禁漫天堂资源中文www| 中文字幕亚洲精品专区| 久久免费观看电影| 午夜91福利影院| 少妇的逼水好多| 日韩不卡一区二区三区视频在线| 亚洲三区欧美一区| 国产成人精品久久久久久| 黄色怎么调成土黄色| 亚洲国产毛片av蜜桃av| 99香蕉大伊视频| 女性被躁到高潮视频| 精品人妻熟女毛片av久久网站| 亚洲国产欧美在线一区| 熟妇人妻不卡中文字幕| freevideosex欧美| 国产在视频线精品| 捣出白浆h1v1| 久久久精品94久久精品| 精品亚洲成a人片在线观看| 亚洲精品,欧美精品| 日韩免费高清中文字幕av| 国产成人一区二区在线| 熟女少妇亚洲综合色aaa.| 女人被躁到高潮嗷嗷叫费观| 久久99精品国语久久久| 亚洲精品日本国产第一区| 成年人免费黄色播放视频| 日本猛色少妇xxxxx猛交久久| 成年av动漫网址| 麻豆av在线久日| 精品卡一卡二卡四卡免费| 亚洲人成网站在线观看播放| 亚洲成国产人片在线观看| 亚洲精品中文字幕在线视频| 91午夜精品亚洲一区二区三区| 99久久综合免费| 亚洲视频免费观看视频| 日韩三级伦理在线观看| 熟妇人妻不卡中文字幕| 国产精品成人在线| 欧美日韩精品网址| 日本免费在线观看一区| 在线观看免费高清a一片| 免费看不卡的av| 国产97色在线日韩免费| 另类精品久久| 亚洲精品日韩在线中文字幕| 久久久久精品性色| 人体艺术视频欧美日本| 一区二区三区四区激情视频| 国产毛片在线视频| 欧美另类一区| 国产精品香港三级国产av潘金莲 | 亚洲美女黄色视频免费看| 亚洲av免费高清在线观看| 伦精品一区二区三区| 一级片'在线观看视频| 九九爱精品视频在线观看| 亚洲一区二区三区欧美精品| 丰满少妇做爰视频| 精品人妻熟女毛片av久久网站| 亚洲欧洲日产国产| 日本欧美视频一区| 亚洲av男天堂| 久久久久久人妻| 汤姆久久久久久久影院中文字幕| 97人妻天天添夜夜摸| 香蕉丝袜av| 一区二区三区四区激情视频| 亚洲成av片中文字幕在线观看 | 少妇精品久久久久久久| 国产一区二区在线观看av| 亚洲精品中文字幕在线视频| 2022亚洲国产成人精品| 亚洲欧美一区二区三区国产| 交换朋友夫妻互换小说| 婷婷色综合www| 欧美激情极品国产一区二区三区| 我的亚洲天堂| 欧美激情 高清一区二区三区| 777米奇影视久久| 亚洲伊人色综图| 国产 一区精品| 精品少妇一区二区三区视频日本电影 | 久久久精品区二区三区| 亚洲激情五月婷婷啪啪| a 毛片基地| 在线天堂最新版资源| 国产成人免费无遮挡视频| 1024香蕉在线观看| 极品人妻少妇av视频| 女性被躁到高潮视频| 99re6热这里在线精品视频| 国产精品免费大片| 国产xxxxx性猛交| 国产成人免费观看mmmm| 国产深夜福利视频在线观看| 一区在线观看完整版| 一级黄片播放器| 男女下面插进去视频免费观看| av电影中文网址| 岛国毛片在线播放| 久久热在线av| 亚洲,欧美精品.| 男女边吃奶边做爰视频| 成年女人在线观看亚洲视频| 日韩av免费高清视频| 日本91视频免费播放| 日韩不卡一区二区三区视频在线| 国产爽快片一区二区三区| 日本av免费视频播放| 亚洲五月色婷婷综合| 高清黄色对白视频在线免费看| 咕卡用的链子| 欧美激情高清一区二区三区 | 久热这里只有精品99| 最近最新中文字幕免费大全7| 久久久久久人妻| 哪个播放器可以免费观看大片| 亚洲精品成人av观看孕妇| 亚洲欧美精品自产自拍| 高清欧美精品videossex| 黑丝袜美女国产一区| 亚洲精品av麻豆狂野| 欧美+日韩+精品| 国产高清国产精品国产三级| 日韩在线高清观看一区二区三区| 搡老乐熟女国产| 少妇的逼水好多| 有码 亚洲区| 在线观看一区二区三区激情| xxx大片免费视频| 国产亚洲精品第一综合不卡| 国产97色在线日韩免费| 久久精品国产a三级三级三级| 日韩伦理黄色片| 国产成人欧美| 欧美日韩一级在线毛片| 欧美日韩亚洲高清精品| 最新中文字幕久久久久| 久久久久久人妻| 日韩伦理黄色片| 亚洲国产成人一精品久久久| 免费在线观看黄色视频的| 国产精品蜜桃在线观看| 91国产中文字幕| 亚洲三区欧美一区| 国产精品无大码| 精品国产超薄肉色丝袜足j| 国产极品天堂在线| 久久久久久免费高清国产稀缺| 这个男人来自地球电影免费观看 | 18禁国产床啪视频网站| 久久综合国产亚洲精品| 亚洲精品美女久久av网站| 黄片播放在线免费| 老汉色av国产亚洲站长工具| 国产免费视频播放在线视频| 老司机影院毛片| 99久久综合免费| 亚洲av电影在线进入| 男男h啪啪无遮挡| 国产一区二区激情短视频 | 在线看a的网站| www日本在线高清视频| 久久精品国产亚洲av高清一级| 99久国产av精品国产电影| 黄频高清免费视频| 国产片特级美女逼逼视频| 满18在线观看网站| 中文乱码字字幕精品一区二区三区| 国产福利在线免费观看视频| 国产黄色视频一区二区在线观看| 2022亚洲国产成人精品| 蜜桃国产av成人99| 成人国产麻豆网| 九九爱精品视频在线观看| 一区二区av电影网| 一个人免费看片子| 欧美日韩亚洲国产一区二区在线观看 | 两个人看的免费小视频| 亚洲久久久国产精品| 日本av手机在线免费观看| 欧美日本中文国产一区发布| 午夜日韩欧美国产| 狠狠婷婷综合久久久久久88av| 欧美人与性动交α欧美软件| 国产女主播在线喷水免费视频网站| 亚洲国产色片| 亚洲精品自拍成人| 久久久久久免费高清国产稀缺| 波多野结衣一区麻豆| 亚洲国产成人一精品久久久| 亚洲综合色惰| 精品一区二区三卡| 91久久精品国产一区二区三区| 亚洲欧美成人精品一区二区| 黄色怎么调成土黄色| 久久久久国产网址| 一级毛片我不卡| 狠狠婷婷综合久久久久久88av| 在现免费观看毛片| 亚洲一区二区三区欧美精品| 国产精品久久久久久精品古装| 男人操女人黄网站| 亚洲伊人色综图| 一本—道久久a久久精品蜜桃钙片| 久久久久久久亚洲中文字幕| 黑人猛操日本美女一级片| 18+在线观看网站| 国产日韩一区二区三区精品不卡| 中文字幕人妻丝袜一区二区 | 丝瓜视频免费看黄片| 国产97色在线日韩免费| 午夜福利乱码中文字幕| 黄片播放在线免费| 成人亚洲精品一区在线观看| 亚洲欧美清纯卡通| av网站免费在线观看视频| 少妇 在线观看| 国产黄频视频在线观看| 狠狠精品人妻久久久久久综合| 人人妻人人澡人人看| 精品人妻一区二区三区麻豆| 免费久久久久久久精品成人欧美视频| av视频免费观看在线观看| 丰满迷人的少妇在线观看| 国产在线视频一区二区| 日韩欧美精品免费久久| 国产成人精品久久二区二区91 | 自线自在国产av| 一级毛片 在线播放| 国产xxxxx性猛交| 精品一区在线观看国产| 黄色配什么色好看| 亚洲欧美清纯卡通| 亚洲一区二区三区欧美精品| 国产片特级美女逼逼视频| 国产熟女欧美一区二区| 嫩草影院入口| 老司机影院毛片| 热99国产精品久久久久久7| 久久久亚洲精品成人影院| 啦啦啦啦在线视频资源| 搡女人真爽免费视频火全软件| 精品亚洲成国产av| 国产亚洲午夜精品一区二区久久| 国产精品二区激情视频| 欧美日韩国产mv在线观看视频| 高清黄色对白视频在线免费看| 一级黄片播放器| 欧美日韩精品成人综合77777| 亚洲少妇的诱惑av| 国产精品 欧美亚洲| 亚洲精品国产av蜜桃| 亚洲欧洲国产日韩| 国产日韩欧美在线精品| 国产视频首页在线观看| 免费大片黄手机在线观看| 一本久久精品| 边亲边吃奶的免费视频| 少妇人妻精品综合一区二区| 国产成人精品无人区| 精品久久久精品久久久| 成年人免费黄色播放视频| 国产视频首页在线观看| 成人午夜精彩视频在线观看| 国产男人的电影天堂91| 久久久a久久爽久久v久久| 26uuu在线亚洲综合色| 久久精品久久久久久久性| 日韩人妻精品一区2区三区| www日本在线高清视频| 午夜福利网站1000一区二区三区| 亚洲欧美精品综合一区二区三区 | xxxhd国产人妻xxx| 女性生殖器流出的白浆| 美国免费a级毛片| 成人毛片60女人毛片免费| 狂野欧美激情性bbbbbb| 久久99蜜桃精品久久| 老司机影院成人| 熟女电影av网| 在线观看免费高清a一片| 精品国产一区二区三区久久久樱花| 热99国产精品久久久久久7| a级毛片在线看网站| 青春草国产在线视频| 观看av在线不卡| 国产在线一区二区三区精| 大片电影免费在线观看免费| 亚洲视频免费观看视频| 久久综合国产亚洲精品| 伦理电影大哥的女人| 超碰成人久久| 视频在线观看一区二区三区| 国产乱来视频区| av福利片在线| 精品亚洲乱码少妇综合久久| 亚洲成国产人片在线观看| 一区二区日韩欧美中文字幕| 日本av手机在线免费观看| 大陆偷拍与自拍| 国产黄色视频一区二区在线观看| 黄片无遮挡物在线观看| 日韩制服丝袜自拍偷拍| 久久久久国产一级毛片高清牌| 少妇被粗大的猛进出69影院| 久久综合国产亚洲精品| 亚洲经典国产精华液单| 9191精品国产免费久久| 肉色欧美久久久久久久蜜桃| 免费日韩欧美在线观看| 超色免费av| 乱人伦中国视频| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲视频免费观看视频| 久久久久久久精品精品| 精品国产乱码久久久久久男人| 久久久久精品性色| 如日韩欧美国产精品一区二区三区| 三上悠亚av全集在线观看| 成人国语在线视频| 日韩电影二区| 国产成人av激情在线播放| 国产成人精品在线电影| 日本欧美国产在线视频| freevideosex欧美| 欧美日韩亚洲高清精品| 日韩中文字幕视频在线看片| 亚洲视频免费观看视频| 久久 成人 亚洲| 寂寞人妻少妇视频99o| 天堂俺去俺来也www色官网| 一区二区av电影网| 日韩电影二区| 综合色丁香网| 国产av码专区亚洲av| 女性被躁到高潮视频| 亚洲色图 男人天堂 中文字幕| 十八禁高潮呻吟视频| 多毛熟女@视频| 不卡视频在线观看欧美| 人妻人人澡人人爽人人| 午夜日本视频在线| 自线自在国产av| 2018国产大陆天天弄谢| 不卡视频在线观看欧美| 高清欧美精品videossex| 婷婷色麻豆天堂久久| 亚洲成人一二三区av| 99久久中文字幕三级久久日本| 80岁老熟妇乱子伦牲交| 久久 成人 亚洲| 高清黄色对白视频在线免费看| 免费看av在线观看网站| 精品一区二区三区四区五区乱码 | 青春草亚洲视频在线观看| 国产av码专区亚洲av| 超碰97精品在线观看| 在线观看三级黄色| 五月开心婷婷网| 精品少妇一区二区三区视频日本电影 | 成人手机av| 日本欧美视频一区| 日日爽夜夜爽网站| 亚洲av中文av极速乱| 在线观看人妻少妇| 看免费av毛片| 天天躁夜夜躁狠狠久久av| 国产乱人偷精品视频| 久久综合国产亚洲精品| 国产成人精品久久二区二区91 | 精品一品国产午夜福利视频| 精品少妇内射三级| 91国产中文字幕| 精品一区二区三卡| 99久久综合免费| 国产精品成人在线| 日韩电影二区| 9191精品国产免费久久| 欧美精品一区二区免费开放| 精品一区二区三卡| 王馨瑶露胸无遮挡在线观看| 青草久久国产| 欧美97在线视频| 国产精品久久久久久av不卡| 一区二区三区精品91| 精品一区二区三卡| 精品卡一卡二卡四卡免费| 中文欧美无线码| 在线观看免费视频网站a站| 亚洲国产欧美网| 午夜av观看不卡| 久久婷婷青草| 一区二区三区乱码不卡18| 亚洲色图 男人天堂 中文字幕| 日韩一区二区视频免费看| 亚洲国产最新在线播放| 午夜av观看不卡| 80岁老熟妇乱子伦牲交| 国产在视频线精品| 日韩伦理黄色片| 亚洲第一青青草原| av女优亚洲男人天堂| 免费观看a级毛片全部| 男女无遮挡免费网站观看| 久久精品国产综合久久久| 哪个播放器可以免费观看大片| 少妇的逼水好多| 在现免费观看毛片|