海上小井眼脈沖空化射流發(fā)生器射流特性室內(nèi)試驗(yàn)研究

蘭起超,李根生,史懷忠,韋明輝,師帥帥

(油氣資源與探測國家重點(diǎn)實(shí)驗(yàn)室,中國石油大學(xué)(北京),北京 102249)

0 引 言

隨著油氣勘探開發(fā)越來越多的向深部地層、深海、低滲超低滲和低豐度邊際油氣方向發(fā)展，機(jī)械鉆速低、鉆井周期長、鉆井成本升高和井身結(jié)構(gòu)復(fù)雜等問題越來越突出[1-3]。如何有效解決這些問題成為國內(nèi)外鉆井工作者研究的重點(diǎn)[4-6]。已進(jìn)行的現(xiàn)場應(yīng)用表明水力脈沖空化射流技術(shù)是一種有效的提高機(jī)械鉆速的方法[7-8]。近年來，簡單的井身結(jié)構(gòu)——小井眼單通道井也已被證實(shí)可以高效開發(fā)邊際油田，尤其是海上邊際油田[9-11]。研發(fā)與之相配套的小尺寸水力脈沖空化射流工具對進(jìn)一步節(jié)約開發(fā)成本有重要意義。

到目前為止，針對水力脈沖空化射流發(fā)生器射流特性的試驗(yàn)研究，李根生[12]等通過現(xiàn)場試驗(yàn)地面測試立管壓力波動(dòng)曲線間接計(jì)算出了不同試驗(yàn)排量下的脈動(dòng)壓力、頻率和壓耗，得出了脈動(dòng)壓力、頻率和壓耗均隨排量的增大而升高的結(jié)論。馬東軍[13-14]等通過室內(nèi)實(shí)驗(yàn)和數(shù)值模擬研究了入口流量和葉輪葉片結(jié)構(gòu)對水力脈沖空化射流發(fā)生器產(chǎn)生的脈動(dòng)壓力振幅的影響，得出了3葉片葉輪(同4、6葉片相比)調(diào)制出脈沖射流的脈動(dòng)壓力振幅較大的結(jié)果。然而，壓力脈動(dòng)頻率、入出口壓力脈動(dòng)振幅和工具壓耗等對鉆井工具優(yōu)選、提速效果及水力參數(shù)設(shè)計(jì)有重要意義的參數(shù)特性的室內(nèi)測試與規(guī)律分析還未進(jìn)行，針對工具產(chǎn)生的空化射流特性的研究亦屬空白。根據(jù)研制成功的178、198和228mm 3種尺寸工具的現(xiàn)場應(yīng)用[15]可知，通過室內(nèi)試驗(yàn)測定不同排量下小尺寸工具不同射流參數(shù)的數(shù)值將對工具未來的現(xiàn)場應(yīng)用有重要的參考作用。鑒于此，本文通過室內(nèi)試驗(yàn)對小尺寸水力脈沖空化射流發(fā)生器的射流特性進(jìn)行比較全面的研究和分析，并揭示其鉆井提速機(jī)理。

1 室內(nèi)試驗(yàn)

1.1試驗(yàn)工具結(jié)構(gòu)及工作原理

小尺寸水力脈沖空化射流發(fā)生器的結(jié)構(gòu)主要包括殼體和芯子(導(dǎo)流體、葉輪、葉輪座、葉輪軸和自激振蕩噴嘴)，具體結(jié)構(gòu)如圖1所示。流體經(jīng)過裝置上部進(jìn)入導(dǎo)流體，經(jīng)導(dǎo)流體的斜坡通道改變流速和方向從而對葉輪葉片產(chǎn)生切向力使葉輪高速旋轉(zhuǎn)。葉輪通過軸安裝在葉輪座上，在流體對葉片沖擊力作用下，葉輪高速旋轉(zhuǎn)連續(xù)改變流道面積從而產(chǎn)生截?cái)嗍矫}沖射流，當(dāng)流體經(jīng)過自激振蕩噴嘴時(shí)，被自激振蕩腔室調(diào)制產(chǎn)生空化射流，最后兩種射流耦合成脈沖空化復(fù)合射流由噴嘴噴出。

1.2試驗(yàn)裝置

小尺寸(120mm)水力脈沖空化射流發(fā)生器；4個(gè)量程為5MPa(精度0.5級)的壓力傳感器；美國NI公司生產(chǎn)的PCI-6259數(shù)據(jù)采集卡及配套的數(shù)據(jù)采集系統(tǒng)(采樣時(shí)間：0.005～1s)；試驗(yàn)用泵為BQ700型撬裝泵(最大排量650L/min，最大工作壓力60MPa)，試驗(yàn)介質(zhì)為水。試驗(yàn)裝置連接如圖2所示。

(a) 試驗(yàn)工具結(jié)構(gòu)

(b) 試驗(yàn)工具照片

圖2 試驗(yàn)裝置連接簡圖

1.3試驗(yàn)方法

本試驗(yàn)用數(shù)據(jù)采集系統(tǒng)工作原理為：在一定排量下，采集測量點(diǎn)不同時(shí)刻的靜壓力。因此，為真實(shí)并完整描述壓力脈動(dòng)曲線規(guī)律，必須確定試驗(yàn)用采樣周期。由于最大排量下葉輪轉(zhuǎn)速最高，對應(yīng)脈動(dòng)曲線的周期最小，所以本實(shí)驗(yàn)采用最大排量下確定的采樣周期可以滿足小排量下測量數(shù)據(jù)的真實(shí)性。由于試驗(yàn)用泵排量的限制，試驗(yàn)最大排量設(shè)為10L/s。采樣周期的確定方法為：通過改變數(shù)據(jù)采集系統(tǒng)2個(gè)測量時(shí)間點(diǎn)的時(shí)間間隔(采樣時(shí)間)，得到不同采樣時(shí)間下的試驗(yàn)工具出口壓力脈動(dòng)曲線，當(dāng)在相同時(shí)間內(nèi)兩個(gè)脈動(dòng)曲線上的脈動(dòng)周期數(shù)和脈動(dòng)幅值相同時(shí)，認(rèn)為采集的數(shù)據(jù)穩(wěn)定性和精確度滿足要求。為保證數(shù)據(jù)采集系統(tǒng)采集工作穩(wěn)定，本次試驗(yàn)采樣周期采用兩種采樣時(shí)間中的最大者。

采樣周期下，測取不同排量下的數(shù)據(jù)，對工具產(chǎn)生的射流特性進(jìn)行室內(nèi)試驗(yàn)研究(由于采集系統(tǒng)采樣時(shí)間限制，該試驗(yàn)未研究自激振蕩噴嘴結(jié)構(gòu)等影響參數(shù))。

2 試驗(yàn)結(jié)果與分析

2.1試驗(yàn)工具脈動(dòng)頻率分析

井下脈沖射流發(fā)生裝置的使用對于井下儀器特別是帶有脈沖信號傳輸儀器的使用有較大影響，因此分析本試驗(yàn)工具產(chǎn)生的脈動(dòng)頻率對于井下工具的選型有重要意義。

在采樣時(shí)間分別取0.1、0.05、0.01和0.005s，泥漿泵排量設(shè)為10L/s時(shí)的試驗(yàn)條件下進(jìn)行脈沖噴射試驗(yàn)。其中采樣時(shí)間為0.01和0.005s時(shí)試驗(yàn)工具出口壓力P隨時(shí)間的脈動(dòng)關(guān)系如圖3所示。

由圖3可知，采樣時(shí)間為0.01和0.005s時(shí)，在相同的時(shí)間間隔(0.5s)內(nèi)出現(xiàn)了相同數(shù)量(4個(gè))的壓力脈動(dòng)周期，且壓力脈動(dòng)曲線和壓力脈沖幅值都比較近似，因此可以判斷10L/s排量下的采樣周期為0.01s。并且通過在0.5s的時(shí)間間隔內(nèi)出現(xiàn)4個(gè)壓力脈動(dòng)周期可計(jì)算出在10L/s的排量條件下，試驗(yàn)工具—小尺寸水力脈沖空化射流發(fā)生器的脈動(dòng)周期為0.125s，即脈動(dòng)頻率為8Hz。

采樣周期0.01s時(shí)，試驗(yàn)工具脈動(dòng)頻率隨排量變化關(guān)系曲線如圖4所示。

由圖4可知，小尺寸水力脈沖空化射流發(fā)生器的脈動(dòng)頻率與排量存在線性增大關(guān)系。這是因?yàn)樾〕叽缢γ}沖空化射流發(fā)生器內(nèi)的葉輪采用圓周均布的6葉片結(jié)構(gòu)，排量一定時(shí)葉輪基本有一個(gè)相對穩(wěn)定的轉(zhuǎn)速，當(dāng)排量增大一倍時(shí)，單位時(shí)間需要通過葉輪的流體體積增大一倍，則需要葉輪轉(zhuǎn)速相應(yīng)增大一倍，因此葉輪轉(zhuǎn)速與排量存在線性增大關(guān)系。當(dāng)葉輪旋轉(zhuǎn)一周時(shí)，葉輪葉片周期性的改變流體通過流道截面積6次，即測得脈動(dòng)頻率是葉輪轉(zhuǎn)數(shù)的6倍，所以脈動(dòng)頻率與排量也呈線性增大關(guān)系。

(a) 采樣時(shí)間為0.01s

(b) 采樣時(shí)間為0.005s

圖4 試驗(yàn)工具脈動(dòng)頻率隨排量變化關(guān)系曲線

2.2試驗(yàn)工具入口、出口壓力脈動(dòng)振幅分析

水力脈沖發(fā)生裝置調(diào)制的脈動(dòng)壓力能有效改善井底巖屑的壓持效應(yīng)，提高井底清巖效率進(jìn)而提高鉆速。而壓力脈動(dòng)振幅是影響清巖效率的重要參數(shù)，因此分析試驗(yàn)工具入口和出口壓力脈動(dòng)振幅具有重要意義。

采樣周期0.01s時(shí)，試驗(yàn)工具入口和出口壓力脈動(dòng)振幅pp1和pp2與排量關(guān)系曲線如圖5所示，排量為10L/s時(shí)試驗(yàn)工具入口和出口的壓力脈動(dòng)曲線如圖6所示。

由圖5可以看出，小尺寸水力脈沖空化射流發(fā)生器的入口和出口壓力脈動(dòng)振幅與排量之間均呈增大關(guān)系，且出口壓力脈動(dòng)振幅大于入口壓力脈動(dòng)振幅。這是因?yàn)殡S排量增大，輸送的流體動(dòng)能也隨之增大，導(dǎo)流體出口處流體速度和動(dòng)能也隨之增大，當(dāng)葉輪旋轉(zhuǎn)時(shí)，有更多的流體動(dòng)能不斷地轉(zhuǎn)化成流體壓能，導(dǎo)致工具出口脈動(dòng)壓力振幅也隨入口流量呈增大趨勢。由于自激振蕩噴嘴的自激振蕩腔室對流體的脈動(dòng)信號有放大的作用，所以產(chǎn)生出口壓力脈動(dòng)振幅大于入口壓力脈動(dòng)振幅的結(jié)果。

圖5 工具入口、出口壓力脈動(dòng)振幅與排量關(guān)系曲線

圖6 工具入口、出口壓力隨時(shí)間變化曲線

由圖6可以看出小尺寸水力脈沖空化射流發(fā)生器的入口和出口壓力脈動(dòng)曲線的趨勢是相同的，這是由于試驗(yàn)流體(水)的不可壓縮性導(dǎo)致的。

2.3試驗(yàn)工具壓耗分析

鉆井工具的壓耗對于鉆井水力參數(shù)設(shè)計(jì)有重要意義，而小尺寸水力脈沖空化射流發(fā)生器是在鉆進(jìn)單通道小井眼時(shí)采用，因此壓耗的測試更為關(guān)鍵。

采樣周期0.01s時(shí)，試驗(yàn)工具裝配芯子前后的壓耗pc與排量的關(guān)系曲線如圖7所示。

由圖7可以看出，裝配芯子的小尺寸水力脈沖空化射流發(fā)生器的壓耗與排量存在2次增大關(guān)系。在鉆井過程中，鉆井液流體通過噴嘴等流道收縮段時(shí)會(huì)由于阻力作用損耗一部分能量，產(chǎn)生壓耗，該壓耗可以用噴嘴壓降計(jì)算公式(1)計(jì)算得出。小尺寸水力脈沖空化射流發(fā)生器的壓耗等于以下2部分壓降之和：

(1) 工具自激振蕩噴嘴壓降的計(jì)算公式為：

(1)

式中：Δp為壓降，MPa；ρ為泥漿密度，g/cm3；Q為排量，L/s；C為噴嘴流量系數(shù)，取0.95；d為噴嘴當(dāng)量直徑，cm。

圖7 試驗(yàn)工具壓耗與排量關(guān)系曲線

(2) 葉輪間隙壓降：由于葉輪間隙是矩形，將其轉(zhuǎn)化成噴嘴圓形的當(dāng)量直徑模式，再用工具噴嘴壓降計(jì)算公式(1)計(jì)算。

因此，小尺寸水力脈沖空化射流發(fā)生器的壓耗與排量的關(guān)系為pc∝Q2，即存在2次增大關(guān)系。

2.4自激振蕩噴嘴射流特性

為了進(jìn)一步研究小尺寸水力脈沖空化射流發(fā)生器的工作原理，測試了只保留自激振蕩噴嘴(無導(dǎo)流體和葉輪總成)時(shí)工具出口的射流特性。

采樣周期0.01s，排量10L/s時(shí)，測得裝配自激振蕩噴嘴前后的工具出口壓力脈動(dòng)隨時(shí)間變化關(guān)系曲線如圖8所示，流體經(jīng)過有、無自激振蕩噴嘴試驗(yàn)工具時(shí)出口射流對比如圖9所示。

圖8 裝配自激振蕩噴嘴前后工具出口脈動(dòng)特性曲線對比

(a) 未安裝自激振蕩噴嘴出口射流 (b) 安裝自激振蕩噴嘴出口射流

由圖8可知，未裝配自激振蕩噴嘴的工具出口產(chǎn)生了正常的圓管近出口端的小幅正壓力(數(shù)值為測量點(diǎn)至圓管出口壓耗)，而裝配自激振蕩噴嘴的工具出口產(chǎn)生了比較明顯的負(fù)壓。這是因?yàn)樽约ふ袷幥皇业慕Y(jié)構(gòu)對流體的自激共振作用使流體內(nèi)部局部壓力降低，而工具出口端為大氣壓力，所以測得壓力表現(xiàn)為負(fù)壓力。

由圖9中安裝自激振蕩噴嘴前后射流特性對比可知，流體通過自激振蕩噴嘴后呈“發(fā)散狀”噴射，產(chǎn)生了較明顯的空化現(xiàn)象。這是由于自激振蕩腔室內(nèi)的負(fù)壓促生了流體內(nèi)部空泡的形成與發(fā)展，當(dāng)流體從出口噴射接觸大氣后，流體內(nèi)的大量空泡破裂使流體發(fā)散噴射，表現(xiàn)為空化現(xiàn)象。

3 提速機(jī)理分析

3.1脈沖射流提高清巖效應(yīng)

脈沖射流主要由葉輪旋轉(zhuǎn)周期性的改變流體流道面積(截?cái)嗍?產(chǎn)生，在鉆井過程中，與連續(xù)射流相比，井底的脈動(dòng)射流可以改變井底流場的壓力狀況，迫使巖屑翻轉(zhuǎn)、啟動(dòng)并及時(shí)脫離井底，提高井底的凈化程度，減少由于壓持效應(yīng)引起的巖石重復(fù)切削，進(jìn)而提高清巖效率。其力學(xué)分析如下文所述。

井底巖屑的壓持效應(yīng)模型如圖10所示，井底巖屑受鉆井液循環(huán)壓耗Δpa和靜液柱壓力pb共同作用，表示為pb+Δpa；獨(dú)立巖屑與巖石母體間的孔隙壓力pp。當(dāng)pb+Δpa>pp時(shí)，巖屑脫離井底必須受到水平推力的作用。該水平推力可近似由摩擦定律表示如下：

τi=μ·(pb+Δpa-pp+W)

(2)

式中：pb為靜液柱壓力，MPa；Δpa為鉆井液循環(huán)壓耗，MPa；pp為巖屑與巖石母體間的孔隙壓力，MPa；W為巖屑重量，kg；μ為摩擦系數(shù)。

圖10 井底巖屑的壓持效應(yīng)模型

當(dāng)脈沖射流使井底出現(xiàn)負(fù)壓脈沖Δpb時(shí)，井底巖屑所受的總壓持力(pb+Δpa-pp+Δpb+W)降低，由公式(2)可知巖屑脫離井底所需的橫向推力也降低。

井底負(fù)壓脈沖Δpb的絕對值與本試驗(yàn)工具產(chǎn)生的出口壓力脈動(dòng)振幅成正相關(guān)。因此，由公式(2)可知，出口壓力脈動(dòng)振幅越大，巖屑脫離井底所需的橫向推力越低，清巖效率也越高。

3.2空化射流提高破巖效率

空化射流主要由基于風(fēng)琴管諧振腔原理設(shè)計(jì)的自激振蕩噴嘴通過流體共振調(diào)制產(chǎn)生。對于空化現(xiàn)象，現(xiàn)在還難以給出一個(gè)簡明而嚴(yán)格的定義。一般把液體內(nèi)部局部壓力降低時(shí)，液體內(nèi)部或液固交界面上蒸汽或氣體的空穴(空泡)的形成、發(fā)展和潰滅過程稱為空化[16]。當(dāng)攜帶大量空泡的空化射流從鉆頭孔眼噴出到井底接觸巖石表面時(shí)，空泡會(huì)發(fā)生潰滅過程。由于空泡潰滅過程發(fā)生于瞬間(微秒級)，在局部產(chǎn)生極高的瞬時(shí)壓力作用在巖石表面附近，水流中不斷潰滅的空泡所產(chǎn)生的極高壓力引起的反復(fù)沖擊作用(也稱空蝕作用)，使巖石表面產(chǎn)生破壞，提高破巖效率。

本文所采用試驗(yàn)工具不僅巧妙地結(jié)合了脈沖射流與空化射流的技術(shù)優(yōu)勢，而且使兩種射流實(shí)現(xiàn)了互生互利的技術(shù)效果：截?cái)嗍矫}沖射流的負(fù)壓可以促進(jìn)空化的產(chǎn)生，而自激振蕩腔室能夠?qū)財(cái)嗍矫}沖射流產(chǎn)生放大作用，最終最大限度的釋放了流體能量，提高了破巖與清巖效率，進(jìn)而提高鉆井的機(jī)械鉆速。

4 結(jié) 論

利用室內(nèi)試驗(yàn)研究了小尺寸水力脈沖空化射流發(fā)生器的射流特性，并分析了其提速機(jī)理，得到如下結(jié)論：

(1) 小尺寸水力脈沖空化射流發(fā)生器能將連續(xù)射流調(diào)制成脈沖和空化復(fù)合射流，最大限度地釋放了流體能量，能有效提高鉆井的機(jī)械鉆速。

(2) 脈沖空化射流的脈動(dòng)頻率與排量呈線性增大關(guān)系，入口和出口壓力脈動(dòng)振幅隨排量增大而呈增大趨勢，工具壓耗與排量呈2次增大關(guān)系。

(3) 自激振蕩噴嘴的主要作用是產(chǎn)生負(fù)壓和空化，負(fù)壓有利于空穴的形成與發(fā)展，促生空化現(xiàn)象。

參考文獻(xiàn):

[1]沈忠厚.現(xiàn)代鉆井技術(shù)發(fā)展趨勢[J].石油勘探與開發(fā),2005,32(1): 89-91 .

Shen Zhonghou.Development trend of the modern drilling technology[J].Petroleum Exploration and Development,2005,32(1): 89-91 .

[2]管志川,李春山,周廣陳,等.深井和超深井鉆井井身結(jié)構(gòu)設(shè)計(jì)方法[J].石油大學(xué)學(xué)報(bào)(自然科學(xué)版),2002,25(6): 42-44.

Guan Zhichuan,Li Chunshan,Zhou Guangchen,et al.A method for designing casing program in deep and superdeep wells[J].Journal of the University of Petroleum (Edition of Natural Science),2002,25(6): 42-44.

[3]李根生,沈忠厚.高壓水射流理論及其在石油工程中應(yīng)用研究進(jìn)展[J].石油勘探與開發(fā),2005,32(1): 96-99.

Li Gensheng,Shen Zhonghou.Advances in researches and applications of water jet theory in petroleum engineering[J].Petroleum Exploration and Development,2005,32(1): 96-99.

[4]Cohen J H,Deskins G,Rogers J.High pressure jet kerf drilling shows fignificant potential to increase ROP.SPE 96557,2005.

[5]汪海閣,鄭新權(quán).中石油深井鉆井技術(shù)現(xiàn)狀與面臨的挑戰(zhàn)[J].石油鉆采工藝,2005,27(2):4-8.

Wang Haige,Zheng Xinquan.Status quo and faced challenges of deep well drilling techniques of PetroChina[J].Oil Drilling & Production Technology,2005,27(2):4-8.

[6]Kuru E,Okunsebor O M,Li Y.Hydraulic optimization of foam drilling for maximum drilling rate.SPE 91610,2004.

[7]史懷忠,李根生,沈忠厚,等.水力脈沖空化射流發(fā)生器配合Power V 防斜打快技術(shù)[J].石油鉆探技術(shù),2009,37(1):14-17.

Shi Huaizhong,Li Gensheng,She Zhonghou,et al.Application of hydraulic pulsation-cavitating generator with power V system in well fengcheng 1[J].Production Drilling Technology,2009,37(1): 14-17.

[8]王學(xué)杰,李根生,康延軍,等.利用水力脈沖空化射流復(fù)合鉆井技術(shù)提高鉆速[J].石油學(xué)報(bào),2009,30(1): 117-120.

Wang Xuejie,Li Gensheng,Kang yanjun,et al.Improvement of penetration rate with hydraulic pulsating-cavitation jet compound drilling technology[J].Acta Petrolei Sinica,2009,30(1):117-120.

[9]周建良,劉書杰,耿亞楠,等.單通道井鉆完井技術(shù)在邊際油田應(yīng)用前景分析[J].石油科技論壇,2012,37(3):13-17.

Zhou Jianliang,Liu Shujie,Geng Yanan,et al.Prospect analysis about application of monobore drilling & completion and in marginal oilfields[J].Petroleum Science and Technology Forum,2012,37(3):13-17.

[10]Duhen T,Sagot A,Kerbart Y.Deep offshore slim hole drilling[C]//Offshore Technology Conference,3-8May,1997,Texas.

[11]Fultz D J,Plttard J F,etc.Slim-Hole drilling in harsh environments[J].SPE Driling Engineering,1991: 225-229.

[12]李根生,史懷忠,沈忠厚,等.水力脈沖空化射流鉆井機(jī)理與試驗(yàn)[J].石油勘探與開發(fā),2008,35(2): 239-243.

Li Gensheng,Shi Huaizhong,Shen Zhonghou,et al.Mechanisms and tests for hydraulic pulsed cavitating jet assisted drilling[J].Petroleum Exploration and Development,2008,35(2): 239-243.

[13]馬東軍,李根生,史懷忠,等.水力脈沖空化射流發(fā)生器參數(shù)優(yōu)化試驗(yàn)研究[J].石油機(jī)械,2009,35(2): 239-243.

Ma Dongjun,Li Gensheng,Shi Huaizhong,et al.An experimental study of the parametric optimization of the hydro-pulse cavitation jet generator CPM[J].China Petroleum Machinery,2009,35(2): 239-243.

[14]馬東軍,李根生,史懷忠,等.水力脈沖射流特性室內(nèi)實(shí)驗(yàn)和數(shù)值模擬[J].實(shí)驗(yàn)流體力學(xué),2011,37(12): 9-11.

Ma Dongjun,Li Gensheng,Shi Huaizhong,et al.Study on hydraulic pulsed jet characteristic by in-house experiment and numerical simulation[J].Journal of Experiments in Fluid Mechanics,2011,37(12): 9-11.

[15]Fu J,Li G,Shi H,et al.A novel tool to improve the rate of penetration hydraulic pulsed cavitating-jet generator[J].SPE Drilling and Completion,2012,27(3): 355-362.

[16]易燦,李根生,陳日吉.圍壓下自振空化射流沖蝕性能實(shí)驗(yàn)研究[J].實(shí)驗(yàn)力學(xué),2005,20(2): 291-296.

Yi Can,Li Gensheng,Chen Riji,et al.Experimental study on self-oscillating cavitating jet erosion ability under ambient pressure[J].Journal of Experiments Mechanics,2005,20(2): 291-296.

作者簡介:

蘭起超(1989-)，男，黑龍江大慶人，碩士研究生。研究方向：油氣井流體力學(xué)與工程。通信地址：北京市昌平區(qū)中國石油大學(xué)(北京)石油工程教育部重點(diǎn)實(shí)驗(yàn)室(102249)。E-mail：lanqc0910@163.com