基于流速冪級數(shù)解的明渠縱向離散系數(shù)計算方法

張文俊,胡 煜,任華堂,夏建新

(中央民族大學(xué)生命與環(huán)境科學(xué)學(xué)院,北京 100081)

基于流速冪級數(shù)解的明渠縱向離散系數(shù)計算方法

張文俊,胡 煜,任華堂,夏建新

(中央民族大學(xué)生命與環(huán)境科學(xué)學(xué)院,北京 100081)

針對縱向離散系數(shù)計算多直接采用經(jīng)驗公式或在流速分布的求解中引入經(jīng)驗公式,不能完全反映離散物理機制的問題,基于矩形斷面一維水流控制方程,利用冪級數(shù)求解流速分布的理論解,提出了矩形斷面明渠縱向離散系數(shù)的計算方法,并進一步將其推廣至天然河流。利用該方法計算了104條天然河流的縱向離散系數(shù),結(jié)果表明計算值與實測值處于同一量級,計算結(jié)果在可接受范圍內(nèi)。

縱向離散系數(shù);明渠;冪級數(shù);流速分布

近10年來,河流突發(fā)水污染惡性事件時有發(fā)生,不僅對河流水生生態(tài)系統(tǒng)造成嚴(yán)重破壞,更為嚴(yán)重的是對下游大中城市供水安全造成威脅,如2005年松花江特大突發(fā)苯污染事故、2010年福建紫金礦業(yè)銅酸水滲漏事故、2012年廣西龍江鎘污染事故等。如何模擬污染物云團在河流中的遷移過程,準(zhǔn)確預(yù)測污染物到達(dá)各取水口時間以及恢復(fù)供水時間,是制定污染事故應(yīng)急預(yù)案的基礎(chǔ)。由于斷面流速分布的不均勻性,排放到河流中污染物形成的云團形狀將發(fā)生變化,即為河流的縱向離散作用。離散作用的強弱直接決定了污染云團遷移過程,在天然河流污染物預(yù)測中,由于計算量小、所需資料少等優(yōu)點,多采用一維縱向遷移離散方程求解,將流速分布不均勻?qū)е碌碾x散作用概化為離散項,而離散項準(zhǔn)確求解主要依賴于縱向離散系數(shù),因此,準(zhǔn)確確定縱向離散系數(shù)對于精準(zhǔn)預(yù)測河流突發(fā)污染的影響范圍和時間具有重要意義。

天然河流縱向離散系數(shù)為河流水環(huán)境模型的關(guān)鍵參數(shù),一直是環(huán)境水力學(xué)研究中的熱點問題。1954年Taylor[1]根據(jù)圓管試驗基于徑向擴散和縱向離散平衡首先提出了縱向離散系數(shù)的概念,Elder[2]將縱向離散系數(shù)的計算從圓管拓展到不考慮橫向差異的明渠中,基于垂向擴散和縱向離散的平衡得到離散系數(shù)計算公式。對于天然河流,Fischer等[3]指出明渠縱向離散主要是由于縱向流速的橫向分布不均勻引起的,而流速垂向差異的影響可忽略不計,并基于橫向擴散和縱向離散的平衡推導(dǎo)出三重積分形式的縱向離散系數(shù)計算公式:

(1)

(2)

在天然河流中,由于流體運動的非線性特性,無法得到流速的分布資料。在離散系數(shù)的確定上,除了示蹤試驗[4]方法之外,仍然以經(jīng)驗公式為主,通過量綱分析等方法將離散系數(shù)表示為河流特征參數(shù)平均流速、底部摩阻流速、水深、河寬等變量的函數(shù)[5],通過對實測數(shù)據(jù)的擬合分析得到了一系列公式,如Seo等[6]、Kashefipour等[7]、Sahay等[8]的公式,但這些公式物理意義不明確,且主要是用于特定類型的河流,不易推廣。

隨著對河流水動力特性認(rèn)識的不斷深入,基于合理假定得到流速分布特性,使得利用式(1)計算河流縱向離散系數(shù)具備可行性。鄧志強等[9]以拋物型斷面河道為對象,基于謝才公式得到橫斷面流速分布,繼而導(dǎo)出縱向離散系數(shù)理論公式;陳永燦等[10]基于最大熵值原理提出確定梯形斷面縱向流速分布的方法,得到梯形斷面明渠縱向離散系數(shù)計算公式; Wang 等[11]則通過矩形斷面的水流控制方程,采用傅里葉級數(shù)形式表達(dá)的流速分布,求解得到縱向離散系數(shù)計算公式。

雖然相關(guān)方法在計算離散系數(shù)中都取得了一定的效果,但是在流速分布求解中均不同程度地引入了經(jīng)驗公式或數(shù)理原理等,而采用流體動力學(xué)方程進行流速求解再得到離散系數(shù)的公式尚不多見。本文以常見的矩形斷面形態(tài)河流為對象,基于恒定均勻流條件下河流控制方程,采用冪級數(shù)形式求解得到河流流速分布,將流速分布代入式(1)求解可得矩形斷面河流縱向離散系數(shù)值,并將其推廣到天然河流縱向離散系數(shù)的計算。

1 縱向離散系數(shù)的冪級數(shù)求解方法

1.1 矩形斷面流速分布的求解

對于圖1所示的均勻恒定明渠流動,流速分布滿足以下動力學(xué)方程:

(3)

(4)

圖1 矩形河流斷面示意圖

令y′=2y/B,對河寬進行歸一化,式(4)化為

(5)

(6)

(7)

當(dāng)斷面最大流速umax已知時,則系數(shù)a0=umax,再代入式(6)和式(7),即可得到不同冪次的系數(shù);當(dāng)umax未知時,利用岸邊界位置流速為零的條件,得到

(8)

聯(lián)立式(6)(7)(8),可以得到不同冪次項的系數(shù)。在實際應(yīng)用求解中,冪級數(shù)項數(shù)一般可取6～10項,而對于寬深比較大的河流,所取項數(shù)則相應(yīng)增加。

將計算工況設(shè)置為H=1.0 m,改變寬深比B/H=10、50、100、150、200,Am=1.0 m2/s,cd=0.001,J=0.000 75,計算結(jié)果如圖2所示。由圖可2見,隨寬深比增加,斷面流速更加均勻:當(dāng)深度一定時,隨著河寬的增加,水流受岸邊摩擦力的影響減弱,流速分布更加趨于均勻,符合水流的運動規(guī)律,說明所得到的冪級數(shù)解是合理的。

圖2 矩形斷面縱向流速分布

1.2 矩形斷面縱向離散系數(shù)計算

(9)

將其代入式(1),可得:

(10)

(11)

其中,三重積分

(12)

2 天然河流縱向離散系數(shù)計算

2.1 計算條件

2.2 結(jié)果與討論

為分析本文計算方法計算結(jié)果的精度,采用對數(shù)比偏差KDR、對數(shù)平均絕對誤差Me和對數(shù)均方根Rms這3個統(tǒng)計指標(biāo)進行評價,定義如下:

(13)

(14)

(15)

式中:Dlp為本文公式計算值;Dlm為實測值。

對數(shù)比偏差反映計算值與實測值的接近程度,當(dāng)計算值與實測值之比越接近1時,則KDR越接近0,而KDR∈(-0.3,0.3)時,則計算值與實測值比值為0.5～2。對數(shù)平均絕對誤差反映出計算值與實測值的整體偏離情況，對數(shù)均方根反映計算值與實測值的離散情況，對數(shù)平均絕對誤差和對數(shù)均方根值越接近于0,說明計算值與實測值越接近,方程預(yù)測效果越好。

將本文公式計算值與實測值比較(圖3),可見本文公式計算值較為均勻地分布在實測值兩側(cè),但也存在一定的偏離,這主要是由于本文流速的冪級數(shù)解是基于矩形斷面均勻流假定得到的,而天然河流的斷面形式多種多樣,將其統(tǒng)一概化為水面寬度與平均水深的矩形斷面,會出現(xiàn)一定程度的失真,縱向離散系數(shù)的計算值和實測值之間存在一定的誤差。

圖3 縱向離散系數(shù)計算值與實測值的比較

將本文公式計算值與理論性較強的Wang等[11]、Deng等[13]的公式計算值進行了比較,結(jié)果見圖3。由圖3可看出,本文公式計算精度優(yōu)于Wang等[11]的公式,但和Deng等[13]的公式相比則偏低。這是由于Deng等[13]的公式在離散系數(shù)的推導(dǎo)中流速分布由謝才公式計算,而謝才公式基于天然河流大量數(shù)據(jù)得到,對于自然河流具有更好的適用性。本文的流速冪級數(shù)解是基于矩形斷面均勻流條件下的動力學(xué)方程得到的,故在將自然河流概化為矩形斷面過程中出現(xiàn)了一定程度的失真。盡管如此,由于流速分布公式的物理意義相較謝才公式更為明確且計算精度總體不低于其他公式,本文計算方法仍然具有一定的意義。

表1 不同縱向離散系數(shù)經(jīng)驗公式對比分析

為進一步分析KDR的誤差與寬深比的關(guān)系,圖4給出了不同寬深比下的誤差分布情況。當(dāng)寬深比較小時,本文公式計算值和Deng等[13]的公式計算值與實測值偏差處在相對合理的區(qū)間,而當(dāng)寬深比大于100時,則所有理論公式計算結(jié)果的誤差都會增大。

圖4 對數(shù)比偏差與寬深比的關(guān)系

本文公式計算值與前人公式的計算值對比分析見表1,本文公式對于104條天然河流計算值的KDR分布在(-1,1)區(qū)間的比例為90%,即計算值與實測值基本在一個量級,與其他傳統(tǒng)公式在同一水平;對數(shù)平均絕對誤差Me為0.46,即本文公式整體計算值在實測值的0.34～2.8倍之間;對數(shù)均方根Rms為0.054,計算結(jié)果波動性不大,與各傳統(tǒng)公式計算值誤差相比相差不大,處于居中水平。

3 結(jié) 語

基于一維均勻恒定流控制方程,利用冪級數(shù)解法得到斷面流速分布表達(dá)式,該解法能夠較好地反映流體運動的力學(xué)機制,且形式簡單,易于計算,與天然河流流動規(guī)律相吻合。通過將流速分布的冪級數(shù)解代入Fischer理論公式推導(dǎo)了矩形明渠縱向離散系數(shù)計算方法,該方法物理機制明確,可根據(jù)精度要求確定求解項數(shù)。將該縱向離散系數(shù)計算方法推廣到天然河流中,并利用104條天然河流的實測資料進行了驗證,結(jié)果顯示計算值和實測值量級一致,能夠反映天然河流的離散特性。

[1] TAYLOR G. The dispersion of matter in turbulent flow through a pipe[J]. Proceedings of the Royal Society of London, 1954, 223(1155):446-468.

[2] ELDER J W.The dispersion of marked fluid in turbulent shear flow[J].Journal of Fluid Mechanics,1959,5(4):544-560.

[3] FISCHER H B,LIST J,KOHR C Y,et al.Mixing in inland and coastal water [M].New York: Academic Press,1979:81-147.

[4] 李駿旻,黃平,李毅能,等.近似擬合法確定離散系數(shù)的適用性和選點策略[J].水利水電科技進展,2008,28(5):10-12.(LI Junmin,HUANG Ping,LI Yineng,et al.Applicability of approximate fitting method in estimation of longitudinal diffusivity coefficient[J].Advances in Science and Technology of Water Resources,2008,28(5):10-12.(in Chinese))

[5] 顧莉,華祖林.天然河流縱向離散系數(shù)確定方法的研究進展[J].水利水電科技進展,2007,27(2):85-89.(GU Li,HUA Zulin.Advances in determination of longitudinal dispersion coefficient of natural rivers[J].Advances in Science and Technology of Water Resources,2007,27(2):85-89.(in Chinese))

[6] SEO I W,CHEONG T S.Predicting longitudinal dispersion coefficient in natural streams[J].Journal of Hydraulic Engineering,1998,124(1):25-32.

[7] KASHEFIPOUR S M,FALCONER R A.Longitudinal dispersion coefficients in natural channels[J].Water Research, 2002,36(6):1596-608.

[8] SAHAY R R,DUTTA S.Prediction of longitudinal dispersion coefficients in natural rivers using genetic algorithm[J]. Hydrology Research,2009,40(6):544-552.

[9] 鄧志強,褚君達(dá).河流縱向分散系數(shù)研究[J].水科學(xué)進展,2001,12(2):137-142.(DENG Zhiqiang,CHU Junda.Longitudinal dispersion coefficient in natural rivers [J].Advances in Water Science,2001,12(2):137-142.(in Chinese))

[10] 陳永燦,朱德軍.梯形斷面明渠中縱向離散系數(shù)研究[J].水科學(xué)進展,2005,16(4):511-517.(CHEN Yongcan,ZHU Dejun.Study on longitudinal dispersion coefficient in open channel trapezoidal cross-section [J].Advances in Water Science,2005,16(4):511-517.(in Chinese))

[11] WANG Y, HUAI W. Estimating the longitudinal dispersion coefficient in straight natural rivers[J]. Journal of Hydraulic Engineering, 2016, 142(11):04016048.

[12] ZENG Yuhong,HUAI Wenxing.Estimation of longitudinal dispersion coefficient in rivers[J].Journal of Hydro-Environment Research,2014,8(1):2-8.

[13] DENG Zhiqiang,SINGGH V P,BENGTSSON L.Longitudinal dispersion coefficient in straight rivers[J].Journal of Hydraulic Engineering,2001,127(11):919-927.

[14] LIU H.Predicting dispersion coefficient of streams[J].Journal of the Environmental Engineering Division,1977,103(1): 59-69.

Computationoflongitudinaldispersioncoefficientinopenchannelsbasedonpowerseriessolutiontovelocitydistribution//

ZHANG Wenjun,HU Yu,REN Huatang,XIA Jianxin

(CollegeofLifeandEnvironmentalSciences,MinzuUniversityofChina,Beijing100081,China)

The longitudinal dispersion coefficient is usually computed directlyby empirical formulas or indirectly through introducing empirical formulas during the solution tovelocity distributions,which cannot fully reflect the physical mechanisms of dispersion.Based on the one-dimensional governing equation of river dynamics, a method to calculate longitudinal dispersion coefficient involving a power series solution to velocity distributionwas proposed.This formulation has been extended to the applications in 104 natural rivers with simplified rectanglecross-sections.The computed values are in the same order of magnitudes with the measured ones, indicating acceptable results.

longitudinal dispersion coefficient; open channel; power series; velocity distributions

國家自然科學(xué)基金(51479218)

張文俊(1992—),男,碩士研究生,主要從事環(huán)境水力學(xué)研究。E-mail:s151093@muc.edu.cn

任華堂(1976—),男,副教授,博士,主要從事環(huán)境水力學(xué)研究。E-mail:renhuatang@muc.edu.cn

10.3880/j.issn.1006-7647.2017.06.008

TV133

A

1006-7647(2017)06-0044-04

2016-12-11 編輯：熊水斌)