白玉湖，陳桂華，徐兵祥，陳 嶺.

(中海油研究總院，北京 100028)

?

頁巖油氣產(chǎn)量遞減典型曲線預(yù)測(cè)推薦做法

白玉湖，陳桂華，徐兵祥，陳 嶺.

(中海油研究總院，北京 100028)

針對(duì)頁巖油氣產(chǎn)量遞減典型曲線預(yù)測(cè)因預(yù)測(cè)方法不一致而導(dǎo)致預(yù)測(cè)結(jié)果因人而異的問題，建立預(yù)測(cè)產(chǎn)量遞減典型曲線的推薦做法及基本流程。對(duì)目前典型曲線模型進(jìn)行分析和對(duì)比，給出典型曲線模型選擇建議。本文在典型曲線關(guān)鍵參數(shù)確定方法的基礎(chǔ)上，建立了單井和區(qū)塊確定性典型曲線預(yù)測(cè)流程；針對(duì)單井和區(qū)塊，在有、無生產(chǎn)數(shù)據(jù)條件下，建立了相對(duì)應(yīng)的不確定性典型曲線預(yù)測(cè)方法。本文方法可為頁巖油氣產(chǎn)量遞減典型曲線預(yù)測(cè)提供規(guī)范性做法。

頁巖油氣；產(chǎn)量遞減；典型曲線；推薦做法

頁巖油氣藏作為一種非常規(guī)油氣藏，其儲(chǔ)層非均質(zhì)性極強(qiáng)，具納米級(jí)孔隙尺度、納達(dá)西級(jí)滲透率，因此其商業(yè)性開發(fā)必須依賴于對(duì)頁巖儲(chǔ)層的體積改造。長水平井多級(jí)壓裂技術(shù)是頁巖油氣得以商業(yè)開發(fā)的關(guān)鍵技術(shù)之一，體積改造進(jìn)一步加劇了儲(chǔ)層的非均質(zhì)性，這導(dǎo)致了頁巖油氣產(chǎn)能評(píng)價(jià)技術(shù)與常規(guī)天然氣具有很大的差別。

頁巖氣產(chǎn)能評(píng)價(jià)方法大體上可以分為3類：第一類是基于生產(chǎn)動(dòng)態(tài)數(shù)據(jù)的典型曲線方法[1]；第二類是基于基質(zhì)和裂縫耦合的氣體滲流機(jī)理的簡化解析方法[2-3]；第三類是考慮儲(chǔ)層和流體復(fù)雜因素及滲流、解吸附等機(jī)理的數(shù)值模擬方法[4-7]。目前工程實(shí)踐中應(yīng)用最廣泛的就是典型曲線方法。但在應(yīng)用典型曲線預(yù)測(cè)遞減的過程中仍舊存在一系列問題，缺乏統(tǒng)一的流程和方法，比如與典型曲線相關(guān)的關(guān)鍵參數(shù)的確定方法、典型曲線的預(yù)測(cè)流程、如何進(jìn)行不確定性遞減預(yù)測(cè)等。本文以實(shí)踐為基礎(chǔ)，結(jié)合理論研究成果，給出了在頁巖油氣實(shí)際生產(chǎn)中確定產(chǎn)量遞減典型曲線的推薦做法及基本流程，以期為頁巖油氣產(chǎn)能評(píng)價(jià)工作提供參考和指導(dǎo)。

1 頁巖油氣產(chǎn)量遞減典型曲線模型

1.1 基于Arps模型的典型曲線模型

Arps遞減模型根據(jù)遞減指數(shù)取值不同，可為雙曲遞減、指數(shù)遞減和調(diào)和遞減3種形式。在頁巖油氣遞減預(yù)測(cè)中，為了擬合瞬態(tài)流動(dòng)，雙曲遞減曲線的遞減指數(shù)常常大于1。由于其瞬時(shí)遞減率隨時(shí)間增加而降低，因此會(huì)造成樂觀的后期產(chǎn)量預(yù)測(cè)。

1.2 修改的雙曲遞減模型

為解決瞬時(shí)遞減率隨時(shí)間增加而減小容易導(dǎo)致高估最終可采儲(chǔ)量(EUR)的問題，頁巖油氣常借用Robertson等引入的修改的雙曲遞減方法；但卻帶來了新的困難，即如何確定雙曲遞減何時(shí)轉(zhuǎn)換為指數(shù)遞減，這一問題目前沒有公認(rèn)的方法。

1.3 冪律指數(shù)模型

為解決雙曲遞減后期預(yù)測(cè)偏于樂觀的問題，Ilk等提出了冪律指數(shù)遞減模型，以此作為對(duì)傳統(tǒng)指數(shù)遞減方法的修改[8]；但其參數(shù)較多，且參數(shù)物理意義不明確。

1.4 混合典型曲線模型

Ambrose等假定裂縫是無限導(dǎo)流能力和定井底流壓，基于此提出了混合典型曲線的預(yù)測(cè)方法，即采用分析法和經(jīng)驗(yàn)法的混合[9]。在線性流動(dòng)階段采用分析法，擬壓力差和流量的比值與時(shí)間的平方根在對(duì)數(shù)圖上呈線性關(guān)系。預(yù)測(cè)邊界控制流動(dòng)階段時(shí)采用Arps雙曲遞減。

1.5 Duong 模型

Duong模型是建立在微裂縫發(fā)育的前提條件下，認(rèn)為在生產(chǎn)過程中，隨著頁巖中壓力的降低，由于應(yīng)力影響會(huì)將一部分原先閉合的微裂縫及小斷層激活，從而使裂縫滲透率增加，流動(dòng)能力增加[10]。

此外，還有一些其他的典型曲線模型[11]，但其他模型沒有上述模型應(yīng)用廣泛，在此不作贅述。

針對(duì)不同模型的對(duì)比情況，Marie等對(duì)比了不同典型曲線模型預(yù)測(cè)的EUR，結(jié)果表明，雙曲遞減模型最高，其次為修改的雙曲遞減、Duong模型、冪律指數(shù)模型，最后是指數(shù)模型[12]。Charles等認(rèn)為如果存在流動(dòng)階段變化，冪律指數(shù)模型是一個(gè)較好的選擇，認(rèn)為雙曲遞減預(yù)測(cè)結(jié)果要高于冪律指數(shù)遞減[13]。Seshadri 等認(rèn)為冪律指數(shù)模型和修改的雙曲遞減模型在多數(shù)情況差別不大，但由于冪律指數(shù)模型很復(fù)雜，因此推薦修改的雙曲遞減模型[14]。

筆者在多年的研究工作基礎(chǔ)上，基于對(duì)頁巖油氣生產(chǎn)動(dòng)態(tài)數(shù)據(jù)的認(rèn)識(shí)和分析認(rèn)為，如果頁巖油氣井的生產(chǎn)歷史時(shí)間長、遞減規(guī)律明顯，可以分別采用多個(gè)模型進(jìn)行預(yù)測(cè)、對(duì)比，根據(jù)儲(chǔ)層實(shí)際情況確定推薦的典型曲線模型及其對(duì)應(yīng)的預(yù)測(cè)結(jié)果。如果生產(chǎn)歷史較短，則推薦采用修改的雙曲遞減模型進(jìn)行頁巖油氣典型曲線預(yù)測(cè)。

2 確定性典型曲線預(yù)測(cè)方法

2.1 與典型曲線相關(guān)的關(guān)鍵參數(shù)確定方法

2.1.1 日產(chǎn)數(shù)據(jù)和月產(chǎn)數(shù)據(jù)對(duì)典型曲線的影響

考慮到頁巖油氣井生產(chǎn)動(dòng)態(tài)的震蕩性，采用日產(chǎn)和月產(chǎn)兩種類型生產(chǎn)數(shù)據(jù)所獲取的典型曲線會(huì)有一定差異，推薦采用日產(chǎn)數(shù)據(jù)為基礎(chǔ)。

2.1.2 典型曲線初始產(chǎn)量的確定方法

頁巖油氣初始產(chǎn)量是一個(gè)非常重要的參數(shù)，是表示頁巖油氣井產(chǎn)能的一個(gè)重要參數(shù)。目前如何計(jì)算頁巖油氣井初始產(chǎn)量還沒有統(tǒng)一的計(jì)算標(biāo)準(zhǔn)。推薦確定方法：以初期產(chǎn)量中的最大值為基點(diǎn)，依次選取比最大產(chǎn)量值略小的29個(gè)點(diǎn)，然后對(duì)該30個(gè)值進(jìn)行求和、平均，作為初始日產(chǎn)量。

2.1.3 生產(chǎn)時(shí)率的處理方法

在頁巖油氣井生產(chǎn)過程中，常常由于各種原因?qū)е玛P(guān)井一段時(shí)間后再開井繼續(xù)生產(chǎn)，生產(chǎn)動(dòng)態(tài)曲線上就會(huì)出現(xiàn)產(chǎn)量為零的數(shù)據(jù)點(diǎn)，而動(dòng)態(tài)曲線的形態(tài)也會(huì)發(fā)生變化。推薦處理方法：扣除掉產(chǎn)量為零的點(diǎn)，把關(guān)井時(shí)間段忽略掉，然后把再次開井的時(shí)間點(diǎn)前移至關(guān)井時(shí)間點(diǎn)，從而保證時(shí)間點(diǎn)是連續(xù)的。

2.1.4 生產(chǎn)歷史對(duì)產(chǎn)量遞減典型曲線的影響

研究發(fā)現(xiàn)，在針對(duì)同一口頁巖油氣井進(jìn)行產(chǎn)量遞減分析時(shí)，生產(chǎn)歷史長短會(huì)影響典型曲線，從而最終影響到EUR。建議在對(duì)頁巖油氣區(qū)塊進(jìn)行產(chǎn)量遞減典型曲線分析時(shí)，在選擇有代表性井的基礎(chǔ)上，盡量選擇生產(chǎn)歷史超過一年的井作為研究對(duì)象進(jìn)行預(yù)測(cè)。

2.1.5 壓裂級(jí)數(shù)對(duì)產(chǎn)量遞減典型曲線的影響

不同的頁巖油氣井，壓裂級(jí)數(shù)一般會(huì)有差別。同時(shí)，隨著技術(shù)進(jìn)步，水平段長度和壓裂級(jí)數(shù)也在逐漸增加。目前現(xiàn)場作業(yè)的壓裂級(jí)數(shù)可多達(dá)30～40級(jí)，水平段長度可達(dá)4000 m。當(dāng)一個(gè)區(qū)塊內(nèi)各井壓裂級(jí)數(shù)相差較大時(shí)，推薦采用單級(jí)折算產(chǎn)量或者每單位水平段長度折算產(chǎn)量作為產(chǎn)量遞減典型曲線分析數(shù)據(jù)，所得結(jié)果更能反映典型井的生產(chǎn)特征。

上述5個(gè)方面的詳細(xì)論證見文獻(xiàn)[15]。

2.2 確定性典型曲線的評(píng)價(jià)流程

2.2.1 單井確定性典型曲線評(píng)價(jià)流程

(1)有生產(chǎn)數(shù)據(jù)的單井典型曲線評(píng)價(jià)流程。

對(duì)于有一定生產(chǎn)歷史，且生產(chǎn)數(shù)據(jù)呈現(xiàn)較好遞減規(guī)律時(shí)，可以通過對(duì)生產(chǎn)井?dāng)?shù)據(jù)進(jìn)行擬合而獲得典型曲線。具體流程如下：①選生產(chǎn)歷史長的井：一般而言，如果條件許可，選擇有效生產(chǎn)時(shí)間在一年左右的生產(chǎn)井；如果生產(chǎn)時(shí)間不足，也可以選擇生產(chǎn)時(shí)間超過半年但遞減趨勢(shì)較為明顯的井進(jìn)行分析。②原始生產(chǎn)動(dòng)態(tài)數(shù)據(jù)處理：在選定生產(chǎn)井的基礎(chǔ)上，對(duì)原始生產(chǎn)動(dòng)態(tài)數(shù)據(jù)進(jìn)行篩選，消除明顯偏離數(shù)據(jù)總體遞減趨勢(shì)的異常數(shù)據(jù)點(diǎn)。③處理關(guān)井階段數(shù)據(jù)：去掉關(guān)井時(shí)間段，把關(guān)井后時(shí)間依次前提，保證生產(chǎn)時(shí)間的連續(xù)性。④去掉數(shù)據(jù)異常的點(diǎn)：對(duì)處理關(guān)井階段之后的原始數(shù)據(jù)進(jìn)行分析，消除數(shù)據(jù)明顯偏離總體規(guī)律的數(shù)據(jù)點(diǎn)。⑤動(dòng)態(tài)數(shù)據(jù)規(guī)整化：根據(jù)水平段長度或者壓裂級(jí)數(shù)對(duì)生產(chǎn)數(shù)據(jù)進(jìn)行歸整化，即用產(chǎn)量除以水平段長度或者壓裂級(jí)數(shù)，以消除不同水平段長度或者壓裂級(jí)數(shù)對(duì)預(yù)測(cè)的影響。⑥初始產(chǎn)量計(jì)算：按照推薦的初始產(chǎn)量計(jì)算方法對(duì)初始產(chǎn)量進(jìn)行計(jì)算。⑦典型曲線預(yù)測(cè)：對(duì)處理好的生產(chǎn)動(dòng)態(tài)數(shù)據(jù)進(jìn)行擬合，獲取初始產(chǎn)量、遞減率、遞減指數(shù)、EUR等參數(shù)。

對(duì)于由于生產(chǎn)時(shí)間較短而無法進(jìn)行典型曲線預(yù)測(cè)的，可以采用類比方法，即采用本井計(jì)算出的初始產(chǎn)量，并借鑒其周邊地質(zhì)和壓裂工程參數(shù)相近的井的典型曲線參數(shù)，進(jìn)行典型曲線預(yù)測(cè)。

(2)無生產(chǎn)數(shù)據(jù)的單井典型曲線評(píng)價(jià)流程。

針對(duì)新井，在沒有生產(chǎn)動(dòng)態(tài)數(shù)據(jù)時(shí)，可以采用解析法或數(shù)值模擬方法，利用已有的或者估算的地質(zhì)、油藏參數(shù)，結(jié)合完井壓裂設(shè)計(jì)參數(shù)進(jìn)行典型曲線的初步評(píng)估。由于參數(shù)未經(jīng)過校正，預(yù)測(cè)結(jié)果僅供參考。

2.2.2 區(qū)塊確定性典型曲線評(píng)價(jià)流程

(1)有生產(chǎn)數(shù)據(jù)的區(qū)塊典型曲線評(píng)價(jià)流程。

有時(shí)需要對(duì)一個(gè)區(qū)塊進(jìn)行典型曲線預(yù)測(cè)，獲取該區(qū)塊具有代表性意義的典型曲線。如果該區(qū)塊具有一定數(shù)量的生產(chǎn)井，可以對(duì)這些生產(chǎn)井進(jìn)行分析、預(yù)測(cè)，獲取該區(qū)塊的典型曲線。具體做法是：首先優(yōu)選出該區(qū)塊內(nèi)具有代表性的生產(chǎn)井，按照“2.2.1有生產(chǎn)數(shù)據(jù)的單井典型曲線評(píng)價(jià)流程”對(duì)單井的生產(chǎn)數(shù)據(jù)進(jìn)行處理，完成生產(chǎn)動(dòng)態(tài)數(shù)據(jù)的歸整化處理。然后對(duì)該區(qū)塊內(nèi)的所有生產(chǎn)井歸整化后的動(dòng)態(tài)數(shù)據(jù)進(jìn)行算術(shù)平均，得到一條平均產(chǎn)量曲線，以此為基礎(chǔ)，計(jì)算初始產(chǎn)量，預(yù)測(cè)典型曲線參數(shù)。

(2)無生產(chǎn)數(shù)據(jù)的區(qū)塊典型曲線評(píng)價(jià)流程。

針對(duì)新區(qū)塊，在沒有生產(chǎn)動(dòng)態(tài)數(shù)據(jù)時(shí)，可以采用解析法或數(shù)值模擬方法，利用已有的或者估算的區(qū)塊地質(zhì)、油藏參數(shù)，結(jié)合完井壓裂設(shè)計(jì)參數(shù)進(jìn)行典型曲線的評(píng)估。同樣由于參數(shù)未經(jīng)過校正，預(yù)測(cè)結(jié)果僅供參考。

3 頁巖油氣產(chǎn)量遞減典型曲線的不確定性分析方法

由于頁巖非均質(zhì)性強(qiáng)，單井產(chǎn)量差異也較大，對(duì)于探尋能夠代表一口井或者一個(gè)區(qū)域性的典型曲線，如果采用單一的典型曲線，會(huì)有一定風(fēng)險(xiǎn)，因此，需要采用不確定性分析方法確定頁巖油氣產(chǎn)量遞減典型曲線。

3.1 單井典型曲線的不確定性分析方法

3.1.1 有生產(chǎn)數(shù)據(jù)的單井典型曲線不確定性分析方法

(1)修改雙曲遞減方法。

對(duì)于有一定生產(chǎn)歷史的頁巖油氣井，可以采用修改的雙曲遞減方法進(jìn)行不確定性遞減分析。具體方法：對(duì)已有生產(chǎn)數(shù)據(jù)采用雙曲遞減進(jìn)行擬合，對(duì)指數(shù)部分的遞減率進(jìn)行概率假設(shè)，一般假定為正態(tài)分布，從而獲得指數(shù)遞減率不同概率條件下的典型曲線，獲得P10、P50、P90的典型曲線參數(shù)。圖1為以某一頁巖氣井為例，預(yù)測(cè)的P10、P50、P90的典型曲線與生產(chǎn)數(shù)據(jù)的對(duì)比。

(2)分段預(yù)測(cè)方法。

在頁巖油氣的生產(chǎn)周期內(nèi)，油氣在儲(chǔ)層中的流動(dòng)形態(tài)是不同的。一般而言，線性流動(dòng)階段和擬穩(wěn)態(tài)流動(dòng)階段是兩個(gè)最主要的流動(dòng)階段，據(jù)此可以采用分段典型曲線來預(yù)測(cè)油氣產(chǎn)量，在線性流動(dòng)階段和擬穩(wěn)態(tài)流動(dòng)階段分別采用不同的典型曲線。在線性流動(dòng)階段采用歸整化壓力和時(shí)間平方根的線性函數(shù)來預(yù)測(cè)遞減規(guī)律，對(duì)于擬穩(wěn)態(tài)流動(dòng)階段采用Arps遞減，而且擬穩(wěn)態(tài)流動(dòng)階段也符合Arps遞減，樂觀估計(jì)為調(diào)和遞減，保守估計(jì)為指數(shù)遞減[16]。

圖1 預(yù)測(cè)的P10、P50、P90典型曲線與生產(chǎn)數(shù)據(jù)的對(duì)比Fig.1 Comparison of the predict typical curves of P10, P50, P90 with the production data注：①M(fèi)cfd即千英尺3/日，1 Mcfd=28.32 m3/d。

3.1.2 無生產(chǎn)數(shù)據(jù)的單井典型曲線不確定性分析方法

在對(duì)單個(gè)頁巖油氣井進(jìn)行產(chǎn)量預(yù)測(cè)時(shí)，如果該井沒有生產(chǎn)數(shù)據(jù)，可采用解析法、數(shù)值模擬法等確定的典型曲線預(yù)測(cè)方法，由于沒有參數(shù)的校正，預(yù)測(cè)結(jié)果不確定性較大，因此，需要采用不確定分析方法。根據(jù)對(duì)儲(chǔ)層參數(shù)的認(rèn)識(shí)、井的初步設(shè)計(jì)，獲得該井的儲(chǔ)層厚度、含烴飽和度、基質(zhì)滲透率、流體物性、孔隙度、預(yù)測(cè)壓裂裂縫半長、預(yù)計(jì)的裂縫條數(shù)等參數(shù)的概率分布情況，對(duì)于認(rèn)識(shí)程度較高、把握較大的參數(shù)，可以直接給定該參數(shù)的值；對(duì)于認(rèn)識(shí)程度不高、把握不大的參數(shù)，則給出參數(shù)的范圍及概率分布函數(shù)，然后采用解析方法，獲取該井的產(chǎn)量概率分布。

3.2 區(qū)塊典型曲線的不確定性分析

3.2.1 有生產(chǎn)數(shù)據(jù)的區(qū)塊典型曲線不確定性分析方法

基本思路是：針對(duì)頁巖油氣區(qū)塊，篩選生產(chǎn)動(dòng)態(tài)規(guī)律較好、生產(chǎn)歷史較長的生產(chǎn)井為研究對(duì)象，按照“2.2.1單井確定性典型曲線評(píng)價(jià)流程”進(jìn)行典型曲線預(yù)測(cè)，獲取每口井的初始產(chǎn)量、遞減率、遞減指數(shù)、EUR等典型曲線參數(shù)，從而確定該區(qū)塊每個(gè)典型曲線參數(shù)的概率分布，并對(duì)概率分布進(jìn)行擬合獲取概率分布函數(shù)；然后用蒙特卡洛方法對(duì)典型曲線參數(shù)進(jìn)行隨機(jī)抽樣，利用典型曲線模型進(jìn)行EUR計(jì)算，分析EUR的概率分布，從而獲取各個(gè)概率下的EUR，計(jì)算出 P10、P50、P90的EUR作為該區(qū)塊的推薦EUR。圖2為以某個(gè)具有54口頁巖氣井的區(qū)塊為例，進(jìn)行100000次隨機(jī)抽樣計(jì)算得到的該區(qū)塊的EUR概率分布。

圖2 進(jìn)行100000次隨機(jī)抽樣計(jì)算得到的EUR概率分布Fig.2 The probability distribution of the EUR probability calculated by 100000 random samples

3.2.2 無生產(chǎn)數(shù)據(jù)的區(qū)塊典型曲線不確定性分析方法

在對(duì)頁巖油氣區(qū)塊進(jìn)行評(píng)價(jià)時(shí)，有時(shí)候沒有任何井資料，或者在只有先導(dǎo)井的條件下就需要對(duì)該區(qū)塊的產(chǎn)量進(jìn)行預(yù)測(cè)，為該區(qū)塊的開發(fā)提出相應(yīng)建議。此時(shí)，需要根據(jù)地質(zhì)認(rèn)識(shí)，結(jié)合導(dǎo)眼井測(cè)試化驗(yàn)認(rèn)識(shí)，獲得該區(qū)塊儲(chǔ)層厚度、含烴飽和度、基質(zhì)滲透率、流體物性、孔隙度、預(yù)測(cè)壓裂裂縫半長、預(yù)計(jì)的裂縫條數(shù)等參數(shù)的范圍及概率分布函數(shù)，然后采用解析方法，獲取該區(qū)塊的產(chǎn)量概率分布。在區(qū)塊典型曲線預(yù)測(cè)時(shí)，也可以根據(jù)地質(zhì)、流體的認(rèn)識(shí)，對(duì)區(qū)塊進(jìn)行分區(qū)，針對(duì)每一個(gè)分區(qū)開展上述產(chǎn)量遞減的概率分析。具體流程與方法和“3.1.2無生產(chǎn)數(shù)據(jù)的單井典型曲線不確定性分析方法”一致，此處不再贅述。

4 結(jié)論

得益于技術(shù)進(jìn)步，頁巖油氣作為一種儲(chǔ)量巨大的非常規(guī)油氣資源，在北美已經(jīng)實(shí)現(xiàn)商業(yè)性開發(fā)。在頁巖油氣開發(fā)中，如何清楚地認(rèn)識(shí)產(chǎn)量遞減規(guī)律對(duì)頁巖油氣規(guī)模開發(fā)具有重要的指導(dǎo)意義。

(1)在生產(chǎn)實(shí)踐工作基礎(chǔ)上，本文給出了預(yù)測(cè)頁巖油氣產(chǎn)量遞減典型曲線的推薦做法以及基本流程，以期為頁巖油氣的產(chǎn)量估算提供參考。

(2)建立了確定性典型曲線的評(píng)價(jià)流程，包括典型曲線模型的選擇、與典型曲線相關(guān)的關(guān)鍵參數(shù)的確定方法、單井和區(qū)塊確定性典型曲線評(píng)價(jià)流程等。

(3)建立了在有、無生產(chǎn)動(dòng)態(tài)數(shù)據(jù)條件下，單井和區(qū)塊不確定性典型曲線的預(yù)測(cè)方法，推薦采用P10、P50、P90表示典型曲線的可能情況。

[1] ILK D. Exponential vs hyperbolic decline in tight gas sands-understanding the origin and implication for reserve estimates using Arps decline curves[C]. Denver: SPE Annual Technical Conference and Exhibition, 2008.

[2] LUO S, LANE N, PATHMAN A, et al. Flow region analysis of multi-stage hydraulically-fractured horizontal wells with reciprocal rate derivative function: Bakken case study[C]. CSUG/SPE 137514, 2010.

[3] ORKHAN S, HASAN S A, WATTENBARGER R A. A semi-analytic method for history matching fractured shale gas reservoir[C]. SPE 144583, 2011.

[4] THOMPSON J M, NOBAKHT M, ANDERSON D M. Modeling well performance data from overpressure shale gas reservoirs[C]. CSUG/SPE 137755, 2010.

[5] FAN L, THOMPSON J W, ROBINSON J R. Understanding gas production mechanism and effectiveness of well stimulation in the Haynesville shale through reservoir simulation[C]. CSUG/SPE 136696, 2010.

[6] FREEMAN C M, MORIDIS G, IIK D, et al. A numerical study of transport and storage effects for tight gas and shale gas reservoir systems[C]. SPE 131583, 2010.

[7] MORIDIS G J, BLASINGAME T A, FREEMAN C M. Analysis of mechanisms of flow in fractured tight-gas and shale-gas reservoirs[C]. SPE 139250, 2011.

[8] ROBERTSON S. Generalized hyperbolic equation[C]. Paper SPE 18731, 1988.

[9] AMBROSE R J, CLARKSON C R, YOUNGBLOOD J, et al. Life cycle decline curve estimation for tight/Shale gas reservoirs[C]. SPE 140519, 2011.

[10] DUONG A N. An unconventional rate decline approach for tight and fracture-dominated gas wells[C]. SPE 137748, 2011.

[11] CLARK A J, LAKE L W, PATEK T W. Production forecasting with logistic growth models[C]. SPE 144970, 2011.

[12] MARIE M, RITEJA D, CHRIS B. Comparison of decline curve analysis methods with analytical models in unconventional plays[C]. SPE 166365, 2013.

[13] CHARLES R V. Production decline analysis lessons from classic shale gas wells[C]. SPE 166205, 2013.

[14] SESHADRI J, MATTAR L. Comparison of power law and modified hyperbolic decline methods[C]. CSUG/SPE 137320, 2010.

[15] 白玉湖,楊皓,陳桂華,等.頁巖氣產(chǎn)量遞減典型曲線應(yīng)用分析[J].可再生能源,2013,31(5):115-119.

[16] 白玉湖,王小林,楊皓,等.流動(dòng)階段的頁巖油氣典型曲線分段預(yù)測(cè)方法[J].特種油氣藏,2014,21(6):112-114.

Recommended Method for Production Decline TypicalCurve Prediction of Shale Oil and Gas

Bai Yuhu, Chen Guihua, Xu Bingxiang, Chen Ling

(CNOOC Research Institute, Beijing 100028, China)

For shale oil and gas production decline analysis, the prediction results are usually different for different researchers because the prediction methods differ from each other. To solve this problem, a series of methods and procedures of production decline typical curve analysis were recommended. Based on the comparison and analysis of present typical curve models, the recommended models were proposed. The typical curve prediction procedures for single well and block were established based on the key parameters analysis of typical curve. The uncertainty analysis procedures of typical curve for single well and block were established under the conditions of with and without production data. The proposed method and procedure in this paper can provide a norm for shale oil and gas typical curve prediction.

shale oil and gas; production decline; typical curve; recommended method

中海油綜合科研“海外頁巖油氣產(chǎn)能評(píng)價(jià)技術(shù)與方法研究”(YXKY-2015-ZY-03)資助。

白玉湖(1976—)，男，遼寧鞍山人，2002年畢業(yè)于石油大學(xué)(華東)石油工程專業(yè)，2006年獲中國科學(xué)院力學(xué)研究所流體力學(xué)博士學(xué)位，高級(jí)工程師，主要從事頁巖油氣、天然氣水合物等非常規(guī)資源開發(fā)方面的研究。郵箱：byh_2002@163.com.

TE328

A