壓裂參數對頁巖氣井產量遞減典型曲線影響分析

白玉湖 楊 皓 陳桂華 馮汝勇 徐兵祥

1.中海油研究總院，北京 100027；2.新興重工有限公司，北京 100070

0 前言

長水平井多級壓裂技術是頁巖氣得以有效開發的關鍵技術之一。 但正是由于這種技術對儲層的強烈改造，使產量遞減預測非常困難。 前人已經開展了包括理論分析方法［1-2］、數值模擬方法［3-4］以及典型曲線方法［5］等大量的預測頁巖氣產量的探索工作。 各種方法中基于生產動態數據的典型曲線方法是應用最廣泛的一種預測頁巖氣產量遞減規律的方法。 本文以美國某一頁巖油氣區塊27 口頁巖氣井生產動態數據為基礎，分析了典型曲線關鍵參數和壓裂參數之間的關系，以期明晰壓裂參數對頁巖氣井產量遞減規律的影響。

1 區塊內典型曲線影響因素

頁巖油氣由于存儲在以納米尺度為主的復雜孔隙空間中，存儲方式除自由氣之外還具備一定的吸附氣含量。 頁巖油氣的成功開發需要對頁巖儲層進行大規模的長水井多級壓裂，盡最大可能使得頁巖油氣暴露、溝通，增大泄油泄氣面積。 頁巖油氣在納米尺度孔隙中微觀流動機理目前尚不清楚，加之大型壓裂導致的儲層破裂使得頁巖儲層非均質性非常嚴重，存在著大規模的跨尺度流動， 這些都使精確計算頁巖油氣的產能變得非常困難。 因此，基于生產動態數據的典型曲線方法是目前估算頁巖油氣產能的重要手段之一，也是應用最廣泛的手段。

目前，描述頁巖油氣產量遞減的典型曲線主要是雙曲遞減曲線：

式中：qi為初始產量，m3/d；n 為遞減指數，無量綱；D 為遞減率，d-1；t 為時間，d。

可見，典型曲線有三個關鍵參數，即初始產量、遞減率及遞減指數，但影響典型曲線參數的關鍵因素分析未見報道，而這恰恰是頁巖氣產量預測、壓裂參數優化的重要研究內容之一。 為了分析壓裂參數對典型曲線參數的影響， 本文選取美國某一頁巖氣區塊作為分析對象，由于該區塊地質情況比較簡單， 地質參數差別不大，可以排除地質因素對產量的影響，把生產井之間產量變化歸結為壓裂參數的不同。 選取該區塊內生產歷史超過1a的頁巖氣井27 口，對每一口井進行典型曲線預測，獲取每口井的典型曲線參數，即初始產氣量、遞減率及遞減指數，然后分析重要的完井壓裂參數，如水平段長度、壓裂級數、總支撐劑用量、每級支撐劑用量、總壓裂液用量、每級壓裂液用量、總射孔數、總簇數等對典型曲線的初始產氣量、遞減指數、遞減率影響。

2 初始產氣量和壓裂參數的關系

初始產氣量是表征頁巖氣產能的重要參數，一般而言，初始產氣量越高，頁巖氣井的產能越大。 針對本研究區塊， 圖1～10 給出了初始產氣量和壓裂參數之間的關系， 其中圖1～6 分別給出了初始產氣量和水平段長度、壓裂級數、總射孔數、總簇數、總/每級支撐劑用量和總/每級壓裂液用量等之間的關系。 從圖1 可見， 隨水平段長度增加，初始產氣量增加，但當水平段增加到2 300 m時，初始產氣量隨水平段的增加出現降低趨勢，因此，針對該研究區塊，最優的水平段長度應為2 300 m，究其原因，本區塊頁巖氣中凝析油含量較高，隨著水平井段的增加、井筒摩阻損失增加而導致初始產氣量降低。從圖2可見，在壓裂級數低于20 級時，初始產氣量隨壓裂級數增加而增加，而當壓裂級數高于20 級時，初始產氣量增加緩慢。 因此，本區塊最優的壓裂級數為20 級。 從圖3～6 可見，初始產氣量隨著總射孔數、總簇數、總支撐劑和總壓裂液用量的增加而增加，但增加幅度有所降低。圖7給出了每簇的初始產氣量和總簇數之間的關系， 可見，兩者之間呈較好的線性關系，隨著簇數的增加，每簇對初始產氣量的貢獻逐漸降低。 圖8 給出了每米初始產氣量和水平段之間的關系，可見兩者之間具有較好的線性關系， 隨著水平段的增加每米的貢獻量逐漸降低。 圖9給出了每級產氣量和每級支撐劑用量關系，可見，隨著每級支撐劑用量的增加， 其對初始產氣量的貢獻率降低。 圖10 給出了每級壓裂液用量和每級產氣量之間的關系，可見，隨著每級壓裂液用量的增加，其對初始產氣量的貢獻率增加。 因此，從圖1～10 可得出：針對本區塊，在地質特征參數差別不大的前提下，水平段長度、壓裂級數、射孔數、簇數、支撐劑和壓裂液用量等都有一個最優值，當然，水平段長度、壓裂級數、射孔數、簇數這幾個參數密切相關，根據現有工藝條件，水平段長度確定之后，即可大致確定壓裂級數，而每級之內的簇數和每簇內的射孔數則可在一定范圍內進行優化，這些參數的優化，可以結合每級、每簇、單位支撐劑用量和壓裂液用量對初始產氣量的貢獻來確定。

圖1 初始產氣量和水平段長度關系

圖2 初始產氣量和壓裂級數關系

圖3 初始產氣量和總射孔數關系

圖4 初始產氣量和總簇數關系

圖5 初始產氣量和支撐劑總量關系

圖6 初始產氣量和總壓裂液用量關系

圖7 每簇初始產氣量和總簇數關系

圖8 每米初始產氣量和水平段長度關系

圖9 每級初始產氣量和每級支撐劑量關系

圖10 每級初始產氣量和每級壓裂液用量關系

3 遞減指數和壓裂參數的關系

遞減指數表示頁巖油氣在經歷快速遞減之后， 在保持較穩定遞減程度時所對應的產量高低，遞減指數越大，穩定遞減階段的產量就相應越高。 為了分析影響遞減指數的主要因素，圖11～18 給出了遞減指數和水平段長度、壓裂級數、總射孔數、總簇數、總/每級支撐劑用量、總/每級壓裂液用量之間的關系。從圖11～18 可見，遞減指數隨著這些參數的增加均呈現增加的趨勢， 這是因為隨著水平段的增加，壓裂級數相應增加，則頁巖儲層被改造的范圍也就越大，在一定程度上增加了泄油泄氣面積，保證了頁巖油氣流動的物質基礎。總簇數和總射孔數的增加，意味著產生裂縫的幾率增加，在地層中會產生更多的裂縫，孔隙及微裂縫之間的連通性增加，油氣流動能力增強。總支撐劑用量、總壓裂液用量以及每級支撐劑用量、壓裂液用量的增加，意味著壓裂時頁巖中的造縫能力增加，且生產過程中保持裂縫有效開啟的能力及有效泄油泄氣面積增加。泄油泄氣面積的增加，改造程度的增加以及保持裂縫有效能力的增強都會增加油氣的供給能力， 保持較穩定的生產能力，從而使遞減指數增大。

圖11 遞減指數和水平段長度關系

圖12 遞減指數和壓裂級數關系

圖13 遞減指數和總射孔數關系

圖14 遞減指數和總簇數關系

圖15 遞減指數和總支撐劑用量關系

圖16 遞減指數和總壓裂液用量關系

圖17 遞減指數和每級支撐劑用量關系

圖18 遞減指數和每級壓裂液用量關系

4 遞減率和壓裂參數的關系

遞減率是表示在初始生產時刻頁巖油氣產量的遞減快慢，遞減率越大，表示初始時刻產量遞減越快。 為了分析影響遞減率的主要因素， 圖19～26 給出了遞減率和水平段長度、壓裂級數、總射孔數、總簇數、總/每級支撐劑用量、總/每級壓裂液用量之間的關系。 從圖19～26 可見，隨著水平段的增加，遞減率基本上呈現先快速降低，然后降低程度減慢，甚至基本不再降低。隨著壓裂級數的增加遞減率快速降低，這是因為水平段長度增加，壓裂級數隨之增加，泄油泄氣體積增大，從而可有效地延緩遞減。 隨著總射孔數和總簇數的增加，遞減率呈現下降的趨勢，但隨著總簇數的增加，遞減率的下降幅度更大，由此說明，相比總射孔數而言，總簇數是影響遞減率更大的參數，這是因為增加簇數能有效地促進裂縫的形成和溝通， 有效降低遞減率。 隨著總壓裂液用量和每級壓裂液用量的增加，遞減率呈現先快速下降，然后下降趨勢減緩并趨于相對穩定的數值。 遞減率隨著總支撐劑用量和每級支撐劑用量的增加呈現遞減趨勢， 這是因為壓裂液用量和支撐劑用量的增加能夠增加有效裂縫的范圍， 并增加裂縫的導流能力，有效降低遞減率。 由此可見，從遞減率角度而言，壓裂參數均存在最優值。

圖19 遞減率和水平段長度關系

圖20 遞減率和壓裂級數關系

圖21 遞減率和總射孔數關系

圖22 遞減率和總簇數關系

圖23 遞減率和總壓裂液用量關系

圖24 遞減率和總支撐劑用量關系

圖25 遞減率和每級壓裂液用量關系

圖26 遞減率和每級支撐劑用量關系

5 結論

a）初始產氣量隨著水平段長度、壓裂級數、壓裂液用量、支撐劑用量、總射孔數的增加而增加，但增加幅度逐漸減小。

b）每簇、每米、每級初始產氣量隨簇數、水平段長度、每級支撐劑用量和每級壓裂液用量的增加而降低。 遞減指數隨著水平段長度、壓裂級數、總簇數、總射孔數、總/每級支撐劑用量和壓裂液用量的增加變大。

c）遞減率隨水平段長度、壓裂級數、總簇數、總射孔數、總/每級壓裂液用量的增加先快速降低后緩慢降低。

［1］ Moghadams S, Mattat L, Pooladi-darvish. Dual Porosity Typecurves for Shale Reservoir ［C］. Paper CSUG/SPE 137535 Presented at the Canadian Unconventional Resources &International Petroleum Conference,19-21 October 2010,Calgary,Alberta,Canada.New York:SPE,2010.

［2］ Seshardri J, Mattar L. Comparison of Power Law and Modified Hyperbolic Decline Methods ［C］. Paper CSUG/SPE 137320 Presented at the Canadian Unconventional Resources &International Petroleum Conference,19-21 October 2010,Calgary,Alberta,Canada.New York:SPE,2010.

［3］ Thompson J M, Nobakht M, Anderson D M. Modeling Well Performance Data from Overpressure Shale Gas Reservoirs ［C］.Paper CSUG/SPE 137755 Presented at the Canadian Unconventional Resources&International Petroleum Conference,19-21 October 2010,Calgary,Alberta,Canada.New York:SPE,2010.

［4］ 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］.Paper CSUG/SPE 136696 Presented at the Canadian Unconventional Resources& International Petroleum Conference, 19-21 October 2010,Calgary,Alberta,Canada.New York:SPE,2010.

［5］ Ilk D, Rushing J A, Perego A D, et al. Exponential vs. Hyperbolic Decline in Tight Gas Sands-understanding the Origin and Implication for Reserve Estimates Using Arps`Decline Curves［C］.Paper 116731 Presented at SPE Annual Technical Conference and Exhibition,21-24 September 2008,Denver,Colorado.New York:SPE,2008.