余 果 劉海峰 李海濤 鄭 姝 劉林清 敬 雷

1.中國(guó)石油西南油氣田公司規(guī)劃計(jì)劃處 2.中國(guó)石油西南油氣田公司勘探開發(fā)研究院 3.四川川港燃?xì)庥邢挢?zé)任公司

0 引言

石炭系氣藏主要分布于四川盆地東部，區(qū)域構(gòu)造位置隸屬于川東高陡構(gòu)造帶，受印支運(yùn)動(dòng)影響，高陡構(gòu)造成排成帶分布[1]。由于構(gòu)造幅度較高、氣水分異徹底，氣藏的氣水分布多為受構(gòu)造控制的邊水賦存模式[2]。川東石炭系氣藏總體具有埋藏較深、儲(chǔ)層連通性好、地層水活躍、單井控制儲(chǔ)量大、單井穩(wěn)產(chǎn)能力強(qiáng)等特點(diǎn)[3]。

1977年，位于相國(guó)寺構(gòu)造的相18井獲氣，發(fā)現(xiàn)石炭系氣藏。隨后相繼在臥龍河、張家場(chǎng)、福成寨等構(gòu)造鉆獲多個(gè)石炭系氣藏，石炭系成為四川盆地天然氣的主力產(chǎn)層。從1977年到2022年，川東石炭系氣藏已走過(guò)45年的勘探開發(fā)歷程，多數(shù)氣藏生產(chǎn)年限在15年以上，歷經(jīng)了上產(chǎn)、穩(wěn)產(chǎn)、遞減的開發(fā)階段。目前，川東石炭系氣藏累計(jì)探明地質(zhì)儲(chǔ)量約2 400×108m3，氣藏的地質(zhì)儲(chǔ)量采出程度普遍高于50%，處于開發(fā)中后期。

未來(lái)，川東石炭系氣藏仍具有深化勘探、滾動(dòng)擴(kuò)邊和提高采收率的潛力[4]。選取臥龍河、張家場(chǎng)、福成寨、云和寨、五百梯、沙坪場(chǎng)等25個(gè)石炭系氣藏，對(duì)已開發(fā)氣藏的產(chǎn)量開展了風(fēng)險(xiǎn)量化研究，為該區(qū)氣藏調(diào)整開發(fā)技術(shù)對(duì)策和提高采收率提供了理論依據(jù)，對(duì)四川盆地油氣產(chǎn)量風(fēng)險(xiǎn)量化評(píng)價(jià)具有重要的借鑒意義。

行業(yè)內(nèi)產(chǎn)量風(fēng)險(xiǎn)量化評(píng)價(jià)工作長(zhǎng)期采用定性分析方式，而定量化的風(fēng)險(xiǎn)評(píng)價(jià)尚屬技術(shù)盲點(diǎn)。本文采用蒙特卡洛法創(chuàng)新建立川東石炭系產(chǎn)量風(fēng)險(xiǎn)量化模型，通過(guò)篩選關(guān)鍵風(fēng)險(xiǎn)因素，構(gòu)建各因素與產(chǎn)量的概率分布曲線，明確已開發(fā)氣藏產(chǎn)量風(fēng)險(xiǎn)主控因素；建立穩(wěn)產(chǎn)期產(chǎn)量計(jì)算公式作為蒙特卡洛模擬的目標(biāo)函數(shù)，計(jì)算25個(gè)石炭系氣藏穩(wěn)產(chǎn)期年產(chǎn)量最大概率范圍，實(shí)現(xiàn)了風(fēng)險(xiǎn)評(píng)價(jià)從定性到定量研究的突破，可為新區(qū)氣田的產(chǎn)量風(fēng)險(xiǎn)評(píng)價(jià)提供技術(shù)參考。

1 風(fēng)險(xiǎn)量化評(píng)價(jià)方法

常用風(fēng)險(xiǎn)定量評(píng)價(jià)方法有專家調(diào)查法、樹型分析法、模糊綜合法、蒙特卡洛法、組合優(yōu)化法等[5]（表1）。從優(yōu)缺點(diǎn)考慮，5種方法可分為兩類（A類和B類）：A類方法計(jì)算過(guò)程簡(jiǎn)單，但指標(biāo)量化時(shí)人為影響大，爭(zhēng)議多，使用這類評(píng)價(jià)方法必須具備一定的風(fēng)險(xiǎn)評(píng)價(jià)經(jīng)驗(yàn)；B類方法計(jì)算過(guò)程復(fù)雜，需要較多數(shù)據(jù)，但方法的物理意義明確，能直接回答實(shí)現(xiàn)目標(biāo)的概率，便于接受，例如蒙特卡洛法[6]。

表1 五種風(fēng)險(xiǎn)量化評(píng)價(jià)方法對(duì)比表

比較以上5種方法，蒙特卡洛法在精細(xì)量化風(fēng)險(xiǎn)評(píng)價(jià)中占統(tǒng)治性地位。它的應(yīng)用條件是“假設(shè)未來(lái)與過(guò)去具有相同的統(tǒng)計(jì)規(guī)律”。在此前提下，通過(guò)對(duì)歷史數(shù)據(jù)的統(tǒng)計(jì)規(guī)律分析，從而預(yù)測(cè)未來(lái)的可能情況。專家調(diào)查法、模糊綜合法這兩種方法與蒙特卡洛法不同，它們主要依靠人力來(lái)總結(jié)事物規(guī)律、劃分風(fēng)險(xiǎn)等級(jí)標(biāo)準(zhǔn)、評(píng)價(jià)判斷風(fēng)險(xiǎn)大小[7]。樹型分析法并不算獨(dú)立的一種方法，它常常與蒙特卡洛法同時(shí)使用，進(jìn)行方案決策[8]。組合優(yōu)化法主要應(yīng)用于篩選多方案最優(yōu)化組合，而每套方案是相對(duì)確定或無(wú)風(fēng)險(xiǎn)的，即該方法不對(duì)底層風(fēng)險(xiǎn)進(jìn)行分析，而是轉(zhuǎn)為對(duì)確定問(wèn)題（方案）的組合風(fēng)險(xiǎn)研究[9]。

由于氣藏產(chǎn)量風(fēng)險(xiǎn)量化需要綜合評(píng)價(jià)較多的風(fēng)險(xiǎn)因素，那么需要對(duì)每個(gè)風(fēng)險(xiǎn)因素進(jìn)行充分的數(shù)值分析，尋找其數(shù)學(xué)規(guī)律。依據(jù)每個(gè)風(fēng)險(xiǎn)因素的數(shù)值變化規(guī)律，用隨機(jī)數(shù)產(chǎn)生器來(lái)生成大量具有相同分布概率的風(fēng)險(xiǎn)因素?cái)?shù)值。只有通過(guò)這種方式，綜合所有風(fēng)險(xiǎn)因素的概率分布，預(yù)測(cè)天然氣產(chǎn)量的概率分布情況，才能對(duì)產(chǎn)量風(fēng)險(xiǎn)量化系統(tǒng)進(jìn)行客觀、準(zhǔn)確、全面的評(píng)價(jià)。這正是蒙特卡洛法所長(zhǎng)。因此，該方法適用于天然氣產(chǎn)量風(fēng)險(xiǎn)量化研究。

2 蒙特卡洛模擬

蒙特卡洛法，即蒙特卡洛模擬（Monte Carlo Simulation），又稱為隨機(jī)抽樣法，常用于求解某個(gè)目標(biāo)問(wèn)題的概率或期望。該方法的基本思想：依據(jù)已有或建立的目標(biāo)變量的數(shù)學(xué)模型，將具有一定概率分布的自變量參數(shù)代入目標(biāo)變量的數(shù)學(xué)模型，得到目標(biāo)變量的概率分布函數(shù)，從而研究多重參數(shù)對(duì)目標(biāo)變量的影響程度。

當(dāng)隨機(jī)變量f(x)的分布密度函數(shù)為時(shí)，該變量的數(shù)學(xué)期望表述為：

式中E為隨機(jī)變量的分布密度。

同理，對(duì)于多重變量的目標(biāo)函數(shù)，分別根據(jù)各變量的概率分布密度函數(shù)依次隨機(jī)抽取變量值，通過(guò)對(duì)各變量的數(shù)值進(jìn)行大量反復(fù)的獨(dú)立模擬，即可獲得目標(biāo)函數(shù)的概率密度分布情況，蒙特卡洛模擬可以實(shí)現(xiàn)多變量隨機(jī)取樣的計(jì)算過(guò)程[10-12]。

蒙特卡洛模擬的主要步驟如下：①對(duì)需要求解的問(wèn)題建立數(shù)學(xué)模型；②確定數(shù)學(xué)模型中變量的概率分布密度函數(shù)；③依據(jù)概率分布密度函數(shù)，對(duì)模型變量進(jìn)行隨機(jī)抽樣；④將抽取的隨機(jī)變量代入數(shù)學(xué)模型中，求取問(wèn)題的解。

蒙特卡洛模擬的關(guān)鍵是要確定變量的概率分布密度函數(shù)。而確定變量的概率分布密度函數(shù)，需要依據(jù)變量數(shù)值的變化規(guī)律，改進(jìn)已有的概率分布密度函數(shù)，才能更加精確地求取各變量的分布函數(shù)。

蒙特卡洛法中常用的概率分布密度函數(shù)包括Gauss、正弦、多項(xiàng)式、傅里葉等分布函數(shù)。其中，Gauss分布是蒙特卡洛法中較為常用的概率分布函數(shù)。

Gauss模型是基于生命旋回得到的生長(zhǎng)曲線[13]，在油氣領(lǐng)域研究工作中常被用于油氣產(chǎn)量等參數(shù)的概率分布研究[14]。Gauss函數(shù)的形態(tài)對(duì)稱，且曲線趨于平緩[15]。曲線開始先平緩增加，然后迅速上升，再緩慢增加，到達(dá)高峰之后以對(duì)稱的形態(tài)開始下降[16]。依據(jù)峰值數(shù)量的不同，該模型被分為單旋回及多旋回Gauss模型概率分布[17-18]。單旋回Gauss模型概率分布可表示為：

式中N為概率分布；Nm為概率分布峰值；t為概率分布對(duì)應(yīng)的參數(shù)值；tm為概率分布峰值對(duì)應(yīng)的參數(shù)值；s為模型參數(shù)，可以一定程度代表峰的波動(dòng)大小，但不完全等同于標(biāo)準(zhǔn)差。

那么，多旋回Gauss模型概率分布可表示為：

式中k為總旋回?cái)?shù)；i為旋回個(gè)數(shù)，表示某個(gè)旋回?cái)?shù)的參數(shù)下標(biāo)。

在油氣領(lǐng)域的概率分布研究中，多旋回情況出現(xiàn)較多。因此，多旋回Gauss模型概率分布常被用于各種參數(shù)的概率分布研究工作中。

3 石炭系產(chǎn)量風(fēng)險(xiǎn)量化模型

選取臥龍河、張家場(chǎng)、福成寨、云和寨、五百梯、沙坪場(chǎng)等25個(gè)石炭系氣藏，開展蒙特卡洛模擬，計(jì)算風(fēng)險(xiǎn)因素概率分布及產(chǎn)量實(shí)現(xiàn)概率情況。

氣藏開發(fā)階段可依次劃分為上產(chǎn)期、穩(wěn)產(chǎn)期、遞減期3個(gè)階段。其中，上產(chǎn)期是資源勘探程度較低的時(shí)期，對(duì)區(qū)域資源分布情況的認(rèn)識(shí)程度還有待提高，此時(shí)尚不能發(fā)揮勘探區(qū)域的最大開發(fā)潛力[19]；遞減期是指從產(chǎn)量開始遞減至遞減到開發(fā)方案設(shè)計(jì)規(guī)模20%的階段，代表該區(qū)域生產(chǎn)能力受采出程度增加和地層能量下降等因素影響，產(chǎn)量逐年下降；而居于上產(chǎn)期之后、遞減期之前的穩(wěn)產(chǎn)期，是指從產(chǎn)量達(dá)到開發(fā)方案設(shè)計(jì)規(guī)模并穩(wěn)定生產(chǎn)的階段，主要通過(guò)合理配產(chǎn)、加密井網(wǎng)、老井挖潛、區(qū)塊和井間接替、立體開發(fā)、適時(shí)增壓、排水采氣、防砂、治水等工作保障并延長(zhǎng)穩(wěn)產(chǎn)期，穩(wěn)產(chǎn)期末可采儲(chǔ)量采出程度一般在50%～60%。穩(wěn)產(chǎn)期代表資源開發(fā)較為穩(wěn)定的時(shí)期，因此，選擇穩(wěn)產(chǎn)期的年均產(chǎn)氣量作為產(chǎn)量風(fēng)險(xiǎn)量化的目標(biāo)函數(shù)。

3.1 產(chǎn)量概率模型構(gòu)建

通過(guò)層次分析法分析并選取川東石炭系天然氣產(chǎn)量目標(biāo)風(fēng)險(xiǎn)因素。建立石炭系氣藏產(chǎn)量風(fēng)險(xiǎn)量化層次分析圖（圖1），分為以下3個(gè)層次。

圖1 石炭系氣藏產(chǎn)量風(fēng)險(xiǎn)量化層次分析圖

第一層：目標(biāo)層。確定“石炭系氣藏穩(wěn)產(chǎn)期的年均產(chǎn)氣量”總目標(biāo)。

第二層：準(zhǔn)則層。將天然氣產(chǎn)量風(fēng)險(xiǎn)分為4個(gè)大類風(fēng)險(xiǎn)：資源規(guī)模風(fēng)險(xiǎn)、開發(fā)地質(zhì)風(fēng)險(xiǎn)、規(guī)劃部署風(fēng)險(xiǎn)、技術(shù)水平風(fēng)險(xiǎn)[20]。其中，資源規(guī)模風(fēng)險(xiǎn)、開發(fā)地質(zhì)風(fēng)險(xiǎn)是客觀風(fēng)險(xiǎn)因素，只能認(rèn)識(shí)，不能人為改變；規(guī)劃部署風(fēng)險(xiǎn)、技術(shù)水平風(fēng)險(xiǎn)是決策風(fēng)險(xiǎn)因素，既能夠認(rèn)識(shí)也能人為改變。在產(chǎn)量風(fēng)險(xiǎn)量化研究中，往往選取決策風(fēng)險(xiǎn)因素進(jìn)行量化分析，依據(jù)結(jié)果改變決策風(fēng)險(xiǎn)因素，從而起到提高氣藏采收率的目的。

第三層：風(fēng)險(xiǎn)指標(biāo)層。依據(jù)石炭系氣藏產(chǎn)量主控因素，將準(zhǔn)則層的大類風(fēng)險(xiǎn)因素細(xì)化成眾多的小類風(fēng)險(xiǎn)因素。由于風(fēng)險(xiǎn)因素?cái)?shù)量較多，圖1中僅列舉14個(gè)風(fēng)險(xiǎn)因素。

在勘探開發(fā)過(guò)程中，規(guī)劃部署受技術(shù)水平風(fēng)險(xiǎn)的影響較大，這兩個(gè)風(fēng)險(xiǎn)的相關(guān)性較高，獨(dú)立性較低。主要考慮從技術(shù)水平風(fēng)險(xiǎn)中選取決策風(fēng)險(xiǎn)因素，穩(wěn)產(chǎn)期采速表征地質(zhì)資源轉(zhuǎn)化為現(xiàn)實(shí)產(chǎn)量的速率；穩(wěn)產(chǎn)年限表征生產(chǎn)規(guī)模的合理性；穩(wěn)產(chǎn)期末采出程度表征氣藏靜態(tài)特征與動(dòng)態(tài)生產(chǎn)組織的耦合關(guān)系；千米井深井均產(chǎn)能是評(píng)價(jià)氣藏效益開發(fā)的主要指標(biāo)。選取技術(shù)水平風(fēng)險(xiǎn)的4個(gè)決策風(fēng)險(xiǎn)指標(biāo)因素建立穩(wěn)產(chǎn)期產(chǎn)量計(jì)算公式，作為蒙特卡洛產(chǎn)量模擬的目標(biāo)函數(shù)，從而研究多重參數(shù)對(duì)穩(wěn)產(chǎn)期產(chǎn)量這個(gè)目標(biāo)變量的影響程度。建立蒙特卡洛產(chǎn)量模擬的目標(biāo)函數(shù)為：

式中Q為穩(wěn)產(chǎn)期年產(chǎn)量，108m3；A為穩(wěn)產(chǎn)期采速，百分?jǐn)?shù)；B為穩(wěn)產(chǎn)年限；C為穩(wěn)產(chǎn)期末采出程度，百分?jǐn)?shù)；D為千米井深井均產(chǎn)能，104m3/（km·口）；J為穩(wěn)產(chǎn)階段生產(chǎn)井?dāng)?shù)，口；L為生產(chǎn)井平均埋深，km。

在建立的產(chǎn)量目標(biāo)函數(shù)中，穩(wěn)產(chǎn)階段生產(chǎn)井?dāng)?shù)（J）、生產(chǎn)井平均埋深（L）為客觀風(fēng)險(xiǎn)因素；其余4個(gè)因素：穩(wěn)產(chǎn)期采速（A）、穩(wěn)產(chǎn)年限（B）、穩(wěn)產(chǎn)期末采出程度（C）、千米井深井均產(chǎn)能（D）均為決策風(fēng)險(xiǎn)因素，可以調(diào)整參數(shù)取值，從而改變穩(wěn)產(chǎn)期產(chǎn)量規(guī)模。因此，選取技術(shù)水平風(fēng)險(xiǎn)中的4個(gè)決策風(fēng)險(xiǎn)因素，開展基于蒙特卡洛法的產(chǎn)量風(fēng)險(xiǎn)模擬研究。

3.2 風(fēng)險(xiǎn)因素概率分布

蒙特卡洛法產(chǎn)量精準(zhǔn)計(jì)算的關(guān)鍵在于準(zhǔn)確地選取參數(shù)的統(tǒng)計(jì)模型，確定各參數(shù)的概率分布曲線。選取Gauss、正弦、多項(xiàng)式、傅里葉4種概率分布函數(shù)，對(duì)25個(gè)石炭系氣藏4個(gè)決策風(fēng)險(xiǎn)因素的原始數(shù)據(jù)樣本點(diǎn)進(jìn)行蒙特卡洛概率分布計(jì)算，得到4種連續(xù)性的概率分布圖（圖2），分別對(duì)每種風(fēng)險(xiǎn)因素的4種計(jì)算結(jié)果和原始數(shù)據(jù)進(jìn)行相關(guān)性分析，并利用相關(guān)性分析選取最佳的概率分布類型。

圖2 4個(gè)決策風(fēng)險(xiǎn)因素概率分布模擬圖

結(jié)果表明，Gauss概率分布函數(shù)計(jì)算結(jié)果與原始數(shù)據(jù)的相關(guān)性最高（表2），且原始數(shù)據(jù)的概率分布為多旋回。

表2 概率分布函數(shù)相關(guān)性分析結(jié)果表

因此，根據(jù)生產(chǎn)實(shí)際數(shù)據(jù)的統(tǒng)計(jì)規(guī)律，利用多旋回Gauss概率分布密度函數(shù)對(duì)川東石炭系產(chǎn)量風(fēng)險(xiǎn)因素的樣點(diǎn)數(shù)進(jìn)行特征分析及概率分布擬合，計(jì)算出穩(wěn)產(chǎn)期采速（A）、穩(wěn)產(chǎn)年限（B）、穩(wěn)產(chǎn)期末采出程度（C）、千米井深井均產(chǎn)能（D）的分布概率密度函數(shù)，如式（6）～式（9）所示，并得到4個(gè)產(chǎn)量決策風(fēng)險(xiǎn)因素的分布概率與累積概率圖（圖3）。從式（6）～式（9）及圖3可以看出，不同產(chǎn)量風(fēng)險(xiǎn)因素的分布概率密度方程擁有不同的多旋回?cái)?shù)。其曲線為累積概率密度函數(shù)，由分布密度函數(shù)一次累加得到，即

圖3 產(chǎn)量風(fēng)險(xiǎn)因素多旋回Gauss概率分布擬合圖

3.3 穩(wěn)產(chǎn)期年產(chǎn)量概率計(jì)算

設(shè)定模擬次數(shù)10 000次，根據(jù)產(chǎn)量風(fēng)險(xiǎn)因素的分布概率密度函數(shù)循環(huán)抽取隨機(jī)數(shù)，代入產(chǎn)量計(jì)算方程中，計(jì)算25個(gè)川東石炭系氣藏的穩(wěn)產(chǎn)期年均產(chǎn)量，模擬得到穩(wěn)產(chǎn)期年均產(chǎn)量概率分布曲線（圖4）。

圖4 25個(gè)石炭系氣藏穩(wěn)產(chǎn)期年均產(chǎn)量蒙特卡洛模擬結(jié)果圖

模擬結(jié)果表明，穩(wěn)產(chǎn)期年產(chǎn)量概率分布規(guī)律與4個(gè)風(fēng)險(xiǎn)因素概率分布規(guī)律相近，由分布概率與累積概率組成。25個(gè)石炭系氣藏穩(wěn)產(chǎn)期年產(chǎn)量對(duì)應(yīng)的最大概率區(qū)間范圍是43.43×108～126.35×108m3，區(qū)間的最小分布概率是1.80%，最大分布概率為4.56%，平均分布概率為3.65%，累積概率區(qū)間為14.59%～92.88%。

在穩(wěn)產(chǎn)期年均產(chǎn)量概率曲線中，累積概率代表產(chǎn)量大于相應(yīng)規(guī)模的可能性，如“P90=48.1×108m3”表示模型預(yù)測(cè)25個(gè)石炭系氣藏穩(wěn)產(chǎn)期年產(chǎn)量大于48.1×108m3的概率為90%，累積概率P90對(duì)應(yīng)的產(chǎn)量為證實(shí)產(chǎn)量；“P50=86.8×108m3”表示模型預(yù)測(cè)25個(gè)石炭系氣藏穩(wěn)產(chǎn)期年產(chǎn)量大于86.8×108m3的概率為50%，累積概率P50對(duì)應(yīng)的產(chǎn)量為基準(zhǔn)產(chǎn)量。即25個(gè)石炭系氣藏穩(wěn)產(chǎn)期最穩(wěn)妥的產(chǎn)量是48.1×108m3，可實(shí)現(xiàn)的產(chǎn)量是86.8×108m3。從石炭系的生產(chǎn)歷史來(lái)看，雖然氣藏投產(chǎn)時(shí)間不同，但預(yù)測(cè)結(jié)果仍與生產(chǎn)實(shí)際吻合度較高，表明該方法較為可靠。產(chǎn)量模擬結(jié)果可以得到25個(gè)石炭系氣藏穩(wěn)產(chǎn)期年產(chǎn)量實(shí)現(xiàn)的概率，為氣藏開發(fā)技術(shù)對(duì)策調(diào)整和提高采收率提供了理論依據(jù)。

4 結(jié)論及認(rèn)識(shí)

1）深入分析風(fēng)險(xiǎn)量化方法優(yōu)缺點(diǎn)，明確蒙特卡洛法可作為產(chǎn)量風(fēng)險(xiǎn)量化推薦方法。產(chǎn)量目標(biāo)風(fēng)險(xiǎn)評(píng)價(jià)涉及上中下游各個(gè)領(lǐng)域，是一個(gè)龐大復(fù)雜的系統(tǒng)工程，要對(duì)不同的風(fēng)險(xiǎn)因素進(jìn)行客觀、準(zhǔn)確的評(píng)價(jià)，必然要生成及計(jì)算處理大量的風(fēng)險(xiǎn)因素?cái)?shù)據(jù)，蒙特卡洛法比其他方法更為適用。

2）通過(guò)層次分析法篩選出穩(wěn)產(chǎn)期采速、穩(wěn)產(chǎn)年限、穩(wěn)產(chǎn)期末采出程度、千米井深井均產(chǎn)能等4個(gè)決策風(fēng)險(xiǎn)因素。采用多旋回Gauss概率分布密度函數(shù)，構(gòu)建各因素與產(chǎn)量的概率分布曲線，明確已開發(fā)氣藏產(chǎn)量風(fēng)險(xiǎn)主控因素。

3）建立穩(wěn)產(chǎn)期產(chǎn)量計(jì)算公式作為蒙特卡洛模擬的目標(biāo)函數(shù)，計(jì)算得到25個(gè)石炭系氣藏年產(chǎn)量最大概率范圍為43.43×108～126.35×108m3，概率區(qū)間為14.59%～92.88%，實(shí)現(xiàn)了風(fēng)險(xiǎn)評(píng)價(jià)從定性到定量研究的突破。