李玉梅， 呂　煒， 宋　杰， 李　軍， 楊宏偉， 于麗維

(1.中國石油大學(xué)(北京)石油工程學(xué)院，北京 102249；2.中國石油玉門油田分公司勘探開發(fā)研究院，甘肅酒泉 735019；3.中國石油新疆油田分公司工程技術(shù)研究院，新疆克拉瑪依 834000)

?

層理性頁巖氣儲(chǔ)層復(fù)雜網(wǎng)絡(luò)裂縫數(shù)值模擬研究

李玉梅1， 呂煒2， 宋杰2， 李軍1， 楊宏偉1， 于麗維3

(1.中國石油大學(xué)(北京)石油工程學(xué)院，北京 102249；2.中國石油玉門油田分公司勘探開發(fā)研究院，甘肅酒泉 735019；3.中國石油新疆油田分公司工程技術(shù)研究院，新疆克拉瑪依 834000)

在考慮天然裂縫的條件下，為了更好地對(duì)層理性頁巖氣儲(chǔ)層的復(fù)雜網(wǎng)絡(luò)裂縫進(jìn)行數(shù)值模擬，采用模擬非連續(xù)介質(zhì)的通用離散元程序(UDEC)，基于滲流-應(yīng)力耦合數(shù)值算法，建立了人工水力裂縫與天然裂縫相互作用的網(wǎng)絡(luò)裂縫數(shù)值計(jì)算模型，并利用該模型分析了水力裂縫長度、天然裂縫傾角、內(nèi)摩擦角及施工凈壓力對(duì)縫網(wǎng)擴(kuò)展的影響。結(jié)果表明：水力裂縫從近井筒處裂縫尖端起裂擴(kuò)展，并沿著天然裂縫的走向發(fā)生剪切破壞，且隨水力裂縫長度增長，天然裂縫網(wǎng)絡(luò)連通面積增大；天然裂縫傾角較大，形成復(fù)雜縫或網(wǎng)絡(luò)縫的概率也相對(duì)較大；天然裂縫內(nèi)摩擦角越小，天然裂縫連通面積越大，越易形成復(fù)雜網(wǎng)絡(luò)裂縫；水平地應(yīng)力差在一定范圍內(nèi)，凈壓力系數(shù)越大，裂縫的擴(kuò)展形態(tài)越復(fù)雜，相鄰裂縫的尖端越易連通形成網(wǎng)絡(luò)裂縫。數(shù)值模擬研究結(jié)果可為進(jìn)一步認(rèn)識(shí)遠(yuǎn)井地帶頁巖氣壓裂裂縫擴(kuò)展機(jī)理提供指導(dǎo)。

頁巖氣；壓裂；數(shù)值模擬；離散元；幾何模型；天然裂縫；內(nèi)摩擦角

頁巖富含天然裂隙和薄頁狀層理，不同的天然層理產(chǎn)狀組合可形成錯(cuò)綜復(fù)雜的離散性裂縫體系，增大裂縫表面與基質(zhì)體的接觸面積，可為頁巖氣藏“體積壓裂”提供地質(zhì)條件。體積壓裂縫網(wǎng)的形成是水力壓裂作業(yè)、天然裂縫及非均勻地應(yīng)力共同作用的結(jié)果。傳統(tǒng)壓裂理論認(rèn)為水力裂縫為對(duì)稱于井筒的雙翼裂縫且沿著最大水平主應(yīng)力方向擴(kuò)展，但事實(shí)上，受天然裂縫的干擾，裂縫性儲(chǔ)層壓裂后的水力裂縫網(wǎng)絡(luò)結(jié)構(gòu)呈明顯的離散性和多樣性[1-4]。網(wǎng)絡(luò)裂縫的形成與天然裂縫的發(fā)育程度相關(guān)，頁巖氣能否有效產(chǎn)出取決于水力裂縫和壓裂過程中誘導(dǎo)天然裂縫開啟而形成的相互交錯(cuò)的網(wǎng)絡(luò)裂縫的面積[5-7]。

近年來，國內(nèi)外學(xué)者對(duì)結(jié)構(gòu)面等非連續(xù)體與水力裂縫相互耦合的關(guān)系進(jìn)行了較多研究，并在天然裂縫建模和研究方法等方面均取得了重要進(jìn)展。T.L.Blanton[8]基于天然裂縫的走向和傾角對(duì)水力裂縫開啟的影響，研究了天然裂縫的受力狀態(tài)和破壞方式。W.S.Dershowitz等人[9]根據(jù)現(xiàn)場(chǎng)獲取的微地震裂縫數(shù)據(jù)，利用FracMan裂縫建模軟件建立了三維離散裂縫網(wǎng)絡(luò)模型，研究了天然裂縫存在條件下水力裂縫的擴(kuò)展形態(tài)。A.Riahi等人[10]利用離散元數(shù)值模擬方法建立了流固耦合離散裂縫網(wǎng)絡(luò)(DFN)數(shù)值模型，研究了壓裂液注入裂縫性儲(chǔ)層后對(duì)裂縫連通性、裂縫分布、裂縫開度以及滑移等的影響規(guī)律。M.Cottrell等人[11]同時(shí)考慮天然裂縫幾何參數(shù)、壓裂液注入水力參數(shù)和水力裂縫與天然裂縫的相互作用，建立了離散裂縫網(wǎng)絡(luò)(DFN)數(shù)值模型，并與有限元數(shù)值法進(jìn)行了對(duì)比。M.McClure等人[12]在假設(shè)裂縫高度一定的條件下，考慮水力裂縫與天然裂縫相互連通的作用，利用CFRAC計(jì)算程序建立了離散裂縫數(shù)值模型，并與PKN模型進(jìn)行了比較，分析了裂縫長度、裂縫開度、縫內(nèi)凈壓力等參數(shù)隨時(shí)間的變化規(guī)律。但是，以上研究都沒有考慮天然裂縫隨機(jī)分布特征對(duì)天然裂縫開啟的影響。N.B.Nagel等人[13-14]利用離散元數(shù)值計(jì)算方法建立了水力裂縫與天然裂縫相互作用的數(shù)值模型，并對(duì)天然裂縫的的剪切破壞行為進(jìn)行了分析。筆者在N.B.Nagel等人的研究基礎(chǔ)上，考慮了天然裂縫的隨機(jī)分布特征，利用通用離散元程序(UDEC)中獨(dú)有的天然裂縫建模技術(shù)建立了水力裂縫與天然裂縫相互干擾的離散裂縫網(wǎng)絡(luò)(DFN)數(shù)值模型，對(duì)層理性頁巖氣儲(chǔ)層的復(fù)雜網(wǎng)絡(luò)裂縫進(jìn)行了數(shù)值模擬。

1　天然裂縫剪切破壞機(jī)制

(1)

(2)

圖1　單裂縫變形受力示意Fig.1　Schematic diagram of the stress and deformation around a single fracture

剪切應(yīng)力為：

(3)

(4)

式中：Δσs為剪切應(yīng)力，Pa;ks為剪切剛度，Pa/m；Δus為剪切位移，m；F為表征剪切位移增量的相關(guān)系數(shù)；as和es為模型參數(shù)。

隨著裂縫開度增大以及裂縫的滑移，裂縫剪切應(yīng)力超過了剪切強(qiáng)度(σc)時(shí)[17]，裂縫就會(huì)發(fā)生剪切膨脹。裂縫開度的變化與剪切位移變化之間的關(guān)系為：

(5)

式中：Δun-d為裂縫開度的變化，m；φd為膨脹角，rad。

2　離散元數(shù)值模型的建立

2.1幾何模型的建立

頁巖具有薄頁狀或薄片層狀的節(jié)理，不同的節(jié)理產(chǎn)狀組合形成錯(cuò)綜復(fù)雜的裂縫網(wǎng)絡(luò)，根據(jù)現(xiàn)場(chǎng)的頁巖露頭可以觀察到天然裂縫傾角排布具有一定的規(guī)律性。筆者基于地震源破碎平面的方向和傾角確定天然裂縫的方向，根據(jù)微地震事件的離散統(tǒng)計(jì)確定天然裂縫走向，以威201-H1井儲(chǔ)層微地震事件的離散性隨機(jī)裂縫統(tǒng)計(jì)結(jié)果作為建模依據(jù)。

依據(jù)天然裂縫在地質(zhì)學(xué)上的合理分布，不考慮頁巖氣儲(chǔ)層中的大量不規(guī)則裂縫，只考慮主要天然裂縫。利用UDEC軟件獨(dú)有的天然離散裂縫建模技術(shù)，將天然裂縫隨機(jī)分布函數(shù)內(nèi)嵌至裂縫生成器中，建立了人工裂縫與天然裂縫相互干擾的離散裂縫網(wǎng)絡(luò)(DFN)數(shù)值模型，模型尺寸為600 m×600 m。圖2為網(wǎng)絡(luò)裂縫數(shù)值計(jì)算幾何模型，水平井筒沿x軸方向(最小水平主應(yīng)力方向σh)，水力裂縫垂直于井筒方向并沿著y軸方向擴(kuò)展(最大水平主應(yīng)力方向σH)。模型中預(yù)設(shè)了不同傾角組合的天然裂縫，相鄰天然裂縫均為閉環(huán)式連通。所建DFN數(shù)值模型為規(guī)則模型，假設(shè)水力裂縫為垂直裂縫，所以借助模型兩側(cè)邊線進(jìn)行創(chuàng)建，兩條邊線間的距離即為裂縫間距。邊界條件設(shè)置過程中，對(duì)模型左右兩側(cè)邊線和上下兩側(cè)邊線分別進(jìn)行x方向和y方向的位移約束，以約束數(shù)值模型位移。假設(shè)壓裂液從井筒注入后沿水力裂縫進(jìn)入，所以在水力裂縫內(nèi)表面單元施加凈壓力，模擬隨著縫內(nèi)液體壓力的逐漸升高，水力裂縫擴(kuò)展過程中產(chǎn)生的誘導(dǎo)應(yīng)力使其周圍天然裂縫的剪切應(yīng)力超過剪切強(qiáng)度而發(fā)生破壞，從而通過水力裂縫溝通儲(chǔ)層中的天然裂縫而形成縫網(wǎng)。

圖2　網(wǎng)絡(luò)裂縫數(shù)值計(jì)算幾何模型Fig.2　Geometry model for the numerical simulation of network fractures

2.2力學(xué)參數(shù)設(shè)置

網(wǎng)絡(luò)裂縫數(shù)值模型建立后，參考N.B.Nagel等人建立的模型中設(shè)定的頁巖氣儲(chǔ)層的地層參數(shù)、天然裂縫幾何參數(shù)、頁巖巖塊及天然裂縫基本力學(xué)參數(shù)[18]，設(shè)定主要力學(xué)參數(shù)為：地層流體密度為1.0 g/cm3，垂直地應(yīng)力為94.5 MPa，最大水平主應(yīng)力為90.5 MPa，最小水平主應(yīng)力為86.2 MPa，地層初始孔隙壓力為77.91 MPa。

天然裂縫產(chǎn)狀數(shù)據(jù)是數(shù)值計(jì)算分析的前提和必要條件，天然裂縫組合不同，裂縫網(wǎng)絡(luò)模型不同，為此，建立了組合A(0°/70°))和組合B(45°/135°)兩組離散裂縫網(wǎng)絡(luò)模型。每一組天然裂縫網(wǎng)絡(luò)模型都設(shè)定了合理的天然裂縫密度、長度、間距等參數(shù)(見表1)。頁巖巖樣及天然裂縫基本力學(xué)參數(shù)為：頁巖巖塊的密度為2 440 kg/m3，體積模量為17.1 GPa，剪切模量為16.4 GPa，內(nèi)聚力為5 MPa，內(nèi)摩擦角為10°～40°，抗拉強(qiáng)度為2 MPa；裂縫的法向剛度為20 GPa/m，切向剛度為20 GPa/m，內(nèi)聚力為5 MPa，內(nèi)摩擦角為19°，滲透率為0.001 mD。

表1　天然裂縫幾何參數(shù)

3　數(shù)值模擬結(jié)果分析

3.1水力裂縫縫長對(duì)縫網(wǎng)擴(kuò)展的影響

預(yù)置水力裂縫半長為50，100，200和300 m，分析天然裂縫傾角組合A和組合B條件下水力裂縫縫長對(duì)縫網(wǎng)擴(kuò)展的影響規(guī)律，結(jié)果見圖3—圖5。由圖3可以看出，應(yīng)力集中區(qū)分布在水力裂縫尖端，有利于裂縫向遠(yuǎn)井地帶擴(kuò)展。由圖4和圖5可以看出：水力裂縫從裂縫尖端開始擴(kuò)展并沿著天然裂縫的走向發(fā)生剪切破壞，且有向更深遠(yuǎn)地層延伸的趨勢(shì)；隨著縫長增長，剪切破壞區(qū)域增大。

利用二維離散裂縫數(shù)值模型，計(jì)算了組合A和組合B條件下，不同水力裂縫縫長度下的剪切破環(huán)面積，結(jié)果見圖6。由圖6可以看出：水力裂縫縫長相同時(shí)，組合B條件下的剪切破壞面積要大于組合A條件下的剪切破壞面積；隨著水力裂縫增長，組合B與組合A的剪切破壞面積差逐漸增大。分析其原因?yàn)椋禾烊涣芽p傾角較大條件下，最小水平主應(yīng)力方向(垂直水力裂縫方向)的誘導(dǎo)應(yīng)力相對(duì)變化較大，干擾強(qiáng)度增大，形成復(fù)雜縫或網(wǎng)絡(luò)縫的概率相對(duì)較大。

3.2天然裂縫內(nèi)摩擦角對(duì)網(wǎng)絡(luò)裂縫擴(kuò)展的影響

理論研究認(rèn)為，裂縫性頁巖氣儲(chǔ)層水力裂縫擴(kuò)展過程中對(duì)天然裂縫形成干擾，遵循天然裂縫張開和剪切破壞機(jī)制。當(dāng)天然裂縫在水力作用下剪切應(yīng)力超過裂縫剪切強(qiáng)度時(shí)，認(rèn)為閉合的天然裂縫開啟或者沿天然裂隙面發(fā)生剪切滑移。假設(shè)天然裂縫開啟的區(qū)域面積為網(wǎng)絡(luò)裂縫的連通面積，選取開啟天然裂縫擴(kuò)展區(qū)域邊界節(jié)點(diǎn)位置坐標(biāo)的方法計(jì)算網(wǎng)絡(luò)裂縫面積。數(shù)值模型預(yù)置水力裂縫半長為200 m，計(jì)算天然裂縫傾角組合A和組合B條件下，不同裂縫內(nèi)摩擦角下的天然裂縫網(wǎng)絡(luò)連通面積，結(jié)果見圖7。

圖3　水力裂縫尖端應(yīng)力集中區(qū)Fig.3　Stress concentration area of hydraulic fracture tips

圖4　天然裂縫傾角組合A(0°/70°)條件下水力裂縫縫長對(duì)縫網(wǎng)擴(kuò)展的影響Fig.4　Effect of hydraulic fracture length on fracture network extension under natural fracture dip angles combination A (0°/70°)

由圖7可以看出：天然裂縫內(nèi)摩擦角相同時(shí)，組合B條件下天然裂縫網(wǎng)絡(luò)連通面積要大于組合A條件下的天然裂縫網(wǎng)絡(luò)連通面積；隨天然裂縫內(nèi)摩擦角增大，天然裂縫網(wǎng)絡(luò)連通面積減小。分析其原因?yàn)椋禾烊涣芽p傾角較大條件下，最小水平主應(yīng)力方向(垂直水力裂縫方向)誘導(dǎo)應(yīng)力的相對(duì)變化量增大，干擾強(qiáng)度增大，形成復(fù)雜縫或網(wǎng)絡(luò)縫的概率相對(duì)較大；天然裂縫內(nèi)摩擦角越小，天然裂縫網(wǎng)絡(luò)連通面積越大，這是由于較小的內(nèi)摩擦角導(dǎo)致巖石抗剪強(qiáng)度減弱，裂縫容易發(fā)生剪切破壞，儲(chǔ)層改造易于形成復(fù)雜網(wǎng)絡(luò)裂縫。

3.3施工凈壓力對(duì)縫網(wǎng)擴(kuò)展的影響

水平井多級(jí)分簇射孔壓裂可以形成水力裂縫與天然裂縫相互交織的復(fù)雜網(wǎng)絡(luò)裂縫。針對(duì)裂縫性儲(chǔ)層，J.E.Olson等人提出采用凈壓力系數(shù)Rn來表征施工凈壓力對(duì)裂縫延伸的影響[19]，認(rèn)為水力裂縫與天然裂縫之間的應(yīng)力干擾是相對(duì)凈壓力的函數(shù)，在裂縫延伸過程中凈壓力系數(shù)與水平地應(yīng)力差成反比關(guān)系：

圖5　天然裂縫傾角組合B(45°/135°)條件下水力裂縫縫長對(duì)縫網(wǎng)擴(kuò)展的影響Fig.5　Effect of hydraulic fracture length on fracture network extension under natural fracture dip angle combination B(45°/135°)

圖6　水力裂縫縫長與剪切破壞面積的關(guān)系Fig.6　Relationship between hydraulic fracture length and shear failure area

(6)

式中：pf為裂縫內(nèi)的流體壓力，MPa；σH和σh分別為最大水平主應(yīng)力和最小水平主應(yīng)力，MPa。

以天然裂縫傾角組合A為例，建立水力裂縫與天然裂縫相互作用的數(shù)值模型，假設(shè)相鄰兩條水力裂縫半長均為200 m，裂縫間距為100 m。當(dāng)初始水平地應(yīng)力差為2～10 MPa時(shí)，分析凈壓力系數(shù)Rn對(duì)相鄰裂縫間網(wǎng)絡(luò)裂縫擴(kuò)展的影響，結(jié)果見圖8。

由圖8可以看出，裂縫性頁巖氣儲(chǔ)層在壓裂過程中形成的水力裂縫不是一條單一的平面縫，而是

圖7　裂縫內(nèi)摩擦角與天然裂縫網(wǎng)絡(luò)連通面積的關(guān)系Fig.7　Relationship between internal friction angles and connected area of natural fracture network

一種復(fù)雜的網(wǎng)狀裂縫。這是由于水力裂縫擴(kuò)展過程中誘導(dǎo)應(yīng)力場(chǎng)的產(chǎn)生促使裂縫剪切應(yīng)力超過裂縫剪切強(qiáng)度，使天然裂縫開啟并與主裂縫連通。主裂縫表現(xiàn)出的這種動(dòng)態(tài)分叉效應(yīng)具有隨機(jī)性和條件性，也是整體裂縫形態(tài)變得復(fù)雜的主要因素。由圖8還可以看出，水平地應(yīng)力差在一定范圍內(nèi)，凈壓力系數(shù)越大，裂縫的擴(kuò)展形態(tài)越復(fù)雜，相鄰裂縫的尖端越易連通形成網(wǎng)絡(luò)裂縫。造成這種現(xiàn)象的主要原因是，較小的初始水平地應(yīng)力差導(dǎo)致剪切裂縫轉(zhuǎn)向最大水平主應(yīng)力方向，在遠(yuǎn)井地帶更容易與天然裂縫溝通形成復(fù)雜的網(wǎng)絡(luò)裂縫，增加儲(chǔ)層改造面積。因此，在水力壓裂過程中，低地應(yīng)力差可以增大形成復(fù)雜網(wǎng)絡(luò)裂縫的概率，而高應(yīng)力差則會(huì)降低形成復(fù)雜網(wǎng)絡(luò)裂縫的概率。但是，考慮到工程地質(zhì)和作業(yè)情況，水平應(yīng)力差較低或壓裂液濾失量增大均可能導(dǎo)致近井筒處過早形成網(wǎng)絡(luò)裂縫而形成砂堵，導(dǎo)致裂縫無法向遠(yuǎn)井地帶擴(kuò)展。因此，制定壓裂作業(yè)方案時(shí)要綜合考慮多因素對(duì)水力壓裂網(wǎng)絡(luò)裂縫擴(kuò)展的影響。

圖8　凈壓力系數(shù)對(duì)相鄰兩條水力裂縫間網(wǎng)絡(luò)裂縫擴(kuò)展的影響Fig.8　The effect of net-pressure coefficients on network fracture extension between two adjacent fractures

4　結(jié)　論

1) 水力裂縫從近井筒處裂縫尖端起裂擴(kuò)展并沿著天然裂縫的走向發(fā)生剪切破壞。在天然裂縫傾角較大條件下，最小水平主應(yīng)力方向(垂直水力裂縫方向)的誘導(dǎo)應(yīng)力相對(duì)變化較大，干擾強(qiáng)度增大，形成復(fù)雜縫或網(wǎng)絡(luò)縫的概率相對(duì)較大。

2) 裂縫內(nèi)摩擦角越小，天然裂縫網(wǎng)絡(luò)連通面積越大，進(jìn)行壓裂改造時(shí)越易形成復(fù)雜的網(wǎng)絡(luò)裂縫。

3) 裂縫性頁巖氣儲(chǔ)層在壓裂過程中形成的水力裂縫不是一條單一的平面縫，而是一種復(fù)雜的網(wǎng)狀裂縫。水平地應(yīng)力差在一定范圍內(nèi)，凈壓力系數(shù)越大，裂縫的擴(kuò)展形態(tài)越復(fù)雜，相鄰裂縫的尖端越容易連通形成網(wǎng)絡(luò)裂縫。

References

[1]吳奇，胥云，王騰飛，等.增產(chǎn)改造理念的重大變革：體積改造技術(shù)概論[J].天然氣工業(yè)，2011，31(4)：7-12,16.

WU Qi,XU Yun,WANG Tengfei,et al.The revolution of reservoir stimulation:an introduction of volume fracturing[J].Natural Gas Industry,2011,31(4):7-12,16.

[2]JARIPATKE O A,CHONG K K,GRIESER W V,et al.A completions roadmap to shale-play development:a review of successful approaches toward shale-play stimulation in the last two decades[R].SPE 130369,2010.

[3]王歡，廖新維，趙曉亮，等.非常規(guī)油氣藏儲(chǔ)層體積改造模擬技術(shù)研究進(jìn)展[J].特種油氣藏，2014,21(2)：8-15.

WANG Huan,LIAO Xinwei,ZHAO Xiaoliang,et al.The progress of reservoir stimulation technology in unconventional oil and gas reservoir[J].Special Oil & Gas Reservoirs,2014,21(2):8-15.

[4]衡帥，楊春和，郭印同，等.層理對(duì)頁巖水力裂縫擴(kuò)展的影響研究[J].巖石力學(xué)與工程學(xué)報(bào)，2015，34(2)：228-237.

HENG Shuai,YANG Chunhe,GUO Yintong,et al.Influence of bedding planes on hydraulic fracture propagation in shale formations[J].Chinese Journal of Rock Mechanics and Engineering,2015,34(2):228-237.

[5]蔣廷學(xué)，賈長貴，王海濤，等.頁巖氣網(wǎng)絡(luò)壓裂設(shè)計(jì)方法研究[J].石油鉆探技術(shù)，2011，39(3)：36-40.

JIANG Tingxue,JIA Changgui,WANG Haitao,et al.Study on network fracturing design method in shale gas[J].Petroleum Drilling Techniques,2011,39(3):36-40.

[6]劉雨，艾池.多級(jí)壓裂誘導(dǎo)應(yīng)力作用下天然裂縫開啟規(guī)律研究[J].石油鉆探技術(shù)，2015，43(1)：20-26.

LIU Yu,AI Chi.Opening of natural fractures under induced stress inmulti-stage fracturing[J].Petroleum Drilling Techniques,2015,43(1):20-26.

[7]蔣廷學(xué).頁巖油氣水平井壓裂裂縫復(fù)雜性指數(shù)研究及應(yīng)用展望[J].石油鉆探技術(shù)，2013，41(2)：7-12.

JIANG Tingxue.The fracture complexity index of horizontal wells in shale oil and gas reservoirs[J].Petroleum Drilling Techniques,2013,41(2):7-12.

[8]BLANTON T L.Propagation of hydraulically and dynamically induced fractures in naturally fractured reservoirs[R].SPE 15261,1986.

[9]DERSHOWITZ W S,COTTRELL M G,LIM D H,et al.A discrete fracture network approach for evaluation of hydraulic fracture stimulation of naturally fractured reservoirs[R].ARMA10-475,2010.

[10]RIAHI A,DAMJANAC B.Numerical study of the interaction between injection and the discrete fracture network in enhanced geothermal reservoirs[R].ARMA-2013-333,2013.

[11]COTTRELL M,HOSSEINPOUR H,DERSHOWITZ W.Rapid discrete fracture analysis of hydraulic fracture development in naturally fractured reservoirs[R].SPE 168843,2013.

[12]McCLURE M,BABAZADEH M,SHIOZAWA S,et al.Fully coupled hydromechanical simulation of hydraulic fracturing in three-dimensional discrete fracture networks[R].SPE 173354,2015.

[13]NAGEL N B,SANCHEZ-NAGEL M A,GARCIA X,et al.A numerical evaluation of the geomechanical interactions between a hydraulic fracture stimulation and a natural fracture system[R].ARMA-2012-287,2012.

[14]NAGEL N B,SANCHEZ-NAGEL M,LEE B,et al.Hydraulic fracturing optimization for unconventional reservoirs-the critical role of the mechanical properties of the natural fracture network[R].SPE 161934,2012.

[15]Itasca Consulting Group Inc.Universal distinct element code user’s guide[M].4th ed.Minneapolis:Itasca Consulting Group Inc,2000.

[16]CUNDALL P A,HART R D.Analysis of block test No.1.inelastic rock mass behavior:phase 2-a characterization of joint behavior (final report)[R].[S.l.]:[s.n.],1984.

[17]BARTON N,BANDIS S,BAKHTAR K.Strength,deformation and conductivity coupling of rock joints[J].International Journal of Rock Mechanics and Mining Science & Geomechanics Abstracts,1985,22(3):121-140.

[18]NAGEL N B,SANCHEZ M A,LEE B.Gas shale hydraulic fracturing:a numerical evaluation of the effect of geomechanical parameters[R].SPE 152192,2012.

[19]OLSON J E,TALEGHANI A D.Modeling simultaneous growth of multiple hydraulic fractures and their interaction with natural fractures[R].SPE 119739,2009.

[編輯劉文臣]

Numerical Simulation Study on the Complex Network Fractures of Stratified Shale Gas Reservoirs

LI Yumei1, LYU Wei2, SONG Jie2, LI Jun1, YANG Hongwei1, YU Liwei3

(1.CollegeofPetroleumEngineering,ChinaUniversityofPetroleum(Beijing),Beijing，102249，China; 2.ExplorationandDevelopmentInstitute,PetroChinaYumenOilfieldCompany,Jiuquan,Gansu， 735019,China; 3.EngineeringTechnologyInstitute,PetroChinaXinjiangOilfieldCompany,Karamay,Xinjiang，834000，China)

To perform the numerical simulation of complex network fractures in stratified shale gas reservoirs with natural fractures, a numerical model for network fractures was set up to determine interactions between hydraulic fractures and natural fractures by using universal distinct element code (UDEC) of non-continuous medium simulation based on numerical algorithms of seepage-stress coupling. The model was used to analyze the effect of the hydraulic fracture length, the natural fracture dip angle, the internal friction angle and the net treatment pressures on fracture network extension. The results showed that hydraulic fractures extend from the fracture tip near the well bore, and they experience shear failure along the strike of natural fractures. In the process, the connected area of natural fracture network would increase with the increases of hydraulic fractures length. It is more likely to form complex fractures or network fractures in the case of relatively large natural fracture dip angles. When the internal friction angle of the natural fractures is reduced, the connected area of natural fractures will increase and it is easier to form complex network fractures. Under a certain range of horizontal stress, fracture growth patterns would become more complex and it would be easier for tips near fractures to form network fractures as net pressure coefficients increase. Numerical simulation results could provide guidance for further understanding on shale gas fracturing fracture extension mechanism in areas away from wells.

shale gas; fracturing; numerical simulation; discrete element; geometric model; natural fracture; internal friction angle

2015-07-31；改回日期：2016-04-18。

李玉梅(1981—)，女，山東鄆城人，2007年畢業(yè)于大慶石油學(xué)院金屬材料工程專業(yè)，2010年獲中國石油大學(xué)(北京)油氣井工程專業(yè)碩士學(xué)位，在讀博士研究生，主要從事油氣井巖石力學(xué)與工程以及數(shù)值模擬計(jì)算方法研究。E-mail:liyumei3680238@163.com。

中國石油大學(xué)(北京)前瞻導(dǎo)向項(xiàng)目專項(xiàng)基金“各向異性頁巖氣藏分簇射孔參數(shù)優(yōu)化方法研究”(ZX20150039)資助。

doi:10.11911/syztjs.201604019

TE377

A

1001-0890(2016)04-0108-06

?油氣開發(fā)?