李 卓

（中國石化中原油田普光分公司，四川 達(dá)州 635000）

0 引言

異常高壓氣藏壓力系數(shù)大于1.3，在中國和世界分布廣泛，是極具開發(fā)價值的一種氣藏類型。國內(nèi)外資料統(tǒng)計(jì)表明，異常高壓氣藏普遍存在水體，60%以上的氣井在開采初期出現(xiàn)產(chǎn)水現(xiàn)象［1］。由于氣藏壓力系數(shù)高，氣藏的滲流規(guī)律和開發(fā)動態(tài)較正常壓力系數(shù)氣藏［2］有較大不同，再疊加水侵影響因素，使得氣藏開發(fā)困難、水侵預(yù)測難度大，特別是研究區(qū)長興組異常高壓氣藏礁體底水發(fā)育，氣水關(guān)系復(fù)雜，常規(guī)氣藏開發(fā)特征及水侵動態(tài)無法準(zhǔn)確地預(yù)測水侵動態(tài)，研究區(qū)缺少合理配產(chǎn)及高效開發(fā)方案［3］。因此，筆者綜合影響氣藏開發(fā)的各項(xiàng)因素，采用數(shù)值模擬方法研究水侵規(guī)律［4］，以期有針對性地制訂合理控水治水開發(fā)方案提供重要依據(jù)。

1 底水錐進(jìn)主控因素分析

1.1 見水時間影響因素分析

為了明確水侵規(guī)律，采用氣藏工程與數(shù)值模擬技術(shù)綜合考慮井型、采氣速度、儲層厚度、滲透率級差、水體倍數(shù)、打開程度等6項(xiàng)參數(shù)［5］，設(shè)計(jì)120套預(yù)測方案，研究異常高壓影響下底水錐進(jìn)的主控因素。

通過建立不同滲透率極差的底水模型，分別預(yù)測直井、水平井見水情況（圖1、圖2），選取井控儲量為35×108m3、滲透率為0.5×10-3μm2、采氣速度為5%，儲層厚度為60 m、距離底水距離為20 m的直井、水平井模型，在水體倍數(shù)分別為0.2、0.5、1、3、5時的錐進(jìn)速度［6］（表1）。

表1 影響水侵各因素的機(jī)理數(shù)模模型設(shè)計(jì)方案表

圖1 不同滲透率級差直井底水錐進(jìn)結(jié)果圖

圖2 不同滲透率級差水平井底水錐進(jìn)結(jié)果圖

1.2 見水時間主控因素分析

灰色關(guān)聯(lián)度分析法是將研究對象及影響因素的因子值視為一條線上的點(diǎn)，與待識別對象及影響因素的因子值所繪制的曲線進(jìn)行比較，比較它們之間的貼近度，并分別量化，計(jì)算出研究對象與待識別對象各影響因素之間的貼近程度的關(guān)聯(lián)度，通過比較各關(guān)聯(lián)度的大小來判斷待識別對象對研究對象的影響程度。

利用灰色關(guān)聯(lián)分析明確影響見水時間的主控因素：采用均質(zhì)化方法做無量綱處理后，進(jìn)行關(guān)聯(lián)度分析［7］（表2）。

表2 影響水侵各因素的機(jī)理數(shù)模模型設(shè)計(jì)方案表

因素間的關(guān)聯(lián)程度，主要是用關(guān)聯(lián)度的大小次序描述，而不僅是關(guān)聯(lián)度的大?。?］。將m個子序列對同一母序列的關(guān)聯(lián)度按大小順序排列起來，便組成了關(guān)聯(lián)序，記為｛x｝，它反映了對于母序列來說各子序列的“優(yōu)劣”關(guān)系。若，則稱｛x｝i對于同一母序列｛x0｝優(yōu)于｛xj｝，記為｛xi｝＞｛xj｝；表示第i個子序列對母數(shù)列特征值［9-10］。

通過對比直井、水平井不同水體規(guī)模的氣井水侵速度，可以很明顯地看到，水平井水侵速度受水體倍數(shù)變化敏感性強(qiáng)于直井。在水體倍數(shù)大于3倍時，水體大小對水侵的影響變?。▓D3）。

圖3 水平井、直井水體倍數(shù)與見水時間關(guān)系圖

通過對比直井、水平井不同儲層厚度的水侵速度，可以很明顯地看到，水平井水侵速度受儲層厚度的變化敏感性略強(qiáng)于直井，相同儲層厚度時，水平井見水時間晚于直井（圖4）。

圖4 水平井、直井儲層厚度與見水時間關(guān)系圖

通過對比直井、水平井不同滲透率級差的水侵速度，可以很明顯的看到，水平井水侵速度水平井水侵速度受滲透率級差的變化敏感性略強(qiáng)于直井，相同滲透率級差時，水平井見水時間晚于直井（圖5）。

圖5 水平井、直井滲透率極差與見水時間關(guān)系圖

通過對比直井、水平井不同儲量打開程度的水侵速度，可以很明顯的看到，隨著打開程度的增加，氣井越易見水；相同打開程度時，水平井見水時間長于直井，在打開程度大于70%時，水平井和直井見水時間接近；水平井對打開程度比直井敏感（圖6）。

圖6 不同打開程度直井、水平井見水時間關(guān)系圖

通過對比直井、水平井不同采氣速度的見水時間，可以很明顯的看到，隨著采氣速度的變化，見水時間的隨采氣速度的變化幅度基本接近，即水平井和直井的見水時間受采氣速度變化敏感性接近，在采氣速度大于4%時，采氣速度對水侵的影響變小（圖7）。

圖7 不同采氣速度直井、水平井見水時間關(guān)系圖

對于底水氣藏的直井來說，根據(jù)關(guān)聯(lián)系數(shù)求關(guān)聯(lián)度得：r4＞r3＞r5＞r2＞r1。

對于底水氣藏的水平井來說，根據(jù)關(guān)聯(lián)系數(shù)求關(guān)聯(lián)度得：r4＞r3＞r5＞r2＞r1。

因此，對于直井和水平井來說，儲層打開程度、滲透率級差均是底水錐進(jìn)的主控因素［11］。

對于直井，采氣速度對見水時間的影響大于儲層厚度；對于水平井，儲層厚度對見水時間的影響大于采氣速度［12］。

建立底水氣藏?cái)?shù)值模擬機(jī)理模型，分析水體倍數(shù)、打開程度、采氣速度等3個參數(shù)與底水推進(jìn)速度呈正相關(guān)；滲透率級差、儲層厚度兩個參數(shù)與底水推進(jìn)速度呈負(fù)相關(guān)。直井水侵呈“錐形”，水平井呈“凸形”，水平井控水效果優(yōu)于直井［13-14］。

2 異常高壓底水氣藏水侵預(yù)測模型研究

2.1 見水時間預(yù)測

根據(jù)前期研究成果，直井和水平井的見水時間與水體規(guī)模、儲層打開程度、采氣速度成反比，與滲透率級差、儲層厚度成正比。

令：

見水時間二元回歸公式為：

式中，H為儲層厚度，m；Kmax為最大滲透率，10-3μm2；Kavg為平均滲 透 率10-3μm2；Wpv為水體倍數(shù)；P為儲層打開程度，%；R為采氣速度，%；th為見水時間，a。

2.2 液氣比預(yù)測

為使多因素水侵預(yù)測間預(yù)測模型有較好的適應(yīng)性，將滲透率級差、采氣速度、水體倍數(shù)等水侵影響因素?zé)o因次化，在多因素?cái)?shù)值模擬實(shí)驗(yàn)基礎(chǔ)上，回歸分析建立多因素水侵模型［15］（表3）。

表3 參數(shù)變化表

不同滲透率級差下，液氣比與采出程度符合對數(shù)關(guān)系（圖8）：

圖8 液氣比隨采出程度的變化圖（Kmax／Kavg＝10）

式中，Kmax為最大滲透率，10-3μm2；Kavg為平均滲透率10-3μm2；Gp為 累積 采氣 量，108m3；G為 地質(zhì) 儲量，108m3；WGR為液氣比，103／104m3。

按以上方法，同樣可得到采氣速度、水體倍數(shù)、平均滲透率等單因素及系數(shù)A，B對應(yīng)的函數(shù)關(guān)系（圖9、10）。

圖9 系數(shù)A隨采氣速度的變化圖

圖10 系數(shù)B隨采氣速度的變化圖

采氣速度影響的液氣比方程：

水體倍數(shù)影響的液氣比方程：

式中，WGR為液氣比，103／104m2；Vg為卸載流速，m/s；Wpv為水體倍數(shù)；Gp為累積采氣量，108m3；G為地質(zhì)儲量，108m3。

針對不同因素分析得到的液氣比方程系數(shù)A，作多元回歸分析，確定滲透率級差、采氣速度、水體大小等影響因素的回歸系數(shù)，得到系數(shù)A的多因素方程表征。

式中，Kmax為最大滲透率，10-3μm3；Kavg為平均滲透率10-3μm2；Vg為卸載流速，m／s；Vwpv為過水體積倍數(shù)；H為儲層厚度，m；Hl為有效厚度，m。

此外，對液氣比方程系數(shù)B也作多元回歸分析，并確定了滲透率級差、采氣速度、水體大小等因素的回歸系數(shù)，得到系數(shù)B的多因素方程表征。

式中，Kmax為最大滲透率，10-3μm3；Kavg為平均滲透率，10-3μm2；Vg為卸載流速，m／s；Vwpv為過水體積倍數(shù)；H為儲層厚度，m；Hl為有效厚度，m。

考慮多因素影響的液氣比預(yù)測方程的建立及檢驗(yàn)。

直井液氣比定量預(yù)測模型：

式中，WGR為液氣比，103／104m2；Kmax為最大滲透率，10-3μm2；Kavg為平均滲透率10-3μm2；Vg為卸載流速，m／s；Gp為累積采氣量，108m3；G為地質(zhì)儲量，108m3；H為儲層厚度，m；Hl為有效厚度，m；Vwpv為過水體積倍數(shù)。

采用相同的方法建立水平井考慮多因素影響的液氣比預(yù)測方程：

水平井液氣比定量預(yù)測模型：

式中，WGR為液氣比，103／104m3；Kmax為最大滲透率，10-3μm2；Kavg為平均滲透率，10-3μm2；Vg為卸載流速，m／s；Gp為累積采氣量，108m3；G為地質(zhì)儲量，108m3；H為儲層厚度，m；Hl為有效厚度，m；Vwpv為過水體積倍數(shù)。

3 現(xiàn)場應(yīng)用情況

分別應(yīng)用見水時間預(yù)測模型和液氣比預(yù)測模型，預(yù)測異常高壓FX-1井見水時間及液氣比變化情況，結(jié)果反映了預(yù)測方法準(zhǔn)確性好。

3.1 見水時間預(yù)測

以異常高壓FX-1井參數(shù)進(jìn)行預(yù)測，預(yù)測見水時間為1.8年，實(shí)際見水時間1.6年，與氣井實(shí)際見水時間相比，誤差小于3個月。

直井見水時間與水體規(guī)模、儲層打開程度、采氣速度成反比，與滲透率級差、儲層厚度成正比。

將5個變量分為兩組，做二元線性回歸得到見水時間預(yù)測公式：

式中，th為見水時間，年。

根據(jù)氣藏各項(xiàng)參數(shù)，對FX-1井見水時間進(jìn)行預(yù)測，與氣井實(shí)際見水時間相比，平均誤差在三個月以內(nèi)。

3.2 液氣比預(yù)測

異常高壓FX-1井儲層滲透率級差約為11，水侵層厚度分別為13 m，與測井解釋高滲層厚度接近，液氣比預(yù)測模型與實(shí)測液氣比誤差在1.5%以內(nèi)（圖11），反映了筆者所研究的方法的準(zhǔn)確性較好。

直井液氣比定量預(yù)測模型：

4 結(jié)束語

通過對異常高壓底水氣藏水侵規(guī)律研究，得到下面3點(diǎn)結(jié)論：

1）水體倍數(shù)、儲層打開程度、采氣速度與見水時間呈正相關(guān)，當(dāng)水體倍數(shù)大于3.0倍、儲層打開程度大于70%、采氣速度大于4%，水侵影響的變化趨勢變緩。

2）滲透率級差、儲層厚度與底水推進(jìn)速度呈負(fù)相關(guān)，滲透率級差大于10水侵影響的變化趨勢變緩，儲層厚度越厚見水越晚。

3）基于異常高壓底水氣藏水侵單因素研究，建立見水時間預(yù)測公式和液氣比預(yù)測模型，見水時間預(yù)測誤差在3個月內(nèi)，液氣比預(yù)測相對誤差小于1.5%。