王志云 滕 凱

（1.齊齊哈爾市水利工程質(zhì)量與安全監(jiān)督站，黑龍江齊齊哈爾 161006；2.齊齊哈爾市水務(wù)局，黑龍江齊齊哈爾 161006）

由于拋物線形斷面具有良好的水流條件及力學(xué)性能，因此在城市供排水及水利水電工程中，工程設(shè)計(jì)部門常根據(jù)輸水流量及地質(zhì)條件，將渠道斷面設(shè)計(jì)為半立方拋物線、二次拋物線或三次拋物線形式之一，因此相關(guān)拋物線形式的水力計(jì)算問(wèn)題也逐漸引起了有關(guān)學(xué)者的重視[1～5]，其中關(guān)于半立方拋物線及二次拋物線形渠道斷面收縮水深的計(jì)算已有相關(guān)研究成果[6～9]，而三次拋物線形斷面渠道收縮水深的計(jì)算研究則相對(duì)較少。

文獻(xiàn)[10]針對(duì)三次拋物線形渠道斷面收縮水深計(jì)算涉及高次方程的求解問(wèn)題，棄用常規(guī)的計(jì)算過(guò)程繁復(fù)且成果精度不高的試算法及圖解法，從三次拋物線形斷面渠道收縮水深的基本方程入手，通過(guò)對(duì)初始迭代值的優(yōu)化擬合，獲得了可直接完成求解的近似計(jì)算公式，具有一定的實(shí)際意義。但因該公式需分別完成中間參數(shù)及初值計(jì)算，求解過(guò)程仍顯繁瑣。

為了進(jìn)一步簡(jiǎn)化三次拋物線形斷面收縮水深的計(jì)算過(guò)程，本文采用優(yōu)化擬合的方法，以標(biāo)準(zhǔn)剩余差最小為目標(biāo)函數(shù)，獲得了一種表達(dá)式較為簡(jiǎn)捷、計(jì)算精度較高的近似公式。

1 收縮水深的基本計(jì)算公式

收縮水深的基本方程為[11]：

式中：E0為以收縮斷面底部為基準(zhǔn)面的過(guò)水建筑物上游總水頭，m；hc為收縮斷面處的水深，m；Q為過(guò)水流量，m3/s；g為重力加速度，m/s2；φ為流速系數(shù)；Ac為三次拋物線形斷面面積，m2。

三次拋物線形斷面的曲線方程為：

其過(guò)水?dāng)嗝婷娣e為：

將式（3）、（4）代人式（1），并設(shè)：

經(jīng)整理即可獲得計(jì)算三次拋物線形斷面收縮水深的計(jì)算公式為：

2 近似計(jì)算公式的建立及精度分析

2.1 公式建立

式（6）為高次函數(shù)的超越方程，無(wú)法直接獲解。為避免利用式（6）的超越方程求解問(wèn)題，現(xiàn)通過(guò)以下方法尋求其替代函數(shù)：

(1) 根據(jù)式（6）函數(shù)，展繪α～k關(guān)系曲線。

(2) 經(jīng)對(duì)α～k關(guān)系曲線的線形分析，初步擬定以下函數(shù)為備選替代函數(shù)。

式中：A、B、C、x、β分別為待定系數(shù)及指數(shù)。

(3) 在工程實(shí)用范圍內(nèi)[10]（即 0.01≤α≤0.5），采用優(yōu)化擬合的方法，以標(biāo)準(zhǔn)剩余差最小為目標(biāo)函數(shù)[12]即：

經(jīng)對(duì)式（7）～（10）逐次逼近擬合[13]可知，式（10）的標(biāo)準(zhǔn)剩余差S最小，其替代函數(shù)式為

2.2 精度分析及比較

為比較式（11）與式（6）的擬合精度，考慮在工程實(shí)用范圍內(nèi)（即 0.01≤α≤0.5， 2.59≤k≤100.38），取不同的iα值即可由式（6）分別計(jì)算出與之相對(duì)應(yīng)的ki，再將ki代入式（11）求得與之相對(duì)應(yīng)的αi′，并由下式完成式（11）替代式（6）的擬合相對(duì)誤差，結(jié)果見(jiàn)表1。

表1 式（11）替代式（6）誤差計(jì)算結(jié)果

由表1可見(jiàn)，在工程實(shí)用范圍內(nèi)。

用式（11）替代式（6）的最大擬合相對(duì)誤差：

式（11）具有較好的擬合替代精度。

表2 拋物線均勻流水深公式形式、最大相對(duì)誤差比較

通過(guò)對(duì)公式形式比較可見(jiàn)，本文公式較文獻(xiàn)[5]公式更加簡(jiǎn)單，計(jì)算過(guò)程也更加簡(jiǎn)捷。具體比較結(jié)果見(jiàn)表2所示。

3 應(yīng)用舉例

選文獻(xiàn)[5]計(jì)算實(shí)例：已知閘前斷面總水頭E0=15m，通過(guò)流量Q =162 m3/s，流速系數(shù)φ=0.95，若采用三次拋物線形斷面渠道，其斷面曲線方程為

求閘后斷面收縮水深hc。

將已知參數(shù)代人式（5）可求得

將k =5.4306代人式（11）即可求得

本例收縮水深的精確解為hc=3.00m，本文公式計(jì)算成果的相對(duì)誤差為0。

4 結(jié) 語(yǔ)

通過(guò)優(yōu)化擬合的方法完成了對(duì)三次拋物線形渠道斷面收縮水深計(jì)算公式的擬合，獲得了表達(dá)形式簡(jiǎn)單、求解精度較高的計(jì)算公式，經(jīng)誤差分析及實(shí)例計(jì)算表明，在工程實(shí)用范圍內(nèi)，本文公式的最大相對(duì)誤差為0.46% 。

[1]魏文禮,楊國(guó)麗.立方拋物線渠道水力最優(yōu)斷面的計(jì)算[J].武漢大學(xué)學(xué)報(bào)：工學(xué)版,2006,(3)∶49-51

[2]張志昌,劉亞菲,劉松艦.拋物線形渠道水力最優(yōu)斷面的計(jì)算[J].西安理工大學(xué)學(xué)報(bào),2002,18（3）：235-237

[3]明萬(wàn)才,黃開路,張曉蓮.立方拋物線形斷面明渠水力計(jì)算探討[J].水利科技與經(jīng)濟(jì),2002,8（2）：74

[4]文輝,李鳳玲.立方拋物線形渠道水力計(jì)算的顯式計(jì)算式[J].人民黃河,2010，（1）：75-76

[5]趙延風(fēng),王中正,方興等.半立方拋物線形渠道正常水深算法[J].排灌機(jī)械工程學(xué)報(bào),2011,2（3）：241-245

[6]文輝,李風(fēng)玲.立方拋物線斷面渠道收縮水深的直接計(jì)算方法[J].人民長(zhǎng)江,2009,40（13）：58-59

[7]文輝,李風(fēng)玲.拋物線形斷面渠道收縮水深的解析解[J].長(zhǎng)江科學(xué)院院報(bào),2009,26（9）：32-33

[8]蘆琴,王正中,任武剛.拋物線形渠道收縮水深簡(jiǎn)捷計(jì)算公式[J].干旱地區(qū)農(nóng)業(yè)研究,2007,25（2）∶134-136

[9]趙延風(fēng),宋松柏,孟秦倩.拋物線形斷面渠道收縮水深的直接計(jì)算方法[J].水利水電技術(shù),2008,39（3）∶36-37

[10]冷暢儉，王中正.三次拋物線形渠道斷面收縮水深的計(jì)算公式[J].長(zhǎng)江科學(xué)院院報(bào),2011,28（4）∶29-31,35

[11]清華大學(xué)主編.水力學(xué)（修訂本）上冊(cè)[M].北京∶清華大學(xué)出版社,1990

[12]王慧文.偏最小二乘回歸法及其應(yīng)用[M].北京國(guó)防工業(yè)出版社,1999

[13]閻鳳文.測(cè)量數(shù)據(jù)處理方法[M].原子能出版社,1988