張返成

(吉安市螺灘水利水電中心，江西 吉安 343060)

滲透系數(shù)可反應(yīng)心墻堆石壩內(nèi)部的滲流情況，該因素對(duì)大壩的力學(xué)性能及穩(wěn)定性影響顯著，近年來，許多專家學(xué)者針對(duì)心墻內(nèi)部滲流情況開展相關(guān)研究。

張超萍等人[1]基于監(jiān)測(cè)數(shù)據(jù)，對(duì)某心墻堆石壩的滲流情況進(jìn)行監(jiān)測(cè)，分析不同防滲效果對(duì)其滲流情況的影響，結(jié)果表明，防滲墻的防滲效果最好。許云鵬等人[2]基于土壤水分特征曲線，建立心墻堆石壩的滲流模型，分析該模型預(yù)測(cè)其滲漏系數(shù)的精度，結(jié)果表明，該模型可準(zhǔn)確反應(yīng)大壩內(nèi)部的滲流情況，準(zhǔn)確性較高。丁艷輝等人[3]對(duì)某水庫(kù)大壩的滲流及變形規(guī)律進(jìn)行監(jiān)測(cè)，分析大壩頂部發(fā)生不均勻沉降的原因，結(jié)果表明，大壩內(nèi)部滲流會(huì)導(dǎo)致壩頂變形，但是大壩總體穩(wěn)定性及安全性較高。王芳等人[4]結(jié)合正交試驗(yàn)與數(shù)值模擬，分析庫(kù)水位變化條件下，某心墻堆石壩的變形及滲流規(guī)律，結(jié)果表明，大壩變形與心墻內(nèi)部的滲流情況有關(guān)，且與滲透系數(shù)呈正比。楊啟貴等人[5]以某瀝青混凝土心墻堆石壩為研究對(duì)象，對(duì)其壩體的滲流及穩(wěn)定性進(jìn)行計(jì)算，并分析動(dòng)、靜荷載下壩體的變形情況，并依據(jù)計(jì)算結(jié)果，對(duì)大壩進(jìn)行設(shè)計(jì)，為后續(xù)滲流控制提供相關(guān)參考依據(jù)。

本研究以某水利工程的心墻堆石壩為研究對(duì)象，基于鯨魚優(yōu)化算法，構(gòu)建滲透系數(shù)預(yù)測(cè)模型，分析各種因素對(duì)滲透系數(shù)的影響。

1 工程概況

本研究以某水利工程為研究對(duì)象，工程由大壩樞紐、供水工程兩部分組成，大壩為面板堆石壩，防滲體為礫石土直心墻，最大壩高63.0m，壩頂高程1668.0m，壩軸線長(zhǎng)121.84m。流域地處第二地形階梯帶，地貌屬高原中山溝谷地貌類型，地下水由大氣降雨及地表徑流的補(bǔ)給，通過脈狀流和巖溶管道流形式往低處排泄，最終呈分散狀或集中形式排出地表補(bǔ)給下游河流。庫(kù)區(qū)為溝谷地形，出露地層為峨眉山玄武巖組(P2β)的碎屑巖地層及第四系覆蓋層，該地區(qū)不良物理地質(zhì)現(xiàn)象總體發(fā)育較弱，主要表現(xiàn)為風(fēng)化、岸坡卸荷、第四系堆積體等，分布在河流、溝谷兩岸緩坡及坡腳地帶。

2 滲透系數(shù)預(yù)測(cè)模型構(gòu)建

本研究基于鯨魚優(yōu)化算法，構(gòu)建滲透系數(shù)預(yù)測(cè)模型，分析各種因素對(duì)滲透系數(shù)(k)的影響，其目標(biāo)函數(shù)如式(1)所示。

k=f(Cu，Cc，D5，D0.075，w，N0，N1，N2，H)

(1)

式中，Cu—不均勻系數(shù)；Cc—曲率系數(shù)；D5—粒徑大于5mm顆粒質(zhì)量分?jǐn)?shù)，%；、D0.075—粒徑小于0.075mm顆粒質(zhì)量分?jǐn)?shù)，%；w—含水率，%；Ni—碾遍次數(shù)，次；H—壓實(shí)厚度，m。

采用最大信息系數(shù)(MIC)分析式(1)中參數(shù)對(duì)滲透系數(shù)的影響，其計(jì)算公式如式(2)所示。

(2)

式中，p—聯(lián)合概率；I—互信息。

為分析不同參數(shù)對(duì)心墻礫石土滲透系數(shù)的影響，分別選取9組數(shù)據(jù)，對(duì)滲透系數(shù)進(jìn)行預(yù)測(cè)，其輸入?yún)?shù)見表1。

表1 輸入?yún)?shù)

3 結(jié)果分析

BPNN模型和WOA-SVM模型是鯨魚優(yōu)化算法中兩種常見的訓(xùn)練算法，為選取合適的算法模型對(duì)滲透系數(shù)預(yù)測(cè)模型進(jìn)行訓(xùn)練，采用BPNN模型和WOA-SVM模型對(duì)滲透系數(shù)預(yù)測(cè)模型進(jìn)行訓(xùn)練，并將其預(yù)測(cè)值與實(shí)際值進(jìn)行對(duì)比，采用BPNN模型得出的測(cè)試樣本編號(hào)-滲透系數(shù)曲線如圖1所示。

圖1 BPNN模型的測(cè)試樣本編號(hào)-滲透系數(shù)曲線

由圖1可知，采用BPNN模型得出的滲透系數(shù)波動(dòng)范圍范圍較大，其值集中于0～0.8×10-5cm/s間，滲透系數(shù)實(shí)際值的波動(dòng)范圍在0～1.0×10-5cm/s間，二者間的最大誤差為0.7×10-5cm/s，最小誤差小于0.01×10-5cm/s，說明采用該模型對(duì)滲透系數(shù)預(yù)測(cè)模型進(jìn)行訓(xùn)練的效果較好，但預(yù)測(cè)值與實(shí)際值間仍存在一定的差異，預(yù)測(cè)精度存在改進(jìn)的空間。

采用WOA-SVM模型得出的測(cè)試樣本編號(hào)-滲透系數(shù)曲線如圖2所示。

圖2 WOA-SVM模型的測(cè)試樣本編號(hào)-滲透系數(shù)曲線

由圖2可知，采用WOA-SVM模型得出的滲透系數(shù)預(yù)測(cè)值波動(dòng)范圍較小，其滲透系數(shù)集中于0.2×10-5cm/s～0.4×10-5cm/s，采用WOA-SVM模型得出的滲透系數(shù)大多大于滲透系數(shù)實(shí)際值，二者間的滲透系數(shù)最大差值為0.4×10-5cm/s，滲透系數(shù)最小差值小于0.02×10-5cm/s，說明采用WOA-SVM模型對(duì)滲透系數(shù)進(jìn)行預(yù)測(cè)的準(zhǔn)確性良好。綜合以上分析可得，采用WOA-SVM模型訓(xùn)練得出的滲透系數(shù)較為穩(wěn)定，其波動(dòng)范圍較小，但是其預(yù)測(cè)精度劣于BPNN模型；采用BPNN模型訓(xùn)練得出的滲透系數(shù)波動(dòng)性較大，但是其滲透系數(shù)預(yù)測(cè)結(jié)果與實(shí)際值間的差值較??；以上兩種模型對(duì)于心墻礫石土的滲透系數(shù)的預(yù)測(cè)效果及精度良好，且各有優(yōu)勢(shì)，但是其預(yù)測(cè)精度仍存在可以改進(jìn)的空間，為進(jìn)一步提高滲透系數(shù)預(yù)測(cè)效果和預(yù)測(cè)精度，可以結(jié)合BPNN模型和WOA-SVM模型對(duì)滲透系數(shù)預(yù)測(cè)模型進(jìn)行訓(xùn)練。

根據(jù)以上分析可得，結(jié)合BPNN模型和WOA-SVM模型可提高滲透系數(shù)預(yù)測(cè)的準(zhǔn)確性，為進(jìn)一步提高滲透系數(shù)預(yù)測(cè)效果和預(yù)測(cè)精度，對(duì)比分析BPNN模型和WOA-SVM模型結(jié)合得出的滲透系數(shù)預(yù)測(cè)結(jié)果與滲透系數(shù)實(shí)際值間的差異，測(cè)試樣本編號(hào)-滲透系數(shù)曲線如圖3所示。

由圖3可知，采用BPNN-WOA-SVM模型得出的滲透系數(shù)預(yù)測(cè)值與實(shí)際值間的差值較小，且二者間的測(cè)試樣本編號(hào)-滲透系數(shù)曲線變化趨勢(shì)具有一致性，預(yù)測(cè)值與實(shí)際值間的最大差值為0.15×10-5cm/s，對(duì)比僅使用BPNN模型和WOA-SVM模型訓(xùn)練的滲透系數(shù)預(yù)測(cè)結(jié)果可得，結(jié)合2種模型進(jìn)行預(yù)測(cè)誤差減小了0.25×10-5cm/s～0.55×10-5cm/s，且預(yù)測(cè)結(jié)果較為穩(wěn)定，BPNN-WOA-SVM模型對(duì)滲透系數(shù)預(yù)測(cè)的準(zhǔn)確性提升效果顯著。

為直觀反映BPNN-WOA-SVM模型對(duì)滲透系數(shù)的預(yù)測(cè)效果，對(duì)采用不同模型預(yù)測(cè)的結(jié)果進(jìn)行誤差分析，誤差分析結(jié)果見表2。

表2 誤差分析結(jié)果

由表2可知，采用BPNN-WOA-SVM模型預(yù)測(cè)得出的滲透系數(shù)平均值與實(shí)測(cè)值間的差值最小，二者間的差值僅為0.5×10-7cm/s。通過不同模型的誤差分析可得，BPNN-WOA-SVM模型的均方誤差與平均絕對(duì)誤差有最小值，說明采用該模型的預(yù)測(cè)準(zhǔn)確性較高。采用3種模型得出的相對(duì)分析誤差均大于2，說明以上模型的穩(wěn)定性好，能夠達(dá)到實(shí)際應(yīng)用標(biāo)準(zhǔn)，其中，BPNN-WOA-SVM模型的相對(duì)分析誤差最大，該模型的預(yù)測(cè)效果最好。

根據(jù)以上分析可得，采用BPNN-WOA-SVM模型對(duì)心墻礫石土的滲透系數(shù)進(jìn)行預(yù)測(cè)的效果較好，為分析碾壓施工參數(shù)對(duì)其滲透系數(shù)的影響，以表1所示的參數(shù)為研究對(duì)象，對(duì)滲透系數(shù)進(jìn)行預(yù)測(cè)，將預(yù)測(cè)結(jié)果與不考慮碾壓施工參數(shù)的預(yù)測(cè)值進(jìn)行對(duì)比，其樣本序號(hào)-滲透系數(shù)曲線如圖4所示。

圖4 樣本序號(hào)-滲透系數(shù)曲線

由圖4可知，不同碾壓參數(shù)下的心墻礫石土滲透系數(shù)具有一定的差異性，其中，樣本5的預(yù)測(cè)值最大，當(dāng)考慮碾壓參數(shù)預(yù)測(cè)值時(shí)，滲透系數(shù)為0.61×10-5cm/s，樣本4的預(yù)測(cè)值最小，當(dāng)考慮碾壓參數(shù)預(yù)測(cè)值時(shí)，其滲透系數(shù)為0.17×10-5cm/s。對(duì)比是否考慮碾壓參數(shù)下的心墻礫石土滲透系數(shù)可得，二者得出的滲透系數(shù)具有一定的差異性，其中，考慮碾壓參數(shù)的滲透系數(shù)預(yù)測(cè)值與實(shí)測(cè)值間的差值較小，說明考慮碾壓參數(shù)對(duì)心墻礫石土滲透系數(shù)進(jìn)行預(yù)測(cè)更符合實(shí)際工程的情況，預(yù)測(cè)結(jié)果更加準(zhǔn)確。

分析碾壓參數(shù)對(duì)各倉(cāng)面的滲透系數(shù)影響，見表3。

表3 碾壓參數(shù)影響結(jié)果

由表3可知，考慮碾壓參數(shù)的滲透系數(shù)與實(shí)測(cè)值間的差值較小，其中，倉(cāng)面1與倉(cāng)面3的預(yù)測(cè)精度最高，二者間的滲透系數(shù)平均值差值均為0.002×10-5cm/s，不考慮碾壓參數(shù)的預(yù)測(cè)精度較差，與實(shí)測(cè)值間的滲透系數(shù)差值均大于0.2×10-5cm/s，說明碾壓參數(shù)對(duì)心墻礫石土滲透系數(shù)的影響較大，為提高預(yù)測(cè)結(jié)果的準(zhǔn)確性，應(yīng)考慮碾壓施工參數(shù)對(duì)滲透系數(shù)的影響[6]。

以表1中的參數(shù)為研究對(duì)象，對(duì)考慮碾壓參數(shù)下的心墻礫石土滲透系數(shù)進(jìn)行預(yù)測(cè)，對(duì)比不同模型下的預(yù)測(cè)結(jié)果，其樣本序號(hào)-滲透系數(shù)曲線如圖5所示。

圖5 樣本序號(hào)-滲透系數(shù)曲線

由圖5可知，采用WOA-SVM模型得出的預(yù)測(cè)結(jié)果與實(shí)測(cè)值間的差值最大，說明采用該模型對(duì)心墻礫石土滲透系數(shù)進(jìn)行預(yù)測(cè)的精度最低，采用BPNN-WOA-SVM模型得出的樣本序號(hào)-滲透系數(shù)曲線與實(shí)測(cè)值曲線的吻合度最高，說明采用組合模型可顯著提升心墻礫石土的滲透系數(shù)的預(yù)測(cè)精度。

4 結(jié)論

本文以某心墻堆石壩為研究對(duì)象，基于鯨魚優(yōu)化算法，構(gòu)建滲透系數(shù)預(yù)測(cè)模型，分析各種因素對(duì)滲透系數(shù)的影響，對(duì)比僅使用BPNN模型和WOA-SVM模型訓(xùn)練的滲透系數(shù)預(yù)測(cè)結(jié)果表明，單一模型與組合模型的穩(wěn)定性好，能夠達(dá)到實(shí)際應(yīng)用標(biāo)準(zhǔn)，采用WOA-SVM模型得出的預(yù)測(cè)結(jié)果與實(shí)測(cè)值間的差值最大，精度最低，采用BPNN-WOA-SVM模型得出的樣本序號(hào)-滲透系數(shù)曲線與實(shí)測(cè)值曲線的吻合度最高，模型預(yù)測(cè)結(jié)果較為穩(wěn)定，得出的滲透系數(shù)預(yù)測(cè)值與實(shí)際值間的差值較小，對(duì)滲透系數(shù)預(yù)測(cè)的準(zhǔn)確性提升效果顯著。施工參數(shù)對(duì)心墻礫石土滲透系數(shù)的影響較大，為提高預(yù)測(cè)結(jié)果的準(zhǔn)確性，應(yīng)考慮碾壓施工參數(shù)對(duì)滲透系數(shù)的影響。