閘后擴散式消力池池深的簡化計算法

滕　凱

(齊齊哈爾市水務局，黑龍江 齊齊哈爾 161006)

?

閘后擴散式消力池池深的簡化計算法

滕凱

(齊齊哈爾市水務局，黑龍江 齊齊哈爾 161006)

摘要：閘后擴散式消力池池深計算需聯立求解高次方程，常規的數學方法不能直接獲得。現有的圖解法及試算法或是求解過程依賴圖表精度差，或是計算過程繁瑣工作效率低，而利用微機編程獲解又不便基層工程技術人員應用。經對池深水力計算相關公式的整理變換，采用優化擬合的方法，完成了對無量綱求解方程的擬合替代，在工程實用參數范圍內，獲得了可直接完成擴散式消力池池深求解的簡化計算公式，計算過程簡捷直觀、成果精度滿足工程設計要求。為此，文章以實例加以論證。

關鍵詞：擴散式消力池；池深；優化擬合；簡捷計算；實例

池式消能工是閘后消能設計選取的主要結構形式，因此《水閘設計規范》通過附錄的形式專門給出了池式消能工關于池深的計算公式。由于該種消能工的水力計算涉及求解繁復的超越方程[1-3]，無法通過常規的數學方法完成求解，為了簡化求解過程，相關學者開展了大量的研究工作，并獲得了許多有益的研究成果。文獻[4-10]分別給出了當b1/b2=1(b1、b2分別為消力池首及末端的寬度)情況下借助圖表試算、直接試算、簡化計算及智能求解的池深計算方法。但在實際水閘工程中，由于閘室的總布置寬度常常小于閘后渠(或河)道的底寬，為了獲得更好的消能效果及與下游渠(或河)道的水流銜接，工程上多將閘后消力池設計為擴散式，而有關擴散式消力池池深計算方法的研究則相對較少[11]，為此，文獻[11]通過建立動量方程給出了擴散式消力池池深的試算獲解方法，計算過程較繁復，不便應用。

1基本計算公式

底流消能的挖深式消力池如圖1所示，當在池中發生完整的淹沒式水躍，并與下游渠道(或河道)水面形成自然銜接時，池深d的確定需聯立求解下列方程：

圖1　水閘消力池消能簡圖

(1)

(2)

(3)

(4)

將式(4)代入式(1)減式(3)經整理可得：

(5)

(6)

式中：x為中間變量；Fr為弗汝德數。

由式(6)可得：

(7)

(8)

式中：η為已知參數。

將式(7)、(8)代入式(5)整理得：

(9)

在式(9)中，設：

(10)

根據水力學原理可得：

(11)

則式(10)為：

(12)

式中：hk為臨界水深，m；t0為下游渠道(或河道)護底以上的總水頭，m；E為已知綜合參數。

將式(12)代入式(9)可得：

(13)

在式(13)中，因當過閘流量Q為已知時，H0、q即可根據閘前渠道(或河道)的斷面參數及閘室寬度求得，進而t0、hk及E可求，則式(10)即為含有1個未知數x的超越方程，可通過采用試算或迭代(整理出迭代形式)法可完成求解。

2簡化公式的建立及精度分析

2.1公式的建立

因式(13)為較繁雜的超越方程，無論采用試算法還是迭代法均需借助計算機編程完成，不便實際應用，為此，本文提出以下簡化計算法。

y=x2-1

(14)

將式(14)代入式(13)得：

(15)

式中：y為引入的中間變量。

(a)η=1

(b)η=0.7

由圖2可見，當φ、η為定值情況下，E隨y的增大呈非線性增加，構成的Ei～yi曲線類型具有較好的拋物線關系。為此，提出以下擬合替代備選函數：

Ay2α+Byα+C=E

(16)

Ay2α+Byα+C=InE

(17)

Ay2α+Byα+C=eE

(18)

采用優化擬合的方法，以標準剩余差最小為目標函數，即為：

(19)

式中：n為數值擬合分析計算的數組數。

經逐次逼近擬合比選[15-16〗可知，式(16)擬合效果最好，并可獲得式(15)的最優替代式為：

Ay2/3.37+By1/3.37+C=E

(20)

其中：A=0.131413φ-2-0.025598η-0.67+0.039922；B=1.615939η-0.23-0.259066φ-2-1.956872；C=0.335885φ-2-0.136459η-2.55+0385221。

式中：A、B及C均為中間參數。

由式(19)即可求得y為：

(21)

2.2精度分析

AI2.0將推動人工智能邁向新的階段，即從人工知識表達技術到大數據驅動知識學習技術，人工智能的實現方法正轉變為從大數據中進行知識發現和學習，使得機器學習從表象深入到綜合推理；從聚焦研究“個體智能”到基于互聯網絡的群體智能，形成群智能力的互聯網服務創新（如共享經濟）體系；從只能處理類型單一的數據到能夠綜合視覺、聽覺、文字等多種媒體的語義，人工智能正在邁向跨媒體認知、學習和推理的新高度；從追求“智能機器”和高水平的人機與腦機交互技術，走向人機混合的增強智能；從機器人到自主智能系統。

圖3　計算y值的誤差包絡線

(22)

(23)

(24)

(25)

當z=-0.94%時

(26)

3計算實例

已知某堰高為P=3m(相對于下游渠底)的溢流壩，堰頂部設有閘門控制，過水單寬流量為q=m3/(s·m)，φ=0.90。下游為底寬16m、邊坡比為1∶2.5的梯形斷面渠道，相應于過閘流量的水深為hs=2.75m。閘后設b1=12m、b2=16m的矩形消力池，試計算當閘前水頭4.0m時的消力池深d。

由上述已知條件可求得：H0=7.0374m，hk=1.5425m，t0=2.8168m，E=1.3682m。根據式(20)～(22)可分別求得：A=0.1752985，B=-0.6338129，C=0.6359646，則由式(23)可得：

(27)

(28)

根據式(1)～(4)利用計算機編程，通過逐次逼近即可求得池深的精確解為d=0.634m，本文方法求解結果的相對誤差為0.32%。由此可見，本文方法具有較好的求解精度。

4結語

消力池深計算是水閘工程底流消能設計的重要環節，因此，開展有關簡化計算方法的研究對提高設計工作效率及精度具有重要意義。本文通過對整理后的無量綱求解方程的擬合替代，獲得了可直接完成相關參數求解的簡化計算公式，僅借助計算器即可非常簡捷地完成池深求解，計算過程大大簡化。誤差分析及實例計算表明，在工程適用參數范圍內，求解精度完全滿足工程計算要求，具有較好的實用推廣價值。

參考文獻：

[1]華東水利學院主編.水工設計手冊[M]. 北京：水利電力出版社，1983.

[2]武漢水利水電學院水力學教研室主編.水力計算手冊[M]. 北京：水利電力出版社，1983.

[3]SL265—2001水閘設計規范[S].北京：中國水利水電出版社，2001.

[4]程應昌.介紹一種消力池水力計算的圖解法[J].水利科技，1980，(2)：11-15.

[5]辛孝明.設計綜合式消力池的簡便計算[J].山西水利科技，2007，(1)：56-57，,6.

[6]張志軍，陶衛.計算消力池深度的一種新方法[J].中國農村水利水電，2000，(5)：16-17.

[7]滕凱.消力池深的簡化計算法[J].人民長江,2012,43(15)：77-79.

[8]張志昌，李若冰.基于動量方程的挖深式消力池深度計算[J].西北農林科技大學學報(自然科學版)，2012，40(12)：214-218.

[9]劉仲桂，陳振坤，王樂秀. 消力池的水力計算[J].廣西科學院學報，1986，(1)38-44.

[10]葉培聰，許小健.遺傳算法工具箱在消力池水力計算中的應用[J].科學技術與工程，2008，(7)：262-263,267.

[11]張志昌，傅銘煥，趙瑩，等.平底漸擴式消力池深度的計算[J].武漢大學學報(工學版)，2013，46(3)：295-299.

[12]謝成玉,滕凱.三次拋物線形渠道斷面收縮水深的簡化計算[J].南水北調與水利科技,2012,10(2)：148-150.

[13]謝成玉,滕凱.拋物線形斷面渠道均勻流水深的近似計算公式[J].水電能源科學,2012,30(7):94～95,172.

[14]滕凱.半立方拋物線形渠道正常水深的近似計算公式[J].長江科學院院報，2012,29(12)：30-32.

[15]王慧文.偏最小二乘回歸法及其應用[M].北京國防工業出版社，1999.

[16]閻鳳文.測量數據處理方法[M].原子能出版社，1988.

Simplified Calculation Method of Diffusive Stilling Pool Depth behind Gate

TENG Kai

(Qiqihaer Urban Water Affairs Bureau，Qiqihaer 161006，China)

Abstract：To calculate the depth of diffusive stilling pool behind the gate needs to be simultaneous with high order equations，can't be gotten by the conventional mathematical method directly. The existing graphic method and trial method is low in accuracy due to depending on the figures and charts during the solution process，or low in working efficiency because complex calculation process，but solutions obtained by computer process are not convenient for technique personnel working in-sit to apply. Through regulation and transmission for the formulas relevant with pool depth hydraulic calculation，the optional fit method was adopted to replace the fit of dimensionless equation，within the scope of practical project parameters，the simplified calculation formulas to calculate directly the depth of diffusive stilling pool were gotten，the calculation process is simple and intuitive，the result accuracy may meet the engineering design demands. Therefore，the paper demonstrates in example herein.

Key words：diffusive stilling pool；depth of pool；optional fit；simple calculation；example

文章編號：1007-7596(2016)04-0046-05

[收稿日期]2016-03-23

[作者簡介]滕凱(1957-)，男，齊齊哈爾人，高級工程師，現任單位總工程師，從事水利防災減災及工程優化設計研究。先后獲齊齊哈爾市政府市長特別獎一項，齊齊哈爾市政府科技進步一等及三等獎各一項。先后被中共黑龍江省委評為“創業、創新、創優”優秀共產黨員，被省水利廳先后授予“全省水利科技工作先進個人”、“全省水利工程建設先進個人”。截至目前，已在《水利學報》、《水利水電科技進展》、《長江科學院院報》、《水文》、《人民長江》、《人民黃河》、《南水北調與水利科技》、《水電能源科學》、《華北水利水電大學學報》等省部級以上刊物發表論文140余篇。

中圖分類號：TV131.4

文獻標識碼：A