馮 曦 彭 先 李 騫 趙曉亮 張 平 潘 登

1.中國石油西南油氣田公司勘探開發(fā)研究院 2.中國石油大學(北京)石油工程學院3.中國石油川慶鉆探工程公司工程技術處 4.中國石油川慶鉆探工程公司鉆采工程技術研究院

0 引言

氣井產能評價是天然氣勘探開發(fā)研究的重要工作之一。其技術內涵主要包括：計算氣井絕對無阻流量，預測無阻流量隨開井時間增加而變化的規(guī)律，分析控制初始無阻流量及其變化特征的主要地質因素，確定氣井的合理配產。在測試資料或生產數(shù)據(jù)充足的條件下，依靠現(xiàn)有氣藏工程方法準確評價氣井產能并非難事。然而，隨著勘探開發(fā)一體化模式的普及，氣田開發(fā)評價及產能建設周期縮短，相關管理和研究工作者經常面臨動態(tài)資料不充分但又期望在早期能夠獲得準確認識的情況，在新井試氣工作中，這一需求尤為突出。為了避免在短時測試情況下嚴重誤判氣井的持續(xù)供給能力，對進一步完善氣井產能評價方法提出了新的要求。為此，筆者結合近年四川盆地氣井試氣分析的疑難問題，重點圍繞定壓生產試井模型的求解與計算，從試井分析理論層面研究了影響氣井不穩(wěn)定產能特征的主控因素，并且面向應用改進了分析預測氣井產能變化規(guī)律的方法，以期豐富中、高滲透氣藏氣井產能評價的技術手段，同時促進低滲透、強非均質性等復雜氣藏氣井產能早期評價問題的解決。

1 現(xiàn)有方法評價不穩(wěn)定產能的可行性

1.1 以穩(wěn)定滲流分析為主的方法

1935年Rawlins等[1]提出評價氣井產能的穩(wěn)定試井方法，1959年Katz等[2]奠定修正等時試井分析的基礎，這些方法一直沿用至今。針對礦場簡化分析的需求，以及特殊類型氣藏的復雜地質情況，后續(xù)研究又形成了一些補充方法[3-6]。該類方法在應用時至少需要一個穩(wěn)定生產制度的壓力、產量數(shù)據(jù)，評價結果代表氣井當時的“穩(wěn)定”產能。若儲層低滲、強非均質性導致氣井產能遞減顯著，方法適用性將變差。盡管可采用不同時間點穩(wěn)定滲流分析結果來表征氣井產能的變化[7]，但需要較多時間點的測試資料，在早期預測氣井產能不穩(wěn)定特征的能力相對較弱。

1.2 產量遞減分析方法

在理論核心層面，生產數(shù)據(jù)遞減分析方法與氣井無阻流量遞減分析方法有一定相似性，針對前者已有較多研究成果。1945年Arps[8]提出指數(shù)遞減、雙曲遞減及調和遞減分析模型，1980年Fetkovich[9]在前人研究定壓生產試井模型成果的基礎上形成了適用范圍更廣的產量遞減分析圖版，1991年Blasingame等[10]考慮產量、壓力同時遞減的普遍現(xiàn)象重新制作了分析圖版，1998年Agarwal等[11]針對有限導流和無限導流垂直裂縫井繪制了遞減分析圖版。該類方法主要用于評價生產井有效控制范圍和動用儲量，雖然在一定程度上也能預測氣井產能的變化，但需要較長時間的生產數(shù)據(jù)，不太適合早期預測。

1.3 不穩(wěn)定試井分析方法

不穩(wěn)定試井分析理論的核心是求解不穩(wěn)定滲流模型，雖然在此基礎上建立的現(xiàn)代產量遞減分析方法不適用于氣井產能的早期預測，但仍有可供借鑒的研究成果。

基于簡單試井模型尋求簡化計算公式的研究起源較早，1951年Horner[12]提出的分析方法屬于這類成果趨于成熟的典型代表。在解決復雜試井模型研究的關鍵問題方面，1949年Van Everdingen等[13]采用拉普拉斯變換方法求解不穩(wěn)定滲流模型，1973年Gringarten等[14]應用格林函數(shù)法改進上述模型解的計算，之后Rosa等[15]將Stehfest[16]研究形成的拉普拉斯數(shù)值反演算法引入到試井分析。在此基礎上，經過數(shù)十年發(fā)展，形成了現(xiàn)代試井分析的技術框架，以及眾多具有不同針對性的分析方法，其中，考慮在定產條件下壓力變化的分析方法占主導地位。

在氣井不穩(wěn)定產能評價方面，考慮定壓生產的試井模型具有獨特的優(yōu)點。雖然過去已開展了大量研究，形成了較多成果[17-18]，但在面向應用的早期預判、簡化計算和可靠性保障方面，仍存在疑難問題，這是長期以來定壓試井模型遠不如定產模型應用廣泛的主要原因。

筆者此次重點關注直井試氣分析，探索研究了滿足分析預測氣井不穩(wěn)定產能應用需求的改進方法。

2 技術關鍵點聚焦

2.1 不穩(wěn)定試井二項式產能方程應用及潛在問題

考慮均質儲層直井平面徑向不穩(wěn)定滲流，氣井擬壓力變化滿足以下關系式[19]：

式中Ψ表示擬壓力，MPa2/(mPags)；pi、pwf、psc分別表示開井前的地層壓力、開井時井底流動壓力、標準大氣壓，MPa；K表示儲層滲透率，mD；h表示儲層有效厚度，m；T、Tsc分別表示地層溫度、標準溫度，K；t表示開井時間，h；φ表示孔隙度；μ表示地層條件下天然氣黏度，mPags；Ct表示地層流體與巖石的綜合壓縮系數(shù)，MPa－1；rw表示井眼半徑，m；Sc表示表征井底污染或改善的凈表皮系數(shù)，無量綱；q表示氣產量，104m3/d；D表示紊流效應影響的非達西流動系數(shù)，(104m3/d)－1。公式系數(shù)對應于滲透率單位10－3μm2，在采用mD單位時近似相同。

式（1）為一次項系數(shù)隨時間變化的二項式產能方程，可用于計算不同時刻的無阻流量，前提是測試時間處于常規(guī)試井分析的徑向流直線特征階段，不滿足這一條件時使用式（1）計算無阻流量容易產生嚴重偏差，而二項式方程求根的計算方式往往使這種偏差難以被察覺。

2.2 影響常規(guī)分析方法適應性的敏感因素

借助定壓生產試井模型解的公式，能了解計算氣井無阻流量產生偏差的原因。前人已研究確認定壓、定產試井模型的長時間漸近解形式相同[20]，有

式中qAOF表示氣井無阻流量，104m3/d；S表示井底污染或改善以及紊流效應綜合影響狀況下的總表皮系數(shù)（以下簡稱表皮系數(shù)），無量綱。

在儲層滲透率低、表皮系數(shù)為負且絕對值較大的極端情況下，如果測試時間明顯小于徑向流直線段起點，則式（2）的計算結果可能為負數(shù)；若接近極端狀態(tài)，即使計算結果為正數(shù)，往往也存在較大誤差。

2.3 異常情況的根源及其影響程度

式（1）、（2）源自均質地層直井平面徑向滲流模型的長時間漸近解。考慮化簡前的指數(shù)積分表達式，采用有效井徑表征表皮效應，無阻流量計算公式為：

式中rwa表示井底污染或改善以及紊流效應綜合影響狀況下折算的井眼有效半徑，m。

式（3）比式（2）的適用范圍廣，已沒有負數(shù)問題，但x增大時，-Ei(-x)急劇變小[21]，如-Ei(-5)約等于0.001 15、-Ei(-10)約等于4.16h10－6，這導致在開井很短時刻計算的無阻流量異常偏高。另一方面，如果忽略紊流表皮效應隨產量變化的特點，將較低產量條件下的表皮系數(shù)用于式（3），也會使無阻流量計算結果偏大。

以四川盆地產能穩(wěn)定性較好、有可靠測試分析結果的4口高產氣井為例，說明應用式（3）產生的問題。LJ6、D4、LG1、MX10等井通過產能試井（穩(wěn)定試井或修正等時試井）分析得到的“穩(wěn)定”無阻流量依次為267.3h104m3/d、758.2h104m3/d、367.9h104m3/d、1 200.6h104m3/d，而采用試井分析商業(yè)軟件預測無阻流量的變化，初期計算值遠遠超過1 000h104m3/d（圖1），誤差很大。

圖1 氣井初期qAOF變化曲線圖（不穩(wěn)定滲流方法）

定壓生產試井模型理論研究早已成形，由于未完全解決上述異常問題，導致長期以來應用較少。

3 氣井絕對無阻流量遞減公式推導

3.1 定壓生產試井模型的無阻流量解析解

采用標準SI單位制進行理論分析，均質儲層邊界影響出現(xiàn)前氣體平面徑向滲流方程的通解為[22]：

式中r表示到井眼的距離，m；η表示地層導壓系數(shù)，m2/s；ξ表示公式推導的中間變量，m/s－0.5。

以井底流動壓力等于大氣壓為內邊界條件，求解得到井壁以外區(qū)域擬壓力的計算式為：

無阻流量的計算式則為：

式中Bg表示地層條件下天然氣體積系數(shù)，無量綱。

3.2 與經典試井理論長時間漸近解的差異

采用國內石油天然氣行業(yè)標準單位制，考慮有效井徑與表皮系數(shù)的關系，式（6）轉變?yōu)椋?/p>

式（7）未做近似性簡化處理，適用范圍比式（3）更廣。若開井時間較長，式（7）中指數(shù)項趨于1，過去研究長時間漸近解忽略該項，得到式（3）。在開井初期計算無阻流量時，這種做法容易引起較大計算誤差，由此產生如圖1所示的異常情況。

筆者根據(jù)試井模型精確解，結合指數(shù)積分高精度算法[23]，解決了計算公式僅適用于測試時間處于徑向流直線時間段的問題，增強了在早期評價氣井產能的適應性，并且在長時間條件下計算的準確性完全不受影響，能顯著抑制前述異常情況的出現(xiàn)。

4 無阻流量的控制因素分析

4.1 非線性影響關系圖版建立及特征分析

由式（7），定義非線性影響因子（σ）為：

根據(jù)式（8）計算、繪制得到σ與關系曲線圖版（圖2），由此，可以通過圖版查詢來確定σ，再結合式（7），用于計算、比較不同σ值情況下氣井產能的差異。

圖2 σ與關系曲線圖版

為了深入認識主要敏感參數(shù)對氣井產能的影響規(guī)律，拆分后再進行計算分析。設定導壓系數(shù)與時間的乘積（ηt）變化范圍介于（大致涵蓋滲透率介于0.01～1 000.00 mD的情況），表皮系數(shù)變化范圍介于－6～10，根據(jù)式（8）計算得到不同參數(shù)組合條件下的σ值（表 1）。為簡便起見，井眼半徑近似取值為0.1 m，該取值不影響對規(guī)律的認識。

如表1所示，在ηt小、負表皮系數(shù)絕對值大的情況下，計算不收斂，這間接說明了當井底存在大裂縫并且儲層滲透率很低時，短時間內難以達到平面徑向流供給狀態(tài)，并定量顯示出了平面徑向流試井模型完全不適用的范圍；同時，當ηt較?。ù篌w對應于滲透率較小、時間較短的情況）時，氣井產能對儲層改善程度的敏感性較強，而當ηt較大時該敏感性明顯降低。

表1 不同S、ηt下σ數(shù)據(jù)表

圖3 σ與S、ηt關系圖版

4.2 無阻流量計算中圖版的使用方法

筆者此次重點關注試氣評價，這種情況通常有壓力恢復測試數(shù)據(jù)，由此可以分析獲得滲透率、表皮系數(shù)，在此基礎上首先采用圖2或圖3確定所關注時間點對應的σ值，再用式（7）計算氣井無阻流量。在其他的應用情形下，對基礎數(shù)據(jù)的需求和分析方法與之類似。

根據(jù)式（7）計算qAOF需要估算在井底壓力為大氣壓的理論假設條件下，由于產量較高導致紊流效應影響嚴重時的S，這是傳統(tǒng)試井分析未涉及的問題，主要解決辦法包含以下兩種：①當有多個產量制度的試井資料時，分析不同產量對應的S，根據(jù)式（2）包含的S—Sc關系式，采用直線回歸法確定Sc和D；②當已知某一時刻qAOF時，按式（7）計算與之對應的S，結合試井測試時的產量和由試井解釋得到的S，確定Sc和D。若有產能發(fā)揮程度較高條件下的生產壓差測試數(shù)據(jù)，可以省去第②種方法中無阻流量、試井解釋得到的表皮系數(shù)兩個條件之一。

在獲得Sc和D之后，對任意時間點，先估算qAOF初值，根據(jù)S—Sc關系式計算與qAOF對應的S，再按式（7）重新計算qAOF，當相鄰兩次qAOF計算結果的差值超出允許誤差范圍時，以最近一次qAOF計算結果為新起點重復上述過程進行迭代計算，直至最終獲得相應時刻無阻流量的準確值。

4.3 井底裂縫的影響

在裂縫較長的情況下，平面徑向流有效井徑等效模型（簡稱平面徑向流等效模型）與無限導流垂直裂縫井模型（簡稱垂直裂縫井模型）相似性較低，從嚴格意義上講文中方法不適用。但是，若氣井產能較高，因大產量條件下表皮系數(shù)較大，這時即使凈表皮系數(shù)為負且絕對值較大，新方法也能算出結果，是否近似可用則根據(jù)具體情況來定。

以定產量生產試井模型來進行計算分析，結果顯示在實際工作常見的裂縫長度、儲層滲透率數(shù)值區(qū)間，兩種模型計算結果最大差距不超過15%（圖4），對于氣井產能評價而言尚可接受。裂縫越長、儲層滲透率越低，平面徑向流等效模型在長時間情況下的計算誤差越大，這時應注意避免誤用。

圖4 平面徑向流等效模型與垂直裂縫井模型計算生產壓差之比變化曲線圖

4.4 井筒儲集效應的影響

氣井井筒儲集效應顯著，前面的試井理論公式推導未考慮該因素，這里增加相關分析。

井筒儲集系數(shù)對產量的影響式為[24]：

式中qwh表示井口產量，104m3/d；qwb表示井底產量，104m3/d；pw表示井底壓力，MPa；pwh表示井口壓力，MPa；C表示井筒儲集系數(shù)，m3/MPa；表示井筒平均壓力、平均溫度條件下的天然氣體積系數(shù)，無量綱。

在定井底壓力生產的條件下，井筒儲集效應影響時間段結束之后，井底與井口壓力差值近似不變，此時式（9）導數(shù)項為零。這表明只要不涉及開井早期井筒儲集效應影響階段的分析，即使理論模型未考慮該因素，對產能評價的影響也不大。

5 復雜情況下無阻流量的變化特征

5.1 定量分析與定性預判搭配的策略

筆者基于均質無邊界地層直井開展研究，形成了快速簡便的分析方法，可手工預測氣井不穩(wěn)定產能。針對大型壓裂井、水平井、雙重介質及復合儲層、封閉邊界等復雜情況，前述新方法不適用。相關研究涉及拉普拉斯變換求解、無窮級數(shù)或積分等計算[25-26]，雖然這方面早已有較多成果[27-29]，并且近年有一些新的探索研究取得成功[30-31]，但基本上都高度依賴計算機、不太適合直觀理解與簡便分析。

早期預測氣井產能，并分析其不穩(wěn)定特征屬長期性疑難問題，諸多研究有待深入。在儲層非均質性不強、直井試氣的簡單情況下可以直接采用前述新方法，在復雜情況下不具備精細定量分析條件時，用新方法也可以提供定性的類比參考，這是現(xiàn)階段滿足實際工作需求的較好方式。

5.2 針對應用關注點的理論分析

以均質儲層平面徑向流模型為參照，針對大型壓裂井、裂縫—孔隙型儲層、遠井區(qū)物性變差的典型情況另做計算研究。將式（7）等號右邊與σ無關的部分移至左邊，然后定義無量綱無阻流量，再將計算結果進行對比，從而深化對不同類型氣井不穩(wěn)定產能特征存在差異的認識。這里側重于研究早期階段氣井的產能特征，沒有考慮供給邊界的影響，圖5中設定的較長時間跨度主要用于適配不同滲透率條件下無量綱時間分布區(qū)間的不同，而不是為了研究超長時間階段氣井產能的變化規(guī)律。

圖5 典型井初期無阻流量遞減曲線圖

5.3 主要規(guī)律性認識

如圖5所示，當無量綱時間值較小時，負表皮或穿過井底的大裂縫將使氣井初期產能大幅度升高，但早期遞減快；對于裂縫—孔隙型儲層，當裂縫滲透率高但與孔隙搭配關系不好、基質孔隙系統(tǒng)供給嚴重滯后時，初期產能主要受裂縫滲透率影響，在裂縫主導供給的早期階段結束后，將出現(xiàn)產能遞減較快的階段；若氣藏宏觀非均質性較強、遠井區(qū)儲層物性比近井區(qū)差，當生產井壓力擾動完全傳播到遠井區(qū)之后，將在一定時間階段表現(xiàn)出產能遞減加快的現(xiàn)象。

以上3種情況屬于氣井初期產能與長期穩(wěn)產能力差異較大的典型情況。以理論研究獲得的認識為指導，在實際工作中針對性地強化地質與動態(tài)分析研究，研判是否屬于上述典型情況，能在短時測試條件下有效減少對氣井產能不穩(wěn)定特征的誤判。

6 新方法的應用

6.1 氣井不穩(wěn)定產能評價的簡化參照時間點

不同類型氣井無阻流量變化特征差異較大，若未采用統(tǒng)一的時間基準，則不同類型氣井產能評價結果的可比性較差，不利于后續(xù)研究與應用。四川盆地氣井試氣多在白天開井測試，連續(xù)開井測試一般不超過10 h，在大多數(shù)情況下試氣期間關井壓力恢復測試時間接近或超過3 d。結合實際情況，建議以10 h、3 d、30 d為參照時間點，分別評價開井初期、相對穩(wěn)定狀態(tài)及較長時間狀態(tài)下的無阻流量，依次對應于根據(jù)開井實測資料直接分析或適度外推預測、經關井壓力恢復實測數(shù)據(jù)驗證的試井模型預測、超出實測數(shù)據(jù)檢驗范圍具有一定不確定性的長時間預測結果。這樣做既減少分析工作量，又便于統(tǒng)一對比不同類型氣井產能特征的差異、強化整體認識。

6.2 新方法分析結果的含義及其與類似方法的差異

新方法預測絕對無阻（井底壓力為大氣壓的理論假設）條件下產量變化規(guī)律，偏重于揭示不受人為因素影響的本質特征，唯一性較強。過去經常采用變通性處理方法，即根據(jù)氣井實際生產條件下的壓力變化趨勢預測后續(xù)井底流動壓力，再采用產能方程或“一點法”計算相應時刻的無阻流量，側重于評價特定產量條件下無阻流量的變化規(guī)律。在儲層滲透率低或單井供給范圍小的情況下，其分析結果會因配產不同而明顯不一致。

與上述變通性處理方法類似，采用定產量生產試井模型預測井底流動壓力的變化，可以分析多個產量條件下后續(xù)時刻的無阻流量，但是，由于較難估算與虛擬產量對應的表皮系數(shù)，該預測方法的應用受限。

6.3 新方法在GS1井應用

GS1井是四川盆地某特大型碳酸鹽巖氣藏的發(fā)現(xiàn)井，該井于2011年7月在震旦系燈二段試氣獲工業(yè)氣流，測試氣產量為102.15h104m3/d，測試過程中以穩(wěn)定產量開井2.4 h，之后關井壓力恢復48 h。與新方法應用相關的GS1井基礎數(shù)據(jù)、試氣測試與分析數(shù)據(jù)如表2所示。

表2 GS1井數(shù)據(jù)表

GS1井試氣解釋儲層滲透率為0.27 mD；井底存在裂縫發(fā)育區(qū)，等效于半長為27.7 m的無限導流垂直裂縫；轉換為有效井徑近似處理模型，計算凈表皮系數(shù)為－6.078。從圖5知這屬于初期產能高但遞減較快的類型。因此應重點關注該井產能不穩(wěn)定特征。

GS1井試氣時生產壓差達到地層壓力的20%，屬于評價氣井產能時采用常規(guī)一點法（經驗系數(shù)α等于0.25）適用性較好的范疇[3,32]，相關研究表明在大生產壓差條件下常規(guī)一點法計算無阻流量的誤差不超過10%，采用該方法計算得到GS1井在開井2.4 h時對應的“穩(wěn)定”無阻流量為187.77h104m3/d。

根據(jù)上述短時開井時刻的qAOF值，用式（7）計算相應條件下的S，其迭代計算方法見下文。該井測試產量達到無阻流量的54.4%，紊流效應的影響反應較充分，根據(jù)測試產量及其相應的生產壓差，用式（7）迭代計算確定與之對應的S，也能用于后續(xù)分析。通過試井分析獲得了Sc，計算了初期qAOF，實測q與qAOF的比值較高，滿足這三個條件之二，即可計算出非達西流動系數(shù)（表3）。

表3 新方法計算GS1井S及D統(tǒng)計表

按前述10 h、3 d、30 d時間點，取表3中第1種方式對應的D值，以一點法計算的“穩(wěn)定”qAOF為初值，用式（7）估算絕對無阻條件下的S，在此基礎上由圖2或圖3查詢得到σ并代入式（7），重新計算相應時刻的qAOF，迭代計算過程如表4所示。

表4 預測GS1井qAOF變化特征的迭代計算數(shù)據(jù)表

筆者采用制作圖版的原始數(shù)據(jù)手工計算，當qAOF穩(wěn)定至整數(shù)部分的個位（四舍五入后）不變時，終止迭代；若采用查圖版手工計算，可將S穩(wěn)定至小數(shù)第一位（四舍五入后）不變作為迭代計算的終止條件。從表4數(shù)據(jù)看，時間越長迭代收斂越快；按上述第二種終止條件，查圖版迭代計算大體不超過5次即可獲得最終結果，該條件下qAOF計算值與按第一種迭代終止條件所計算的qAOF相比誤差小于5%。

GS1井（燈二段為主產層）2012年9月投產后，配產介于5h104～10h104m3/d時生產較穩(wěn)定；在試采過程中開展過氣產量為35h104m3/d生產的試驗，6天內氣產量遞減到12h104m3/d。該井近井區(qū)存在裂縫使得氣井初期產能較高，遠井區(qū)滲透率低則引起產能快速下降。表4中的最終計算結果簡要描述了該井無阻流量的遞減規(guī)律，可以看出，即使是試氣期高產這種容易使人們對氣井持續(xù)穩(wěn)產能力產生經驗性誤判的情況，新方法仍然能準確預測產能不穩(wěn)定本質特征，并且預測結果與后續(xù)生產情況相符。

6.4 新方法的其他應用示例

下面以近十余年四川盆地天然氣開發(fā)熱點區(qū)塊的部分重點井為例，來驗證新方法的應用效果。

圖1中涉及的4口氣井均是探井，且做過專項試井，有高質量壓力恢復試井資料，以及穩(wěn)定試井或修正等時試井資料。試井解釋得到了對應于短時開井時刻的“穩(wěn)定”qAOF，以及多個測試產量分別對應的S，由此計算了D，結果見表5。采用文中介紹的新方法另外計算D，并推算后續(xù)時刻的qAOF，表6定量揭示了穩(wěn)產能力較強的氣井不易被發(fā)現(xiàn)的產能遞減特征及存在的差異，有助于深化氣井產能評價的認識。

表5 四川盆地4口高產探井產能試井計算qAOF和D統(tǒng)計表

表6 四川盆地4口高產探井不同開井時間點qAOF預測數(shù)據(jù)表（新方法）

新方法為獲取紊流效應影響下的非達西流動系數(shù)提供了新途徑，在此基礎上對比不同氣井的非達西流動系數(shù)，能間接認識儲層微觀儲集空間特征存在的差異。表5、6中數(shù)據(jù)顯示新方法分析的非達西流動系數(shù)比試井解釋結果略大，且產能越大的氣井這種差異越明顯?；诶碚摲治?，結合部分氣井產量增大后表皮系數(shù)與產量線性關系發(fā)生變化的實際情況，認為在遠低于無阻流量的開井產量狀態(tài)下試井，解釋得到的非達西流動系數(shù)偏小。

采用新方法預測氣井無阻流量的變化需要儲層孔隙度和滲透率、流體黏度、綜合壓縮系數(shù)以及初期無阻流量參數(shù)值；另外需要凈表皮系數(shù)、非達西流動系數(shù)二者之一。研究發(fā)現(xiàn)對一些高產氣井而言，因在絕對無阻條件下紊流效應對其無阻流量的影響遠遠大于凈表皮系數(shù)的影響。因此即使凈表皮系數(shù)未知，無阻流量預測也相當準確（表7）。

表7 Sc誤差較大情況下高產氣井qAOF預測數(shù)據(jù)表

大多數(shù)新井試氣時都未進行產能試井，應用新方法時通常需要依靠“一點法”提供參照基準。對于生產壓差較小的氣井，常規(guī)“一點法”誤差較大，宜采用參數(shù)約束改進“一點法”[32]。

以上列舉了測試資料較全的高產探井典型應用案例，另有部分探井及上試新層老井的分析實例如表8所示。新方法與傳統(tǒng)方法配合，能夠顯著提高氣井產能評價的質量，獲得更豐富的認識，為天然氣勘探與開發(fā)評價、產能建設可行性論證及氣藏開發(fā)設計提供重要技術依據(jù)。

7 結論

1）獲得均質儲層平面徑向流模型氣井無阻流量隨時間變化精確解，其適用范圍拓展至徑向流直線段出現(xiàn)之前。

2）確認表皮系數(shù)和地層導壓系數(shù)是氣井不穩(wěn)定產能特征多樣化的主要控制因素，分析揭示了其影響關系。

表8 四川盆地部分探井及上試新層氣井qAOF預測數(shù)據(jù)表

3）形成氣井不穩(wěn)定產能非線性影響因子分析圖版、無阻流量遞減特征預測的迭代算法，達到實用化程度。

4）提出計算非達西流動系數(shù)新方法，既可用于分析氣井不穩(wěn)定產能，也可用于間接評價儲層微觀儲集特征。

5）對于均質、視均質或與之相似儲層中的直井，新方法分析準確，具有較強實用性。