徐 珂，戴俊生，付曉龍，任啟強(qiáng)，劉叢寧，趙 惲

（1.中國(guó)石油大學(xué)地球科學(xué)與技術(shù)學(xué)院，山東青島 266580；2.中國(guó)石油冀東油田勘探開(kāi)發(fā)研究院，河北唐山 063004；

3.廣州大學(xué)教育學(xué)院，廣州 510006）

基于有限元法的層狀巖體破裂規(guī)律探討

徐 珂1，戴俊生1，付曉龍1，任啟強(qiáng)1，劉叢寧2，趙 惲3

（1.中國(guó)石油大學(xué)地球科學(xué)與技術(shù)學(xué)院，山東青島 266580；2.中國(guó)石油冀東油田勘探開(kāi)發(fā)研究院，河北唐山 063004；

3.廣州大學(xué)教育學(xué)院，廣州 510006）

層狀巖體的非均質(zhì)性及各向異性導(dǎo)致其破裂方式及規(guī)律與均質(zhì)巖體有顯著不同。對(duì)層狀巖體分別進(jìn)行不同方式的單軸、雙軸、三軸試驗(yàn)，分析應(yīng)力-應(yīng)變曲線(xiàn)特征；再利用ANSYS有限元軟件進(jìn)行數(shù)值模擬，觀(guān)察應(yīng)力、應(yīng)變?cè)趲r體上的分布，通過(guò)曲線(xiàn)和圖件的對(duì)比分析，并結(jié)合巖石破裂理論，總結(jié)不同應(yīng)力狀態(tài)下層狀巖體的破裂方式、順序及規(guī)律；最后以富臺(tái)地區(qū)為例，對(duì)分析結(jié)果進(jìn)行驗(yàn)證。研究結(jié)果表明，不同受力方式對(duì)層狀巖體破裂的影響體現(xiàn)在施加的載荷及約束與層面的方位。當(dāng)應(yīng)力方向與巖層面平行時(shí)，強(qiáng)度大的石灰?guī)r巖體發(fā)生集中應(yīng)力，首先破裂；而應(yīng)力與巖層面垂直時(shí)，強(qiáng)度小的泥巖巖體首先破裂。巖石試驗(yàn)、數(shù)值模擬結(jié)果以及實(shí)例均成功驗(yàn)證了這個(gè)規(guī)律。

層狀巖體；巖石試驗(yàn)；數(shù)值模擬；破裂規(guī)律

0 引言

層狀巖體普遍分布于自然界中。與均質(zhì)塊狀巖體相比，層狀巖體的變形和強(qiáng)度特性具有明顯的各向異性，因此其穩(wěn)定性和破壞條件也有特殊表現(xiàn)，特別是裂縫的產(chǎn)生和延伸變得比較復(fù)雜，現(xiàn)有的破裂準(zhǔn)則也難以精準(zhǔn)描述它們的破裂規(guī)律。目前研究過(guò)層狀巖體力學(xué)性質(zhì)和變形方式的學(xué)者不在少數(shù)，胡明研等［1］對(duì)層狀復(fù)合巖石做了壓縮試驗(yàn)，分析了不同條件下裂縫在巖層上的發(fā)育和展布；黎立云等［2］通過(guò)巴西圓盤(pán)劈裂試驗(yàn)及單軸壓縮試驗(yàn)討論了層狀巖體破裂的特殊現(xiàn)象及成因機(jī)制；劉立［3］對(duì)復(fù)合巖體的特殊性質(zhì)進(jìn)行了深入研究，在損傷演化和斷裂破壞規(guī)律方面有了實(shí)質(zhì)進(jìn)展；許寶田等［4］對(duì)泥巖進(jìn)行了三軸壓縮實(shí)驗(yàn)，分析了泥巖應(yīng)力-應(yīng)變特性；戴俊生等［5］提出了適合脆性低滲透砂巖的破裂準(zhǔn)則，盧虎勝等［6］則在戴俊生的基礎(chǔ)上建立了用于砂泥巖間互地層的破裂準(zhǔn)則；雷宇［7］通過(guò)砂泥巖壓縮、回彈、破裂實(shí)驗(yàn)對(duì)裂縫成因進(jìn)行了初步分析；楊仕教等［8］、李永盛［9］和尹小濤等［10］分別研究了不同加載速率對(duì)巖石力學(xué)效應(yīng)的影響。

利用ANSYS軟件進(jìn)行有限元數(shù)值模擬，可以對(duì)不連續(xù)的非均勻介質(zhì)組成的巖體進(jìn)行應(yīng)力特征研究。將整個(gè)巖體分為數(shù)目有限的單元，通過(guò)分析計(jì)算每個(gè)單元的應(yīng)力場(chǎng)進(jìn)而綜合所有單元來(lái)研究整體特征。單元的位移、應(yīng)變及應(yīng)力都可以經(jīng)過(guò)計(jì)算而得出，因此巖體在受力狀態(tài)下的應(yīng)力、應(yīng)變及位移分布可以直觀(guān)地表現(xiàn)出來(lái)。

本文以巖石試驗(yàn)為基礎(chǔ)，并用ANSYS軟件建立與實(shí)體1：1的模型，約束與載荷按照真實(shí)情景加載，分析層狀巖體在不同應(yīng)力狀態(tài)下的應(yīng)力、應(yīng)變及強(qiáng)度分布，據(jù)此推斷和解釋巖體的破裂規(guī)律。研究成果有助于判斷裂縫的產(chǎn)生、發(fā)育與分布，對(duì)儲(chǔ)層裂縫帶的預(yù)測(cè)具有指導(dǎo)意義。

1 巖石力學(xué)試驗(yàn)

巖石受到不同構(gòu)造應(yīng)力作用時(shí)，將發(fā)生不同形式的變形。當(dāng)其所承受的應(yīng)力超過(guò)極限強(qiáng)度，就會(huì)發(fā)生破裂并產(chǎn)生裂縫。由于不同巖石本身性質(zhì)及所受應(yīng)力條件不同，破裂的方式也不同。對(duì)于層狀巖石而言，各分層不同的力學(xué)性質(zhì)必然導(dǎo)致破裂規(guī)律的特殊性，而不同的應(yīng)力狀態(tài)也會(huì)產(chǎn)生不同的破裂方式。例如，裂紋在均質(zhì)花崗巖上的擴(kuò)展比較平直（見(jiàn)圖1a），而在層狀花崗巖上，裂紋越過(guò)層面時(shí)擴(kuò)展路徑發(fā)生變化，呈鋸齒狀（見(jiàn)圖1b）［2］。

圖1 均質(zhì)巖石與層狀巖石的裂紋擴(kuò)展Fig.1 The crack growth in homogeneous rock and layered rock

本文根據(jù)相似原理［1］，利用與泥巖、砂巖和石灰?guī)r3種常見(jiàn)且性質(zhì)差異明顯的材料，制成復(fù)合成層狀巖體，分別為邊長(zhǎng)100 mm的正方體及直徑50 mm、高100 mm的圓柱體。通過(guò)對(duì)單一巖體進(jìn)行三軸試驗(yàn)確定其彈性模量和泊松比。計(jì)算公式如下：

式中：E為巖石彈性模量，GPa；μ為巖石泊松比；σ1、σ2分別為應(yīng)力-應(yīng)變曲線(xiàn)直線(xiàn)段開(kāi)始和結(jié)束時(shí)的應(yīng)力，MPa；εz1、εz2分別為應(yīng)力-應(yīng)變曲線(xiàn)直線(xiàn)段開(kāi)始和結(jié)束時(shí)的軸向應(yīng)變［11］，mm；εx1、εx2分別為應(yīng)力-應(yīng)變曲線(xiàn)直線(xiàn)段開(kāi)始和結(jié)束時(shí)的側(cè)向應(yīng)變，mm。

試驗(yàn)得到的參數(shù)如表1。

表1 巖石力學(xué)參數(shù)Table 1 Rock mechanical parameters

試驗(yàn)用的每件樣品應(yīng)完好無(wú)損，分別進(jìn)行單軸、雙軸、三軸壓縮試驗(yàn)，其中單軸和雙軸試驗(yàn)再分為巖層水平和豎直兩種情況。試驗(yàn)得到的應(yīng)力-應(yīng)變曲線(xiàn)如圖2。除三軸壓縮外，其余4曲線(xiàn)雖大致趨勢(shì)相似，但形態(tài)不盡相同，意味著不同受力方式對(duì)巖石破裂有重要影響。

圖2 不同應(yīng)力狀態(tài)下的應(yīng)力-應(yīng)變曲線(xiàn)Fig.2 The stress-strain curves in different load states

2 數(shù)值模擬與分析

模擬依據(jù)巖石力學(xué)試驗(yàn)的5種方式進(jìn)行，即巖層水平的單軸壓縮，巖層豎直的單軸壓縮，巖層水平的雙軸壓縮，巖層豎直的雙軸壓縮，三軸壓縮試驗(yàn)（見(jiàn)圖3）。

圖3 5種受力方式示意圖Fig.3 The sketch map of five stress modes

2.1 巖層水平的單軸壓縮

首先研究分層排列順序?qū)ζ屏训挠绊?。進(jìn)行3次試驗(yàn)，各層自上而下的排列順序分別為：①砂巖、泥巖、石灰?guī)r；②泥巖、砂巖、石灰?guī)r；③砂巖、石灰?guī)r、泥巖。

以試驗(yàn)①為例，隨著壓力增大，強(qiáng)度最弱的泥巖先出現(xiàn)裂紋，但不能很快貫通界面；繼續(xù)增加荷載，當(dāng)裂縫的擴(kuò)張力大于層間的束縛力時(shí)，砂巖和石灰?guī)r也開(kāi)始出現(xiàn)裂紋，且很快與泥巖的裂紋形成了貫通縫（見(jiàn)圖4）。

圖4 試驗(yàn)①裂紋發(fā)展與巖體破壞示意圖Fig.4 The crack development and rock failure process in Test①

試驗(yàn)②、③的結(jié)果與試驗(yàn)①類(lèi)似，總是強(qiáng)度最弱的泥巖首先破裂。因此，受力條件相同的情況下，破裂與巖層排列順序無(wú)明顯關(guān)系。

水平巖層單軸壓縮的應(yīng)力-應(yīng)變曲線(xiàn)（見(jiàn)圖2a）大體分4個(gè)階段：孔隙壓密階段，表現(xiàn)為曲線(xiàn)“上凹”；彈性變形階段，表現(xiàn)為近直線(xiàn)；破裂發(fā)展階段，其上界為峰值強(qiáng)度；破裂后階段，曲線(xiàn)迅速下降。

使用ANSYS 15.0軟件進(jìn)行數(shù)值模擬。模型按實(shí)際試件1∶1比例建立，層間為Glue方式膠結(jié)。采用Solid45單元進(jìn)行網(wǎng)格劃分，每個(gè)單元為8節(jié)點(diǎn)，z方向施加擠壓力。為了看出應(yīng)力-應(yīng)變分布的顯著變化，擠壓力設(shè)定為15 MPa。模擬計(jì)算結(jié)果見(jiàn)圖5。

圖5 巖層水平單軸壓縮數(shù)值模擬結(jié)果Fig.5 The numerical simulation results of uniaxial compression experiment with horizontal layers

3個(gè)主應(yīng)力中，最小主應(yīng)力和中間主應(yīng)力均平行于層面，最大主應(yīng)力豎直向下。應(yīng)力分布特征：泥巖和石灰?guī)r的最小主應(yīng)力和中間主應(yīng)力為擠壓應(yīng)力，泥巖應(yīng)力值低于石灰?guī)r，而砂巖的最小主應(yīng)力和中間主應(yīng)力呈現(xiàn)部分拉張應(yīng)力。這是由于泥巖和砂巖交界處的水平變形不一致，而砂巖的彈性模量大于泥巖，且泊松比小于泥巖，導(dǎo)致泥巖的水平變形大于砂巖。在巖體受壓至達(dá)到峰值強(qiáng)度的過(guò)程中，由于變形特性的差異及變形連續(xù)靜力平衡條件，交界面附近的砂巖和泥巖在水平方向上必然發(fā)生應(yīng)力突變，就會(huì)在界面上產(chǎn)生方向相反的應(yīng)力［12］。正是這個(gè)原因，在交界處附近的石灰?guī)r層上也出現(xiàn)局部拉張力。

應(yīng)力強(qiáng)度分布圖說(shuō)明在單軸擠壓力的作用下，石灰?guī)r承受的壓力最大，強(qiáng)度最大，而泥巖承擔(dān)的壓力最小，強(qiáng)度低。巖石物理試驗(yàn)結(jié)果證明了這一點(diǎn)。

巖體的最小應(yīng)變和中間應(yīng)變的方向與應(yīng)力方向一致，為平行于層面的水平方向。可以看出，泥巖在水平方向的應(yīng)變呈現(xiàn)出最大，砂巖次之，石灰?guī)r最小，即在單軸壓縮下，泥巖發(fā)生的水平膨脹最明顯，z方向泥巖也表現(xiàn)出最大的變形，而砂巖的變形依然較石灰?guī)r大。這是由于泥巖的強(qiáng)度最小，泊松比大，在力的作用下，相對(duì)砂巖和石灰?guī)r更容易發(fā)生水平膨脹和縱向擠壓變形。應(yīng)變強(qiáng)度分布圖同樣表現(xiàn)出泥巖應(yīng)變強(qiáng)度高于砂巖及石灰?guī)r，說(shuō)明泥巖最容易發(fā)生形變，易于首先產(chǎn)生裂紋。

2.2 巖層豎直單軸壓縮

巖層豎直單軸壓縮試驗(yàn)的壓力方向與層面平行，其最大主應(yīng)力方向平行于壓力方向，即z軸方向；最小主應(yīng)力和中間主應(yīng)力垂直于最大主應(yīng)力，即x，y所在平面?？梢钥闯瞿鄮r在y方向的應(yīng)變較大，而石灰?guī)r最小；z分量應(yīng)變與最大應(yīng)變方向平行，其應(yīng)變分布規(guī)律也基本一致：泥巖在z方向的應(yīng)變?nèi)匀蛔畲?，砂巖次之，石灰?guī)r最小。應(yīng)變強(qiáng)度分布圖直觀(guān)地表明了各巖體應(yīng)變強(qiáng)度的順序?yàn)槟鄮r＞砂巖＞石灰?guī)r（見(jiàn)圖6）。

圖6 巖層豎直的單軸壓縮的數(shù)值模擬結(jié)果Fig.6 The numerical simulation results of uniaxial compression experiment with vertical layers

應(yīng)力在巖體上的分布有如下特征：石灰?guī)r和砂巖上的最大主應(yīng)力值較大，而泥巖的最大主應(yīng)力值較小，應(yīng)力強(qiáng)度的分布也顯示此特征。這是由于豎直方向施加壓力時(shí)，三者同時(shí)受力，且位移相同，故彈性模量和強(qiáng)度比較大的石灰?guī)r和砂巖集中受力，承受的壓力較大。若增大載荷，石灰?guī)r首先出現(xiàn)裂紋，繼續(xù)施壓，裂紋擴(kuò)大，砂巖、泥巖緊接著出現(xiàn)裂紋，隨后破裂貫穿石灰?guī)r分層，發(fā)生完全破壞。該情況下的應(yīng)力-應(yīng)變曲線(xiàn)如圖2b，與水平巖層單軸試驗(yàn)曲線(xiàn)相比，此時(shí)曲線(xiàn)的第二階段具有較高的斜率，當(dāng)達(dá)到峰值強(qiáng)度后，出現(xiàn)突變式下落，代表完全破壞的時(shí)刻，之后恢復(fù)平緩，是由于砂巖和泥巖很快承受了壓力。

2.3 巖層水平雙軸壓縮

雙向加載試驗(yàn)與單向加載試驗(yàn)最大的不同是，雙向約束時(shí)只有一個(gè)自由面，相當(dāng)于在豎向單軸試驗(yàn)的基礎(chǔ)上增加了一組約束（見(jiàn)圖3c），因此受力時(shí)產(chǎn)生垂直于分界面向外的拉力。

巖層水平雙軸試驗(yàn)的模擬結(jié)果如圖7。從圖7中可見(jiàn)，泥巖在z方向應(yīng)力與應(yīng)變均小于石灰?guī)r和砂巖的應(yīng)力與應(yīng)變，而最小應(yīng)力和應(yīng)變的方向與z軸平行，故泥巖最小應(yīng)力、最小應(yīng)變分布與z向應(yīng)力、應(yīng)變基本一致。應(yīng)力強(qiáng)度仍表現(xiàn)為泥巖較低，而石灰?guī)r較高。此外，泥巖和砂巖的應(yīng)變強(qiáng)度高于石灰?guī)r。

圖7 巖層水平的雙軸壓縮的數(shù)值模擬結(jié)果Fig.7 The numerical simulation results of biaxial compression experiment with horizontal layers

該情況下z方向受力產(chǎn)生拉張，由于x和y方向均受約束，所以位移一致，故應(yīng)力集中于石灰?guī)r上；隨著擠壓增大，石灰?guī)r會(huì)首先破裂，其在z方向的高應(yīng)變量說(shuō)明了石灰?guī)r已有產(chǎn)生裂紋的趨勢(shì)。應(yīng)力-應(yīng)變曲線(xiàn)圖（見(jiàn)圖2c）中，波峰為石灰?guī)r破裂的時(shí)刻，隨后的似階梯狀變化代表石灰?guī)r破裂后，壓力由砂巖和泥巖承受，以致其逐漸開(kāi)裂、破壞。

2.4 巖層豎直雙軸壓縮

巖層豎直的雙軸壓縮試驗(yàn)是在橫向單軸試驗(yàn)的基礎(chǔ)上增加了一組約束，受力時(shí)會(huì)產(chǎn)生平行于界面的壓力。其最小應(yīng)力與應(yīng)變方向平行于z方向，在該方向上泥巖的應(yīng)力值高于石灰?guī)r和砂巖，應(yīng)變值小于砂巖和石灰?guī)r（見(jiàn)圖8）。這是因?yàn)樵诋?dāng)前擠壓力（15 MPa）下，石灰?guī)r和砂巖先后發(fā)生破裂，應(yīng)力釋放，此時(shí)泥巖為應(yīng)力的主要承受者，又由于泥巖具較高的泊松比，發(fā)生了更大的形變。

2.5 三軸試驗(yàn)

三軸試驗(yàn)最大應(yīng)力方向?yàn)閳A柱的軸向，從圖9可以看出泥巖在軸向的應(yīng)變最大，石灰?guī)r的應(yīng)變量最小，砂巖的應(yīng)變則介于兩者之間。至于應(yīng)力的分布，沒(méi)有明顯的石灰?guī)r＞砂巖＞泥巖的現(xiàn)象，這是由于圍壓的存在，增大了巖體的抗壓強(qiáng)度。巖石力學(xué)試驗(yàn)結(jié)果也表明，三軸擠壓的情況下巖體沒(méi)有明顯的破裂，泥巖的膨脹現(xiàn)象也并不突出，巖體的塑性大大增強(qiáng)。應(yīng)力-應(yīng)變曲線(xiàn)（見(jiàn)圖2e）也表現(xiàn)出了這個(gè)特征，三軸壓縮的曲線(xiàn)沒(méi)有明顯的抗壓極限，壓力增大到一定程度保持平緩，沒(méi)有急劇下降的階段，即巖體未完全破壞，塑性特征逐漸呈現(xiàn)出來(lái)。

受力方式對(duì)破裂的影響主要體現(xiàn)在作用于巖體的約束及載荷的方位上。當(dāng)力與巖層面平行時(shí)，約束面則和層面垂直，像石灰?guī)r這樣強(qiáng)度大的“強(qiáng)硬”巖體受到集中應(yīng)力，故首先破裂；而當(dāng)力與巖層面垂直時(shí)，垂直層面的方向沒(méi)有約束，泥巖這樣強(qiáng)度小的“軟弱”巖體首先破裂。

單軸壓縮的自由面最多，變形、破裂的空間最廣，因此巖體的應(yīng)變量和位移量最大，變形也以彈性為主；而三軸壓縮的約束較緊，應(yīng)變量較小，當(dāng)載荷超過(guò)一定限度變形則轉(zhuǎn)化為塑性變形。

圖8 巖層豎直雙軸壓縮的數(shù)值模擬結(jié)果Fig.8 The numerical simulation results of biaxial compression experiment with vertical layers

圖9 三軸試驗(yàn)的數(shù)值模擬結(jié)果Fig.9 The numerical simulation results of triaxial compression experiment

3 研究意義與實(shí)例驗(yàn)證

探討層狀巖體的破裂規(guī)律對(duì)儲(chǔ)層裂縫研究有直接的意義。根據(jù)不同地質(zhì)體裂縫的開(kāi)裂、延伸、擴(kuò)展特征，可以預(yù)測(cè)裂縫富集帶，特別是在致密低滲透巖體里，裂縫是重要的油氣儲(chǔ)集空間。商琳等［11］對(duì)富臺(tái)碳酸鹽潛山各期巖層做了深入的巖石試驗(yàn)及裂縫數(shù)值模擬研究。從該區(qū)燕山期地層的力學(xué)參數(shù)及裂縫參數(shù)特征（見(jiàn)表2）可以看出，彈性模量較大、泊松比較小的鳳山組和八陡組巖層發(fā)育的裂縫開(kāi)度較大，為2.9～4.7 mm，而在相對(duì)較“軟”的冶里—亮甲山組巖層發(fā)育的裂縫開(kāi)度較小，為2.5～2.9 mm，表明較“強(qiáng)硬”的巖層容易破裂，且裂縫較寬。

表2 富臺(tái)潛山燕山期力學(xué)參數(shù)及裂縫參數(shù)特征Table 2 The mechanical parameters and fracture features of Futai Buried Hill in Yanshanian

富臺(tái)地區(qū)在燕山期主要受北東東—南西西向擠壓作用，各組巖層基本水平，類(lèi)似于水平巖層的雙軸壓縮或豎直巖層的單軸壓縮，該受力狀態(tài)下層狀巖體中的“強(qiáng)硬”巖層易于破裂。故巖石試驗(yàn)、數(shù)值模擬結(jié)果與富臺(tái)地區(qū)實(shí)例吻合，為本研究的正確性和合理性提供了支持。

4 結(jié)論

不同受力方式顯著影響層狀巖體的破裂順序，相同受力方式下的破壞順序與巖層排列順序無(wú)明顯關(guān)系。巖層水平和豎直兩種情況下的單軸壓縮試驗(yàn)結(jié)果具有明顯差異，巖層水平時(shí)首先在強(qiáng)度小的“軟弱”巖體出現(xiàn)裂紋，最終貫穿整體；而巖層豎直時(shí)，強(qiáng)度大的“強(qiáng)硬”巖體首先破裂。

不同約束方式影響層狀巖體的破裂順序。約束面與巖層面平行時(shí)，“軟弱”巖體首先變形、破壞，約束面垂直層面時(shí)，“強(qiáng)硬”巖體先破裂。

巖體的彈性和塑性并非其固有屬性，受外界條件的影響。圍壓的存在大大提高了巖體的強(qiáng)度，強(qiáng)圍壓能使彈性向塑性轉(zhuǎn)化。

ANSYS數(shù)值模擬法在本次研究中得到了合理而正確的應(yīng)用，結(jié)果滿(mǎn)足試驗(yàn)的要求并符合客觀(guān)的情況，為本次研究提供了保障。

［1］ 胡明研.壓縮條件下層狀復(fù)合巖石力學(xué)特性研究［D］.焦作：河南理工大學(xué)，2009.

HU Yan-ming.Mechanical characteristic study of stratified composite rock under compression［D］.Henan polytechnic university，2009.

［2］ 黎立云，寧海龍，劉志寶，等.層狀巖體斷裂破壞特殊現(xiàn)象及機(jī)制分析［J］.巖石力學(xué)與工程學(xué)報(bào)，2006，25（2）：3933～3938.

LI Li-yun，NING Hai-long，LIU Zhi-bao，et al.Special phenomenon of fracture and mechanism analysis of layered rock mass［J］.Chinese Journal of Rock Mechanics and Engineering，2006，25（2）：3933～3938.

［3］ 劉立，邱賢德，黃木坤，等.層狀復(fù)合巖石損傷破壞的實(shí)驗(yàn)研究［J］.重慶大學(xué)學(xué)報(bào)：自然科學(xué)版，1999，22（4）：28～32.

LIU Li，QIU Xian-de，HUANG Mu-kun，et al.Experimental study on damage and failure of layered composite rock［J］. Journal of Chongqing University：Natural Science Edition，1999，22（4）：28～32.

［4］ 許寶田，閻長(zhǎng)虹，許宏發(fā).三軸試驗(yàn)?zāi)鄮r應(yīng)力-應(yīng)變特性分析［J］.巖土工程學(xué)報(bào)，2005，26（6）：863～865.

XU Bao-tian，YAN Chang-hong，XU Hong-fa.Triaxial tests on stress-strain of mud stone［J］.Chinese Journal of Geotechnical Engineering，2005，26（6）：863～865.

［5］ 戴俊生，汪必峰，馬占榮.脆性低滲透砂巖破裂準(zhǔn)則研究［J］.新疆石油地質(zhì)，2007，28（4）：393～395.

DAI Jun-sheng，WANG Bi-feng，MA Zhan-rong.Research on Cracking Principles of Brittle Low-Permeability Sands［J］. Xinjiang Petroleum Geology，2007，28（4）：393～395.

［6］ 盧虎勝，張俊峰，李世銀，等.砂泥巖間互地層破裂準(zhǔn)則選取及在裂縫穿透性評(píng)價(jià)中的應(yīng)用［J］.中國(guó)石油大學(xué)學(xué)報(bào)：自然科學(xué)版，2012，36（3）：14～19.

LU Hu-sheng，ZHANG Jun-feng，LI Shi-yin，et al.Best selection among crack criterions to sandstone-mudstone formation and its application to fracture penetration ability evaluation［J］.Journal of China University of Petroleum，2012，36（3）：14～19.

［7］ 雷宇.砂、泥巖壓縮、回彈、破裂實(shí)驗(yàn)與裂縫成因的初步分析［D］.西安：西北大學(xué)，2011.

LEI Yu.The experiment of compression，rebound and fracturing of sandstone and mudstone and preliminary analysis on origin of cracks［D］.Northwest University，2011.

［8］ 楊仕教，曾晟，王和龍.加載速率對(duì)石灰?guī)r力學(xué)效應(yīng)的試驗(yàn)研究［J］.巖土工程學(xué)報(bào)，2005，27（7）：786～788.

YANG Shi-jiao，ZENG Sheng，WANG He-long.Experimental analysis on mechanical effects of loading rates on limestone［J］.Chinese Journal of Geotechnical Engineering，2005，27（7）：786～788.

［9］ 李永盛.加載速率對(duì)紅砂巖力學(xué)效應(yīng)的試驗(yàn)研究［J］.同濟(jì)大學(xué)學(xué)報(bào)：自然科學(xué)版，1995，23（3）：265～269.

LI Yong-sheng.Experimental analysis on the mechanical effects of loading rates on red sandstone［J］.Journal of Tongji University，1995，23（3）：265～269.

［10］ 尹小濤，葛修潤(rùn)，李春光，等.加載速率對(duì)巖石材料力學(xué)行為的影響［J］.巖石力學(xué)與工程學(xué)報(bào)，2010，29（S1）：2610～2615.

YIN Xiao-tao，GE Xiu-run，LI Chun-guang，et al.Influences of loading rates on mechanical behaviors of rock materials［J］.Chinese Journal of Rock Mechanics and Engineering，2010，29（S1）：2610～2615.

［11］ 商琳，戴俊生，賈開(kāi)富，等.碳酸鹽巖潛山不同級(jí)別構(gòu)造裂縫分布規(guī)律數(shù)值模擬——以渤海灣盆地富臺(tái)油田為例［J］.天然氣地球科學(xué)，2013，24（6）：1260～1267.

SHANG Lin，DAI Jun-sheng，JIA Kai-fu，et al.Numerical simulation for the distribution of different levels of tectonicfractures in carbonate buried hills：Taking Futai Oilfield in Bohai Bay Basin as an example［J］.Natural Gas Geoscience，2013，24（6）：1260～1267.

［12］ 徐素國(guó)，梁衛(wèi)國(guó)，莫江，等.軟弱泥巖夾層對(duì)層狀鹽巖體力學(xué)特性影響研究［J］.地下空間與工程學(xué)報(bào)，2009，5（5）：878～883.

XU Su-guo，LIANG Wei-guo，MO Jiang，et al.Influence of weak mudstone intercalated layer on mechanical properties of laminated salt rock［J］.Chinese Journal of Underground Space and Engineering，2009，5（5）：878～883.

［13］ 姜德義，任濤，陳結(jié)，等.含軟弱夾層鹽巖型鹽力學(xué)特性試驗(yàn)研究［J］.巖石力學(xué)與工程學(xué)報(bào)，2012，31（9）：1797～1803.

JIANG De-yi，REN Tao，CHEN Jie，et al.Experiment study of mechanical characteristics of molded salt rock with weak interlayer［J］.Chinese Journal of Rock Mechanics and Engineering，2012，31（9）：1797～1803.

DISCUSSION ON THE FRACTURE OF LAYERED ROCK MASS BASED ON THE FINITE ELEMENT METHOD

XU Ke1，DAI Jun-sheng1，F(xiàn)U Xiao-long1，REN Qi-qiang1，LIU Cong-ning2，ZHAO Yun3

（1.School of Geoscience，China University of Petroleum，Qingdao 266580，Shangdong，China；
2.Jidong Oil Exploration and Development Institute，CNPC，Tangshan 063004，Hebei，China；
3.School of education，Guangzhou University，Guangzhou 510006，China）

The failure mode and laws of layered rock mass are unique because of its heterogeneity and anisotropy.By conducting uniaxial，biaxial and three axis compression experiments on layered rock mass and analyzing their stress strain curves，different characteristics have been shown.And then with the software ANSYS，finite element numerical simulation have been done.Comparing stress strain curves and maps of simulation as well as combining the theory of fracture，failure modes and rules of layered rock mass in different stress conditions have been summarized.Finally Futai area was taken as an example to validate the result.This research suggests that the orientations of loads and constraints as well as the rock level are the key factors to rock failure.Rock with strong intensity like limestone would be cracked in priority if stress direction is paralleled to rock level，for stress concentration occurs to the strong rock.While rock with weak intensity like mudstone would be cracked first if stress direction is perpendicular to the rock level.Rock compression experiment，numerical simulation and the example of Futai area have verified this law successfully.

layered rock mass；rock compression experiment；numerical simulation；failure laws

TD315

A

1006-6616（2015）03-0330-11

2014-12-26

國(guó)家科技重大專(zhuān)項(xiàng)“復(fù)雜裂縫性碳酸鹽巖油藏開(kāi)發(fā)關(guān)鍵技術(shù)”（2011ZX05014-004）；中國(guó)石油大學(xué)（華東）研究生創(chuàng)新工程（YCX2015009）

徐珂（1991-），男，碩士研究生，主要從事構(gòu)造地質(zhì)學(xué)方面的研究。E-mail：232437658＠qq.com