唐騰飛，黃耀英，周宜紅，郝 鵬，王 飛

(1.中國(guó)電建集團(tuán)貴陽(yáng)勘測(cè)設(shè)計(jì)研究院有限公司，貴州貴陽(yáng)550081；2.三峽大學(xué)水利與環(huán)境學(xué)院，湖北宜昌443002)

高拱壩施工期垂直向溫度分布智能動(dòng)態(tài)預(yù)測(cè)

唐騰飛1，黃耀英2，周宜紅2，郝 鵬1，王 飛1

(1.中國(guó)電建集團(tuán)貴陽(yáng)勘測(cè)設(shè)計(jì)研究院有限公司，貴州貴陽(yáng)550081；2.三峽大學(xué)水利與環(huán)境學(xué)院，湖北宜昌443002)

針對(duì)當(dāng)前特高拱壩施工期混凝土溫控采取分區(qū)分時(shí)段通水冷卻難以實(shí)時(shí)掌握并準(zhǔn)確預(yù)測(cè)垂直向溫度分布的問(wèn)題，利用分布式光纖實(shí)時(shí)在線監(jiān)測(cè)采集數(shù)據(jù)量大的優(yōu)勢(shì)，發(fā)掘其線空間內(nèi)的數(shù)據(jù)關(guān)系，基于廣義回歸神經(jīng)網(wǎng)絡(luò)建立隱式垂直向溫度分布模型，并對(duì)某處于建設(shè)期的高拱壩垂直向溫度分布進(jìn)行動(dòng)態(tài)預(yù)測(cè)分析。結(jié)果表明，建立的垂直向溫度預(yù)測(cè)模型的預(yù)測(cè)溫度與實(shí)測(cè)溫度之間存在良好的相關(guān)性，能很好地反映各溫控因素與監(jiān)測(cè)數(shù)據(jù)的非線性關(guān)系，能夠準(zhǔn)確、合理地預(yù)測(cè)垂直向溫度。

高拱壩；施工期；廣義回歸神經(jīng)網(wǎng)絡(luò)；垂直向溫度分布；智能動(dòng)態(tài)預(yù)測(cè)

0 引 言

目前特高拱壩的典型溫降過(guò)程和分區(qū)冷卻方案的最大特點(diǎn)就是，在溫度控制上按“小溫差，早冷卻，緩慢冷卻”[1]指導(dǎo)思想將溫度峰值出現(xiàn)后的降溫分為若干段降溫和控溫過(guò)程，在高度方向分別設(shè)擬灌區(qū)、同冷區(qū)、過(guò)渡區(qū)和蓋重區(qū)，自上而下形成動(dòng)態(tài)溫度梯度。計(jì)算結(jié)果表明，二期水冷時(shí)的溫差大于14 ℃時(shí)，會(huì)出現(xiàn)超標(biāo)拉應(yīng)力。因此，二期水冷時(shí)的降溫幅度要嚴(yán)格控制，上下層溫差應(yīng)控制在14 ℃以下[2]，上下層溫差問(wèn)題引起的裂縫則可通過(guò)合理控制溫降過(guò)程和設(shè)置上下層溫度梯度來(lái)解決[3]。在高拱壩施工期安全監(jiān)測(cè)中，垂直向溫度分布是大壩整體施工進(jìn)度、仿真分析、溫控措施調(diào)整優(yōu)化的基本點(diǎn)。因此，對(duì)大壩施工期的垂直向溫度過(guò)程進(jìn)行動(dòng)態(tài)監(jiān)測(cè)十分必要。

目前國(guó)內(nèi)外大壩安全監(jiān)測(cè)中現(xiàn)有監(jiān)測(cè)模型一是基于統(tǒng)計(jì)方法建立的回歸分析模型，二是基于模糊數(shù)學(xué)或灰色系統(tǒng)理論建立的監(jiān)控模型，三是基于優(yōu)化算法和人工智能技術(shù)建立的智能監(jiān)控模型[4]。由于大壩施工期工作條件復(fù)雜，溫度影響因素眾多，無(wú)法用確切的定量關(guān)系描述，因此很難建立考慮各因素的非線性數(shù)學(xué)監(jiān)測(cè)模型。模糊數(shù)學(xué)模型需要對(duì)各項(xiàng)指標(biāo)賦予不同的權(quán)值，而權(quán)值的選取不可避免地帶有一些主觀性和隨意性，容易導(dǎo)致結(jié)果失真。王利杰[5]等建立了施工期混凝土壩垂直向溫度分布模型，可以采用逐步回歸分析法回歸確定式中的各系數(shù)，該模型本質(zhì)上來(lái)說(shuō)是統(tǒng)計(jì)模型，雖然簡(jiǎn)單易行，但存在著非線性映射不強(qiáng)等問(wèn)題，不能充分反映各影響因素序列中的非線性特性，有時(shí)甚至不能保證預(yù)測(cè)精度；此外統(tǒng)計(jì)模型需要做一定的假設(shè)和事先已知確定因子。而神經(jīng)網(wǎng)絡(luò)模型具有處理高度非線性問(wèn)題的能力，且采用隱性表達(dá)式，不必像統(tǒng)計(jì)模型需要假設(shè)和事先已知確定因子，在理論上可以實(shí)現(xiàn)任意函數(shù)的逼近[6- 7]。本文采用人工神經(jīng)網(wǎng)絡(luò)建立各輸入量與垂直向溫度分布之間的非線性關(guān)系，動(dòng)態(tài)預(yù)測(cè)未來(lái)N天的溫度信息。

1 基于廣義回歸神經(jīng)網(wǎng)絡(luò)垂直向溫度分布模型

1.1 廣義回歸網(wǎng)絡(luò)理論[8]

廣義回歸神經(jīng)網(wǎng)絡(luò)由輸入層、模式層、求和層、輸出層4層組成。

(1) 輸入層神經(jīng)元的數(shù)目等于學(xué)習(xí)樣本中輸入向量的維數(shù)k，各神經(jīng)元直接將輸入變量傳遞給模式層。

(2) 模式層神經(jīng)元數(shù)目等于學(xué)習(xí)樣本的數(shù)目，各神經(jīng)元對(duì)應(yīng)不同的樣本，模式層神經(jīng)元的傳遞函數(shù)為

(1)

式中，X為網(wǎng)絡(luò)輸入變量；Xi為第i個(gè)神經(jīng)元對(duì)應(yīng)的學(xué)習(xí)樣本；σ為光滑因子。

1.2 垂直向溫度分布的隱式模型

高拱壩垂直向溫度分布的預(yù)測(cè)建模問(wèn)題可歸結(jié)為建立垂直向溫度與其影響因素之間非線性輸入輸出關(guān)系，所以，高拱壩垂直向溫度分布預(yù)測(cè)問(wèn)題可概括為建立多輸入、多層的網(wǎng)絡(luò)結(jié)構(gòu)?；炷翝仓}(cāng)垂直向溫度分布是一個(gè)復(fù)雜多因素問(wèn)題，涉及到混凝土澆筑倉(cāng)的空間布置、澆筑時(shí)機(jī)、混凝土水泥水化反應(yīng)、通水冷卻、表面保溫等方面。從場(chǎng)論出發(fā)，應(yīng)用冪級(jí)數(shù)的處理方法，并結(jié)合大體積混凝土的熱傳導(dǎo)理論，在建立神經(jīng)網(wǎng)絡(luò)模型時(shí)，采用4個(gè)周期項(xiàng)sin(2πt/365)、cos(2πt/365)、sin(4πt/365)、cos(4πt/365)作為環(huán)境氣溫的影響因子，兩個(gè)指數(shù)函數(shù)1-e-0.318t，1-e-0.295t作為水化熱溫升和通水冷卻的影響因子，垂直向高程Z構(gòu)成的處理樣本數(shù)為網(wǎng)絡(luò)的輸入，其相應(yīng)的溫度T作為單一的輸出建立網(wǎng)絡(luò)模型，構(gòu)建混凝土澆筑倉(cāng)垂直向溫度分布的隱式網(wǎng)絡(luò)模型如圖1所示。

圖1 混凝土垂直向溫度分布的隱式模型

1.3 混凝土澆筑倉(cāng)垂直向溫度分布的智能動(dòng)態(tài)預(yù)測(cè)模型

根據(jù)當(dāng)前垂直向溫度分布信息，通過(guò)隱式模型訓(xùn)練當(dāng)前溫控信息數(shù)據(jù)并預(yù)測(cè)未來(lái)N天的垂直向溫度分布，通過(guò)分析預(yù)測(cè)的未來(lái)N天的垂直向溫度數(shù)據(jù)，評(píng)價(jià)當(dāng)前溫控措施下未來(lái)時(shí)刻的垂直向溫度分布是否滿足設(shè)計(jì)條件(即上下層容許溫差，冷卻水與混凝土容許溫差等)，如滿足，可繼續(xù)實(shí)施當(dāng)前溫控措施；如不滿足，可根據(jù)分析結(jié)果，調(diào)整當(dāng)前溫控措施，使其未來(lái)時(shí)間段內(nèi)的垂直向溫度分布向設(shè)計(jì)值靠攏直至滿足要求。待監(jiān)測(cè)到未來(lái)N天垂直向溫度分布信息實(shí)測(cè)值后，更新預(yù)測(cè)模型輸入輸出數(shù)據(jù)，重復(fù)上述流程，即可達(dá)到動(dòng)態(tài)更新歷史數(shù)據(jù)、動(dòng)態(tài)訓(xùn)練隱式網(wǎng)絡(luò)模型、動(dòng)態(tài)預(yù)測(cè)未來(lái)N天垂直向溫度分布、動(dòng)態(tài)評(píng)價(jià)未來(lái)N天垂直向溫度分布、動(dòng)態(tài)調(diào)控當(dāng)前溫控措施的效果?；炷翝仓}(cāng)垂直向溫度分布的動(dòng)態(tài)預(yù)測(cè)模型見(jiàn)圖2。

圖2 混凝土澆筑倉(cāng)垂直向溫度分布的動(dòng)態(tài)預(yù)測(cè)模型

2 工程實(shí)例

2.1 工程背景

以西南某建設(shè)中的混凝土雙曲拱壩為例，在5、15、16、23號(hào)4個(gè)典型壩段埋設(shè)分布式光纖監(jiān)測(cè)混凝土澆筑塊溫度場(chǎng)，實(shí)時(shí)在線采集了大量的壩體混凝土內(nèi)部溫度數(shù)據(jù)。本文基于這些實(shí)測(cè)溫度數(shù)據(jù)分析壩體垂直向溫度分布的相關(guān)規(guī)律。

表1 隱式網(wǎng)絡(luò)模型部分訓(xùn)練樣本

sin2πt365coσ2πt365sin4πt365cos4πt3651-e-0.318t1-e-0.295t高程/m溫度/℃0.03440.99940.06880.99760.47060.4457441.312.00000.03440.99940.06880.99760.47060.4457444.311.91600.03440.99940.06880.99760.47060.4457447.311.54100.03440.99940.06880.99760.47060.4457450.311.95200.03440.99940.06880.99760.47060.4457453.311.72300.03440.99940.06880.99760.47060.4457489.319.64200.03440.99940.06880.99760.47060.4457492.319.88500.03440.99940.06880.99760.47060.4457495.320.59900.03440.99940.06880.99760.47060.4457498.321.8770

2.2 光纖布置和樣本數(shù)據(jù)選擇

該拱壩典型壩段(15號(hào)壩段)在基礎(chǔ)約束區(qū)的澆筑倉(cāng)厚1.5 m，脫離基礎(chǔ)約束區(qū)后，澆筑倉(cāng)厚為3.0 m。為監(jiān)測(cè)該壩段溫度狀態(tài)，從該壩段澆筑第一倉(cāng)混凝土開(kāi)始，就在每個(gè)澆筑倉(cāng)中間層均埋設(shè)分布式光纖實(shí)時(shí)在線測(cè)溫。采用分布式光纖測(cè)溫主機(jī)進(jìn)行監(jiān)測(cè)時(shí)，將測(cè)溫系統(tǒng)的空間分辨率設(shè)置為1.0 m，采樣頻率為2 h/次，采樣空間間隔設(shè)置為1.0 m，這樣可保證沿分布式光纖傳感網(wǎng)絡(luò)每米即獲得一個(gè)溫度測(cè)值。將光纖測(cè)溫?cái)?shù)據(jù)上傳至大壩信息管理系統(tǒng)時(shí)，從外到內(nèi)依次選取10個(gè)光纖測(cè)點(diǎn)的溫度作為該混凝土澆筑倉(cāng)的溫度數(shù)據(jù)。多個(gè)光纖測(cè)點(diǎn)的平均溫度作為該澆筑倉(cāng)混凝土的當(dāng)前實(shí)測(cè)溫度。

選取15號(hào)壩段440～500 m高程范圍內(nèi)20個(gè)混凝土澆筑倉(cāng)2011- 11- 20～2012- 01- 13分布式光纖實(shí)測(cè)溫度作為垂直向溫度分布信息，在此時(shí)間段的初始時(shí)刻48～51號(hào)倉(cāng)處于二期控溫階段，52～56號(hào)倉(cāng)處于二期冷卻階段，57～62號(hào)倉(cāng)處于中期冷卻階段，63～67號(hào)倉(cāng)處于一期冷卻階段，此后這些混凝土澆筑倉(cāng)溫度由于溫控措施的影響發(fā)生動(dòng)態(tài)變化，隨之其溫控階段也不斷調(diào)整變化。選2011- 11- 20～2012- 12- 31實(shí)測(cè)數(shù)據(jù)作為隱式模型的訓(xùn)練值，2012- 1- 1～2012- 01- 14數(shù)據(jù)用來(lái)預(yù)測(cè)檢驗(yàn)。隱式網(wǎng)絡(luò)模型部分訓(xùn)練樣本見(jiàn)表1。

2.3 訓(xùn)練預(yù)測(cè)和結(jié)果分析

以2011- 11- 20～2012- 12- 31分布式光纖實(shí)測(cè)溫度信息作為訓(xùn)練樣本預(yù)測(cè)此后7d的溫度信息，即以前39個(gè)時(shí)刻的溫度數(shù)據(jù)作為樣本，預(yù)測(cè)接下來(lái)7d共140個(gè)數(shù)據(jù)點(diǎn)(2011- 12- 31～2012- 01- 06)的垂直向溫度；隱式網(wǎng)絡(luò)擬合效果見(jiàn)圖3a，預(yù)測(cè)效果見(jiàn)圖3b。以2011- 11- 20～2012- 01- 06分布式光纖實(shí)測(cè)溫度信息作為訓(xùn)練樣本預(yù)測(cè)此后7d的溫度信息，即以前46個(gè)時(shí)刻的溫度數(shù)據(jù)作為樣本，預(yù)測(cè)接下來(lái)7d共140個(gè)數(shù)據(jù)點(diǎn)(2012- 01- 07～2012- 01- 13)的垂直向溫度；隱式網(wǎng)絡(luò)擬合效果見(jiàn)圖4a，預(yù)測(cè)效果見(jiàn)圖4b。

圖3 混凝土澆筑倉(cāng)垂直向溫度分布隱式模型的擬合預(yù)測(cè)效果

圖4 混凝土澆筑倉(cāng)垂直向溫度分布隱式模型的動(dòng)態(tài)訓(xùn)練預(yù)測(cè)

(1)從圖3可以看出，垂直向溫度分布廣義回歸隱式網(wǎng)絡(luò)模型能夠很好地完成對(duì)未來(lái)溫度的預(yù)測(cè)，未來(lái)7d預(yù)測(cè)值最大誤差為0.19 ℃，模型擬合數(shù)據(jù)、預(yù)測(cè)數(shù)據(jù)與實(shí)測(cè)溫度數(shù)據(jù)的復(fù)相關(guān)系數(shù)均在0.99以上。計(jì)算結(jié)果表明建立的垂直向溫度分布隱式模型預(yù)測(cè)精度能夠滿足工程實(shí)際要求。

(2)對(duì)比圖3、4可知，將隱式模型預(yù)測(cè)時(shí)刻(2011- 12- 31～2012- 01- 06)的溫度實(shí)測(cè)值作為動(dòng)態(tài)預(yù)測(cè)模型的歷史數(shù)據(jù)，更新輸入樣本，構(gòu)建隱式動(dòng)態(tài)預(yù)測(cè)模型，進(jìn)預(yù)測(cè)未來(lái)7 d(2012- 01- 07～2012- 01- 13)垂直向溫度，預(yù)測(cè)溫度最大誤差為0.24 ℃，滿足預(yù)測(cè)精度的要求，且能為動(dòng)態(tài)分析評(píng)價(jià)線空間內(nèi)的溫度分布提供數(shù)據(jù)支撐。

3 結(jié) 論

(1)在建立神經(jīng)網(wǎng)絡(luò)模型時(shí)，從場(chǎng)論出發(fā)，應(yīng)用冪級(jí)數(shù)的處理方法，并結(jié)合大體積混凝土的熱傳導(dǎo)理論，采用4個(gè)正余弦周期項(xiàng)作為環(huán)境氣溫的影響因子，兩個(gè)指數(shù)函數(shù)作為水化熱溫升和通水冷卻的影響因子，垂直向高程Z坐標(biāo)構(gòu)成的數(shù)據(jù)樣本作為隱式網(wǎng)絡(luò)模型的輸入，實(shí)測(cè)溫度作為輸出。實(shí)踐表明，建立的隱式模型預(yù)測(cè)精度較高，能夠滿足工程實(shí)際需要。

(2)鑒于處于施工期的高拱壩各混凝土澆筑倉(cāng)在不同的溫控階段采取不同的溫控措施，導(dǎo)致垂直向分布時(shí)刻處于變化中，結(jié)合動(dòng)態(tài)控制的理念，在建立的隱式模型基礎(chǔ)上，構(gòu)建了施工期高拱壩垂直向溫度分布的智能動(dòng)態(tài)預(yù)測(cè)模型，能為動(dòng)態(tài)分析評(píng)價(jià)線空間內(nèi)的溫度分布提供數(shù)據(jù)支撐。

[1]朱伯芳. 小溫差早冷卻緩慢冷卻是混凝土壩水管冷卻的新方向[J]. 水利水電技術(shù). 2009， 40(1)： 44- 50．

[2]張國(guó)新， 趙文光. 特高拱壩溫度應(yīng)力仿真與溫度控制的幾個(gè)問(wèn)題探討[J]. 水利水電技術(shù). 2008， 39(10)： 36- 40．

[3]劉有志， 張國(guó)新， 劉毅. 特高拱壩施工期裂縫成因分析與溫控防裂方法新思路[C]∥北京力學(xué)會(huì)第17屆學(xué)術(shù)年會(huì)論文集， 2011．

[4]楊杰， 吳中如， 顧沖時(shí). 大壩變形監(jiān)測(cè)的BP網(wǎng)絡(luò)模型與預(yù)報(bào)研究[J]. 西安理工大學(xué)學(xué)報(bào)， 2001， 17(1)： 25- 29．

[5]WANG Lijie， HUANG Yaoying， ZHOU Yihong， et al. Primary vertical temperature distribution model of a concrete dam during construction based on distributed optical fiber measuring temperature[J]. WIT Transactions on the Built Environment， 2013(140)： 353- 359 ．

[6]韓力群. 人工神經(jīng)網(wǎng)絡(luò)教程[M]. 北京： 北京郵電大學(xué)出版社， 2006： 39- 41．

[7]戴奎. 神經(jīng)網(wǎng)絡(luò)實(shí)現(xiàn)技術(shù)[M]. 北京： 國(guó)防科技大學(xué)出版社， 1998： 303- 309．

[8]史峰， 王小川. 郁磊. matlab神經(jīng)網(wǎng)絡(luò)30個(gè)案例分析[M]. 北京： 北京航空航天大學(xué)出版社， 2010．

(責(zé)任編輯焦雪梅)

Intelligent Dynamic Prediction about the Vertical Temperature Distribution of High Arch Dam during Construction Period

TANG Tengfei1, HUANG Yaoying2, ZHOU Yihong2, HAO Peng1, WANG Fei1

(1. PowerChina Guiyang Engineering Corporation Limited, Guiyang 550081, Guizhou, China;2. College of Hydraulic & Environmental Engineering, Three Gorges University, Yichang 443002, Hubei, China)

The concrete temperature control in construction period of high arch dam always adopts water cooling in different time and space, but it is difficult to grasp and predict the vertical temperature distribution. By using abundant online monitoring data obtained from distributed optical fiber, the relationship of data in line space is explored based on generalized regression neural network to set up an implicit expression model about the vertical temperature distribution. The model is used to predict the vertical temperature distribution of a high arch dam in construction period. The results show that the predicted temperature of model has a good relation with actual online monitor temperature and it can reflect the non-linear relationship between temperature control factors and monitoring data. The model can accurately and reasonably predict vertical temperature distribution．

high arch dam; construction period; generalized regression neural network; vertical temperature distribution; intelligent dynamic prediction

2015- 01- 26

國(guó)家自然科學(xué)基金(51079079，51209124)

唐騰飛(1987—)，男，湖北荊州人，碩士，從事水利水電工程設(shè)計(jì)工作．

TV698.1;TV642.4

A

0559- 9342(2015)12- 0054- 04