太湖流域淀山湖水位代表站選擇分析

程公德

摘要：為準(zhǔn)確計(jì)算淀山湖水位，根據(jù)環(huán)淀山湖水文資料，分別采用算術(shù)平均法和加權(quán)平均法，計(jì)算了環(huán)湖13個(gè)斷面的平均值，作為淀山湖平均水位；選擇部分?jǐn)嗝娼M合的水位平均值作為淀山湖代表站水位，并分析比較各個(gè)組合方案的離差程度，以確定最優(yōu)方案。結(jié)合淀山湖實(shí)際情況和測(cè)站現(xiàn)狀，選擇2021年全年小時(shí)水位數(shù)據(jù)進(jìn)行驗(yàn)證，結(jié)果表明：千燈浦、商榻、淀浦河西閘（外）3個(gè)站加權(quán)平均水位精度最高，水位數(shù)據(jù)權(quán)重分別為0.16，0.46和0.38。因此，建議將這3站的水位加權(quán)平均值作為淀山湖代表站水位值。

關(guān)鍵詞：平均水位；水位代表站；算術(shù)平均；加權(quán)平均；淀山湖

中圖法分類號(hào)：P332.3

文獻(xiàn)標(biāo)志碼：A

DOI：10.15974/j.cnki.slsdkb.2023.06.015

文章編號(hào)：1006-0081（2023）06-0076-04

0 引 言

湖泊水位是流域和區(qū)域防洪、供水安全、水生態(tài)及水資源調(diào)度的重要指標(biāo)，一般采用多個(gè)水位站的水位算術(shù)平均值或加權(quán)平均值表示。目前，采用大浦口、望亭（太）、西山、小梅口、夾浦等沿湖5站的水位算術(shù)平均值來(lái)代表太湖水位［1］。淀山湖作為長(zhǎng)三角生態(tài)綠色一體化發(fā)展示范區(qū)內(nèi)的重要湖泊，是太湖流域第二大省際湖泊，也屬于全國(guó)重點(diǎn)河湖生態(tài)流量監(jiān)測(cè)與信息報(bào)送斷面名錄內(nèi)容。淀山湖湖區(qū)水域面積不大，但受潮汐、臺(tái)風(fēng)、出入流、蓄水變化等因素影響，環(huán)湖各站水位相差較大，尤其是大潮、風(fēng)暴潮、暴雨、洪水期間的上下游站點(diǎn)。淀山湖水位是衡量淀山湖水域水情最重要的指標(biāo)。確定淀山湖湖區(qū)水位時(shí)，常需依賴個(gè)別水位站站點(diǎn)的水位，如西部入湖河道急水港上的商榻站水位。但單一站點(diǎn)的水位難以反映整個(gè)湖區(qū)的水位變化情況，因此，需對(duì)淀山湖水位代表站進(jìn)一步研究。

1 淀山湖概況

淀山湖位于太湖流域蘇滬邊界地區(qū)，分屬上海市青浦區(qū)和江蘇省蘇州市，湖區(qū)面積59.2 km2，屬吞吐型弱感潮平原淺水湖泊，是太湖流域重要的洪水下泄通道、水資源蓄集涵養(yǎng)地以及上海市重要濱水生態(tài)水域［2-3］。環(huán)湖河網(wǎng)密布，進(jìn)、出湖的河道多達(dá)60余條，主要有商榻、珠砂港、攔路港、淀浦河、元蕩、千燈浦等。目前，淀山湖湖區(qū)及環(huán)湖沿線共布設(shè)水位站8處，湖區(qū)為淀中站，環(huán)湖沿線為商榻、千燈浦、珠砂港、淀浦河西閘（外）、葉水路泵閘、石塘港、河祝，平均每7.44 km2有1處水位站。淀山湖湖區(qū)水位因受降水、徑流、潮汐、水利工程調(diào)度、水面蒸發(fā)、風(fēng)速風(fēng)向等因素影響，其水位值在時(shí)間和空間上不斷變化［4-5］，每日水位一般有明顯的兩次漲落變化過(guò)程。

2 研究方法

2.1 水位平均法

在實(shí)際計(jì)算中，湖泊平均水位常采用多站算術(shù)平均或加權(quán)平均計(jì)算，并在滿足計(jì)算精度的前提下盡可能減少水位站點(diǎn)的數(shù)量，這就要求選擇合理的水位代表站。若環(huán)湖水位站點(diǎn)分布比較均勻，可以采用算術(shù)平均方法計(jì)算湖泊平均水位，計(jì)算公式如下：

Z—=1n∑ni=1Zi （1）

式中：Z—為平均水位，m；Zi為第i個(gè)站點(diǎn)水位，m；i=1，2，3，…，n。

若水位站點(diǎn)分布不均勻，可以采用加權(quán)平均方法計(jì)算湖泊平均水位，即各站水位值與對(duì)應(yīng)面積權(quán)重系數(shù)的乘積和［6］，計(jì)算公式如下：

Z—=∑ni=1αiZi（2）

式中：αi為第i個(gè)站點(diǎn)水位的權(quán)重系數(shù)，可采用垂直平分法劃分各站點(diǎn)的控制水域面積，再與湖區(qū)水域總面積計(jì)算得到面積權(quán)重系數(shù)。

考慮到數(shù)據(jù)精度、測(cè)站管理、控制要求等，一般情況下水位站布置在出入湖主要河道、重要工程等附近，以準(zhǔn)確反映河道與湖泊水情規(guī)律。分析計(jì)算中，可以選取多個(gè)環(huán)湖站點(diǎn)的水位數(shù)據(jù)的算術(shù)或加權(quán)平均值作為湖泊平均水位，選擇部分環(huán)湖站點(diǎn)的算術(shù)或加權(quán)平均值作為代表站水位值，比較分析平均水位和代表值之間的誤差，作為衡量水位代表站方案優(yōu)劣的標(biāo)準(zhǔn)［1，7］。

2.2 誤差分析法

環(huán)湖站點(diǎn)可以組合成多種代表站方案，一般情況下站點(diǎn)越多，代表站水位越接近湖泊平均水位。結(jié)合實(shí)際應(yīng)用要求和湖區(qū)水域面積大小，選擇的站點(diǎn)應(yīng)該數(shù)量合適、位置分散分布，且自動(dòng)監(jiān)測(cè)精度較高。判斷水位代表站組合方案的優(yōu)劣一般采用統(tǒng)計(jì)分析方法，即統(tǒng)計(jì)分析代表站水位和平均水位之間的離差情況［8］。常用的指標(biāo)有相關(guān)系數(shù)、標(biāo)準(zhǔn)離差、平均絕對(duì)離差、均值離差、最大離差等［9］。

r=cov（Zs—，Z—）varZs—var（Z—）（3）

ez=（Zs—-Z—）2（4）

ea=1n∑ni=1Zs—-Z—（5）

eo=1n∑ni=1（Zs—-Z—）（6）

Δzm=Max{（Zsi-Zi）}（7）

式中：r為相關(guān)系數(shù)；Zs—為代表站平均水位，m；Z—為平均水位，m；cov（Zs—，Z—）為代表站水位與平均水位的協(xié)方差；var為對(duì)應(yīng)代表站水位和平均水位的方差；ez為標(biāo)準(zhǔn)離差，m；ea為平均絕對(duì)離差，m；eo為均值離差，m；Δzm為最大離差，m。

3 代表站選擇分析

3.1 數(shù)據(jù)資料

數(shù)據(jù)資料原則上應(yīng)選擇淀山湖環(huán)湖各個(gè)口門的長(zhǎng)系列觀測(cè)數(shù)據(jù)，但考慮到自動(dòng)站較少，人工觀測(cè)任務(wù)較重，現(xiàn)有歷史數(shù)據(jù)資料較少，故分析資料選擇短序列和長(zhǎng)序列相結(jié)合方式，其中短序列資料的口門覆蓋了環(huán)湖所有主要出入湖口門，口門數(shù)量較多且分散，站點(diǎn)為自動(dòng)監(jiān)測(cè)站和臨時(shí)布置站，數(shù)據(jù)來(lái)源為水文同步調(diào)查，主要用于關(guān)系分析；長(zhǎng)序列資料的口門僅覆蓋了環(huán)湖最主要出入湖口門，站點(diǎn)均為自動(dòng)站，歷史數(shù)據(jù)齊全，監(jiān)測(cè)頻率較高，主要用于驗(yàn)證分析。

短系列同步觀測(cè)資料主要分析環(huán)湖站點(diǎn)較多情況下的水文站代表關(guān)系，資料為環(huán)淀山湖水文同步調(diào)查的監(jiān)測(cè)成果，涉及主要出入湖河道13個(gè)站點(diǎn)（自動(dòng)站6個(gè)、臨時(shí)布置站7個(gè)），其中入湖站點(diǎn)9個(gè)，商榻、珠砂港、千燈浦、白石磯等；出湖站點(diǎn)4個(gè)，河祝、石塘港、西旺港、淀浦河西閘（外），見圖1。同步調(diào)查時(shí)間2019年5月21～30日，連續(xù)16個(gè)潮流期，水位觀測(cè)15 min一次。

長(zhǎng)系列數(shù)據(jù)資料主要驗(yàn)證代表站在長(zhǎng)時(shí)間水情規(guī)律中的精度情況，資料為環(huán)淀山湖已建自動(dòng)水位站的監(jiān)測(cè)數(shù)據(jù)，涉及主要出入湖河道6個(gè)站點(diǎn)，其中入湖站點(diǎn)3個(gè)，商榻、珠砂港和千燈浦；出湖站點(diǎn)3個(gè)，河祝、石塘港和淀浦河西閘（外）。數(shù)據(jù)為2021年全年小時(shí)水位數(shù)據(jù)。

3.2 統(tǒng)計(jì)分析

根據(jù)短系列同步觀測(cè)資料，采用13個(gè)斷面水位數(shù)據(jù)的算術(shù)或加權(quán)平均值作為淀山湖平均水位，選擇部分?jǐn)嗝娼M合水位數(shù)據(jù)的算術(shù)或加權(quán)平均值作為代表站水位，比較分析平均水位和代表站水位值之間的離差程度，其中代表站組合的選擇有4個(gè)原則：① 代表站位置分散，不能距離太近；② 湖泊有較多出入河道口門，這是水量出入的關(guān)鍵，需選擇出入河道水量較大的站點(diǎn)；③ 站點(diǎn)數(shù)量不能太多，需兼顧計(jì)算精度和便利性；④ 優(yōu)先選擇自動(dòng)監(jiān)測(cè)站點(diǎn)，以實(shí)現(xiàn)實(shí)時(shí)自動(dòng)監(jiān)測(cè)。13個(gè)站點(diǎn)根據(jù)序號(hào)依次為珠砂港（1）、尤車港（2）、蓮湖港（3）、商榻（4）、白石磯（5）、西旺港（6）、石塘港（7）、河祝（8）、淀浦河西閘（外）（9）、斜路港（10）、千燈浦（11）、毛柴涇（12）、后港（13）。選擇滿足要求的站點(diǎn)組合進(jìn)行計(jì)算，得到表1和表2。

3.2.1 算術(shù)平均法計(jì)算結(jié)果

選擇算術(shù)平均法計(jì)算平均水位和代表站水位，并比較離差程度，結(jié)果見表1。算術(shù)平均計(jì)算中，方案1～9，17～18結(jié)果較好，各項(xiàng)離差指標(biāo)均明顯優(yōu)于其他方案。

3.2.2 加權(quán)平均法計(jì)算結(jié)果

若水位站點(diǎn)分布不均勻，可以采用加權(quán)平均法計(jì)算湖泊平均水位，前提是需要計(jì)算各站水位的控制面積，確定相應(yīng)的權(quán)重系數(shù)。選取自動(dòng)監(jiān)測(cè)站點(diǎn)作為組合斷面，利用ArcGIS軟件通過(guò)垂直平分法劃分計(jì)算各站水位權(quán)重系數(shù)，再根據(jù)權(quán)重系數(shù)計(jì)算加權(quán)平均情況，并比較平均水位和代表站水位離差程度，結(jié)果見表2。結(jié)果表明方案1，4，6，8結(jié)果較好，各項(xiàng)離差指標(biāo)均明顯優(yōu)于其他方案。

3.3 結(jié)果分析

采用算數(shù)平均法和加權(quán)平均法計(jì)算淀山湖平均水位（環(huán)湖13個(gè)斷面的平均值），并比較兩組數(shù)據(jù)誤差，結(jié)果表明：相關(guān)系數(shù)0.993、標(biāo)準(zhǔn)離差0.4 cm、平均離差0.3 cm、均值離差-0.3 cm、最大離差1.7 cm、離差大于1 cm組數(shù)7組，說(shuō)明淀山湖湖區(qū)平均水位采用算術(shù)和加權(quán)平均得到的結(jié)果相差不大。

由離差統(tǒng)計(jì)，一般情況下選取站點(diǎn)越多相關(guān)程度越高、標(biāo)準(zhǔn)離差等指標(biāo)越小，數(shù)值與淀山湖平均水位越接近。若僅采用1個(gè)水位站來(lái)代表淀山湖水位，難以反映整個(gè)湖區(qū)水位的變化?？紤]到淀山湖湖區(qū)面積和上下游水情特點(diǎn)，為便于水位代表站實(shí)際應(yīng)用和管理，優(yōu)先選擇自動(dòng)監(jiān)測(cè)站點(diǎn)，在滿足精度要求的同時(shí)盡可能減少站點(diǎn)數(shù)，每個(gè)方案中至少包含2個(gè)水位站，計(jì)算結(jié)果的離差指標(biāo)相對(duì)較小。根據(jù)算術(shù)和加權(quán)平均法的離差統(tǒng)計(jì)結(jié)果，2個(gè)站精度較高組合方案是千燈浦和商榻，3個(gè)站精度較高組合方案是千燈浦、商榻、淀浦河西閘（外）。

3.4 數(shù)據(jù)驗(yàn)證

為驗(yàn)證短系列同步觀測(cè)資料分析得到的代表站平均水位的精度，采用長(zhǎng)系列數(shù)據(jù)資料（全年小時(shí)水位，共計(jì)8 761組數(shù)據(jù)）計(jì)算所選兩種方案的算術(shù)和加權(quán)平均水位，分析千燈浦和商榻2個(gè)站，千燈浦、商榻和淀浦河西閘（外）3個(gè)站代表水位與湖區(qū)平均水位離差情況，得到結(jié)果見表3所示。

根據(jù)長(zhǎng)系列數(shù)據(jù)資料分析情況，千燈浦、商榻和淀浦河西閘（外）3個(gè)站加權(quán)平均水位精度最高，建議選取千燈浦、商榻和淀浦河西閘（外）3個(gè)站加權(quán)平均水位作為淀山湖水位代表站，其中3站權(quán)重分別為0.16，0.46和0.38。

4 結(jié) 論

本文基于淀山湖水情特點(diǎn)，采用算術(shù)平均法和加權(quán)平均法分別計(jì)算了不同代表站組合下的水位，并選取歷史數(shù)據(jù)進(jìn)行驗(yàn)證，得到精度較高且站點(diǎn)數(shù)量合適的水位代表站組合。鑒于淀山湖的東部和南部部分地區(qū)無(wú)口門或者口門較小，考慮到水位站點(diǎn)布置較少，建議在相關(guān)區(qū)域適當(dāng)增設(shè)水位站點(diǎn)，進(jìn)一步分析驗(yàn)證淀山湖水位代表站的關(guān)系。

參考文獻(xiàn)：

［1］ 徐向陽(yáng)，浦寶元.太湖水位分析和代表站選擇［J］.河

海大學(xué)學(xué)報(bào)，1992，20（6）：26-31.

［2］ 易文林，方正杰.淀山湖水量交換分析［J］.上海水務(wù)，2014，30（4）：46-47.

［3］ 王壽兵，馬小雪，張韋倩，等.上海淀山湖水環(huán)境容量評(píng)估［J］.中國(guó)環(huán)境科學(xué)，2013，33（6）：1137-1140.

［4］ 章杭惠，劉曙光，季同德.太湖流域降雨與太湖水位關(guān)系分析［J］.水文，2021，41（1）：85-89.

［5］ 曹中華，夏士淳，劉哲.青年水庫(kù)水位代表性分析及動(dòng)庫(kù)容對(duì)調(diào)洪演算影響的探討［J］.黑龍江水利科技，1999（2）：55-56.

［6］ 王忠紅，趙賢產(chǎn)，余菲.義烏降水氣候分析中的面雨量計(jì)算及其應(yīng)用比較［J］.浙江水利水電學(xué)院學(xué)報(bào)，2010，22（1）：47-50.

［7］ 李衛(wèi)東，閔勇，施順成，等.太湖流域平原水網(wǎng)城市生態(tài)補(bǔ)水措施研究［J］.人民長(zhǎng)江，2022，53（8）：75-81.

［8］ 孫逢立，崔海濤，齊云婷，等.南四湖湖區(qū)水位代表性分析［J］.海洋湖沼通報(bào)，2010（2）：96-100.

［9］ 袁玉，魯軍，胡挺，等.三峽水庫(kù)典型運(yùn)用過(guò)程庫(kù)區(qū)水面線變化特點(diǎn)研究［J］.水利水電快報(bào)，2022，43（7）：89-94.

（編輯：江 文）

Analysis on selection of water level representative stations in Dianshan Lake

of Taihu Basin

CHENG Gongde

（Hydrological Survey Team of Shanghai Qingpu District，Shanghai 200434，China）

Abstract：

In order to accurately determine the water level of Dianshan Lake，on the basis of the hydrological synchronous survey data of Dianshan lake，we calculate the average value of 13 sections as the average water level of Dianshan Lake by arithmetic average method and weighted average method.Also，we selected the average water levels of some sections as the water level of the representative station，and compared the degree of each combination scheme to determine the optimal scheme.Then，we selected 2021 annual hourly water level data to verify the calculation results.The results showed that the weighted average water level accuracy of Qiandengpu Station，Shangta Station and Dianpuhexi gate （outside） Station were the highest，the water level data were 0.16，0.46 and 0.38 respectively.It is suggested to take the weighted average of the water level of the three stations as the representative water level of Dianshan Lake．

Key words：

average water level； representative water level station； arithmetic mean； weighted average； Dianshan Lake