吳立新 孟小敏 韓 陽(yáng)

（江蘇省灌溉總渠管理處 淮安 223200）

1 理論依據(jù)

在水力學(xué)和水文資料整編中，堰閘沉溺式孔流的流量一般用下式進(jìn)行計(jì)算：

Q=MBe △Z0.5

假設(shè)閘門處在即將提出水面而未出水（如距出水相差0.1cm）的臨界狀態(tài)，這時(shí)閘門對(duì)水流的阻礙作用微乎其微，所以此時(shí)的流態(tài)既可以把它看作是孔流，也可以把它視為堰流。因水流具有連續(xù)性，當(dāng)閘門緩慢提升在這兩種流態(tài)之間過(guò)渡時(shí)，流量大小不會(huì)發(fā)生突變，此時(shí)用孔流和堰流公式計(jì)算出的流量應(yīng)該相等。用孔流公式來(lái)推求堰流流量是完全可以的。也可以這樣認(rèn)為，堰流是孔流的一種特殊狀態(tài)，是孔流狀態(tài)下閘門開啟高度最大時(shí)的情形，換句話說(shuō)，就是所有的堰流都可以當(dāng)作孔流來(lái)處理。

2 應(yīng)用實(shí)例

下面以一個(gè)堰閘站為例，根據(jù)江蘇省蘇北灌溉總渠六垛南閘水文站2003—2005年堰流流態(tài)下實(shí)測(cè)流量76 次（見表1），用孔流公式M=Q/（Be△Z0.5）計(jì)算流量系數(shù)。該閘共7 孔，弧形閘門，單寬9.05m，閘底為實(shí)用堰型，閘底高程-0.35m，其中開啟高度e 用上游水頭（上游水位與閘底高程的差值）代替。以e/△Z 為縱坐標(biāo)、M 為橫坐標(biāo)，將76 個(gè)點(diǎn)據(jù)（見表2）點(diǎn)繪在直角坐標(biāo)圖上，如圖1所示。從圖1可以看出，點(diǎn)群自上而下呈帶狀分布，如果從點(diǎn)群中心畫一條垂直于橫坐標(biāo)的直線，則此直線與橫坐標(biāo)交點(diǎn)的M 值為6.48。

圖1 六垛南閘沉溺式堰流e/△Z 的關(guān)系圖

表1 六垛南閘2003—2005年實(shí)測(cè)流量表

表2 六垛南閘流量系數(shù)與水力因子的關(guān)系表

3 定線檢驗(yàn)和均方差計(jì)算

為了驗(yàn)證測(cè)點(diǎn)與關(guān)系線的擬合情況，根據(jù)《水文資料整編規(guī)范》（以下簡(jiǎn)稱《規(guī)范》）要求，對(duì)關(guān)系線分別進(jìn)行符號(hào)檢驗(yàn)、適線檢驗(yàn)和偏離數(shù)值檢驗(yàn)。

3.1 符號(hào)檢驗(yàn)

可用下式進(jìn)行符號(hào)檢驗(yàn)：

u=（k-0.5n-0.5）/（0.5n0.5）

式中：u—符號(hào)檢驗(yàn)統(tǒng)計(jì)量；k—正號(hào)數(shù)；n—測(cè)點(diǎn)總數(shù)。由表3 可知，k=39，n=76。所以：

表3 六垛南閘流量關(guān)系線符號(hào)檢驗(yàn)計(jì)算表

u=（39-0.5×76-0.5）÷（0.5×760.5）≈0.11

當(dāng)α=0.25 時(shí)，由《規(guī)范》查表得u1-α/2=1.15，因u＜u1-α/2，故符號(hào)檢驗(yàn)合格。

為解決一個(gè)問(wèn)題而讀，獲得的是一時(shí)的參考，廣博而有計(jì)劃的系統(tǒng)閱讀是一種長(zhǎng)遠(yuǎn)的眼光，是建立在普遍聯(lián)系與宏觀綜合的視野之上的自我錘煉。完善的知識(shí)結(jié)構(gòu)將最終決定我們?cè)诮逃@條路上能夠走多遠(yuǎn)。

3.2 適線檢驗(yàn)

可用下式進(jìn)行適線檢驗(yàn)：

u=[0.5×（n-1）-k-0.5]÷[0.5×（n-1）0.5]

式中：u—適線檢驗(yàn)統(tǒng)計(jì)量；k—變換符號(hào)次數(shù)；n—測(cè)點(diǎn)總數(shù)。

將e/△Z 從小到大排序（表4），自上而下統(tǒng)計(jì)表中的相對(duì)偏差變號(hào)次數(shù)k=30。

表4 六垛南閘流量關(guān)系線適線檢驗(yàn)計(jì)算表

u=[0.5×（76-1）-30-0.5]÷[0.5×（76-1）0.5]≈1.62

當(dāng)α=0.05 時(shí)，從《規(guī)范》查表得u1-α=1.64。因u ＜u1-α,故適線檢驗(yàn)合格。

3.3 偏離數(shù)值檢驗(yàn)

平均相對(duì)偏離值：

P平均=∑（pi-6.48）÷n=2.2÷76 ≈0.029

式中：pi—測(cè)點(diǎn)與關(guān)系線的相對(duì)偏差；n—測(cè)點(diǎn)總數(shù)；如表5所示。

平均相對(duì)偏離值的標(biāo)準(zhǔn)差：

Sp=｛∑（pi-6.48）2÷[n（n-1）]｝0.5

=｛3622.98÷[76（76-1）]｝0.5≈0.80

式中相對(duì)偏差平方的累計(jì)計(jì)算如表5所示。

偏離數(shù)值檢驗(yàn)公式：

t=P平均÷Sp=0.029÷0.80 ≈0.036

式中：t—偏離數(shù)值統(tǒng)計(jì)量。

當(dāng)α=0.20，n=76 時(shí)，由《規(guī)范》查表得t1-α/2=1.28。因t ＜t1-α/2，故偏離數(shù)值檢驗(yàn)合格。

3.4 均方差計(jì)算

Se=[∑（pi-6.48）2/（n-1）]0.5

=[3622.98/（76-1）]0.5≈6.95%

4 結(jié)論

（1）由3 項(xiàng)檢驗(yàn)結(jié)果可知，推算結(jié)果完全滿足《規(guī)范》定線要求，實(shí)測(cè)流量與計(jì)算流量的均方差為6.95%，達(dá)到二類精度標(biāo)準(zhǔn)。

（2）六垛南閘堰流2003—2005年采用堰流公式定線，由于閘下河床泥沙淤積嚴(yán)重，河床變化較大，關(guān)系線每年都在變化，如使用孔流公式，關(guān)系線較為穩(wěn)定，可達(dá)到較好效果。

（3）用孔流公式推算堰流流量的優(yōu)點(diǎn)是：首先是定線簡(jiǎn)單，平均M 值可直接取各點(diǎn)據(jù)M 值的算術(shù)平均值，而以往堰流流量關(guān)系線為曲線，定線時(shí)有很大的任意性。其次是對(duì)一般的堰閘，可用公式Q=MBe△Z0.5快速地進(jìn)行流量估算，對(duì)平底閘M 值一般取4～6，對(duì)實(shí)用堰閘M 取值要略大些。

（4）堰流狀態(tài)下，閘門是提出水面的，上、下游水位差一般都比較小，小的時(shí)候往往只有3～6cm，目前水位的測(cè)量技術(shù)由于觀測(cè)誤差和水尺零高毫米位的取舍，正常誤差都在1～2cm，這會(huì)造成小水位差點(diǎn)據(jù)散亂較難定線，所以與使用堰流公式定線一樣，選一個(gè)水位差比較大的流量實(shí)測(cè)資料來(lái)分析，或者提高水位差觀測(cè)精度，效果會(huì)更好■