蔡 華 鹿克峰 何賢科 王 理
(中海石油(中國)有限公司上海分公司勘探開發(fā)研究院 上海 200335)
東海盆地西湖凹陷低滲氣主要分布在3 500 m以下的砂巖地層中,資源量占到總資源量的80%以上[1],但受制于地層出水問題尚未實(shí)現(xiàn)規(guī)模性開發(fā),具體表現(xiàn)為:西次凹、中央反轉(zhuǎn)構(gòu)造帶及裙邊3 800 m以下的低幅度構(gòu)造特低滲氣藏,自然測(cè)試基本無產(chǎn)出,加砂壓裂測(cè)試日產(chǎn)氣量為(1~2)×104m3,水氣比可達(dá)到10~300 m3/104m3;3 500~3 800 m深的低滲氣藏,投產(chǎn)起即表現(xiàn)為氣水同出狀態(tài),日產(chǎn)氣量在10×104m3左右,水氣比在1 m3/104m3左右,后續(xù)生產(chǎn)水氣比未見明顯下降;3 500 m以上的常規(guī)氣,投產(chǎn)初期基本不產(chǎn)水。初步定性判斷氣井投產(chǎn)初期水氣比的差異源自氣藏物性引起的氣水過渡帶差異,但要達(dá)到澄清問題、指導(dǎo)生產(chǎn)的目的,需要系統(tǒng)評(píng)價(jià)儲(chǔ)層物性與氣水過渡帶高度以及生產(chǎn)水氣比的定量關(guān)系。
現(xiàn)有文獻(xiàn)主要基于3類方法研究低滲、特低滲氣藏地層出水問題。方法1:基于特殊巖心分析的機(jī)理研究。借助微觀模型、離心毛管壓力、核磁共振、非穩(wěn)態(tài)氣驅(qū)水等組合實(shí)驗(yàn),明確水在孔隙中的賦存狀態(tài)及其可流動(dòng)性,認(rèn)為氣藏開發(fā)過程中巖石壓實(shí)作用[2]、天然氣降壓膨脹作用[3]、驅(qū)替壓差增大作用[4-5]會(huì)將部分易動(dòng)束縛水轉(zhuǎn)化為可動(dòng)水而產(chǎn)出。這類文獻(xiàn)從機(jī)理上解釋了“純氣藏”少量出水的原因,但對(duì)低滲氣井大量產(chǎn)水的情況討論較少。方法2:基于測(cè)井解釋與核磁共振實(shí)驗(yàn)的可動(dòng)水分布研究。首先采用巖心核磁共振實(shí)驗(yàn)確定束縛水飽和度,進(jìn)而統(tǒng)計(jì)巖心束縛水飽和度與巖心物性關(guān)系式,建立井點(diǎn)束縛水飽和度剖面,最終結(jié)合測(cè)井解釋的總含水飽和度剖面,實(shí)現(xiàn)井點(diǎn)可動(dòng)水分布的定量解釋[6-7]。這類文獻(xiàn)實(shí)現(xiàn)了井點(diǎn)可動(dòng)水定量表征,并通過生產(chǎn)井出水狀況建立可動(dòng)水飽和度與產(chǎn)水率的經(jīng)驗(yàn)統(tǒng)計(jì)關(guān)系,但經(jīng)驗(yàn)關(guān)系式存在普適性的問題。方法3:基于壓汞和非穩(wěn)態(tài)相滲實(shí)驗(yàn)的可動(dòng)水表征及產(chǎn)水率計(jì)算。采用進(jìn)汞毛管壓力曲線計(jì)算垂向含水飽和度分布;同時(shí)采用毛管壓力經(jīng)驗(yàn)截取值確定束縛水飽和度;二者結(jié)合確定可動(dòng)水垂向分布,最終采用非穩(wěn)態(tài)相對(duì)滲透率曲線計(jì)算不同飽和度條件下的含水率[8-10]。這類文獻(xiàn)實(shí)現(xiàn)了可動(dòng)水分布表征與產(chǎn)水率預(yù)測(cè),從思路和方法上值得借鑒,但應(yīng)用于西湖凹陷低滲氣藏還須做進(jìn)一步完善:①研究對(duì)象以油藏為主,適用于氣藏的定量評(píng)價(jià)方法有待建立;②研究目標(biāo)局限于具體氣藏,并未形成系統(tǒng)的、規(guī)律性認(rèn)識(shí);③采用壓汞曲線確定飽和度分布[11-13],采用非穩(wěn)態(tài)法相滲曲線計(jì)算產(chǎn)水率[14-16],與現(xiàn)行標(biāo)準(zhǔn)實(shí)驗(yàn)方法相比會(huì)存在一定偏差,實(shí)驗(yàn)方法有待進(jìn)一步優(yōu)選。
為揭示西湖凹陷不同物性氣藏出水差異的原因,本文以低滲氣藏為對(duì)象,以現(xiàn)行標(biāo)準(zhǔn)實(shí)驗(yàn)數(shù)據(jù)為基礎(chǔ),形成了定量的、系統(tǒng)評(píng)價(jià)初始產(chǎn)水狀況的思路和方法:①以半滲透隔板曲線代替?zhèn)鹘y(tǒng)的壓汞曲線,經(jīng)J函數(shù)處理將有限的實(shí)驗(yàn)數(shù)據(jù)拓展到任意物性條件,實(shí)現(xiàn)區(qū)域含水飽和度垂向表征;②以穩(wěn)態(tài)法相滲曲線代替?zhèn)鹘y(tǒng)的非穩(wěn)態(tài)法相滲曲線,經(jīng)實(shí)驗(yàn)曲線端點(diǎn)規(guī)律統(tǒng)計(jì)拓展到任意物性條件,實(shí)現(xiàn)區(qū)域氣、水相對(duì)滲透率預(yù)測(cè),進(jìn)而結(jié)合含水飽和度垂向表征,實(shí)現(xiàn)初始生產(chǎn)水氣比計(jì)算;③綜合不同物性、不同氣柱高度條件計(jì)算的含水飽和度、初始生產(chǎn)水氣比,形成定量的、系統(tǒng)的區(qū)域初始產(chǎn)水評(píng)價(jià)圖版,簡(jiǎn)單、快速評(píng)價(jià)低滲氣藏的出水狀況。
已有文獻(xiàn)指出[11-13],壓汞實(shí)驗(yàn)中汞并不是潤濕相,即使進(jìn)行了實(shí)驗(yàn)條件和油氣藏條件的轉(zhuǎn)換,轉(zhuǎn)換后的曲線并不能反映油氣藏條件下毛管壓力與含水飽和度的關(guān)系;而半滲透隔板法實(shí)驗(yàn)條件比較接近油氣藏條件,是檢驗(yàn)其他實(shí)驗(yàn)方法的標(biāo)準(zhǔn)測(cè)定方法。為了更精準(zhǔn)表征氣藏垂向飽和度的分布,實(shí)驗(yàn)選擇半滲透隔板法測(cè)定巖心氣水毛管壓力曲線,由美國某巖心公司負(fù)責(zé)實(shí)驗(yàn)測(cè)定,共計(jì)完成西湖凹陷16條氣水毛管壓力曲線(圖1a),總體表現(xiàn)為滲透率越低對(duì)應(yīng)的最終含水飽和度越高。
圖1 西湖凹陷半滲透隔板法實(shí)驗(yàn)毛管壓力曲線及J函數(shù)曲線Fig .1 Experimental capillary pressure curve and J function curve of porous diaphragm method in Xihu sag
如圖1a所示,各巖樣含水飽和度均隨毛細(xì)管壓力增大而降低。在實(shí)驗(yàn)達(dá)到最高毛細(xì)管壓力(0.83 MPa)時(shí),含水飽和度仍處于小幅降低中,由此可見,從實(shí)驗(yàn)數(shù)據(jù)中無法直接獲取臨界水飽和度數(shù)值。據(jù)調(diào)研,多次離心+核磁共振組合實(shí)驗(yàn)是目前確定臨界水飽和度的主流方法,但臨界水飽和度對(duì)應(yīng)的氣水毛細(xì)管壓力值在不同文獻(xiàn)中存在差異,典型值有0.69 MPa[17]、1.38 MPa[18-19]、2.87 MPa[20-21]。據(jù)2014年西湖凹陷多次離心+核磁共振組合實(shí)驗(yàn)測(cè)定結(jié)果,臨界水飽和度對(duì)應(yīng)的氣水毛細(xì)管壓力為2.07 MPa[22],遠(yuǎn)高于本次半滲透隔板法實(shí)驗(yàn)最高毛細(xì)管壓力,需要基于本次半滲透隔板法實(shí)驗(yàn)數(shù)據(jù)確定毛細(xì)管壓力等于2.07 MPa時(shí)的臨界水飽和度。采用二元回歸方法對(duì)實(shí)驗(yàn)數(shù)據(jù)進(jìn)行處理,得到毛管壓力為2.07 MPa下的臨界水飽和度相關(guān)關(guān)系式,即
Swc=-0.071 7ln(K/φ)+0.765 1
(R2=0.961 8)
(1)
對(duì)16塊巖樣測(cè)定的毛管壓力分別進(jìn)行J函數(shù)計(jì)算(計(jì)算中,氣藏氣水界面張力按照呂延防 等[23]建立的經(jīng)驗(yàn)關(guān)系式計(jì)算,在地層溫度140 ℃時(shí)計(jì)算氣水界面張力為25 mN/m;地層天然氣密度取160 kg/m3,地層水密度取1 000 kg/m3),相應(yīng)的含水飽和度分別進(jìn)行規(guī)格化處理,得到平均J函數(shù)曲線(圖1b)。
J(SW)=0.003SW-1.2
(2)
式(2)中:SW=(Sw-Swc)/(1-Swc)。
結(jié)合毛管壓力定義式、J函數(shù)定義式可建立氣柱高度與含水飽和度之間的關(guān)系式[24],即
-ρg)g]
(3)
將式(2)及相關(guān)參數(shù)代入式(3),得
(4)
式(4)即為依據(jù)有限的巖石實(shí)驗(yàn)數(shù)據(jù)經(jīng)J函數(shù)處理后,擴(kuò)展到區(qū)域含水飽和度垂向分布預(yù)測(cè)關(guān)系式。在已知?dú)獠乜紫抖群蜐B透率的情況下,即可結(jié)合式(4)和式(1)計(jì)算含水飽和度、可動(dòng)水飽和度的垂向變化。
不同含水飽和度條件下的氣、水兩相相對(duì)滲透率一般表達(dá)式[25]為
(5)
(6)
通過有限的實(shí)驗(yàn)數(shù)據(jù)獲取區(qū)域氣水相對(duì)滲透率預(yù)測(cè)關(guān)系式,需要確定出式(5)、(6)中包含的相對(duì)滲透率端點(diǎn)值[Krw(Sgr)、Krg(Swc)]、飽和度端點(diǎn)值(Swc、Sgr)、水相與氣相相對(duì)滲透率曲線常數(shù)(nw、ng)的統(tǒng)計(jì)規(guī)律。
用于氣藏開發(fā)過程的相對(duì)滲透率曲線有非穩(wěn)態(tài)法滲吸型、穩(wěn)態(tài)法滲吸型2種,其中穩(wěn)態(tài)法是相滲曲線測(cè)定的標(biāo)準(zhǔn)方法,但因?qū)嶒?yàn)時(shí)間較長,實(shí)際應(yīng)用較少。但與穩(wěn)態(tài)法相比,非穩(wěn)態(tài)法測(cè)取的氣相相對(duì)滲透率明顯偏高,殘余氣飽和度明顯偏低[14-16]。為了更準(zhǔn)確地反映低滲氣藏的滲流規(guī)律,選擇穩(wěn)態(tài)法,由中-加天然氣實(shí)驗(yàn)中心負(fù)責(zé)實(shí)驗(yàn)測(cè)定,完成了7條低滲巖樣穩(wěn)態(tài)法滲吸型相滲曲線的測(cè)定(圖2,表1)。
圖2 西湖凹陷穩(wěn)態(tài)法滲吸型氣水相對(duì)滲透率曲線Fig .2 Experimental gas-water relative permeability imbibition type curve of steadystate flow method in Xihu sag
表1 西湖凹陷穩(wěn)態(tài)法測(cè)定的氣水相對(duì)滲透率數(shù)據(jù)Table 1 Data of gas-water relative permeability measured by steady-state flow method in Xihu sag
1) 相對(duì)滲透率端點(diǎn)的求取。通過參數(shù)間的相關(guān)性分析,發(fā)現(xiàn)Krw(Sgr)與K存在較好的半對(duì)數(shù)關(guān)系,Krg(Swc)與Swc存在較好的冪函數(shù)關(guān)系,關(guān)系式分別為
Krw(Sgr)=0.034 6lnK+0.076 3 (R2=0.969 5)
(7)
(8)
2) 飽和度端點(diǎn)的求取。7塊巖樣的Sgr非常接近,介于0.296~0.324,平均值為0.307 6。
3) 水相與氣相相對(duì)滲透率曲線常數(shù)的求取。采用式(4)、(5)分別對(duì)7塊巖樣的實(shí)驗(yàn)數(shù)據(jù)進(jìn)行回歸,得到各巖樣的nw與ng(表1),多數(shù)數(shù)據(jù)點(diǎn)分布集中,取nw平均值為2.871 9,ng平均值為2.300 4。
將相對(duì)滲透率端點(diǎn)表達(dá)式(7)和(8)、飽和度端點(diǎn)平均值、相對(duì)滲透率曲線常數(shù)一并代入式(5)、(6),建立起區(qū)域氣藏氣水相對(duì)滲透率預(yù)測(cè)關(guān)系式如下:
Krw=(0.034 6lnK+0.076 3)×
(9)
(10)
依據(jù)式(9)、(10)并結(jié)合式(1),在已知?dú)獠乜紫抖群蜐B透率的情況下,即可計(jì)算任一孔隙度和滲透率的情況下氣、水相對(duì)滲透率與含水飽和度關(guān)系曲線(即氣水相對(duì)滲透率曲線)。
不考慮凝析水的情況下,氣藏中任一氣柱高度下生產(chǎn)水氣比由氣水兩相達(dá)西公式導(dǎo)出[26],即
WGR=(μgBgKrw)/(μwBwKrg)
(11)
由二項(xiàng)式產(chǎn)能方程確定氣井無阻流量[27],即
(12)
已知?dú)獠乜紫抖?、滲透率及自由水面,由式(1)、(4)、(9)、(10),可依次計(jì)算氣藏臨界水飽和度、原生水飽和度垂向分布、水相和氣相的相對(duì)滲透率垂向分布,進(jìn)而采用式(11)計(jì)算初始生產(chǎn)水氣比的垂向分布,采用式(12)計(jì)算氣藏產(chǎn)能的垂向分布。式(11)、(12)計(jì)算涉及的區(qū)域基本參數(shù)列于表2。
表2 西湖凹陷低滲氣藏初始生產(chǎn)水氣比與 產(chǎn)能計(jì)算基本參數(shù)Table 2 Basic parameters of production water-gas ratio and production capacity calculation of low permeability reservoirs in Xihu sag
取9個(gè)不同儲(chǔ)層物性,空氣滲透率分別為0.19、0.28、0.63、1.41、3.15、7.06、15.81、35.41、76.02 mD,將計(jì)算得到的不同物性條件下氣柱高度與含水飽和度的關(guān)系曲線、不同生產(chǎn)水氣比界限值對(duì)應(yīng)的氣柱高度與含水飽和度的關(guān)系曲線,無阻流量界限值對(duì)應(yīng)的氣柱高度與含水飽和度的關(guān)系曲線分別繪制在圖3、4中。
圖3 西湖凹陷不同儲(chǔ)層物性氣藏垂向產(chǎn)水區(qū)帶劃分圖Fig .3 Vertical water production zones division chart of different physical properties gas reservoirs in Xihu sag
圖4 西湖凹陷不同儲(chǔ)層物性氣藏垂向產(chǎn)能區(qū)帶劃分圖Fig .4 Vertical productivity zones division chart of different physical properties gas reservoirs in Xihu sag
圖3、4中黑色水平實(shí)線為以毛管壓力2.07 MPa為標(biāo)準(zhǔn)計(jì)算的純氣底界線,以西湖凹陷地層流體密度計(jì)算對(duì)應(yīng)的氣柱高度為251 m,可見假定成藏動(dòng)力相同時(shí)純氣底線只取決于流體密度,與儲(chǔ)層物性無關(guān)。
圖3、4中9條“S”形虛線為計(jì)算的9個(gè)不同滲透率氣藏氣柱高度與含水飽和度關(guān)系曲線,可見臨界水飽和度隨滲透率降低而增大;含水飽和度、可動(dòng)水飽和度(含水飽和度減去臨界水飽和度)隨氣柱高度減小而增大;距氣水界面較近時(shí),可動(dòng)水飽和度隨滲透率降低而減小,反之則隨滲透率降低而增大。
以西湖凹陷開發(fā)經(jīng)驗(yàn),在生產(chǎn)水氣比達(dá)到0.4 m3/104m3時(shí)(含有約0.2 m3/104m3左右的凝析水),氣井產(chǎn)能會(huì)出現(xiàn)快速降低;在生產(chǎn)水氣比達(dá)到10.0 m3/104m3時(shí),氣井將失去自噴能力。將2個(gè)重要的生產(chǎn)水氣比界限值繪制在圖3中(紅色實(shí)線為計(jì)算的水氣比為0.2 m3/104m3的等值線,藍(lán)色實(shí)線為生產(chǎn)水氣比為10.0 m3/104m3的等值線),連同純氣底界線(Sw=Swc)和自由水面(Sw=1),將不同儲(chǔ)層物性氣藏垂向上劃分為5個(gè)產(chǎn)水區(qū)帶(定義WGRi為初始生產(chǎn)水氣比,單位m3/104m3):
Ⅰ純氣區(qū):指Sw≤Swc,WGRi=0,只產(chǎn)氣不產(chǎn)水的氣藏區(qū)域;
Ⅱ近氣區(qū):指Sw>Swc,且0 Ⅲ氣水同產(chǎn)區(qū):指0.2 Ⅳ近水區(qū):指Sw<1,且WGRi>10,以產(chǎn)水為主,少量產(chǎn)氣的氣藏區(qū)域; Ⅴ純水區(qū),指Sw=1,只產(chǎn)水不產(chǎn)氣的氣藏區(qū)域。 表3列出了圖3中兩條水氣比等值線對(duì)應(yīng)的滲透率和氣柱高度值,在滲透率由76.02 mD降低至0.19 mD,水氣比為0.2 m3/104m3的等值線(氣水同產(chǎn)區(qū)頂界深度線)對(duì)應(yīng)的氣柱高度由2.35 m增加至251 m,水氣比為10 m3/104m3的等值線(氣水同產(chǎn)區(qū)底界深度線)對(duì)應(yīng)的氣柱高度由0.87 m增加至155.85 m,即隨著儲(chǔ)層滲透率降低,氣水同產(chǎn)區(qū)對(duì)應(yīng)的頂界深度和底界深度急劇上升,這揭示出了低滲氣藏易出水的根本原因。 表3 西湖凹陷低滲氣藏儲(chǔ)層滲透率與指標(biāo)界限 對(duì)應(yīng)的氣柱高度Table 3 Gas column height corresponding to permeability and index limit of low permeability gas reservoirs in Xihu sag 將圖3中水氣比為0.2 m3/104m3的等值線對(duì)應(yīng)的氣柱高度與滲透率建立相關(guān)關(guān)系,即可估計(jì)任意滲透率條件下近氣區(qū)底界深度 ΔH=59.169 6K-0.766 8(R2=0.994 1) (13) 要精確計(jì)算所需滲透率條件的氣柱高度(或改變水氣比界限值),可按本文方法和步驟完成。 圖4中綠色實(shí)線為計(jì)算的氣井無阻流量為6×104m3/d的等值線,將不同儲(chǔ)層物性氣藏垂向上劃分為2個(gè)產(chǎn)能區(qū)帶:自然產(chǎn)能區(qū)、非自然產(chǎn)能區(qū)。隨著滲透率降低,自然產(chǎn)能區(qū)所需的氣柱高度越大,這反映出含水飽和度對(duì)氣藏有效滲透率和產(chǎn)能的影響。 圖3、4中黑色實(shí)心點(diǎn)為西湖凹陷深層10個(gè)特低滲氣藏實(shí)際數(shù)據(jù)點(diǎn),自然測(cè)試基本無產(chǎn)出,加砂壓裂測(cè)試日產(chǎn)氣量為(1~2)×104m3,水氣比可達(dá)到10~300 m3/104m3,數(shù)據(jù)點(diǎn)落在近水區(qū)、氣水同產(chǎn)區(qū)底部及非自然產(chǎn)能區(qū);藍(lán)色實(shí)點(diǎn)為西湖凹陷中深層的低滲氣,投產(chǎn)起即表現(xiàn)為氣水同出狀態(tài),日產(chǎn)氣量在10×104m3左右,水氣比在1 m3/104m3左右,數(shù)據(jù)點(diǎn)落在氣水同產(chǎn)區(qū)及自然產(chǎn)能區(qū);紅色實(shí)點(diǎn)為西湖凹陷中—淺層常規(guī)氣,生產(chǎn)初期基本不產(chǎn)水,無阻流量(20~600)×104m3/d不等,數(shù)據(jù)點(diǎn)落在近氣區(qū)及自然產(chǎn)能區(qū)。由此可見,西湖凹陷各類氣井實(shí)際測(cè)試、生產(chǎn)情況與圖版預(yù)測(cè)結(jié)果吻合。 西湖凹陷低滲氣藏實(shí)際動(dòng)態(tài)表現(xiàn)和圖版相互吻合,揭示出西湖凹陷低滲氣井出水規(guī)律:儲(chǔ)層物性和氣柱高度決定了氣藏的出水狀況,有效產(chǎn)氣區(qū)(水氣比為10 m3/104m3的等值線)深度下限隨滲透率降低而急劇升高,盡管常規(guī)氣藏與特低滲氣藏同處于氣柱高度相近的氣水過渡區(qū)內(nèi),但常規(guī)氣藏以產(chǎn)氣為主,特低滲氣藏則以產(chǎn)水為主。 需要說明的是,本文中含水飽和度的計(jì)算是建立在氣水毛細(xì)管壓力平衡理論基礎(chǔ)上,適用于運(yùn)移成藏、存在氣水重力分異且靜止的氣藏中,從本文研究結(jié)果看,要實(shí)現(xiàn)該類低滲氣藏的有效開發(fā),需要相當(dāng)?shù)臍庵叨然蛐枰獙ふ矣行А疤瘘c(diǎn)”。在自生自儲(chǔ)或近源成藏條件下,由于不存在氣水界面,該類低滲氣藏含水飽和度可低于臨界水飽和度,形成毛管壓力欠平衡狀態(tài),是有利的開發(fā)目標(biāo),但本文提出的方法及建立的評(píng)價(jià)圖版均不適用。 1) 綜合不同物性、不同氣柱高度條件計(jì)算的含水飽和度、初始生產(chǎn)水氣比,形成了定量的、系統(tǒng)的區(qū)域初始產(chǎn)水評(píng)價(jià)圖版,西湖凹陷各類氣井實(shí)際測(cè)試、生產(chǎn)情況與圖版預(yù)測(cè)結(jié)果相吻合,驗(yàn)證了評(píng)價(jià)圖版的可靠性和實(shí)用性。 2) 儲(chǔ)層物性和氣柱高度決定了氣藏的出水狀況,有效產(chǎn)氣區(qū)(水氣比為10 m3/104m3的等值線)深度下限隨滲透率降低而急劇升高,盡管常規(guī)氣藏與特低滲氣藏同處于氣柱高度相近的氣水過渡區(qū)內(nèi),但常規(guī)氣藏以產(chǎn)氣為主,特低滲氣藏則以產(chǎn)水為主。 符號(hào)注釋 ΔH—?dú)庵叨龋琺; J(SW)—J函數(shù),無量綱; δgw—?dú)馑缑鎻埩Γ琺N/m; K—空氣滲透率,mD; φ—孔隙度,f; ρw—地層水密度,kg/m3; ρg—地層氣密度,kg/m3; SW—規(guī)格化含水飽和度,f; Sw—含水飽和度,f; Swc—臨界水飽和度,f; Krw—水相相對(duì)滲透率,f; Krw(Sgr)—?dú)堄鄽怙柡投认滤嘞鄬?duì)滲透率,f; Sgr—?dú)堄鄽怙柡投?,f; nw—水相相對(duì)滲透率曲線常數(shù),f; Krg—?dú)庀嘞鄬?duì)滲透率,f; Krg(Swc)—臨界水飽和度下氣相相對(duì)滲透率,f; ng—?dú)庀嘞鄬?duì)滲透率曲線常數(shù),f; WGR—生產(chǎn)水氣比,m3/m3; H—取樣深度,m; μg—地層氣相黏度,mPa·s; Bg—天然氣體積系數(shù),m3/m3; μw—地層水相黏度,mPa·s; Bw—地層水體積系數(shù),m3/m3; pR—地層壓力,MPa; pwf—井底流壓,MPa; qg—產(chǎn)氣量,m3/d; T—地層溫度,K; Z—天然氣偏差因子,f; re—泄氣半徑,m; rw—井腔半徑,m; S—污染表皮系數(shù),f; h—?dú)鈱雍穸?,m; D—紊流系數(shù); γg—天然氣相對(duì)密度,f; Qaof—?dú)饩疅o阻流量。4 結(jié)論