李 揚(yáng)， 鄧金根， 蔚寶華， 劉 偉， 陳建國

(油氣資源與探測國家重點(diǎn)實(shí)驗(yàn)室(中國石油大學(xué)(北京))，北京 102249)

水力壓裂在非常規(guī)油氣資源開發(fā)中具有重要作用，但該技術(shù)的成功率并不高，其中一個重要原因是最大施工壓力預(yù)測不準(zhǔn)確，缺少對隔層中壓裂縫高的有效控制方法。國內(nèi)外許多學(xué)者就地層彈性模量、地應(yīng)力、斷裂韌性以及界面性質(zhì)等因素對壓裂縫高的影響進(jìn)行了研究：Liu He等人[1]研究了孔隙度和泥質(zhì)含量對水力壓裂縫高的影響；Gu Hongren等人[2]通過三維模擬試驗(yàn)得出，除地應(yīng)力之外，彈性模量對壓裂縫高的影響也較為明顯，當(dāng)裂縫未進(jìn)入隔層時，彈性模量大的隔層對縫高擴(kuò)展起阻礙作用，當(dāng)裂縫已擴(kuò)展進(jìn)入隔層時，彈性模量小的隔層對縫高限制效果較好；M.B.Smith等人[3]認(rèn)為僅憑儲/隔層彈性模量之間的差異不足以影響裂縫縫高的擴(kuò)展；A.A.Daneshy和R.D.Barree等人[4-7]認(rèn)為儲/隔層界面之間的剪切滑移是影響縫高的重要因素之一；G.Anderson[8]通過試驗(yàn)證實(shí)了界面之間存在剪切滑移的可能；黃榮樽[9]指出界面抗剪強(qiáng)度一方面可能源于層間膠結(jié)物的粘合力，另一方面可能來源于抵抗變形的摩擦力產(chǎn)生的機(jī)械聯(lián)結(jié)力；周文高等人[10]研究了巖石物質(zhì)特性、施工參數(shù)、地應(yīng)力差和裂縫上下端阻抗對縫高的影響；M.K.Fisher等人[11]根據(jù)實(shí)際壓裂作業(yè)的微地震數(shù)據(jù)，從地應(yīng)力、巖性、弱面、斷層、壓裂液梯度和巖石滲透性等方面討論了壓裂縫縫長和縫高的擴(kuò)展；陳治喜等人[12]認(rèn)為在一定的地應(yīng)力條件下，裂縫能否向隔層擴(kuò)展以及擴(kuò)展范圍的大小取決于作業(yè)壓力的大?。煌蹂热薣13]利用二維平面應(yīng)變非線性流固耦合水力壓裂模型對裂縫高度的擴(kuò)展進(jìn)行了模擬，極大地減小了計(jì)算量?？傮w而言，現(xiàn)有研究成果主要從巖石力學(xué)特性和地應(yīng)力差兩方面來評價(jià)隔層控制縫高的能力，但是就部分巖石力學(xué)性質(zhì)對縫高的影響而言還存在爭議，而且影響裂縫縫高的因素特別多，各參數(shù)的影響程度并不明確，儲/隔層界面性質(zhì)對裂縫縫高的影響缺乏有效的評價(jià)方法。鑒于此，筆者基于有限元軟件ABAQUS研究了儲/隔層巖石彈性模量、地應(yīng)力和抗拉強(qiáng)度等因素對水力裂縫縫高擴(kuò)展的影響，并著重分析了儲/隔層界面性質(zhì)的影響。

1 Cohesive單元損傷模型

有限元計(jì)算平臺ABAQUS中的Cohesive單元基于損傷力學(xué)原理、采用剛度衰減的方法模擬介質(zhì)中裂縫的起裂與擴(kuò)展，并且可以模擬裂縫擴(kuò)展過程中壓裂液的流動與濾失，比較適合用來模擬水力壓裂裂縫的擴(kuò)展。

1.1 Cohesive單元損傷模式

圖1 Cohesive單元損傷T-S準(zhǔn)則Fig.1 The traction-separation law of Cohesive element damage

該T-S準(zhǔn)則是以Cohesive單元承受的拉應(yīng)力為損傷判據(jù)，單元在損傷之前其剛度非常大，所受應(yīng)力與位移成正比，卸載后可以恢復(fù)，但當(dāng)拉應(yīng)力超過材料抗拉強(qiáng)度之后，其能承受的應(yīng)力隨著位移的增加呈線性衰減，并且不可恢復(fù)。該準(zhǔn)則常用來模擬巖石、混凝土等材料中裂縫的擴(kuò)展，比較適合于模擬水力壓裂中裂縫的擴(kuò)展過程。

1.2 Cohesive單元的起裂與擴(kuò)展準(zhǔn)則

筆者所采用的Cohesive單元的起裂準(zhǔn)則為二次應(yīng)力準(zhǔn)則，其表達(dá)式為[15]：

(1)

Cohesive單元采用剛度退化來描述單元的損傷演化過程，其表達(dá)式為：

(2)

損傷因子的計(jì)算公式為[16]：

(3)

1.3 Cohesive單元損壞區(qū)的流體流動

Cohesive單元損壞區(qū)的流體流動分為沿Cohesive單元的切向流動和垂直于Cohesive單元上下表面的法向?yàn)V失，如圖2所示。

圖2 Cohesive單元損壞區(qū)流體流動示意Fig.2 The fluid flow schematic within a damaged Cohesive zone

假設(shè)Cohesive單元內(nèi)部流體為不可壓縮的牛頓流體，其切向流動的計(jì)算式為[17]：

(4)

Cohesive單元上下表面法向方向的濾失可描述為[18]：

(5)

式中：pt，pb分別為裂縫上下表面處的孔隙壓力，Pa；pf為裂縫內(nèi)的流體壓力，Pa；ct，cb分別為上下表面的濾失系數(shù)，m3/(Pa·s)；qt，qb分別為上下表面的法向體積流量，m3/s。

Cohesive單元的流體質(zhì)量守恒方程為[19]：

(6)

式中：Q(t)為壓裂流體注入速度，m3/s。

2 計(jì)算模型及結(jié)果分析

為了研究各因素對縫高的影響，根據(jù)界面的膠結(jié)程度分了2種情況：一種界面膠結(jié)完好，界面不產(chǎn)生滑移；另一種界面膠結(jié)不完全，儲/隔層之間可能會產(chǎn)生滑移。假定井為垂直井，地層為水平層狀均質(zhì)地層，建立了如圖3所示的計(jì)算模型。該模型中，豎直方向設(shè)置Cohesive單元模擬水力裂縫的垂向擴(kuò)展，水平方向設(shè)置Cohesive單元層間界面性質(zhì)，水平Cohesive單元可以模擬界面的拉伸破壞和剪切破壞。

圖3 計(jì)算模型示意Fig.3 Sketch of computational model

2.1 界面膠結(jié)完好條件下的情況

假定產(chǎn)層與隔層之間膠結(jié)完好，界面不會產(chǎn)生滑移，討論隔層和儲層的彈性模量、地應(yīng)力以及抗拉強(qiáng)度等參數(shù)對縫高的影響(計(jì)算模型見圖3)。模型儲層厚度5 m，隔層厚度30 m，寬度50 m。左右兩側(cè)分層施加非均勻最小水平主應(yīng)力，上邊界施加上覆巖層壓力，下邊界設(shè)置為關(guān)于y軸的對稱邊界條件，左右和上邊界保持恒定的孔隙壓力(34 MPa)，初始孔隙壓力為34 MPa，初始孔隙度為0.15。由于模型為平面應(yīng)變，因此模型中的泵注排量根據(jù)現(xiàn)場排量2～3 m3/min進(jìn)行相應(yīng)的折減，設(shè)置為0.012 m3/min，持續(xù)600 s，壓裂液的黏度為100 mPa·s，其他的主要計(jì)算參數(shù)初始值見表1和表2。

表1 地層巖石力學(xué)參數(shù)Table 1 The rock mechanics parameters of the stratum

表2 Cohesive單元材料參數(shù)Table 2 The material parameters of Cohesive elements

2.1.1 彈性模量對壓裂縫高的影響

儲層巖石的彈性模量取30 GPa，隔層的彈性模量分別取20、30、40和50 GPa，其他材料參數(shù)值見表1和表2，研究隔層巖石彈性模量對裂縫縫高的影響，結(jié)果如圖4所示。從圖4可以看出，隔層彈性模量越大，裂縫在隔層擴(kuò)展的高度越大。分析認(rèn)為，這是因?yàn)閺椥阅Ａ看髸沽芽p的寬度變窄，在注入流體體積和濾失量一定的情況下，裂縫高度就會相應(yīng)增大。同時，巖石的彈性模量越大，越容易在裂縫尖端形成高應(yīng)力，在同樣的抗拉強(qiáng)度條件下巖石更容易向前擴(kuò)展。因此，彈性模量小的隔層對控制壓裂縫高有利。

圖4 不同隔層彈性模量下的裂縫形態(tài)Fig.4 Fracture shapes with different elastic modulus of insulating layers

圖5為壓裂過程中井底壓力隨時間的變化曲線，并標(biāo)出了裂縫擴(kuò)展穿過儲/隔層界面時的井底壓力。

圖5 不同隔層彈性模量下的井底壓力Fig.5 Bottom hole pressure with different elastic modulus of insulating layers

從圖5可以看出，隨隔層彈性模量增大，裂縫從儲層進(jìn)入隔層時刻的井底壓力會略微增大，說明隔層彈性模量大對裂縫擴(kuò)展進(jìn)入隔層確實(shí)有阻礙作用，但是阻礙效果不明顯，和現(xiàn)場施工過程中的壓力波動相比幾乎可以忽略。

綜合分析圖4、圖5可知，當(dāng)泵入排量不足以使裂縫擴(kuò)展至分界面時，隔層彈性模量大對縫高延伸起限制作用；而一旦裂縫已經(jīng)進(jìn)入隔層，彈性模量大的隔層其縫寬相對較小，由壓裂液體積守恒可知，縫高會急劇增大，壓穿隔層進(jìn)入水層的風(fēng)險(xiǎn)增大。

2.1.2 地應(yīng)力對壓裂縫高的影響

計(jì)算參數(shù)值取自表1和表2，逐漸改變隔層巖石的最小水平主應(yīng)力(分別為50、52、54、56和58 MPa)，研究隔層最小水平主應(yīng)力對裂縫縫高的影響，結(jié)果如圖6所示。從圖6可以看出，隔層的最小水平主應(yīng)力越大，隔層中裂縫的寬度和高度越小，儲層中裂縫的寬度越大。分析認(rèn)為，這是因?yàn)榱芽p寬度與凈壓力相關(guān)，隔層的最小水平主應(yīng)力越大，那么儲層中的凈壓力就相對越大，壓裂液就能更多地停留在儲層中。

圖6 不同隔層地應(yīng)力下的裂縫形態(tài)Fig.6 Fracture shapes with different in-situ stress of insulating layers

2.1.3 抗拉強(qiáng)度差對壓裂縫高的影響

計(jì)算參數(shù)值取自表1和表2，逐漸改變隔層巖石的抗拉強(qiáng)度(分別為4、7、10、13和16 MPa)，研究隔層抗拉強(qiáng)度對裂縫縫高的影響，結(jié)果如圖7所示。

圖7 不同隔層抗拉強(qiáng)度下的裂縫形態(tài)Fig.7 Fracture shapes with different tensile strength of insulating layers

從圖7可以看出，隔層的抗拉強(qiáng)度越大，裂縫的高度就越小。分析認(rèn)為，這是因?yàn)榱芽p必須克服巖石的抗拉強(qiáng)度才能繼續(xù)向前擴(kuò)展。因此，高抗拉強(qiáng)度的隔層也會限制裂縫縫高的擴(kuò)展。

通過分析界面完全膠結(jié)時儲/隔層彈性模量、地應(yīng)力和抗拉強(qiáng)度可知，彈性模量對水力裂縫的起裂影響較小，但是對裂縫擴(kuò)展的影響較大，地應(yīng)力和抗拉強(qiáng)度對裂縫起裂與擴(kuò)展的影響均較大。總體來說，地應(yīng)力大、抗拉強(qiáng)度高和彈性模量小的隔層有利于控制縫高。

2.2 界面不完全膠結(jié)條件下的情況

假設(shè)儲層和隔層界面不完全膠結(jié)，分層界面之間存在著一定量的膠結(jié)物，并且在上覆巖層壓力作用下有一定程度的壓實(shí)，能夠承受一定的拉應(yīng)力和剪應(yīng)力。在該情況下裂縫能否穿過界面進(jìn)入隔層受到多種因素的影響，其中包括層間地應(yīng)力、力學(xué)參數(shù)以及界面自身的性質(zhì)。關(guān)于層間地應(yīng)力和力學(xué)參數(shù)對裂縫縫高的影響，已經(jīng)在前文討論，以下主要討論層間界面性質(zhì)對壓裂縫高的影響。

2.2.1 界面抗剪強(qiáng)度對裂縫高度的影響

儲層和隔層界面抗剪強(qiáng)度應(yīng)廣義地理解為2部分：一是界面膠結(jié)物本身的抗剪強(qiáng)度，二是在上覆巖層壓力作用下由于接觸界面比較粗糙引起的摩擦強(qiáng)度。對于儲/隔層界面的廣義抗剪強(qiáng)度，可以依據(jù)弱面理論表示為[20]：

τ=Cw+σtanφw

(7)

式中：τ為結(jié)構(gòu)面的抗剪強(qiáng)度，Pa；Cw為結(jié)構(gòu)面內(nèi)聚力，Pa；σ為結(jié)構(gòu)面上的法向應(yīng)力，Pa；φw為結(jié)構(gòu)面有效摩擦角，(°)。

文獻(xiàn)[21]給出了各類軟弱結(jié)構(gòu)面的抗剪強(qiáng)度，根據(jù)測井解釋的巖層性質(zhì)可以查出對應(yīng)結(jié)構(gòu)面的內(nèi)聚力和有效摩擦角，從而能確定界面抗剪強(qiáng)度。

水力壓裂過程中通常有以下2種情況：一是層間界面抗剪強(qiáng)度小于隔層抗拉強(qiáng)度與最小水平主應(yīng)力之和，此時水力裂縫將會完全被限制在產(chǎn)層里面，界面處產(chǎn)生滑移；二是當(dāng)層間界面抗剪強(qiáng)度足夠大，界面處會發(fā)生部分滑移，但還是會擴(kuò)展進(jìn)入隔層。

計(jì)算模型見圖3，計(jì)算參數(shù)見表3。

表3 地層巖石力學(xué)參數(shù)Table 3 The rock mechanics parameters of the stratum

分別研究了界面抗剪強(qiáng)度為50、58和66 MPa時不同注入時間下裂縫的形態(tài)演化情況，并以裂紋中心為(0，0)點(diǎn)，繪制了裂縫縫高的幾何形態(tài)，見圖8—11。

圖8 界面抗剪強(qiáng)度為50 MPa時不同注入時間下裂縫形態(tài)的演化Fig.8 Evolution of fracture shapes with 50 MPa of interface shear strength at different injection times

圖9 界面抗剪強(qiáng)度為58 MPa時不同注入時間下裂縫形態(tài)的演化Fig.9 Evolution of fracture shapes with 58 MPa of interface shear strength at different injection times

圖10 界面抗剪強(qiáng)度為66 MPa時不同注入時間下裂縫形態(tài)的演化Fig.10 Evolution of fracture shapes with 66 MPa of interface shear strength at different injection times

從圖8可以看出，當(dāng)層間界面抗剪強(qiáng)度為50 MPa時，產(chǎn)層在層間界面處完全滑移，隔層不會起裂。從圖9可以看出，當(dāng)層間界面抗剪強(qiáng)度為58 MPa時，大于隔層最小水平主應(yīng)力與抗拉強(qiáng)度之和，產(chǎn)層在層間界面處部分產(chǎn)生滑移，同時裂縫也擴(kuò)展至隔層。

圖11 不同界面抗剪強(qiáng)度下的裂縫形態(tài)Fig.11 Fracture shapes with different interface shear strength

從圖10可以看出，當(dāng)層間界面抗剪強(qiáng)度為66 MPa時，產(chǎn)層在層間界面處產(chǎn)生部分滑移，裂縫也擴(kuò)展至隔層。比較注入時間為4和6 min時的裂縫形態(tài)可以看出，當(dāng)裂縫進(jìn)入隔層之后，儲層的裂縫寬度會有一個減小的過程，這是因?yàn)榱芽p擴(kuò)展進(jìn)入高應(yīng)力的隔層后，儲層內(nèi)壓裂液壓力會有一個釋放的過程。

從圖11可以看出，層間界面抗剪強(qiáng)度越大，水力裂縫進(jìn)入隔層的可能性越大，層間滑移越小，同時裂縫上竄的高度越大。

水力裂縫擴(kuò)展至儲/隔層交界面后，通常情況下只有當(dāng)界面抗剪強(qiáng)度大于隔層最小水平主應(yīng)力與隔層抗拉強(qiáng)度之和這一臨界值時，裂縫才會擴(kuò)展至隔層，但是水力裂縫的縫長方向通常延伸得較遠(yuǎn)，難免有某些地方的隔層本身存在天然微裂隙，此時若有壓裂液滲入其中，即使界面的抗剪強(qiáng)度很低，水力裂縫也有可能擴(kuò)展進(jìn)入隔層。

2.2.2 界面上覆巖層壓力對裂縫高度的影響

當(dāng)上覆巖層壓力遠(yuǎn)大于最大和最小水平主應(yīng)力時，難以形成水平裂縫或T形裂縫。但對于逆斷層發(fā)育和淺層等構(gòu)造運(yùn)動相對強(qiáng)烈的地區(qū)，上覆巖層壓力小于或相近于水平地應(yīng)力，則非常容易形成水平裂縫。因此，對于深部地層，可能存在層間水平地應(yīng)力相差較大的情況，比如上覆巖層壓力在儲層是最大主應(yīng)力，但是在隔層可能是中間主應(yīng)力。這種情況下，當(dāng)層間界面膠結(jié)強(qiáng)度較低時，裂縫容易在界面處轉(zhuǎn)向，沿著界面擴(kuò)展，形成T形裂縫，限制裂縫高度擴(kuò)展。可以用如下關(guān)系式來判定能否形成此類T形裂縫：

σhb+Stb>σv+St>σhp+Stp

(8)

式中：σhb為隔層最小水平主應(yīng)力，MPa；Stb為隔層抗拉強(qiáng)度，MPa；σhp為儲層最小水平主應(yīng)力，MPa；Stp為儲層抗拉強(qiáng)度，MPa；σv為界面處的上覆巖層壓力，MPa；St為界面抗拉強(qiáng)度，MPa。

當(dāng)σv+St<σhp+Stp時，會直接在儲層產(chǎn)生水平裂縫；當(dāng)σv+St>σhb+Stb時，裂縫會直接穿透界面進(jìn)入隔層或者在儲/隔層界面產(chǎn)生滑移；當(dāng)σhb+Stb>σv+St>σhp+Stp時，在儲層形成垂直裂縫，并沿層間界面轉(zhuǎn)向形成水平裂縫。

為驗(yàn)證此準(zhǔn)則的正確性，利用圖3所示模型進(jìn)行了2組數(shù)值模擬試驗(yàn)，計(jì)算參數(shù)值見表4，其他參數(shù)值與表3相同。

表4 數(shù)值試驗(yàn)參數(shù)Table 4 The parameters of numerical experiments

數(shù)值試驗(yàn)1和2分別滿足關(guān)系式(9)和(10):

σv+St=58>σhb+Stb=57.8

(9)

σv+St=57<σhb+Stb=57.8

(10)

對上述2種情況分別進(jìn)行數(shù)值計(jì)算，結(jié)果如圖12和圖13所示。

圖12 數(shù)值模擬試驗(yàn)1結(jié)果Fig.12 The result of the numerical experiment 1

圖13 數(shù)值模擬試驗(yàn)2結(jié)果Fig.13 The result of the numerical experiment 2

由圖12和圖13可知，當(dāng)?shù)貙訔l件滿足式(8)時，水力裂縫擴(kuò)展至儲/隔層界面后會發(fā)生轉(zhuǎn)向形成T形縫。

3 結(jié) 論

1) 彈性模量小、地應(yīng)力大和抗拉強(qiáng)度高的隔層對縫高能夠起到一定的限制作用，其中地應(yīng)力差的限制作用最為明顯。在裂縫進(jìn)入隔層之前，彈性模量大的隔層可以對裂縫進(jìn)入隔層起到一定的阻礙作用，但裂縫進(jìn)入隔層后彈性模量小的隔層對縫高有較好的限制作用。彈性模量對裂縫的起裂影響較小，僅憑儲/隔層之間彈性模量的差異難以將裂縫完全限制在儲層內(nèi)。

2) 除彈性模量、地應(yīng)力外，儲/隔層界面滑移也是影響壓裂縫高的主要因素之一，當(dāng)界面抗剪強(qiáng)度低于某一臨界值時，裂縫尖端可能在界面處產(chǎn)生滑移，使裂縫上端變鈍，能有效限制裂縫沿垂向擴(kuò)展。界面的抗剪強(qiáng)度是水力壓裂作業(yè)裂縫能否進(jìn)入隔層的重要判斷指標(biāo)之一，其中儲/隔層界面的抗剪強(qiáng)度可以用弱面準(zhǔn)則來估算。

3) 上覆巖層壓力對裂縫的形態(tài)也存在一定的影響，在特定上覆巖層壓力等地層條件下可能產(chǎn)生T形縫。給出了在逆斷層控制區(qū)域產(chǎn)生T形縫的判斷準(zhǔn)則，并用數(shù)值模擬試驗(yàn)對其進(jìn)行了驗(yàn)證。

參考文獻(xiàn)
References

[1] Liu He,Wang Han,Wu Heng’an,et al.Effect of reservoir porosity and clay content on hydraulic fracture height containment[R].IPTC 16415,2013.

[2] Gu Hongren,Siebrits E.Effect of formation modulus contrast on hydraulic fracture height containment[R].SPE 103822,2006.

[3] Smith M B,Bale A B,Britt L K,et al.Layered modulus effects on fracture propagation,proppant placement,and fracture modeling[R].SPE 71654,2001.

[4] Daneshy A A.Factors controlling the vertical growth of hydraulic fractures[R].SPE 118789,2009.

[5] Daneshy A A.Hydraulic fracture propagation in layered formations[J].SPE Journal,1978,18(1):33-41.

[6] Daneshy A A.Hydraulic fracture propagation in the presence of planes of weakness[R].SPE 4852,1974.

[7] Barree R D,Winterfeld P H.Effects of shear planes and interfacial slippage on fracture growth and treating pressures[R].SPE 48926,1998.

[8] Anderson G.Effects of friction on hydraulic fracture growth near unbonded interfaces in rocks[J].SPE Journal,1981,21(1):21-29.

[9] 黃榮樽.水力壓裂裂縫的起裂和擴(kuò)展[J].石油勘探與開發(fā)，1981，8(5):62-74.

Huang Rongzun.The initiation and propagation of hydraulic fracture[J].Petroleum Exploration and Development，1981，8

(5):62-74.

[10] 周文高，胡永全，趙金洲，等.控制壓裂縫高技術(shù)研究及影響因素分析[J].斷塊油氣田，2006，13(4):70-72.

Zhou Wengao,Hu Yongquan,Zhao Jinzhou,et al.Research on fracture height containment technology and analysis of influence factors[J].Fault-Block Oil & Gas Field，2006,13(4):70-72.

[11] Fisher M K,Warpinski N R.Hydraulic-fracture-height growth:real data[J].SPE Production & Operations,2012,27(1):8-19.

[12] 陳治喜,陳勉,黃榮樽，等.層狀介質(zhì)中水力裂縫的垂向擴(kuò)展[J].石油大學(xué)學(xué)報(bào):自然科學(xué)版,1997，21(4):23-26，32.

Chen Zhixi,Chen Mian，Huang Rongzun,et al.Vertical growth of hydraulic fracture in layered formations[J].Journal of the University of Petroleum,China：Edition of Natural Science,1997,21(4):23-26,32.

[13] 王瀚,劉合,張勁,等.水力裂縫的縫高控制參數(shù)影響數(shù)值模擬研究[J].中國科學(xué)技術(shù)大學(xué)學(xué)報(bào),2011，41(9):820-825.

Wang Han,Liu He,Zhang Jin,et al.Numerical simulation of hydraulic fracture height control with different parameters[J].Journal of University of Science and Technology of China，2011,41(9):820-825.

[14] Tomar V,Zhai Jun,Zhou Min.Bounds for element size in a variable stiffness cohesive finite element model[J].International Journal for Numerical Methods in Engineering,2004,61(11):1894-1920.

[15] Camanho P P,Dávila C G.Mixed-mode decohesion finite elements for the simulation of delamination in composite materials[R].NASA/TM-2002-211737,2002.

[16] Turon A,Camanho P P,Costa J,et al.A damage model for the simulation of delamination in advanced composites under variable-mode loading[J].Mechanics of Materials,2006,38(11):1072-1089.

[17] Economides M J,Nolte K G,Ahmed U.Reservoir stimulation[M].Chichester:Wiley，2000:5.15-5.16.

[18] Hagoort J,Weatherill B,Settari A.Modeling the propagation of waterflood-induced hydraulic fractures[J].SPE Journal,1980,20(4):293-303.

[19] Peirce A,Detournay E.An implicit level set method for modeling hydraulically driven fractures[J].Computer Methods in Applied Mechanics and Engineering，2008,197(33-40):2858-2885.

[20] Fjar E,Holt R M,Raaen A M,et al.Petroleum related rock mechanics[M].2nd ed.Amsterdam:Elsevier,2008:60-64.

[21] 賈喜榮.巖石力學(xué)與巖層控制[M].徐州：中國礦業(yè)大學(xué)出版社，2010:111-114.

Jia Xirong.Rock mechanics and strata control[M].Xuzhou:China University of Mining and Technology Press,2010:111-114.