郭艷東

(1.中國(guó)石油化工股份有限公司石油勘探開(kāi)發(fā)研究院，北京 100083；2.頁(yè)巖油氣富集機(jī)理與有效開(kāi)發(fā)國(guó)家重點(diǎn)實(shí)驗(yàn)室，北京 100083)

0 引 言

頁(yè)巖氣是指主體位于暗色泥頁(yè)巖或高碳泥頁(yè)巖中的天然氣聚集[1]。頁(yè)巖氣藏屬于非常規(guī)氣藏，主要以游離氣和吸附氣兩種形式存在，游離氣主要賦存在頁(yè)巖的各種孔隙和各級(jí)裂縫中，吸附氣主要吸附在頁(yè)巖中有機(jī)質(zhì)和黏土礦物表面[2]。頁(yè)巖氣藏地質(zhì)特征復(fù)雜，儲(chǔ)集空間具有多尺度特征，存在納米孔隙、微米孔隙、微裂隙、裂縫等，頁(yè)巖的孔隙度低，滲透率極低。頁(yè)巖氣井一般沒(méi)有自然產(chǎn)能，需采用水平井和大規(guī)模水力壓裂的方式進(jìn)行開(kāi)發(fā)。因此頁(yè)巖氣井的產(chǎn)能評(píng)價(jià)難度較大，國(guó)內(nèi)外學(xué)者主要從遞減曲線分析法、產(chǎn)能解析公式和數(shù)值模擬法三個(gè)方面進(jìn)行了相關(guān)研究。

在遞減曲線分析方面，ARPS[3]提出了產(chǎn)量遞減分析雙曲模型；DILHAN[4]提出了修正的ARPS模型;FETKOVICH等[5]指出在擬穩(wěn)定流動(dòng)階段產(chǎn)量遞減符合雙曲模型；DUONG[6]建立了產(chǎn)量指數(shù)遞減方程。遞減曲線方法都是根據(jù)氣井的遞減趨勢(shì)進(jìn)行擬合，然后預(yù)測(cè)氣井產(chǎn)能，這就需要?dú)饩M(jìn)入遞減階段生產(chǎn)一定的時(shí)間，如果氣井還處于定產(chǎn)生產(chǎn)階段或者放噴生產(chǎn)時(shí)間較短，均不適用于遞減分析方法。在產(chǎn)能解析公式方面，有學(xué)者綜合考慮頁(yè)巖氣解吸、擴(kuò)散等滲流特征，建立并求解頁(yè)巖氣藏不穩(wěn)定滲流數(shù)學(xué)模型，劃分了流動(dòng)階段，開(kāi)展了產(chǎn)能影響因素分析，預(yù)測(cè)氣井產(chǎn)能[7-9]；尹洪軍等[10]考慮頁(yè)巖大型壓裂改造特征將儲(chǔ)層分為五個(gè)區(qū)，建立了五區(qū)復(fù)合產(chǎn)能模型；顧岱鴻等[11]將頁(yè)巖儲(chǔ)層看成三重介質(zhì)，建立了頁(yè)巖氣井產(chǎn)能預(yù)測(cè)模型；田冷等[12]、姜瑞忠等[13]建立了考慮頁(yè)巖的儲(chǔ)層應(yīng)力敏感特征的產(chǎn)能預(yù)測(cè)方法。解析公式法需要確定相關(guān)地質(zhì)和工程參數(shù)，建立較為精確的地質(zhì)模型，才能進(jìn)行產(chǎn)能預(yù)測(cè)。在數(shù)值模擬法方面，樊冬艷等[14]、趙金洲等[15]采用數(shù)值模擬的方法研究了考慮頁(yè)巖氣微觀滲流機(jī)理的壓裂井產(chǎn)能預(yù)測(cè)。采用數(shù)值模擬法進(jìn)行產(chǎn)能預(yù)測(cè)時(shí)，需要建立頁(yè)巖氣藏三維地質(zhì)模型。

上述頁(yè)巖氣井產(chǎn)能預(yù)測(cè)方法均沒(méi)有涉及頁(yè)巖氣開(kāi)發(fā)的初期，利用物質(zhì)平衡方程結(jié)合產(chǎn)能方程對(duì)頁(yè)巖氣井產(chǎn)能進(jìn)行快速預(yù)測(cè)的相關(guān)內(nèi)容。為此，本文在頁(yè)巖氣物質(zhì)平衡方程和產(chǎn)能方程研究的基礎(chǔ)上，建立了基于物質(zhì)平衡方程的頁(yè)巖氣井產(chǎn)能預(yù)測(cè)方法。該方法需要參數(shù)較少，可在開(kāi)發(fā)初期沒(méi)有建立地質(zhì)模型時(shí)對(duì)頁(yè)巖氣井產(chǎn)能進(jìn)行快速預(yù)測(cè)。同時(shí)，該方法綜合考慮了頁(yè)巖氣吸附氣解吸、異常高壓氣藏巖石彈性能量的影響，產(chǎn)能預(yù)測(cè)結(jié)果準(zhǔn)確、可靠。

1 頁(yè)巖氣物質(zhì)平衡方程

1.1 剩余游離氣儲(chǔ)量

頁(yè)巖氣藏開(kāi)采初期，氣藏地層壓力隨著天然氣的產(chǎn)出不斷下降，隨著地層壓力的下降，必將引起天然氣發(fā)生膨脹作用、儲(chǔ)氣層的壓實(shí)和巖石顆粒的彈性膨脹作用、地層束縛水的彈性膨脹作用[16]。巖石的壓實(shí)和顆粒膨脹以及地層束縛水的膨脹會(huì)占據(jù)巖石地下孔隙體積，使得游離氣占據(jù)的地下孔隙體積減少。

假設(shè)氣藏原始游離氣儲(chǔ)量為Gf，當(dāng)?shù)貙訅毫τ蓀i下降到p時(shí)，由于流體壓力下降及有效應(yīng)力升高，游離氣占據(jù)地下孔隙體積減少量表達(dá)為式(1)。

(1)

式中：ΔV為游離氣占據(jù)地下孔隙體積減少量，104m3；Gf為氣藏原始游離氣儲(chǔ)量，104m3；Bgi為初始地層狀態(tài)下的氣體體積系數(shù)；Sgi為氣藏原始含氣飽和度，1；Swi為氣藏原始含水飽和度，1；pi為原始地層壓力，MPa；p為地層壓力，MPa;Cf為地層巖石壓縮系數(shù)，MPa-1；Cw為地層水壓縮系數(shù)，MPa-1。

對(duì)于異常高壓頁(yè)巖氣藏，與常規(guī)異常高壓氣藏一樣，巖石壓縮系數(shù)是地層有效應(yīng)力的函數(shù)[17-18]，參考常規(guī)異常高壓氣藏，巖石壓縮系數(shù)與有效應(yīng)力可以用式(2)擬合。

式中：pob為上覆巖層壓力，MPa；a0、a1、a2、a3為擬合系數(shù)。

將式(2)代入式(1)，然后對(duì)式(1)積分項(xiàng)進(jìn)行積分、化簡(jiǎn)可以得到式(3)。

式(3)即為考慮巖石壓縮系數(shù)隨地層壓力變化，異常高壓頁(yè)巖氣藏由于地層壓力下降引起的游離氣占據(jù)地下孔隙體積變化量。因此，當(dāng)?shù)貙訅毫ο陆档絧時(shí)，根據(jù)式(3)，此時(shí)剩余游離氣儲(chǔ)量N表達(dá)為式(4)。

式中,Bg為地層壓力p時(shí)氣體體積系數(shù)，1。

1.2 剩余吸附氣儲(chǔ)量

頁(yè)巖氣藏游離氣和吸附氣共存，在裂縫孔隙和基質(zhì)孔隙中存在的主要是游離氣，吸附氣則吸附在基質(zhì)孔隙的內(nèi)表面[19]。隨著地層壓力下降，達(dá)到臨界解吸壓力后，吸附在基質(zhì)表面的吸附氣將發(fā)生解吸，通過(guò)Langmuir等溫吸附方程來(lái)描述頁(yè)巖氣藏的吸附氣量[20-21]，見(jiàn)式(5)。

(5)

式中：V(p)為地層壓力p時(shí)頁(yè)巖飽和吸附氣含量，m3/t；VL為蘭氏體積，表示當(dāng)?shù)貙訅毫呌跓o(wú)窮大時(shí)頁(yè)巖的最大理論飽和吸附氣量，m3/t；pL為蘭氏壓力，表示蘭氏曲線中吸附氣量50%對(duì)應(yīng)的壓力，MPa；蘭氏壓力越低，吸附氣在開(kāi)采過(guò)程中越不容易解吸。

根據(jù)式(5)，當(dāng)?shù)貙訅毫閜時(shí)頁(yè)巖氣藏中吸附氣儲(chǔ)量Nxf可以表示為式(6)。

(6)

式中：Nxf為吸附氣儲(chǔ)量，104m3；ρB為頁(yè)巖密度，t/m3；φ為頁(yè)巖孔隙度，1。

1.3 頁(yè)巖氣物質(zhì)平衡方程建立

根據(jù)質(zhì)量守恒原理，可以建立綜合考慮吸附氣解吸和異常高壓影響的頁(yè)巖氣藏物質(zhì)平衡方程,推導(dǎo)過(guò)程如下所述。

在原始地層壓力為pi時(shí)，頁(yè)巖氣藏總儲(chǔ)量G由原始游離氣儲(chǔ)量和原始吸附氣儲(chǔ)量構(gòu)成，可以表示為式(7)。

(7)

天然氣體積系數(shù)定義為式(8)。

式中：psc為標(biāo)準(zhǔn)狀態(tài)壓力，MPa；Tsc為標(biāo)準(zhǔn)狀態(tài)溫度，K；Zsc為標(biāo)準(zhǔn)狀態(tài)下天然氣的偏差系數(shù)；Z為天然氣偏差系數(shù)；T為地層溫度，K。

根據(jù)式(7)和式(8)推導(dǎo)可得原始條件下游離氣儲(chǔ)量Gf與總儲(chǔ)量G之間關(guān)系為式(9)。

(9)

當(dāng)頁(yè)巖氣產(chǎn)出量為Gp，地層壓力降至p時(shí)，根據(jù)物質(zhì)守恒定律，原始游離氣儲(chǔ)量+原始吸附氣儲(chǔ)量=剩余游離氣儲(chǔ)量+剩余吸附氣儲(chǔ)量+累計(jì)產(chǎn)氣量，可以表示為式(10)。

(10)

式中,Gp為累計(jì)產(chǎn)氣量，104m3。將式(4)、式(6)代入式(10)，可以得到式(11)。

胃切除術(shù)引起的體內(nèi)神經(jīng)-內(nèi)分泌紊亂可導(dǎo)致繼發(fā)性PEI[8]。胃切除患者體內(nèi)碳酸氫鹽和脂肪酶分泌顯著降低。胃部分切除患者 PEI的發(fā)生率約70%，全胃切除患者PEI發(fā)生率高達(dá)100%[9]。

(12)

Z*p=

如果不考慮巖石壓縮系數(shù)隨地層有效應(yīng)力的變化，將巖石壓縮系數(shù)看成常數(shù)，則Z*(p)可以簡(jiǎn)化為式(14)。

Z*p=

在氣井生產(chǎn)過(guò)程中，通過(guò)實(shí)測(cè)地層壓力，利用新建立的頁(yè)巖氣物質(zhì)平衡方程式(12)可以計(jì)算出氣藏的總儲(chǔ)量G；當(dāng)氣井繼續(xù)生產(chǎn)至某一時(shí)刻時(shí)，就可以根據(jù)氣井累計(jì)產(chǎn)氣量Gp代入式(12)，計(jì)算該時(shí)刻氣藏的平均地層壓力p。

2 頁(yè)巖氣井二項(xiàng)式產(chǎn)能方程

對(duì)于頁(yè)巖氣井，根據(jù)頁(yè)巖氣地質(zhì)特征和多段壓裂水平井的特點(diǎn)，假設(shè)單條裂縫為一個(gè)流動(dòng)單元，將地層向裂縫的流動(dòng)考慮為變質(zhì)量流，當(dāng)流動(dòng)階段處在擬穩(wěn)態(tài)階段時(shí)，以單條裂縫滲流方程為基礎(chǔ)，通過(guò)疊加原理建立頁(yè)巖氣多段壓裂水平井產(chǎn)能方程，頁(yè)巖氣壓裂水平井二項(xiàng)式產(chǎn)能方程表達(dá)為式(15)[22]；式中A、B分別表示為式(16)、式(17)；式(16)和式(17)中的c表示為式(18)。

(15)

式中：qgsc為頁(yè)巖氣井標(biāo)準(zhǔn)狀態(tài)下日產(chǎn)氣量，104m3/d；T為氣井溫度，K；h為氣層有效厚度，m；km為頁(yè)巖基質(zhì)滲透率，mD；xf為裂縫半長(zhǎng)，m；n為裂縫條數(shù)；kf為裂縫滲透率，mD；w為裂縫寬度，m；ze為裂縫面到阻流邊界距離，m；psc為標(biāo)準(zhǔn)壓力，MPa；Tsc為標(biāo)準(zhǔn)溫度，K；S為表皮系數(shù)；D為高速非達(dá)西系數(shù)。

通過(guò)式(15)可以看出，當(dāng)頁(yè)巖氣井進(jìn)入擬穩(wěn)態(tài)階段時(shí)，頁(yè)巖氣多段壓裂水平井的產(chǎn)能方程滿足二項(xiàng)式產(chǎn)能方程形式，只是方程系數(shù)A和B的表達(dá)式與常規(guī)氣井二項(xiàng)式產(chǎn)能方程不同，可以通過(guò)多工作制度試氣資料，確定頁(yè)巖氣井的二項(xiàng)式產(chǎn)能方程[22]。因此，在頁(yè)巖氣井生產(chǎn)過(guò)程中，只要知道某一時(shí)刻的地層壓力和日產(chǎn)氣量，就可通過(guò)二項(xiàng)式產(chǎn)能方程計(jì)算氣井在該時(shí)刻的井底流壓。

3 基于物質(zhì)平衡方程的產(chǎn)能預(yù)測(cè)方法

根據(jù)頁(yè)巖氣物質(zhì)平衡方程可以計(jì)算頁(yè)巖氣藏不同階段的地層壓力，得到地層壓力后，就可以根據(jù)頁(yè)巖氣二項(xiàng)式產(chǎn)能方程對(duì)不同工作制度的頁(yè)巖氣井產(chǎn)氣量和井底流壓進(jìn)行計(jì)算，因此將頁(yè)巖氣物質(zhì)平衡方程和二項(xiàng)式產(chǎn)能方程結(jié)合起來(lái)就可以對(duì)頁(yè)巖氣井的產(chǎn)能進(jìn)行預(yù)測(cè)?；谝陨纤悸罚疚慕⒘嘶谖镔|(zhì)平衡方程的頁(yè)巖氣井產(chǎn)能預(yù)測(cè)方法，具體方法和流程如下所述，示意圖見(jiàn)圖1。

圖1 頁(yè)巖氣井產(chǎn)能預(yù)測(cè)方法流程示意圖

根據(jù)靜壓測(cè)試資料，采用建立的頁(yè)巖氣物質(zhì)平衡方程式(12)計(jì)算氣井動(dòng)態(tài)儲(chǔ)量G；根據(jù)系統(tǒng)測(cè)試資料，采用式(15)計(jì)算頁(yè)巖氣井的二項(xiàng)式產(chǎn)能方程；按照給定的工作制度進(jìn)行產(chǎn)能預(yù)測(cè)。

1) 定壓生產(chǎn)。流程如下：①根據(jù)氣井二項(xiàng)式產(chǎn)能方程，利用t0時(shí)刻給定的井底流壓pwf計(jì)算t0時(shí)刻日產(chǎn)氣量qg；②根據(jù)日產(chǎn)氣量qg計(jì)算t+dt時(shí)刻的累產(chǎn)氣量Gp，根據(jù)建立的頁(yè)巖氣物質(zhì)平衡方程式(12)，利用氣井動(dòng)態(tài)儲(chǔ)量G和累產(chǎn)氣量Gp計(jì)算t+dt時(shí)刻地層壓力p；③利用t+dt時(shí)刻地層壓力p和給定的井底流壓pwf，根據(jù)氣井二項(xiàng)式產(chǎn)能方程計(jì)算t+dt時(shí)刻日產(chǎn)氣量qg；④下一個(gè)時(shí)間重復(fù)②～③，如果t+dt大于或者等于給定的預(yù)測(cè)時(shí)間tmax，則停止計(jì)算；⑤輸出預(yù)測(cè)期內(nèi)的日產(chǎn)氣量、井底流壓、累產(chǎn)氣量等預(yù)測(cè)結(jié)果。

2) 定產(chǎn)生產(chǎn)。流程如下：①根據(jù)t0時(shí)刻給定的產(chǎn)量qg計(jì)算t+dt時(shí)刻的累產(chǎn)氣量Gp，根據(jù)新建立的頁(yè)巖氣物質(zhì)平衡方程式(12)，利用氣井動(dòng)態(tài)儲(chǔ)量G和累產(chǎn)氣量Gp計(jì)算t+dt時(shí)刻地層壓力p；②根據(jù)氣井二項(xiàng)式產(chǎn)能方程，利用t+dt時(shí)刻地層壓力p計(jì)算t+dt時(shí)刻井底流壓pwf；③將t+dt時(shí)刻井底流壓pwf與給定的極限井底流壓pwfmin比較，如果pwf>pwfmin則繼續(xù)重復(fù)①～③步，如果pwf >=pwfmin則轉(zhuǎn)為定壓生產(chǎn)，令pwf=pwfmin；④如果t+dt大于等于給定的預(yù)測(cè)時(shí)間tmax，則停止計(jì)算；⑤輸出預(yù)測(cè)期內(nèi)的日產(chǎn)氣量、井底流壓、累產(chǎn)氣量等預(yù)測(cè)結(jié)果。

以上計(jì)算方法和流程，可以通過(guò)計(jì)算機(jī)編程來(lái)實(shí)現(xiàn)。

4 實(shí)例計(jì)算

以涪陵地區(qū)龍馬溪組頁(yè)巖氣井X1井地質(zhì)參數(shù)為基礎(chǔ)參數(shù)進(jìn)行實(shí)例計(jì)算，預(yù)測(cè)該井產(chǎn)能。該井產(chǎn)氣層位為五峰-龍馬溪組，埋深2 460 m，地層壓力系數(shù)1.55，為一個(gè)高壓頁(yè)巖氣藏，其他相關(guān)參數(shù)見(jiàn)表1。在開(kāi)采初期，該井進(jìn)行了系統(tǒng)試井和靜壓測(cè)試。

表1 X1井基礎(chǔ)參數(shù)

通過(guò)系統(tǒng)試井測(cè)試確定X1井的產(chǎn)能方程表達(dá)為式(19)。

(19)

通過(guò)不同階段靜壓測(cè)試，采用頁(yè)巖氣物質(zhì)平衡方程評(píng)價(jià)該井動(dòng)態(tài)儲(chǔ)量，目前X1井總動(dòng)態(tài)儲(chǔ)量為1.91×108m3，其中游離氣儲(chǔ)量為1.2×108m3，見(jiàn)圖2。

采用本文建立的基于物質(zhì)平衡方程的頁(yè)巖氣井產(chǎn)能預(yù)測(cè)方法，對(duì)該井分別按照定壓生產(chǎn)和定產(chǎn)生產(chǎn)兩種方式，考慮吸附氣和不考慮吸附氣兩種情況進(jìn)行了產(chǎn)能預(yù)測(cè)，同時(shí)也采用數(shù)值模擬方法對(duì)該井進(jìn)行了產(chǎn)能預(yù)測(cè)，具體計(jì)算結(jié)果見(jiàn)表2和圖3、圖4。為了研究異常高壓頁(yè)巖氣藏巖石彈性能量對(duì)產(chǎn)能預(yù)測(cè)的影響， 將巖石壓縮系數(shù)Cf取為常數(shù)， 分別對(duì)考慮和不考慮巖石壓縮系數(shù)Cf兩種情況產(chǎn)能進(jìn)行了預(yù)測(cè)，結(jié)果見(jiàn)表3和圖5。

圖2 頁(yè)巖氣物質(zhì)平衡方程動(dòng)態(tài)儲(chǔ)量評(píng)價(jià)曲線

評(píng)價(jià)方法生產(chǎn)方式是否考慮吸附氣初期日產(chǎn)氣量/(104m3/d)穩(wěn)產(chǎn)時(shí)間/a20年累產(chǎn)氣量/(104m3)基于物質(zhì)平衡方程的產(chǎn)能預(yù)測(cè)法定壓不考慮18.6不穩(wěn)產(chǎn)9 720Pwf=7 MPa考慮18.6不穩(wěn)產(chǎn)11 744定產(chǎn)不考慮63.09 718Pwf min=7 MPa考慮63.211 690數(shù)值模擬法定產(chǎn)Pwf min=7 MPa考慮63.312 800

表3 X1井物質(zhì)平衡產(chǎn)能預(yù)測(cè)法開(kāi)發(fā)指標(biāo)預(yù)測(cè)結(jié)果

圖3 X1井定壓生產(chǎn)開(kāi)發(fā)指標(biāo)曲線

圖4 X1井定產(chǎn)生產(chǎn)開(kāi)發(fā)指標(biāo)曲線

圖5 X1井定產(chǎn)生產(chǎn)不同Cf情況開(kāi)發(fā)指標(biāo)曲線

X1井試采期間按照6×104m3/d左右的產(chǎn)量定產(chǎn)生產(chǎn)，穩(wěn)產(chǎn)3.34年，穩(wěn)產(chǎn)期末累產(chǎn)氣7 715×104m3，已進(jìn)入遞減階段。根據(jù)表2計(jì)算結(jié)果可知，采用本文建立的頁(yè)巖氣物質(zhì)平衡快速產(chǎn)能評(píng)價(jià)方法在定產(chǎn)生產(chǎn)考慮吸附氣、考慮巖石壓縮系數(shù)影響情況下，預(yù)測(cè)X1井穩(wěn)產(chǎn)期為3.2年，20年末累產(chǎn)氣量為11 690×104m3；采用數(shù)值模擬法預(yù)測(cè)X1井穩(wěn)產(chǎn)期為3.3年，20年末累產(chǎn)氣量為12 800×104m3；筆者建立的方法與實(shí)際試采數(shù)據(jù)、數(shù)值模擬方法相比穩(wěn)產(chǎn)期接近，因?yàn)樵摲椒](méi)有考慮外圍補(bǔ)給，預(yù)測(cè)累產(chǎn)氣量略低于數(shù)值模擬的預(yù)測(cè)結(jié)果，20年末累產(chǎn)氣量誤差僅為8.6%，符合率91.4%。因此在開(kāi)發(fā)初期沒(méi)有建立地質(zhì)模型時(shí)，可以采用頁(yè)巖氣物質(zhì)平衡快速產(chǎn)能評(píng)價(jià)方法快速對(duì)頁(yè)巖氣井的產(chǎn)能進(jìn)行預(yù)測(cè)，計(jì)算結(jié)果可靠。

根據(jù)X1井不同情況產(chǎn)能預(yù)測(cè)結(jié)果對(duì)生產(chǎn)方式、吸附氣和巖石壓縮系數(shù)的影響進(jìn)行研究。

根據(jù)表2的計(jì)算結(jié)果可知，X1井定壓生產(chǎn)與定產(chǎn)生產(chǎn)相比，20年末累產(chǎn)氣量相差不大，定壓生產(chǎn)略高于定產(chǎn)生產(chǎn)。在不考慮應(yīng)力敏感和井筒影響的情況下，生產(chǎn)方式對(duì)氣井累產(chǎn)氣量影響不大。

根據(jù)表2計(jì)算結(jié)果可知：X1井定壓生產(chǎn)，考慮吸附氣較不考慮吸附氣，20年末累產(chǎn)氣量增加2 024×104m3，占總累產(chǎn)氣量的17%，通過(guò)對(duì)比累產(chǎn)氣量曲線可以看出，吸附氣的貢獻(xiàn)在2.5年逐漸體現(xiàn)出來(lái)，見(jiàn)圖3；X1井定產(chǎn)生產(chǎn)，考慮吸附氣與不考慮吸附氣，20年末累產(chǎn)氣量增加1 972×104m3，占總累產(chǎn)氣量的16.9%，穩(wěn)產(chǎn)期增加0.2年，吸附氣對(duì)穩(wěn)產(chǎn)期的影響不大，見(jiàn)圖4；定壓生產(chǎn)初期產(chǎn)量高，生產(chǎn)壓差大，初期地層壓力下降快，而采取定產(chǎn)方式生產(chǎn)，初期產(chǎn)量得到控制，地層壓力下降慢，因此吸附氣的貢獻(xiàn)在生產(chǎn)4年逐漸體現(xiàn)出來(lái)，比定壓生產(chǎn)要晚。

根據(jù)表3的計(jì)算結(jié)果可知，X1井定產(chǎn)6×104m3/d生產(chǎn)時(shí)，考慮巖石壓縮系數(shù)Cf影響較不考慮，穩(wěn)產(chǎn)時(shí)間增加0.2年，20年末累產(chǎn)氣量增加172×104m3，見(jiàn)圖5。對(duì)于異常高壓頁(yè)巖氣藏，考慮巖石壓縮系數(shù)影響，頁(yè)巖的彈性能量的貢獻(xiàn)使X1井穩(wěn)產(chǎn)時(shí)間和累產(chǎn)氣有一定的增加，但增加不明顯。

5 結(jié) 論

1) 建立了基于物質(zhì)平衡方程的頁(yè)巖氣井產(chǎn)能預(yù)測(cè)方法，該方法可以在開(kāi)發(fā)初期對(duì)頁(yè)巖氣井的產(chǎn)能進(jìn)行快速預(yù)測(cè)，能夠考慮吸附氣和異常高壓氣藏巖石彈性能量的影響，適用于定產(chǎn)和定壓兩種生產(chǎn)方式，預(yù)測(cè)結(jié)果準(zhǔn)確、可靠。

2) 吸附氣對(duì)X1井累產(chǎn)氣影響明顯，吸附氣占總累產(chǎn)氣量的約17%，定產(chǎn)生產(chǎn)時(shí)吸附氣對(duì)氣井穩(wěn)產(chǎn)期影響不大；吸附氣對(duì)頁(yè)巖氣井的貢獻(xiàn)在2.5～4年逐漸體現(xiàn)出來(lái)，定壓放噴生產(chǎn)初期產(chǎn)量高，生產(chǎn)壓差大，地層壓力下降快，而采取定產(chǎn)方式生產(chǎn)，初期產(chǎn)量受到控制，地層壓力下降慢，定壓放噴生產(chǎn)吸附氣的貢獻(xiàn)比定產(chǎn)生產(chǎn)要早。

3) X1井進(jìn)行產(chǎn)能預(yù)測(cè)時(shí)，考慮巖石壓縮系數(shù)的影響，穩(wěn)產(chǎn)時(shí)間和累產(chǎn)氣有所增加，但不明顯。

[1] 張金川,金之鈞,袁明生．頁(yè)巖氣成藏機(jī)理和分布[J].天然氣工業(yè),2004,24(7):15-18．

[2] KING G R.Material-balance techniques for coal-seam Devonian shale gas reservoirs with limited water influx[J].SPE20730,1993:67-72.doi:https:∥doi.org/10.2118/20730-MS.

[3] ARPS J J.Analysis of decline curves[J].Transactions of the Aime,1945,160(1):228-247.

[4] DILHAN IIK.Exponential vs.hyperbolic decline in tight gas sands:understanding the origin and implications for reserve estimates using Arps decline curves[C]∥SPE Annual Technical Conference and Exhibition,21-24 September,Denver,Colorado,USA.doi:https:∥doi.org/10.2118/116731-MS.

[5] FETKOVICH M J,FETKOVICH E J,Fetkovich M D.Useful concepts for decline curve forecasting,reserve estimation and analysis[J].SPE Reservoir Engineering,1996,11(1):13-22.

[6] DUONG Anh N.An unconventional rate decline approach for tight and fracture dominated gas wells[C]∥Canadian Unconventional Resources and International Petroleum Conference,19-21 October,Calgary,Alberta,Canada.doi:https:∥doi.org/10.2118/137748-MS.

[7] 王偉,樊冬艷,孫海.頁(yè)巖氣藏分段多簇壓裂水平井產(chǎn)能模型及影響因素分析[J].科學(xué)技術(shù)與工程,2015,15(31):36-41.

[8] 謝亞雄,劉啟國(guó),王衛(wèi)紅.頁(yè)巖氣藏多段壓裂水平井產(chǎn)能預(yù)測(cè)模型[J].大慶石油地質(zhì)與開(kāi)發(fā),2016,35(5):163-169.

[9] 任俊杰,郭平,王德龍.頁(yè)巖氣藏壓裂水平井產(chǎn)能模型及影響因素[J].東北石油大學(xué)學(xué)報(bào),2012,36(6):76-81.

[10] 尹洪軍,趙二猛,付京.頁(yè)巖氣藏壓裂水平井五區(qū)復(fù)合模型產(chǎn)能分析[J].西南石油大學(xué)學(xué)報(bào):自然科學(xué)版,2015(3):9-16.

[11] 顧岱鴻,丁道權(quán),劉軍.三重介質(zhì)頁(yè)巖氣藏分段壓裂水平井產(chǎn)能預(yù)測(cè)模型[J].大慶石油地質(zhì)與開(kāi)發(fā),2016,35(1):158-165.

[12] 田冷,肖聰,顧岱鴻.考慮應(yīng)力敏感與非達(dá)西效應(yīng)的頁(yè)巖氣產(chǎn)能模型[J].天然氣工業(yè),2014,34(12):70-75.

[13] 姜瑞忠,韓光偉,汪洋.低滲氣藏特殊滲流機(jī)理及穩(wěn)定產(chǎn)能預(yù)測(cè)方法[J].石油鉆采工藝,2015(4):67-71.

[14] 樊冬艷,姚軍,孫海.考慮多重運(yùn)移機(jī)制耦合頁(yè)巖氣藏壓裂水平井?dāng)?shù)值模擬[J].力學(xué)學(xué)報(bào),2015,47(6):906-915.

[15] 趙金洲,李志強(qiáng),胡永全.考慮頁(yè)巖儲(chǔ)層微觀滲流的壓裂產(chǎn)能數(shù)值模擬[J].天然氣工業(yè),2015,35(6):53-58.

[16] 陳元千,胡建國(guó).確定異常高壓氣藏地質(zhì)儲(chǔ)量和有效壓縮系數(shù)的新方法[J].天然氣工業(yè),1993,13(1):53-58.

[17] 劉傳喜,劉華,王衛(wèi)紅.基于物質(zhì)平衡法確定超高壓氣藏動(dòng)態(tài)儲(chǔ)量的新方法[J].天然氣工業(yè),2009,29(12):68-70.

[18] 王衛(wèi)紅,郭艷東,劉華.超高壓氣藏單井控制動(dòng)態(tài)儲(chǔ)量評(píng)價(jià)方法[J].新疆石油地質(zhì),2011,32(3):291-294.

[19] 賈承造,鄭民,張永峰.中國(guó)非常規(guī)油氣資源與勘探開(kāi)發(fā)前景[J].石油勘探與開(kāi)發(fā),2012,39(2):129-136.

[20] 曹廷寬,段永剛,王榮.考慮多組分吸附的頁(yè)巖氣儲(chǔ)量計(jì)算[J].油氣藏評(píng)價(jià)與開(kāi)發(fā),2011,1(6):68-70.

[21] LANGMUIR I.The adsorption of gases on plane surfaces of Glass,Mica and Platinum[J].Journal of American Chemical Society,1918,40(9):1361-1370.

[22] 劉華,胡小虎,王衛(wèi)紅.頁(yè)巖氣壓裂水平井?dāng)M穩(wěn)態(tài)階段產(chǎn)能評(píng)價(jià)方法研究[J].西安石油大學(xué)學(xué)報(bào):自然科學(xué)版,2016,31(2):76-81.