屈平 申瑞臣

中國石油鉆井工程技術(shù)研究院

煤層氣鉆井井壁穩(wěn)定機(jī)理及鉆井液密度窗口的確定

屈平 申瑞臣

中國石油鉆井工程技術(shù)研究院

在煤層中鉆井時,保持井壁的穩(wěn)定具有重要意義。為此,利用有限元軟件分析了近井壁裂紋尖端的受力狀況,比較了尖端與井壁兩處應(yīng)力集中的情況,結(jié)果表明,裂紋尖端受力高于井壁的應(yīng)力集中約一個數(shù)量級。從斷裂力學(xué)的角度,研究了近井壁裂紋的擴(kuò)展,分析了裂紋擴(kuò)展在節(jié)理煤巖井壁失穩(wěn)中所起的作用,討論了裂紋表面閉合、孔隙壓力及井壁滲透性等因素對井壁失穩(wěn)的影響。利用經(jīng)典斷裂力學(xué)理論即應(yīng)力強(qiáng)度因子理論和最大周向應(yīng)力準(zhǔn)則,對近井壁的裂紋進(jìn)行了抽象建模,建立了煤巖井壁失穩(wěn)的判據(jù)。根據(jù)失穩(wěn)判據(jù),提出拉剪破壞壓力和壓剪破壞壓力的概念,兩者分別為節(jié)理煤層鉆井液密度窗口的上下臨界值,兩者之間為井壁穩(wěn)定,兩者之外便是井壁失穩(wěn)。利用山西沁水盆地的巖心數(shù)據(jù)及現(xiàn)場鉆井資料進(jìn)行了實(shí)例計算,基于常規(guī)巖石力學(xué)方法求出的破裂壓力和坍塌壓力分別為19.1 MPa、6.4 M Pa,而用所建立的方法得到的上下臨界壓力為18.1 M Pa、7.6 M Pa。由此證明了所建立的方法得到的密度窗口小于前者,與實(shí)際情況更加相符。

斷裂力學(xué) 節(jié)理煤層 井壁穩(wěn)定 應(yīng)力強(qiáng)度因子 鉆井液密度窗口

在煤層中鉆井時,保持井壁的穩(wěn)定具有重要意義。由于煤巖及煤層具有強(qiáng)度低、割理多等特征,導(dǎo)致在煤層中鉆井時,特別是在煤層氣進(jìn)行多分支水平井時,井壁失穩(wěn)時常發(fā)生。據(jù)不完全統(tǒng)計,2006—2007年在煤層多分支水平井段內(nèi)發(fā)生垮塌埋卡鉆具13套,造成嚴(yán)重的經(jīng)濟(jì)損失。目前,已有學(xué)者初步探討了斷裂力學(xué)理論在節(jié)理煤層井壁穩(wěn)定方面的應(yīng)用[1-2],但都沒能提出具體的判據(jù)或方法。因此,無法定量的分析井壁穩(wěn)定。筆者試圖建立井壁失穩(wěn)的判據(jù),從而為確定鉆井液密度窗口、井眼軌跡與分支走向等鉆井施工參數(shù)提供定量的計算手段。

1 井壁穩(wěn)定力學(xué)機(jī)理

當(dāng)前井壁穩(wěn)定力學(xué)機(jī)理的研究主要從巖石力學(xué)角度進(jìn)行開展。然而對于節(jié)理地層,巖石力學(xué)方法不再適用,主要體現(xiàn)在:①無法準(zhǔn)確確定巖石的強(qiáng)度,由于裂紋的存在,導(dǎo)致巖石的強(qiáng)度存在尺寸效應(yīng);②節(jié)理地層不屬于連續(xù)性介質(zhì);③預(yù)測的破裂壓力和坍塌壓力與實(shí)際不符。針對節(jié)理性煤巖的特點(diǎn),筆者嘗試?yán)脭嗔蚜W(xué)的方法來研究節(jié)理煤巖的井壁失穩(wěn)問題。

1.1 基于斷裂力學(xué)的機(jī)理分析

在漫長的地質(zhì)演化條件下,巖石在各種各樣的構(gòu)造作用下產(chǎn)生了多種性質(zhì)的缺陷如斷層、節(jié)理和裂隙等,而煤巖的此種現(xiàn)象尤為典型。依據(jù)斷裂力學(xué)分析,巖石材料中的缺陷對巖石力學(xué)性質(zhì)有很大的影響,巖石變形破壞過程的實(shí)質(zhì)就是巖石材料中缺陷的萌生、擴(kuò)展、相互作用和貫通的過程。當(dāng)鉆遇節(jié)理地層時,一方面圍巖的地應(yīng)力將重新分布,使得近井壁應(yīng)力集中;另一方面彌散分布的微裂紋表面受到的應(yīng)力也將重新分布,裂紋尖端將出現(xiàn)應(yīng)力集中。因此,一方面當(dāng)井壁應(yīng)力超過圍巖的強(qiáng)度時,井壁圍巖破壞,進(jìn)而使得井壁失穩(wěn);另一方面裂紋尖端應(yīng)力超過巖石的斷裂韌度時,裂紋會進(jìn)一步擴(kuò)展失穩(wěn),達(dá)到一定條件時形成井壁失穩(wěn)坍塌。在上述兩種失穩(wěn)現(xiàn)象中,若占主導(dǎo)地位的是圍巖的強(qiáng)度破壞,則可用彈塑性理論推導(dǎo)的強(qiáng)度準(zhǔn)則來判斷井壁失穩(wěn);若是裂紋失穩(wěn)則需要用斷裂力學(xué)的理論來解決[3-7]。

筆者利用有限元軟件ANSYS對近井壁有裂紋存在時井壁的應(yīng)力分布狀態(tài)進(jìn)行了研究。研究表明:近井壁裂紋的尖端應(yīng)力集中遠(yuǎn)高于井壁的應(yīng)力集中,約高一個數(shù)量級,相比原始地應(yīng)力約高兩個數(shù)量級。另外,巖石斷裂力學(xué)研究也表明在遠(yuǎn)低于圍巖強(qiáng)度的條件下,裂紋就會失穩(wěn)擴(kuò)展。材料斷裂力學(xué)研究表明脆性材料的斷裂破壞多屬于低應(yīng)力破壞。因此,可以假定裂紋失穩(wěn)擴(kuò)展是造成含裂紋地層井壁破壞的主因。以斷裂力學(xué)的觀點(diǎn),對于不含彌散分布裂紋的地層如砂巖等,井壁失穩(wěn)破壞可分為4個環(huán)節(jié):①在應(yīng)力的作用下,巖石的局部缺陷產(chǎn)生微裂紋;②若微裂紋尖端的應(yīng)力增加到一定程度,彌散分布的微裂紋就會失穩(wěn)擴(kuò)展;③若應(yīng)力加大,則微裂紋就會串接為宏觀裂紋;④若應(yīng)力再大則宏觀裂紋就會演化至災(zāi)難性失穩(wěn)裂紋,此時就會造成井壁的失穩(wěn)(圖1)。對于含彌散分布裂紋的地層,井壁失穩(wěn)破壞就只有后3個環(huán)節(jié)。對于節(jié)理煤層而言,相比常規(guī)的含裂紋地層如碳酸鹽巖地層,裂紋更加發(fā)育。煤層中存在大量相互正交的面割理和端割理(統(tǒng)稱為微裂紋),微裂紋之間基本處于貫通狀態(tài)。因此,煤層井壁失穩(wěn)破壞只有后兩個環(huán)節(jié)[8]。

圖1 不含裂紋的地層的井壁失穩(wěn)示意圖

1.2 力學(xué)機(jī)理影響因素

1.2.1 裂紋表面閉合分析

唐立強(qiáng)[2]認(rèn)為井壁附近的巖石中,微裂紋表面始終閉合,只能產(chǎn)生裂紋面的滑動擴(kuò)展。然而,當(dāng)一個塊體被放在粗糙的面上時,則該塊體與該粗糙面實(shí)際接觸的表面積幾乎等于零,這時,全部接觸力是被幾個或更多的接觸點(diǎn)支持著,即點(diǎn)接觸[3]。在增加法向荷載的情況下,點(diǎn)的接觸可能因彈性變形、壓碎和張裂紋等而增大[3]。事實(shí)上由于雜質(zhì)的充填、地層水的存在和點(diǎn)接觸的影響,地層的裂紋不可能完全閉合。特別是對于煤層中的割理,通過試井和測井等方法,已經(jīng)得出割理具有很高的滲透率。煤巖中的地層水和煤層氣正是通過沒有閉合的裂紋而不斷涌入井筒排到地面的。近井壁裂紋主要為壓剪或拉壓復(fù)合型裂紋。裂紋表面是否閉合對井壁失穩(wěn)機(jī)理有很大的影響,主要體現(xiàn)在以下幾點(diǎn):①若裂紋表面完全閉合則裂紋只能產(chǎn)生裂紋面的滑動擴(kuò)展,此時裂紋的Ⅰ型應(yīng)力強(qiáng)度因子為零,在計算裂紋的擴(kuò)展條件時只需考慮Ⅱ型應(yīng)力強(qiáng)度因子[8];②若裂紋表面不閉合則裂紋會受到滑動、拉裂或壓裂的復(fù)合作用,近井壁裂紋的失穩(wěn)破壞機(jī)理相對復(fù)雜;③裂紋面不閉合時還需要考慮孔隙壓力、鉆井液的滲漏對井壁失穩(wěn)的影響。考慮到巖石的破壞主要為剪應(yīng)力破壞。因此,在用有限元模擬井壁失穩(wěn)時可假定裂紋完全閉合。而在考慮孔隙壓力、鉆井液的滲漏以及應(yīng)力敏感性時需要認(rèn)定裂紋未完全閉合。

1.2.2 孔隙壓力

地層中一般會有孔隙水存在,特別是煤層中存在孔隙水的可能性較大。在上覆巖層的壓力下,孔隙水會形成孔隙壓力,又稱地層壓力。存在孔隙壓力時,孔隙壓力會改變裂紋上下兩表面所受的正應(yīng)力,但不影響剪應(yīng)力。下面根據(jù)裂紋表面的受力情況來分別討論。若裂紋表面受到的是拉應(yīng)力,當(dāng)孔隙壓力存在時,裂紋表面的正應(yīng)力為:

式中σ為裂紋表面的正應(yīng)力,M Pa;σ1為不考慮孔隙壓力時裂紋表面的正應(yīng)力,M Pa;p為孔隙壓力。

顯然孔隙壓力會加重裂紋的受力狀況,增加裂紋的Ⅰ型應(yīng)力強(qiáng)度因子,使裂紋更容易擴(kuò)展破壞。若裂紋表面受到的是壓應(yīng)力,當(dāng)孔隙壓力存在時,裂紋表面的正應(yīng)力為:

當(dāng) p≤σ1時,孔隙壓力會緩解裂紋的受力狀況,減小裂紋的Ⅰ型應(yīng)力強(qiáng)度因子,裂紋趨于穩(wěn)定。當(dāng) p> σ1時,孔隙壓力將使裂紋從受壓轉(zhuǎn)為受拉,轉(zhuǎn)而使Ⅰ型應(yīng)力強(qiáng)度因子增加。在地層中裂紋大多受壓剪作用,孔隙壓力將緩解裂紋的受壓作用。因此,當(dāng)煤層氣井的開采井排水采氣時,近井壁的孔隙壓力將逐漸降低,孔隙壓力對裂紋受壓的緩解作用將逐漸減小,裂紋可能由穩(wěn)定轉(zhuǎn)向失穩(wěn)。這可以解釋在排水采氣過程中為什么會有大量的煤粉出現(xiàn),以及滲透率越來越低的現(xiàn)象。

1.2.3 井壁滲透性影響

井壁有無滲透性對近井壁的孔隙壓力有很大的影響。當(dāng)井壁無滲透性時只需考慮原始孔隙壓力即可。當(dāng)井壁有滲透性時,無論是在欠平衡還是在過平衡狀態(tài)下,近井壁地層的孔隙壓力都會受到井底壓力的干擾而發(fā)生改變,由此對井壁的失穩(wěn)產(chǎn)生影響。下面從欠平衡和過平衡兩個方面來分析。欠平衡鉆井時,近井筒的地層流體將向井筒中滲流。由于滲流速度差異,隨著時間的推移,近井筒的孔隙壓力逐漸降落到井筒中的壓力水平[4]。過平衡鉆井時,井筒中的鉆井液將向近井筒的地層滲流。同樣由于滲流速度的差異,隨時間的推移,近井筒的地層會逐漸上升到井底壓力水平。根據(jù)井底壓力和孔隙壓力,利用達(dá)西滲流公式可以求出隨時間變化的近井筒壓力的變化趨勢。把改變后的孔隙壓力代入式(1)和式(2)中就可以分析滲透性對井壁的穩(wěn)定性,具體過程略去。據(jù)此可得:過平衡鉆井時,鉆井液的滲漏會改變近井壁的孔隙壓力,從而間接影響裂紋的擴(kuò)展與穩(wěn)定。由于鉆井液滲漏具有時間延遲效應(yīng)。因此,使得節(jié)理地層井壁失穩(wěn)具有時間效應(yīng)。

2 井壁失穩(wěn)判據(jù)

2.1 假設(shè)條件

假設(shè):近井壁裂紋擴(kuò)展的開始是井壁失穩(wěn)的臨界點(diǎn)。即如果存在井壁裂紋的擴(kuò)展條件,就認(rèn)為井壁會發(fā)生失穩(wěn)坍塌現(xiàn)象,并且近井壁裂紋擴(kuò)展的開始,就是井壁失穩(wěn)的開始。根據(jù)此假設(shè),就可以把定量確定裂紋擴(kuò)展的判據(jù)應(yīng)用于確定井壁失穩(wěn)。

2.2 裂紋擴(kuò)展判據(jù)

國內(nèi)外學(xué)者對裂紋的失穩(wěn)擴(kuò)展研究取得了豐碩的成果,建立了相應(yīng)的斷裂準(zhǔn)則[6-7]。通過適用性分析,筆者選用了基于巖石破壞準(zhǔn)則的巖石斷裂判據(jù)[6]。筆者利用D-P破壞準(zhǔn)則,并結(jié)合裂紋尖端的應(yīng)力場、應(yīng)力強(qiáng)度因子和斷裂韌度提出了該判據(jù),具體推導(dǎo)過程見本文參考文獻(xiàn)[6]。

通過推導(dǎo)可得Ⅰ—Ⅱ壓剪裂紋擴(kuò)展判據(jù):

2.3 應(yīng)力強(qiáng)度因子計算

應(yīng)力強(qiáng)度因子是上述壓剪判據(jù)中最重要的兩個參數(shù)。因此,需要單獨(dú)討論其求解方法。應(yīng)力強(qiáng)度因子為斷裂力學(xué)中用來表征裂紋應(yīng)力集中程度的參數(shù),應(yīng)力強(qiáng)度因子越大,表明裂紋尖端的應(yīng)力集中越大,反之亦然。

2.3.1 水平井中的應(yīng)力強(qiáng)度因子計算

本文參考文獻(xiàn)[8]中對水平井近井壁裂紋的應(yīng)力強(qiáng)度因子計算進(jìn)行了詳盡分析,本文略去了推導(dǎo)過程。

2.3.2 直井中的應(yīng)力強(qiáng)度因子計算

為了能計算出直井近井壁裂紋的應(yīng)力強(qiáng)度因子,提出兩條假設(shè):①裂紋在井筒的某個過井軸的截面內(nèi);②假定裂紋為貫穿性裂紋。

圖2 直井近井壁裂紋應(yīng)力強(qiáng)度因子計算簡圖

把式(4)代入井周應(yīng)力計算式求得σr、σφ和τrφ,再利用σr、σφ和τrφ求得裂紋面的正應(yīng)力σ和剪應(yīng)力τ,再把σ和τ代入應(yīng)力強(qiáng)度因子計算式[7-8]就可求得應(yīng)力強(qiáng)度因子 KⅠ、KⅡ。

3 鉆井液密度窗口的確定

確定鉆井液密度窗口是井壁穩(wěn)定研究的一項重要任務(wù)。確定鉆井液密度窗口即確定了鉆井液的最大或最小密度,也就確定了保持井壁穩(wěn)定的最大或最小井底壓力。

利用巖石力學(xué)方法可以確定出鉆井液密度窗口即確定坍塌壓力(pc)和破裂壓力(pf)。與之對應(yīng),從斷裂力學(xué)角度分析,也可以確定出兩個井底壓力值(圖3)。但與前者無論在意義上或計算方法上都大不同。當(dāng)井底壓力過小時,裂紋受到壓應(yīng)力和剪應(yīng)力的共同作用,導(dǎo)致裂紋失穩(wěn)擴(kuò)展,由此使井壁破壞,此時的井底壓力可記為 pfc,名為壓剪破壞壓力。當(dāng)井底壓力過大時,裂紋主要受到拉應(yīng)力和剪應(yīng)力的作用,導(dǎo)致裂紋擴(kuò)展,進(jìn)而井壁破壞,此時的井底壓力可記為 pff,名為拉剪破壞壓力。顯然,若破裂壓力 pf>pff,則地層不存在破裂壓力,井壁失穩(wěn)時的上臨界壓力為 pff。若破裂壓力 pf≤pff,則地層存在破裂壓力,井壁失穩(wěn)時的上臨界壓力即為 pf。其實(shí),有裂紋的地層中往往有 pf>pff,這是導(dǎo)致壓裂失效的重要原因。若坍塌壓力 pc>pfc,則井壁在裂紋尚未破壞時已失穩(wěn),此時井壁失穩(wěn)的下臨界壓力為 pc。若坍塌壓力 pc≤pfc,則井壁的破壞由裂紋造成,此時井壁失穩(wěn)的下臨界壓力為 pfc。因此,當(dāng)從巖石力學(xué)和斷裂力學(xué)角度共同考慮時,鉆井液密度窗口為[max(pc,pfc)、min(pf,pff)]。如此確定的密度窗口更加準(zhǔn)確也全面。

圖3 從斷裂力學(xué)角度確定井壁失穩(wěn)上下臨界值方法示意圖

3.1 斷裂力學(xué)的臨界壓力計算方法

當(dāng)近井壁只有一條裂紋,利用第3節(jié)中的井壁失穩(wěn)判據(jù),通過迭代的方法就可以求出[pfc、pff],見圖3。

近井壁存在多條裂紋時,鉆井液密度窗口的計算方法非常復(fù)雜,由于篇幅在此略去。但其主要的研究思路等同于近井壁只有單裂紋的情況。

3.2 實(shí)例分析

其他參數(shù):孔隙度 f=3.5%;抗拉強(qiáng)度St=1.15 M Pa;有效應(yīng)力系數(shù)為1;非線性修正系數(shù)為1。

巖石力學(xué)方法求出的破裂壓力和坍塌壓力[10-11]分別為 pf=19.1 M Pa、pc=6.4 M Pa。利用本文中的斷裂力學(xué)方法計算求得的上下臨界壓力分別為 pff= 18.1 M Pa、pfc=7.6 M Pa,計算結(jié)果見圖4。綜合兩種方法,鉆井液的密度窗口為[7.6 M Pa、18.1 M Pa]。

圖4 鉆井液上下臨界密度計算結(jié)果圖

4 結(jié)論

1)裂紋的失穩(wěn)擴(kuò)展是導(dǎo)致節(jié)理地層井壁失穩(wěn)的主因,在節(jié)理煤層中可把裂紋擴(kuò)展可視為井壁失穩(wěn)的充分條件。

2)孔隙壓力會改變裂紋的受力狀況,影響裂紋的擴(kuò)展。鉆井液的滲漏會改變近井壁的孔隙壓力,從而間接影響裂紋的擴(kuò)展與穩(wěn)定。由于鉆井液滲漏具有時間延遲效應(yīng)。因此,使得節(jié)理地層井壁失穩(wěn)具有時間效應(yīng)。

3)基于斷裂力學(xué)的井壁穩(wěn)定模型同樣可確定地層的密度窗口。在密度窗口內(nèi),裂紋不擴(kuò)展,井壁保持穩(wěn)定;在密度窗口外,裂紋擴(kuò)展,井壁失穩(wěn)。

[1]鮮保安.煤層氣田多分支井優(yōu)化設(shè)計研究[D].北京:中國石油大學(xué),2007.

[2]唐立強(qiáng),楊敬源,王勇,等.井壁穩(wěn)定性的斷裂損傷力學(xué)分析[J].哈爾濱工程大學(xué)學(xué)報,2007,28(6):642-646.

[3]黃醒春,陶連金,曹文貴.巖石力學(xué)[M].北京:高等教育出版社,2005.

[4]屈平,申瑞臣,袁進(jìn)平,等.煤儲層的應(yīng)力敏感性理論研究[J].石油鉆探技術(shù),2007,35(5):68-71.

[5]李建林,孫志宏.節(jié)理巖體壓剪斷裂及其強(qiáng)度研究[J].巖石力學(xué)與工程學(xué)報,2000,19(4):444-448.

[6]周家文,徐衛(wèi)亞,石崇.基于破壞準(zhǔn)則的巖石壓剪斷裂判據(jù)研究[J].巖石力學(xué)與工程學(xué)報,2007,26(6):1194-1201.

[7]程靳,趙樹山.斷裂力學(xué)[M].北京:科學(xué)出版社,2006.

[8]屈平,申瑞臣,楊恒林,等.節(jié)理煤層井壁穩(wěn)定性的評價模型[J].石油學(xué)報,2009,30(3):455-459.

[9]朱煥春,BRUMM ER RICHARD,ANDRIEUX PA TRICK.節(jié)理巖體數(shù)值計算方法及其應(yīng)用(一):方法與討論[J].巖石力學(xué)與工程學(xué)報,2004,23(20):3444-3449.

[10]鄧金根,張洪生.鉆井工程中井壁失穩(wěn)的力學(xué)機(jī)理[M].北京:石油工業(yè)出版社,1998.

[11]呂開河,孫明波,邱正松.塔里木盆地依奇克里克區(qū)塊煤層鉆井技術(shù)研究[J].石油學(xué)報,2006,27(5):108-111.

Mechan ism of wellbore stability and determ ination of drilling fluid density window in coalbed methane drilling

Qu Ping,Shen Ruichen
(D rilling Engineering Technology Research Institute,CN PC,Beijing 100195,China)

NATUR.GAS IND.VOLUM E 30,ISSUE 10,pp.64-68,10/25/2010.(ISSN 1000-0976;In Chinese)

It is impo rtant to keep the wellbore stable during coalbed methane drilling.Here in this paper the finite element analysis softw are is adop ted to study the stress condition at the tipsof the near-wellbo re cracks,w hich is compared w ith that at the wellbo re. The result show s that the fo rmer stress is higher than the latter one by amagnitude order.From the perspective of the fracturemechanics,the development of the near-wellbore cracks and its bad effects on wellbore instability in coal seam are both studied;and many other influencing facto rson wellbo re instability like surface closed crack,po re p ressure,and permeability of a wellwall,etc., are also discussed.Mo reover,the classical fracture mechanics theo ries of stress intensity facto r and the maximum circum ferential stress p rinciple are used in abstractmodelling for near-wellbo re cracks,and then the criteria for wellbo re instability in coal seam is developed,based on w hich the concep tsof pull-shear collapsing p ressure and comp ression-shear collapsing p ressure are thus put forward.The two concep ts are the upper and lower critical limitsof the drilling fluid density w indow of the coal seam.If between these two values,the wellbo re w ill be stable;if not,the wellbo re w ill be unstable.Calculation example ismade by using the co re data and the field data taken from the Qinshui Basin of Shanxi.The fracture p ressure and the collapse p ressure obtained by the conventional rock mechanicsmethod are 19.1 and 6.4 M Pa respectively,w hile those obtained by the above-mentioned method are 18.1 and 7.6 M Pa respectively.This demonstrates that the density w indow obtained by this new method is smaller,show ingmore agreementw ith the actual situation.

fracturemechanics,coal seam,wellbore stabilizing,stress intensity factor,drilling fluid density w indow

屈平等.煤層氣鉆井井壁穩(wěn)定機(jī)理及鉆井液密度窗口的確定.天然氣工業(yè),2010,30(10):64-68.

DO I:10.3787/j.issn.1000-0976.2010.10.016

國家科技重大專項“煤層氣鉆井工程技術(shù)及裝備研制”(編號:2008ZX05000-036-001)。

屈平,1982年生,工程師,中國石油勘探開發(fā)研究院博士研究生;主要從事煤層氣鉆井技術(shù)研究。地址:(100195)北京市海淀區(qū)四季青鎮(zhèn)北塢村路甲25號靜芯園小區(qū)K座302室。電話:(010)52781873。E-mail:qp1012@163.com

2010-05-29 編輯 鐘水清)

DO I:10.3787/j.issn.1000-0976.2010.10.016

Qu Ping,engineer,born in 1982,is studying fo r a Ph.D degree at the Exp lo ration&Development Research Institute,CNPC,doing research on coalbed methane drilling technology.

Add:Room 302,Building K,Jingxinyuan Residential A rea,No.Jia 25,Beiw ucun Rd.,Sijiqing,Haidian District,Beijing 100195, P.R.China

Tel:+86-10-5278 1873 E-mail:qp1012@163.com