層狀巖體破裂規律探析范文

本站小編為你精心準備了層狀巖體破裂規律探析參考范文，愿這些范文能點燃您思維的火花，激發您的寫作靈感。歡迎深入閱讀并收藏。

《地質力學學報》2016年第2期

摘要:

層狀巖體的非均質性及各向異性導致其破裂方式及規律與均質巖體有顯著不同。對層狀巖體分別進行不同方式的單軸、雙軸、三軸試驗，分析應力-應變曲線特征;再利用ANSYS有限元軟件進行數值模擬，觀察應力、應變在巖體上的分布，通過曲線和圖件的對比分析，并結合巖石破裂理論，總結不同應力狀態下層狀巖體的破裂方式、順序及規律;最后以富臺地區為例，對分析結果進行驗證。研究結果表明，不同受力方式對層狀巖體破裂的影響體現在施加的載荷及約束與層面的方位。當應力方向與巖層面平行時，強度大的石灰巖巖體發生集中應力，首先破裂;而應力與巖層面垂直時，強度小的泥巖巖體首先破裂。巖石試驗、數值模擬結果以及實例均成功驗證了這個規律。

關鍵詞:

層狀巖體;巖石試驗;數值模擬;破裂規律

0引言

層狀巖體普遍分布于自然界中。與均質塊狀巖體相比，層狀巖體的變形和強度特性具有明顯的各向異性，因此其穩定性和破壞條件也有特殊表現，特別是裂縫的產生和延伸變得比較復雜，現有的破裂準則也難以精準描述它們的破裂規律。目前研究過層狀巖體力學性質和變形方式的學者不在少數，胡明研等［1］對層狀復合巖石做了壓縮試驗，分析了不同條件下裂縫在巖層上的發育和展布;黎立云等［2］通過巴西圓盤劈裂試驗及單軸壓縮試驗討論了層狀巖體破裂的特殊現象及成因機制;劉立［3］對復合巖體的特殊性質進行了深入研究，在損傷演化和斷裂破壞規律方面有了實質進展;許寶田等［4］對泥巖進行了三軸壓縮實驗，分析了泥巖應力-應變特性;戴俊生等［5］提出了適合脆性低滲透砂巖的破裂準則，盧虎勝等［6］則在戴俊生的基礎上建立了用于砂泥巖間互地層的破裂準則;雷宇［7］通過砂泥巖壓縮、回彈、破裂實驗對裂縫成因進行了初步分析;楊仕教等［8］、李永盛［9］和尹小濤等［10］分別研究了不同加載速率對巖石力學效應的影響。利用ANSYS軟件進行有限元數值模擬，可以對不連續的非均勻介質組成的巖體進行應力特征研究。將整個巖體分為數目有限的單元，通過分析計算每個單元的應力場進而綜合所有單元來研究整體特征。單元的位移、應變及應力都可以經過計算而得出，因此巖體在受力狀態下的應力、應變及位移分布可以直觀地表現出來。本文以巖石試驗為基礎，并用ANSYS軟件建立與實體1:1的模型，約束與載荷按照真實情景加載，分析層狀巖體在不同應力狀態下的應力、應變及強度分布，據此推斷和解釋巖體的破裂規律。研究成果有助于判斷裂縫的產生、發育與分布，對儲層裂縫帶的預測具有指導意義。

1巖石力學試驗

巖石受到不同構造應力作用時，將發生不同形式的變形。當其所承受的應力超過極限強度，就會發生破裂并產生裂縫。由于不同巖石本身性質及所受應力條件不同，破裂的方式也不同。對于層狀巖石而言，各分層不同的力學性質必然導致破裂規律的特殊性，而不同的應力狀態也會產生不同的破裂方式。例如，裂紋在均質花崗巖上的擴展比較平直(見圖1a)，而在層狀花崗巖上，裂紋越過層面時擴展路徑發生變化，呈鋸齒狀(見圖1b)［2］。本文根據相似原理［1］，利用與泥巖、砂巖和石灰巖3種常見且性質差異明顯的材料，制成復合成層狀巖體，分別為邊長100mm的正方體及直徑50mm、高100mm的圓柱體。通過對單一巖體進行三軸試驗確定其彈性模量和泊松比。計算公式如下:E=σ2－σ1εz2－εz1(1)μ=εx1－εx2εz1－εz2(2)式中:E為巖石彈性模量，GPa;μ為巖石泊松比;σ1、σ2分別為應力-應變曲線直線段開始和結束時的應力，MPa;εz1、εz2分別為應力-應變曲線直線段開始和結束時的軸向應變［11］，mm;εx1、εx2分別為應力-應變曲線直線段開始和結束時的側向應變，mm。試驗得到的參數如表1。試驗用的每件樣品應完好無損，分別進行單軸、雙軸、三軸壓縮試驗，其中單軸和雙軸試驗再分為巖層水平和豎直兩種情況。試驗得到的應力-應變曲線如圖2。除三軸壓縮外，其余4曲線雖大致趨勢相似，但形態不盡相同，意味著不同受力方式對巖石破裂有重要影響。

2數值模擬與分析

模擬依據巖石力學試驗的5種方式進行，即巖層水平的單軸壓縮，巖層豎直的單軸壓縮，巖層水平的雙軸壓縮，巖層豎直的雙軸壓縮，三軸壓縮試驗(見圖3)。

2.1巖層水平的單軸壓縮

首先研究分層排列順序對破裂的影響。進行3次試驗，各層自上而下的排列順序分別為:①砂巖、泥巖、石灰巖;②泥巖、砂巖、石灰巖;③砂巖、石灰巖、泥巖。以試驗①為例，隨著壓力增大，強度最弱的泥巖先出現裂紋，但不能很快貫通界面;繼續增加荷載，當裂縫的擴張力大于層間的束縛力時，砂巖和石灰巖也開始出現裂紋，且很快與泥巖的裂紋形成了貫通縫(見圖4)。試驗②、③的結果與試驗①類似，總是強度最弱的泥巖首先破裂。因此，受力條件相同的情況下，破裂與巖層排列順序無明顯關系。水平巖層單軸壓縮的應力-應變曲線(見圖2a)大體分4個階段:孔隙壓密階段，表現為曲線“上凹”;彈性變形階段，表現為近直線;破裂發展階段，其上界為峰值強度;破裂后階段，曲線迅速下降。使用ANSYS15.0軟件進行數值模擬。模型按實際試件1∶1比例建立，層間為Glue方式膠結。采用Solid45單元進行網格劃分，每個單元為8節點，z方向施加擠壓力。為了看出應力-應變分布的顯著變化，擠壓力設定為15MPa。模擬計算結果見圖5。3個主應力中，最小主應力和中間主應力均平行于層面，最大主應力豎直向下。應力分布特征:泥巖和石灰巖的最小主應力和中間主應力為擠壓應力，泥巖應力值低于石灰巖，而砂巖的最小主應力和中間主應力呈現部分拉張應力。這是由于泥巖和砂巖交界處的水平變形不一致，而砂巖的彈性模量大于泥巖，且泊松比小于泥巖，導致泥巖的水平變形大于砂巖。在巖體受壓至達到峰值強度的過程中，由于變形特性的差異及變形連續靜力平衡條件，交界面附近的砂巖和泥巖在水平方向上必然發生應力突變，就會在界面上產生方向相反的應力［12］。正是這個原因，在交界處附近的石灰巖層上也出現局部拉張力。應力強度分布圖說明在單軸擠壓力的作用下，石灰巖承受的壓力最大，強度最大，而泥巖承擔的壓力最小，強度低。巖石物理試驗結果證明了這一點。巖體的最小應變和中間應變的方向與應力方向一致，為平行于層面的水平方向。可以看出，泥巖在水平方向的應變呈現出最大，砂巖次之，石灰巖最小，即在單軸壓縮下，泥巖發生的水平膨脹最明顯，z方向泥巖也表現出最大的變形，而砂巖的變形依然較石灰巖大。這是由于泥巖的強度最小，泊松比大，在力的作用下，相對砂巖和石灰巖更容易發生水平膨脹和縱向擠壓變形。應變強度分布圖同樣表現出泥巖應變強度高于砂巖及石灰巖，說明泥巖最容易發生形變，易于首先產生裂紋。

2.2巖層豎直單軸壓縮

巖層豎直單軸壓縮試驗的壓力方向與層面平行，其最大主應力方向平行于壓力方向，即z軸方向;最小主應力和中間主應力垂直于最大主應力，即x，y所在平面。可以看出泥巖在y方向的應變較大，而石灰巖最小;z分量應變與最大應變方向平行，其應變分布規律也基本一致:泥巖在z方向的應變仍然最大，砂巖次之，石灰巖最小。應變強度分布圖直觀地表明了各巖體應變強度的順序為泥巖＞砂巖＞石灰巖(見圖6)。應力在巖體上的分布有如下特征:石灰巖和砂巖上的最大主應力值較大，而泥巖的最大主應力值較小，應力強度的分布也顯示此特征。這是由于豎直方向施加壓力時，三者同時受力，且位移相同，故彈性模量和強度比較大的石灰巖和砂巖集中受力，承受的壓力較大。若增大載荷，石灰巖首先出現裂紋，繼續施壓，裂紋擴大，砂巖、泥巖緊接著出現裂紋，隨后破裂貫穿石灰巖分層，發生完全破壞。該情況下的應力-應變曲線如圖2b，與水平巖層單軸試驗曲線相比，此時曲線的第二階段具有較高的斜率，當達到峰值強度后，出現突變式下落，代表完全破壞的時刻，之后恢復平緩，是由于砂巖和泥巖很快承受了壓力。

2.3巖層水平雙軸壓縮

雙向加載試驗與單向加載試驗最大的不同是，雙向約束時只有一個自由面，相當于在豎向單軸試驗的基礎上增加了一組約束(見圖3c)，因此受力時產生垂直于分界面向外的拉力。巖層水平雙軸試驗的模擬結果如圖7。從圖7中可見，泥巖在z方向應力與應變均小于石灰巖和砂巖的應力與應變，而最小應力和應變的方向與z軸平行，故泥巖最小應力、最小應變分布與z向應力、應變基本一致。應力強度仍表現為泥巖較低，而石灰巖較高。此外，泥巖和砂巖的應變強度高于石灰巖。該情況下z方向受力產生拉張，由于x和y方向均受約束，所以位移一致，故應力集中于石灰巖上;隨著擠壓增大，石灰巖會首先破裂，其在z方向的高應變量說明了石灰巖已有產生裂紋的趨勢。應力-應變曲線圖(見圖2c)中，波峰為石灰巖破裂的時刻，隨后的似階梯狀變化代表石灰巖破裂后，壓力由砂巖和泥巖承受，以致其逐漸開裂、破壞。

2.4巖層豎直雙軸壓縮

巖層豎直的雙軸壓縮試驗是在橫向單軸試驗的基礎上增加了一組約束，受力時會產生平行于界面的壓力。其最小應力與應變方向平行于z方向，在該方向上泥巖的應力值高于石灰巖和砂巖，應變值小于砂巖和石灰巖(見圖8)。這是因為在當前擠壓力(15MPa)下，石灰巖和砂巖先后發生破裂，應力釋放，此時泥巖為應力的主要承受者，又由于泥巖具較高的泊松比，發生了更大的形變。

2.5三軸試驗

三軸試驗最大應力方向為圓柱的軸向，從圖9可以看出泥巖在軸向的應變最大，石灰巖的應變量最小，砂巖的應變則介于兩者之間。至于應力的分布，沒有明顯的石灰巖＞砂巖＞泥巖的現象，這是由于圍壓的存在，增大了巖體的抗壓強度。巖石力學試驗結果也表明，三軸擠壓的情況下巖體沒有明顯的破裂，泥巖的膨脹現象也并不突出，巖體的塑性大大增強。應力-應變曲線(見圖2e)也表現出了這個特征，三軸壓縮的曲線沒有明顯的抗壓極限，壓力增大到一定程度保持平緩，沒有急劇下降的階段，即巖體未完全破壞，塑性特征逐漸呈現出來。受力方式對破裂的影響主要體現在作用于巖體的約束及載荷的方位上。當力與巖層面平行時，約束面則和層面垂直，像石灰巖這樣強度大的“強硬”巖體受到集中應力，故首先破裂;而當力與巖層面垂直時，垂直層面的方向沒有約束，泥巖這樣強度小的“軟弱”巖體首先破裂。單軸壓縮的自由面最多，變形、破裂的空間最廣，因此巖體的應變量和位移量最大，變形也以彈性為主;而三軸壓縮的約束較緊，應變量較小，當載荷超過一定限度變形則轉化為塑性變形。

3研究意義與實例驗證探討

層狀巖體的破裂規律對儲層裂縫研究有直接的意義。根據不同地質體裂縫的開裂、延伸、擴展特征，可以預測裂縫富集帶，特別是在致密低滲透巖體里，裂縫是重要的油氣儲集空間。商琳等［11］對富臺碳酸鹽潛山各期巖層做了深入的巖石試驗及裂縫數值模擬研究。從該區燕山期地層的力學參數及裂縫參數特征(見表2)可以看出，彈性模量較大、泊松比較小的鳳山組和八陡組巖層發育的裂縫開度較大，為2.9～4.7mm，而在相對較“軟”的冶里—亮甲山組巖層發育的裂縫開度較小，為2.5～2.9mm，表明較“強硬”的巖層容易破裂，且裂縫較寬。富臺地區在燕山期主要受北東東—南西西向擠壓作用，各組巖層基本水平，類似于水平巖層的雙軸壓縮或豎直巖層的單軸壓縮，該受力狀態下層狀巖體中的“強硬”巖層易于破裂。故巖石試驗、數值模擬結果與富臺地區實例吻合，為本研究的正確性和合理性提供了支持。

4結論

不同受力方式顯著影響層狀巖體的破裂順序，相同受力方式下的破壞順序與巖層排列順序無明顯關系。巖層水平和豎直兩種情況下的單軸壓縮試驗結果具有明顯差異，巖層水平時首先在強度小的“軟弱”巖體出現裂紋，最終貫穿整體;而巖層豎直時，強度大的“強硬”巖體首先破裂。不同約束方式影響層狀巖體的破裂順序。約束面與巖層面平行時，“軟弱”巖體首先變形、破壞，約束面垂直層面時，“強硬”巖體先破裂。巖體的彈性和塑性并非其固有屬性，受外界條件的影響。圍壓的存在大大提高了巖體的強度，強圍壓能使彈性向塑性轉化。ANSYS數值模擬法在本次研究中得到了合理而正確的應用，結果滿足試驗的要求并符合客觀的情況，為本次研究提供了保障。

參考文獻：

［1］胡明研．壓縮條件下層狀復合巖石力學特性研究［D］．焦作:河南理工大學，2009.

［2］黎立云，寧海龍，劉志寶，等．層狀巖體斷裂破壞特殊現象及機制分析［J］．巖石力學與工程學報，2006，25(2):3933～3938.

［3］劉立，邱賢德，黃木坤，等．層狀復合巖石損傷破壞的實驗研究［J］．重慶大學學報:自然科學版，1999，22(4):28～32.

［4］許寶田，閻長虹，許宏發．三軸試驗泥巖應力-應變特性分析［J］．巖土工程學報，2005，26(6):863～865.

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

［6］盧虎勝，張俊峰，李世銀，等．砂泥巖間互地層破裂準則選取及在裂縫穿透性評價中的應用［J］．中國石油大學學報:自然科學版，2012，36(3):14～19.

［7］雷宇．砂、泥巖壓縮、回彈、破裂實驗與裂縫成因的初步分析［D］．西安:西北大學，2011.

［8］楊仕教，曾晟，王和龍．加載速率對石灰巖力學效應的試驗研究［J］．巖土工程學報，2005，27(7):786～788.

［9］李永盛．加載速率對紅砂巖力學效應的試驗研究［J］．同濟大學學報:自然科學版，1995，23(3):265～269.

［10］尹小濤，葛修潤，李春光，等．加載速率對巖石材料力學行為的影響［J］．巖石力學與工程學報，2010，29(S1):2610～2615.

［11］商琳，戴俊生，賈開富，等．碳酸鹽巖潛山不同級別構造裂縫分布規律數值模擬———以渤海灣盆地富臺油田為例［J］．天然氣地球科學，2013，24(6)

［12］徐素國，梁衛國，莫江，等．軟弱泥巖夾層對層狀鹽巖體力學特性影響研究［J］．地下空間與工程學報，2009，5(5):878～883.

［13］姜德義，任濤，陳結，等．含軟弱夾層鹽巖型鹽力學特性試驗研究［J］．巖石力學與工程學報，2012，31(9):1797～1803.

作者:徐珂 戴俊生 付曉龍 任啟強劉叢寧 趙惲 單位:中國石油大學地球科學與技術學院 中國石油冀東油田勘探開發研究院 廣州大學教育學院