姜宇航 戴妙林 劉曉青

摘 要：某土工膜堆石壩的壩基防滲體系包括趾板及壩基帷幕，帷幕后實(shí)測水位偏高，與上游水位連通性較強(qiáng)，需要對其成因進(jìn)行分析?；跐B流觀測資料建立回歸模型進(jìn)行初步分析，運(yùn)用正交試驗(yàn)設(shè)計(jì)和滲流有限元計(jì)算生成訓(xùn)練樣本，對各影響因素進(jìn)行敏感性分析。借助支持向量機(jī)在小樣本中的高度非線性映射能力，建立滲流參數(shù)與滲壓水位的對應(yīng)關(guān)系。采用遺傳算法對支持向量機(jī)模型進(jìn)行參數(shù)優(yōu)化，以支持向量機(jī)預(yù)測值與實(shí)測值誤差作為適應(yīng)度值，對壩基地層和防滲體各區(qū)域滲流參數(shù)進(jìn)行優(yōu)化搜索，并將反演結(jié)果進(jìn)行反饋計(jì)算驗(yàn)證。結(jié)果表明，基于支持向量機(jī)-遺傳算法反演滲流參數(shù)是可行的，帷幕后滲壓水位偏高由趾板與帷幕間的裂縫及壩基滲流各向異性綜合引起，可針對上述薄弱部位進(jìn)行工程處理以降低幕后水位。

關(guān)鍵詞：土工膜堆石壩;回歸模型;正交試驗(yàn)設(shè)計(jì);敏感性分析;支持向量機(jī);遺傳算法;滲流參數(shù)反演

中圖分類號：TV698.1 ? 文獻(xiàn)標(biāo)志碼：A

doi：10.3969/j.issn.1000-1379.2020.07.021

Abstract：The anti-seepage system of a geomembrane rockfill dam includes a toe board and a dam foundation curtain. The measured water level behind the curtain is higher and it has stronger connectivity with the upstream water level， and its cause needs to be analyzed. A regression model was established for preliminary analysis based on the seepage observation data. It used orthogonal experimental design and seepage finite element calculation to generate training samples and conduct sensitivity analysis of each influencing factor. With the support of the highly non-linear mapping capability of support vector machines in small samples， the correspondence relationship of seepage parameters and seepage pressure levels were established. The genetic algorithm was used to optimize the parameters of the support vector machine model and the error between the predicted value and the measured value of the support vector machine was used as the fitness value. The seepage parameters of the dam base layer and the impervious body were optimized and searched. The results were verified by feedback calculations. The results show that inversion of seepage parameters based on support vector machine-genetic algorithm is feasible， and the high level of seepage pressure behind the curtain is caused by the cracks between the toe plate and the curtain and the anisotropy of dam foundation seepage. The above weak parts are engineered to lower the water level behind the scenes.

Key words： geomembrane rockfill dam; regression model; orthogonal test design; sensitivity analysis; support vector machine; genetic algorithm; inversion of seepage parameters

在大壩蓄水和運(yùn)行過程中，受各種因素的共同作用，部分區(qū)域滲透系數(shù)偏大、滲壓水位偏高。為了對其成因進(jìn)行分析，需要結(jié)合壩基地質(zhì)條件和防滲體可能存在的缺陷進(jìn)行研究[1]。張兵等[2]運(yùn)用正交試驗(yàn)設(shè)計(jì)，對云駕嶺礦區(qū)巖層力學(xué)參數(shù)進(jìn)行了敏感性分析與反演;Harrouni等[3]提出了運(yùn)用優(yōu)化理論識別含水層參數(shù)的方法;鄧祥輝[4]采用遺傳算法（GA）和有限元法進(jìn)行了高壩工程滲流反問題研究;魏進(jìn)兵等[5]運(yùn)用遺傳算法優(yōu)化BP神經(jīng)網(wǎng)絡(luò)反演得到了滑坡體飽和滲透系數(shù);萬文新等[6]對壩基防滲墻可能出現(xiàn)的裂縫進(jìn)行了滲流計(jì)算與研究。

土工膜堆石壩的趾板與壩基帷幕的接觸部位相對薄弱，可能存在接觸面裂縫，此外壩基巖體存在的滲流各向異性特性，呈現(xiàn)出與各向同性不同的特點(diǎn)。針對上述因素，筆者基于某土工膜堆石壩實(shí)測資料和數(shù)值計(jì)算等方法，進(jìn)行影響規(guī)律和滲流參數(shù)的反演研究。近年興起的支持向量機(jī)（Support Vector Machine，SVM）是基于數(shù)理統(tǒng)計(jì)理論建立起來的新型學(xué)習(xí)機(jī)器，具有適應(yīng)性強(qiáng)、所需樣本容量小和運(yùn)行速度快等特點(diǎn)，且具有極強(qiáng)的非線性映射能力。通過對壩基滲流回歸模型進(jìn)行初步分析，基于正交設(shè)計(jì)方案與滲流有限元計(jì)算得到的樣本，并借助支持向量機(jī)建立滲壓水位與壩基滲流參數(shù)之間的預(yù)測模型，可以為反演滲流參數(shù)提供快速有效的途徑。

1 工程概況

某土工膜堆石壩由復(fù)合土工膜、擠壓邊墻、墊層、過渡層、主堆石區(qū)、次堆石區(qū)和排水體等組成。擠壓式邊墻混凝土共202層，斷面為梯形，邊墻外側(cè)坡比1∶1.6，內(nèi)側(cè)坡比8∶1，頂寬10 cm，底寬79 cm，墻高40 cm。上游沿壩面擠壓邊墻每隔6 m設(shè)置錨固帶，錨固帶由擠壓邊墻固定，土工膜與錨固帶相連接。壩基設(shè)置防滲帷幕，土工膜與壩基面、帷幕之間采用混凝土趾板連接。壩頂高程515 m，趾板建基面最低高程430 m，最大壩高85 m，壩頂全長362.09 m，壩頂寬度為8.00 m。壩頂以上設(shè)4.2 m高防浪墻，上游壩坡1∶1.6，下游壩坡1∶1.8。

大壩最大剖面滲壓計(jì)DB-A-P-01～10觀測布置見圖1，其中位于土工膜后的滲壓計(jì)DB-A-P-08 10滲壓實(shí)測值為0，說明壩體浸潤面低于儀器所在高程，DB-A-P-06儀器損壞無數(shù)據(jù)。

2 實(shí)測滲壓水位分析

帷幕后從高到低依次布置滲壓計(jì)DB-A-P-03、DB-A-P-02、DB-A-P-01，用于監(jiān)測帷幕后滲壓水位。帷幕下游沿壩基面依次布置滲壓計(jì)DB-A-P-04、DB-A-P-05、DB-A-P-07，用于監(jiān)測壩基面滲壓水位，實(shí)測滲壓水位如圖2所示，帷幕后3個(gè)測點(diǎn)2016年和2017年最高水位、最低水位和年變幅見表1。從實(shí)測滲壓水位來看，DB-A-P-01～03滲壓水位與上游水位相關(guān)性較為密切，而DB-A-P-04～05和DB-A-P-07受上游水位變化程度影響不明顯。DB-A-P-01～03測值和下游水位的差值與上下游水位差值之比分別為0.72、0.69、0.86，數(shù)值偏大，表明這3個(gè)測點(diǎn)水位偏高。

由表2可知，DB-A-P-01 03的上游水位分量較大，且受上游水位變化影響程度從高到低為DB-A-P-03、DB-A-P-01、DB-A-P-02，與測點(diǎn)水位高低順序是一致的。以上分析表明，DB-A-P-03測點(diǎn)與上游庫水連通性相對較強(qiáng)。

4 數(shù)值計(jì)算分析

4.1 計(jì)算模型和參數(shù)

Geo-studio為加拿大GEO-SLOPE公司開發(fā)的二維有限元軟件，其中SEEP/W模塊可用來分析多孔隙介質(zhì)滲流狀態(tài)與孔隙水壓力消散過程。本文涉及土工膜、防滲帷幕的模擬，由于該軟件可模擬各向異性滲流，且計(jì)算結(jié)果較為合理準(zhǔn)確，因此采用其進(jìn)行計(jì)算和研究。假設(shè)壩體滲流為層流并滿足達(dá)西定律，建立有限元計(jì)算模型，模型范圍壩基面上下游側(cè)各取200 m，深度取180 m，并包含壩基防滲帷幕，地基兩側(cè)及底部施加邊界約束。根據(jù)設(shè)計(jì)資料對不同區(qū)域賦予滲透系數(shù)，各材料設(shè)計(jì)滲透系數(shù)見表3。采用2017年1月6日水位工況即上游水位508.25 m、下游水位440.32 m進(jìn)行計(jì)算。鑒于壩體上游鋪設(shè)的土工膜厚度僅為3.5 mm，若以此作為邊長來剖分網(wǎng)格單元，不僅計(jì)算精度無法保證，而且會因單元畸形而造成計(jì)算溢出。為此，在實(shí)際有限元計(jì)算中依據(jù)滲流量等效原理，適當(dāng)增大復(fù)合土工膜的計(jì)算厚度與滲透系數(shù)，使計(jì)算精度不受影響[7]。

4.2 基于正交試驗(yàn)設(shè)計(jì)的滲流參數(shù)敏感性分析

4.2.1 壩基地層與防滲體參數(shù)選取

趾板與土工膜和壩基防滲帷幕連接，其自身的滲透系數(shù)較小，考慮到趾板與帷幕連接部位有可能存在缺陷，故在趾板與帷幕連接處設(shè)置寬度為10 cm的等效裂縫，該值為設(shè)計(jì)軟件所允許的最小寬度，并與上游庫水連通。裂縫滲透系數(shù)采用裂隙滲流立方定律進(jìn)行估算，為6.29×10-3 m/s，故裂縫滲透系數(shù)范圍取1×10-3～1×10-2 m/s。考慮到壩基地層為砂質(zhì)板巖，可能存在各向異性的滲流特性[8]，故將壩基垂直層向滲透系數(shù)Kn、順層向滲透系數(shù)Ks與垂直層向滲透系數(shù)Kn的比值Ks/Kn以及Kn與水平向夾角θ作為影響因素。防滲帷幕深度為30 m，滲透系數(shù)取值范圍為1×10-7～9×10-7 m/s。

4.2.2 基于正交試驗(yàn)方案的敏感性分析

正交試驗(yàn)設(shè)計(jì)是依據(jù)正交性原則，從所有試驗(yàn)中挑選出最具有代表性的試驗(yàn)組進(jìn)行試驗(yàn)，這些代表性的試驗(yàn)組具有均勻分散、齊整可比的特點(diǎn)，極大地簡化了多因素多水平情況下的設(shè)計(jì)試驗(yàn)組數(shù)，是一種高效、快速、經(jīng)濟(jì)的試驗(yàn)方法[9]。將Kn、Ks/Kn、θ、帷幕的滲透系數(shù)K帷幕以及裂縫滲透系數(shù)K裂縫作為影響因素，由此設(shè)計(jì)出25組試驗(yàn)的5因素5水平正交表，壩基滲透系數(shù)分組參數(shù)水平設(shè)置見表4。運(yùn)用有限元軟件對每一組試驗(yàn)進(jìn)行計(jì)算，得到25組數(shù)據(jù)的正交試驗(yàn)表。

通過計(jì)算極差Rj進(jìn)行極差分析可以對正交試驗(yàn)各因素的敏感性進(jìn)行評價(jià)[10]，極差Rj越大，說明因素水平改變對結(jié)果的影響越大，即敏感性越大;極差Rj越小，說明因素水平改變對結(jié)果的影響性越小。正交試驗(yàn)數(shù)據(jù)極差分析結(jié)果見表5，其中R1～R7分別為DB-A-P-01～07測點(diǎn)各因素極差。結(jié)果顯示，因素水平改變對DB-A-P-01～03測值影響程度較大，對DB-A-P-04～05和DB-A-P-07測值影響較小;對DB-A-P-01測值影響程度較大的因素為Ks/Kn、θ、K帷幕;對DB-A-P-02測值影響程度較大的因素為Ks/Kn、K帷幕、θ;對DB-A-P-03測值影響程度較大的因素為Ks/Kn、K裂縫、K帷幕。

5 滲流參數(shù)反演

5.1 構(gòu)建支持向量機(jī)預(yù)測樣本

因?yàn)橛邢拊浖礼EO-slope不支持其他程序調(diào)用，所以在參數(shù)反演優(yōu)化計(jì)算過程中無法直接調(diào)用此軟件計(jì)算測點(diǎn)水位。為解決這個(gè)問題，在反演滲流參數(shù)之前，首先建立因變量y（各監(jiān)測點(diǎn)水位）與自變量x（待反演滲流參數(shù)）之間關(guān)系的預(yù)測模型。借助預(yù)測模型，可以計(jì)算得到待反演滲流參數(shù)任意取值組合的各測點(diǎn)水位預(yù)測值，用于反演分析中預(yù)測值與實(shí)測值誤差的計(jì)算。當(dāng)因變量y與自變量x之間無確定性解析表達(dá)式時(shí)，建立預(yù)測模型的方法有神經(jīng)網(wǎng)絡(luò)和支持向量機(jī)等智能學(xué)習(xí)方法。支持向量機(jī)由Vapnik首先提出，可用于模式分類和非線性回歸，它的理論基礎(chǔ)是統(tǒng)計(jì)學(xué)習(xí)理論，具有解決問題的有效性高、計(jì)算方法簡單、運(yùn)行速度快及理論完善等優(yōu)點(diǎn)。在非線性回歸中核函數(shù)K的選擇至關(guān)重要，其種類主要有線性核函數(shù)、多項(xiàng)式核函數(shù)、徑向基核函數(shù)（RBF）和sigmoid核函數(shù)。徑向基核函數(shù)作為局部性核函數(shù)，只對測試點(diǎn)附近領(lǐng)域的數(shù)據(jù)點(diǎn)有影響，故適合作為通用模式的核函數(shù)[11]，本文選擇徑向基核函數(shù)進(jìn)行非線性回歸分析。