周敏，王涵，張娜，龐宇來(lái)，陳黎，歐陽(yáng)沐鯤，張偉

1.中國(guó)石油西南油氣田分公司重慶氣礦地質(zhì)研究所，重慶 400021 2.重慶天然氣儲(chǔ)運(yùn)有限公司，重慶 400021

氣井的絕對(duì)無(wú)阻流量定義為井底壓力等于大氣壓力時(shí)的產(chǎn)氣量，它反映了氣井生產(chǎn)的一種最大能力。因此，獲取氣井絕對(duì)無(wú)阻流量都是氣井產(chǎn)能分析最重要的任務(wù)之一。目前，對(duì)已投產(chǎn)井的產(chǎn)能評(píng)價(jià)方法很多，其中穩(wěn)定試井和“一點(diǎn)法”產(chǎn)能試井是獲取氣井絕對(duì)無(wú)阻流量較為可靠的方法[1-4]。兩種方法均需要?dú)饩懂a(chǎn)，并獲得較為準(zhǔn)確的測(cè)試資料，因此，如何在氣井鉆井之前對(duì)氣井絕對(duì)無(wú)阻流量進(jìn)行預(yù)測(cè)是目前產(chǎn)能評(píng)價(jià)的最大的問(wèn)題[5,6]。筆者以四川盆地川東地區(qū)碳酸鹽巖氣田為研究目標(biāo)，在系統(tǒng)分析川東地區(qū)已投產(chǎn)井絕對(duì)無(wú)阻流量與氣井各主要?jiǎng)?、靜態(tài)參數(shù)之間關(guān)系的基礎(chǔ)上，應(yīng)用回歸分析研究氣井絕對(duì)無(wú)阻流量與其他變量之間的關(guān)系，使用氣井在未鉆井之前容易獲得的參數(shù)預(yù)測(cè)建立了氣井絕對(duì)無(wú)阻流量的方程，以指導(dǎo)開(kāi)發(fā)井的鉆井決策。

1 地質(zhì)背景

研究樣本均來(lái)自于川東地區(qū)碳酸鹽巖氣田實(shí)際生產(chǎn)測(cè)試數(shù)據(jù)。川東地區(qū)碳酸鹽巖氣藏開(kāi)發(fā)始于20世紀(jì)30年代石油溝氣藏，目前川東地區(qū)在縱向上已發(fā)現(xiàn)十余個(gè)含油氣層系。川東氣田主力生產(chǎn)層系為石炭系氣藏，長(zhǎng)興組氣藏是主要的接替氣藏，研究的樣本多來(lái)自這2個(gè)生產(chǎn)層系。

1.1 主要地質(zhì)特征

川東地區(qū)石炭系沉積相模式為相對(duì)局限的海灣向正常開(kāi)闊的淺演化規(guī)律和潮汐作用，主要產(chǎn)層石炭系和礁灘在川東地區(qū)廣泛分布。

川東地層石炭系氣藏普遍為受構(gòu)造控制的，具有邊水的整裝氣藏，部分氣藏為巖性和構(gòu)造控制的巖性-構(gòu)造復(fù)合氣藏，還有部分氣藏具有高陡構(gòu)造特征。礁灘氣藏主要受巖性控制，氣藏氣水關(guān)系較復(fù)雜，但生產(chǎn)中地層水對(duì)生產(chǎn)影響較小。

川東地區(qū)石炭系儲(chǔ)層主要分布于黃龍組二段。儲(chǔ)集巖主要為白云巖，儲(chǔ)集空間總體可分為孔隙、洞穴、裂縫3大類(lèi)。儲(chǔ)層孔隙形態(tài)極不規(guī)則，孔喉體系分選差，歪度細(xì)，峰態(tài)平緩，連通性較好，多具有中孔小喉特征。儲(chǔ)層總體上屬于低孔低滲儲(chǔ)層，橫向分布較為連續(xù)，局部裂縫較發(fā)育，儲(chǔ)層非均質(zhì)性強(qiáng)。根據(jù)川東地區(qū)多年來(lái)的勘探開(kāi)發(fā)實(shí)踐石炭系主要是裂縫-孔隙型。

川東地區(qū)另一主產(chǎn)層長(zhǎng)興組氣藏儲(chǔ)集巖主要也為白云巖。儲(chǔ)集空間集孔、洞、縫三類(lèi)于一體，起主要儲(chǔ)集作用的是各類(lèi)孔隙和溶洞。整體上表現(xiàn)為中低孔、中低滲的裂縫-孔隙型儲(chǔ)層特征，局部存在高孔高滲儲(chǔ)層。

1.2 主要生產(chǎn)動(dòng)態(tài)特征

1)低滲透氣藏生產(chǎn)特征：①多數(shù)氣井原始絕對(duì)無(wú)阻流量小于20×104m3/d，單井穩(wěn)定產(chǎn)量大都在3×104m3/d以下；②多數(shù)氣井生產(chǎn)壓差大于10MPa，要求對(duì)氣井進(jìn)行精心配產(chǎn)；③多數(shù)氣井開(kāi)采初期壓力下降較快，但中后期低壓下，壓力下降緩慢，生產(chǎn)時(shí)間很長(zhǎng)；④氣藏連通關(guān)系到開(kāi)采中后期才能搞清楚；⑤關(guān)井壓力表現(xiàn)出恢復(fù)緩慢的特點(diǎn)，往往1個(gè)月以上；⑥多數(shù)氣井壓降儲(chǔ)量小于5×108m3。

2)中高滲氣藏生產(chǎn)特征：①多數(shù)氣井原始絕對(duì)無(wú)阻流量大于20×104m3/d，單井穩(wěn)定產(chǎn)量大都在3×104m3/d以上；②多數(shù)氣井生產(chǎn)壓差小于3MPa，調(diào)整配產(chǎn)對(duì)壓力影響較??；③多數(shù)氣井開(kāi)采初期壓力下降較慢，氣藏穩(wěn)產(chǎn)時(shí)間長(zhǎng)，穩(wěn)產(chǎn)期采出程度高；④氣藏連通關(guān)系在開(kāi)采早期就能搞清楚，開(kāi)采中各井壓力下降均衡，氣藏連通性好，生產(chǎn)中一般沒(méi)有大的壓降漏斗；⑤關(guān)井壓力表現(xiàn)出恢復(fù)很快的特點(diǎn)，許多井1天就能夠恢復(fù)；⑥多數(shù)氣井壓降儲(chǔ)量大于5×108m3，部分井壓降儲(chǔ)量甚至超過(guò)30×108m3。

從生產(chǎn)特征可以知道，在川東地區(qū)絕對(duì)無(wú)阻流量越大，氣井越能在較高產(chǎn)量下穩(wěn)定生產(chǎn)，同時(shí)單井動(dòng)態(tài)儲(chǔ)量也越大。絕對(duì)無(wú)阻流量幾乎成為后期開(kāi)發(fā)井是否實(shí)施的關(guān)鍵決策指標(biāo)。

2 影響絕對(duì)無(wú)阻流量的主要因素

假設(shè)氣藏?zé)o窮大，當(dāng)測(cè)試井生產(chǎn)所引起的壓力變化未波及到邊界時(shí)，氣井產(chǎn)能的二項(xiàng)式方程為[1,2,7]：

(1)

其中：

(2)

(3)

(4)

根據(jù)二項(xiàng)式產(chǎn)能方程式(1)，使用絕對(duì)壓力pwf=0.101325MPa時(shí)，氣井絕對(duì)無(wú)阻流量計(jì)算式：

(5)

從式(5)可以看出，影響氣井絕對(duì)無(wú)阻流量是A、B和pe，且絕對(duì)無(wú)阻流量Qaof是隨A的增大而減小的[7]。根據(jù)式(2)影響A的主要是地層系數(shù)(kh)，此外天然氣黏度、氣藏溫度及表皮系數(shù)對(duì)氣井的產(chǎn)能也具有一定的影響[8-10]，但同一區(qū)塊天然氣黏度、氣藏溫度差異較小，且對(duì)于新鉆井來(lái)說(shuō)，表皮系數(shù)一般較小，因此天然氣黏度、氣藏溫度及表皮系數(shù)對(duì)氣井產(chǎn)能影響可以忽略不計(jì)。此外，鄭超等[11]對(duì)川東石炭系氣藏氣井有效孔隙度及儲(chǔ)層厚度對(duì)產(chǎn)能的影響進(jìn)行了分析研究，認(rèn)為隨著儲(chǔ)層厚度及孔隙度降低，氣井產(chǎn)能降低。

根據(jù)上述理論分析，結(jié)合川東地區(qū)天然氣生產(chǎn)實(shí)際數(shù)據(jù)，筆者對(duì)川東地區(qū)現(xiàn)場(chǎng)穩(wěn)定試井獲得的絕對(duì)無(wú)阻流量與主要?jiǎng)屿o態(tài)參數(shù)(kh、pe、φ)和動(dòng)態(tài)儲(chǔ)量相關(guān)性進(jìn)行了分析。結(jié)果表明，絕對(duì)無(wú)阻流量與kh、pe、φ和動(dòng)態(tài)儲(chǔ)量均有一定的正相關(guān)性(見(jiàn)圖1)。在天然氣生產(chǎn)中，對(duì)于已開(kāi)發(fā)氣藏的未鉆開(kāi)發(fā)井更容易獲得較可靠的kh、pe等動(dòng)靜態(tài)參數(shù)，而動(dòng)態(tài)儲(chǔ)量在鉆井前無(wú)法獲得，因此研究這些參數(shù)與絕對(duì)無(wú)阻流量定量關(guān)系就能夠獲得開(kāi)發(fā)井的絕對(duì)無(wú)阻流量預(yù)測(cè)方程。

3 新鉆開(kāi)發(fā)井絕對(duì)無(wú)阻流量預(yù)測(cè)方程的建立

3.1 基本思想

回歸分析方法是常用的數(shù)理統(tǒng)計(jì)方法，主要用于研究因變量y與多個(gè)自變量x1、x2,…,xn之間的數(shù)量依存關(guān)系，并利用這種統(tǒng)計(jì)關(guān)系在一定置信度下由各因素的取值去預(yù)測(cè)響應(yīng)值，根據(jù)預(yù)報(bào)變量的給定值來(lái)估計(jì)和預(yù)測(cè)精度，在生產(chǎn)實(shí)踐的各個(gè)方面都具有廣泛的應(yīng)用[12-16]。

影響氣井絕對(duì)無(wú)阻流量的因素非常多，通過(guò)回歸分析方法來(lái)獲取少量有關(guān)鍵意義的因素。其基本思想是將變量一個(gè)一個(gè)引入，并對(duì)引入的變量逐個(gè)進(jìn)行顯著性檢驗(yàn)，當(dāng)原引入的變量由于后面變量的引入而變得不再顯著時(shí)，要將其剔除。整個(gè)過(guò)程反復(fù)進(jìn)行，直到既無(wú)顯著的自變量選入，也無(wú)不顯著自變量從回歸方程中剔除為止，這樣就保證了最后所得的回歸子集是“最優(yōu)”回歸子集[17-20]。

3.2 預(yù)測(cè)方程

1)樣本井的選取。根據(jù)川東地區(qū)的地質(zhì)特征，川東地區(qū)主要產(chǎn)層以裂縫-孔隙型氣藏為主，裂縫-孔隙型氣藏投產(chǎn)井占?xì)馓锶客懂a(chǎn)井的90%以上。系統(tǒng)分析川東氣田裂縫-孔隙型氣藏已投產(chǎn)井資料錄取情況，隨機(jī)從資料錄取較全、數(shù)據(jù)解釋可靠的氣井中選取30口氣井作為分析樣本[21，22]。筆者選取地層系數(shù)(kh)、地層壓力(pe)、孔隙度(φ)作為多元線性回歸分析的主要參數(shù)(見(jiàn)表1)。

2)絕對(duì)無(wú)阻流量預(yù)測(cè)的線性回歸方程。以氣井絕對(duì)無(wú)阻流量Qaof為因變量，以kh、pe、φ和kh、pe為自變量，采用多元線性回歸分析，綜合研究因變量與其他自變量的關(guān)系(見(jiàn)表2)，得到氣井產(chǎn)能與地層系數(shù)(kh)、孔隙度(φ)、地層壓力(pe)相關(guān)系數(shù)矩陣，從表2中可以看出，氣井絕對(duì)無(wú)阻流量與kh、φ、pe成正相關(guān)關(guān)系，其中與kh、pe具有較好的相關(guān)性，各因變量之間的相關(guān)系數(shù)較低，各因變量之間不存在多重共線性。

以氣井絕對(duì)無(wú)阻流量為因變量，分別以kh、pe、φ和kh、pe為自變量，且以常數(shù)是否為零建立4種多元線性回歸模型(見(jiàn)表3)，從表3中可以看出，當(dāng)常數(shù)項(xiàng)為零時(shí)，回歸方程的相關(guān)系數(shù)、測(cè)定系數(shù)即修正測(cè)定系數(shù)均高于常數(shù)不為零時(shí)；即常數(shù)項(xiàng)為零時(shí)，回歸方程的擬合程度較好。

根據(jù)數(shù)理統(tǒng)計(jì)理論，4種模型均是多元線性回歸方程，方程是否成立需要進(jìn)行顯著性檢驗(yàn)。①當(dāng)回歸方程得到的F統(tǒng)計(jì)量(F統(tǒng)計(jì)量計(jì)算方法見(jiàn)文獻(xiàn)[12,13])大于F臨界值時(shí)，方程為真。即在給定顯著性水平α之下。如F>Fα(n,n-2)(n為樣本數(shù))，則拒絕假設(shè)，認(rèn)為回歸方程整體是顯著的；②F統(tǒng)計(jì)量對(duì)應(yīng)的P值表示當(dāng)原假設(shè)為真時(shí)，比所得到的樣本觀察結(jié)果更極端的結(jié)果出現(xiàn)的概率，即為棄真概率，所以P值越小越好，其計(jì)算方法見(jiàn)文獻(xiàn)[12，13]。

對(duì)于給定的顯著性水平α=0.05，查F分布表[12,13]得到4種模型的F臨界值分別為：F0.05(3,28)=2.98、F0.05(2,28)=3.35、F0.05(3,28)=2.96、F0.05(2,28)=3.34，而4種模型的F統(tǒng)計(jì)量值分別為39.03、59.98、102.77、157.24，遠(yuǎn)大于相應(yīng)的F臨界值；同時(shí)對(duì)應(yīng)的P值均小于0.0001，即所得到的樣本觀察結(jié)果更極端的結(jié)果出現(xiàn)的概率很小，說(shuō)明4種模型的置信度達(dá)到99.99%以上。由此可知4種模型的線性回歸均為真的概率大。

表3 不同模型線性回歸分析結(jié)果

在不同模型線性回歸分析后，通過(guò)回歸系數(shù)檢驗(yàn)，確定最優(yōu)模型。對(duì)方程線性系數(shù)采取t分布檢驗(yàn)，該檢驗(yàn)主要用來(lái)判斷因變量與自變量之間的顯著程度，主要的參數(shù)為t統(tǒng)計(jì)量及對(duì)應(yīng)的P值[12,13]，兩個(gè)參數(shù)通過(guò)概率學(xué)t分布表查詢和計(jì)算[12,13]，一般情況下t>2則說(shuō)明這個(gè)參數(shù)顯著，也就是說(shuō)對(duì)模型的貢獻(xiàn)量比較大，是不可以剔除的參數(shù)，P值與上述F分布檢驗(yàn)中參數(shù)含義一致，為棄真概率[12,13]。通過(guò)線性回歸得到絕對(duì)無(wú)阻流量預(yù)測(cè)方程的系數(shù)如表4，查t分布表得在置信度為95%時(shí)，t0.05(28)=2.0484，從4種模型線性回歸結(jié)果來(lái)看(見(jiàn)表4)，模型4的t檢驗(yàn)值均大于t0.05(28)=2.0484，且對(duì)應(yīng)的P值均遠(yuǎn)小于0.0001，方程為真的概率高于99.99%，故說(shuō)明模型4為最優(yōu)回歸方程：

Qaof=1.04pe+0.47kh

(6)

4 應(yīng)用實(shí)例

以川東地區(qū)DMP長(zhǎng)興組氣藏為例，該氣藏為典型的生物礁氣藏。氣藏已完鉆井5口，均為氣井。其中，YA012-2井于2009年11月最先投產(chǎn)，該井測(cè)試產(chǎn)量101.78×104m3/d；2011年又增加2口滾動(dòng)勘探開(kāi)發(fā)井(YA012-6、012-X7井)，分別在2010～2013年間相繼投產(chǎn)；2017年部署開(kāi)發(fā)井YA012-12井。根據(jù)氣藏生產(chǎn)情況，該氣藏具有以下特征：①早期投產(chǎn)井表明，該氣藏為統(tǒng)一壓力系統(tǒng)，氣藏連通性好；②氣藏滲透性好，氣井產(chǎn)能高，YA012-1井及YA012-X7井井區(qū)滲透率最高，分別為7.12、10.23mD，并且早期投產(chǎn)4口井單井日均產(chǎn)能達(dá)到了211.8×104m3/d。

根據(jù)YA012-6井的壓力恢復(fù)試井資料，結(jié)合地震等靜態(tài)地質(zhì)認(rèn)識(shí)認(rèn)為該氣藏南端滲透性是逐漸變差，據(jù)此繪制了地層壓力分布和氣藏kh圖(見(jiàn)圖2、3)，獲得擬部署井YA012-12井kh為157.2mD·m，地層壓力為35.1MPa，根據(jù)式(6)計(jì)算得到擬鉆井YA012-12井絕對(duì)無(wú)阻流量110.38×104m3/d。2017年YA012-12井投產(chǎn)，投產(chǎn)后進(jìn)行了穩(wěn)定試井獲得絕對(duì)無(wú)阻流量104.31×104m3/d，與筆者方法預(yù)測(cè)值比較，相對(duì)誤差為5.50%，結(jié)果比較可靠。

圖2 DMP長(zhǎng)興組氣藏2017年地層壓力分布圖 圖3 DMP長(zhǎng)興組kh值分布圖 Fig.2 Distribution of formation pressure of Changxing Formation of DMP in 2017 Fig.3 Distribution of kh of Changxing Formation of DMP

TD69、TD18、W89、W58、QL41、QL9均為川東氣田裂縫-孔隙型氣藏生產(chǎn)井，氣井投產(chǎn)初期均進(jìn)行過(guò)穩(wěn)定試井，以穩(wěn)定試井絕對(duì)無(wú)阻流量作為標(biāo)準(zhǔn)，使用上述產(chǎn)能預(yù)測(cè)方程計(jì)算氣井絕對(duì)無(wú)阻流量，并與穩(wěn)定試井絕對(duì)無(wú)阻流量進(jìn)行對(duì)比，從而驗(yàn)證預(yù)測(cè)方程的可靠性。從隨機(jī)選取的6口氣井和ZG4井的計(jì)算結(jié)果來(lái)看(見(jiàn)表5)，產(chǎn)能預(yù)測(cè)方程計(jì)算結(jié)果與穩(wěn)定試井絕對(duì)無(wú)阻流量相對(duì)誤差均在10%以下，計(jì)算誤差較為合理。

表5 氣井產(chǎn)能計(jì)算結(jié)果對(duì)比表

5 結(jié)論與建議

1)根據(jù)氣井二項(xiàng)式產(chǎn)能方程及前人研究成果分析，影響氣井絕對(duì)無(wú)阻流量的主要因素為氣井地層系數(shù)(kh)、地層壓力(pe)、儲(chǔ)層孔隙度(φ)和動(dòng)態(tài)儲(chǔ)量。

2)通過(guò)對(duì)氣井實(shí)際生產(chǎn)數(shù)據(jù)的分析，氣井絕對(duì)無(wú)阻流量與地層系數(shù)(kh)、地層壓力(pe)、儲(chǔ)層孔隙度(φ)及動(dòng)態(tài)儲(chǔ)量均呈正相關(guān)關(guān)系，其中與地層系數(shù)(kh)及地層壓力(pe)相關(guān)系數(shù)較高，說(shuō)明地層系數(shù)與地層壓力對(duì)氣井絕對(duì)無(wú)阻流量影響較為顯著；與儲(chǔ)層孔隙度相關(guān)關(guān)系較小，說(shuō)明儲(chǔ)層孔隙度對(duì)氣井產(chǎn)能影響較小。

3)通過(guò)多元線性回歸與逐步回歸相結(jié)合的方法，得到了氣井的絕對(duì)無(wú)阻流量預(yù)測(cè)方程，將新建立的方程在實(shí)際生產(chǎn)中進(jìn)行驗(yàn)證應(yīng)用，發(fā)現(xiàn)新方程計(jì)算誤差較小，該方程的建立可用于未鉆開(kāi)發(fā)井的絕對(duì)無(wú)阻流量預(yù)測(cè)，對(duì)于指導(dǎo)鉆井決策具有一定的指導(dǎo)意義，對(duì)將數(shù)學(xué)地質(zhì)應(yīng)用在天然氣井產(chǎn)能評(píng)價(jià)上具有一定的參考價(jià)值。

4)建立的預(yù)測(cè)方程是針對(duì)川東氣田裂縫-孔隙型氣藏，在應(yīng)用方程時(shí)應(yīng)注意氣藏類(lèi)型。由于該方程需要獲取地層系數(shù)(kh)和地層壓力(pe)，因此，在應(yīng)用時(shí)對(duì)于儲(chǔ)層橫向分布穩(wěn)定，氣藏連通關(guān)系比較清楚氣藏的開(kāi)發(fā)補(bǔ)充井預(yù)測(cè)精度較高。