張士誠 牟松茹 崔勇

1.中國石油大學(xué)（北京）石油工程教育部重點實驗室 2.中國石油海外勘探開發(fā)公司

頁巖氣壓裂數(shù)值模型分析

張士誠1牟松茹1崔勇2

1.中國石油大學(xué)（北京）石油工程教育部重點實驗室 2.中國石油海外勘探開發(fā)公司

水力壓裂和水平井開采是頁巖氣開發(fā)的主要技術(shù)，在我國尚處在工業(yè)試驗階段，存在很多技術(shù)瓶頸。在總結(jié)分析了頁巖氣壓裂的特點基礎(chǔ)上，探討了網(wǎng)狀裂縫形成的主控因素及裂縫擴(kuò)展模型、產(chǎn)能預(yù)測模型的類型以及優(yōu)缺點。結(jié)果認(rèn)為，特殊的賦存生產(chǎn)機(jī)理、復(fù)雜的裂縫形態(tài)和多尺度的滲流模式是頁巖氣壓裂的主要特點，其目的是形成網(wǎng)狀裂縫，擴(kuò)大儲層改造體積；網(wǎng)狀裂縫的形成主要受天然裂縫與人工裂縫的夾角、水平主應(yīng)力差和巖石的脆性等因素的控制。頁巖氣壓裂產(chǎn)能預(yù)測模型面臨的主要問題是裂縫形態(tài)的模擬和氣體流態(tài)的描述，主要有非常規(guī)裂縫模型、離散裂縫模型和雙重介質(zhì)模型等，這些模型和方法在一定程度上表征了頁巖氣壓裂裂縫形態(tài)和滲流特點，但沒有考慮不規(guī)則的裂縫形態(tài)等。

頁巖氣 開發(fā) 壓裂（巖石） 裂縫擴(kuò)展模型 產(chǎn)能預(yù)測模型 滲流 特點

1 頁巖氣藏的特點

1.1 特殊的賦存生產(chǎn)機(jī)理

頁巖既是烴源巖又是儲集層，就近賦存是頁巖氣成藏的特點。頁巖氣的賦存方式多樣，游離方式、吸附狀態(tài)和溶解狀態(tài)并存?？傮w上主要以游離氣和吸附氣為主，吸附狀態(tài)天然氣的含量變化介于20%～85%。

目前認(rèn)為頁巖氣的產(chǎn)出分為3個階段：①在壓降的作用下，基質(zhì)系統(tǒng)中的頁巖氣在基質(zhì)表面進(jìn)行解吸附；②在濃度差的作用下，頁巖氣由基質(zhì)系統(tǒng)向裂縫系統(tǒng)進(jìn)行擴(kuò)散；③在流動勢的作用下，頁巖氣通過裂縫系統(tǒng)流向生產(chǎn)井筒。由于裂縫空間的有限性，因此早期以游離氣為主的天然氣產(chǎn)量快速下降并且達(dá)到穩(wěn)定，穩(wěn)定期的產(chǎn)量主要是基質(zhì)孔隙里的游離氣和解吸氣。水力壓裂可以增大裂縫空間和連通性，使更多的吸附氣發(fā)生解吸附而向裂縫聚集。

1.2 納米級的微觀孔隙結(jié)構(gòu)

通過掃描電鏡等成像技術(shù)和脈沖法等測試技術(shù)研究表明，納米級的有機(jī)質(zhì)孔隙是頁巖的主要儲集空間［1］和孔隙類型，其形成與分布與有機(jī)質(zhì)的豐度密切相關(guān)；巖心觀察表明天然裂縫較為發(fā)育，但絕大部分被礦物充填處于閉合狀態(tài)?？紫逗秃鸬赖某叽鐬榧{米級別，孔隙、吼道配置關(guān)系復(fù)雜；基質(zhì)滲透率為納達(dá)西級別，孔隙度一般小于7%。

頁巖儲層納米級的微觀孔隙結(jié)構(gòu)，與相同孔隙度的微米級孔隙相比提供了更大的比表面積，為氣體的吸附提供了條件。但是也相應(yīng)引發(fā)如下的問題：①納米級儲層的物性特征參數(shù)難以用常規(guī)的方法測量和計算；②氣體在納米級孔隙中的滲流復(fù)雜多變，流動規(guī)律目前尚不明確；③需采用如水平井多級壓裂等特殊的開發(fā)方式才能獲得經(jīng)濟(jì)產(chǎn)量，且增產(chǎn)的機(jī)理也與常規(guī)壓裂不同。

1.3 水力壓裂形成復(fù)雜的裂縫形態(tài)

常規(guī)壓裂形成的裂縫一般呈雙翼對稱裂縫的形式。但頁巖氣壓裂中微地震監(jiān)測的結(jié)果表明，裂縫的形態(tài)復(fù)雜多變，如圖1所示［2］。1口水平井壓裂后微地震監(jiān)測結(jié)果表明，第1、2段壓裂施工形成了垂直于水平井段的平面縫，第3、4段施工形成了網(wǎng)狀裂縫。目前認(rèn)為頁巖氣壓裂目的，就是要建立一個獨立于傳統(tǒng)意義裂縫半長的更加龐大的裂縫網(wǎng)絡(luò)系統(tǒng)，實現(xiàn)更大規(guī)模的儲層改造波及體積；生產(chǎn)實踐也證明了儲層改造體積越大壓后增產(chǎn)效果越好。

1.4 多尺度的流動狀態(tài)

圖1 頁巖水平井壓裂微地震監(jiān)測結(jié)果圖

頁巖儲層壓裂后形成了多尺度的流動空間，包括納米級的微孔隙和吼道，微米級的大孔隙或天然裂縫以及厘米級別的人工裂縫等，不能用唯一的流動狀態(tài)進(jìn)行描述?？刹捎肒nudsen系數(shù)進(jìn)行流態(tài)的劃分，在不同的流動狀態(tài)下，對應(yīng)以不同的滲流方程進(jìn)行描述。頁巖氣體的流動狀態(tài)主要包括：在有機(jī)質(zhì)表面的解吸附，可采用Langmuir等溫吸附模型進(jìn)行表征；在納米級孔隙中的自由分析運動狀態(tài)，近似地認(rèn)為氣體遵循擴(kuò)散定律運動；在微裂縫中可能出現(xiàn)滑脫現(xiàn)象；在人工裂縫中遵循線性滲流規(guī)律。其中在納米級孔隙中的微觀滲流機(jī)理目前尚不明確，有學(xué)者采用流體力學(xué)格子Boltzmann方法結(jié)合數(shù)字成像技術(shù)從介觀層次的角度進(jìn)行研究，提出了一個新的研究方向。

2 網(wǎng)狀裂縫的形成及主控因素分析

Jon等［3］采用邊界元法對壓裂時多裂縫的同時延伸和它們與天然裂縫之間的相互作用進(jìn)行了研究，認(rèn)為在天然裂縫發(fā)育的條件下，天然裂縫與人工裂縫的夾角、擬凈壓力系數(shù)是影響網(wǎng)狀裂縫形成的主要因素，水平井中人工裂縫和天然裂縫夾角越大，值越大則越容易形成網(wǎng)狀裂縫；直井條件下不容易形成網(wǎng)狀裂縫。

式中Rn為擬凈壓力系數(shù)；pfrac為凈壓力，MPa；Sh，max為水平最大主應(yīng)力，MPa；Sh，min為水平最小主應(yīng)力，MPa。

Gu等［4］建立了天然裂縫和人工裂縫相互作用的判斷準(zhǔn)則，考慮了非正交的模式，采用UFM模型進(jìn)行計算，結(jié)果表明水平面主應(yīng)力差、縫內(nèi)凈壓力、天然裂縫密度以及巖性是影響裂縫形態(tài)的主要因素；水平主應(yīng)力差越小、天然裂縫發(fā)育程度越高，易于形成網(wǎng)狀裂縫；楊氏模量越大，巖石脆性越強(qiáng)，易于形成粗糙節(jié)理并保持裂縫開啟；天然裂縫密度、基質(zhì)滲透率等參數(shù)也對裂縫的形態(tài)和規(guī)模有一定的影響。

生產(chǎn)實踐表明，頁巖中石英的含量越多巖石的脆性越大，越容易形成網(wǎng)狀裂縫，獲得較好的增產(chǎn)效果。提出了脆度指數(shù)的概念，即巖石中石英的含量占全部礦物成分的百分比。在施工設(shè)計中，脆度指數(shù)越大（大于50%），則采用清水壓裂，大液量、少支撐劑量的方法，形成網(wǎng)狀裂縫獲得較好的效果；若脆度指數(shù)小（低于30%），則采用常規(guī)壓裂的方法。

3 復(fù)雜裂縫模型

3.1 線網(wǎng)模型

Xu等［5］提出了表示復(fù)雜裂縫形態(tài)的線網(wǎng)模型（wire－mesh），認(rèn)為頁巖氣藏水平井壓裂產(chǎn)生的裂縫網(wǎng)絡(luò)是沿井筒對稱的橢球體，通過將該橢球體劃分為數(shù)條正交的水平、垂直均勻截面來描述高滲裂縫。線網(wǎng)模型計算時應(yīng)用巖石力學(xué)方法考慮了壓裂過程中裂縫橢球體的實時擴(kuò)展，考慮了施工參數(shù)的影響，并計算了支撐劑在裂縫中的分布情況。線網(wǎng)的不足之處在于：①它必須將油藏改造區(qū)域近似為沿井筒對稱的橢球體，不能模擬不規(guī)則的裂縫形態(tài)；②沒有建立判斷準(zhǔn)則，直接地認(rèn)為天然裂縫與人工裂縫相連接；③沒有考慮人工裂縫之間的相互干擾；④裂縫間距和改造體積由微地震監(jiān)測結(jié)果確定，僅限于本段壓裂施工模擬，計算結(jié)果不具有普遍適用性。

線網(wǎng)模型為半解析模型（圖2），泄流區(qū)域呈橢圓形，有兩組垂直正交的裂縫分別平行于橢圓平面上兩個主應(yīng)力（σh、σH），為避免裂縫間距的主觀性，dx、dy依據(jù)實際測試數(shù)據(jù)確定，橢圓的長短軸a和b以及高度h由微地震結(jié)果給出。

求解方程如下：

圖2 線網(wǎng)模型圖

式中q為壓裂液量，m3／d；ti為壓裂液排量，m3／min；Lxi、Lyj分別為x和y方向裂縫段長度，m；Wxi、Wyi分別為x和y方向裂縫段寬度，m；Kfx、Kfy分別為x和y方向裂縫滲透率，mD；r為橢圓的寬高比；B為橢圓積分。

3.2 非常規(guī)模型

Weng等［6］提出了非常規(guī)裂縫擴(kuò)展模型（UFM－unconventional fracture model）。該模型為數(shù)值模型（圖3），能夠模擬天然裂縫和人工裂縫之間的相互作用，建立了裂縫端部擴(kuò)展準(zhǔn)則，考慮了壓裂液的一維流動、支撐劑的輸送和裂縫寬度的彈性變形，通過計算應(yīng)力陰影考慮相鄰人工裂縫間的相互作用；采用三層模型模擬支撐劑輸送，從底至上依次為支撐劑，混砂液和壓裂液。非常規(guī)模型比線網(wǎng)模型提供了一種更為精確預(yù)測裂縫分布、幾何形態(tài)和支撐劑分布的方法，充分考慮了儲層巖石力學(xué)性質(zhì)以及不規(guī)則的裂縫形態(tài)，有助于更好地了解天然裂縫對于預(yù)想的人工裂縫幾何形狀和對產(chǎn)能的影響，模型的計算結(jié)果可以通過微地震監(jiān)測進(jìn)行校正。非常規(guī)模型的主要問題在于：①天然裂縫的分布依賴于離散裂縫地質(zhì)建模的結(jié)果；②非常規(guī)復(fù)雜裂縫模型對輸入?yún)?shù)的精確性要求較高。

圖3 非常規(guī)復(fù)雜裂縫模型圖

4 產(chǎn)能預(yù)測模型

4.1 離散裂縫模型

離散裂縫模型建立的基礎(chǔ)是Fisher等［7］提出的復(fù)雜裂縫形態(tài)理論，認(rèn)為頁巖氣壓裂后形成了復(fù)雜的裂縫形態(tài)，簡化為多裂縫或交錯分布的形態(tài)；離散裂縫模型包括三維的線網(wǎng)模型、二維的離散裂縫模型以及隨機(jī)分布的多裂縫模型（圖4）。

圖4 產(chǎn)能預(yù)測數(shù)值模型圖

Bruce等［8］建立了線網(wǎng)模型，采用正交的橢球面模擬網(wǎng)狀裂縫面，裂縫在空間上呈三維分布互相正交；考慮了穩(wěn)態(tài)和擬穩(wěn)態(tài)流動的情況，計算出了產(chǎn)能預(yù)測圖版。Cipolla等［9］建立了等距分布正交分布的網(wǎng)狀裂縫產(chǎn)能預(yù)測模型，采用SRV方法計算改造體積。模擬計算中考慮了支撐劑的分布形態(tài)（支撐裂縫和未支撐裂縫），滲流模式（解吸附、達(dá)西滲流和非達(dá)西滲流）以及裂縫參數(shù)對產(chǎn)能的影響。Williams等［10］通過裂縫識別技術(shù)確認(rèn)裂縫面和裂縫數(shù)量；結(jié)合成像測井確定的裂縫方位，建立裂縫分布模型后簡化為氣藏數(shù)值模擬模型，通過歷史擬合進(jìn)行模型校正和參數(shù)優(yōu)化。

4.2 雙重介質(zhì)模型

雙重介質(zhì)模型建立的基礎(chǔ)是Warren和Root提出的雙重孔隙介質(zhì)模型，該模型強(qiáng)調(diào)了裂縫性油藏的雙孔隙的本質(zhì)，簡化了裂縫性油藏中的連通性和與模型規(guī)模有關(guān)的非均質(zhì)性問題。

Zhang等［11］采用了微地震解釋獲得的潛在裂縫分布狀態(tài)結(jié)果，得到模擬主要橫切裂縫的雙重介質(zhì)油藏模型。計算中考慮了裂縫間距、非達(dá)西效應(yīng)等對產(chǎn)能的影響，表明產(chǎn)能對竄流因子參數(shù)最為敏感。

Changan等［12］建立了離散裂縫網(wǎng)絡(luò)模型，通過ESV方法確定壓裂有效體積，作為進(jìn)一步約束雙重介質(zhì)參數(shù)的條件。該方法利用巖石力學(xué)參數(shù)和施工參數(shù)估算平均裂縫寬度，建立裂縫強(qiáng)度、裂縫寬度和施工規(guī)模的關(guān)系式，通過和微地震監(jiān)測結(jié)果對比，形成超大規(guī)模的雙重介質(zhì)模擬。模擬計算中考慮了等溫吸附、不穩(wěn)定擴(kuò)散和竄流等問題。

Schepers等［13］建立了三重介質(zhì)模型，針對頁巖吸附氣和游離氣并存的特點，考慮了氣體的解吸附、擴(kuò)散和達(dá)西流動的滲流模式以及氣水兩相滲流規(guī)律，認(rèn)為水的存在反應(yīng)了部分頁巖氣藏生產(chǎn)初期產(chǎn)水的實際情況，只有基質(zhì)中氣體向裂縫中流入時才有含氣飽和度的存在，故產(chǎn)氣時間實際上出現(xiàn)了延緩，有利產(chǎn)能的準(zhǔn)確預(yù)測。

5 結(jié)論

由于頁巖氣藏特殊的賦存運移機(jī)理、復(fù)雜的滲流模式和裂縫形態(tài)等特殊性，需要建立針對頁巖氣壓裂特點的裂縫擴(kuò)展和產(chǎn)能預(yù)測模型。產(chǎn)能預(yù)測模型研究面臨的主要問題是復(fù)雜的裂縫形態(tài)和多變的滲流模式，目前的模型在裂縫形態(tài)上沒有考慮不規(guī)則的壓裂改造區(qū)域、網(wǎng)狀裂縫不規(guī)則分布的情況；裂縫擴(kuò)展模型沒能考慮天然裂縫和人工裂縫任意角度相交的問題；在滲流模式方面，未建立針對頁巖氣吸附或解吸附特點的吸附模型和納米級孔隙中微觀滲流模型，上述問題亟待進(jìn)行深入的探討和研究。

［1］RAY J，ROBERT C，ERY DIAZ－CAMPOS，et al.New pore－scale considerations for shale gas in place calculations［C］∥paper 131772presented at the SPE Unconventional Gas Conference，23－25February 2010，Pittsburgh，Pennsylvania，USA.New York：SPE，2010.

［2］MAYERHOFER M，LOLON E，WARPINSKI N，et al.What is stimulated reservoir volume（SRV）？［C］∥paper 119890presented at the SPE Shale Gas Production Conference，16－18November 2008，F(xiàn)ort Worth，Texas，USA.New York：SPE，2008.

［3］JON E，ARASH DAHI－TALEGHANI.Modeling simultaneous growth of multiple hydraulic fractures and their interaction with natural fractures［C］∥paper 119739presented at the SPE Hydraulic Fracturing Technology Conference，19－21January 2009，The Woodlands，Texas，USA.New York：SPE，2009.

［4］GU H，WENG X，LUND J，et al.Hydraulic fracture crossing natural fracture at non－orthogonal angles，a criterion，its validation and applications［C］∥paper 139984presented at the SPE Hydraulic Fracturing Technology Conference and Exhibition，24－26January 2011，The Woodlands，Texas，USA.New York：SPE，2011.

［5］XU Wenyue，THIERCELIN M，WALTON I.Characterization of hydraulically－induced shale fracture network using an analytical／semi－analytical model［C］∥paper 124697 presented at the SPE Annual Technical Conference and Exhibition，4－7October 2009，New Orleans，Louisiana，USA.New York：SPE，2009.

［6］WENG X，KRESSE O，COHEN C，et al.Modeling of hydraulic fracture network propagation in a naturally fractured formation［C］∥paper 140253presented at the SPE Hydraulic Fracturing Technology Conference and Exhibition，24－26January 2011，The Woodlands，Texas，USA.New York：SPE，2011.

［7］FISHER M，WRIGHT C，DAVIDSON B，et al.Integrating fracture mapping technologies to optimize stimulations in the Barnett Shale［C］∥paper 77441presented at the SPE Annual Technical Conference and Exhibition，29September－2October 2002，San Antonio，Texas，USA.New York：SPE，2002.

［8］BRUCE R.MEYER，LUCAS W.BAZAN.A discrete fracture network model for hydraulically induced fractures：theory，parametric and case studies［C］∥paper 140514presented at the SPE Hydraulic Fracturing Technology Conference and Exhibition，24－26January 2011，The Woodlands，Texas，USA.New York：SPE，2011.

［9］CIPOLLA C，LOLON E，MAYER HOFER M，et al.Reservoir modeling and production evaluation in shale－gas reservoirs［C］∥paper 13185presented at the International Petroleum Technology Conference，7－9December 2009，Doha.Qatar：IPTC，2009.

［10］WILLIAMS M J，KHADHRAOUI B，BRADFORD I.

Quantitative interpretation of major planes from microseismic event locations with application in production prediction［J］.SEG Expanded Abstracts 29.2085（2010）.［11］ZHANG X，DU C，DEIMBACHER F，et al.Sensitivity studies of horizontal wells with hydraulic fractures in shale gas reservoirs［C］∥paper 13388presented at the International Petroleum Technology Conference，7－9December 2009，Doha.Qatar：IPTC，2009.

［12］CHANGAN M，XU ZHANG，LANG ZHAN，et al.Modeling hydraulic fracturing induced fracture networks in shale gas reservoirs as a dual porosity system［C］∥paper 132180presented at the CPS／SPE International Oil ＆Gas Conference and Exhibition，8－10June 2010，Beijing，China.New York：SPE，2010.

［13］SCHEPERS K，GONZALEZ R，KOPERNA G，et al.Reservoir modeling in support of shale gas exploration［C］∥paper 123057presented at the SPE Latin American and Caribbean Petroleum Engineering Conference，31May－3June，2009，Cartagena，Colombia.New York：SPE，2009.

2011－09－27 編輯 韓曉渝）

DOI：10.3787／j.issn.1000－0976.2011.12.014

Zhang Shicheng，professor，born in 1963，is mainly engaged in teaching and research of theories，technologies and their application in oil production engineering，and oil ＆gas seepage flows.

Add：No.18，F(xiàn)uxue Rd.，Changping District，Beijing 102249，P.R.China

Tel：＋86－10－8973 3047 E－mail：zhangsc＠cup.edu.cn

Numerical simulation models with hydraulic fracturing in shale gas reservoirs

Zhang Shicheng1，Mou Songru1，Cui Yong2
（1.Key Laboratory of Petroleum Engineering，Ministry of Education∥China University of Petroleum，Beijing102249，China；2.Overseas Exploration ＆Development Company，PetroChina，Beijing 100034，China）

NATUR.GAS IND.VOLUME31，ISSUE12，pp.81－84，12／25／2011.（ISSN1000－0976；In Chinese）

Hydraulic fracturing treatment and horizontal well exploitation are major techniques in the shale gas development.Through a summary analysis of the characteristics of the hydraulic fracturing in shale gas reservoirs，the dominating factors resulting in the formation of net fractures are first studied，then the extended models of fractures and the types of fluid rate forecast models including their pros and cons are discussed.This study indicates that the major characteristics of hydraulic fracturing treatment in shale gas reservoirs include special mechanism of absorption，complex fractures，and multi－scale effect of seepage flows，with an aim of forming net fractures and enlarging the volume of reservoir stimulation.Furthermore，the formation of net fractures is mainly controlled by the factors such as the intersection angle of natural and hydraulic fractures，the maximum and the minimum horizontal principal stresses，and the brittleness of shales.The major problems facing the fluids rate forecast of hydraulic fracturing in shale gas reservoirs are the simulation of fracture morphology and the description of gas flow patterns，for which there are models for unconventional fractures and discrete fractures，and a dual porosity model.These models and methods represent，to some degree，the characteristics of the fracture morphology and seepage flows in the shale gas fracturing treatment but without considering the irregular fracture morphology.These research findings and cognitions give guidance to further shale gas development.

shale gas，development；hydraulic fracturing，complex fracture model，fluid rate forecast model，seepage flow，characteristics

張士誠，1963年生，教授，博士生導(dǎo)師，本刊第七屆編委會委員；長期從事采油工程理論與技術(shù)、油氣滲流理論與應(yīng)用的教學(xué)與研究工作。地址：（102249）北京市昌平區(qū)府學(xué)路18號。電話：（010）89733047。E－mail：zhangsc＠cup.edu.cn

張士誠等.頁巖氣壓裂數(shù)值模型分析.天然氣工業(yè)，2011，31（12）：81－84.

10.3787／j.issn.1000－0976.2011.12.014