高溫高壓氣井、凝析氣井井筒及近井地帶的壓力變化規(guī)律

常志強 肖香姣 朱忠謙 王小培 王海應

中國石油塔里木油田公司勘探開發(fā)研究院

高溫高壓氣井、凝析氣井井筒及近井地帶的壓力變化規(guī)律

常志強 肖香姣 朱忠謙 王小培 王海應

中國石油塔里木油田公司勘探開發(fā)研究院

高溫高壓氣井的井筒溫度與壓力的分布是氣藏開發(fā)方案制定和調整的重要數據。為此，在井筒摩阻室內模擬實驗及高溫井筒溫度室內模擬實驗的基礎上，優(yōu)選了高溫高壓天然氣PVT物性參數計算方法，開展了高溫高壓氣井、凝析氣井井筒及近井帶溫度、壓力分布規(guī)律的理論研究，解決了該類氣井的動態(tài)監(jiān)測難題；建立了考慮井流物組成、流體相態(tài)、含水量、井身結構、動能損耗、重烴含量、環(huán)境溫度等多種因素影響的井筒壓力溫度耦合計算模型；編制了高溫高壓氣井、凝析氣井井筒動力學軟件，可實時獲取氣井壓力、產能變化規(guī)律及儲層動態(tài)參數信息。現場應用結果表明，這為實時了解氣田及單井的生產動態(tài)、及時進行平面產氣結構優(yōu)化、制定合理開發(fā)技術對策提供了技術手段。

塔里木盆地 高溫高壓氣井 凝析氣井 井筒 滲流 相態(tài) 凝析水 變化 數學模型

1 氣井井筒壓力、溫度計算模型

筆者綜合利用Ramey的假設條件，考慮環(huán)境溫度和井身結構的影響建立了分段的井筒溫度計算模型。同時，結合能量守恒方程與狀態(tài)方程，基于Cullender－Smith方法［1－3］，建立了考慮流體相態(tài)、動能損耗、動態(tài)摩阻等因素影響的井筒壓力計算模型。將井筒壓力與溫度計算模型通過程序耦合，確保獲取準確的井筒壓力、溫度變化規(guī)律如圖1所示。

1.1 井筒壓力計算模型

假設：①氣體在氣井中的流動為穩(wěn)定流動；②在氣體流動過程中，沒有功的輸入和輸出。根據質量守恒定律，推導得到考慮動能損耗的氣井井筒壓力計算數學模型為：

其中，相關參數的計算限于篇幅，此處不再詳細介紹，可詳見本文參考文獻［3］。

1.2 井筒溫度計算模型

在地層溫度分布研究過程中，要進行相關假設。根據熱量守恒定律，考慮第二接觸面（水泥環(huán)和地層的交界面）向周圍地層傳遞的熱量服從Ramey無因次時間函數，且井筒向第二接觸面?zhèn)鬟f的熱量與第二接觸面向周圍地層傳遞的熱相等，可推導得到高溫高壓氣井、凝析氣井井筒溫度計算模型為：

圖1 氣井井筒壓力、溫度計算流程示意圖

1.3 相關參數確定

對于凝析氣井，由于井筒內壓力和溫度的變化較大，當壓力和溫度沿著井筒向上不斷下降時，井筒中的液體則不斷增加，氣、液組成則不斷發(fā)生變化，特別是對于富凝析氣井，然而常規(guī)方法并未考慮井筒相態(tài)變化。因此要考慮兩相偏差因子、熱物性參數以及動態(tài)摩阻對井筒壓力計算模型的影響。由于篇幅有限，具體計算模型見本文參考文獻［4－7］。

2 氣井近井地帶壓力、溫度計算模型

氣井開采時，氣體從高壓到低壓膨脹吸熱，使得產層的溫度會緩慢的降低，其溫度的傳播規(guī)律同樣滿足熱傳導方程，但其擴散的速度明顯要小于壓力擴散的速度。尤其對高溫高壓氣井，開關井瞬間，近井帶溫度及壓力隨時間會發(fā)生明顯變化。如果不考慮溫度變化的影響，將嚴重影響測試資料準確性，難以獲取精度較高的測試資料及儲層動態(tài)參數信息［8－12］。

2.1 模型建立

2.1.1 近井帶壓力計算模型

在地層壓力分布研究過程中，主要基于以下假設：①地層為均勻各向同性地層；②地層中的流體為定常黏度的牛頓流體，孔隙度是常數；③流體在地層中的流動為層流狀態(tài)，遵從達西定律；④流體是單相氣體，流動是一個等溫過程；⑤忽略重力作用，不考慮其他物理化學的影響。

可得無量綱壓力的控制方程為：

初始條件：

內邊界條件：

外邊界條件，即無限大邊界條件：

定壓邊界條件：

封閉邊界：

2.1.2 近井帶溫度計算模型

在地層溫度分布研究過程中，要進行相關假設。

所得到無量綱溫度控制方程為：

初始條件：

內邊界條件：

外邊界條件：

2.2 模型求解

原設計中，超前管棚每2.0 m高打一環(huán)，環(huán)向間距0.5 m，每環(huán)共57根，每根長3.0 m，共需1 160根，長3 480 m。實際施工中，采用本施工方案從地面自上而下一次性鉆孔打進，深入地下38.4 m，間距0.5 m，共57個鉆孔需無縫鋼管約2 050 m，比原設計少1 274 m長，節(jié)省超前管棚支護（灌漿）投資約35萬元。本出渣方案省去調壓井口門機、鋼絲繩和吊罐等吊裝設備，減少門機司機、信號員等特種作業(yè)人員，減少施工資源投入。

采用有限元方法求解，例如溫度模型求解，利用伽遼金加權余量法，令其權函數為差值函數：Ni＝ai＋bix＋ciy，其中i＝1，2，3，然后將方程弱化，最終得到方程組為：

求解系統(tǒng)方程組，就可得任意網格點（x，y）上n＋1時刻溫度值T（x，y）。

3 軟件編制及實例分析

3.1 軟件編制

在以上井筒及近井帶溫度、壓力分析理論與方法基礎上，采用VC＋＋編程語言結合OLE數據庫的方法，編制了高溫高壓氣井、凝析氣井井筒及近井地帶壓力、溫度計算軟件。結合實測數據擬合得到地層的熱力參數，實時獲取井筒及井筒附近壓力、溫度分布及變化規(guī)律，解決了高溫、高壓氣井動態(tài)監(jiān)測難題，為高溫、高壓氣藏開發(fā)調控與生產管理提供了科學的技術支撐和有效的輔助工具。軟件基本界面如圖2所示。

3.2 實例分析

3.2.1 模型精度分析

為驗證軟件的可靠性，筆者利用克拉2、牙哈、吉拉克等區(qū)塊30余井次的實測數據進行了對比分析。從誤差分析結果表可以看出（表1），溫度最大相對誤差不超過2%，壓力最大相對誤差不超過1.8%，滿足工程計算精度要求。KL205、KL2－8井計算壓力剖面與實測壓力數據擬合的較好，說明所建立的模型實用性較強。

圖2 高溫氣井及凝析氣井壓力、溫度計算軟件基本界面圖

表1 JLK、YH等井8組實測數據與計算數據誤差分析表

3.2.2 溫度、壓力計算敏感參數分析

3.2.2.1 井身結構對溫度壓力分布的影響

從圖3可以看出，考慮井身結構時，相當于減小了井底處水泥環(huán)的半徑，而水泥環(huán)的導熱系數比地層小，井筒向近井地帶熱量傳遞更快，由此溫度會降低。從圖4可以看出井身結構對壓力幾乎沒有影響。

3.2.2.2 含水量對井筒壓力溫度分布的影響

從圖5可以看出，隨著含水量增加，從井底算到井口的溫度越大。這是由于含水量增大相當于增加了氣井的流量，流量越大，平均流速越快，這樣氣井與周圍地層的熱交換越不充分，熱量損失越少，溫度越高。從圖6中可以看出，含水量越大，井底算到井口的壓力越低。這是由于在建立擬單向流的壓力計算模型時，加入了含水校正系數（Fw），含水量越大，Fw越大，計算的壓力越低。另一方面，從圖6中可以看出，含水量越高，溫度越低，井筒中的密度越小，因此壓力越小。

圖3 井身結構對溫度的影響圖

圖4 井身結構對壓力的影響圖

圖5 含水量對溫度的影響圖

圖6 含水量對壓力的影響圖

3.2.2.3 重烴組分對壓力溫度分布的影響

從圖7可以看出重烴含量越高，從井底計算到井口的溫度越大。從圖8可以看出重烴組分含量越高，從井底算到井口的壓力越小。

圖7 重烴組分對溫度的影響圖

圖8 重烴組分對壓力的影響圖

4 結論

1）首次建立了綜合考慮井流物組成、流體相態(tài)、含水量、井身結構、動能損耗、重烴含量、環(huán)境溫度等因素影響的井筒壓力溫度耦合計算模型，初步形成了高溫高壓氣井井筒溫度、壓力計算理論分析方法。通過20井次實測數據驗證，模型計算精度較高。利用建立的井筒壓力溫度耦合計算模型，完成了15井次高溫高壓氣井的井筒溫度、壓力剖面預測，獲取了準確的壓力、產能及儲層動態(tài)參數信息，解決了該類氣井動態(tài)監(jiān)測難題，為克拉2氣田調整方案編制、迪那2、大北氣田開發(fā)方案編制提供了科學依據。

2）建立了井筒附近地層溫度、壓力變化的理論模型，獲取了高溫高壓氣井近井帶壓力和溫度分布及其變化規(guī)律。提出了動井筒溫度的概念，針對于產層與非產層的區(qū)別，給出了產層及非產層溫度的計算方法，獲取了準確的近井帶地層溫度剖面。

3）編制了高溫高壓氣井井筒動力學軟件，利用軟件完成了35井次實例分析，計算的溫度、壓力誤差均低于2%，能夠滿足工程要求，軟件具有實用性。

4）高溫高壓氣井近井帶及井筒壓力溫度耦合計算模型的建立為實時了解氣田及單井生產動態(tài)、及時進行平面產氣結構優(yōu)化、制定合理開發(fā)技術對策提供了重要的技術手段，為高溫高壓氣藏開發(fā)調控與生產管理提供了科學的技術支撐和有效的輔助工具。該配套技術目前已在塔里木油田克拉2、迪那2、大北等氣田進行了應用，推廣應用前景非常廣闊。

符 號 說 明

pD為無因次壓力；CD為無因次等溫壓縮系數；S為表皮系數；ww為水的質量流量，kg／m3；qg為氣體流量，m3／d；qo為油流量，m3／d；Zmix為復合氣體偏差因子；T為體系平衡溫度，K；d為油管內徑，m；f為Moody摩阻系數，用Jain公式計算；xD、yD分別為無因次距離；pwf為井底流動壓力，MPa；ptf為井口流動壓力，MPa；p為體系平衡壓力，MPa；Fw為含水校正系數；wg為氣體的質量流量，kg／m3；γmix為復合氣體相對密度；γg為干氣的相對密度；γL為液體的相對密度。

［1］李士倫.天然氣工程［M］.北京：石油工業(yè)出版社，2000：38－80.

［2］常志強.考慮多相流體復雜相態(tài)變化的凝析氣井生產動態(tài)分析［D］.南充：西南石油學院，2003.

［3］肖香嬌，姜漢橋，劉曰武.高壓氣井井筒附近地層溫度分布計算方法［J］.科技導報，2011，29（32）：20－24.

［4］何麗萍，毛美麗，廖紅梅，等.蘇里格氣田壓力監(jiān)測方法和遠程試井技術［J］.天然氣工業(yè)，2011，31（1）：59－61.

［5］朱德武，何漢平.凝析氣井井筒溫度分布計算［J］.天然氣工業(yè)，1998，18（1）：60－63.

［6］黃煒，楊蔚.高氣水比氣井井筒壓力的計算方法［J］.天然氣工業(yè)，2002，22（1）：64－67.

［7］張曉輝，張同義，劉鷹.利用氣井動態(tài)摩阻系數預測井口流出動態(tài)［J］.天然氣工業(yè)，2009，29（3）：80－83.

［8］BEGGS H D，BRILL J P.A study of two－phase flow in inclined pipes［J］.Journal of Petroleum Technology，1973，25（5）：607－617.

［9］NELSON P A.Rapid phase determinations in multiple phase flash calculation［J］.Computers ＆Chemical Engineering，1987，11（6）：581－591.

［10］趙智勇，劉道杰，蔣華，等.動態(tài)監(jiān)測指標在預判區(qū)塊開發(fā)中的應用研究［J］.西南石油大學學報：自然科學版，2009，31（1）：116－120.

［11］童景山.流體熱物理性質計算［M］.北京：清華大學出版社，1983.12.

［12］張允，姚軍，王子勝.測試井數據資料的數據清洗技術研究［J］.西南石油大學學報，2007，29（6）：27－30.

Pressure change law at the wellbore and well－adjacent zone of a HTHP or condensate gas well

Chang Zhiqiang，Xiao Xiangjiao，Zhu Zhongqian，Wang Xiaopei，Wang Haiying
（Exploration ＆Development Research Institute of Tarim Oilfield Company，PetroChina，Korla，Xinjiang 841000，China）

NATUR.GAS IND.VOLUME 32，ISSUE 10，pp.47－51，10／25／2012.（ISSN 1000－0976；In Chinese）

The distribution of shaft temperature and pressure of HTHP gas wells is significant for preparing and adjusting gas reservoir development plans.Therefore，based on the indoor experiments of wellbore friction simulation and temperature simulation in a high－temperature bore hole，the method of calculating PVT physical property parameters for HTHP natural gas is selected for theoretical research on the law of pressure and temperature distribution of the wellbore and well－adjacent zone of a HTHP or condensate gas well，which helps solve the problem in dynamic monitoring of such gas wells.A shaft pressure and temperature coupling calculation model is established considering such influencing factors as well fluid composition，fluid phase state，water cut，well bore structure，energy consumption，heavy hydrocarbon content，and ambient temperature.Moreover，the dynamic software is prepared for the wellbore of a HTHP or condensate gas well，which is capable of obtaining in real time the pressure of gas wells，the production change law as well as dynamic parameters of gas reservoirs.Field application indicates that this study provides technical approaches to promptly mastering the performance of gas fields and single wells，optimizing the plane structure of gas production，and preparing feasible development technologies.

Tarim Basin，HTHP gas well，condensate gas well，shaft，seepage，phase state，condensate，change，mathematical model

常志強等.高溫高壓氣井、凝析氣井井筒及近井地帶的壓力變化規(guī)律.天然氣工業(yè)，2012，32（10）：47－51.

10.3787／j.issn.1000－0976.2012.10.011

國家科技重大專項“塔里木盆地庫車前陸沖斷帶油氣開發(fā)示范工程”（編號：2011ZX05046）。

常志強，1980年生，高級工程師，博士；現從事油氣田開發(fā)研究工作。地址：（841000）新疆維吾爾自治區(qū)庫爾勒市石化大道26號。電話：（0996）2174787。E－mail：changzhq－tlm＠petrochina.com.cn

2012－06－21 編輯 韓曉渝）

DOI：10.3787／j.issn.1000－0976.2012.10.011

Chang Zhiqiang，senior engineer，born in 1980，holds a Ph.D degree and is currently engaged in development and research of oil ＆gas fields.

Add：No.26，Shihua Avenue，Korla，Xinjiang 841000，P.R.China

E－mail：changzhq－tlm＠petrochina.com