李丞　蔡立明　張偉鋒　侯建香　董凱　蘇立君

摘 要：渣土場(chǎng)是容納棄土和棄渣的主要場(chǎng)所，由于堆填材料和壓實(shí)方法的差異性，研究初始滲流場(chǎng)對(duì)計(jì)算渣土場(chǎng)邊坡的降雨入滲特征及其穩(wěn)定性的影響具有現(xiàn)實(shí)意義。以深圳市部九窩渣土場(chǎng)邊坡為研究對(duì)象，運(yùn)用非飽和滲流原理和非飽和土強(qiáng)度理論建立數(shù)值計(jì)算模型，設(shè)計(jì)了不同初始滲流條件和不同降雨強(qiáng)度的計(jì)算方案，揭示了渣土場(chǎng)邊坡在降雨作用下的滲流特征及其安全系數(shù)變化規(guī)律。結(jié)果表明：初始滲流場(chǎng)和降雨強(qiáng)度對(duì)渣土場(chǎng)邊坡的降雨入滲特征和穩(wěn)定性影響顯著;豎直方向流速峰值和孔隙水壓力上升速率與初始基質(zhì)吸力呈正相關(guān)關(guān)系;當(dāng)降雨強(qiáng)度為3.5 mm/h，監(jiān)測(cè)點(diǎn)1在初始基質(zhì)吸力為-25、-50、-75 kPa時(shí)達(dá)到飽和的時(shí)間分別為10.9、20.2、25.5 h;豎直方向流速、體積含水率上升速度和孔隙水壓力上升速率與降雨強(qiáng)度呈正相關(guān)關(guān)系;在初始基質(zhì)吸力為-25 kPa和降雨強(qiáng)度為3.5 mm/h時(shí)，渣土場(chǎng)邊坡安全系數(shù)下降明顯，需要采取必要的工程措施。

關(guān)鍵詞：渣土場(chǎng);邊坡;滲流場(chǎng);滲流特征;邊坡穩(wěn)定性

中圖分類號(hào)：TU431? ?文獻(xiàn)標(biāo)志碼：A ??文章編號(hào)：2096-6717（2021）02-0001-9

Abstract： The municipal solid waste landfill is the main place to hold the waste soil and slag.Due to the difference? in the filling material and the compaction method， it is of practical significance to study the influence of initial seepage field on the calculation of rainfall infiltration characteristics and stability of the municipal solid waste landfill.Therefore， taking the municipal solid waste landfill in Bujiuwo， Shenzhen as the research object， a numerical calculation model was established based on the unsaturated seepage principle and unsaturated soil strength theory.The calculation schemes of different initial seepage conditions and rainfall intensities are designed. Then the seepage characteristics and the variation law of safety factor of the slope under rainfall are revealed. The results show that the initial seepage field and rainfall intensity have significant influence on the rainfall infiltration characteristics and therefore stability of the municipal solid waste landfill. The peak value of vertical velocity and the increase rate of pore water pressure are positively correlated with the initial matrix suction. When the rainfall intensity is 3.5 mm/h， the saturation time of monitoring point 1 is 10.9， 20.2， 25.5 h? and the corresponding initial matrix suction? equals to -25， -50， -75 kPa， respectively.Rainfall intensity，velocity of flow and volume water content? are positively correlated with pore water pressure in vertical direction.When the initial matrix suction is -25 kPa and the rainfall intensity is 3.5 mm/h， the safety factor of the waste dump slope decreases obviously，? and necessary engineering measures should be taken.

Keywords：municipal solid waste landfill; slope; seepage field; seepage characteristics; slope stability

隨著城市化進(jìn)程的加快，產(chǎn)生了以大量廢土和廢渣為主的余泥渣土，對(duì)城市固體廢棄物的管理和城市規(guī)劃提出了挑戰(zhàn)[1]，而填埋是最常見(jiàn)的渣土處理方法[2]。對(duì)于渣土場(chǎng)邊坡而言，其穩(wěn)定性受填料性質(zhì)、施工工藝、降雨、地下水位等影響[3-4]，在眾多因素中，降雨入滲將導(dǎo)致渣土場(chǎng)邊坡含水率增加、基質(zhì)吸力降低和地下水位上升，這是渣土邊坡變形與穩(wěn)定性下降的主要原因之一[5-7]。

降雨作用下渣土場(chǎng)邊坡的滲流特征和穩(wěn)定性的研究主要集中在幾個(gè)方面：1）降雨強(qiáng)度和初始含水率對(duì)滲濾液產(chǎn)量的影響[8-9];2）降雨模式和降雨強(qiáng)度對(duì)邊坡水分運(yùn)移及其穩(wěn)定性的影響[10-11];3）降雨模式和降雨強(qiáng)度對(duì)填埋場(chǎng)覆蓋層滲流特性的影響[12-13];4）降雨強(qiáng)度對(duì)渣土場(chǎng)邊坡變形破壞特征的影響[14-15]。以上成果對(duì)于理解渣土場(chǎng)邊坡的水分運(yùn)移規(guī)律和穩(wěn)定性計(jì)算起到了關(guān)鍵的作用。然而，上述研究鮮有涉及初始滲流場(chǎng)對(duì)渣土場(chǎng)邊坡降雨入滲機(jī)制的影響。值得注意的是，在常規(guī)邊坡的研究中，學(xué)者們分析了初始滲流場(chǎng)對(duì)降雨作用下邊坡滲流特征和穩(wěn)定性的影響[16-18]。其中，唐棟等[17]研究了不同初始滲流條件對(duì)不同土體邊坡穩(wěn)定性的影響，結(jié)果表明，初始滲流條件對(duì)邊坡安全系數(shù)影響很大，并建議將多年平均降雨量對(duì)應(yīng)的穩(wěn)態(tài)滲流場(chǎng)作為初始滲流條件。何忠明等[19]比較了4類初始滲流場(chǎng)對(duì)邊坡降雨入滲過(guò)程的影響，發(fā)現(xiàn)初始滲流場(chǎng)分布形式對(duì)降雨過(guò)程中邊坡表面孔隙水壓力的影響較小，且初始基質(zhì)吸力越大，受到降雨入滲的影響程度則越小。然而，由于填筑材料的差異性，渣土場(chǎng)邊坡的降雨滲流特征和穩(wěn)定性分析較常規(guī)邊坡復(fù)雜[14，20-21]，因此，初始滲流場(chǎng)對(duì)渣土場(chǎng)邊坡降雨入滲和穩(wěn)定性計(jì)算的影響值得深入研究。

筆者以深圳市部九窩渣土場(chǎng)邊坡為研究對(duì)象，結(jié)合室內(nèi)土工試驗(yàn)及實(shí)際降雨資料，采用數(shù)值分析方法，在探討不同初始滲流場(chǎng)和不同降雨強(qiáng)度對(duì)渣土場(chǎng)邊坡的滲流特征影響的基礎(chǔ)上，研究了渣土場(chǎng)邊坡的安全系數(shù)變化規(guī)律。

1 計(jì)算原理

1.1 非飽和滲流原理

邊坡降雨入滲過(guò)程可認(rèn)為是典型的飽和非飽和滲流過(guò)程，通常認(rèn)為，在非飽和狀態(tài)下達(dá)西定律依然適用，介質(zhì)水力傳導(dǎo)度是含水率的函數(shù)，Richards等[22]建立了土體中液體瞬時(shí)流動(dòng)的二維控制方程

2 有限元模型及計(jì)算方案

2.1 工程實(shí)例概況

選取深圳市部九窩人工渣土場(chǎng)邊坡作為研究對(duì)象，依據(jù)壓實(shí)度和施工工藝不同，可劃分為3個(gè)區(qū)：渣土區(qū)（壓實(shí)度為80%，坡比為1∶1.275）、填埋嚴(yán)控區(qū)（壓實(shí)度為85%，坡比為1∶1.25）和壩體區(qū)（壓實(shí)度為90%，坡比為1∶1.25），圖1為渣土場(chǎng)邊坡2017年12月的平面圖。研究場(chǎng)地屬于亞熱帶季風(fēng)氣候區(qū)，年平均溫度22 ℃，年平均降雨量1 900 mm，降水量集中在三月到九月。

2.2 邊界條件及計(jì)算模型

根據(jù)工程地質(zhì)調(diào)查確定圖1中的AB斷面為計(jì)算采用的剖面，剖面圖如圖2所示。采用理想彈塑性模型和Mohr-Coulomb屈服準(zhǔn)則在Geo-studio中建立數(shù)值分析模型[24]，模型共劃分13 232個(gè)節(jié)點(diǎn)和13 250個(gè)單元。在模型中設(shè)置監(jiān)測(cè)截面1、2、3，每個(gè)截面的高度為8.8 m，其中，在監(jiān)測(cè)截面1和2坡面處各有一個(gè)監(jiān)測(cè)點(diǎn)，在監(jiān)測(cè)截面3坡面處設(shè)置監(jiān)測(cè)點(diǎn)3，在離監(jiān)測(cè)點(diǎn)3的豎直方向4.4 m處設(shè)置監(jiān)測(cè)點(diǎn)4，各監(jiān)測(cè)截面和監(jiān)測(cè)點(diǎn)用來(lái)監(jiān)測(cè)邊坡在降雨過(guò)程中的滲流特征，監(jiān)測(cè)方案如圖2所示。初始地下水位位于巖土分界面，地下水位線以上土體處于飽和非飽和狀態(tài)。在降雨滲流計(jì)算中，定義坡面為單位流量邊界以模擬降雨入滲過(guò)程，同時(shí)，保證兩側(cè)和底部為不透水邊界。

根據(jù)室內(nèi)土工試驗(yàn)所得結(jié)果，測(cè)得各巖土體的飽和滲透系數(shù)值和飽和體積含水率，并用式（2）和式（3）所示VG模型擬合得到巖土體非飽和滲透系數(shù)與基質(zhì)吸力的關(guān)系曲線，如圖3所示。此外，當(dāng)土體的體積含水率達(dá)到其自身飽和體積含水率的90%時(shí)，則認(rèn)為土體已經(jīng)達(dá)到飽和狀態(tài)[19]，因此，沒(méi)有特別說(shuō)明時(shí)，體積含水率大于或等于飽和含水率的90%的區(qū)域均視為飽和區(qū)范圍，即暫態(tài)飽和區(qū)（Transient Saturated Zone，簡(jiǎn)稱TSZ）。巖土體物理力學(xué)參數(shù)見(jiàn)表1。

2.3 計(jì)算方案

結(jié)合深圳地區(qū)區(qū)域氣象條件，設(shè)計(jì)了如表2所示的3種工況。

由于直接測(cè)定坡體內(nèi)的初始基質(zhì)吸力分布比較困難，而且費(fèi)時(shí)費(fèi)力[16]，為了研究初始滲流條件對(duì)邊坡滲透特性和穩(wěn)定性的影響，參考前人研究[16-19]，假設(shè)地下水位以上土體基質(zhì)吸力呈線性增加，達(dá)到吸力上限后保持常數(shù)。假設(shè)坡體內(nèi)基質(zhì)吸力上限分別為-25、-50、-75 kPa，并通過(guò)微小降雨量形成穩(wěn)態(tài)滲流達(dá)到上限值[25]，將穩(wěn)態(tài)滲流形成的基質(zhì)吸力作為初始滲流條件，如圖4所示。3種工況在進(jìn)行降雨入滲瞬態(tài)分析時(shí)均設(shè)置3種降雨強(qiáng)度q1、q2和q3，以此研究不同降雨強(qiáng)度下初始滲流場(chǎng)對(duì)邊坡降雨入滲特征和穩(wěn)定性的影響。

3 結(jié)果與分析

3.1 邊坡監(jiān)測(cè)點(diǎn)流速變化規(guī)律

當(dāng)降雨強(qiáng)度為q1～q3（1.5～3.5 mm/h）和不同初始滲流場(chǎng)時(shí)，監(jiān)測(cè)點(diǎn)1、2、3處Y方向流速隨降雨持時(shí)的變化關(guān)系如圖5所示。由圖5可知，監(jiān)測(cè)點(diǎn)1、2、3的流速均是先增大后減小，最后略大于降雨強(qiáng)度，且同一情況下，監(jiān)測(cè)點(diǎn)2的流速最大，監(jiān)測(cè)點(diǎn)3的流速次之，監(jiān)測(cè)點(diǎn)1的流速最小;對(duì)于同一初始基質(zhì)吸力和監(jiān)測(cè)點(diǎn)，降雨強(qiáng)度越大，流速的峰值越大，且達(dá)到峰值的時(shí)間越短;對(duì)于同一降雨強(qiáng)度和監(jiān)測(cè)點(diǎn)，初始基質(zhì)吸力越大，流速的峰值越大，且達(dá)到峰值的時(shí)間越長(zhǎng)。此外，當(dāng)初始基質(zhì)吸力為-25 kPa時(shí)，降雨強(qiáng)度越小，各監(jiān)測(cè)點(diǎn)間流速的數(shù)值差距越小，而當(dāng)初始基質(zhì)吸力減小至-75 kPa后，各監(jiān)測(cè)點(diǎn)間流速的數(shù)值差距越大。

3.2 邊坡監(jiān)測(cè)點(diǎn)暫態(tài)飽和區(qū)形成時(shí)間規(guī)律

當(dāng)降雨強(qiáng)度為q1～q3（1.5～3.5 mm/h）和不同初始滲流場(chǎng)時(shí)，監(jiān)測(cè)點(diǎn)1、2、3處的體積含水率隨降雨持時(shí)的變化關(guān)系如圖6所示。由圖6可知，對(duì)于同一初始基質(zhì)吸力和監(jiān)測(cè)點(diǎn)，降雨強(qiáng)度越大，體積含水率上升速度越快，達(dá)到飽和的時(shí)間越短，降雨強(qiáng)度越小，暫態(tài)飽和區(qū)形成的時(shí)間越長(zhǎng)，甚至不會(huì)出現(xiàn)暫態(tài)飽和區(qū)，以圖6（a）中的監(jiān)測(cè)點(diǎn)1為例，監(jiān)測(cè)點(diǎn)1在q2和q3時(shí)達(dá)到飽和的時(shí)間分別為29、10.9 h，而在q1時(shí)沒(méi)有達(dá)到飽和;對(duì)于同一降雨強(qiáng)度和監(jiān)測(cè)點(diǎn)，初始基質(zhì)吸力越大，體積含水率上升速度越慢，達(dá)到飽和的時(shí)間越長(zhǎng)，以圖6中的監(jiān)測(cè)點(diǎn)1為例（降雨強(qiáng)度q3），監(jiān)測(cè)點(diǎn)1在初始基質(zhì)吸力為-25、-50、-75 kPa時(shí)達(dá)到飽和的時(shí)間分別為10.9、20.2、25.5 h。

圖7為當(dāng)降雨強(qiáng)度為q3時(shí)，不同深度監(jiān)測(cè)點(diǎn)的體積含水率和Y方向流速的變化。由圖7可知，斜坡表面監(jiān)測(cè)點(diǎn)（監(jiān)測(cè)點(diǎn)3）的流速與體積含水率變化要快于坡體內(nèi)部的監(jiān)測(cè)點(diǎn)（監(jiān)測(cè)點(diǎn)4），表明暫態(tài)飽和區(qū)從坡表向坡內(nèi)擴(kuò)展;對(duì)監(jiān)測(cè)點(diǎn)3而言，體積含水率隨著流速的增加而緩緩增大，當(dāng)流速達(dá)到峰值時(shí)，體積含水率也達(dá)到絕對(duì)飽和（28%），隨后流速開(kāi)始下降，體積含水率開(kāi)始保持不變。值得注意的是，監(jiān)測(cè)點(diǎn)3在120 h時(shí)的流速和體積含水率發(fā)生了突變（A點(diǎn)和B點(diǎn)），這是由于降雨導(dǎo)致了地下水位大幅度上升，且地下水已經(jīng)淹沒(méi)了監(jiān)測(cè)點(diǎn)3。

3.3 邊坡暫態(tài)飽和區(qū)形成深度規(guī)律

當(dāng)降雨強(qiáng)度在q1～q3（1.5～3.5 mm/h）和不同初始滲流場(chǎng)時(shí)，監(jiān)測(cè)截面1、2、3處的暫態(tài)飽和區(qū)深度隨降雨持時(shí)的變化關(guān)系如圖8所示。由圖8可知，對(duì)于同一初始基質(zhì)吸力和監(jiān)測(cè)截面，降雨強(qiáng)度越大，暫態(tài)飽和區(qū)深度越大;對(duì)于同一降雨強(qiáng)度和監(jiān)測(cè)截面，初始基質(zhì)吸力越小，入滲深度越大，截面處暫態(tài)飽和區(qū)出現(xiàn)的范圍越大，體積含水率變化范圍越小，體積含水率在同一位置增大的幅度也越大。此外，當(dāng)降雨強(qiáng)度為q1時(shí)，截面1和2都沒(méi)有產(chǎn)生暫態(tài)飽和區(qū)，截面3的暫態(tài)飽和區(qū)深度也較小，最大值僅為1.15 m，而當(dāng)降雨強(qiáng)度為q3時(shí)，截面3的暫態(tài)飽和區(qū)深度最大值卻達(dá)8.79 m，由此可見(jiàn)，降雨強(qiáng)度會(huì)顯著影響邊坡是否會(huì)出現(xiàn)暫態(tài)飽和區(qū)，而初始基質(zhì)吸力則會(huì)對(duì)暫態(tài)飽和區(qū)深度的數(shù)值產(chǎn)生影響。

3.4 邊坡孔隙水壓力變化規(guī)律

當(dāng)降雨強(qiáng)度為q1～q3（1.5～3.5 mm/h）和不同初始滲流場(chǎng)時(shí)，監(jiān)測(cè)點(diǎn)1、2、3處的孔隙水壓力隨降雨持時(shí)的變化關(guān)系如圖9所示。由圖9可知，隨著降雨時(shí)間的持續(xù)，坡面孔隙水壓力出現(xiàn)了升高的趨勢(shì)。具體表現(xiàn)為：對(duì)于同一初始基質(zhì)吸力和監(jiān)測(cè)點(diǎn)處，降雨強(qiáng)度越大，入滲深度越大，孔隙水壓力上升速率越快，降雨結(jié)束后的孔隙水壓力也越小;對(duì)于同一降雨強(qiáng)度和監(jiān)測(cè)點(diǎn)處，初始基質(zhì)吸力較大的邊坡孔隙水壓力上升速率快于初始基質(zhì)吸力較低的邊坡，同時(shí)，坡面孔隙水壓力升高速率整體上逐漸減小并趨于水平，就降雨持續(xù)時(shí)間內(nèi)孔隙水壓力的數(shù)值而言，初始滲流條件對(duì)各情況下降雨結(jié)束后的孔隙水壓力的影響較小。由此可見(jiàn)，降雨強(qiáng)度會(huì)顯著影響孔隙水壓力增長(zhǎng)的快慢，而初始基質(zhì)吸力影響較小。

3.5 降雨條件下邊坡穩(wěn)定性分析

由前述可知，在不同初始滲流場(chǎng)和降雨強(qiáng)度的影響下，降雨過(guò)程中邊坡的滲流特性和孔隙水壓力變化規(guī)律有所差異，因此，針對(duì)圖2所示邊坡，探討其邊坡穩(wěn)定性的變化規(guī)律。當(dāng)不考慮降雨時(shí)，對(duì)邊坡初始安全系數(shù)進(jìn)行計(jì)算，最危險(xiǎn)滑動(dòng)面通過(guò)SLOPE/W程序自動(dòng)搜索確定，采用有限單元應(yīng)力法求解安全系數(shù)，以初始基質(zhì)吸力為-25 kPa時(shí)為例，邊坡安全系數(shù)為1.623，最危險(xiǎn)滑動(dòng)面如圖10所示。由圖10可知，邊坡最危險(xiǎn)滑動(dòng)面同時(shí)位于壩體區(qū)和填埋嚴(yán)控區(qū)內(nèi)，考慮到壩體潰決的可能性，為簡(jiǎn)單起見(jiàn)，在降雨持續(xù)過(guò)程中，分析壩體區(qū)內(nèi)（即監(jiān)測(cè)截面3）的暫態(tài)飽和區(qū)厚度的變化規(guī)律對(duì)邊坡穩(wěn)定性的影響。

降雨強(qiáng)度為q1～q3（1.5～3.5 mm/h）和不同初始滲流場(chǎng)時(shí)的邊坡安全系數(shù)變化如圖11所示。由圖11可知，暫態(tài)飽和區(qū)厚度對(duì)邊坡穩(wěn)定性影響顯著，且邊坡初始基質(zhì)吸力為-50、-75 kPa時(shí)，邊坡安全系數(shù)下降趨勢(shì)和暫態(tài)飽和區(qū)厚度增加趨勢(shì)均呈線性。此外，由圖11（a）可知，當(dāng)降雨強(qiáng)度為q1時(shí)，邊坡沒(méi)有產(chǎn)生暫態(tài)飽和區(qū)，且初始基質(zhì)吸力為-25 kPa時(shí)的安全系數(shù)下降速率最大;由圖11（b）可知，降雨強(qiáng)度為q2時(shí)暫態(tài)飽和區(qū)的最大厚度分別為6.00、2.15、1.34 m，降雨結(jié)束后，邊坡安全系數(shù)分別為1.439、1.653、1.693;由圖11（c）可知，當(dāng)降雨強(qiáng)度為q3時(shí)，初始基質(zhì)吸力為-25 kPa的邊坡安全系數(shù)下降呈非線性趨勢(shì)，最小值為1.241，相應(yīng)地，其暫態(tài)飽和區(qū)厚度增加也呈非線性趨勢(shì)，最大值為8.79 m。綜合圖11可知，在相同初始基質(zhì)吸力條件下，暫態(tài)飽和區(qū)厚度為8.79 m時(shí)的安全系數(shù)與沒(méi)有暫態(tài)飽和區(qū)時(shí)的安全系數(shù)最大相差0.272，降幅約21.91%;在相同降雨強(qiáng)度下，暫態(tài)飽和區(qū)厚度為8.79 m時(shí)的安全系數(shù)與暫態(tài)飽和區(qū)厚度為2.27 m時(shí)的安全系數(shù)最大相差0.429，降幅約35.37%。

當(dāng)初始基質(zhì)吸力為-25 kPa時(shí)，邊坡初始安全系數(shù)為1.623，具備很好的安全儲(chǔ)備，而經(jīng)歷120 h的強(qiáng)降雨后（q3=3.5 mm/h），安全系數(shù)只有1.241，且結(jié)合監(jiān)測(cè)點(diǎn)暫態(tài)飽和區(qū)形成時(shí)間規(guī)律可知，在降雨作用下，壩體區(qū)的地下水位大幅度上升，并已達(dá)到壩體區(qū)坡頂。目前，中國(guó)關(guān)于渣土場(chǎng)邊坡相關(guān)的設(shè)計(jì)規(guī)范還未成熟，若借鑒碾壓土石壩的設(shè)計(jì)標(biāo)準(zhǔn)：壩坡正常條件下的最小安全系數(shù)為1.25（簡(jiǎn)化畢曉普法），非正常運(yùn)用條件Ⅰ時(shí)的最小安全系數(shù)為1.15（簡(jiǎn)化畢曉普法）[26]。由此可知，在初始基質(zhì)吸力更小或遭遇強(qiáng)度更大的降雨時(shí)，若沒(méi)有相應(yīng)的工程措施，渣土場(chǎng)邊坡具有失穩(wěn)風(fēng)險(xiǎn)。

4 討論

基于非飽和滲流原理和非飽和土強(qiáng)度理論建立數(shù)值模型，探討了初始滲流場(chǎng)和降雨強(qiáng)度對(duì)渣土場(chǎng)邊坡的降雨入滲特征及穩(wěn)定性的影響，得到的結(jié)論為渣土場(chǎng)的長(zhǎng)期運(yùn)營(yíng)及排水防滲工作提供了參考依據(jù)。

一般而言，在數(shù)值模擬研究中，常規(guī)邊坡多為單級(jí)邊坡，且構(gòu)成物質(zhì)單一，其在降雨作用下的滲流特征和穩(wěn)定性分析已經(jīng)比較明確。然而，渣土場(chǎng)邊坡基本為多級(jí)邊坡，其填筑過(guò)程分階段、分區(qū)域進(jìn)行，這使得不同區(qū)域的巖土體力學(xué)參數(shù)存在差異，從而改變了渣土場(chǎng)邊坡的滲流特征，并顯著影響邊坡穩(wěn)定性。具體體現(xiàn)為：需要綜合分析渣土場(chǎng)邊坡的整體穩(wěn)定性和渣及其各自區(qū)域內(nèi)單極邊坡的穩(wěn)定性。從文獻(xiàn)[15]可知，在降雨作用下，渣土場(chǎng)邊坡下部的安全系數(shù)要小于其上部和中部的安全系數(shù)，為簡(jiǎn)單起見(jiàn)，只計(jì)算了渣土場(chǎng)邊坡下部的安全系數(shù)，分析了渣土場(chǎng)邊坡發(fā)生“潰壩”的失穩(wěn)隱患。綜合分析渣土場(chǎng)邊坡的多種失效模式將是下一步研究的工作重點(diǎn)。

另一方面，考慮初始滲流場(chǎng)影響的渣土場(chǎng)邊坡降雨入滲過(guò)程分析是從數(shù)值模擬的角度進(jìn)行的，一定程度上可以為實(shí)際場(chǎng)地的相關(guān)分析提供理論支撐，在后續(xù)研究中應(yīng)通過(guò)現(xiàn)場(chǎng)監(jiān)測(cè)和模型試驗(yàn)等手段對(duì)結(jié)果進(jìn)行檢驗(yàn)與修正。

5 結(jié)論

1）初始滲流場(chǎng)和降雨強(qiáng)度對(duì)渣土場(chǎng)邊坡的降雨入滲特征有顯著影響。初始基質(zhì)吸力越大，Y方向流速的峰值越大，達(dá)到峰值的時(shí)間越長(zhǎng)，且體積含水率上升速度越慢，達(dá)到飽和的時(shí)間越長(zhǎng)。降雨強(qiáng)度越大，Y方向流速的峰值越大，達(dá)到峰值的時(shí)間越短，且體積含水率上升速度越快，達(dá)到飽和的時(shí)間越短。

2）渣土場(chǎng)邊坡暫態(tài)飽和區(qū)的滲流特征可描述為：壩體區(qū)坡面以下土體迅速飽和，暫態(tài)飽和區(qū)范圍由壩體區(qū)坡腳附近加速向上部土體擴(kuò)散，并逐漸覆蓋整個(gè)壩體區(qū)，同時(shí)由壩體區(qū)向填埋嚴(yán)控區(qū)漫延。

3）初始基質(zhì)吸力較大的邊坡孔隙水壓力上升速率快于初始基質(zhì)吸力較低的邊坡。渣土區(qū)的孔隙水壓力上升速率最快，填埋嚴(yán)控區(qū)和壩體區(qū)的孔隙水壓力上升速率差別不大。降雨強(qiáng)度會(huì)顯著影響孔隙水壓力增長(zhǎng)的快慢，而初始基質(zhì)吸力影響較小。

4）初始滲流場(chǎng)對(duì)渣土場(chǎng)邊坡安全系數(shù)有顯著影響。邊坡安全系數(shù)與暫態(tài)飽和區(qū)厚度的關(guān)系為負(fù)相關(guān)。在初始基質(zhì)吸力為-25 kPa時(shí)，經(jīng)歷了120 h降雨后（3.5 mm/h），渣土場(chǎng)邊坡安全系數(shù)下降明顯，具有一定的失穩(wěn)風(fēng)險(xiǎn)，因此，在渣土場(chǎng)進(jìn)行填筑的過(guò)程中，為確保壩坡長(zhǎng)期穩(wěn)定，應(yīng)充分考慮區(qū)域氣象條件對(duì)渣土場(chǎng)運(yùn)營(yíng)期穩(wěn)定性的影響，做好渣土場(chǎng)邊坡排水和防滲設(shè)施，防止發(fā)生壩體區(qū)潰壩。

參考文獻(xiàn)：

[1] HUANG Y， CHENG H L. A simplified analytical model for run-out prediction of flow slides in municipal solid waste landfills [J]. Landslides， 2017， 14（1）： 99-107.

[2] BRUNNER P H， FELLNER J. Setting priorities for waste management strategies in developing countries [J]. Waste Management & Research， 2007， 25（3）： 234-240.

[3] XU Q， PENG D L， LI W L， et al. The catastrophic landfill flowslide at Hongao dumpsite on 20 December 2015 in Shenzhen， China [J]. Natural Hazards and Earth System Sciences， 2017， 17（2）： 277-290.

[4] YIN Y P， LI B， WANG W P， et al. Mechanism of the December 2015 catastrophic landslide at the Shenzhen landfill and controlling geotechnical risks of urbanization [J]. Engineering， 2016， 2（2）： 230-249.

[5] MOHAN Y W， DARLINGTON J. Effects of monsoonal rainfall on waste dump stability and respective geo-environmental issues： A case study [J]. Environmental Earth Sciences， 2011， 63（6）： 1169-1177

[6] BLIGHT G E， FOURIE A B. Catastrophe revisited - disastrous flow failures of mine and municipal solid waste [J]. Geotechnical and Geological Engineering， 2005， 23（3）： 219-248.

[7] CHO Y C， SONG Y S. Deformation measurements and a stability analysis of the slope at a coal mine waste dump [J]. Ecological Engineering， 2014， 68：189-199.

[8] 蘭吉武， 詹良通， 李育超， 等. 填埋垃圾初始含水率對(duì)滲濾液產(chǎn)量的影響及修正滲濾液產(chǎn)量計(jì)算公式[J]. 環(huán)境科學(xué)， 2012， 33（4）： 1389-1396.

LAN J W， ZHAN L T， LI Y C， et al. Impacts of initial moisture content of MSW waste on leachate generation and modified formula for predicting leachate generation [J]. Environmental Science， 2012， 33（4）： 1389-1396. （in Chinese）

[9] 柯瀚， 陳曉哲， 陳云敏， 等. 基于滲流–壓縮耦合作用的填埋場(chǎng)滲瀝液導(dǎo)排量分析[J]. 巖土工程學(xué)報(bào)， 2013， 35（9）： 1634-1641.

KE H， CHEN X Z， CHEN Y M， et al. Calculation of leachate drainage volume of landfills based on seepage-compression coupling action [J]. Chinese Journal of Geotechnical Engineering， 2013， 35（9）： 1634-1641.（in Chinese）

[10] 邱戰(zhàn)洪， 何春木， 朱兵見(jiàn)， 等. 不同降雨模式下山谷型垃圾填埋場(chǎng)水分運(yùn)移及其穩(wěn)定性研究[J]. 巖土力學(xué)， 2012， 33（10）： 3151-3155， 3170.

QIU Z H， HE C M， ZHU B J， et al. Investigations of water transport in valley-type MSW landfills and their stabilities subjected to various rainfall patterns [J]. Rock and Soil Mechanics， 2012， 33（10）： 3151-3155， 3170. （in Chinese）

[11] 黃剛海. 強(qiáng)降雨入滲下復(fù)雜地形排土場(chǎng)穩(wěn)定性分析[J]. 巖土工程學(xué)報(bào)， 2013， 35（Sup2）： 292-295.

HUANG G H. Stability analysis of waste dump with complex terrain under heavy rainfall infiltration [J]. Chinese Journal of Geotechnical Engineering， 2013， 35（Sup2）： 292-295. （in Chinese）

[12] 馮世進(jìn)， 李夕林， 高麗亞. 不同降雨模式條件下填埋場(chǎng)封頂系統(tǒng)最大飽和深度[J]. 巖土工程學(xué)報(bào)， 2012， 34（5）： 924-930.

FENG S J， LI X L， GAO L Y. Maximum saturated depth overbarrier in landfill cover system under various rainfall patterns [J]. Chinese Journal of Geotechnical Engineering， 2012， 34（5）： 924-930. （in Chinese）

[13] 賈官偉. 固廢堆場(chǎng)終場(chǎng)土質(zhì)覆蓋層中水分運(yùn)移規(guī)律及調(diào)控方法研究[D]. 杭州： 浙江大學(xué)， 2010.

JIA G W. Study on the water transport in the landfill earthen final cover and its controlling method [D]. Hangzhou： Zhejiang University， 2010. （in Chinese）

[14] 李宏達(dá). 降雨影響下排土場(chǎng)邊坡滲流特征及穩(wěn)定性分析[D]. 福州： 福州大學(xué)， 2015.

LI H D. The analysis of dump field slopes seepage characteristics and stability effect by rainfall [D]. Fuzhou： Fuzhou University， 2015. （in Chinese）

[15] 秦宏楠. 紫金山金銅礦排土場(chǎng)滑坡誘發(fā)機(jī)理及監(jiān)測(cè)預(yù)警技術(shù)研究[D]. 北京： 北京科技大學(xué)， 2016.

QIN H N. Study on waste dump landslide inducement mechanism and monitoring and early warning technology of Zijinshan gold and copper mine [D]. Beijing： University of Science and Technology Beijing， 2016. （in Chinese）

[16] 汪華斌， 李建梅， 金怡軒， 等. 降雨誘發(fā)邊坡破壞數(shù)值模擬兩個(gè)關(guān)鍵問(wèn)題的解決方法[J]. 巖土力學(xué)， 2019， 40（2）： 777-784.

WANG H B， LI J M， JIN Y X， et al. The numerical methods for two key problems in rainfall-induced slope failure [J]. Rock and Soil Mechanics， 2019， 40（2）： 777-784. （in Chinese）

[17] 唐棟， 李典慶， 周創(chuàng)兵， 等. 考慮前期降雨過(guò)程的邊坡穩(wěn)定性分析[J]. 巖土力學(xué)， 2013， 34（11）： 3239-3248.

TANG D， LI D Q， ZHOU C B， et al. Slope stability analysis considering antecedent rainfall process [J]. Rock and Soil Mechanics， 2013， 34（11）： 3239-3248. （in Chinese）

[18] 何忠明， 邱祥， 卞漢兵， 等. 考慮初始滲流場(chǎng)影響的坡積土邊坡降雨入滲過(guò)程[J]. 長(zhǎng)安大學(xué)學(xué)報(bào)（自然科學(xué)版）， 2017， 37（5）： 39-48.

HE Z M， QIU X， BIAN H B， et al. Process of rainfall infiltration of residual slope considering the effect on initial seepage field [J]. Journal of Changan University （Natural Science Edition）， 2017， 37（5）： 39-48. （in Chinese）

[19] RAHARDJO H， ONG T H， REZAUR R B， et al. Factors controlling instability of homogeneous soil slopes under rainfall [J]. Journal of Geotechnical and Geoenvironmental Engineering， 2007， 133（12）： 1532-1543.

[20] 蔣中明， 伍忠才， 馮樹(shù)榮， 等. 考慮暫態(tài)飽和區(qū)的邊坡穩(wěn)定性極限平衡分析方法[J]. 水利學(xué)報(bào)， 2015， 46（7）： 773-782.

JIANG Z M， WU Z C， FENG S R， et al. Limit equilibrium method for stability analysis of slope with transient saturated area [J]. Journal of Hydraulic Engineering， 2015， 46（7）： 773-782. （in Chinese）

[21] 汪海濱， 李小春， 米子軍， 等. 排土場(chǎng)空間效應(yīng)及其穩(wěn)定性評(píng)價(jià)方法研究[J]. 巖石力學(xué)與工程學(xué)報(bào)， 2011， 30（10）： 2103-2111.

WANG H B， LI X C， MI Z J， etal. Researches on space effect of waste dumps and its stability evaluation method [J]. Chinese Journal of Rock Mechanics and Engineering， 2011， 30（10）： 2103-2111. （in Chinese）

[22] RICHARDS L A. Capielary conduction of liquids through porous mediums[J]. Physics， 1931， 1（5）： 318-333.

[23] VAN GENUCHTEN M T. A closed-form equation for predicting the hydraulic conductivity of unsaturated soils [J]. Soil Science Society of America Journal， 1980， 44（5）： 892-898.

[24] GEO-SLOPE International Ltd. Seepage modeling with SEEP/W 2007 [M]. Calgary， Alberta， Canada： Geo-Slope International Ltd.， 2010.

[25] 楊攀， 楊軍. 考慮前期降雨的邊坡穩(wěn)定降雨閾值曲面[J]. 巖土力學(xué)， 2015， 36（Sup1）： 169-174.

YANG P， YANG J. Rainfall threshold surface for slopes stability considering antecedent rainfall [J]. Rock and Soil Mechanics， 2015， 36（Sup1）： 169-174. （in Chinese）

[26] 小型水利水電工程碾壓式土石壩設(shè)計(jì)規(guī)范： SL 189—2013 [S]. 北京： 中國(guó)水利水電出版社， 2014.

Design code for rolled earth-rockfill dams in small size water resources and hydroelectric engineering： SL 189-2013 [S]. Beijing： China Water Power Press， 2014. （in Chinese）

（編輯 王秀玲）