摘 要

摘　要：對低滲透頁巖儲層進行體積壓裂改造以形成復雜裂縫網(wǎng)絡是獲得頁巖氣經(jīng)濟產(chǎn)能的關(guān)鍵，壓裂改造體積和縫網(wǎng)導流能力是評價體積壓裂施工效果的關(guān)鍵指標，同時對壓裂優(yōu)化設計、

摘　要：對低滲透頁巖儲層進行體積壓裂改造以形成復雜裂縫網(wǎng)絡是獲得頁巖氣經(jīng)濟產(chǎn)能的關(guān)鍵，壓裂改造體積和縫網(wǎng)導流能力是評價體積壓裂施工效果的關(guān)鍵指標，同時對壓裂優(yōu)化設計、壓后產(chǎn)能預測及經(jīng)濟評價也具有重要意義。為此，在分析頁巖氣體積壓裂特點的基礎上，對兩種主要頁巖氣體積壓裂縫網(wǎng)模型的假設、數(shù)學方程及參數(shù)優(yōu)化方法進行了比較分析，并結(jié)合美國Marcellus頁巖區(qū)塊現(xiàn)場參數(shù)對頁巖儲層壓裂方案進行了優(yōu)選。結(jié)果表明：離散化縫網(wǎng)模型及線網(wǎng)模型均能有效表征復雜縫網(wǎng)幾何特征，模擬縫網(wǎng)的擴展規(guī)律和縫網(wǎng)中壓裂液流動及支撐劑運移，獲得縫網(wǎng)幾何形態(tài)參數(shù)，可優(yōu)選壓裂施工方案；天然裂縫發(fā)育的頁巖層是體積壓裂改造的重點，水平地應力差越小則越易形成復雜縫網(wǎng)，施工排量越大，壓裂液泵入總量越大，則儲層改造體積范圍越大，縫網(wǎng)導流能力越高，頁巖氣產(chǎn)能就越高。

關(guān)鍵詞：頁巖氣  體積壓裂  離散化縫網(wǎng)模型  線網(wǎng)模型  排量  壓裂改造體積  縫網(wǎng)導流能力

Analysis and application of fracture network models of volume fracturing in shale gas reservoirs

Abstract：Volume fracturing in low permeability shale gas reservoirs which aims to induce a complicated fracture network is a key measure to acquire commercial productivity of shale gas．Stimulated reservoir volume and fracture network conductivity are important criteria of evaluating volume fracturing treatment effect and essential for fracturing optimization，production prediction after fracturing，and economic evaluation．Based on the analysis of volume fracturing characteristics of shale gas，a comparative analysis was made of the two dominant shale gas volume fracturing network models in terms of modeling assumptions，mathematical equations，and parameter optimization methods．In addition，we also selected an optimal shale reservoir fracturing plan in combin5tion with field parameters of the Marcellus Shale Block in the USA．It shows that both discrete fracture network model and wiremesh model can not only describe geometrical characteristics of a complicated fracture network，but simulate the extending rule of a fracture network and fracturing fluid flow patterns and proppant transport methodologies in a fracture network，thereby to acquire geometrical parameters of a fracture network and optimize a fracturing scheme．In terms of a shale gas reservoir with a natural。fracture system，it is a key factor of volume fracturing and easy to induce a complicated fracture network with a small horizontal crustal stress difference．Besides，the larger the pumping rate and fracturing fluid volume，the greater the stimulated reservoir volume，fracture network conductivity as well as shale gas productivity．

Key words：shale gas，volume fracturing，discrete fracture network model，wiremesh model，pumping rate，stimulated reservoir volume，fracture network conductivity

隨著北美地區(qū)頁巖氣勘探的巨大成功，頁巖氣資源的勘探開發(fā)在全球范圍內(nèi)持續(xù)升溫^[1-2]。美國頁巖氣商業(yè)開發(fā)的經(jīng)驗表明，水平井技術(shù)和水力壓裂技術(shù)是獲得頁巖氣工業(yè)氣流的關(guān)鍵技術(shù)^[3]。以體積壓裂^[4]技術(shù)為代表的增產(chǎn)技術(shù)將是未來開采頁巖一類特低滲透非常規(guī)油氣藏的核心技術(shù)，由于形成了復雜的裂縫網(wǎng)絡系統(tǒng)，增大了裂縫壁面與頁巖基質(zhì)的接觸面積，進而提高儲層整體滲透率，實現(xiàn)了對頁巖層整體上的三維壓裂改造。

在頁巖層進行體積壓裂時，由于頁巖特殊物理性質(zhì)及其內(nèi)部天然裂縫的影響，會產(chǎn)生一個水力裂縫與天然裂縫相互連通的復雜縫網(wǎng)系統(tǒng)^[5-7]。而傳統(tǒng)水力壓裂模型(二維、擬三維、全三維模型)都是基于雙翼對稱裂縫理論，假設裂縫為單一形態(tài)裂縫，不適用于天然裂縫及層理發(fā)育、各向異性突出的頁巖氣體積壓裂縫網(wǎng)系統(tǒng)的分析，因此需要建立專門的縫網(wǎng)壓裂模型來模擬頁巖層縫網(wǎng)幾何形態(tài)及其擴展規(guī)律^[8]。

筆者給出了離散化縫網(wǎng)模型和線網(wǎng)模型兩種頁巖氣體積壓裂縫網(wǎng)模型的原理、縫網(wǎng)幾何形態(tài)表征方法以及相應數(shù)學方程，分析了兩種模型的優(yōu)缺點，并利用兩種模型進行了敏感性因素分析，深化了對頁巖氣縫網(wǎng)壓裂優(yōu)化設計方法的理解。

1　體積壓裂縫網(wǎng)模型

1．1　離散化縫網(wǎng)模型

離散化縫網(wǎng)模型(DFN)最早由Meyer^[9-10]等人提出。該模型基于自相似原理及Warren和Root的雙重介質(zhì)模型，利用網(wǎng)格系統(tǒng)模擬解釋裂縫在3個主平面上的擬三維離散化擴展和支撐劑在縫網(wǎng)中的運移及鋪砂方式，通過連續(xù)性原理及網(wǎng)格計算方法獲得壓裂后縫網(wǎng)幾何形態(tài)。

DFN模型基本假設如下：①壓裂改造體積為2a×2b×h的橢球體，由直角坐標系XYZ表征，X軸平行于最大水平主應力(sH)方向，Y軸平行于最小水平主應力(sh)方向，Z軸平行于垂向應力(sv)方向；②包含一條主裂縫及多條次生裂縫，主裂縫垂直于sh方向，在X-Z平面內(nèi)擴展，次生裂縫分別垂直于X、Y、Z軸，縫間距分別為dx、dy、dz；③考慮縫間十擾及壓裂液濾失；④地層及流體不可壓縮?；谝陨霞僭O，作出DFN模型幾何模型的示意圖(圖1)。

 

DFN模型主要數(shù)學方程如下所示。

1)連續(xù)性方程

在考慮濾失的情況下，壓裂液泵入體積與濾失體積之差等于縫網(wǎng)中所含裂縫的總體積。即

 

式中q為壓裂液流量，m³／min；V1為濾失量，m³；Vsp為初濾失量，m³；Vf為總裂縫體積，m³。

2)流體流動方程

假設壓裂液在裂縫中的流動為層流，遵循冪率流體流動規(guī)律，其流動方程為：

 

式中P為縫內(nèi)流體壓力，MPa；n¢為流態(tài)指數(shù)，無因次；k¢為稠度系數(shù)，Pa·sⁿ；a、b分布為橢圓長軸半長及短軸半長，m；F( n¢)為積分函數(shù)，無因次。

3)縫寬方程

主裂縫縫寬方程為：

 

假設所有垂直于z軸的次生縫縫寬相同，與主裂縫縫寬之比為lz，則次生裂縫縫寬方程為：

wz=lzwx　　　　　　　　　　　　　　　　(4)

式中wx，為主裂縫縫寬，mm；Gw為功能函數(shù)；E為彈性模量，MPa；u為泊松比，無因次；sh。為最小水平主應力，MPa；Dsxx，為縫間干擾應力，MPa；z代表x、y、z；wz為垂直于z軸的次生縫縫寬，mm；lz為垂直于z軸的次生縫縫寬與主裂縫縫寬之比，無因次。

應用離散化縫網(wǎng)模型進行壓裂優(yōu)化設計時，需要首先設定次生裂縫縫寬、縫高、縫長等參數(shù)與主裂縫相應參數(shù)的關(guān)系，假設次生裂縫幾何分布參數(shù)；然后按設計支撐劑的沉降速度以及鋪砂方式，將地層物性，施工條件等參數(shù)代入以上數(shù)學模型，通過數(shù)值分析方法求得主裂縫的幾何形態(tài)和次生裂縫幾何形態(tài)；最后得到壓裂改造后的復雜縫網(wǎng)幾何形態(tài)。

DFN模型是目前模擬頁巖氣體積壓裂復雜縫網(wǎng)的成熟模型之一，特別是考慮了縫間干擾和壓裂液濾失間題后，更能夠準確描述縫網(wǎng)幾何形態(tài)及其內(nèi)部壓裂液流動規(guī)律，對縫網(wǎng)優(yōu)化設計具有重要意義。其不足之處在于需要人為設定次生裂縫與主裂縫的關(guān)系，主觀性強，約束條件差，且本質(zhì)上仍是擬三維模型。

1．2　線網(wǎng)模型

線網(wǎng)模型又稱HFN模型，首先由Xu等人^[11-14]提出，該模型基于流體滲流方程及連續(xù)性方程，同時考慮了流體與裂縫及裂縫之間的相互作用。

HFN模型基本假設如下：①壓裂改造體積為沿井軸對稱2a×2b×h的橢柱體，由直角坐標系XYZ表征，X軸平行于sH方向，Y軸平行于sh方向，Z軸平行于sv方向；②將縫網(wǎng)等效成兩簇分別垂直于X軸、Y軸的縫寬、縫高均恒定的裂縫，縫間距分別為dx、dy；③考慮流體與裂縫以及裂縫之間的相瓦作用；④不考慮壓裂液濾失。基于以上假設，做出HFN模型的幾何模型示意圖(圖2)。

 

HFN數(shù)學模型如下所示。

1)連續(xù)性方程

在不考慮壓裂液濾失的情況下，泵入壓裂液的體積與所形成裂縫的總體積相等。即

 

式中q為壓裂液流量，m³／min；ti為施工時間，min；h為裂縫縫高，m；Nx、Ny，分別為垂直于X軸、Y軸的裂縫的條數(shù)，無因次；Lxi、Lyi分別為垂直于于X軸的第i條裂縫和平行于Y軸的第j條裂縫的長度，m；wx、wy，分別為垂直于X軸裂縫和垂直于Y軸裂縫的縫寬，mm。

2)流體滲流方程

由于縫網(wǎng)內(nèi)裂縫寬度很小，因此可以假設流體在裂縫中的流動遵循流體滲流方程，則橢圓滲流方程為：

 

式中B為第二類橢圓積分，無因次；g為橢圓縱橫比，g=b／a；Kfx、Kfy，分別為垂直于X軸、Y軸方向的裂縫滲透率，mD；m為壓裂液黏度，mPa·s；j為裂縫孔隙度，無因次。

3)縫寬方程

假設相互平行的裂縫縫寬相同，則垂直于z方向的裂縫縫寬方程為：

 

式中z代表x、y；wz是垂直于z方向的裂縫縫寬，mm；E為彈性模量，MPa，u為泊松比，無因次；p為縫內(nèi)流體壓力，MPa；sc代表垂直于z方向的水平主應力sh或sH，MPa；Dszz為縫間干擾應力，MPa。

將方程(8)代入方程(5)(6)(7)，并聯(lián)立方程(5)(6)(7)可以獲得線網(wǎng)模型的方程組：

 

其中　　　　　　　　　　　Ds ＝ sH - sh

方程組(9)含有l2個變量，在已知其中9個變量的前提下，即可通過求解該方程組求得另外3個變量。因此在利用HFN模型進行縫網(wǎng)壓裂優(yōu)化設計時，可先通過微地震監(jiān)測獲得縫網(wǎng)幾何形態(tài)參數(shù)(h、a、b)，通過巖石力學常規(guī)三軸實驗獲得巖石物性參數(shù)(E、u)，根據(jù)壓裂施工方案獲得施工參數(shù)(p、q、t、m)，然后利用半解析法求解方程組(11)，獲得縫網(wǎng)分布參數(shù)(dx—dy)及差應力(Ds)。反之，若已知縫網(wǎng)分布及差應力，則通過計算HFN模型，可以獲得縫網(wǎng)幾何形態(tài)參數(shù)(h、a、b)。

HFN模型考慮了壓裂過程中改造體積的實時擴展以及施工參數(shù)的影響，能夠?qū)σ淹瓿蓧毫堰M行縫網(wǎng)分析，同時可以基于該分析對之后的壓裂改造方案進行二次優(yōu)化設計。其不足之處在于模擬縫網(wǎng)幾何形態(tài)較為簡單，需借助于地球物理技術(shù)的幫助獲取部分參數(shù)，同時由于不能模擬水平裂縫的起裂及擴展問題，及忽略了濾失問題，所以使用時具有較大的局限性。

2　縫網(wǎng)模型應用

Marcellus頁巖區(qū)塊是美國成功進行商業(yè)化開采的頁巖區(qū)塊之一，Henry等人給出了該區(qū)塊壓裂的基礎參數(shù)^[15]。因為巖石物性、施工參數(shù)及縫網(wǎng)分布參數(shù)均能對縫網(wǎng)幾何形態(tài)參數(shù)(h、a、b)、壓裂改造體積以及縫網(wǎng)導流能力產(chǎn)生影響，因此這里以Marcellus頁巖區(qū)塊壓裂水平井A為例，分別利用DFN模型和HFN模型重點分析縱橫比、施工排量以及壓裂液泵入總量的影響規(guī)律，進一步優(yōu)化壓裂施工參數(shù)。

A井垂深2380m，水平井段長l828m，目標頁巖層在2370～2397m，層厚27m，基質(zhì)滲透率為0.002mD，各層巖石物性參數(shù)如表1所示；本井共進表1巖石物性參數(shù)表行了七級清水壓裂，每一級有5簇射孔段，每一簇含8個射孔孔眼；選擇100目石英砂作為支撐劑，假設支撐劑在所有裂縫中均勻分布，沉降速率為1.27mm／min。兩個模型均假設縫網(wǎng)由兩簇分別垂直于最大、最小水平主應力的垂直縫組成，不含水平縫。

 

2．1　縱橫比影響規(guī)律

縱橫比b／a定義了縫網(wǎng)在Y方向和X方向延伸的最大長度之比，可通過微地震直接得到。假設壓裂施工參數(shù)恒定，排量為l6m³／min，總泵入液量為1600m³，縫間距dx=dy=l5m；改變縱橫比，其變化范圍為0.1至1，分別利用DFN模型和HFN模型計算縫網(wǎng)的幾何形態(tài)參數(shù)(h、a、b)，得到壓裂改造體積(VSRV)，然后再確定此縫網(wǎng)的導流能力，最終獲得縱橫比對縫網(wǎng)幾何形態(tài)參數(shù)、壓裂改造體積及縫網(wǎng)導流能力的影響規(guī)律。

DFN模型和HFN模型計算所得縫網(wǎng)最大縫高分別為27.9m和27.2m，不隨縱橫比的變化而變化。在相同壓裂施工條件下，隨著縱橫比的增大，a逐漸降低，b逐漸增大，壓裂改造體積逐漸增大；對于DFN模型，a由500m降低至l76m，b由50m增大至l76m，壓裂改造體積由l47×10⁴m³增大至179×l0⁴m³；對于HFN模型，a由435m降低至158m，b由43.5m增大至158m，壓裂改造體積由214×10⁴m³增大至264×10⁴m³。在相同縱橫比下，利用DFN模型計算的縫網(wǎng)幾何參數(shù)比HFN計算值大(圖3-a)；而利用DFN模型計算的壓裂改造體積比HFN模型計算值小(圖3-b)。

 

縫網(wǎng)導流能力為縫網(wǎng)平均縫寬與縫網(wǎng)平均滲透率之積，反映了流體在縫網(wǎng)中的流動能力。由圖4可知：隨著縱橫比的增大，縫網(wǎng)導流能力降低；對于DFN模型，縫網(wǎng)導流能力從18lmD·m降低至117mD·m；對于HFN模型，縫網(wǎng)導流能力從l67mD·m降低至126mD·m。當縱橫比小于0.6時，DFN模型計算得到的縫網(wǎng)導流能力偏大；當縱橫比大于0.6時，HFN模型計算值偏大。

 

縱橫比主要受水平地應力差的控制，水平地應力差越小，縱橫比越大，壓裂改造體積越大而縫網(wǎng)導流能力越小。因此在優(yōu)選壓裂層位時，應該綜合考慮目標壓裂改造體積與目標縫網(wǎng)導流能力，優(yōu)選具有合適水平地應力差的地層。

2．2　排量影響規(guī)律

壓裂過程中，排量為可控施工參數(shù)。假設壓裂液總泵入量為1600m³，dx＝dy＝l5m，g＝0.2；改變壓裂液排量，其變化范圍為2m³／min到20m³／min，分別利用DFN模型和HFN模型計算縫網(wǎng)的幾何形態(tài)參數(shù)，得到壓裂改造體積，然后再確定此縫網(wǎng)的導流能力，最終得到排量對縫網(wǎng)幾何形態(tài)參數(shù)、壓裂改造體積及縫網(wǎng)導流能力的影響規(guī)律。

DFN模型和HFN模型計算所得縫網(wǎng)最大縫高分別為27.9m和27.2m，不隨排量的變化而變化。在其他條件相同的情況下，隨著排量的增大，a和b逐漸增大，壓裂改造體積逐漸增大；對于DFN模型，“由243m增加至380m，b由49m增大至76m，壓裂改造體積由68×10⁴m³增大至169×10⁴m³；對于HFN模型，以由218m增加至325m，b由44m增大至65m，壓裂改造體積由82×10⁴m³增大至181×10⁴m³。在相同排量下，利用DFN模型預測的縫網(wǎng)幾何參數(shù)比HFN模型預測值大(圖5-a)；而利用DFN模型計算的壓裂改造體積比HFN模型計算值小(圖5-b)。

 

由圖6可知，隨著排量的增大，縫網(wǎng)導流能力增大；對于DFN模型，縫網(wǎng)導流能力從ll4mD·m增大至l63mD·m；對于HFN模型，縫網(wǎng)導流能力從106mD·m增大至169mD·m。當流量低于l4m³／min時，DFN模型計算的縫網(wǎng)導流能力偏大；當流量高于l4m³／min時，HFN模型計算值偏大。

 

因此，在施工條件允許范圍內(nèi)，增大排量，有利于形成具有高導流能力和大壓裂改造體積的復雜縫網(wǎng)。

2．3　壓裂液泵入總量影響規(guī)律

假設壓裂施工排量為16m³／min，dx＝dy＝l5m，g＝0.2；改變壓裂液泵入總量，變化范圍為200～2000m³，分別利用DFN模型和HFN模型計算縫網(wǎng)的幾何形態(tài)參數(shù)，得到壓裂改造體積，然后再確定此縫網(wǎng)的導流能力，最終獲得壓裂液泵入總量對縫網(wǎng)幾何形態(tài)參數(shù)、壓裂改造體積及縫網(wǎng)導流能力的影響規(guī)律。

DFN模型和HFN模型計算所得縫網(wǎng)最大縫高分別為27.9m和27.2m，不隨壓裂液泵入總量的變化而變化。在其他條件相同的情況下，隨著壓裂液泵入總量的增大，a和b逐漸增大，壓裂改造體積逐漸增大；對于DFN模型，a由l71m增加至395m，b由34m增大至79m，壓裂改造體積由34×10⁴m³增大至183×10⁴m³；對于HFN模型，以由153m增加至340m，b由31m增大至68m，壓裂改造體積由40×10⁴m³增大至l97×l0⁴m³；在泵入相同體積壓裂液的情況下，利用DFN模型計算的縫網(wǎng)幾何參數(shù)比HFN模型計算值大(圖7-a)；而利用DFN模型計算的壓裂改造體積比HFN模型計算值小(圖7-b)。

 

由圖8可知，隨著壓裂液泵入總量的增大，縫網(wǎng)導流能力增大；對于DFN模型，縫網(wǎng)導流能力從l45mD·m增大至l58mD·m，且當壓裂液泵入總量超過l200m³后，縫網(wǎng)導流能力基本保持不變；對于HFN模型，縫網(wǎng)導流能力從l34mD·m增大至l63mD·m。當壓裂液泵入總量低于l400m³時，DFN模型計算得到的縫網(wǎng)導流能力偏大；當壓裂液泵入總量高于1400m³時，HFN模型計算值偏大。

 

因此，在施工條件允許范圍內(nèi)，適當增大壓裂液泵入總量有利于形成具有高導流能力和大壓裂改造體積的復雜縫網(wǎng)。

3　結(jié)論

1)離散化縫網(wǎng)模型及線網(wǎng)模型均能有效表征復雜縫網(wǎng)幾何特征，模擬縫網(wǎng)中壓裂液的流動及支撐劑的運移，獲得縫網(wǎng)擴展規(guī)律及幾何形態(tài)參數(shù)，優(yōu)選壓裂施工方案。

2)離散化縫網(wǎng)模型基于縫內(nèi)流體流動理論，考慮濾失及支撐劑運移，人為設定次生裂縫分布及與主裂縫幾何參數(shù)關(guān)系，可以對水平井及直井進行參數(shù)優(yōu)化，但人為假設次生縫主觀性較強；線網(wǎng)模型基于流體滲流理淪，能夠模擬裂縫的實時擴展，但需要微地震進行約束，另外由于忽略了濾失問題，且僅能對垂直裂縫進行有效模擬，所以有一定局限性。

3)離散化縫網(wǎng)模型假設壓裂改造體積為橢球體，線網(wǎng)模型假設壓裂改造體積為橢柱體。對于縫網(wǎng)幾何參數(shù)(h、a、b)，離散化縫網(wǎng)模型計算值比線網(wǎng)模型計算值大；而對于壓裂改造體積，離散化縫網(wǎng)模型汁算值比線網(wǎng)模型計算值小。

4)地應力差越小，縱橫比越大，縫網(wǎng)壓裂改造體積越大，導流能力越低。因此在壓裂層位優(yōu)選時，需綜合考慮目標壓裂改造體積和目標縫網(wǎng)導流能力，選擇具有合適水平地應力差的頁巖層。

5)在施工條件允許范圍內(nèi)，增大壓裂施工排量和壓裂液泵入總量，有和于形成具有大壓裂改造體積及高導流能力的復雜裂縫網(wǎng)絡系統(tǒng)。

 

參考文獻

[1]蒲泊伶，蔣有錄，王毅，等．四川盆地志留統(tǒng)龍馬溪組頁巖氣成藏條件及有利地區(qū)分析[J]．石油學報，2010，3l(2)：225-231．

PU Boling，JIANG Youlu，WANG Yi，et al．Analysis of the reservoir-forming conditions of shale gas potential in Sichuan Basin[J]．Aeta Petrolei Sinica，2010，3l(2)：225-231．

[2]鄒才能，董大忠，楊樺，等．中國頁巖氣形成條件及勘探實踐[J]．天然氣工業(yè)，2011，31(12)：26-39．

ZOU Caineng，DONG Dazhong，YANG Hua，et al．Conditions of shale gas accumulation and exploration practices in China[J]．Natural Gas industry，2011，31(12)：26-39．

[3]張抗．頁巖氣發(fā)展的中國之鑒一致密油氣和煤層氣[J]．天然氣工業(yè)，2013，33(4)：18-25．

ZHANG Kang．Lessons for shale oil&gas development from that of tight oil&gas and coalbed methane gas in China[J]．Natural Gas Industry，2013，33(4)：18-25．

[4]王海慶，王勤．體積壓裂在超低滲油藏中的開發(fā)應用[J]．中國石油和化工標準與質(zhì)量，2012，32(2)：143．

WANG Haiqing，WANG Qin．Application of volume fracturing in ultra low permeability reservoirs[J]．China Petroleum and Chemical Standard and Quality，2012，32(2)：143．

[5]杜林麟，春蘭，王玉艷，等．頁巖儲層水力壓裂優(yōu)化設計[J]．石油鉆采工藝，2010，32(增刊l)：l30-132．

DU Linlin，CHUN Lan，WANG Yuyan，et al．Hydraulic fracturing optimization for shale reservoirs LJ]．Oil Drilling &Production Technology，2010，32(S1)：130-132．

[6]雷群，胥云，蔣廷學，等．用于提高低  特低滲油氣藏改造效果的縫網(wǎng)壓裂技術(shù)[J]．石油學報，2009，30(2)：237-241．

LEI Qun，XU Yun，JIANG Tingxue，et al．“Fracture net work”fracturing technique for improving post-fracturing performance of low and ultra low permeability reservoirs[J]．Aeta Petrolei Sinica，2009，30(2)：237-241．

[7]張小濤，吳建發(fā)，馮曦，等．頁巖氣藏水平井分段壓裂滲流特征數(shù)值模擬[J]．天然氣工業(yè)，2013，33(3)：47-52．

ZHANG Xiaotao，WU Jianfa，FENG Xi，et al．Numerical simulation of seepage flow characteristics of multi-stage fracturing(MSF)in horizontal shale gas wells[J]．Natural Gas Industry，2013，33(3)：47-52．

[8]張士誠，牟松茹，崔勇．頁巖氣壓裂數(shù)值模型分析[J]．天然氣工業(yè)，20ll，31(12)：81-84．

ZHANG Shicheng，MU Songru，CUI Yong．Numerical simulation models with hydraulic[racturing in shale gas reservoirs[J]．Natural Gas Industry，2011，3l(12)：81-84．

[9]MEYER B R，BAZAN L W．A discrete fracture network model for hydraulically induced fraclures：Theory，parametric and case studies[C]//paper l40514-MS presented at the SPE Hydraulic Fracturing Technology Conference and Exhibition，24-26 January 2011，the Woodlands，Texas，USA．New York：SPE，2011．

[10] MEYER B R，BAZAN L W，JACOT R H，et al．Optimi zation of muhiple transverse hydraulic fractures in horizontal wellbores[C]//paper l31732 MS presented at the SPE Unconventional Gas Conference，23-25 February 2010，Pittsburgh，Pennsylvania，USA．New York：SPE，2009．

[11]XU W X，THIERCELIN M，GANGUl．Y U．Wiremesh：A novel shale fracturing simulator[C]//paper l40514-MS presented at the CPS／SPE International Oil&Gas Conference and Exhibition，8-10 June 2010，Beijing，China．New York：SPE，2010．

[12]XU W X，CALVEZJ L，THIERCELIN M．Characterization of hydraulically induced fracture network using treatment and microseismic data in a tight gas sand formation： A geomechanical approach[C]//paper l25237-MS presented at the SPE Tight Gas Completions Conference，15-17 June 2009，San Antoni0，Texas，USA．New York：SPE，2009．

[13]XU W X，THIERCELIN M，WALTON I．Characterization of hydraulically induced shale fracture network using an analytical／semi-analytical model[C]//paper l24697-MS presented at the SPE Annual Technical Conference and Exhibition，4-7 0ctober 2009，New Orleans．I．ouisiana，USA．New York：SPE，2009．

[14]XU W X，THIERCELIN M，CALVEZ J L，et al．Fracture network development and proppant placement during slickwater fracturing treatment of Barnett Shale laterals[C]//paper l35484-MS presented at the SPE Annual Technical Conference and Exhibition，19-22 September 2010，Florence，Italy．New York：SPE，2010．

[15]JACOT R H，BAZAN L W，MEYER B R．Technology in tegration：A methodology to enhance production and maximize economic in horizontal Marcellus shale wells[C]//paper l35262-MS presented at the SPE Annual Technical Conference and Exhibition,19-22 September 2010．Florence，Italy．New York：SPE，2010．

 

本文作者：程遠方  李友志  時賢  吳百烈  王欣  鄧文彪

作者單位：中國石油大學(華東)石油工程學院

  中國石油勘探開發(fā)研究院廊坊分院