頁巖氣藏采收率計算方法探討

杜殿發(fā)，付金剛，周志海，張 婧，王妍妍

(1.中國石油大學(xué)，山東 青島 266580;2.中國石化石油勘探開發(fā)研究院，北京 100083)

頁巖氣藏采收率計算方法探討

杜殿發(fā)1，付金剛1，周志海1，張 婧1，王妍妍2

(1.中國石油大學(xué)，山東 青島 266580;2.中國石化石油勘探開發(fā)研究院，北京 100083)

中國頁巖氣資源豐富，開發(fā)潛力巨大，但勘探開發(fā)時間較短，現(xiàn)階段頁巖氣采收率資料的獲得難度仍然較大。通過借鑒煤層氣藏和低滲透氣田采收率的標(biāo)定方法，對頁巖氣藏采收率各種計算方法進(jìn)行分析探討，認(rèn)為現(xiàn)有方法在一定程度上均存在其局限性。結(jié)合頁巖氣藏特殊地質(zhì)情況，運用物質(zhì)平衡原理，利用氣體狀態(tài)方程和等溫吸附曲線，對傳統(tǒng)容積法進(jìn)行修正，提出一種適用于頁巖氣藏采收率計算的新方法，達(dá)到了既操作簡單，又能夠?qū)搸r氣藏采收率進(jìn)行精確計算和預(yù)測的目的。

頁巖氣藏;采收率;容積法;等溫吸附;數(shù)值模擬

0 引言

頁巖氣采收率是指在現(xiàn)有經(jīng)濟(jì)技術(shù)條件下，頁巖氣藏開發(fā)過程中可采儲量與原始地質(zhì)儲量的比值。當(dāng)?shù)貙訅毫ο陆抵翉U棄壓力前，其采出儲量即為可采儲量。頁巖氣藏采收率是頁巖氣開發(fā)者投資決策、實施開發(fā)規(guī)劃、調(diào)整開發(fā)方案的重要依據(jù)，是投資者最關(guān)注的指標(biāo)，采收率的高低直接影響著頁巖氣藏的開發(fā)方式和開發(fā)工藝。

中國對頁巖氣的勘探開發(fā)起步較晚，目前頁巖氣采收率資料匱乏。預(yù)測采收率的方法通常借鑒煤層氣藏和低滲透氣藏的采收率標(biāo)定方法，此次研究主要對類比法、解析法、容積法、產(chǎn)量遞減法以及數(shù)值模擬法等常用方法對頁巖氣藏的適用性進(jìn)行了分析對比，并針對頁巖氣藏的特殊性，對容積法(物質(zhì)平衡法)修正優(yōu)化［1-7］，使之能夠更精確地計算頁巖氣藏采收率，使預(yù)測結(jié)果更加合理、科學(xué)和準(zhǔn)確。

1 采收率方法對比

分別對類比法、解析法、容積法、等溫吸附曲線法、產(chǎn)量遞減法和數(shù)值模擬法進(jìn)行對比分析。

(1)類比法是一種相對較簡單的采收率預(yù)測方法，是通過與地質(zhì)條件相似的區(qū)塊進(jìn)行比較，從而得到目標(biāo)區(qū)塊內(nèi)頁巖氣藏的采收率［8］。目前從美國油氣雜志報道的數(shù)據(jù)可以看出，頁巖氣的采收率變化范圍較大(5%～60%)，且初期采收率較低(5%～10%)。但隨著鉆、完井技術(shù)和儲層改造技術(shù)的進(jìn)步，頁巖氣藏的采收率有望進(jìn)一步提高。由于不同頁巖氣藏的地質(zhì)條件各不相同，且開發(fā)方式和技術(shù)手段也不相同，因此，利用類比法確定的采收率局限性較大，預(yù)測的可靠性不高。

(2)解析法是通過對儲層的巖心解吸實驗來確定采收率的方法。理論上認(rèn)為，解吸氣和損失氣在低壓下能夠解吸并沿著裂縫或孔隙發(fā)生運移，最終被開采到地面;相反，殘余氣則認(rèn)為不可解吸也不能被開采出來［2，5］。因此，可采儲量可以認(rèn)為是解吸氣與損失氣之和，而在實際上，解吸氣和損失氣無法完全采出，其采出的程度也受儲層地質(zhì)狀況、生產(chǎn)制度以及工藝技術(shù)等諸多因素的影響。因此，該方法估算出的采收率結(jié)果可靠性較差，只能為勘探開發(fā)初期提供一定的數(shù)據(jù)參考。

(3)容積法又稱物質(zhì)平衡法，是在物質(zhì)平衡方程基礎(chǔ)上進(jìn)行氣藏動態(tài)地質(zhì)儲量計算進(jìn)而預(yù)測采收率的方法［9］。傳統(tǒng)容積法是以定容封閉氣藏的物質(zhì)平衡原理為基礎(chǔ)，近似地計算采收率。根據(jù)氣藏工程理論，推導(dǎo)出氣藏地質(zhì)儲量為:

式中:GR為氣藏地質(zhì)儲量，m3;A為油層面積，m2;h為油層厚度，m;pa為油層廢棄壓力，MPa;pi為油層初始壓力，MPa;psc為標(biāo)準(zhǔn)狀況下壓力，MPa;T為油層溫度，K;Tsc為標(biāo)準(zhǔn)狀況下溫度，K;Sgi為油層初始含氣飽和度;Za為廢棄壓力下氣體壓縮系數(shù); Zi為原始地層壓力下氣體壓縮因子;φ為油層孔隙度。

針對無水侵氣藏，束縛水及巖石膨脹影響一般較小，可以近似地認(rèn)為氣藏在開發(fā)過程中含氣飽和度保持不變，則氣藏采收率ER為:

在用該方法計算頁巖氣藏采收率時，忽略了頁巖氣藏中吸附氣和游離氣不同賦存狀態(tài)的影響，計算結(jié)果存在誤差。

(4)等溫吸附曲線法簡單方便，是煤層氣地質(zhì)工作者比較常用的方法［10］。煤層氣在開采時需要排水降壓，使吸附氣解吸來達(dá)到開采目的，在計算采收率時，可忽略其中所含有的少量游離氣。但由于頁巖氣藏同時存在吸附狀態(tài)的氣體和游離狀態(tài)的氣體，且游離氣的含量往往比較大。在開采過程中，隨著游離態(tài)氣體的采出，可達(dá)到自然降壓的目的，從而使吸附態(tài)氣體解吸出來，達(dá)到開采的目的，因此游離氣對采收率的影響不可忽略。由于等溫吸附方法只考慮了吸附氣的存在卻忽略了游離氣的影響，用該方法計算的采收率往往偏低。

(5)產(chǎn)量遞減法是根據(jù)氣藏實際生產(chǎn)歷史資料的生產(chǎn)規(guī)律和開發(fā)趨勢，利用產(chǎn)量遞減曲線間接估算采收率的方法。該方法是常規(guī)氣藏進(jìn)行采收率計算的常用方法［2］。對于頁巖氣藏來說，氣藏初期的遞減率較大，后期遞減率變化趨于穩(wěn)定。該方法一般應(yīng)用于產(chǎn)氣高峰期后產(chǎn)量遞減并穩(wěn)產(chǎn)時段，其局限性是需要的生產(chǎn)時間較長，同時當(dāng)對儲層進(jìn)行多次改造時，遞減規(guī)律會因儲層改造而發(fā)生改變，此種情況下不適合用該方法來計算采收率。

(6)數(shù)值模擬法是通過歷史擬合，并用氣藏初始條件來預(yù)測未來產(chǎn)量，從而確定頁巖氣藏采收率的方法。數(shù)值模擬方法適用于已經(jīng)勘探開發(fā)過的試驗區(qū)塊，同時擁有較全面的試采數(shù)據(jù)。雖然數(shù)值模擬方法預(yù)測精度較高，但是頁巖氣藏數(shù)值模擬需要大量的儲層參數(shù)和合適的數(shù)模軟件作為支持。目前中國頁巖氣藏勘探開發(fā)起步較晚，規(guī)模化的生產(chǎn)區(qū)塊較少且缺乏有效的生產(chǎn)數(shù)據(jù)，還沒有一種專門針對頁巖氣藏的數(shù)值模擬軟件，只能依賴于常規(guī)的商業(yè)數(shù)模軟件。

通過對各種方法的對比可以看出(表1)，在計算頁巖氣藏采收率時可采用的方法有類比法、產(chǎn)量遞減法和數(shù)值模擬法，而其他方法由于可靠性差不能夠滿足頁巖氣藏采收率計算的需要。

表1 頁巖氣藏采收率計算方法對比

由于現(xiàn)有方法均存在一定的局限性，可考慮通過對現(xiàn)有的容積法進(jìn)行修正，使之更加符合頁巖氣藏的特殊性，從而使頁巖氣藏的采收率預(yù)測結(jié)果更加合理、科學(xué)和準(zhǔn)確。

2 修正的容積法

對容積法的修正借鑒了煤層氣采收率計算中的等溫吸附曲線法［11-14］，同時又考慮到頁巖氣藏開發(fā)中游離氣含量不可忽略的特點。在計算采收率時，基于物質(zhì)平衡原理并將吸附氣和游離氣分開考慮，具體步驟如下。

(1)確定廢棄壓力。氣藏廢棄壓力是指氣井具有工業(yè)開采價值的最低壓力，廢棄壓力受到氣藏地質(zhì)狀況、開采工藝和經(jīng)濟(jì)指標(biāo)等因素影響。由于中國目前沒有枯竭的頁巖氣藏，因此，無法獲得頁巖氣藏的準(zhǔn)確廢棄壓力資料。目前計算氣藏廢棄壓力的方法主要有2種，一種是經(jīng)驗公式法，另一種是穩(wěn)定產(chǎn)能方程法。利用經(jīng)驗公式計算廢棄壓力的方法比較簡單，但考慮因素比較單一。穩(wěn)定產(chǎn)能法計算的氣藏因素較全面，但影響氣藏廢棄壓力的因素很多，同時受動態(tài)資料匱乏的制約。因此，要在深入認(rèn)識氣藏地質(zhì)狀況的基礎(chǔ)上，選擇合適的計算方法來確定廢棄壓力。

(2)計算游離氣含氣量變化值。利用氣體狀態(tài)方程計算單位質(zhì)量的頁巖儲層中游離氣在由原始地層壓力降到廢棄壓力的過程中，由于頁巖氣藏中氣體體積膨脹而開采出的游離氣的體積。

由氣體狀態(tài)方程可以得到:

式中:Vi為原始地層下游離氣含氣量，m3/t;Vme(f)為初始游離氣含氣量，m3/t;Vre(f)為廢棄壓力下游離氣殘余含氣量，m3/t;Ti為原始地層溫度，K;pu為廢棄壓力，MPa;Tu為廢棄壓力下地層溫度，K; Zsc為標(biāo)準(zhǔn)狀況下氣體壓縮因子。

由式(3)、(4)可得:

則游離氣含氣量變化為:

式中:ΔV(f)為單位質(zhì)量頁巖中游離氣含氣量變化，m3/t。

(3)計算吸附氣含氣量變化值。研究表明，頁巖層對氣體的吸附規(guī)律符合Langmuir方程，且頁巖氣解吸和吸附過程是可逆的。根據(jù)頁巖等溫吸附實驗，隨著壓力的增大，頁巖對氣體的吸附能力增強;當(dāng)頁巖吸附一定數(shù)量的氣體后，壓力減小，氣體解吸，吸附量減少。利用實驗得到的等溫吸附數(shù)據(jù)可以繪制頁巖氣藏等溫吸附曲線，通過等溫吸附曲線可計算廢棄壓力下吸附氣的殘余含氣量，進(jìn)而得到吸附氣含氣量的變化。

在等溫吸附曲線上，過廢棄壓力作垂線與吸附等溫線相交于一點，過該點做水平線得殘余吸附氣含氣量(圖1)。

利用等溫吸附曲線法可得到單位質(zhì)量頁巖中吸附氣的含氣量變化:

式中:Vme(a)為原始地層壓力下氣藏中吸附氣含氣量，m3/t;Vre(a)為廢棄壓力下氣藏中殘余吸附氣含氣量，m3/t;ΔV(a)為單位質(zhì)量頁巖中吸附氣含氣量變化，m3/t。

(4)確定采收率。根據(jù)上述計算得到的游離氣和吸附氣的含氣量變化量，可計算頁巖氣藏的采收率ER(%)為:

圖1 頁巖氣藏等溫吸附曲線示意圖

由于該方法基于定容封閉氣藏物質(zhì)平衡原理的基礎(chǔ)上，同時也考慮了頁巖氣藏中氣體的不同賦存狀態(tài)，具有比較強的理論性，估算的采收率結(jié)果可靠性較高。與其他計算方法相比，該方法既適用于未開發(fā)的區(qū)塊，又適用于已經(jīng)勘探開發(fā)的區(qū)塊，同時操作簡單易行，預(yù)測結(jié)果精確度高，是一種比較理想的頁巖氣藏采收率預(yù)測方法。但該方法忽略了地質(zhì)狀況和開發(fā)方式的復(fù)雜性，單獨應(yīng)用該法計算頁巖氣采收率存在一定的誤差，在誤差范圍內(nèi)可綜合其他方法確定最終采收率。

3 實例分析

某頁巖區(qū)塊屬五峰組—龍馬溪組，儲層孔隙度為4.87%，滲透率為0.25×10-3μm2，平均含水飽和度為33.8%，含氣面積為41.50 km2，氣藏總厚度為87 m，頁巖密度為2.571 6×103kg/m3，原始地層壓力為34.7 MPa。由單井含氣量實測結(jié)果可知，目的層總含氣量為0.44～5.19 m3/t，平均為1.97m3/t，主要以損失氣與解吸氣為主，殘余氣含量低。損失氣含量為0.11～3.90 m3/t，平均為1.14 m3/t;解吸氣含量為0.31～1.40 m3/t，平均為0.79 m3/t;殘余氣含量為0.01～0.07 m3/t，平均為0.04 m3/t。

分別采用類比法、產(chǎn)量遞減法、數(shù)值模擬法和修正后的容積法計算采收率(表2)。

通過不同方法計算的采收率結(jié)果來看，利用修正前的容積法所得到的采收率計算結(jié)果明顯偏高，不符合頁巖氣藏開采的實際情況。而采用其他4種方法得到的結(jié)果均在合理范圍之內(nèi)，其中修正后的容積法計算結(jié)果與數(shù)值模擬法得到的結(jié)果偏差較小，說明修正后的容積法可靠性較高。考慮到頁巖氣田開發(fā)過程中的復(fù)雜性及風(fēng)險性，可取類比法、產(chǎn)量遞減法、數(shù)值模擬法和修正后的容積法所得結(jié)果的平均值作為最終該頁巖氣藏的采收率，供氣田開發(fā)設(shè)計參考。

表2 不同方法計算采收率結(jié)果統(tǒng)計

4 結(jié)論及建議

(1)通過分析采收率計算方法可知，現(xiàn)有的采收率計算方法均有其適用條件，其中類比法、產(chǎn)量遞減法和數(shù)值模擬法可對頁巖氣藏采收率計算提供一定的參考。

(2)考慮到現(xiàn)有采收率計算方法的局限性和頁巖氣藏的特殊性，提出基于容積法的修正方法。該方法操作簡單易行且預(yù)測結(jié)果精確度高，是一種比較理想的頁巖氣藏采收率預(yù)測方法。

(3)據(jù)現(xiàn)有資料難以確定其采收率的頁巖氣藏區(qū)塊，可綜合利用類比法、產(chǎn)量遞減法、數(shù)值模擬法和修正的容積法進(jìn)行頁巖氣藏開發(fā)設(shè)計。

(4)目前對于頁巖氣藏的滲流機(jī)理研究已日漸深入，但由于頁巖氣藏滲流機(jī)理的復(fù)雜性，導(dǎo)致滲流機(jī)理研究上仍存在很大爭議。應(yīng)繼續(xù)加深對滲流機(jī)理的研究，提高采收率預(yù)測的準(zhǔn)確性，從而保障頁巖氣藏高效合理地開發(fā)。

［1］郭平，景莎莎，彭彩珍.氣藏提高采收率技術(shù)及其對策［J］.天然氣工業(yè)，2014，34(2):48-55.

［2］周尚忠.當(dāng)前中國煤層氣采收率估算方法及存在問題［J］.中國煤層氣，2011，8(4):12-13.

［3］姜瑞忠.低滲透油藏彈性采收率計算方法研究［J］.特種油氣藏，2008，15(4):70-72.

［4］李士倫.天然氣工程［M］.北京:石油工業(yè)出版社，2000:111-115.

［5］李明宅，孫晗森.煤層氣采收率預(yù)測技術(shù)［J］.天然氣工業(yè)，2008，28(3):25-29.

［6］伍帥.低滲透氣藏采收率計算方法分析［J］.內(nèi)江科技，2012，(2):16.

［7］郝玉鴻，許敏，徐小蓉.正確計算低滲透氣藏的動態(tài)儲量［J］.石油勘探與開發(fā)，2002，29(5):66-67.

［8］黃炳光，劉蜀知，唐海，等.氣藏工程與動態(tài)分析方法［M］.北京:石油工業(yè)出版社，2004:4-5.

［9］李寶國，石萬里，王淑娟，等.氣田采收率計算過程中的問題探討［J］.海洋石油，2006，26(3):55-59.

［10］陳春琳，林大楊.等溫吸附曲線方法在煤層氣可采資源量估算中的應(yīng)用［J］.中國礦業(yè)大學(xué)學(xué)報，2005，34 (5):679-682.

［11］陳薇.幾種求取氣藏廢棄壓力的方法探討［J］.遼寧化工，2013，42(8):934-935.

［12］郝玉鴻，閻紀(jì)輝.確定氣井廢棄地層壓力的新方法［J］.油氣采收率技術(shù)，1999，6(4):77-82.

［13］李忻洪.蘇里格氣田蘇20井區(qū)廢棄地層壓力研究［J］.油氣井測試，2010，19(2):22-24.

［14］盧濤，蘭義飛，何磊，等.一種計算氣藏廢棄壓力的方法——二項式產(chǎn)能方程法［J］.天然氣勘探與開發(fā)，2009，32(2):18-20.

編輯姜 嶺

TE377

A

1006-6535(2015)05-0074-04

20150601;改回日期:20150807

中國石油科技管理部科技創(chuàng)新基金項目“頁巖氣藏滲流機(jī)制及產(chǎn)能評價研究”(2011D-5006-0210)

杜殿發(fā)(1972-)，男，教授，1993畢業(yè)于石油大學(xué)(華東)油藏工程專業(yè)，2002畢業(yè)于該校油氣田開發(fā)工程專業(yè)，獲博士學(xué)位，現(xiàn)從事油氣田開發(fā)方面的教學(xué)和科研工作。

10.3969/j.issn.1006-6535.2015.05.015