露天邊坡局部失穩(wěn)特征的數(shù)值分析及監(jiān)測

張擁軍 楊登峰 王金濤 陳關(guān)平

(1.青島理工大學土木工程學院，山東 青島 266033；2.中國礦業(yè)大學(北京)力學與建筑工程學院，北京 100083)

·安全與環(huán)保·

露天邊坡局部失穩(wěn)特征的數(shù)值分析及監(jiān)測

張擁軍1楊登峰2王金濤1陳關(guān)平1

(1.青島理工大學土木工程學院，山東 青島 266033；2.中國礦業(yè)大學(北京)力學與建筑工程學院，北京 100083)

以金川集團露天石英石礦邊坡監(jiān)測項目為依托，利用監(jiān)測數(shù)據(jù)成果，以臺階局部滑坡部位作為研究對象，通過分析其臺階周圍及滑坡剖面上臺階監(jiān)測點的累積水平位移值，研究了采動卸荷效應造成的集中應力反復作用下邊坡的變形破壞機制。利用RFPA強度折減法，數(shù)值模擬過程再現(xiàn)了1 768 m平臺開采過程中典型剖面多臺階的漸進失穩(wěn)過程，宏觀裂縫在細觀上的“破裂-應力轉(zhuǎn)移-破裂”的循環(huán)中不斷發(fā)育，最終形成貫通的滑動面。破壞的整個過程中伴隨著由應力集中、釋放、轉(zhuǎn)移過程造成的應變局部化現(xiàn)象。監(jiān)測結(jié)果表明：局部化損傷發(fā)展過程是從采動平臺向上部臺階邊坡巖體轉(zhuǎn)移，各臺階邊坡內(nèi)部7～11 m位置處受荷載轉(zhuǎn)移影響程度較大。

邊坡 變形局部化 強度折減法 監(jiān)測 數(shù)據(jù)分析

邊坡穩(wěn)定性分析的理論與方法的發(fā)展已經(jīng)有了很長的歷史，涌現(xiàn)出了大量經(jīng)典的理論計算分析方法[1-6]，并在實踐中得到了有效的應用。但是由于邊坡工程的復雜性、影響因素的不確定性以及邊坡變形破壞演化的時空效應，單純的理論分析計算已經(jīng)不再適合邊坡穩(wěn)定性分析的需求。近年來，隨著數(shù)值模擬方法及監(jiān)測技術(shù)的發(fā)展，邊坡數(shù)值模擬分析方法及監(jiān)測預警技術(shù)在工程中得到了越來越廣泛的應用。

鄭穎人[4]、唐春安等[5]將強度折減法的基本原理引入到巖土工程的數(shù)值分析方法中，為邊坡的穩(wěn)定性研究提供了一種新的方法。陳國慶[6]利用動態(tài)強度折減法與整體強度折減法相結(jié)合計算邊坡漸進失穩(wěn)過程中的動態(tài)安全系數(shù)，達到對邊坡失穩(wěn)全過程的分析和調(diào)控的目的。許強[7]提出斜坡變形時間演化的三階段規(guī)律和斜坡裂縫空間演化的分期配套特性。何滿潮[8]提出滑體和滑床相對運動的力學監(jiān)測原理，開發(fā)了對滑動力和抗滑力進行監(jiān)測的滑坡地質(zhì)災害遠程監(jiān)測預報系統(tǒng)。張振華[9]提出了基于設(shè)計安全系數(shù)和破壞模式的邊坡開挖過程中動態(tài)變形監(jiān)測預警指標的研究思路和方法。王東[10]將位移速度變量看作是一隨機變量，提出基于Lilliefors檢驗的邊坡臨滑時刻的動態(tài)識別方法。

對巖體破壞現(xiàn)象的研究發(fā)現(xiàn)，不同尺度巖體的破裂都有局部化破壞特征。露天邊坡災害的發(fā)生往往是集中在采動影響的某些部位[11]。因此，注意捕捉巖石或巖體結(jié)構(gòu)變形局部化信息并掌握其發(fā)展趨勢，就有可能提高邊坡災害監(jiān)測預報的可能性。

本研究以金川露天石英石礦邊坡監(jiān)測項目為依托，以臺階局部滑坡部位作為研究對象，通過分析其周圍臺階及與滑坡剖面相鄰的監(jiān)測點的累積水平位移值并結(jié)合數(shù)值模擬的方法來研究采動卸荷效應造成的集中應力作用下邊坡的變形破壞機制，分析邊坡開采過程中的損傷累積導致的變形失穩(wěn)特征。為礦山邊坡的破壞分析提供理論依據(jù)，提高滑坡監(jiān)測預警的可靠性。

1 工程地質(zhì)條件

1.1 工程概況

金川集團露天石英石礦位于甘肅省金昌市。目前，石英石礦露天采場已形成1 856、1 846、1 836、1 826、1 816、1 792、1 780、1 768、1 756 m共8個臺階，東西長1 100 m、南北寬500 m、高近90 m的山坡型露天采場，1 756 m水平以下的1 744、1 732、1 720、1 708、1 696 m等5個水平將進入凹陷開采。目前許多區(qū)段邊坡已產(chǎn)生不同程度的變形和破壞，主要表現(xiàn)：1 780 m臺階中段出現(xiàn)局部坍塌；1816 m臺階西部出現(xiàn)許多的裂縫，且不斷擴大；1 792 m臺階東部出現(xiàn)滾石。隨著采場的延深、雨水作用及地震、爆破震動等因素的影響，邊坡穩(wěn)定問題日益突顯。

1.2 地質(zhì)條件

礦區(qū)揭露地層依次為第四系，石炭—二疊系，下古生界(上、下組)。由于礦區(qū)受區(qū)域構(gòu)造控制和南西、北東向的擠壓，使含礦巖系形成NWW向的緊閉褶皺和斷裂，并有一系列北東和北西向的剪切斷裂伴生，加上礦區(qū)受開采影響，地表巖石普遍松動，裂隙發(fā)育，礦區(qū)巖體破碎。礦區(qū)石炭統(tǒng)—二疊系(C3-P)巖石質(zhì)量指標RQD值在0.4～0.5之間，完整性較好，其他地層RQD值在0～0.1之間,巖體極其破碎，工程地質(zhì)條件極差。礦體RQD值在0.3～0.5之間，完整性較好。礦區(qū)氣候干旱，年降水量122 mm，蒸發(fā)量2 837 mm，且無常年性地表水體。

1.3 監(jiān)測方案設(shè)計

針對礦區(qū)干旱、少雨等特點，經(jīng)綜合論證決定采用固定式測斜儀對多臺階邊坡的變形狀態(tài)進行監(jiān)測，綜合考慮巖體內(nèi)部斷面位置、節(jié)理發(fā)育情況，施工難易程度，安全等多項因素，在采場北部邊幫1 768、1 780、1 792、1 816、1 826 m 5個平臺上拉網(wǎng)式布置21個監(jiān)測孔共84個測斜儀(位置如圖1所示)，測斜儀通過測斜管裝入監(jiān)測孔內(nèi)，根據(jù)鉆孔取芯資料，將測斜管底部嵌入邊坡體穩(wěn)定的基巖內(nèi)，并以此作為基準進行位移計算，通過遠程數(shù)據(jù)采集系統(tǒng)，對邊坡體內(nèi)部變形進行長期監(jiān)測。

圖1 監(jiān)測孔位置分布圖Fig.1 Location map of monitor holes▲—數(shù)據(jù)采集及傳輸系統(tǒng)；○—鉆孔監(jiān)測點

2 數(shù)值模擬分析

為研究采動卸荷作用導致的集中應力釋放后的傳遞路徑及對礦區(qū)各臺階的影響程度，根據(jù)工程實際地質(zhì)條件，滑坡面附近實際存在的開采臺階狀況，建立數(shù)值模型。

2.1 模型建立

在考慮巖土體材料非均勻性的基礎(chǔ)上，基于RFPA細觀尺度上的本構(gòu)模型，將模擬問題視為靜力條件下的二維平面應變問題。以石英石礦區(qū)地質(zhì)條件作為模型設(shè)計的依據(jù)，建立了1 768 m臺階開挖后的數(shù)值模型(圖2所示)。

圖2 1 768 m平臺開挖前后RFPA數(shù)值模型示意圖Fig.2 RFPA numerical model schematic diagram before and after the excavation of 1768 m platform

數(shù)值模型沿水平方向取130 m，垂直方向取95 m，共設(shè)4層巖層，劃分為49 400個單元。模型邊界條件：垂直方向自重加載，無外荷載，模型兩側(cè)為限制水平方向位移的滑動支座，可垂直移動；底部為限制垂直方向位移的滑動支座。采用RFPA強度折減法[12]，強度折減系數(shù)取0.01，加載總步數(shù)為100步。模擬1 768 m臺階礦體開采后邊坡的卸荷變形破壞情況。

根據(jù)金川石英石露天礦實測工程地質(zhì)數(shù)據(jù)，本數(shù)值模擬選取的各地層巖體力學參數(shù)如表1所示。

表1 各地層巖石物理與力學參數(shù)Table 1 Physical and mechanical parameters of rocks at each layer

2.2 數(shù)值模擬結(jié)果分析

圖3(a)是通過數(shù)值模擬得出的1 768 m臺階開挖后邊坡卸荷損傷漸進失穩(wěn)的彈性模量圖，圖3(b)是對應的聲發(fā)射圖。在聲發(fā)射圖中，黑色圓圈表示拉伸破壞，白色圓圈表示剪切破壞。

分析圖3可知，邊坡強度折減過程中，在拉剪應力持續(xù)作用下，損傷逐步累積，隨著損傷累積，小尺度小范圍內(nèi)的巖體破裂持續(xù)發(fā)生，裂隙發(fā)育(Step 24-02所示)，由于坡體內(nèi)的裂隙、結(jié)構(gòu)面及巖石的非均勻性等因素影響，導致坡體中應力呈現(xiàn)出不均勻。折減過程中1 780 m臺階坡腳處出現(xiàn)局域性的應力集中，巖石聲發(fā)射的分布也由無序向有序發(fā)展，逐漸向該區(qū)域集中，當應力集中超過該處巖體的強度時，巖體將不再穩(wěn)定，必然發(fā)生局部破壞(Step 25-06所示)，出現(xiàn)局部失穩(wěn)滑動面，伴有小規(guī)模崩塌。并伴隨發(fā)生應力釋放，隨后應力發(fā)生轉(zhuǎn)移和應力的重新調(diào)整，應力逐步向1 780 m臺階的鄰近區(qū)域的1 792 m臺階擴展，由于新增應力作用，導致1 792 m臺階處的巖體也出現(xiàn)應力集中，使原先沒有超過巖體的強度值轉(zhuǎn)變?yōu)閼χ党^了巖體的強度值，損傷逐步累積，裂紋萌生(Step 25-08所示)，直至巖體發(fā)生破壞。此后重復上述過程，應力再次發(fā)生轉(zhuǎn)移和重新調(diào)整，導致1 816 m臺階的應力集中，應力值超過了巖體的強度值，臺階裂紋貫通，發(fā)生局部破壞(Step 25-11所示)。宏觀裂縫在這個細觀上“破裂—應力轉(zhuǎn)移—破裂”的循環(huán)中不斷發(fā)育。在此過程中，微裂紋逐步沿臺階向上發(fā)生損傷擴展，隨著轉(zhuǎn)移的應力值逐步減小，裂紋擴展的速度和臺階發(fā)生破壞的規(guī)模逐漸減小，體現(xiàn)了應變局部化特征。隨著過程的繼續(xù)，邊坡最終將有2種可能，一種是破裂面完全貫通，損傷巖體沿著滑動面整體運動；另一種是破壞面沒有完全貫通，破壞擴展到某一區(qū)域后停止擴展，其上部的應力應變未超過巖體的強度值。由數(shù)值模擬可知，最下部臺階坡腳處雖有拉剪裂隙產(chǎn)生，但是破壞較小，并沒有形成貫通的破壞區(qū)域，因此可以判斷卸荷巖體應力集中釋放之后其應力轉(zhuǎn)移主要是向上部轉(zhuǎn)移，對下部巖體影響較小。

圖3 1 768 m臺階開挖后的彈性模量動態(tài)變化規(guī)律和聲發(fā)射規(guī)律Fig.3 The dynamic variation rules and acoustic emission rules of elastic modulus after excavation of 1 768 m platform

3 監(jiān)測成果分析

在監(jiān)測系統(tǒng)中，取坡體內(nèi)部指向邊坡臨空面方向的位移為正，反之為負。通過監(jiān)測發(fā)現(xiàn)1 868 m臺階邊坡于2013年9月28日，靠近1-B監(jiān)測點處發(fā)生過一次較大規(guī)模滑坡，取監(jiān)測點1-B、滑坡點附近左側(cè)監(jiān)測點1-C、上部剖面距滑坡面較近的1 780 m臺階監(jiān)測點2-C進行分析，由于監(jiān)測數(shù)據(jù)量較大，取具有代表性的監(jiān)測點在2013-04-01～2014-05-15日期間的監(jiān)測數(shù)據(jù)進行分析，討論邊坡卸荷過程中的巖體破壞變形規(guī)律。

3.1 監(jiān)測數(shù)據(jù)分析

圖4是離滑坡地點最近的監(jiān)測點1-B在2013-04-01～2013-09-28日滑坡前的位移累積曲線，滑坡時滑落的巖土體將監(jiān)測孔破壞，觀察圖4曲線變化規(guī)律可知，在2013-05-30日前曲線變化較平緩，呈緩慢上升趨勢，之后位移量迅速上升，進入加速變形階段，截止到2013-09-15日時位移量達到106.5 mm，此后位移增速加快，幾乎呈指數(shù)形式增長，監(jiān)測到的位移最大值達到了208 mm。經(jīng)綜合分析，于9月25日發(fā)出滑坡預警。9月28日夜，滑坡發(fā)生，所幸無人員傷亡及設(shè)備損失。

圖4 1-B監(jiān)測點7.2 m測斜儀累積位移-時間曲線Fig.4 The cumulative displacement-time curve of 1-B monitoring at 7.2 m by the inclinometer

隨著1 768 m臺階的局部滑坡破壞的發(fā)生，該區(qū)域應力釋放，發(fā)生應力重新調(diào)整和應力轉(zhuǎn)移，為研究卸荷作用導致的局部破壞對鄰近區(qū)域的影響程度，筆者選取了同一臺階的監(jiān)測點1-C及與滑坡部位在同一剖面上的監(jiān)測點2-C的監(jiān)測數(shù)據(jù)進行分析。

圖5是滑坡區(qū)域附近的1-C監(jiān)測點位移累積曲線，分析曲線變化趨勢可知，該監(jiān)測點3.2、7.2 m深處累積位移隨卸荷過程均保持增長，但漲幅較小，2監(jiān)測點出現(xiàn)了位移突變，3.2 m深處測斜儀監(jiān)測到的突變值為4.94 mm，7.2 m深處的測斜儀監(jiān)測到的突變值為3.07 mm，漲幅均較小。之后3.2 m和 7.2 m深度測點處的水平位移雖持續(xù)增加，但變化速率基本穩(wěn)定。說明滑坡區(qū)域的應力轉(zhuǎn)移對該點影響較小。

圖5 1-C監(jiān)測點累積位移-時間曲線Fig.5 Cumulative displacement-time curves of 1-C monitoring◆—3.2 m測斜儀；■—7.2 m測斜儀；▲—11.2 m測斜儀

圖6是滑坡區(qū)域上部剖面2-C監(jiān)測點位移累積值，分析各曲線位移值可知，2-C監(jiān)測點3個測斜儀監(jiān)測到的位移值均有較大增幅，其中3.2m處的位移值增幅最大，且在2013-05-11日之后位移增加速度加快，明顯高于7.2 m和11.2 m兩測斜儀，但是到11月份之后增速明顯放緩，7.2 m和11.2 m兩測斜儀位移增速基本保持穩(wěn)定。位移增幅也比其他監(jiān)測點大，最大時達到了101.77 mm。

由各位移累積曲線可知，邊坡的變形方向基本一致，總體趨勢都是沿垂直邊坡的方向發(fā)生變形，累計位移逐步增大。不同之處在于局部位移曲線受開采的爆破施工擾動作用出現(xiàn)震蕩性變化。在雨季開始前及結(jié)束后邊坡位移變化趨于相對平穩(wěn)狀態(tài)，而雨季監(jiān)測位移出現(xiàn)明顯增大趨勢，由于近年來礦區(qū)雨季降雨偏大，因此加強雨季的邊坡穩(wěn)定性監(jiān)測尤為重要。

圖6 2-C監(jiān)測點累積位移-時間曲線Fig.6 Cumulative displacement-time curves of 2-C monitoring◆—3.2 m測斜儀；■—7.2 m測斜儀；▲—11.2 m測斜儀

3.2 邊坡變形機制分析

通過分析圖4中位移累積變化曲線可知，受采動卸荷作用及現(xiàn)場環(huán)境及工程因素的反復影響作用，巖體內(nèi)發(fā)生了應力重分布，在拉剪應力作用下，坡體出現(xiàn)應力集中，導致另一些方向的加荷作用，從而造成巖體的損傷破壞，邊坡巖體損傷逐步累積，坡體的原生構(gòu)造節(jié)理、裂隙在集中應力作用下開始擴展，2013-05-28日前坡體水平位移保持較小增幅，隨后發(fā)生裂隙擴展狀態(tài)的分叉。只有某些特定位置的節(jié)理和裂隙的失穩(wěn)擴展會導致坡體破壞，其他節(jié)理裂隙對坡體破壞影響較小，形成局部區(qū)域的損傷變形集中現(xiàn)象。隨著應力集中的繼續(xù)，局部化損傷帶相互影響，張拉應力增大，導致巖體破壞加劇，坡體破壞加速。此過程中巖體由均勻變形逐漸發(fā)展到局部化變形直至破裂。在2013-05-29～2013-09-01日期間水平位移量迅速增加，張拉裂縫增多，分布范圍增大，裂縫長度不斷增大，寬度和深度增大，且在地表相互連接，形成臺階后緣的弧形拉裂縫，伴有下座變形，當變形達到一定程度時，坡體滑動面基本貫通，此時坡體狀態(tài)將不再穩(wěn)定，2013-09-01日后坡體水平位移急速增加，繼而發(fā)生滑坡破壞。因此，巖體結(jié)構(gòu)的破壞過程是一個伴隨巖石微破裂特征的漸進破壞誘致突變的過程。

綜合對比分析曲線圖5、圖6的數(shù)據(jù)變化趨勢可知，1 768 m臺階1-B點處的局部化破壞導致其集中應力的釋放和向周圍區(qū)域的轉(zhuǎn)移，受采動影響的破壞區(qū)域的鄰近區(qū)域都受到了不同程度的影響，通過對比所選監(jiān)測點位移累積曲線可知，變形局部化造成的應力轉(zhuǎn)移具有一定的規(guī)律性。

首先通過對比破壞區(qū)域上部的2-C監(jiān)測點和其附近的1-C監(jiān)測點可知，2-C監(jiān)測點中各處的監(jiān)測數(shù)據(jù)的水平位移累積值均比2-C監(jiān)測點的水平位移累積值大，其中3.2 m處的位移累積數(shù)據(jù)差值最大時達到84.6 mm，而7.2 m處的位移累積數(shù)據(jù)差值最大時達到26 mm。2-C監(jiān)測點臺階表面產(chǎn)生了較大的拉剪裂縫。由于較大的集中應力作用，其臺階內(nèi)部的節(jié)理及裂隙發(fā)生了較大的損傷擴展，導致坡體變形加劇。

滑坡區(qū)域下部是1 756 m臺階，由于當時監(jiān)測系統(tǒng)布設(shè)時該臺階尚未采出，未在該臺階布置監(jiān)測點，但從現(xiàn)場可以看出，由于卸荷引起的應力集中作用時間較短，臺階靠近邊緣表面有裂縫產(chǎn)生，但是裂縫寬度較小，最大處不足5 mm，說明變形局部化作用造成的應力集中對該臺階影響較小。

通過觀察各監(jiān)測點位移累積曲線可以發(fā)現(xiàn)，各點位移在增加的同時都出現(xiàn)震蕩性變化，原因在于監(jiān)測過程中的現(xiàn)場的爆破施工、自然及邊坡巖體內(nèi)部因素所致。1-C和2-C監(jiān)測點中各測斜儀的累積位移值從整體上分析，具有隨距離破壞區(qū)域越遠位移值逐漸減小的趨勢，但是根據(jù)各個巖層特性，這個趨勢也不完全一致，應力重新分布對坡體的影響還與坡體巖體的完整性有關(guān)。露天邊坡災害的發(fā)生是集中在采動影響的某些部位。巖石或巖體結(jié)構(gòu)變形的這種局部化現(xiàn)象，與巖體中的異構(gòu)密切相關(guān)。因此，注意捕捉變形局部化信息并掌握其發(fā)展趨勢，就有可能提高邊坡災害監(jiān)測預報的可能性。

巖體結(jié)構(gòu)的破壞過程是一個伴隨巖石微破裂特征的漸進破壞誘致突變的過程。采動作用導致的卸荷效應造成了局域性的應力集中使邊坡巖體原生節(jié)理、裂隙損傷擴展，產(chǎn)生局部變形破壞，并伴隨著集中應力的釋放、轉(zhuǎn)移和重新調(diào)整，在破壞區(qū)域的周圍都會受到影響，通過監(jiān)測數(shù)據(jù)和數(shù)值模擬分析可知，越靠近破壞區(qū)域受影響程度越大，集中應力釋放之后主要是向上部巖體轉(zhuǎn)移，同一水平臺階邊坡處和下部臺階邊坡受影響較小，但是對臺階邊坡巖體的破壞程度還與巖體性質(zhì)有關(guān)，巖體完整性好，原生節(jié)理、裂隙較少的巖層造成破壞就小，越不容易發(fā)生滑坡。通過監(jiān)測累積位移分析可知，集中應力造成的損傷破壞主要發(fā)生在各臺階邊坡地下7～11 m左右位置，因此7～11 m應是滑移帶形成區(qū)域，是監(jiān)測的重點。監(jiān)測過程中位移累積值由于爆破作業(yè)及自然因素的作用出現(xiàn)了震蕩性變化，這也是影響邊坡穩(wěn)定性的因素之一。

4 結(jié) 論

(1)選取了1 768 m臺階邊坡局部變形垮塌部位為研究對象，分析了其周圍1-C臺階和與滑坡體處于剖面相鄰的2-C臺階邊坡作為研究對象，對監(jiān)測過程中的位移累積值進行了分析。結(jié)果表明采動卸荷造成導致局域性的應力應變集中，使臺階邊坡出現(xiàn)變形局部化現(xiàn)象，發(fā)生應力釋放、轉(zhuǎn)移和重新調(diào)整。分析結(jié)果可知，采動卸荷應力主要向上部臺階邊坡巖體轉(zhuǎn)移，下部及同一臺階邊坡變形較小。

(2)利用RFPA強度折減法，再現(xiàn)了1 768 m平臺開采過程中典型剖面多臺階漸進失穩(wěn)過程，隨著巖體力學性質(zhì)的不斷劣化，裂隙首先從1 768 m平臺坡腳處開始萌生，不斷擴展、貫通，宏觀的裂縫在這個細觀上“破裂—應力轉(zhuǎn)移—破裂”的循環(huán)中不斷發(fā)育，最終形成貫通的滑動面。破壞的整個過程中伴隨著應力集中、釋放、轉(zhuǎn)移造成的應變局部化現(xiàn)象。

(3)通過數(shù)據(jù)分析可知，在各臺階邊坡下部7～11 m位置處受轉(zhuǎn)移荷載影響程度較大，水平位移累積值加速增大，另外受到爆破開采及自然因素的影響作用，水平位移累積值出現(xiàn)了震蕩性變化，也是影響邊坡穩(wěn)定的因素之一。雨季的坡體累積位移值明顯大于平時，需加強監(jiān)測力度。

[1] 張永興.邊坡工程學[M].北京：中國建筑工業(yè)出版社，2008. Zhang Yongxing.Slope Engineering[M].Beijing：China Architecture and Building Press，2008.

[2] 楊光華,鐘志輝,張玉成,等.根據(jù)應力場和位移場判斷滑坡的破壞類型及最優(yōu)加固位置確定[J].巖石力學與工程學報，2012，31(9):1880-1887. Yang Guanghua，Zhong Zhihui，Zhang Yucheng，et al.Identification of landslide type and determination of optimal reinforcement site based on stress field and displaeneat field[J].Chinese Journal of Rock Mechanics and Engineering，2012,31(9):1880-1887.

[3] 尹光志,張衛(wèi)中,張東明,等.基于指數(shù)平滑法與回歸分析相結(jié)合的滑坡預測[J].巖土力學,2007,28(8):1725-1728. Yin Guangzhi，Zhang Weizhong，Zhang Dongming，et al.Forecasting of landslide displacement based on exponential smoothing and nonlinear regression analysis[J].Rock and Soil Mechanics,2007,28(8):1725-1728.

[4] 鄭穎人,趙尚毅.有限元強度折減法在土坡與巖坡中的應用[J].巖石力學與工程學報，2004，23(19):3381-3388. Zheng Yingren，Zhao Shangyi.Application of strength reduction FEM to soil and rock slope[J].Chinese Journal of Rock Mechanics and Engineering,2004,23(19):3381-3388.

[5] 唐春安,李連崇,李常文,等.巖土工程穩(wěn)定性分析RFPA強度折減法[J].巖石力學與工程學報,2006,25(8):1522-1530. Tang Chun'an，Li Lianchong，Li Changwen，et al.RFPA strength reduction method for stability analysis of geotechnical engineering[J].Chinese Journal of Rock Mechanics and Engineering，2006,25(8):1522-1530.

[6] 陳國慶,黃潤秋,石豫川,等.基于動態(tài)和整體強度折減法的邊坡穩(wěn)定性分析[J].巖石力學與工程學報，2014，33(2):243-256. Chen Guoqing，Huang Runqiu，Shi Yuchuan，et al.Stability analysis of slope based on dynamic and whole strength reduction methods[J].Chinese Journal of Rock Mechanics and Engineering，2014，33(2):243-256.

[7] 許 強,湯明高,徐開祥,等.滑坡時空演化規(guī)律及預警預報研究[J].巖石力學與工程學報，2008，27(6):1104-1112. Xu Qiang，Tang Minggao，Xu Kaixiang，et al.Research on space-time evolution laws and early warning-prediction of landslides[J].Chinese Journal of Rock Mechanics and Engineering，2008，27(6):1104-1112.

[8] 何滿潮.滑坡地質(zhì)災害遠程監(jiān)測預報系統(tǒng)及其工程應用[J].巖石力學與工程學報，2009，28(6):1081-1090. He Manchao.Real-time remote monitoring and forecasting system for geological disasters of landsides and its engineering application[J].Chinese Journal of Rock Mechanics and Engineering，2009，28(6):1081-1090.

[9] 張振華,馮夏庭,周 輝,等.基于設(shè)計安全系數(shù)及破壞模式的邊坡開挖過程動態(tài)變形監(jiān)測預警方法研究[J].巖土力學，2009，30(3):603-612. Zhang Zhenhua，F(xiàn)eng Xiating，Zhou Hui，et al.Research on dynamic early warning method of slope deformation monitoring during excavation based on designed safety factor and failure mode[J].Rock and Soil Mechanics，2009，30(3):603-612.

[10] 王 東,曹蘭柱,樸春德,等.基于假設(shè)檢驗原理的邊坡臨滑時刻的動態(tài)識別方法[J].巖石力學與工程學報，2012，31(3):577-585. Wang Dong，Cao Lanzhu，Piao Chunde，et al.Dynamic identification method of slope impending landslide moment on hypothetic testing theory[J].Chinese Journal of Rock Mechanics and Engineering，2012，31(3):577-585.

[11] 王庚蓀.邊坡的漸進破壞及穩(wěn)定性分析[J].巖石力學與工程學報,2000,19(1):29-33. Wang Gengsun.The progressive failure of slope and the stability analysis[J].Chinese Journal of Rock Mechanics and Engineering，2000，19(1):29-33.

[12] 李連崇,唐春安,梁正召,等.RFPA 邊坡穩(wěn)定性分析方法及其應用[J].應用基礎(chǔ)與工程科學學報，2007,15(4):501-508. Li Lianchong，Tang Chun'an，Liang Zhengzhao，et al.RFPA method for slope stability analysis and associated application[J].Journal of Basic Science and Engineering，2007，15(4):501-508.

(責任編輯 石海林)

Numerical Analysis and Monitoring of Local Instability Characteristics of the Open Pit Slope

Zhang Yongjun1Yang Dengfeng2Wang Jintao1Chen Guanping1

(1.SchoolofCivilEngineering，QingdaoTechnologicalUniversity，Qingdao266033，China;2.SchoolofMechanics&CivilEngineering，ChinaUniversityofMining&TechnologyBeijing，Beijing100083，China)

Based on slope monitoring project of Open Pit Quartz Mine of Jinchuan Group as background，with the use of the monitoring data，and taking local landslide of the beach as the research object，the deformation and failure mechanism of the slope under the stress concentration action caused by unloading effect is studied through analyze the cumulative horizontal displacement values of the monitoring sites around the beach and on the landslide.With RFPA strength subtraction method，the gradual instability process of typical profile and multi-steps during the mining process of 1 768 mm platform is reproduced.macro-cracks on the meso cycle of “rupture-stress transfer-rupture” is developed constantly，eventually forming a cut-through sliding surface.The strain localization phenomenon caused by stress concentration，release，and transfer process is accompanied by the whole process of failure.The monitoring results show that:the process of localized damage is migrated from mining step to the upper bench slope rock mass，and 7～11 m point of each beach slope is greatly affected by load transfer.

Slope，Strain localization，Strength reduction method，Monitoring，Data analysis

2014-07-28

國家自然科學基金重點項目(編號：51234005)。

張擁軍(1974—)，男，副教授，博士后。

TD854.6

A

1001-1250(2014)-10-142-06