任富強(qiáng) 常 遠(yuǎn) 常來山 陳東偉

(1. 中國(guó)礦業(yè)大學(xué)（北京）力學(xué)與建筑工程學(xué)院，北京100083;2. 深部巖土力學(xué)與地下工程國(guó)家重點(diǎn)實(shí)驗(yàn)室，北京100083;3. 遼寧科技大學(xué)礦業(yè)工程學(xué)院，遼寧 鞍山114051)

以往對(duì)邊坡的可靠性分析中，巖體一般是作為均質(zhì)的隨機(jī)介質(zhì)考慮，認(rèn)為不同點(diǎn)間的巖體強(qiáng)度參數(shù)是一致的，這顯然是不符合實(shí)際的，實(shí)際的巖土介質(zhì)不同點(diǎn)間、整體和局部的強(qiáng)度參數(shù)都有一定的差異性，這就需要考慮巖體的空間變異性，其中Vanmarcke［1-2］在該理論方面進(jìn)行了開創(chuàng)性研究，隨后Chowdhurry［3］利用隨機(jī)場(chǎng)理論描述了巖體性質(zhì)的空間變異性。地質(zhì)統(tǒng)計(jì)學(xué)理論被廣泛應(yīng)用于巖體空間變異性的評(píng)估，其中王家臣等［4-7］用地質(zhì)統(tǒng)計(jì)學(xué)理論中的變差函數(shù)經(jīng)過推導(dǎo)得出巖體性質(zhì)參數(shù)的相關(guān)距離和相關(guān)函數(shù)，說明了巖體性質(zhì)的空間變異性，并進(jìn)行了考慮巖體空間變異性的邊坡可靠性分析，但只考慮了1 種理論變差函數(shù)模型。

本研究結(jié)合野外實(shí)測(cè)點(diǎn)荷載數(shù)據(jù)擬合出了理論變差函數(shù)的球狀模型、高斯模型和指數(shù)模型的相應(yīng)曲線，并利用蒙特卡洛抽樣方法結(jié)合理論變差函數(shù)模型對(duì)單一巖性的邊坡巖體進(jìn)行了可靠性分析。

1 巖體強(qiáng)度的空間變異性

1.1 現(xiàn)場(chǎng)點(diǎn)荷載測(cè)試

為了得到巖體強(qiáng)度的空間分布特征，在某露天鐵礦的某一勘探剖面上垂直巖層走向方向以間距1 m等間距布置46 個(gè)測(cè)點(diǎn)。在每個(gè)測(cè)點(diǎn)上利用SD －1型數(shù)碼點(diǎn)荷載儀對(duì)選取的若干混合巖巖塊進(jìn)行了點(diǎn)荷載實(shí)驗(yàn)，將實(shí)測(cè)結(jié)果按ISRM 法并經(jīng)過尺寸修正計(jì)算得到標(biāo)準(zhǔn)點(diǎn)荷載強(qiáng)度指數(shù)Is(50)，其垂直于巖層走向分布情況如圖1 所示。由圖1 可知，在空間位置上，標(biāo)準(zhǔn)點(diǎn)荷載強(qiáng)度指數(shù)Is(50)有較大的差異，但其整體較為穩(wěn)定，未曾發(fā)生漂移。對(duì)實(shí)測(cè)的數(shù)據(jù)進(jìn)行統(tǒng)計(jì)分析后得出其平均值為11.60 MPa，標(biāo)準(zhǔn)差為6.79 MPa，變異系數(shù)為0.585。

圖1 混合巖Is(50)垂直巖層走向分布情況Fig.1 The vertical strata distribution of migmatite Is(50)

1.2 實(shí)驗(yàn)變差函數(shù)與理論變差函數(shù)

混合巖Is(50)為區(qū)域化的變量且滿足二階平穩(wěn)假設(shè)，故可按地質(zhì)統(tǒng)計(jì)學(xué)中的實(shí)驗(yàn)變差函數(shù)［8-10］來計(jì)算Is(50)的實(shí)驗(yàn)變差函數(shù)值r(h)，即

式中，h 為取樣點(diǎn)的間隔距離，取1 m;N(h)為被h 相隔的數(shù)據(jù)對(duì)的數(shù)目;Z(xi)和Z(xi+h)在此為Is(50)的不同取值。

為了預(yù)估區(qū)域化的變量Is(50)的未知值，需要根據(jù)實(shí)驗(yàn)變差函數(shù)值擬合成相應(yīng)的理論變差函數(shù)模型，其中球狀模型、高斯模型和指數(shù)模型的理論公式分列如下。

(1)球狀模型:

式中，C0為塊金常數(shù);C0+ C 為基臺(tái)值;C 為拱高;a為變程(即隨機(jī)場(chǎng)中的相關(guān)距離D)。

(2)高斯模型:

(3)指數(shù)模型:

式中，a 不是變程，變程為3a，即h =3a 時(shí)，r(h)≈C0+C。

采用上述理論變差函數(shù)模型中的球狀模型、高斯模型和指數(shù)模型，根據(jù)實(shí)驗(yàn)變差函數(shù)曲線利用Visual Basic 進(jìn)行最小二乘擬合，擬合后3 個(gè)理論變差函數(shù)模型的參數(shù)見表1。由表1 可知指數(shù)模型相關(guān)距離最大，高斯模型的最小，由相關(guān)距離的物理含義可知同等長(zhǎng)度內(nèi)指數(shù)模型的數(shù)據(jù)點(diǎn)的相關(guān)性最大，高斯模型最小;此外指數(shù)模型與高斯模型的基臺(tái)值相同，略大于球狀模型的基臺(tái)值，可近似認(rèn)為3 種模型描述的Is(50)在垂直于巖層走向方向上變化幅度大小相同。實(shí)驗(yàn)變差函數(shù)曲線與理論變差函數(shù)曲線如圖2 所示，由圖2 可知3 種模型的擬合曲線的趨勢(shì)一致，且對(duì)應(yīng)的擬合值基本相同。

表1 混合巖Is(50)計(jì)算結(jié)果Table 1 The calculation result of migmatite Is(50)

圖2 混合巖3 種模型實(shí)驗(yàn)與理論變差函數(shù)曲線Fig.2 Curves of three kinds of model experiment and theoretic variation function for migmatite

1.3 巖體強(qiáng)度局部分析

由理論變差函數(shù)，可以求出巖體強(qiáng)度的相關(guān)函數(shù)。根據(jù)巖體強(qiáng)度，相關(guān)函數(shù)的定義為

式中，R(h)為強(qiáng)度相關(guān)函數(shù);r(h)為理論變差函數(shù);Var(Z(x))為Is(50)的方差，混合巖的方差為46.043。

根據(jù)式(5)可以得到混合巖3 種模型對(duì)應(yīng)的R(h)，而在邊坡可靠性計(jì)算中一般要將巖體分條再將每一分條作為均質(zhì)體處理，在每一分條上巖體強(qiáng)度要局部平均化，在二階平穩(wěn)假設(shè)的前提下局部平均化后的期望與整體的期望相同，但方差要發(fā)生變化。根據(jù)局部平均域大小、強(qiáng)度的空間相關(guān)距離、相關(guān)函數(shù)等不同，其方差衰減的程度不同，一般用式(6)表示巖體強(qiáng)度局部平均化后的方差，式中ZT(x)表示平均域長(zhǎng)度為T 的巖體強(qiáng)度函數(shù)，G(T)為方差衰減函數(shù)，可由式(7)計(jì)算。故可以得到混合巖3 種模型對(duì)應(yīng)的方差衰減函數(shù)以及方差函數(shù)。

2 邊坡可靠性分析

為方便對(duì)比分析，探討不同的變差函數(shù)模型對(duì)邊坡可靠性的影響規(guī)律，構(gòu)造了圖3 所示的均質(zhì)混合巖的簡(jiǎn)單邊坡模型，邊坡高度為400 m，邊坡角為45°。

圖3 混合巖簡(jiǎn)單邊坡模型Fig.3 The simple slope model of migmatite

應(yīng)用Monte－Calo 隨機(jī)抽樣技術(shù)［11-12］，通過Visual Basic 編程得到巖體強(qiáng)度參數(shù)的隨機(jī)分布特征，模型的計(jì)算參數(shù)如表2 所示。極限平衡計(jì)算選用畢肖普(BISHOP)法，用單純形變法［13］進(jìn)行優(yōu)化計(jì)算，確定了危險(xiǎn)滑弧及安全系數(shù)，危險(xiǎn)滑弧如圖3 曲線所示，混合巖邊坡的安全系數(shù)F 為1.835。

表2 巖體強(qiáng)度參數(shù)隨機(jī)分布表Table 2 The random distribution of rock strength parameters

根據(jù)不同的分條寬度將球狀模型、高斯模型、指數(shù)模型的可靠性指標(biāo)［14］的部分計(jì)算結(jié)果列于表3 所示。

表3 混合巖邊坡Monte-Calo 法部分計(jì)算結(jié)果Table 3 Parts of calculation result by Monte-Calo method for migmatite slope

由表3 比較混合巖球狀模型、高斯模型、指數(shù)模型的邊坡可靠性指標(biāo)計(jì)算結(jié)果。若以球狀模型為基準(zhǔn)，高斯模型的絕對(duì)差值為0.003 ～0.312，平均為0.053，相對(duì)差值為0.515% ～17.154%，平均為3.383%;指數(shù)模型的絕對(duì)差值為0.046 ～2.568，平均為1.141，相對(duì)差值為11.348% ～83.229%，平均為59.941%;指數(shù)模型相差的較多，可靠性計(jì)算結(jié)果偏高，而高斯模型與球狀模型的可靠性計(jì)算結(jié)果較為接近。將3 種模型的可靠指標(biāo)隨分條寬度變化的關(guān)系曲線繪制如圖4 所示。

圖4 3 種模型可靠性指標(biāo)對(duì)比Fig.4 The reliability index contrast of three kinds of model

由圖4 可知，在分條寬度T =0 即不考慮巖體強(qiáng)度空間變異性時(shí)與T 趨于0 時(shí)有一個(gè)落差，并且差值的大小與塊金常數(shù)C0有關(guān)。不考慮巖體強(qiáng)度的空間變異性時(shí)邊坡的可靠指標(biāo)最小，隨分條寬度的增加可靠指標(biāo)整體呈增大的趨勢(shì)，其中球狀模型與球狀模型趨勢(shì)基本一致都低于指數(shù)模型的增長(zhǎng)趨勢(shì)。當(dāng)分條寬度大于相關(guān)距離時(shí)3 種模型的可靠指標(biāo)均有下降的趨勢(shì)直到達(dá)到一定的分條寬度時(shí)又逐漸隨分條寬度的增大而上升。

3 結(jié) 論

(1)高斯模型的相關(guān)距離最小，指數(shù)模型的最大，3 種模型描述的巖體強(qiáng)度沿垂直于巖層走向方向的變化幅度基本接近。

(2)當(dāng)考慮巖體強(qiáng)度空間變異性時(shí)計(jì)算邊坡的可靠指標(biāo)增加，破壞概率降低。主要是由于強(qiáng)度參數(shù)方差衰減所致，而且在分條寬度達(dá)到相關(guān)距離前隨分條寬度的增加而增大，分條寬度達(dá)到相關(guān)距離之后的可靠指標(biāo)先減小后增大。

(3)球狀模型與高斯模型的可靠指標(biāo)隨分條寬度的變化基本一致，整體值都低于指數(shù)模型的可靠指標(biāo)，增加的趨勢(shì)也較指數(shù)模型平緩。

［1］ Vanmarke E H. Random Fields:Analysis and Synthesis［M］.［S.l.］:The MIT Press，1983.

［2］ Vanmarke E H. Probabilistic stability analysis of earth slopes［J］.Eng Geology，1980，16(1):29-50.

［3］ Chowdhurry R N.Simulation of risk of progressive slope failure［J］.Can Geotech J，1992，29(1):94-102.

［4］ 譚文輝，王家臣，周汝弟.巖體強(qiáng)度參數(shù)空間變異性分析［J］.巖石力學(xué)與工程學(xué)報(bào)，1999，18(5):497-502.

Tan Wenhui，Wang Jiachen，Zhou Rudi.Analysis on spatial variability of strength parameter of rock mass［J］. Chinese Journal of Rock Mechanics and Engineering，1999，18(5):497-502.

［5］ 譚文輝，蔡美峰，周汝弟.巖體空間變異性對(duì)邊坡可靠性的影響研究［J］.工程力學(xué)，2001(S):517-521.

Tan Wenhui，Cai Meifeng，Zhou Rudi.Influence of spatial variability of rock mass on the slope reliability research［J］. Engineering Mechanics，2001(S):517-521.

［6］ 王家臣，王艷青，寧 柯.應(yīng)用地質(zhì)統(tǒng)計(jì)學(xué)研究巖體強(qiáng)度的空間變異性［J］.黃金科學(xué)技術(shù)，1999，7(45):51-54.

Wang Jiacheng，Wang Yanqing，Ning Ke.Application of geostatistics study spatial variation of rock mass strength［J］. Gold Science and Technology，1999，7(45):51-54.

［7］ 王家臣，譚文輝，寧 柯.義馬北露天煤礦工作幫邊坡可靠性研究［J］.巖土工程學(xué)報(bào)，2000，22(1):105-108.

Wang Jiachen，Tan Wenhui，Ning Ke.The reliability research on the working slope of North Surface Coal Mine of Yima Coal Administration［J］. Chinese Journal of Geotechnical Engineering，2000，22(1):105-108.

［8］ 王仁鐸，胡光道.線性地質(zhì)統(tǒng)計(jì)學(xué)［M］.北京:地質(zhì)出版社，1988.

Wang Renduo，Hu Guangdao.Linear Geostatistics［M］.Beijing:Geological Press，1988.

［9］ 張仁鐸.空間變異理論及應(yīng)用［M］.北京:科學(xué)出版社，2006.

Zhang Renduo. Spatial Variation Theory and Application［M］. Beijing:Science Press，2006.

［10］ 李章林，王 平，李冬梅.實(shí)驗(yàn)變差函數(shù)計(jì)算方法的研究與運(yùn)用［J］.國(guó)土資源信息化，2008，8(2):10-14.

Li Zhanglin，Wang Ping，Li Dongmei. Research and application of the method calculating experiment variogram［J］. Land and Resources Informatization，2008，8(2):10-14.

［11］ 李育超，凌道盛，陳云敏，等.蒙特卡洛法與有限元相結(jié)合分析邊坡穩(wěn)定性［J］. 巖石力學(xué)與工程學(xué)報(bào)，2005，24(11):1933-1941.

Li Yuchao，Ling Daosheng，Chen Yunmin，et al. Slope stability analysis using monte carlo technique with FEM［J］.Chinese Journal of Rock Mechanics and Engineering，2005，24(11):1933-1941.

［12］ 李 亮，褚雪松，袁長(zhǎng)豐，等.快速蒙特卡羅法及其在土坡可靠度分析中的應(yīng)用［J］.煤炭學(xué)報(bào)，2013，38(9):1576-1582.

Li Liang，Chu Xuesong，Yuan Changfeng，et al. Fast Monte-Carlo method and its application in the soil slope reliability analysis［J］.Journal of China Coal Society，2013，38(9):1576-1582.

［13］ 燕子宗，費(fèi)浦生，萬(wàn)仲平. 線性規(guī)劃的單純形法及其發(fā)展［J］.計(jì)算數(shù)學(xué)，2007，29(1):121-127.

Yan Zizong，F(xiàn)ei Pusheng，Wan Zhongping.Classical simplex methods for linear programming and their develppments［J］.Mathematica Numerica Sinica，2007，29(1):121-127.

［14］ 張昀青，李維珍，閆靜昌，等.邊坡可靠性指標(biāo)及敏感性因素分析［J］.石家莊鐵道大學(xué)學(xué)報(bào):自然科學(xué)版，2015，28(1):62-65.

Zhang Yunqing，Li Weizhen，Yan Jingchang，et al.The slope reliability index and sensitivity factor analysis［J］.Journal of Shijiazhuang Railway Institute:Natural Science，2015，28(1):62-65.