周 密
( 中國南方電網(wǎng)調(diào)峰調(diào)頻發(fā)電公司 天生橋水力發(fā)電總廠,貴州 興義562400)
我國水能資源豐富,水資源分布卻極不平衡,為了開發(fā)利用水能資源和解決垮流域供水問題,我國自建國以來興建了將近9 萬座大壩,這些大壩不僅在防洪、灌溉、水力發(fā)電、水路航運和城市供水等方面發(fā)揮著重要作用,同時對流域生態(tài)環(huán)境改造等方面也有著重要影響。由于大壩在運行過程中,受到水壓力、溫度、當?shù)丨h(huán)境等多種因素影響會產(chǎn)生變形,該變形在時間上具有明顯的時效性,是一種復(fù)雜的非線性變化。一旦大壩變形超過了設(shè)定范圍,將直接影響到大壩自身安全,進而對周圍工農(nóng)業(yè)生產(chǎn)和人類生活活動造成極大的威脅。因此,開展大壩自動檢測與趨勢預(yù)測,一直都是水利水文工作者的研究熱點之一[1-2]。
近年來,由于自動化水平的提高以及逐步回歸分析、神經(jīng)網(wǎng)絡(luò)、小波分析等函數(shù)逼近與預(yù)測算法的發(fā)展,大壩變形預(yù)測也得到了較快的發(fā)展。比如黃世秀等[3]針對港口灣大壩實測數(shù)據(jù),采用BP 神經(jīng)網(wǎng)絡(luò)模型進行大壩變形預(yù)測,獲得了有效的預(yù)測值,進一步研究了大壩變化規(guī)律; 朱金海等[4]結(jié)合大壩沉降監(jiān)測的具體實例,應(yīng)用時間序列分析理論和方法對大壩變形進行預(yù)報,該預(yù)測方法對短期預(yù)報效果更佳,對中長期預(yù)報效果不佳;李明然等[5]采用改進的灰色馬爾科夫模型進行大壩預(yù)測,該研究表明大壩變形監(jiān)測數(shù)據(jù)呈現(xiàn)隨機性和波動性變化特點,仿真方法的適當與否對大壩變形預(yù)測至關(guān)重要。
由上述分析可知,由于大壩變形影響參數(shù)較多,變化周期不一致,比如,溫度變化是日星期型的,水位的變化則以季節(jié)為周期,這都極大地影響著大壩變形預(yù)測的精度;同時,大部分的預(yù)測分析方法或多或少都有著一些缺陷,比如,逐步回歸分析模型的精度與觀測樣本數(shù)目有關(guān); 人工神經(jīng)網(wǎng)絡(luò)預(yù)測方法收斂速度又比較慢等等,這些又影響這大壩變形預(yù)測分析的經(jīng)濟性。
為此,本文采用多元線性回歸分析方法對大壩變形影響因素進行深入分析,在此基礎(chǔ)上引入神經(jīng)網(wǎng)絡(luò)模型,構(gòu)成大壩變形預(yù)測混合模型,有效提高了大壩變形預(yù)測的精度。
一般認為,影響大壩變形的3 類因素為: 水壓因素H、溫度因素T 和時效因素S。其中,水壓因素與水庫水位密切相關(guān),水庫水位與大壩變形之間呈復(fù)雜非線性關(guān)系,為此,可將水壓因素分解成H、H2、H3、H4、H5等5個分量,即大壩變形與水位及水位的2、3、4、5 次方有關(guān); 溫度因素又可以分為水溫和氣溫,由于水溫變化較小可以省略,而氣溫變化對大壩的影響又具有滯后效應(yīng),為此,選取溫度T 和T1 ( 為上一個記錄時刻的溫度) 作為溫度影響因素。再考慮到時效因素,則影響大壩變形參數(shù)可分為3 類,共8個變量。
多元線性回歸分析是將變量一個一個引入,進而通過檢驗各個變量的偏回歸平方和是否顯著,以決定是否保留該變量,這樣經(jīng)過顯著性判斷之后,所保留下來的變量就是對大壩變形影響較顯著的變量。
假設(shè)隨機變量y 與p個自變量x1,x2,…xp之間存在著線性相關(guān)關(guān)系,實際樣本量為n,則其n 次觀測值可寫為如下形式:
式中: β1,β2,…βp是未知參數(shù); x1,x2,…xp是p個可以精確測量并可控制的一般變量; ε1,ε2,…εn是隨機誤差。假定εi是相互獨立且服從同一正態(tài)分布N(0,σ) 的隨機變量。
若將方程組(1) 用矩陣表示,則有:
多元線性回歸分析的首要任務(wù)就是通過尋求β 的估計值b,建立多元線性回歸方程見式3,來描述多元線性模型見式1:
大壩變形神經(jīng)網(wǎng)絡(luò)預(yù)測結(jié)構(gòu),見圖1 所示,大壩變形神經(jīng)網(wǎng)絡(luò)預(yù)測模型由輸入層、中間層和輸出層組成,其中輸入層節(jié)點數(shù)n 為多元線性分析后保留下的影響參量個數(shù),輸出節(jié)點只有一個,即大壩變形值。中間層節(jié)點數(shù)m個,一般m 可取6 ~15個。在模型預(yù)測的過程中,將大壩影響參量分別送入輸入層的n個節(jié)點,而中間層每個節(jié)點的輸出值為輸入層所有節(jié)點對其輸出的累加( 利用sigmoid 函數(shù)計算) ,整個模型的輸出,則是中間節(jié)點輸出與其連接權(quán)系數(shù)之積的線性和。
圖1 大壩變形神經(jīng)網(wǎng)絡(luò)預(yù)測結(jié)構(gòu)圖
大壩變形混合預(yù)測模型是針對多元線性回歸分析和神經(jīng)網(wǎng)絡(luò)的缺點構(gòu)建起來的,即針對大壩實測數(shù)據(jù),先使用多元線性回歸方法,分析水壓等8個因素在大壩變形中的權(quán)重,保留影響因素較大的幾個參量,舍棄影響因素較小的參量,進而針對保留的參量,引入神經(jīng)網(wǎng)絡(luò)模型進行網(wǎng)絡(luò)訓(xùn)練,最終完成混合預(yù)測模型。這樣一方面克服了多元線性回歸分析在小容量樣本下精度不夠的問題,又解決了神經(jīng)網(wǎng)絡(luò)在變量較多時收斂速度慢,預(yù)測解振蕩大的特點。
針對貴州某大壩2002年5月—2007年7月實測數(shù)據(jù),共44 組,見表1 以前34 組數(shù)據(jù)作為大壩變形預(yù)測模型的訓(xùn)練樣本,以后10 組數(shù)據(jù)作為預(yù)測目標,分別采用多元線性回歸分析、神經(jīng)網(wǎng)絡(luò)和本文所建立的混合模型進行預(yù)測,并比較其精確度。
表1 貴州某大壩實測水平位移及相關(guān)數(shù)據(jù)表
圖2 多種預(yù)測分析結(jié)果比較圖
比較圖2 可知,多元線性回歸分析模型不論是樣本訓(xùn)練精度還是預(yù)測精度都比較差,只能粗略分析各影響因素之間的權(quán)重; 神經(jīng)網(wǎng)絡(luò)的樣本訓(xùn)練較好,但是預(yù)測精度稍差,這與輸入影響因素較多,各參量權(quán)重不明確有關(guān),而混合模型預(yù)測效果最好。
表2 多種預(yù)測分析結(jié)果比較 mm
表2 表明,采用多元線性回歸分析模型的預(yù)測結(jié)果最大誤差為1.5 mm,采用神經(jīng)網(wǎng)絡(luò)模型進行預(yù)測最大誤差為0.71 mm,采用本文提出的混合模型最大誤差為0.18 mm,誤差最小,預(yù)測精度最高,達到了預(yù)測目的。
另外,多元線性回歸分析結(jié)果表明,貴州某大壩變形的關(guān)鍵影響因素為: 水位、水位2 次方、水位3次方、氣溫和時效性,其中水位和氣溫的權(quán)重最大。
1) 采用多元線性回歸分析方法對貴州某大壩進行變形分析,結(jié)果表明:水位、水位2 次方、水位3 次方、氣溫和時效性,其中水位和氣溫的權(quán)重最大。這不僅為大壩變形研究提供了理論依據(jù),也為下一步展開神經(jīng)網(wǎng)絡(luò)預(yù)測分析奠定了基礎(chǔ)。
2) 采用多元線性回歸分析和神經(jīng)網(wǎng)絡(luò)構(gòu)建大壩變形預(yù)測模型是可行的,實例對比分析結(jié)果表明,采用該方法所獲得變形預(yù)測結(jié)果精度最高。
[1]吳中如. 水工建筑物安全監(jiān)控理論及其應(yīng)用[M]. 北京:高等教育出版社,2003.
[2]顧沖時,吳中如. 大壩與壩基安全監(jiān)控理論和方法及其應(yīng)用[M]. 南京:河海大學(xué)出版社,2006.
[3]黃世秀,洪天求,高飛. 基于小波消噪及BP 神經(jīng)網(wǎng)絡(luò)的大壩變形分析[J]. 人民長江,2011,42(09) :90 -93.
[4]朱金海,叢枝鮮,李秀海. 大壩變形的動態(tài)預(yù)報模型研究[J]. 煤炭技術(shù),2009,28(03) ,124 -126.
[5]李明然,田林亞,洪毅. 改進的灰色馬爾科夫模型在大壩變形預(yù)測中的應(yīng)用[J]. 水利與建筑工程學(xué)報,2012,10(03) :65 -67.