游　健

(云南省水利水電勘測設(shè)計研究院， 昆明　650000)

?

基于改進(jìn)遺傳算法優(yōu)化的混凝土壩壩基揚壓力統(tǒng)計模型研究

游健

(云南省水利水電勘測設(shè)計研究院， 昆明650000)

摘要：混凝土壩壩基揚壓力是外界復(fù)合環(huán)境因素綜合作用的結(jié)果，如何建立準(zhǔn)確反應(yīng)大壩實際滲流性態(tài)的統(tǒng)計模型對于掌握大壩安全性態(tài)具有十分重要的意義．本文基于外界環(huán)境量對混凝土壩壩基揚壓力的影響過程分析，將前期庫水位及降雨對壩基揚壓力的非線性影響過程視為“正態(tài)分布”過程，并采用改進(jìn)遺傳算法對相應(yīng)滯后天數(shù)和影響天數(shù)進(jìn)行智能尋優(yōu)，以提高壩基揚壓力統(tǒng)計模型準(zhǔn)確性．最后，以某混凝土壩為例，采用本文方法對其壩基揚壓力實測進(jìn)行擬合和預(yù)測，通過與傳統(tǒng)統(tǒng)計模型對比，表明本文方法具有更高的擬合和預(yù)測精度，具有一定的科學(xué)和實踐應(yīng)用價值．

關(guān)鍵詞：混凝土壩；壩基揚壓力；統(tǒng)計模型；改進(jìn)遺傳算法

0引言

混凝土壩壩基揚壓力安全監(jiān)測是大壩安全監(jiān)測的重要組成部分，研究表明壩高100m左右的重力壩，壩基面上作用的揚壓力大約為壩體重量的20%左右[1]．因此，如何建立反應(yīng)混凝土大壩壩基實際滲流性態(tài)的揚壓力統(tǒng)計模型對于掌握大壩安全工作性態(tài)具有十分重要的意義．姜諳男等基于支持向量機結(jié)構(gòu)風(fēng)險最小化原理，采用粒子群算法快速優(yōu)化支持向量機模型參數(shù)，對大壩滲流監(jiān)測資料進(jìn)行非線性擬合，建立了基于PSO-SVM的人壩滲流監(jiān)測的時間序列非線性預(yù)報模型[2]．為了消除統(tǒng)計模型中各因子相關(guān)性對分析結(jié)果的影響，田偉等將主成分分析法和神經(jīng)網(wǎng)絡(luò)方法相結(jié)合，建立了大壩滲流監(jiān)測數(shù)據(jù)的主成分神經(jīng)網(wǎng)絡(luò)模型[3]．李波等采用最小二乘回歸法對大壩滲流各因子進(jìn)行篩選，提取對因變量影響強的成分，建立了基于PLS-LSSVM的大壩滲流監(jiān)控模型，有效克服了變量間的多重相關(guān)性問題[4]．梅一韜等利用物元可拓理論與模糊數(shù)學(xué)理論構(gòu)造了滲流性態(tài)評價中的模糊物元，引入嫡權(quán)法確定大壩滲流指標(biāo)體系的權(quán)重系數(shù)，建立了基于嫡權(quán)的模糊可拓評價模型[5]．傳統(tǒng)混凝土壩基揚壓力統(tǒng)計模型將壩基揚壓力視為與前期平均庫水位及降雨量呈線性相關(guān)關(guān)系，然而大壩安全資料分析表明上游庫水位及降雨對混凝土壩壩基揚壓力的影響過程是一個先增長后下降的過程，而不是一個“平均”過程，因此需要進(jìn)一步研究混凝土重力壩壩基揚壓力統(tǒng)計模型的因子選擇以提高模型準(zhǔn)確性．

本文基于上游庫水位及降雨對壩基揚壓力的非線性影響過程，采用改進(jìn)遺傳算法以相應(yīng)統(tǒng)計模型復(fù)相關(guān)系數(shù)為適應(yīng)度值，對相應(yīng)滯后天數(shù)和影響天數(shù)進(jìn)行智能優(yōu)化，從而提高統(tǒng)計模型精度．最后為了驗證本文方法的有效性，以某混凝土重力壩為例，采用上述方法對其多年壩基揚壓力進(jìn)行擬合，通過與傳統(tǒng)統(tǒng)計模型對比分析，表明本文方法具有更高的擬合和預(yù)測精度．

1基本理論

基于混凝土壩實測資料分析，本文混凝土重力壩壩基揚壓力統(tǒng)計模型采用如下形式：

(1)

式中，P為壩基揚壓力，PH為水壓分量，PR為降雨分量，Pθ為時效分量，PT是溫度分量．

1)水位分量PH．壩基揚壓力是在庫水位變動情況下非穩(wěn)定滲流的瞬時狀態(tài)結(jié)果，任意時刻的壩基揚壓力是此時刻以前庫水位連續(xù)變化影響的瞬態(tài)結(jié)果．傳統(tǒng)統(tǒng)計模型將前期水位的平均值作為因子進(jìn)行考慮(如圖1所示)，然而，監(jiān)測資料分析表明，庫水位對壩基揚壓力的影響是一個先上升后下降的過程，基本服從“正態(tài)”分布過程(如圖2所示)[6-7]．不同于傳統(tǒng)統(tǒng)計模型，本文考慮庫水位對壩基滲流非線性影響的水壓分量采用如下形式：

(2)

式中，aj為回歸系數(shù)；x1為滯后天數(shù)；x2為影響天數(shù)；Hi(t)為上游水位；Hum為有效水位．

圖1　前x3天平均水位　　　圖2　正態(tài)分布影響曲線(P-影響權(quán))　　　　　　　　　(P-影響權(quán))

其中，x1，x2需要通過試算求得，在求得準(zhǔn)確滯后天數(shù)和影響天數(shù)后，即可提高相應(yīng)統(tǒng)計模型擬合和預(yù)測精度．

2)降雨分量PR．降雨與大壩滲流的關(guān)系比較復(fù)雜，它與入滲條件、雨型、地質(zhì)條件等因素有關(guān)，降雨量與大壩滲流呈現(xiàn)明顯的非線性關(guān)系[8-9]．基于此，采用指數(shù)關(guān)系法，把降雨量做2/5指數(shù)變化，變換后降雨量變化比較平穩(wěn)，綜合考慮降雨對大壩滲流的滯后影響過程，本文采用壩基揚壓力降雨分量如下：

(3)

式中，bj為降雨分量的回歸系數(shù)；x3為降雨分量的滯后天數(shù)；x4為降雨影響權(quán)正態(tài)分布標(biāo)準(zhǔn)差(影響天數(shù))；R(t)為t時刻的單位時段降雨量；Rum為有效降雨量．其中，x3和x4需要通過試算求得．

3)溫度分量PT．混凝土壩基巖裂隙變化受基巖溫度影響顯著，而基巖溫度變化基本呈年周期性變化規(guī)律，本文采用周期函數(shù)來描述溫度分量：

(4)

式中，c1i，c2i為回歸系數(shù)．

4)時效分量Pθ．時效分量是混凝土重力壩壩基揚壓力的重要組成部分，也是評價大壩滲流穩(wěn)定性的重要依據(jù)[10]，將壩體測壓管水位時效分量表示為：

(5)

式中，d1，d2為回歸系數(shù)；θ為監(jiān)測日自初始監(jiān)測日天數(shù)除以100．

因此，本文定義的新型混凝土壩壩基揚壓力統(tǒng)計模型如下：

(6)

2雙鏈量子遺傳算法

傳統(tǒng)遺傳算法(GA)采用二進(jìn)制編碼，該編碼方式存在隨機性，容易造成算法收斂速度慢、容易陷入局部極值等問題[11-13]．基于此，本文采用雙鏈量子遺傳算法對相應(yīng)統(tǒng)計模型的滯后天數(shù)和影響天數(shù)進(jìn)行智能優(yōu)化，其流程如圖3所示．

圖3　雙鏈量子遺傳算法流程圖

該算法采用量子位概率幅值對染色體進(jìn)行編碼，每條染色體分別包括兩個優(yōu)化解，通過量子旋轉(zhuǎn)門與量子非門來增加種群多樣性，從而提高遺傳算法的收斂速度與精度[14-15]．

一個量子位可能處于0態(tài)、1態(tài)、以及|0〉和|1〉之間的任意疊加態(tài)．可以表示為：

(7)

式中，α，β為相應(yīng)量子位的概率幅值．|α|2為量子態(tài)達(dá)到|0〉態(tài)時相應(yīng)的概率，|β〉2為量子態(tài)達(dá)到|1〉態(tài)的概率，并且滿足歸一化條件

(8)

(9)

式中，θij=2π×rnd；rnd∈(0，1)為隨機數(shù)；i=1，2，…，m，j=1，2，…，n，其中m為種群規(guī)模，n為空間維數(shù)．

(10)

即可表示為：

(11)

式中，i=1，2，…，m，j=1，2，…，n．

量子門更新過程可以表示為：

(12)

(13)

所以，可以得出：

相應(yīng)量子轉(zhuǎn)角定義如下：

(14)

(15)

(16)

為了避免算法陷入局部極值，本文同時引入量子非門作為變異算子來增加種群多樣性，相應(yīng)量子染色體變異操作如下：

(17)

式中，[cosθijsinθij]T為編號i量子染色體|pi〉的第j個基因位．

3工程實例

3.1工程簡介

某水電站位于福建省閩清縣境內(nèi)的閩江干流上，為了監(jiān)測大壩滲流狀況，掌握揚壓力的變化規(guī)律，沿壩軸線方向設(shè)置一條縱向觀測斷面．在帷幕下游側(cè)共設(shè)置了38個觀測孔，沿縱向均勻分布．另外，根據(jù)壩基地質(zhì)情況，在11號、17號、25號、32號壩段的橫向觀測廊道中，分別設(shè)4個揚壓力觀測孔[16]．

3.2壩基揚壓力分析

為驗證本文雙鏈量子遺傳算法有效性，本文同時采用傳統(tǒng)遺傳算法進(jìn)行計算，兩種方法最優(yōu)個體適應(yīng)度值進(jìn)化過程如圖4所示．由圖4可知，雙鏈量子遺傳算法最終復(fù)相關(guān)系數(shù)為0.961，而傳統(tǒng)遺傳算法最終復(fù)相關(guān)系數(shù)為0.935，雙鏈量子遺傳算法具有更高的計算效率和收斂精度．

圖4　種群最優(yōu)適應(yīng)度值發(fā)展曲線

為驗證本文方法有效性，本文同時采用多元非線性回歸統(tǒng)計模型[1]對該測點測值進(jìn)行擬合，相應(yīng)統(tǒng)計模型為：

(18)

式中，c0為待求常數(shù)項，H(t)，H(t-τn)分別為時刻t和計算時間前τn時間的水頭，τn為滯后時間，w(t)，w(t-i)分別為計算當(dāng)日的降水量和計算時刻前i天內(nèi)前期降水量的累加值，cij，αi，βi，a，b，di分別為待定系數(shù)．

通過采用多元非線性統(tǒng)計回歸法，剔除影響因素不明顯的因子，可得關(guān)于滲流監(jiān)測效應(yīng)量統(tǒng)計模型表達(dá)式，該方法可以克服傳統(tǒng)統(tǒng)計模型只能針對單測點監(jiān)測的缺點，相應(yīng)兩種方法擬合結(jié)果及兩種方法擬合殘差如圖5～6所示．

圖5　傳統(tǒng)方法擬合結(jié)果

圖6　傳統(tǒng)方法擬合結(jié)果

由圖可知，本文方法擬合結(jié)果與7號壩段壩蹱位置揚壓力實測值接近，能夠準(zhǔn)確反映該測點壩基揚壓力的變化規(guī)律．本文方法所得復(fù)相關(guān)系數(shù)為0.987，遠(yuǎn)大于多元非線性回歸統(tǒng)計模型擬合結(jié)果0.945，具有更高的擬合精度．因此，本文方法能夠準(zhǔn)確反映前期庫水位及降雨對混凝土壩基揚壓力的非線性影響．

為了驗證本文方法的預(yù)測精度，采用上述兩種模型對該測點2008年10月20日至2008年10月31日之間揚壓力測值進(jìn)行預(yù)測，并與實測值對比，如圖7所示．

圖7　壩基揚壓力預(yù)測結(jié)果對比圖

由圖可知，本文方法預(yù)測結(jié)果相比多元非線性統(tǒng)計模型具有更高的預(yù)測精度，此外，由于本文模型考慮了前期庫水位及降雨對壩基揚壓力測值的非線性影響，因此能夠反映出壩基揚壓力測值的非線性變化．

4結(jié)論

本文考慮上游庫水位及降雨對混凝土壩壩基揚壓力的非線性影響，將前期庫水位及降雨對壩基揚壓力的影響視為“正態(tài)分布”過程，并采用雙鏈量子遺傳算法對其進(jìn)行智能優(yōu)化，通過傳統(tǒng)統(tǒng)計模型計算結(jié)果對比分析，可以得出如下結(jié)論：

1)雙鏈量子遺傳算法綜合量子運算和遺傳算法基本原理，采用量子概率幅值對染色體進(jìn)行編碼，利用量子門及量子非門對染色體進(jìn)行實時更新，避免了傳統(tǒng)遺傳算法容易產(chǎn)生局部極值的現(xiàn)象，具有更高的計算效率和精度．

2)通過將上游庫水位及降雨對壩基揚壓力的影響視為“正態(tài)分布”過程，并采用改進(jìn)遺傳算法對相應(yīng)滯后天數(shù)和影響天數(shù)進(jìn)行智能尋優(yōu)，相比傳統(tǒng)統(tǒng)計模型具有更高的擬合精度．

3)由于本文模型能夠準(zhǔn)確反映前期庫水位及降雨對壩基揚壓力的非線性影響，因此本文模型相比傳統(tǒng)統(tǒng)計模型能夠準(zhǔn)確預(yù)測壩基揚壓力的變化趨勢，具有更高的預(yù)測精度．

參考文獻(xiàn)：

[1]顧沖時，吳中如．大壩與壩基安全監(jiān)控理論和方法及其應(yīng)用[M]．南京：河海大學(xué)出版社，2006．

[2]姜諳男，梁冰．基于PSO-SVM的大壩滲流監(jiān)測時間序列非線性預(yù)報模型[J]．水利學(xué)報，2006，37(3)：331-335．

[3]田偉，魏光輝，高強．基于主成分分析與BP神經(jīng)網(wǎng)絡(luò)模型的大壩滲流監(jiān)測資料分析[J]．大壩與安全，2009(5)：29-31．

[4]李波，顧沖時，李智錄，等．基于偏最小二乘回歸和最小二乘支持向量機的大壩滲流監(jiān)控模型[J]．水利學(xué)報，2008，39(12)：1390-1394．

[5]梅一韜，仲云飛．基于熵權(quán)的大壩滲流性態(tài)模糊可拓評價模型[J]．水電能源科學(xué)，2011，29(8)：58-61．

[6]劉耀儒，楊強，覃振朝．基于統(tǒng)計模型的裂隙巖體滲流場的并行數(shù)值模擬[J]．巖石力學(xué)與工程學(xué)報，2008(4)：：736-742．

[7]洪云，鄭東鍵．古田溪一級大壩14號壩段以左的滲壓確定性模型[J]．水電自動化與大壩監(jiān)測，1996(3)：30-31．

[8]顧沖時，李波，徐寶松，等．碾壓混凝土壩滲漏安全監(jiān)控模型建模方法研究[J]．河海大學(xué)學(xué)報：自然科學(xué)版，2008，36(5)：579-585．

[9]胡江，鄭鵬翔．基于滯后效應(yīng)和防洪調(diào)度的大壩滲流預(yù)測模型研究[J]．三峽大學(xué)學(xué)報：自然科學(xué)版，2008，30(6)：16-19．

[10] 畢奇龍，劉星，徐波，等．因子分析法在繞壩滲流監(jiān)測資料分析中的應(yīng)用[J]．水電能源科學(xué)，2010，28(5)：62-65．

[11] 楊建文，李志鵬，喻哲欣．基于改進(jìn)遺傳算法的泵站優(yōu)化運行[J]．中國水利，2015(8)：33-35．

[12] 張尚然，湯亞芳，林俐．基于改進(jìn)遺傳算法的無功優(yōu)化方法的研究[J]．中小企業(yè)管理與科技，2015(7)：162-163．

[13] 邵改革，盧小平，杜耀剛，等．基于優(yōu)化模糊遺傳算法的地理信息個性化推送方法[J]．測繪工程，2015，24(7)：20-23．

[14] 夏俊，凌培亮，虞麗娟，等．基于量子遺傳算法的無線傳感網(wǎng)絡(luò)路由優(yōu)化[J]．同濟大學(xué)學(xué)報：自然科學(xué)版，2015，43(7)：1097-1103．

[15] 沙林秀，賀昱曜，陳延偉．一種變步長雙鏈量子遺傳算法[J]．計算機工程與應(yīng)用，2012(20)：59-63．

[16] 李延芳，王玉潔，趙花城．水口水電站大壩原型觀測設(shè)計[J]．水力發(fā)電，1994(6)：23-25．

[責(zé)任編輯王康平]

DOI：10.13393/j.cnki.issn.1672-948X.2016.03.002

收稿日期：2015-11-17

基金項目：國家自然科學(xué)基金“靜動結(jié)合的高拱壩健康性態(tài)監(jiān)測和診斷方法研究”

通信作者：游健(1962-)，女，高級工程師，主要研究方向為水工設(shè)計．E-mail：yannan_youjian@126.com

中圖分類號：TV64

文獻(xiàn)標(biāo)識碼：A

文章編號：1672-948X(2016)03-0006-05

Research on Statistical Model of Uplift Pressure Based on Improved Genetic Algorithm

You Jian

(Yunnan Provincial Conservancy & Hydropower Survey, Design Institute, Kunmin 650000， China)

AbstractThe variation of uplift pressure on concrete dam foundation is effected by the compound external environmental factors. How to establish statistical model which accurately reflects the real dam seepage state is quite important to grasp the safety state of a dam. Based on analyzing the influential process of ambient variables, the paper considers the nonlinear influential process as “Normal” process and adopts the improved genetic algorithm to optimize the lagged days and the influential days to raise the accuracy of the statistical model. Finally, taking a certain concrete dam for example, the recorded series of the uplift pressure is fitted and predicted with the proposed method. By comparing with the traditional model, it is verified that the proposed method possesses higher fitting and prediction accuracy and is of certain scientific and practical value.

Keywordsconcrete dam;foundation uplift pressure;statistical model;improved genetic algorithm