王　帆, 焦振華

(陜西省水利電力勘測(cè)設(shè)計(jì)研究院,陜西 咸陽(yáng)　712000)

0　引言

目前，中國(guó)水利水電建設(shè)事業(yè)處在一個(gè)高潮時(shí)期，水利水電工程邊坡、水庫(kù)岸坡、河道岸坡等滑坡事故時(shí)有發(fā)生?；虏粌H造成局部環(huán)境破壞，而且還危害人民生命財(cái)產(chǎn)安全及經(jīng)濟(jì)建設(shè)的發(fā)展。

滑坡穩(wěn)定性受各種內(nèi)在因素及外在誘發(fā)因素的影響，如滑坡所處地形地貌、地層巖性、降雨作用、地表水和地下水，以及地震作用等[1]。為了有效地防治必須了解各影響因素對(duì)滑坡穩(wěn)定性的影響程度，通過(guò)滑坡穩(wěn)定性因素敏感性分析找出敏感性最大的因素，對(duì)于滑坡治理及優(yōu)化設(shè)計(jì)有著相當(dāng)大的意義，使滑坡治理能夠達(dá)到安全、經(jīng)濟(jì)和有效的目的。

1　敏感性分析

1.1　敏感性分析及敏感性系數(shù)定義

敏感性分析就是指從定量分析的角度研究有關(guān)因素發(fā)生某種變化對(duì)某一個(gè)或一組關(guān)鍵指標(biāo)影響程度的一種不確定分析技術(shù)。敏感性系數(shù)定義為各屬性對(duì)模型影響程度的大小。若敏感性系數(shù)越大，則這個(gè)屬性對(duì)模型的輸出結(jié)果影響越大。

獲得各屬性敏感性系數(shù)的大小是敏感性分析最為重要的目的。為了提高模型的精度、降低模型復(fù)雜度，同時(shí)減少數(shù)據(jù)處理量，應(yīng)該去除敏感性系數(shù)較小的屬性，主要考慮敏感性系數(shù)較大屬性。人們可根據(jù)各屬性敏感性系數(shù)的大小排序，解決相應(yīng)的問(wèn)題[2]。

1.2　敏感性分析法

根據(jù)不確定性因素每次變動(dòng)數(shù)目的多少，敏感性分析法可以分為局部分析和全局分析。

(1)局部分析只檢驗(yàn)單個(gè)屬性對(duì)模型的影響程度，該方法能夠比較直觀地反映各因素對(duì)基準(zhǔn)值的影響，但是需要假設(shè)前提條件，與實(shí)際情況不符，因而存在一定局限性。

(2)全局分析可以檢驗(yàn)多個(gè)屬性對(duì)模型結(jié)果產(chǎn)生的總體影響,同時(shí)分析屬性之間的相互作用對(duì)模型輸出結(jié)果的影響，此法是在單因素敏感性分析基礎(chǔ)上進(jìn)行，與單因素分析原理大體相同，但多因素分析須進(jìn)一步假定同時(shí)變動(dòng)的幾個(gè)因素都是相互獨(dú)立的，而且各因素發(fā)生變化的概率相同。

由于現(xiàn)有主要的敏感性分析方法采用的分析工具不同，筆者把它們分為兩種類型:基于數(shù)理統(tǒng)計(jì)的和基于神經(jīng)網(wǎng)絡(luò)的。

(1)數(shù)理統(tǒng)計(jì)方法有一定局限性，當(dāng)模型屬性太多或者得到的結(jié)果與屬性之間是一種非線性關(guān)系時(shí)，其結(jié)果不理想,而且精度達(dá)不到要求。

(2)人工神經(jīng)網(wǎng)絡(luò)方法較統(tǒng)計(jì)法有一定優(yōu)勢(shì)，分析模型時(shí)僅僅需要獲得輸入變量數(shù)據(jù)和輸出數(shù)據(jù)，其自身可以對(duì)訓(xùn)練數(shù)據(jù)集進(jìn)行訓(xùn)練和學(xué)習(xí)，運(yùn)用大量簡(jiǎn)單的人工神經(jīng)元模擬數(shù)據(jù)間的非線性關(guān)系,并且能夠自動(dòng)適應(yīng)和調(diào)節(jié)神經(jīng)元之間的連接權(quán)重，以此建立能夠較好反映數(shù)據(jù)真實(shí)情況的網(wǎng)絡(luò)結(jié)構(gòu)。

以前廣泛采用的敏感性分析的一般方法是單因素分析法，這種方法存在明顯的局限性。而對(duì)于多因素敏感性分析的研究，目前已出現(xiàn)一些方法，如采用全微分法、正交設(shè)計(jì)、均勻設(shè)計(jì)、連環(huán)代替法、模糊算法、可靠度分析、BP神經(jīng)網(wǎng)絡(luò)等方法進(jìn)行多因素分析[3-6]。

1.3　滑坡穩(wěn)定性敏感性分析

滑坡穩(wěn)定性敏感性分析，是指定量地分析影響滑坡穩(wěn)定性各因素與滑坡穩(wěn)定系數(shù)之間的相關(guān)性，研究和預(yù)測(cè)這些因素的變動(dòng)對(duì)穩(wěn)定系數(shù)的影響程度。

2　正交設(shè)計(jì)滑坡穩(wěn)定性影響因素敏感性分析方法

正交設(shè)計(jì)法可以通過(guò)分析部分試驗(yàn)的結(jié)果，了解全面試驗(yàn)的情況，體現(xiàn)其正交性、代表性以及綜合可比性；另一方面其計(jì)算過(guò)程簡(jiǎn)便，且次數(shù)較少，因而文中選用正交設(shè)計(jì)法。正交設(shè)計(jì)方法具有一般性，有其理論意義和工程實(shí)用價(jià)值。

2.1　正交設(shè)計(jì)法

正交設(shè)計(jì)方法是一種數(shù)理統(tǒng)計(jì)方法，可以解決多因素、多水平、多指標(biāo)的統(tǒng)計(jì)分析方法，它是利用已規(guī)格化的正交表來(lái)安排試驗(yàn)的一種科學(xué)方法。它是在全部水平組合內(nèi)挑選部分具有代表性的水平組合完成試驗(yàn)，并通過(guò)對(duì)這些試驗(yàn)結(jié)果的分析達(dá)到了解全面試驗(yàn)的情況，從而找出最優(yōu)的水平組合。它可以分析以下幾種結(jié)果：

(1)因素的主次關(guān)系，即各因素對(duì)所考察指標(biāo)影響大小順序。

(2)因素與指標(biāo)的關(guān)系，即每個(gè)因素水平不同時(shí)，指標(biāo)是如何變化的。

(3)得出較優(yōu)水平組合。

(4)對(duì)試驗(yàn)誤差進(jìn)行粗略估計(jì)。

正交設(shè)計(jì)的基本特點(diǎn)是：用部分試驗(yàn)代替全面試驗(yàn)，通過(guò)分析部分試驗(yàn)的結(jié)果，了解全面試驗(yàn)的情況。作為處理多因素試驗(yàn)的方法，只要較少次數(shù)的試驗(yàn)就可以從眾多的全面試驗(yàn)中挑選出具有代表性的組合條件，進(jìn)一步分析試驗(yàn)結(jié)果又能尋找出各因素作用的影響程度。因此它是進(jìn)行多因素敏感性分析的有力工具。

設(shè)A、B、C……表示不同的因素，r為各因素的水平，Ai表示因素A的第i水平(i=1,2……r)，Xij表示因素j的第i水平的值(i=1,2……r；j=A,B……)。在Xij下進(jìn)行試驗(yàn)得到因素j第i水平的試驗(yàn)結(jié)果指標(biāo)Yij，Yij是服從正態(tài)分布的隨機(jī)變量。在Xij下做了n次試驗(yàn)得到n個(gè)試驗(yàn)結(jié)果，分別為Yijk(k=1,2,……,n)。有計(jì)算參數(shù)如下：

(1)

式中：Kij為因素j在i水平下的統(tǒng)計(jì)參數(shù)；n為因素j在i水平下的試驗(yàn)次數(shù)；Yijk為因素j在i水平下第k個(gè)試驗(yàn)結(jié)果指標(biāo)值。

評(píng)價(jià)因素顯著性的參數(shù)為極差Rj，其計(jì)算公式如下：

Rj=max(K1j,K2j,…Krj)-min(K1j,K2j,…Krj)

(2)

極差大小代表該因素的水平改變對(duì)試驗(yàn)結(jié)果影響程度。極差最大的因素也就是最主要的因素，極差較小的因素為較次要的因素[7-8]。

2.2　正交設(shè)計(jì)基本原理

正交設(shè)計(jì)就是從全面試驗(yàn)點(diǎn)(水平組合)中挑選出有代表性的部分試驗(yàn)點(diǎn)(水平組合)來(lái)進(jìn)行試驗(yàn)。

比如3個(gè)因素的全面試驗(yàn)可以用一個(gè)立方體表示，如圖1，3個(gè)因素分別取3個(gè)水平，可把立方體分成成27個(gè)格點(diǎn)，在圖1上就是立方體中27個(gè)點(diǎn)。若27個(gè)網(wǎng)格點(diǎn)都試驗(yàn)，就是全面試驗(yàn)。

圖1中標(biāo)有試驗(yàn)號(hào)的9個(gè)點(diǎn)，就是利用正交表L9(34)從27個(gè)試驗(yàn)點(diǎn)中挑選出來(lái)的9個(gè)試驗(yàn)點(diǎn)。9個(gè)試驗(yàn)點(diǎn)在選優(yōu)區(qū)中分布是均衡的，在立方體的每個(gè)平面上，都恰有3個(gè)試驗(yàn)點(diǎn)；在立方體的每條線上也恰有一個(gè)試驗(yàn)點(diǎn)。9個(gè)試驗(yàn)點(diǎn)均衡地分布于整個(gè)立方體內(nèi)，有很強(qiáng)的代表性，能夠比較全面地反映選優(yōu)區(qū)內(nèi)的基本情況。

圖1　3因素3水平的均衡分散立體圖

2.3　正交表基本性質(zhì)

正交設(shè)計(jì)安排試驗(yàn)和分析試驗(yàn)結(jié)果都需要用正交表，因此，對(duì)正交表作簡(jiǎn)要介紹。記號(hào)為L(zhǎng)a(bc)的正交表中，“L”表示正交表；“a”指試驗(yàn)總次數(shù)表示有a行 ，括號(hào)內(nèi)“b”為因素水平數(shù)，b的指數(shù)“c”指因素個(gè)數(shù)表示有c列。常用的正交表已由數(shù)學(xué)工作者制定出來(lái)，供進(jìn)行正交設(shè)計(jì)時(shí)選用。如2水平正交表除L8(27)外，還有L4(23)、L16(215)等；3水平正交表有L9(34)、L27(213)等。

正交表有三個(gè)基本性質(zhì)，其中正交性是核心和基礎(chǔ)，代表性與綜合可比性是正交性反應(yīng)的必然結(jié)果。

(1)正交性任一列中，各水平都出現(xiàn)，且出現(xiàn)的次數(shù)相等。例如L8(27)一列中不同數(shù)字只有1和2，它們各出現(xiàn)4次；L9(34)一列中不同數(shù)字有1、2和3，它們各出現(xiàn)3次。任兩列間所有水平組合出現(xiàn)次數(shù)相等。

(2)代表性由于正交表具有正交性，正交試驗(yàn)點(diǎn)必然均勻地分布在全面試驗(yàn)點(diǎn)中，體現(xiàn)了代表性。

(3)綜合可比性因素各水平在任一列出現(xiàn)次數(shù)相等；任兩列間所有水平組合出現(xiàn)次數(shù)也相等，使得任一因素各水平的試驗(yàn)條件相同。最大限度地排除了其他因素的干擾。從而可以綜合比較該因素不同水平對(duì)試驗(yàn)指標(biāo)的影響情況。

2.4　正交設(shè)計(jì)分析過(guò)程

正交試驗(yàn)設(shè)計(jì)的基本程序包括試驗(yàn)方案設(shè)計(jì)及試驗(yàn)結(jié)果分析兩部分。

(1)試驗(yàn)方案設(shè)計(jì)根據(jù)實(shí)際情況，選擇一張適合的正交表，記為L(zhǎng)a(bc)，把各因素隨機(jī)地填入表的上方，稱之為表頭設(shè)計(jì)。

(2)按照計(jì)算方案進(jìn)行計(jì)算，記錄計(jì)算結(jié)果。

(3)試驗(yàn)結(jié)果分析分析這些結(jié)果并繪制各因素與指標(biāo)趨勢(shì)圖，得出結(jié)論。分析正交計(jì)算數(shù)據(jù)的方法有兩種，極差分析法和方差分析法。因前者簡(jiǎn)單又方便，故本文選用極差分析法。

3　工程實(shí)例分析

3.1　滑坡工程概況

涇河?xùn)|莊水利樞紐工程位于陜西省禮泉縣與淳化縣交界的涇河下游峽谷河段，擬建壩址右岸屬禮泉縣，距縣城40 km；左岸屬淳化縣，距縣城30 km；壩址距咸陽(yáng)市區(qū)77 km，距西安市105 km。

東莊水庫(kù)庫(kù)區(qū)21#滑坡體前緣剪出口位于蓄水位以下，蓄水后可能產(chǎn)生變形或滑移破壞。由于影響滑坡穩(wěn)定性因素較多，通過(guò)對(duì)滑坡穩(wěn)定性敏感度分析，找到滑坡失穩(wěn)的主導(dǎo)因素，可以為滑坡治理及優(yōu)化設(shè)計(jì)提出針對(duì)性的整治措施和建議。

東莊水庫(kù)庫(kù)區(qū)21#滑坡體主要組成物質(zhì)為黃土及黃土狀土，下部位砂巖及頁(yè)巖，滑坡體形態(tài)如圖2所示。

圖2　東莊庫(kù)區(qū)21#滑坡地質(zhì)剖面

3.2　正交設(shè)計(jì)方案

對(duì)于滑坡穩(wěn)定性計(jì)算，選擇考慮粘聚力、內(nèi)摩擦角、重度、地震等4個(gè)因素進(jìn)行敏感性分析，各影響因素分析均只考慮土體參數(shù)的自相關(guān)性，不考慮其相關(guān)性影響。為了排除誤差干擾，試驗(yàn)中可隨機(jī)進(jìn)行；安排試驗(yàn)方案時(shí)，部分因素的水平可隨機(jī)安排。根據(jù)參數(shù)的取值范圍選取L9(34)正交表，因素水平表列于表1。

表1　因素水平表

根據(jù)正交表設(shè)計(jì)計(jì)算方案，針對(duì)本滑坡采用Janbu法計(jì)算穩(wěn)定性系數(shù)，計(jì)算結(jié)果和極差分析列于表2。

表2　正交設(shè)計(jì)分析表

圖3　各因素與指標(biāo)趨勢(shì)圖

各因素與指標(biāo)趨勢(shì)分析如圖3所示，其橫坐標(biāo)表示因素的分類，水平值按表1中參數(shù)采取從小到大排列，縱坐標(biāo)為滑坡穩(wěn)定系數(shù)Fs平均值，即表2中統(tǒng)計(jì)參數(shù)Kij的平均值。

3.3　結(jié)果分析

東莊庫(kù)區(qū)滑坡實(shí)例通過(guò)正交設(shè)計(jì)分析有如下結(jié)論：

(1)據(jù)表2可知，滑坡穩(wěn)定性影響因素敏感性由大到小依次為：內(nèi)摩擦角(φ)、地震作用、粘聚力(c)、重度(γ)。

(2)根據(jù)圖3可知，相對(duì)內(nèi)摩擦角及地震作用，重度和粘聚力對(duì)滑坡穩(wěn)定系數(shù)影響不大。穩(wěn)定系數(shù)隨重度的增加略微減小，隨粘聚力的增加總體上呈增大趨勢(shì)，隨內(nèi)摩擦角的增加而明顯增大，對(duì)粘土坡穩(wěn)定性的影響最為顯著。隨著地震作用增大，安全系數(shù)呈現(xiàn)減小趨勢(shì)。

3.4　治理建議

(1)抗剪強(qiáng)度參數(shù)內(nèi)摩擦角φ較其他因素的敏感性強(qiáng)。影響滑帶抗剪強(qiáng)度大小的關(guān)鍵因素是巖土體飽水程度，隨著雨水的下滲，滑帶土的飽和面積越大，強(qiáng)度降低，變形破壞越嚴(yán)重，因而降雨入滲會(huì)降低滑帶土的抗剪強(qiáng)度，因此在滑坡治理設(shè)計(jì)中應(yīng)重點(diǎn)做好地表截排水系統(tǒng),滑坡體后部適當(dāng)削坡，并做護(hù)坡。

(2)抗震設(shè)計(jì)指標(biāo)是僅次于內(nèi)摩擦角的重要參數(shù)，滑坡防治應(yīng)慎重考慮地震作用的影響。

4　結(jié)論

多因素正交設(shè)計(jì)在滑坡穩(wěn)定性因素敏感性分析過(guò)程中能夠找出敏感性最強(qiáng)的因素，可使滑坡防治工作有目的性，使滑坡的治理設(shè)計(jì)安全、可靠和經(jīng)濟(jì)。

根據(jù)實(shí)例分析，采用正交分析法得到的結(jié)論與一般工程經(jīng)驗(yàn)相符。正交設(shè)計(jì)可以極大地簡(jiǎn)化試驗(yàn)數(shù)量，對(duì)單因素、多因素試驗(yàn)均能達(dá)到很好的效果。因此該方法具有一定的可靠性、通用性，并具有一定的工程實(shí)用價(jià)值。

參考文獻(xiàn)：

[1]鄭穎人，陳祖煜，王恭先，等．邊坡與滑坡工程治理[M].北京:人民交通出版社,2007.

[2]蔡毅，邢巖，胡丹，等.敏感性分析綜述[J].北京師范大學(xué)學(xué)報(bào):自然科學(xué)版， 2008(1):9-16.

[3]倪恒，劉佑榮，龍治國(guó).正交設(shè)計(jì)在滑坡敏感性分析中的應(yīng)用[J].巖石力學(xué)與工程學(xué)報(bào)，2002，21(7)：989-992.

[4]陳高峰，程圣國(guó)，盧應(yīng)發(fā)，等.基于均勻設(shè)計(jì)的邊坡穩(wěn)定性敏感性分析[J].水利學(xué)報(bào)，2007，37(11)：1397-1401.

[5]夏元友，熊海豐.邊坡穩(wěn)定性影響因素敏感性人工神經(jīng)網(wǎng)絡(luò)分析[J].巖石力學(xué)與工程學(xué)報(bào)，2004，23(16)：2704-2707.

[6]劉漢東.邊坡失穩(wěn)定時(shí)預(yù)報(bào)理論與方法[M].鄭州：黃河水利出版社，1996.

[7]馬顯春，王 雷，趙法鎖.滑坡穩(wěn)定影響因子敏感性分析及治理方案探討[J].地質(zhì)力學(xué)學(xué)報(bào)，2008,14(4):381-388.

[8]趙瑜,李曉紅,盧義玉,等.深埋隧道圍巖應(yīng)變軟化模型參數(shù)的正交設(shè)計(jì)[J].重慶大學(xué)學(xué)報(bào),2008(7):716-719.