張月征 紀洪廣 宋朝陽
(金屬礦山高效開采與安全教育部重點實驗室,北京100083)
地下工程巖體在開挖擾動后發(fā)生的動力破壞與其開挖前的初始穩(wěn)定性相關,巖體的初始穩(wěn)定狀態(tài)是由巖體自身材料屬性和初始地應力場共同決定的。初始地應力場由巖體自重和地質構造運動所引起,與巖體特性、歷史構造現(xiàn)象、正在發(fā)生的構造運動密切相關,是塑造地質構造、誘發(fā)地震的動力源泉,同時也嚴重影響著人類活動進行的地下工程的穩(wěn)定性,在高地應力地區(qū)巖體開挖后常會出現(xiàn)巖爆、沖擊地壓、鉆孔縮徑等地質災害[1-2]。對于巖石材料而言,當其受力狀態(tài)達到某種極限狀態(tài)時才會發(fā)生破壞,受自身材料屬性、應力狀態(tài)、加載方式等影響會產生脆性斷裂、脆性剪切、延性、弱面剪切等破壞方式,應力是造成巖石破裂的最主要因素。無論何種破壞方式應力狀態(tài)都決定著巖石是否接近破壞失穩(wěn)極限,當應力狀態(tài)越接近破壞極限時,破壞失穩(wěn)危險程度越高。張鏡劍等[3]在研究分析巖爆判據(jù)中歸納出了8 種判據(jù),其中國外學者提出了E.Hoek 方法(最大主應力與單軸抗壓強度之比σmax/Rc)和Russense 判據(jù)方法(硐室最大切向應力與單軸抗壓強度之比σθ/Rc)。徐林生等[4]在二郎山隧道巖爆研究實踐中提出了改進的“σθ/Rc判據(jù)法”及分級標準。谷明成建立的綜合巖爆判據(jù)中也提到了應力標準(σθ/Rc≥0.3 ),張宏偉等[3-4]則采用應力集中系數(shù)來描述地質動力特征。應力場是影響巖體初始穩(wěn)定性的主要的影響因素,因此,量化分析評價巖體初始穩(wěn)定性需要從其所處應力狀態(tài)出發(fā),建立評價指標,分析巖體初始穩(wěn)定性與動力災害發(fā)生的相關性,進而為災害發(fā)生提供危險性預測。
本研究從強度理論和初始應力場出發(fā),采用莫爾-庫倫原理推導出了基于應力狀態(tài)分析的巖石初始穩(wěn)定性描述方法和評價指標——“莫-庫應力”,并采用該指標對地下開采前的巖體三維初始地應力場狀態(tài)進行評價,統(tǒng)計分析巖體初始穩(wěn)定性與沖擊地壓動力災害發(fā)生的相關性,進行巖體穩(wěn)定性分級,為地下工程災害的防治提供參考。
巖體初始應力狀態(tài)由自重應力場和構造應力場疊加而成[5]。由地球重力場引起巖體內部垂直應力分量,形成巖體單元的自重應力狀態(tài),見式(1):
式中,γ 為上覆巖體重度;H 為巖體單元深度。
構造應力場由地質構造運動產生,水平方向構造運動如板塊移動碰撞,對巖體構造應力的形成起控制作用,構造應力以水平應力為主,早期學者認為水平應力只與垂直應力有關,見式(2):
式中,λ 為側壓力系數(shù);μ 為泊松比。
隨著試驗儀器的發(fā)展,目前可通過現(xiàn)場測量方法獲得地層巖體中較為準確的地應力分布。
斷層、褶皺這些保留在巖體中的地質構造與其形成時的應力又密切關系,構造應力場方向可根據(jù)地質力學的方法加以判斷[6]。如圖1:正斷層的自重應力為最大主應力,最小主應力與斷層走向正交;對于逆斷層,最大主應力與斷層走向正交;對于平移斷層,自重應力是中間主應力,最大主應力與斷層走向呈30°~45°夾角,最小主應力與最大主應力均為水平方向;褶皺的最大主應力方向正交于褶皺軸部走向。
圖1 地質構造與主應力方向Fig.1 Geological structure and the main stress direction
在地應力場動力特征研究方面,張宏偉等[7-9]以板塊構造學為基礎,根據(jù)地形地貌的基本形態(tài)決定地質構造形式原理,發(fā)展了地質動力區(qū)劃學說,并將其引入到礦井沖擊地壓、瓦斯突出等動力災害的危險預測當中。采用應力集中系數(shù)K 來劃分高應力區(qū)(K>1.2 )、應力梯度區(qū)(0.8 ≤K ≤1.2 )、低應力區(qū)(K <0.8 ),確定可能發(fā)生礦井動力現(xiàn)象的危險區(qū)域,評估沖擊危險程度,對斷層構造發(fā)育的地區(qū)的動力災害防治起到了重要的作用。K 的計算公式為
地下工程中巖體的初始應力狀態(tài)是巖體在開挖后發(fā)生失穩(wěn)破壞的重要影響因素,一般以最大主應力、主應力方向、應力集中程度描述地應力場狀態(tài),判斷可能發(fā)生破壞的危險區(qū)域。但了解這些因素還不能夠完全確定危險區(qū)域是否容易發(fā)生破壞,例如,在高應力區(qū),如果巖體處于三向受壓狀態(tài)未必產生破壞,而在低應力區(qū),如果應力處在最不利的拉壓組合狀態(tài)下也可能易發(fā)生破壞。因此,需要一種方法或指標來描述某一應力狀態(tài)下主應力的組合方式,量化評價巖石在初始應力場狀態(tài)下的穩(wěn)定性,進而確定巖體接近破壞極限的程度。
對于單軸受力情況,巖石的受力狀態(tài)可根據(jù)相對應力,即巖石所承受應力與峰值強度應力之比予以評價,而對于三軸受力狀態(tài)下巖石所處的應力狀態(tài)評價還需要綜合考慮多個應力的影響。而巖石強度理論中以巖石所處應力狀態(tài)為參量研究其與巖石破壞的關系,當其所處應力狀態(tài)越接近極限應力狀態(tài),則其發(fā)生破壞的可能性越高,巖石的穩(wěn)定性也就越差,發(fā)生破壞的危險性也就越高。經(jīng)典的莫爾-庫倫強度理論為三維應力狀態(tài)的評價提供了基礎依據(jù),該理論體系認為巖石是在不同的正應力和剪應力組合作用下喪失承載能力,巖石的強度值與中間主應力的大小無關,如圖2 所示。
圖2 莫爾應力圓與應力狀態(tài)Fig.2 Mohr stress circle and the stress state
莫爾-庫倫強度理論數(shù)學表達式為
式中,τ 為正應力σ 作用下的極限剪應力,MPa;c 為巖石的內聚力,MPa;θ 為巖石的內摩擦角,(°)。
從圖2 中可知,當最小主應力σ3確定的時候,依據(jù)幾何關系可得出必存在唯一的最大主應力峰值σ峰值與最小主應力σ3所組成的極限應力圓與強度線相切。在此,設巖石處在某一應力狀態(tài)(σ1,σ3)下,定義一參數(shù)Wσ來表征應力狀態(tài)與極限應力狀態(tài)在圖2 中的幾何關系:
式中,Wσ為某一應力狀態(tài)下莫爾應力圓與極限狀態(tài)下應力圓的半徑之比,在此定義其為莫爾-庫倫應力;r0為與直線相切極限狀態(tài)的莫爾應力圓半徑;r'為某一應力狀態(tài)下的莫爾應力圓半徑。
由圖2 幾何關系可知,Wσ越小,莫爾圓越偏離強度線,巖石發(fā)生破壞所需的擾動越小;Wσ的物理意義反映的是一種相對應力狀態(tài),即應力狀態(tài)與極限應力狀態(tài)的關系。Wσ值越大,越接近極限狀態(tài),穩(wěn)定性越差。
為量化這一參數(shù),作者在文獻[10]中做了詳細的推導,得出了Wσ數(shù)學計算式[10]為
式中,σ1為最大主應力;σ3為最小主應力;c 為內聚力;θ 為內摩擦角;ξ 為σ峰值系數(shù),ξ = σ峰值/σ理論峰值。
從解析推導過程中發(fā)現(xiàn),基于莫爾-庫倫強度準則建立的參量Wσ表征了某一應力狀態(tài)下主應力之間的組合方式以及這種組合方式下其應力狀態(tài)接近極限狀態(tài)的強度。
(1)單軸壓縮。當σ3= σ2= 0 ,σ1>0 時,巖石處在單軸壓縮情況下,
式中,σ壓峰值等于巖石的單軸抗壓強度。
(2)單軸拉伸。當σ3= σ2= 0 ,σ1<0 時,巖石處在單軸拉伸情況下,
式中,σ拉峰值等于巖石的單軸抗拉強度。
(3)三向受壓。當σ1>σ2>σ3>0 時,巖石處在三向受壓情況下,0 ≤Wσ≤1 ,當Wσ為0 時說明巖石處在三向等力壓縮情況,無論應力值有多高,其狀態(tài)是穩(wěn)定的;當Wσ為1 時,應力狀態(tài)達到極限狀態(tài),這種情況下主應力組合最不利,巖石發(fā)生破壞;當Wσ從0 向1 變化時,巖石從穩(wěn)定狀態(tài)向失穩(wěn)狀態(tài)過渡,這期間巖石經(jīng)歷彈性階段、塑性階段、破壞階段,由于巖石材料的脆性性質,當其越接近1 時,同樣擾動條件下,其發(fā)生破壞的可能性越大。
(4)受壓與拉伸組合。當σ1>σ2>σ3,σ3<0,σ1>0 時,在受壓與拉伸應力組合下,巖石材料最容易破壞。
無論是單軸壓縮、單軸拉伸、三向受壓、受壓與拉伸組合下的應力組合,Wσ都能表現(xiàn)出其接近極限的程度,當Wσ越接近1,巖石越接近破壞極限,穩(wěn)定性越差。由莫爾-庫倫強度理論衍生而來的參量Wσ能夠描述巖石材料在各種應力組合狀態(tài)下其接近破壞極限的程度。
大量試驗資料研究表明,巖石受力破壞過程中一般經(jīng)歷彈性階段—塑性階段—破壞階段,一般情況下材料在應力平均達到70%左右的時候開始產生擴容現(xiàn)象,即巖石內部損傷裂紋開始發(fā)展,材料進入塑性發(fā)展階段時期,內部裂紋進一步貫通,最終導致巖石整體破壞失穩(wěn)[11-13]。遵循巖石破壞這一物理過程,根據(jù)Wσ將巖體穩(wěn)定性劃分為5 個等級,劃分依據(jù)見表1。
表1 基于Wσ 的巖石穩(wěn)定性等級Table 1 the Level of rock stability based on Wσ
穩(wěn)定:巖石處于受力初級階段和變形彈性階段前期,破壞還需要大的外力作用,安全,穩(wěn)定級別最高。
較穩(wěn)定:巖石處于變形彈性階段后期,破壞需要較大的外力作用,偏安全,穩(wěn)定級別較高。
中等穩(wěn)定:巖石處于變形彈性與塑性過渡階段,內部損傷裂紋開始產生,偏危險,穩(wěn)定級別中等。
較低穩(wěn)定:巖石處于變形塑性發(fā)展階段,內部裂紋發(fā)育速度加快,破壞需要較小的外力作用,偏高危險,穩(wěn)定性較低。
低穩(wěn)定:巖石處于變形塑性發(fā)展和臨界破壞階段,內部裂紋快速發(fā)育貫通,破壞需要較小的外力作用,高危險,穩(wěn)定性最差。
深部巖體在開挖前受賦存條件限制完整性良好,巖體初始穩(wěn)定性主要由地應力場控制。Wσ越高的巖體,其達到失穩(wěn)破壞所需要的外部擾動門檻值越小,也就是說其受到外部擾動因素發(fā)生后越容易達到極限破壞狀態(tài),進而越容易引發(fā)相應的地下工程動力災害?;诖耍蓪σ這一新參量引入到初始地應力場巖體穩(wěn)定性的評價當中,表征基于應力狀態(tài)評價的穩(wěn)定性等級的分布情況,進而為人類工程活動提供必要的地質環(huán)境信息并預測工程活動可能產生的地質動力效應。本研究選取典型褶皺構造發(fā)育的甘肅省華亭煤田區(qū)域內地應力場進行評價分析,劃分巖體穩(wěn)定分級區(qū),通過沖擊地壓事件的統(tǒng)計分析,研究巖體初始穩(wěn)定性與沖擊地壓發(fā)生的相關性。
如圖3 所示,華亭煤田區(qū)域內地質環(huán)境復雜,主要受褶曲地質構造控制,向斜呈北北西—南南東向展布的“S”型構造形態(tài),東緩西陡,中間寬緩,南北兩段收斂,形似紡綞形的復式不對稱向斜構造,受區(qū)域構造應力的扭動作用形成。區(qū)域內煤層厚,煤層厚度變化大,受構造影響,煤層起伏,無斷層,呈多個向斜、背斜相互交替連接的空間形態(tài),整個煤田位于構造應力集中區(qū),煤層底板和頂板皆為砂巖。
該區(qū)域內地應力測試結果表明,華亭煤田受水平構造應力控制為主,最大主應力、最小主應力、垂直主應力與深度成線性比例關系,深度越深其值變化越大,最大主應力值受深度影響最大,地應力與深度的關系:
建模過程中,根據(jù)該區(qū)域內地形高程、地層、走向、煤巖層分布、鉆孔、勘探線等多種地質信息,采用Surpac 軟件將這些地質信息數(shù)字化,建立起三維地質模型,見圖4。以地質模型為基礎,建立了數(shù)值計算模型,并根據(jù)地應力測試結果和巖石物理力學參數(shù),見表2,采用FLAC3D數(shù)值計算軟件進行了地應力場反演,獲取了該區(qū)域內三維地應力分布圖,見圖5。
表2 巖石力學參數(shù)Table 2 Rock mechanics parameters
圖3 計算模型區(qū)域(單位:m)Fig.3 The area of calculation model
圖4 基于Surpac 軟件建立的地質模型Fig.4 The geological model based on software Surpac
圖5 基于FLAC3D軟件建立的反演模型(垂直應力場)Fig.5 The inversion model based on FLAC3D software (vertical stress field)
華亭煤田區(qū)域內有華亭、硯北、陳家溝、山寨4個礦井同時進行開采工作,受地質條件和開采條件等眾多因素的影響,2008 年以來,隨著開采量和開采深度的增加礦井內多次出現(xiàn)沖擊地壓現(xiàn)象,尤其以華亭、硯北兩礦井最為嚴重,最大級別沖擊事件誘發(fā)當?shù)?.9 級地震。通過對來壓事件的數(shù)據(jù)分析,區(qū)域內的沖擊地壓多顯現(xiàn)為底板來壓型,能量級別高,破壞嚴重,因此,提取了煤層底板巖層的地應力信息進行分析,繪制了最大主應力場和巖體彈性能分布圖,見圖6 和圖7。
圖6 最大主應力分布Fig.6 The maximum principal stress
圖7 巖體彈性能分布Fig.7 The elastic energy distribution of rock mass
(1)巖體最大主應力與彈性能特征。區(qū)域最大主應力與巖體彈性能分布規(guī)律類似,兩者分布與底板構造形態(tài)及埋深密切相關,向斜軸內最大主應力和能量富集程度高,背斜軸附近富集程度低,在褶皺翼部呈過渡趨勢;埋深大的地區(qū)最大主應力富集程度低,埋深小的地區(qū)富集程度低。
(2)Wσ分布。根據(jù)式(6),計算獲得了底板莫-庫應力分布情況,見圖8。莫-庫應力與底板巖層構造形態(tài)密切相關,在向斜軸和背斜軸附近莫-庫應力值最高,褶皺翼部莫-庫應力值稍低,在地層平緩地區(qū)莫-庫應力值最低。
對比圖6、圖7、圖8,在最大主應力或巖體賦存能量富集程度高的區(qū)域其應力狀態(tài)未必是危險程度最高的,同樣最大主應力或巖體賦存能量富集程度低的區(qū)域應力狀態(tài)未必危險程度最低。從理論上講,應力狀態(tài)的危險程度由應力大小和組合方式兩者共同決定;如果巖石處在三向等力壓縮時,無論力有多大,積聚能量有多高,其應力狀態(tài)都是穩(wěn)定安全的;如果巖石處在最不利的壓縮拉伸組合狀態(tài)時,則無需太大的力和太高的能量即可破壞,應力狀態(tài)是最不穩(wěn)定不安全的。因此單獨靠最大主應力、能量單個參量的變化來判斷危險程度是不夠完備的,而莫-庫應力考慮了應力組合方式在里面,則能夠彌補這一不足。
為獲取巖石初始穩(wěn)定性與沖擊地壓類動力災害發(fā)生的關系,選取2011 年4 月至2012 年12 月發(fā)生在硯北煤礦和華亭煤礦的沖擊事件進行分析,利用GIS 軟件,將沖擊事件以圓點形式投影至華亭煤田巷道模型內,對其發(fā)生特征與基于Wσ劃分的應力危險區(qū)的關系進行分析。如圖9 所示,圖中圓點大小代表沖擊地壓級別大小,圓點越大沖擊級別越高。
硯北煤礦:這段時間內該礦的工作狀態(tài)為250204 工作面回采工作和250203 工作面巷道掘進工作,期間受回采影響250204 工作面內共發(fā)生42 次沖擊事件,受掘進影響250203 工作面內共發(fā)生40 次沖擊事件。
華亭煤礦:工作狀態(tài)為250104 工作面回采工作和250105 工作面巷道掘進工作,期間受回采影響250104 工作面內共發(fā)生9 次沖擊事件,受掘進影響250105 工作面內共發(fā)生31 次沖擊事件。
圖9 和表3 數(shù)據(jù)表明,大部分沖擊事件發(fā)生在低穩(wěn)定區(qū)和較低穩(wěn)定區(qū)范圍內,其中硯北煤礦25204 工作面靠近背斜軸附近的沖擊事件分布位置與低穩(wěn)定區(qū)分布走向高度保持一致。這說明沖擊地壓事件的空間分布與基于Wσ劃分的巖石初始穩(wěn)定性密切相關,呈現(xiàn)以下特征:①沖擊事件多分布在低穩(wěn)定區(qū)和較低穩(wěn)定區(qū);②沖擊事件分布趨勢與穩(wěn)定性低的區(qū)域分布及走向趨勢相一致;③巖石穩(wěn)定性越差的地區(qū),發(fā)生沖擊地壓事件的概率越大,沖擊事件的級別也越高。
表3 來壓事件統(tǒng)計分析Table 3 Statistical analysis on rock burst
如圖10 所示,基于莫爾-庫倫準則的動力災害危險區(qū)劃分方法在分析過程中,蘊含了巖石材料信息(強度、泊松比、彈性模量)、空間地質信息(地形地貌、地層、地質構造等)、實測地應力等多維信息;在此應力場基礎上建立的Wσ則綜合考慮了內聚力、內摩擦角、最大主應力、最小主應力之間的綜合關系。
圖10 基于莫-庫準則的巖體初始穩(wěn)定性分析技術路線Fig.10 Technical roadmap of the initial stability of rocks based on Wσ
通過對采用Wσ劃分巖體初始穩(wěn)定區(qū)及對沖擊地壓事件進行統(tǒng)計分析表明,采用該方法能夠為沖擊地壓這一類地下工程動力災害提供合理的前瞻性預測,莫-庫應力越高的地區(qū),在開采過程中巖體發(fā)生動力失穩(wěn)引發(fā)災害的可能性越高,因此根據(jù)基于莫爾庫倫準則建立的Wσ進行深部巖體初始穩(wěn)定性劃分并進行動力危險區(qū)預測是合理有效可行的。
(1)基于莫爾-庫倫強度理論建立了評價應力狀態(tài)危險性的Wσ指標,物理本質上反映的是巖石破壞前彈塑性階段某一應力狀態(tài)下巖石材料體接近極限破壞強度的程度。
(2)Wσ可以表達單軸壓縮、單軸拉伸、三向受壓、壓縮-拉伸等模式下的巖石應力狀態(tài)接近極限狀態(tài)的程度,Wσ越接近1,巖石越接近破壞極限,巖體的穩(wěn)定性越差。
(3)運用Wσ評價巖體初始地應力場狀態(tài),能夠表達出巖體所處應力狀態(tài)與極限狀態(tài)的接近程度,尤其是地質構造發(fā)育區(qū)域復雜應力場的狀態(tài),有助于判斷初始應力場下的巖體穩(wěn)定性。
(4)Wσ越高的巖體,其達到失穩(wěn)破壞所需要的外部擾動值越小。可依據(jù)Wσ劃分穩(wěn)定、較穩(wěn)定、中等穩(wěn)定、較低穩(wěn)定、低穩(wěn)定5 級穩(wěn)定區(qū),巖體初始穩(wěn)定性越差的區(qū)域發(fā)生沖擊地壓動力災害的概率越高。據(jù)此可基于地質環(huán)境應力信息對巖體動力危險程度進行評價,為預測工程活動所引發(fā)的沖擊地壓等地質動力效應提供參考。
(5)運用Wσ進行原巖穩(wěn)定區(qū)劃分具有實用意義,然而沖擊地壓發(fā)生機制復雜,誘發(fā)與地質環(huán)境、開采活動等多種因素有關,沖擊地壓的精確預測還需與其他影響因素綜合分析。
[1] 蔡美峰,何滿潮,劉東燕.巖石力學與工程[M].北京:科學出版社,2002:129-139.
Cai Meifeng,He Manchao,Liu Dongyan.Rock Mechanics and Engineering[M].Beijing:Science Press,2002:129-139.
[2] 李 曉,李守定,陳 劍,等. 地質災害形成的外動力耦合作用機制[J].巖石力學與工程學報,2008,27(9):1792-1806.
Li Xiao,Li Shouding,Chen Jian,et al.Coupling effect mechanism of endogenic and exogenic geological processes of geological hazards evolution[J]. Chinese Journal of Rock Mechanics and Engineering,2008,27(9):1792-1806.
[3] 張鏡劍,傅冰駿.巖爆及其判據(jù)和防治[J]. 巖石力學與工程學報,2008,27(10):2034-2042.
Zhang Jingjian,F(xiàn)u Bingjun. Rookburst and its criteria and control[J].Chinese Journal of Rock Mechanics and Engineering,2008,27(10):2034-2042.
[4] 徐林生,王蘭生. 二郎山公路隧道巖爆發(fā)生規(guī)律與巖爆預測研究[J].巖土工程學報,1999,2l(5):569-572.
Xu Linsheng,Wang Lansheng.Study on the laws of rockburst and its forecasting in the tunnel of Erlang Mountain road[J].Chinese Journal of Geotechnical Engineering,1999,2l(5):569-572.
[5] 侯明勛,葛修潤.巖體初始地應力場分析方法研究[J]. 巖土力學,2007,28(8):1626-1630.
Hou Mingxun,Ge Xiurun. Study on fitting analysis of initial stress field in rock masses[J]. Rock and Soil Mechanics,2007,28(8):1626-1630.
[6] 沈明榮,陳建峰. 巖體力學[M]. 上海:同濟大學出版社,2006:120-124.
Shen Mingrong,Chen Jianfeng. Rock Mechanics[M]. Shanghai:Tongji University Press,2006:120-124.
[7] 張宏偉,韓 軍,宋衛(wèi)華.地質動力區(qū)劃[M]. 北京:煤炭工業(yè)出版社,2009:198-202.
Zhang Hongwei,Han Jun,Song Weihua. Geodynamics Zoning[M].Beijing:Coal Industry Press,2009:198-202.
[8] 張宏偉.地質動力區(qū)劃方法在煤與瓦斯突出區(qū)域預測中的應用[J].巖石力學與工程學報,2003,22(4):621-624.
Zhang Hongwei.Application of geo-dynamic division method in prediction of coal and gas outburst region[J].Chinese Journal of Rock Mechanics and Engineering,2003,22(4):621-624.
[9] 王志輝,張宏偉,石永生. 地質動力區(qū)劃在新汶潘西礦的應用[J].遼寧工程技術大學學報,2005,24(S):1-4.
Wang Zhihui,Zhang Hongwei,Shi Yongsheng.Application of geological dynamic division of Panxi mine[J].Journal of Liaoning Technical University:2005,24(S):1-4.
[10] 張月征,紀洪廣,侯昭飛.基于莫爾-庫倫強度理論的巖石沖擊危險性判據(jù)[J].金屬礦山,2014(11):138-142.
Zhang Yuezheng,Ji Hongguang,Hou Zhaofei.Instability criterion of rockburst risk based on Mohr-Coulomb strength theory[J]. Metal Mine,2014(11):138-142.
[11] 張建國,張強勇,楊文東,等.大崗山水電站壩區(qū)初始地應力場反演分析[J].巖土力學,2009,30(10):3071-3078.
Zhang Jianguo,Zhang Qiangyong,Yang Wendong,et al.Regression analysis of initial geostress field in dam zone of Dagangshan hydropower station[J].Rock and Soil Mechanics,2009,30(10):3071-3078.
[12] 宋衛(wèi)華,張宏偉.斷層滑動準則和穩(wěn)定性評價及構造應力場的反演[J].巖石力學與工程學報,2008,27(S2):3470-3475.
Song Weihua,Zhang Hongwei. Assessment on faults slide criterion and stability and back-analysis of tectonic stress field[J]. Chinese Journal of Rock Mechanics and Engineering,2008,27(S2):3470-3475.
[13] 王金安,劉 航,李 鐵.臨近斷層開采動力危險區(qū)劃分數(shù)值模擬研究[J].巖石力學與工程學報,2007,26(1):28-35.
Wang Jin'an,Liu Hang,Li Tie. Study on numerical simulation of dynamic risk regionalization during exploitation approaching to faults[J]. Chinese Journal of Rock Mechanics and Engineering,2007,26(1):28-35.