雷丹鳳 ，王莉 ，張曉偉，王南，于榮澤，徐鳳廷

（1.中國石油勘探開發(fā)研究院廊坊分院，河北 廊坊 065007；2.國家能源頁巖氣研發(fā)（實驗）中心，河北 廊坊 065007；3.中國石油非常規(guī)油氣重點實驗室，河北 廊坊 065007；4.中國石油遼河油田公司高升采油廠，遼寧 盤錦 124125）

頁巖氣井?dāng)U展指數(shù)遞減模型研究

雷丹鳳1，2，3，王莉1，2，3，張曉偉1，2，3，王南1，2，3，于榮澤1，2，3，徐鳳廷4

（1.中國石油勘探開發(fā)研究院廊坊分院，河北 廊坊 065007；2.國家能源頁巖氣研發(fā)（實驗）中心，河北 廊坊 065007；3.中國石油非常規(guī)油氣重點實驗室，河北 廊坊 065007；4.中國石油遼河油田公司高升采油廠，遼寧 盤錦 124125）

頁巖氣藏滲透率極低，用傳統(tǒng)產(chǎn)量遞減模型對頁巖氣井進(jìn)行預(yù)測，結(jié)果過于樂觀。為尋求適用于頁巖氣井的產(chǎn)能評價方法，以Arps遞減分析方法為基礎(chǔ)，對擴展指數(shù)遞減模型展開研究。在模型應(yīng)用過程中，將傳統(tǒng)的遞減曲線分析方法與概率預(yù)測框架相結(jié)合，根據(jù)單井初始遞減率進(jìn)行同類井的劃分，在計算單井遞減曲線參數(shù)的基礎(chǔ)上，利用P指數(shù)概率準(zhǔn)則求取具有代表性的井組遞減曲線參數(shù)值，據(jù)此對單井進(jìn)行概率性預(yù)測，從而規(guī)避了頁巖氣藏復(fù)雜性及上產(chǎn)措施對產(chǎn)量遞減分析的影響，提高了頁巖氣藏產(chǎn)量動態(tài)預(yù)測的準(zhǔn)確性。利用美國多口頁巖氣井的實際生產(chǎn)資料，進(jìn)一步驗證了該方法在頁巖儲層應(yīng)用的準(zhǔn)確性與可行性。

頁巖氣藏；產(chǎn)量預(yù)測；Arps遞減模型；擴展指數(shù)遞減模型

0 引言

在油氣井的整個生產(chǎn)周期中，主要包括產(chǎn)量上升、穩(wěn)定和遞減3個階段。在產(chǎn)量遞減階段，遞減模型的選用對于生產(chǎn)動態(tài)的預(yù)測以及最終采收率的評估具有重要影響［1-2］。

油氣田開發(fā)實踐證明，Arps產(chǎn)量遞減曲線分析是一種可靠的統(tǒng)計方法，可用于預(yù)測擬穩(wěn)定狀態(tài)下的未來動態(tài)產(chǎn)量和油氣藏可采儲量，因此，在常規(guī)油氣藏生產(chǎn)中得到了廣泛的應(yīng)用［3-8］。傳統(tǒng)的遞減曲線分析方法均是基于Arps遞減模型得到的。然而，在頁巖氣藏開發(fā)過程中，由于儲層滲透率極小，氣體需要很長時間才能達(dá)到穩(wěn)定流態(tài)，而且由于氣體的吸附與解吸，使其在基質(zhì)中的流動不穩(wěn)定，各層位的流動狀態(tài)不斷發(fā)生變化［9-11］。此時，應(yīng)用常規(guī)遞減曲線分析方法進(jìn)行頁巖氣井產(chǎn)能預(yù)測，將會產(chǎn)生較大誤差［1，12-14］。

為此，本文在Arps遞減模型的基礎(chǔ)上，對擴展指數(shù)遞減模型展開研究，將其應(yīng)用于頁巖氣井產(chǎn)量遞減趨勢預(yù)測中，并借用美國2個頁巖氣藏多口頁巖氣井的實際生產(chǎn)資料，對模型的可靠性進(jìn)行了驗證。

1 模型的建立

擴展指數(shù)遞減模型由指數(shù)遞減函數(shù)積分形成，受指數(shù)遞減中大量體積微元的影響，具有不均一性。與Arps遞減模型相同，擴展指數(shù)遞減模型亦建立于經(jīng)驗分析之上，它通過將傳統(tǒng)遞減曲線分析法與概率預(yù)測框架相結(jié)合，提高了非常規(guī)儲層油氣井產(chǎn)量遞減趨勢預(yù)測的準(zhǔn)確性。

擴展指數(shù)遞減模型與Arps雙曲遞減模型的對比關(guān)系見表1。在擴展指數(shù)遞減模型中，參數(shù)q0，n，τ為經(jīng)驗常數(shù)，其中，n相當(dāng)于Arps遞減模型里的遞減指數(shù)b，τ是特征時間常數(shù)的中值。參數(shù)對（n，τ）決定著特征時間常數(shù)的分布，需通過求解非線性方程組獲取。特征時間常數(shù)存在一個特定的分布，表現(xiàn)為“肥尾”概率分布形式，n越接近0，尾部分布越大。

表1 擴展指數(shù)遞減模型與Arps雙曲遞減模型對比

2 模型的優(yōu)勢

與Arps遞減模型相比，擴展指數(shù)遞減模型具備2個優(yōu)勢：1）對具有天然邊界的任意單井，可進(jìn)行最終可采儲量的有效預(yù)測；2）利用剩余采收率與累計產(chǎn)量之間的直線關(guān)系，可進(jìn)行剩余采收率的推斷［14-15］。

2.1 最終可采儲量的預(yù)測

由表1可以看出，在Arps遞減模型中，隨b值增大，產(chǎn)量遞減率隨之減小，當(dāng)b值大于1時，最終可采儲量Np，∞為無窮大。這顯然與生產(chǎn)實際不符。

擴展指數(shù)遞減模型通過對天然邊界的假設(shè)解決了這一問題。在該模型中，參數(shù)n介于0到1之間，n=1代表指數(shù)遞減。與Arps遞減模型相反，隨著n值的減小，產(chǎn)量遞減率減小。不同n值、b值條件下，繪制2種遞減模型的產(chǎn)量變化曲線（見圖1）?？梢钥闯?，2種遞減模型的指數(shù)遞減曲線是重合的。當(dāng)b為1.4時，Arps遞減模型曲線隨時間延續(xù)而無限延伸，最終可采儲量無窮大；當(dāng)n為0.3時，盡管擴展指數(shù)遞減模型曲線在生產(chǎn)早期與Arps模型類似，但其晚期的遞減規(guī)律揭示了其中天然邊界的存在。

圖1 不同遞減指數(shù)條件下的模型曲線對比

2.2 剩余采收率的推斷

在Arps遞減模型中，由于當(dāng)b大于1時，產(chǎn)量-時間曲線下的面積無限大（見圖1），因此，無法根據(jù)目前的累計產(chǎn)量進(jìn)行采收率預(yù)測。根據(jù)擴展指數(shù)遞減模型中剩余采收率與累計產(chǎn)量之間的關(guān)系式，利用實際生產(chǎn)數(shù)據(jù)繪制相應(yīng)關(guān)系曲線（見圖2），可以看出二者呈現(xiàn)良好的直線關(guān)系。利用這一相關(guān)關(guān)系，可實現(xiàn)剩余采收率的有效預(yù)測。

圖2 擴展指數(shù)遞減模型剩余采收率與累計產(chǎn)量的關(guān)系

3 模型的應(yīng)用

3.1 方法步驟

1）數(shù)據(jù)過濾。當(dāng)早期的生產(chǎn)數(shù)據(jù)存在不規(guī)律的遞減行為時，可將其作為常數(shù)直接加入累計產(chǎn)量中進(jìn)行計算，而只對起主導(dǎo)遞減趨勢作用的數(shù)據(jù)進(jìn)行曲線擬合（見圖3）。

圖3 產(chǎn)量數(shù)據(jù)篩選示意

2）分井組。依據(jù)單井的初始遞減率對井進(jìn)行分組，以獲取遞減規(guī)律相似井井組的遞減模型參數(shù)。

3）參數(shù)設(shè)定。計算井組內(nèi)所有井的參數(shù)n和τ。國外學(xué)者Valko和Lee［15］對參數(shù)n，τ的獲取方法有較為全面的解釋。首先，將生產(chǎn)時間等分為3個階段，分別統(tǒng)計每一階段的累計產(chǎn)量Q1，Q2和Q3；然后，分別求取后2個階段與第1階段累計產(chǎn)量的比值r21和r31，即

最后，通過求解非線性方程組（式（3）、式（4）），獲取參數(shù)n和τ。

獲得各井的參數(shù)n后，用P指數(shù)概率準(zhǔn)則（P10-P50-P90）進(jìn)行概率計算，構(gòu)建累計分布函數(shù)；利用蒙特卡洛隨機模擬方法，求取各井組中n的高、中、低值nP90，nP50，nP10，將這3個值分別代入式（3）或式（4），即可求取相應(yīng)的τ的高、中、低值τP90，τP50，τP10。

4）單井產(chǎn)量遞減模型的建立。參數(shù)nP90，τP90代表整個井組的產(chǎn)量遞減屬性，適用于井組內(nèi)的任一口單井。因此，利用井組的共同參數(shù)nP90和τP90，結(jié)合單井的最大月產(chǎn)量q0，即可得到該井的相應(yīng)模型。

3.2 實例對比

借用美國德克薩斯州中北部 Fort Worth盆地Barnett頁巖氣藏8口井（包括6口水平井、2口直井），以及路易斯安那州Haynesville頁巖氣藏數(shù)口井的生產(chǎn)數(shù)據(jù)，分別應(yīng)用擴展指數(shù)遞減模型和Arps遞減模型進(jìn)行產(chǎn)能預(yù)測，得到擴展指數(shù)遞減模型的參數(shù)n，τ分別為0.294和3.044，Arps遞減模型的參數(shù)b，Di分別為1.350和0.420。

將其中1口水平井在生產(chǎn)3，6，9 a時的預(yù)測產(chǎn)量與實際產(chǎn)量分別進(jìn)行對比（見圖4），可以看出：Arps遞減模型對頁巖氣井的產(chǎn)能評價，結(jié)果過于樂觀，而擴展指數(shù)遞減模型預(yù)測結(jié)果相對更接近氣藏實際開發(fā)狀況；生產(chǎn)時間為3 a時，2種方法計算得到的預(yù)測產(chǎn)量與實際產(chǎn)量的誤差均較大，隨著生產(chǎn)時間的延長，擴展指數(shù)遞減模型的預(yù)測結(jié)果快速趨于實際產(chǎn)量，可見生產(chǎn)歷史長短對產(chǎn)量遞減模型預(yù)測結(jié)果的準(zhǔn)確性具有重要影響。

圖4 遞減模型預(yù)測結(jié)果對比

4 結(jié)論

1）利用擴展指數(shù)遞減模型進(jìn)行頁巖氣井產(chǎn)能預(yù)測，可有效提高生產(chǎn)動態(tài)預(yù)測的準(zhǔn)確性。

2）由于頁巖氣藏地質(zhì)條件的不確定性和開發(fā)過程中流動形態(tài)的復(fù)雜性，以及水平井增產(chǎn)措施的影響，難以確定頁巖氣井生產(chǎn)初期的遞減規(guī)律，因此，需要選擇生產(chǎn)歷史較長的井作為研究對象。

3）擴展指數(shù)遞減模型通過對井進(jìn)行分組，利用井組整體參數(shù)，結(jié)合單井最大月產(chǎn)量，對單井產(chǎn)能進(jìn)行動態(tài)預(yù)測。隨著頁巖氣井?dāng)?shù)目不斷增加，模型獲取的單井生產(chǎn)動態(tài)將更為準(zhǔn)確。

5 符號注釋

q為產(chǎn)量，m3/月；q0為初始產(chǎn)量，m3/月；t為生產(chǎn)時間，月；n，τ分別為擴展指數(shù)遞減模型的遞減指數(shù)和時間參數(shù)；b，Di分別為Arps遞減模型中的遞減指數(shù)和初始遞減率；Q為累計產(chǎn)量，m3；Np，∞為最終可采儲量，m3；RF為剩余采收率；P10，P50，P90分別為小于10%，50%，90%的概率分布。

［1］ 孫贊東，賈承造，李相方，等.非常規(guī)油氣勘探與開發(fā)（下冊）［M］.北京：石油工業(yè)出版社，2011：1048-1053.

［2］ 朱豫川，劉建儀，張廣東，等.氣井產(chǎn)量遞減規(guī)律分析方法對比分析［J］.天然氣勘探與開發(fā)，2009，32（1）：28-31.

［3］ 鐘儀華，張志銀，朱海雙.特高含水期油田產(chǎn)量預(yù)測新方法［J］.斷塊油氣田，2011，18（5）：641-644.

［4］ Baihly J，Altman R，Malpani R，et al.Shale gas production decline trend comparison over timeand basins［R］.SPE 135555，2010.

［5］ 梁波，陳友蓮，楊洪志.現(xiàn)代產(chǎn)量預(yù)測模型應(yīng)用研究［J］.天然氣勘探與開發(fā)，2004，27（2）：54-56．

［6］ 張順利，劉秀婷，王勝義，等.一種產(chǎn)量預(yù)測新方法的研究及探討［J］.斷塊油氣田，2006，13（2）：51-53，92.

［7］ Jhon L P.Gas Reservoir Engineering［M］.郭晶晶，王鳳江，趙玉龍，等譯.北京：石油工業(yè)出版社，2011：254-270.

［8］ Arps J J.Analysis of decline curve［J］.Trans.A.I.M.E.，1945（160）：228-247.

［9］ 王南，鐘太賢，劉興元，等.復(fù)雜條件下頁巖氣藏生產(chǎn)特征及規(guī)律［J］.斷塊油氣田，2012，19（6）：767-770.

［10］姜寶益，李治平，巨亞鋒，等.修正的產(chǎn)量不穩(wěn)定法預(yù)測頁巖氣動態(tài)儲量［J］.石油鉆探技術(shù)，2012，40（2）：66-69.

［11］李亞洲，李勇明，羅攀，等.頁巖氣滲流機理與產(chǎn)能研究［J］.斷塊油氣田，2013，20（2）：186-190.

［12］ValkóPeter P.Assigning value to stimulation in the Barnett Shale：A simultaneous analysis of 7，000 plus production histories and well completion records［R］.SPE 119369，2009.

［13］Mattar L，Gault B，Morad K，etal.Production analysisand forecasting ofshalegas reservoirs：Case history-based approach［R］.SPE 119897，2008.

［14］Can B，Kabir C S.Probabilistic performance forecasting for unconventional reservoirs with stretched exponential production declinemodel［R］.SPE 143666，2012.

［15］ValkóP P，John L W.A better way to forecast production from unconventionalgaswells［R］.SPE 134231，2010.

（編輯 王淑玉）

Stretched exponential declinemodel for shale gaswell

Lei Danfeng1，2，3,W ang Li1，2，3,Zhang Xiaowei1，2，3,W ang Nan1，2，3,Yu Rongze1，2，3,Xu Fengting4
(1.Langfang Branch,Research Institute of Petroleum Exp loration&Development,PetroChina,Langfang 065007,China;2.National Energy Shale Gas R&D(experim ent)Center,Langfang 065007,China;3.Key Laboratory of Unconventional Oil&Gas, PetroChina,Langfang 065007,China;4.Gaosheng Oil Production Plant,Liaohe Oilfield Com pany,PetroChina,Panjin 124125,China)

The permeability of shale gas reservoir is very low.The prediction resultwith traditional production declinemodel is too optimistic for shale gaswell.The stretched exponential declinemodel is developed to seek for an suitablemethod for productivity evaluation of shale gaswell based on Arps production decline analysismethod.In the process ofmodel application,the traditional production decline analysismethod is combined with probabilistic predictionmodel.Based on the parameters calculation of single welldecline curve,the same kind wellsaremarked offaccording to initial decline ratio of singlewell.The representative parameter valuesofwell group decline curve are determined with P index probability criterion.And then,the probability prediction for single well ismade with these values.Thismethod evades the influence of shale gas reservoir complexity and production increasing measurementson production decline analysis.The reliability of dynamic prediction on production is improved obviously.Case study is performed with practical production data of shale gaswells in USA.The accuracy and feasibility of thismethod used in shale gas reservoir is further verified.

shale gas reservoir;production forecast;Arps declinemodel;stretched exponentialdeclinemodel

國家科技重大專項“頁巖氣勘探開發(fā)關(guān)鍵技術(shù)研究”子課題“頁巖氣開發(fā)機理及技術(shù)政策研究”（2011ZX05018-005）

TE37

A

10.6056/dkyqt201401016

2013-06-13；改回日期：2013-11-28。

雷丹鳳，女，1988年生，工程師，2011年本科畢業(yè)于西安石油大學(xué)石油工程專業(yè)，主要從事頁巖氣井產(chǎn)能評價工作。E-mail：danf_lei1117@sina.com。

雷丹鳳，王莉，張曉偉，等.頁巖氣井?dāng)U展指數(shù)遞減模型研究［J］.斷塊油氣田，2014，21（1）：66-68，82.

Lei Danfeng，Wang Li，Zhang Xiaowei，et al.Stretched exponential declinemodel for shale gas well［J］.Fault-Block Oil&Gas Field，2014，21（1）：66-68，82.