王樂 WANG Le；汪洋 WANG Yang

（①鄂北地區(qū)水資源配置建設與管理局（籌），武漢 430062；②孝感市孝南區(qū)應急管理局，孝感 432001）

0 引言

敏感性分析是指當我們對模型進行數(shù)值模擬時，模型參數(shù)的變化對模擬結(jié)果的影響程度，判斷各個模型參數(shù)對模型結(jié)果的重要程度[3]。當前，對邊坡穩(wěn)定性分析中主要考慮的因素有土體容重γ、粘聚力c 及內(nèi)摩擦角φ，如彭小云等[4]基于灰度關聯(lián)分析法分析了影響高邊坡穩(wěn)定的因素，結(jié)果顯示重力及滑動面的強度參數(shù)是影響邊坡穩(wěn)定的主要因素。VG 模型方面國內(nèi)研究主要是VG 模型的應用及參數(shù)的獲取。其中范嚴偉等[5]利用HYDRUS-1D 模型，分析了不同土壤條件下VG 模型中飽和含水量、殘余含水量、系數(shù)及飽和滲透系數(shù)對土壤水分特征曲線的影響，發(fā)現(xiàn)系數(shù)和飽及滲透系數(shù)n 與濕潤鋒運移距離成正相關，而系數(shù)a 及殘余含水量與濕潤鋒成負相關。徐付橋等[6]研究了VG 模型參數(shù)的選取對堿水層二氧化碳封存的影響。楊紅等[7]基于VG 模型研究了有機肥對土壤水分特征曲線的影響，表明VG 模型可以較好擬合土壤水分特征曲線。田東方等[8]研究了尿素濃度對土壤水分特征曲線的影響，建立了可以考慮尿素濃度影響的VG 模型。彭建平等[9]提出了一種用MATLAB 確定VG 模型參數(shù)的方法。

本文利用Carsel and Parrish（1988）[10]文章中砂土、粉土、粉質(zhì)粘土的VG 參數(shù)，然后將原參數(shù)上調(diào)10%、20%以及下調(diào)10%、20%，以GEO-sudio 軟件分析VG 模型中系數(shù)n、α、飽和含水量θs、殘余含水量θr以及飽和滲透系數(shù)Ks，5 個因素對邊坡中土體穩(wěn)定性系數(shù)的影響，并通過正交試驗設計分析進行VG 模型參數(shù)敏感性分析。

1 基本理論

1.1 滲流基本理論

本文以GEO-studio 軟件為主，GeoStuio 軟件是加拿大專業(yè)的巖土分析軟件，由Fredlund 教授20 世紀70年代開始研發(fā)。GEO-studio 中滲流公式為Richards 方程，對于一維非飽和土壤水分運動，Richards 方程常表示為3 種形式，即以壓力水頭、以含水率和以位置坐標為因變量的公式，下面給出以壓力水頭h 表示的非飽和滲流方程[11]：

式中：kij—分別沿水平及垂直方向的滲透系數(shù)，m/s；ρw—水的密度，g/cm3—土壤水特性曲線的斜率；且，其中θw—含水量；uw—水壓力，kPa。

1.2 極限平衡法基本理論

GEO-studio 軟件中SLOPE/W 模塊采用通用剛體極限平衡法。該方法由Fredlund 等人于20 世紀70年代提出。特點是假設每個土條的受力平衡，每個土條安全系數(shù)相等。同時滿足力矩、水平力、豎向力的平衡。通用極限平衡法力矩平衡安全系數(shù)方程：

水平方向靜力平衡安全系數(shù)方程：

式中：c′—有效粘聚力；φ′—有效摩擦角；u—孔隙水壓力；N—土條底部法向力；W—土條重量；D—集中荷載；β、R、x、f、d、w—幾何參數(shù)；a—土條底面傾角。其中N 為土條底部的法向力。

云計算是為了解決用戶的存儲資源不足和共享資源不足問題，將大量的服務器資源通過網(wǎng)絡整合和調(diào)度重新分配給用戶。高校之間、校內(nèi)老師之間由于受到技術、知識產(chǎn)權等客觀因素的影響，課程資源很難實現(xiàn)共享，還可能有課程資源重復建設的現(xiàn)象，基于云計算技術的網(wǎng)絡教學平臺在充分保護知識產(chǎn)權的情況下實現(xiàn)了跨地區(qū)的教學資源同步共享，通過彼此協(xié)商和網(wǎng)絡確認實現(xiàn)相互間教學資源共享和學習交流[2]。

式中：XR、XL—土體兩側(cè)的剪切力。當土為非飽和時，有效粘聚力c 取值如下：

1.3 VG 模型

VG 模型是1980 由Van Genuchten 提出的一種土壤水分特征曲線的模型[1]，到現(xiàn)在已經(jīng)廣泛應用于地下水分析、滲流分析以及各類軟件中。VG 模型表達式[1]：

1.4 正交實驗設計

正交實驗設計的本質(zhì)是一種數(shù)理統(tǒng)計方法，首先確定有代表性的點，然后根據(jù)正交表進行多因素實驗。其正交表表示為：Ln（tm），式中：L—正交表；n—正交試驗次數(shù)；t—影響因子的水平數(shù)；m—正交試驗表列數(shù)。

正交試驗表格完成后，就是進行數(shù)據(jù)分析；正交試驗計算數(shù)據(jù)分析的一般方法有：極差分析法和方差分析法。極差分析方法簡單、快捷、方便，不足之處是精度，只能定性分析；方差分析于極差分析而言精度要好，但相應的計算要復雜一些。本文采用極差分析方法進行計算。

極差分析方法的計算統(tǒng)計參數(shù)公式為[12]：

式中：Kij—影響因子j 在i 水平下每次試驗結(jié)果的平均值；n—影響因子j 在i 水平下的正交試驗次數(shù)；Yk—第k 個正交實驗的值；—所有試驗的平均值。

設極差分析法的敏感性評價標準為Si，則其計算公式為：

Sj越大，則該影響因子越大，反之。

2 VG 模型參數(shù)對邊坡穩(wěn)定參數(shù)的敏感性分析

2.1 模型概況

通過《工程地質(zhì)手冊（第四版）》表3-1-24 得到砂土、粉土和粉質(zhì)黏土的物理參數(shù)，模型參數(shù)如表1，模型如圖1。

圖1 邊坡模型

表1 模型參數(shù)

邊坡模型邊坡底長60m，左高21m，右高6m，坡高10m。其余參數(shù)如圖1 所示。圖中模型全局網(wǎng)絡單元尺寸為1m，網(wǎng)格四邊形與三邊形一共826 個單元，892 個節(jié)點。

2.2 參數(shù)θr、θs、α、n 以及Ks 的擾動對邊坡穩(wěn)定影響

根據(jù)Carsel and Parrish（1988）[5]得砂土、粉土及粉質(zhì)黏土，以砂石為例對其VG 參數(shù)進行微調(diào)得到其正交實驗水平取值如表2。

表2 正交試驗因素水平取值

由于本文是5 因素5 水平，這里選擇L25（56）正交表，使用GEO-studio 軟件計算邊坡的安全系數(shù)，其正交表及計算結(jié)果見表3。

表3 砂土正交設計方案及安全系數(shù)結(jié)果

通過對正交表數(shù)據(jù)分析計算得到極差分析方法數(shù)據(jù)如表4 所示。

表4 砂土極差分析表

通過對極差分析方法表可以看出：在砂土中，VG 模型中飽和含水量θs對邊坡穩(wěn)定性系數(shù)影響最大，其次是系數(shù)n，再次是殘余含水量θr，然后是系數(shù)α，對邊坡穩(wěn)定性系數(shù)影響最小的則是Ks。

同理得到粉土的極差分析表如表5。

表5 粉土極差分析表

通過上面的正交試驗可以知道對于粉土，VG 模型中飽和含水量θs對邊坡穩(wěn)定性系數(shù)影響最大，系數(shù)n 次之，再次是系數(shù)α，飽和滲透系數(shù)Ks和殘余含水量θr對邊坡穩(wěn)定性系數(shù)影響最小。

對于粉質(zhì)黏土由于VG 模型參數(shù)n＞1，所以只有3 組水平，即5 因素3 水平，這里使用L18（37）正交表，其余過程同砂石一樣，算得粉質(zhì)黏土極差分析表如表6。

表6 粉質(zhì)黏土極差分析表

通過極差分析表可以得到：對于粉質(zhì)黏土，殘余含水量θr對邊坡穩(wěn)定性系數(shù)影響最大，然后是飽和含水量θs，再次是系數(shù)n，最后是系數(shù)α，飽和滲透系數(shù)Ks。

對砂土和粉土，VG 模型參數(shù)中的飽和含水量θs和系數(shù)n 的擾動對邊坡的穩(wěn)定性系數(shù)影響較大；系數(shù)α、飽和滲透系數(shù)Ks和殘余含水量θr對邊坡穩(wěn)定性系數(shù)影響較小。對粉質(zhì)黏土，則是殘余含水量θr和飽和含水量θs的擾動對邊坡穩(wěn)定性系數(shù)影響較大；系數(shù)n、最后是系數(shù)α 以及飽和滲透系數(shù)Ks的影響較小。造成這樣結(jié)果的原因可能是：砂土與粉土的粘聚力小，水在下滲的過程中砂土和黏土的粘聚力有一個增大的過程，然后隨著水下滲量的增加砂土與粉土的粘聚力才開始減小；而對于粉質(zhì)黏土其粘聚力較大，因此水在下滲過程中粉質(zhì)黏土的粘聚力直接就開始減小。

VG 模型在這里可以發(fā)現(xiàn)主要是影響土中水的分布以及土中各處的水頭，從而對邊坡的穩(wěn)定性造成影響。本文通過對正交試驗驗證了VG 模型參數(shù)對邊坡穩(wěn)定性系數(shù)的敏感度，具有一定的工程實際意義。

3 結(jié)語

本文選取了砂土、粉土和粉質(zhì)黏土的VG 模型參數(shù)，研究了其VG 模型參數(shù)的擾動對邊坡穩(wěn)定性的影響，通過GEO-sudio 軟件模擬得到不同VG 系數(shù)下邊坡的安全系數(shù)，設計L25（56）及L18（37）正交試驗表分析各個系數(shù)的擾動對邊坡穩(wěn)定性影響，得到敏感度關系，對于砂土是：θs>n>θr>α>Ks；粉土：θs>n>α>Ks>θr；粉質(zhì)黏土：θr>θs>n>α>Ks。即對砂土和粉土，VG 模型參數(shù)中的飽和含水量θs的擾動對邊坡的穩(wěn)定性系數(shù)影響最大；系數(shù)n 次之，系數(shù)α、飽和滲透系數(shù)Ks和殘余含水量θr對砂土和粉土邊坡穩(wěn)定性系數(shù)影響較小。對粉質(zhì)黏土而言，則是殘余含水量θr的擾動對邊坡穩(wěn)定性系數(shù)影響最大；飽和含水量θs次之，再次是系數(shù)n，最后是系數(shù)α 以及飽和滲透系數(shù)Ks的影響較小。