張楊洋

(西北水利水電工程有限責任公司，西安 710065)

?

巖體工程地應力測定方法綜述及展望

張楊洋

(西北水利水電工程有限責任公司，西安 710065)

原始地應力是巖石地下工程及部分露天開挖工程進行穩(wěn)定性評價和設計的基本參數(shù)，同時也是該類工程發(fā)生變形、失穩(wěn)和破壞的根本因素之一。隨著國防、水利、交通、石油、采礦、核能源等領域工程深度及規(guī)模的不斷加大，地應力測試對指導工程科學、安全、經濟開展的作用越發(fā)顯著，而目前中國地應力測量技術水平和應用范圍還難以滿足眾多工程需要，這就對加強地應力測量理論、技術和應用的研究提出了更高要求。文章對地應力成因和分類、影響因素及目前主要的測量方法進行了總結，并對未來發(fā)展做出了展望。

巖體工程；地應力；測定方法；展望

地應力是引發(fā)各類露天及地下巖體工程變形和破壞的根本力源，是圍巖穩(wěn)定性判斷的基礎，同時也是進行巖體工程科學設計和決策的必要前提[1]。隨著，國防、交通、水電、采礦、石油及煤層氣、核能源開發(fā)利用等工程發(fā)展需要的加深，給地應力測量工作提出了新的挑戰(zhàn)，加強地應力測量理論、方法、技術的研究，對促進地球科學及巖體工程相關領域的發(fā)展具有重要意義[2-7]。

1 地應力成因及分類

1.1 地應力成因

地球的各種動力運動是地應力形成的主因，具體包含：板塊間擠壓、地幔熱對流、地球內應力、引力、地球自轉、巖漿體侵入和地殼非均勻擴容等。另外，溫差、水壓梯度以及其他物理、化學變化也能引發(fā)應力場的改變。不同環(huán)境下形成因素不同，但基本以自重應力和構造應力為主導。

1.2 地應力分類

依照形成因素不同地應力可分為：構造、自重應力、殘余應力及變異應力。需要注意的是，殘余應力不具有方向性，但殘余應力引起的高水平應力具有方向性，且二者常相差較大。

2 地應力的影響因素

關于地應力的影響因素，世界各國都開展了一定研究，但研究結果存在一定差別。如何正確認識每種因素對地應力產生和分布的影響，是合理應用地應力場資料的基礎。Brudy[8]和Brown[9]對影響地應力的因素進行了如下歸納，包括：地形、殘余應力、侵蝕、巖漿侵入、構造、破裂和結構面以及溫度場的變化因素等；Amadei[10]和Stephansson[11]把影響地應力的因素歸結為巖性(地層)、各向異性、構造和非均稱性、地形及構造運動、侵蝕作用、超固結作用、以及邊界和時間、地球球面彎曲等；陳宗基[12]認為地應力的產生是巖體自重、板塊運動、剝蝕、地形、及封閉應力引起；王思敬[13]從巖石的地質本質性出發(fā)對地應力受地質作用的影響機理進行了論述?？傊?，地應力的影響因素多樣，各因素的影響效應各不相同，但應當強調的是，總體而言斷裂構造是地應力場的主控因素。

3 地應力測試方法

3.1 地應力測量的基本方法及分類

根據(jù)彈性力學，1點的應力可由6個應力分量來表達，即(σx、σy、σz、τxy、τyz、τzx)，如圖1所示。當該6個值確定后，即可得出過該點的主應力狀態(tài)。

圖1 任一點的應力狀態(tài)圖

早期地應力測定多采用硐室開挖扁千斤頂系數(shù)修正法，該方法存在修正系數(shù)選取經驗性較強并且難以準確還原擾動前原巖應力的問題。目前多采用鉆孔測量法，即在壁面上開鉆小孔，打至原巖應力區(qū)，以保證數(shù)據(jù)理論上的可靠性。測量方法的分類尚無統(tǒng)一定論，根據(jù)測量手段和原理不同可進行分類，見表1。

表1 測試方法分類表

國外一般根據(jù)測量基本原理的不同而分為直接法和間接法，見表2。

直接法的特點是所測量值直接與應力值相關，計算公式較為簡單，無需知道巖體物理力學參數(shù)及本構關系。間接法通過測取與應力有關的參數(shù)如變形、密度、孔隙率、滲透性、導電率、波速等，再通過參數(shù)-應力函數(shù)關系計算原巖應力，此方法須事先知道巖體的相關參數(shù)，從而確定函數(shù)關系。

表2 直接、間接法分類表

1987年國際巖石力學學會試驗專業(yè)委員會在《巖石應力測定的建議方法》中推薦了5種比較成熟并最具代表性的方法，即：扁千斤頂法、水壓致裂法、孔徑變形法、孔壁應變法和空心包體應變法，后2種方法可實現(xiàn)單一鉆孔的三維應力測量。表3為根據(jù)操作特點進行的分類。

表3 地應力測試方法分類表

關于地應力測量的理論和方法國內相關學者如陳宗基、蔡美峰、蘇愷之[14]、葛修潤、王連捷[15]、楊永杰[16]、蘇生瑞[17]等都進行了相關研究。

3.2 直接測量法

3.2.1 扁千斤頂法

扁千斤頂法的基本原理是應力補償，通過扁千斤頂對開槽加壓，使得測柱間距恢復初始值, 此時壓力值等于原先巖體的豎向壓力值，測量原理見圖2。

圖2 扁千斤頂應力測量圖

該方法的優(yōu)點在于操作簡便，無需測算巖體本構關系及物理力學參數(shù)。缺點是單槽只能進行一維應力測量，且要求該方向為主應力方向，雖然可采用多槽測量推算原巖應力，但此方法假設前提較多計算結果往往不準。另外，該方法是基于線彈性巖體的基本假設，而巖體加、卸載具有路徑相關性，加之開槽本身形成了次生應力場，故所測量值可能存在較大偏差。類似該方法，還有剛性包體應力計、壓磁式、液壓式、剛弦式及光彈式應力計等。

3.2.2 水壓致裂法

水壓致裂法測量的關鍵是對測點封隔段注高壓水，使該段巖體產生初始裂隙并擴展到3倍孔徑深度，加壓過程記錄水壓、流量的時間曲線，由此根據(jù)彈性力學理論解算出垂直于鉆孔平面的二維應力場量值，后利用印模器或井下攝影等獲得主應力方向，測量設備見圖3。

圖3 水壓致裂法測量圖

該方法的突出優(yōu)點在于可適用于深部應力測量(已超5 000 m)，并且無需獲得巖體參數(shù)，操作較為簡單；缺點是需假定一個主應力方向與鉆孔軸線平行，即只能測得垂直于鉆孔平面的地應力，當假定主應力方向與實際偏差超15°時，結果偏差較大。當封隔段存在原生節(jié)理時，因裂紋尖端的應力場集中，初始破裂很可能在該節(jié)理處產生，而非理論上的孔壁切應力最小點，因此該方法主要用于脆性完整巖石。

3.2.3 聲發(fā)射法

聲發(fā)射法的基本原理是基于巖石的凱澤效應。當材料受外載時，所儲應變能會釋放并產生彈性波和聲響，即聲發(fā)射。凱澤(1950年)發(fā)現(xiàn)，當材料受力小于先期最大應力時，不會發(fā)生聲發(fā)射；當所受應力接近或超過先期最大應力時聲發(fā)射頻率會明顯增加，此現(xiàn)象稱凱澤效應，測試中將聲發(fā)射頻率突增點稱凱澤點，對應的應力值為歷史最大應力值。試驗中通過對不同方向的巖石試件進行凱澤點測試，進而尋找該點的三維原巖應力。聲發(fā)射方法的優(yōu)點是成本低、操作簡便、效率高、限制條件少。缺點是只適用于較堅硬脆性巖石，測試工作量大，受水影響數(shù)據(jù)離散度較大。

3.3 間接測量法

3.3.1 應力解除法

應力解除法是為減少人為測試擾動的常用實測方法，主要包括表面應力解除法和鉆孔應力解除法。

表面應力解除法一般用于巷道或硐室表面的應力測試，操作時，將傳感器布設在巖壁表面，當解除該點后，巖塊發(fā)生彈性恢復，故可以通過該應變恢復值測得原始地應力。

鉆孔應力解除法屬于完全應力解除法，是目前發(fā)展最成熟的方法之一，示意圖如圖4所示。其將空心包體應變計放入鉆孔并粘合在內壁上，再用套鉆將應變計解除，記錄取芯后巖芯的相關應變值，通過彈性力學公式計算原巖應力。依據(jù)測量原理和部位的不同，套芯應力解除法又可細分為鉆孔孔壁、孔底和孔徑3種應變測量法。

圖4 鉆孔應力解除法示意圖

該法的優(yōu)點是理論成熟、誤差小，單孔即可實現(xiàn)三維應力測量，成本相對較低。缺點是套鉆取芯時難度較大，易發(fā)生斷芯，且需要測定該處巖體的變形參數(shù)，操作相對繁瑣。

3.3.2 局部應力解除法

局部應力解除法只能實現(xiàn)測點的局部應力解除，常用方法如下：

(1) 徑向切槽法。用特定的鉆頭在鉆孔壁上沿半徑方向開槽，并記錄開槽后切向應變值。垂直于鉆孔軸向并過該切槽的平面上再開鑿2個同樣尺寸的切槽，即可測得該面上的應力狀態(tài)。

(2) 全息干涉測量法。在鉆孔內垂直于孔壁的方向打一小鉆孔，小鉆孔因局部應力解除而產生孔壁位移。當孔壁發(fā)生微小位移就會引起全息干涉儀中相干光的全息干涉，據(jù)此可測出孔壁變形，從而導出垂直于孔壁平面的二維應力。

(3) 平行鉆孔法。鉆孔打好后，安裝孔徑應變計，再在鉆孔周圍打1個或數(shù)個鉆孔，并記錄第1個鉆孔的孔徑變化值。由含有1個或幾個圓孔的無限大板的彈性理論解或數(shù)值計算，即可算出垂直于鉆孔平面的應力狀態(tài)。

(4) 中心鉆孔法。先在巖體壁面上打直徑為200～250 mm鉆孔，并在圓周上等間距地布設3對測量柱，之后在孔底中心處打直徑為150 mm的鉆孔，并測量由于打小鉆孔所引起測量柱的間距變化。由此即可推導出巖石表面的二維應力狀態(tài)。

3.4 地球物探法

(1) 超聲波測量法：研究發(fā)現(xiàn)巖體中彈性波特別是縱波的傳播，波速和衰減具有應力相關性。因此可以通過測定巖體彈性縱波速的大小，再由已標定的應力-波速關系推導出巖體的當前應力值。值得注意的是，該方法測得的是平面內的主應力狀態(tài)。

(2) 超聲波譜法：向巖體發(fā)射超聲剪切波時，受巖體應力的影響會發(fā)生雙折射現(xiàn)象，研究表明，雙折射率是應力的函數(shù)。故可以通過事先標定的雙折射率與應力關系獲得原始應力值大小。

實際操作中，由于巖體的非均質不連續(xù)地質特點，測量難度很大。地球物理方法測定地應力的技術尚未成熟，測值的可靠度有待考證，故此類方法只能作為一種參考的技術手段。

3.5 震源機制分析法

震源機制分析法對分析深部地應力比較有效。當研究區(qū)域相對震源較大時，那么該震源可視為點震源。通過寫出一組地震矩陣張量或震源機制解，來推算該研究區(qū)域內的平均構造應力、主應力方向及相對比值。

彈性模型是目前常用震源模型，由此計算出的因地震引起的位移場、應變場和應力場，實質上是相對于震前的相對值，所以理論上只能推出震源地震引起的應力變化、區(qū)域構造應力的方向及相對大小量。另一方面，地震波從震源發(fā)出后，傳播過程中可攜帶路徑中介質的應力信息，所以，有利于用地震波來研究地殼的空間受力。

3.6 深部地應力測量方法

當前，世界范圍內對深部地應力的測量主要采用水壓致裂法以及對其進行的改進方法。蔡美峰[18]利用煤礦的地質勘查測孔，進行了水壓致裂的測量工作，測量深度達1 105 m?？导t普[19]等在針對傳統(tǒng)水壓致裂法的改進上做了相關研究，使其能更方便地適用于較深部環(huán)境的地應力測量。盧云虎[20]等，利用粘滯剩磁分量對深部巖芯進行地理方位標定，為Kaiser地應力測試獲取方位角信息提供了解決方法。

葛修潤[21]等研究出一種三維地應力測量方法(BWSRM)——鉆孔局部壁面應力解除法及配套的測井機器人，在錦屏水電站獲得成功應用，取代了沿鉆孔軸向套取巖芯的傳統(tǒng)方法，降低了取芯難度，也摒棄了水壓致裂法必須假定一個主應力方向和鉆孔軸線一致的前提. 王連捷[22]等針對騰沖科學鉆探孔，進行了非彈性應變恢復法(ASR)的地應力測量，實測結果和震源機制解一致，同時具有較高的經濟實用性，該方法為深部三維應力測量提供了新思路。

4 展 望

(1) 深部地層地應力的測量問題

隨著中國經濟的快速增長，對能源的要求逐步增加，煤炭的開采深度以每年15 m左右的速度向下延伸，核能的逐漸使用涉及到核廢料的儲存問題，采油、煤層氣的開發(fā)，這都涉及到深部地應力的確定問題。所以深部地層地應力的測量顯得尤為重要，地應力確定不了，將嚴重制約深部工程技術的發(fā)展。

(2) 測試方法的推新和優(yōu)化

目前地應力測量的應用還較低，主要用于大型工程和部分采礦工程，這主要是測試費用較高及精度不夠導致。因此，如何對已有設備進行改進或者研制出更有普遍適用意義的新設備，將大大提高地下工程建設的科學性和經濟效益。

(3) 復雜地貌地應力場的分析

地殼作為地質體，具有地質本質性，復雜地貌和斷裂構造等處地應力異常復雜。中國還有相當儲量的資源由于儲藏條件復雜，限于工程技術條件，難以開發(fā)。因此，如何根據(jù)測試數(shù)據(jù)提出復雜地貌下的地應力場的分析方法或是設計出針對局部復雜地質結構的地應力測試設備，這對資源合理開發(fā)和指導工程設計意義重大。

[1] 蔡美峰.地應力測量原理和技術[M].北京:科學出版社,2000.

[2] 張文忠.當金山隧道區(qū)域地應力特征分析與應用[J].鐵道工程學報,2014,31(12):18-22.

[3] 袁風波,劉建,李蒲健,等.拉西瓦工程河谷區(qū)高地應力場反演與形成機理[J].巖土力學,2007,28(04):836-842.

[4] 高魁,劉澤功,劉健.地應力在石門揭構造軟煤誘發(fā)煤與瓦斯突出中的作用[J].巖石力學與工程學報,2015, 34(02):305-312.

[5] 宋連騰,劉忠華,李潮流,等. 基于橫向各向同性模型的致密砂巖地應力測井評價方法[J].石油學報,2015(06):707-714.

[6] 李志強,鮮學福,黃滾.地應力地溫場中煤層氣富集區(qū)高精度定量預測的力學方法[J].煤炭學報, 2012, 37(S2):395-400.

[7] 趙星光,王駒,馬利科,等.高放廢物地質處置庫北山預選區(qū)新場巖體地應力場分布規(guī)律[J].巖石力學與工程學報, 2014,33(z2):3750-3759.

[8] BRUDY M,ZOBACK M D,FUCHS K,et al. Estimation of the complete stress tensor to 8 km depth in the KTB scientific drill holes:implications for crustal strength[J].Journal of Geophysical Research,1997,102(B8):453-475.

[9] BROWN E T, HOEK E.Echnical note trends inrelationships between measured in-situ stress and depth[J]. Int. J. Rock Mech. Min. Sci. & Geomech. Abstr.1978(15): 211-215.

[10] AMADEI B, STEPHANSSON O. Rock stress and its measurement[M]. London:Chapman and Hall，1997:95-476.

[11] Stephansson.Stress estimation in rock:a brief history and review[J]. International Journal of Rock Mechanics and Mining Sciences,2003,40(7/8):957-973.

[12] 陳宗基.地球構造動力學[M].北京:北京大學出版社,1982.

[13] 王思敬.論巖石的地質本質性及其巖石力學演繹[J].巖石力學與工程學報.2009,28(03):433-450.

[14] 蘇愷之.地應力測量方法[M].北京:地震出版社,1985.

[15] 王連捷,潘立宙.地應力測量及其在工程中的應用[M].北京:地震出版社,1991.

[16] 楊永杰.地應力測量及綜放沿空掘巷支護技術[M].北京:煤炭工業(yè)出版社,2010.

[17] 蘇生瑞,黃潤秋,王士天. 斷裂構造對地應力場的影響及其工程應用[M].北京:科學出版社,2002.

[18] 蔡美峰,陳長臻,彭華,等.萬福煤礦深部水壓致裂地應力測量[J].巖石力學與工程學報, 2006,25(05):1069-1074.

[19] 康紅普,林建,張曉.深部礦井地應力測量方法研究與應用[J].巖石力學與工程學報,2007,26(05):929-933.

[20] 盧運虎,陳勉,金衍,等.深層地應力地理方位確定的新方法[J].巖石力學與工程學報,2011，30(02):233-237.

[21] 葛修潤,候明勛.三維地應力 BWSRM 測量新方法及其測井機器人在重大工程中的應用[J].巖石力學與工程學報. 2011,30(11):2161-2180.

[22] 王連捷,崔軍文,孫東生,等. 騰沖科學鉆探孔非彈性應變恢復法三維地應力測量[J].地球學報, 2016(01):111-115.

Description and Prospect of Determination Method for Crustal Stress of Rockmass Engineering

ZHANG Yangyang

(Northwest Water Resources and Hydropower Engineering Co., Ltd., Xi'an 710065,China)

The original crustal stress is the basic parameter for stability assessment and design of the underground engineering in rockmass and partial excavation engineering in open air. Meanwhile, the parameter is also one of essential factors of deformation, instability and failure of the engineering. With the increasing of the engineering degree and scale in the fields such as the national defense, water resources, transportation, petroleum, mining and nuclear energy, etc, function of determination of crustal stress in guiding engineering to develop scientifically, safety and economically is quite obvious. But in China, the technical level and application scope of the crustal stress determination currently cannot satisfy requirements of projects in a large quantity. Therefore, study of theory, technology and application of the crustal determination is required to reinforce. In the paper, causes, classification and impact factors of the crustal stress, and main determination methods at present are summarized. Meanwhile, prospects on the future development of the crustal stress determination are performed. Key words: rockmass engineering; crustal stress; determination method; prospect

1006—2610(2016)05—0001—04

2016-07-12

張楊洋(1991- )，男，陜西省興平縣人，助理工程師，碩士，主要從事巖土工程支護與數(shù)值計算工作.

TU45;TV223.31

A

10.3969/j.issn.1006-2610.2016.05.001