河道一維非恒定流水溫數(shù)學模型研究

1、河道一維非恒定流水溫數(shù)學模型研究     摘要：本文介紹了影響河道水溫的各種因素，并分析各因素對河道水溫的影響程度。根據(jù)引黃入津河道實測資料，確定了縱向離散系數(shù)的經(jīng)驗公式，建立了適用于河道的一維非恒定流水溫數(shù)學模型。通過對實際河道水溫的數(shù)值模擬，結(jié)果表明，該模型具備模擬河道縱向水溫變化的功能。此模型可以用于河道的水溫模擬和預測，為河道實施生態(tài)調(diào)水方案和水溫對環(huán)境的影響提供科學依據(jù)。 關(guān)鍵詞：非恒定流；水溫；數(shù)值模擬；河道 Abstract: this paper introduces the influence of various factors of

2、 river water temperature, and analyzes the factors on the influence degree of the river water temperature. According to the measured data, the Yellow River diversion into the tianjin river, determined the longitudinal dispersion coefficient of experience formula, is established through the channel o

3、f 1-d unsteady flow water temperature mathematical model. Through to the real river water temperature of numerical simulation, the results show that the model has simulation river water temperature changes of longitudinal function. This model can be used to simulate the river water temperature and t

4、he forecast, for the implementation of ecological river water temperature scheme and the influence of the environment to provide the scientific basis. Key words: the unsteady flow; Water temperature; Numerical simulation; river 中圖分類號：TV697 文獻標識碼：A 1河道水溫影響因素 引起河道水溫變化的因素主要有四個方面，即水面熱交換、與河床的能量交換、內(nèi)部產(chǎn)生的熱和

5、人為的加熱和減熱。 1.1 水面熱交換 水面熱交換包括太陽短波輻射、凈蒸發(fā)、熱傳導和凈長波輻射四個方面。 (1) 式中， ， 為太陽短波輻射強度(W/m2)； 為凈蒸發(fā)(W/m2)； 為熱傳導(W/m2)； 為凈長波輻射(W/m2)， ， 為大氣長波輻射， 為水面對長波的反射， 為水體的長波返回輻射。 1.2 與河床的能量交換 水體和河床的熱交換 是通過固體的熱傳導確定的， (2) 為河床溫度， 為河床熱傳導系數(shù)(cal/(ms)。水體與河床的能量交換，對于一般河流而言，因為河床溫度 的梯度很小，即 ，熱量交換很小，故可忽略不記，即 。但是對于封凍河道和潛水河道，特別是水溫晝夜溫差變化較大的河

6、道，水體與河床之間的熱交換則不可忽略。 1.3 內(nèi)部產(chǎn)生的熱 水流內(nèi)部產(chǎn)生的熱產(chǎn)生有兩種因素，一種為水的機械能通過摩擦轉(zhuǎn)化為熱；另一種為水中物質(zhì)的化學能經(jīng)過生物化學作用轉(zhuǎn)換成熱。無論是通過摩擦轉(zhuǎn)化成熱還是通過化學作用轉(zhuǎn)化成熱，熱量產(chǎn)生很小，此項可忽略不計。 1.4 人為的加熱和減熱 若河流沿程無人為的排熱或減熱，此項可忽略不計。 上游來水溫度和氣象條件是對下游河道水溫影響最主要的因素。氣象因素主要包括氣溫、太陽輻射強度、總云量、相對濕度、風速、露點溫度等。根據(jù)靈敏度分析，各影響因素由大到小的排列順序依次為上游來水水溫、空氣溫度、太陽短波輻射。因此運用數(shù)學模型預測壩址下游河道水溫的沿程變化規(guī)律時

7、，應保證這三個參數(shù)有足夠的可靠性和準確性。 2河道水溫數(shù)學模型 2.1 圣維南方程組 圣維南方程組由水流連續(xù)方程和動量方程組成，以Z、Q為變量的圣維南方程組形式如下： 水流連續(xù)方程： (3) 水流動量方程： (4) 式中，Z為水位(m)；Q為流量(m3/s)；A為過水斷面面積(m2)；n為糙率；R為水力半徑(m)；C為謝才系數(shù)。 2.2 溫度控制方程 一維河道溫度控制方程為對流-擴散方程，形式如下： (5) 式中，Ex為縱向彌散系數(shù)(m2/s)；為水的密度(kg/m3)；Cp為水的比熱(cal/(kg)；B為河面寬度(m)；S為單位表面積凈熱交換通量，表示水流與外界(太陽、空氣、河道邊界)之間

8、的熱交換量(W/m2)。 2.3 方程求解 圣維南方程采用Preissmann四點隱式加權(quán)格式，由此可得到系數(shù)為三對角矩陣的非線性方程組，用TDMA(TriDiagonal-Matrix Algorithm)法進行求解。 溫度控制方程，對時間項采用向前差分，空間采用中心差分，可得： 整理可得： (6) 其中 采用擴散格式求解方程(6)。 3參數(shù)率定 縱向彌散系數(shù) 是水溫模型中的一個重要參數(shù)，其取值是否合理直接關(guān)系到水溫模擬的精度。 是隨空間變化的，確定其分布十分困難。一般來說， 應分為漸進值、恒定值和瞬時值。為方便起見，天然河流中 一般都取不隨時間變化的漸進值來代替隨時間變化的實際值。 目前，

9、天然河流的彌散系數(shù) 都是由經(jīng)驗公式確定的，主要有： (1) 式中， 為斷面平均流速(m/s)； 為河寬(m)； 為平均水深(m)； 為摩阻流速(m/s)， ，g為重力加速度(m/s2)， 為水力半徑(m)， 為水力坡降。 (2) MCQuivey and Keefer經(jīng)驗公式： 式中， 為水面坡降， 為水面寬度(m)， 為流量(m3/s)。 (3)Fischer(1973)和Jain(1976)經(jīng)驗公式： 式中，0.001< <0.016， 為平均流速(m/s)； 為水面寬度(m)； 為橫向擴散系數(shù)。 (4)Liu(1977)經(jīng)驗公式： 式中， 為水力半徑(m)； ， 為摩阻流速(

10、m/s)， 為平均流速(m/s)。 經(jīng)過比較以上四個經(jīng)驗公式后可知，根據(jù)不同研究水體的得到的經(jīng)驗公式的差別是比較大的。針對引黃濟津河道的實測資料對上述公式進行校驗，認為經(jīng)驗公式(1)求得的各個斷面的 比較穩(wěn)定，故采用該計算式計算縱向彌散系數(shù)。 4模型驗證 采用某河流的沿程水溫對數(shù)學模型進行驗證。河道長20km，計算斷面間距為1km，共20個計算斷面，計算時間30天。 計算值和實測值如圖1所示。計算值和實測值吻合較好，表明該模型能準確的模擬河道的沿程水溫，對水溫進行預測分析。具有較高的實用性和準確性。 5參考文獻 1 莊春義. 河道一維非恒定流水溫預測模型研究D. 四川: 四川大學. 2005. 2 郝紅升, 李克峰, 莊春義. 關(guān)于河道一維非恒定流