多級壓裂頁巖氣水平井的不穩(wěn)定生產(chǎn)數(shù)據(jù)分析

李 欣 段勝楷 侯大力 孫 揚 李士倫

1．美國 Knowledge Reservoir咨詢公司 2．Chevron Energy Technology Company 3．“油氣藏地質(zhì)及開發(fā)工程”國家重點實驗室·西南石油大學(xué)

1 概念

多級水力壓裂水平井（MFHW）完井技術(shù)是開發(fā)頁巖氣資源的有效技術(shù)［1］。但MFHW產(chǎn)生的壓裂氣體體積（SRV）改變了多孔介質(zhì)的流態(tài)［2－3］。壓裂出的SRV的形態(tài)及其增強的滲流能力，對于井的生產(chǎn)動態(tài)乃至最終采收率都有很大的影響。目前通過改進的具有天然裂縫的致密氣藏模型解析解［3－9］描述 MFHW的不穩(wěn)定直線流。對文中用到的一些概念和MFHW的流動階段［10－12］介紹如下。

SRV是頁巖氣藏中最接近井筒部分，1口 MFHW有多個SRV（圖1）。SRV的邊界是裂縫半長的位置。

XRV是單井泄流面積的外部未被壓裂或有壓裂縫卻未被支撐劑撐開的氣藏體積（圖1）。XRV的大小取決于MFHW在一個地質(zhì)單元里的井距。頁巖氣井的流動先發(fā)生在SRV，隨著SRV中壓力的遞減，XRV中的頁巖氣才開始向SRV中供給，所以頁巖氣多級裂縫水平井流動狀態(tài)分以下5個階段：①SRV直線流；②SRV衰竭流；③復(fù)合直線流（圖2）；④擬徑向流；⑤油氣藏系統(tǒng)的邊界流。

2 生產(chǎn)數(shù)據(jù)不穩(wěn)定流分析方法

多級壓裂井早期的壓力和產(chǎn)量數(shù)據(jù)，主要是基于Wattenbarger［8］等人的工作，經(jīng)過改進后，可用于分析MFHW早期不穩(wěn)定直線流。筆者在后面給出了詳細的公式推導(dǎo)。

2．1 用log－log診斷圖識別直線流

2．2 Δ［m（p）］／qg和分析圖

時間疊加的開方圖是一個更準(zhǔn)確地確定和描述直線流的圖。無論是定產(chǎn)還是定壓，在Δ［m（p）］／qg和圖中直線段就是log－log診斷圖中的1／2直線段。

圖1 SRV和XRV示意圖［9］

圖2 MFHW復(fù)合直線流圖

圖3 MFHW的壓降變化圖

圖4 log－log診斷圖確定定產(chǎn)和定壓條件下的直線流圖

圖5 Δ［m（p）］／qg和分析圖

2．3 SRV氣藏性質(zhì)的計算

該方法的優(yōu)點是在氣藏的φ、Κ、h和A不確定的情況下，也能計算SRV的Vp和OGIP。如果所有的生產(chǎn)歷史數(shù)據(jù)仍然是1／2斜率段或無限導(dǎo)流直線段（即鄰近裂縫的壓力波沒有交叉），采用表1中定產(chǎn)或定壓表達式計算的SRV的Vp和OGIP視為最低值。從上面的介紹中可以看到有3個變量在SRV裂縫系統(tǒng)早期直線流的分析中至關(guān)重要：①SRV直線流段的斜率，珦mCRL或者珦mCPL；②直線流的結(jié)束時間，tef；③Y軸的截距，bCRL或者bCPL。

表1 定產(chǎn)和定壓條件下直線性流段公式表

3 公式推導(dǎo)

3．1 裂縫井直線流的解析解

以下的公式適用于液體的裂縫直線流。如果采用氣體的擬壓力代替pwD，則適用于氣井。下面的兩個方程分別是定產(chǎn)和定井底壓力內(nèi)邊界條件和封閉線性油藏外邊界條件下的解析解。式（1）是一個封閉有垂直裂縫線性流油藏中1口定產(chǎn)的生產(chǎn)井模型的解：

定pwf生產(chǎn)的解：

3．2 早期段裂縫直線流近似解

早期或無限大地層中裂縫直線流對定產(chǎn)和定壓力氣井的無因次近似解是：

單相氣體流動的解可以采用真實氣體擬壓力，即

m（p）是真實氣體擬壓力，即

相對于裂縫半長的無因次時間：

不同于無限大徑向油藏的解析解，無限大線性油藏的定產(chǎn)量和定井底壓力的解相差為π／2。

在頁巖氣藏中，無因次的時間是基于Xmf，即

Xmf是所有裂縫半長的總和［1］。

Δ［m（p）］／qg和圖直線段的斜率珦mCRL是：

Ac是直線流通過裂縫的總表面積，因為裂縫有兩面，Ac＝4Xmfh，故

這個解隱含的假設(shè)是多級壓裂水平井的多級裂縫的n個水平裂縫的半長，相當(dāng)于一個半長是nxf的超長水平裂縫。

3．3 水平裂縫的邊界

隨著生產(chǎn)時間的增加，井底壓力響應(yīng)的一個重要特征是無限導(dǎo)流直線段的結(jié)束，即1／2斜率段的結(jié)束。這說明不穩(wěn)定直線流向鄰間裂縫干擾流過渡。這個時間可用于計算鄰近的兩條水平縫的距離，因為在定產(chǎn)條件下直線流相對于ye的無因次結(jié)束時間是0．5；在定壓條件下因為在定產(chǎn)條件下直線流相對于ye的無因次結(jié)束時間是0．25。即

但是，在通常的情況下都假設(shè)ye是多級水平井每一級長度的一半，然后計算SRV的滲透率：

如果直線流動段尚未結(jié)束，滲透率的估計值可能偏大。

3．4 早單井泄流面積的計算

在Κ或者Κf未知的情況下，仍然可以通過直線流的斜率計算單井的泄流面積。單井的SRV泄流面積A是：

A可以用直線流的斜率和1／2斜率結(jié)束時間計算：

3．5 孔隙體積的計算

式（17）乘以φ和h就是SRV的孔隙體積，即

3．6 原始天然氣儲量

一旦計算出孔隙體積，原始天然氣儲量就很容易計算：

4 討論

用生產(chǎn)數(shù)據(jù)分析裂縫系統(tǒng)的流動狀態(tài)是解釋氣藏水力壓裂縫空間展布的關(guān)鍵。采用Δ［m（p）］／qg對lnt的導(dǎo)數(shù)曲線可以使曲線平滑，更便于分析。許多生產(chǎn)井的壓力和產(chǎn)量都隨時間波動，對于這些情況，在診斷圖和特征圖上，可能看不到穩(wěn)定的特征斜率段，即使井的流動可能正處于直線流段，只是壓力產(chǎn)量的波動掩蓋了1／2斜率段，使得生產(chǎn)分析變得復(fù)雜。因此，在這些情況下，常規(guī)的log－log特征圖只能認為是一個定性的方法，而不是一個明確的定量分析。

5 結(jié)論

筆者介紹了如何分析MFHW井SRV的不穩(wěn)定直線流數(shù)據(jù)來描述壓裂縫的形態(tài)和壓裂區(qū)域的滲流能力，它可用于壓裂效果的評價和配合單井模型歷史擬合。通過以上的討論得到以下結(jié)論。

1）MFHW井SRV解析解模型適用于氣藏單相氣體流動。如果氣藏中有兩相以上的流動相，要做定量分析，模型就要進行修改。

3）如果井的工作制度經(jīng)常改變或長期關(guān)井，可采用時間疊加方法。

4）非零的截距可用于完井質(zhì)量評價。

符 號 說 明

A為井的泄流面積，ft2，1ft＝0．304 8m，下同；Ac為裂縫流動的截面積，ft2；b為地層污染和更早期流動段的影響值，無量綱；bCPL為定壓圖的直線段截距；bCRL為定產(chǎn)圖直線段截距；Bg為氣體體積系數(shù)，無量綱，scf／scf，1scf＝0．028 3m3，下同；ct為總壓縮系數(shù)，psi－1，1psi＝6．894 8kPa，下同；h為氣層凈厚度，ft；Κ為有效SRV滲透率，mD；m（p）為真實氣擬壓力，psi2／cp，1cp＝1mPa·s，下同為定圖直線段斜率；珦mCPL為定圖直線段斜率；n為多級壓裂水平井的級數(shù)；OGIP為原始天然氣儲量，scf；p為壓力，psi；qg為單井日產(chǎn)量，Mscf／d；Sgi為初始含氣飽和度；tDx為相當(dāng)于裂縫半長的無因次時間；tef為直線流的結(jié)束時間，d；T為油藏溫度，，1℃＝33．8，下同；Vp為孔隙體積；xf為平均裂縫半長，ft；Xmf為多級壓裂水平井的總裂縫半長，ft；ye為每級壓裂距離的半長，ft；φ為孔隙度；Δm（p）為擬壓力差，Δm（p）＝m（pi）－m（pwf）；μg為氣體黏度，cp。

下標(biāo)：i為初始狀態(tài)；sup為疊加。

［1］唐穎，張金川，張琴，等．頁巖氣井水力壓裂技術(shù)及其應(yīng)用分析［J］．天然氣工業(yè)，2010，30（10）：33－38．

［2］付永強，馬發(fā)明，曾立新，等．頁巖氣藏儲層壓裂實驗評價關(guān)鍵技術(shù)［J］．天然氣工業(yè)，2011，31（4）：51－54．

［3］RASDI F，CHU Lifu．Diagnosing fracture network pattern and flow regime aids production performance analysis in unconventional oil reservoirs［C］∥paper 151623presented at the SPE／EAGE European Unconventional Resources Conference and Exhibition，20－22March 2012，Vienna，Austria．New York：SPE，2012．

［4］WATTENBARGER R A，RAMEY Jr H J，et al．Well test interpretation of vertically fractured gas wells［J］．Journal of Petroleum Technology，1969，21（5）：625－632．

［5］CINCO－LEY H，SAMANIEGO V F，DOMINGUEZ A N．Transient pressure behavior for a well with a finite conductivity vertical fracture［J］．SPE Journal，1978，18（4）：253－264．

［6］GRINGARTEN A C．Reservoir limit testing for fractured wells［C］∥paper 7452－MS presented at the 53rdSPE Annual Fall Technical Conference and Exhibition，1－3October 1978，Houston，Texas，USA．New York：SPE，1978．

［7］AGARWAL R G，CARTER R D，POLLOCK C B．Evaluation and performance prediction of low－permeability gas wells simulated by massive hydraulic fracturing［J］．Journal of Petroleum Technology，1979，31（3）：362－372．

［8］WATTENBARGER R A，EL－BANBI A H，VILLEGAS M E，et al．Production analysis of linear flow into fractured tight gas wells［C］∥paper 39931presented at the SPE Rocky Mountain Regional／Low－Permeability Reservoirs Symposium and Exhibition，5－8April 1998，Denver，Colorado，USA．New York：SPE，1998．

［9］AREVALO－VILLAGRAN J A，WATTENBARGERR A，SAMANIEGO－VERDUZCO F，et al．Production analysis of long－term linear flow in tight gas reservoirs：case histories［C］∥paper 71516presented at SPE Annual Technical Conference and Exhibition，30September －3October 2001，New Orleans，Louisiana，USA．New York：SPE，2001．

［10］IBRAHIM M，WATTENBARGER R A．Analysis of rate dependence in transient linear flow in tight gas wells［C］∥paper 2005－057presented at Canadian International Petroleum Conference，7－9June 2005，Calgary，Alberta，Canada．New York：SPE，2011．

［11］NAIR H，MILLER M．Northside：Cross－play shale gas performance analysis［C］∥SPE Gulf Coast Section Luncheon，8May 2012，Houston，USA．New York：SPE，2012．

［12］SONG B，EHLIG－ECONOMIDES C A．Rate－normalized pressure analysis for determination of shale gas well performance［C］∥paper 144013presented at North American Unconventional Gas Conference and Exhibition，14－16 June 2011，Woodlands，Texas，USA．New York：SPE，2011．