張 駿，姬戰(zhàn)生，孫映宏

(1.建德市水文水資源監(jiān)測站，浙江 杭州 310016；2.杭州市水文水資源監(jiān)測總站，浙江 杭州 310016)

?

桐鄉(xiāng)市西部平原河網(wǎng)數(shù)值模擬水源取水可行性研究

張 駿1，姬戰(zhàn)生2，孫映宏2

(1.建德市水文水資源監(jiān)測站，浙江 杭州 310016；2.杭州市水文水資源監(jiān)測總站，浙江 杭州 310016)

利用MIKE11軟件對桐鄉(xiāng)西部平原河網(wǎng)的水流進行模擬分析，并利用流域內(nèi)河網(wǎng)實測水位資料進行驗證.結果表明，計算值與實測值吻合較好，具有較好的精度.運用建立的河網(wǎng)模型，對河網(wǎng)內(nèi)部擬設置取水口處的來水情況進行了計算與分析，結果表明，4種不同的工況條件下，取水口處水量、水位均能滿足工程取水的要求.

平原河網(wǎng)；水動力學模型；取水口；可靠性研究

桐鄉(xiāng)市位于杭嘉湖平原，境內(nèi)地勢平坦，河道縱橫，水網(wǎng)密布，屬太湖水系.近年來，隨著經(jīng)濟社會的發(fā)展和城鄉(xiāng)供水一體化的實施，全市公共水廠的供水量逐年增加，目前全市承擔集中供水的果園橋水廠和運河水廠總供水能力為30萬t/日，2010年供水總量為8 012萬t，比上年增長10.6%，日最高供水量突破28萬t，基本處于滿負荷運行狀態(tài)，2011年日最大供水量達到30.3萬t，已經(jīng)啟用崇福應急水廠，水量供需矛盾日益凸顯.并且，桐鄉(xiāng)市地表水污染不斷加劇，大部分河湖塘漾的水質均處于Ⅴ類和劣Ⅴ類，已很難找到符合飲用水原水標準的水源.為此，桐鄉(xiāng)市積極尋求境內(nèi)新的取水水源地，擬在市區(qū)西部河網(wǎng)中設置取水口，開展集中供水水廠水源從運河向西部水域轉移的可行性研究.

本文針對桐鄉(xiāng)西部水域內(nèi)河道復雜，流向順逆不定的特點，利用丹麥DHI公司的Mike11數(shù)值模擬軟件[1]建立了桐鄉(xiāng)西部水域的河網(wǎng)水動力模型，在各種不同工況條件下，通過模型計算，對區(qū)域內(nèi)取水口附近的水位、流量以及來水方向進行了研究.

1 區(qū)域概況

桐鄉(xiāng)市屬太湖流域運河水系，河道稠密，是典型的江南水鄉(xiāng).河道總長2 398.3 km，河道面積46.3 km2，河道密度3.3 km/km2，河網(wǎng)率6.4%.主要河道長440 km有余.京杭古運河斜貫全境，是全市水利、水運的大動脈，東部主要河道以南北向為主，西部河道以東西向為主.

桐鄉(xiāng)西部水域位于洲泉鎮(zhèn)眾安村與德清交界附近的野凌灘漾至桑柴湖漾、白蕩漾一線水域.分析區(qū)域內(nèi)主要河流有六縱八橫，六縱分別為龍溪港(亦稱老龍溪)、十字港、京杭運河、大東港、三洞環(huán)橋港、京杭古運河(桐鄉(xiāng)段)；八橫為漾溪港、含山塘、西圣埭港、橫塘港、大羔羊港、祝香橋港、大紅橋港、京杭古運河(余杭段)，主要湖漾有百畝漾、野凌灘漾、桑柴湖漾.分析區(qū)域周邊水文站有新市水位站、臨平下水位站、塘棲水位站和崇德水位站.區(qū)域主要河道及水文測站分布(見圖1).

圖1 桐鄉(xiāng)西部水域主要河道及水文測站分布圖

2 一維水動力模型

Mike11河流模型軟件是一個結果清晰、界面友好的模擬系統(tǒng)，廣泛應用于模擬河口、河流、河網(wǎng)、灌溉系統(tǒng)的水流、水質、泥沙輸運等一維問題的專業(yè)軟件包，適合于包括復雜平原河網(wǎng)在內(nèi)的一維非恒定流計算[2-3].

2.1 水動力方程

河網(wǎng)水動力模型控制方程為Saint-Venant方程組：

(1)

(2)

式中：x和t—分別為空間坐標和時間坐標；

B—河寬；

Z—水位；

A—過水斷面面積；

Q—斷面流量；

q—旁側入流、出流流量；

R—水力半徑；

g—重力加速度；

A—動量校正系數(shù)；

C—謝才系數(shù).

Mike11在求解上述方程時，利用Abbott六點隱式格式離散上述控制方程組，該離散格式在每一個網(wǎng)格點并不同時計算水位和流量，而是按順序交替計算水位或流量.該格式無條件穩(wěn)定，可以在相當大的Courant數(shù)下保持計算穩(wěn)定，可以取較長的時間步長以節(jié)省計算時間[4-5].

2.2 河網(wǎng)概化

桐鄉(xiāng)西部水域一維河網(wǎng)數(shù)學模型的計算范圍：北以漾溪港和無量橋港為界；西以老龍溪和東大河為界；東、南以京杭古運河為界.

由于區(qū)域內(nèi)河網(wǎng)河道縱橫交織，水系呈網(wǎng)狀，流向復雜，因此必須對河網(wǎng)進行概化.河網(wǎng)概化的原則是：概化所得的河網(wǎng)必須滿足與天然河網(wǎng)輸水能力和調蓄能力相似、水面率相近.本次模擬概化過程中側重考慮對水流輸送占主導作用的骨干河道，也考慮了少數(shù)輸水能力較小但溝通水系的河道.模擬的主要河道有：老龍溪、東大河、十字港、京杭運河、大東港、白米塘港、高橋港、漾溪港、長洛河、德新線、德博港、橫塘港、三洞環(huán)橋港、大有橋港、無量橋港、西圣埭港、大羔羊港、小羔羊港、祝香橋港、大紅橋港、京杭古運河等.概化后的河網(wǎng)(見圖2)

圖2 桐鄉(xiāng)西部水域河網(wǎng)概化示意圖

2.3 邊界條件

桐鄉(xiāng)西部水域一維河網(wǎng)數(shù)學模型共選用了14個邊界，利用2011年1月1日至8月31日新市、塘棲、臨平下、崇福4個測站實測水位來控制河網(wǎng)邊界水位，以模擬研究范圍內(nèi)河網(wǎng)水體流動情況.

2.4 模型率定

模型對計算斷面對應的河段糙率系數(shù)進行了率定，水動力模型河道糙率取值基本上在0.015～0.03之間.

由于桐鄉(xiāng)西部水域研究范圍內(nèi)部沒有水文測站，為滿足數(shù)學模型驗證要求，在河網(wǎng)內(nèi)部德博港上設置劉皇橋臨時水位站(見圖1)，該站于2011年5月16日開始觀測，同年8月31日結束觀測.選擇該時間段內(nèi)劉皇橋逐日平均水位與模型計算結果進行對比，結果見圖3.從逐日平均水位的率定結果可以看出，劉皇橋臨時水位站點模型計算的逐日水位變化曲線與實測逐日水位變化曲線擬合良好，表明模型可以滿足計算要求.

圖3 劉皇橋逐日平均水位率定結果

3 計算結果與分析

3.1 計算工況

根據(jù)桐鄉(xiāng)西部水域2011年1月1日至8月31日的水文資料，選擇高水、中水、低水三種工況對擬取水點處來水進行模擬計算.同時為充分考慮枯水年份對取水點的影響，選擇2003年作為典型年，計算各邊界水位處于最低狀況時的取水點來水情況.

方案1 高水位工況，以2011年6月份為計算時段，選取6月19日作為該方案的代表時段，各邊界站點的水位取同一日的日平均水位，其中新市站2.37m，崇福站2.42m，塘棲站2.64m，臨平下2.59m.

方案2 中水位工況，以2011年8月份為計算時段，選取8月20日作為該方案的代表時段，各邊界站點的水位取同一日的日平均水位，其中新市站1.56m，崇福站1.49m，塘棲站1.68m，臨平下1.62m.

方案3 低水位工況，以2011年4月份為計算時段，選取4月30日作為該方案的代表時段，各邊界站點的水位取同一日的日平均水位，其中新市站0.78m，崇福站0.75m，塘棲站0.90m，臨平下0.90m.

方案4 2003年典型年工況，分析2003年各邊界站最低逐日平均水位，選取2003年6月25日作為該方案的代表時段，各邊界站點的水位取同一日的日平均水位，其中新市站0.76m，崇福站0.77m，塘棲站0.78m，臨平下0.72m.

3.2 四種工況計算成果與分析

根據(jù)四種計算方案和已建立的河網(wǎng)模型，對桐鄉(xiāng)西片水源擬取水點處按取水規(guī)模0萬t/d、15萬t/d、30萬t/d、45萬t/d分別進行模擬計算，分析不同情況下取水口周邊河道水流的流向及流量，結果(見表1).

表1 桐鄉(xiāng)西部水域不同工況下補給水源計算成果

(1)高水位工況

方案1在4種不同取水規(guī)模下取水點處水位基本不變，均處在2.46m左右，河網(wǎng)內(nèi)部東西向河道水流流向自西向東，南北向河道水流流向由南向北.由表1可知，高水位工況下取水點附近水量主要來自兩條河道，分別是德博港和大東港，其中德博港自西向東流，大東港由南向北而來.從水量分析，不同取水規(guī)模下取水點附近總的流量在73.4m3/s～74.1m3/s之間，其中來自大東港水量的比例在83.8%～84.0%之間，德博港水量的比例在16.0%～16.2%.

(2)中水位工況

方案2在4種不同取水規(guī)模下取水點處水位基本不變，均處在1.56m左右，河網(wǎng)內(nèi)部東西向河道水流流向自西向東，南北向河道水流流向由南向北.由表1可知，中水位工況與高水為工況類似，取水點附近水量主要也來德博港和大東港，其中德博港自西向東流，大東港由南向北而來.從水量分析，不同取水規(guī)模下取水點附近總的流量在27.3m3/s～27.8m3/s之間，其中來自大東港水量的比例在76.6%～76.9%之間，德博港水量的比例在23.1%～23.4%之間.

(3)低水位工況

方案3在4種不同取水規(guī)模下取水點處水位基本不變，均處在0.80m左右，河網(wǎng)內(nèi)部東西向河道水流流向自西向東，南北向河道水流流向由南向北.由表1可知，低水工況也前面兩種工況類似，取水點附近水量主要也來德博港和大東港，其中德博港自西向東流，大東港由南向北而來.從水量分析，不同取水規(guī)模下取水點附近總的流量在19.0m3/s～19.5m3/s之間，其中來自大東港水量的比例在82.1%～82.4%之間，德博港水量的比例在17.6%～17.9%之間.

(4)2003年典型年工況

方案4在4種不同取水規(guī)模下取水點處水位基本不變，均處在0.75m左右，與前3種工況相比，典型年工況下河網(wǎng)水流流向相對復雜，南北向河道水流流向基本自北向南，東西向河道的水流均匯入于南北向河道.由表1可知，典型年工況下，當不取水和取水規(guī)模為15萬t/d時，取水點附近水量來自大有橋港和大紅橋港；當取水規(guī)模為30萬噸/日和45萬t/d時，取水點附近水量來自大有橋港、大紅橋港和三環(huán)洞橋港.從流向分析，大有橋港和三環(huán)洞橋港均從北向南流，大紅橋港自東向西流.從水量分析，不同取水規(guī)模下取水點附近總的流量在1.0m3/s～5.2m3/s之間，當取水規(guī)模為30 萬t/d和45 萬t/d時，大有橋港水量的比例在10%～45%之間，大紅橋港水量的比例在55%～90.%之間；當取水規(guī)模為30 萬t/d和45 萬t/d時，三環(huán)洞橋港水量的比例在17.1%～26.9%之間，大有橋港水量的比例在44.2%～45.7%之間，大紅橋港水量的比例在28.9%～37.2%之間.

綜上所述，在四種不同方案條件下，模型計算的結果表明在不同取水規(guī)模下，取水口處的水位變化不大，均可以滿足水廠取水頭埠高程的要求.前三種方案取水口附近水量主要來自南部余杭方向來水，由南向北經(jīng)由大東港到達取水口.2003年典型年工況下，取水口附近水量主要來自北部“引江濟太”太湖來水，由北向南經(jīng)由三環(huán)洞橋港、大有橋港和大紅橋港達到取水口.

4 結 語

桐鄉(xiāng)西部水域河道密布，縱橫交錯，十分復雜，利用Mike11模型軟件，建立了適合該區(qū)域的河網(wǎng)模型，率定結果表明，計算值和實測值吻合良好，計算精度可以滿足實際應用要求.運用建立的河網(wǎng)模型，對河網(wǎng)內(nèi)部擬設置取水口處的來水情況進行了計算與分析，結果表明，4種不同的工況條件下，取水口處水量、水位均能滿足工程取水的要求，其中前3種方案在不同取水規(guī)模下，取水口處的水位變化不大，取水口附近水量主要來自南部余杭方向來水，由南向北經(jīng)由大東港到達取水口處；2003年典型年工況下，取水口附近水量主要來自北部“引江濟太”太湖來水，由北向南經(jīng)由三環(huán)洞橋港和大紅橋港匯集至取水口處.

[1]DanishHydraulicInstitute(DHI).MIKE11:AModellingSystemforRiversandChannelsReferenceManual[R].Compenhagen:DHI.2007.

[2] 徐祖信,李雨芯,廖振良,等.蘇州河綜合整治決策系統(tǒng)的模型系統(tǒng)開發(fā)[J].中國給水排水,2008,24(9):72-74.

[3] 徐海波,蘇 飛,王士武.余姚市河網(wǎng)平枯水期水質模擬研究[J].浙江水利科技,2008(1):27-29,34.

[4] 吳挺峰,周 鍔,崔廣柏,等.河網(wǎng)概化密度對河網(wǎng)水量水質模型的影響研究[J].人民黃河,2006,28(3):46-48.

[5] 楊 洵,梁國華,周惠成.基于MIKE11的太子河觀葠河段水文水動力模型研究[J].水電能源科學,2010,28(11):84-87,171.

Reliability Study on Water Intake in Plain River
Networks Based on Numerical Simulation

ZHANG Jun1, JI Zhan-sheng2, SUN Ying-hong2

(1.Jiande Hydrology and Water Resources Monitoring Station, Hangzhou 310016, China;2.Hangzhou Hydrology and Water Resources Monitoring Central Station, Hangzhou 310016, China)

In this paper, the software named MIKE11 is used to simulate unsteady flow of western plain river networks in Tongxiang, which is verified by the measured water-level data. The results show that the calculated results agree well with the measured values. Using the established model, the inflow of river water intake is calculated and analyzed. The results show that under four different working conditions, the water flow as well as the water level in water intake can meet the requirements of the project.

plain river network; hydrodynamic model; water intake; reliability research

2015-05-10

國家自然科學基金資助項目(51109054)

張 駿(1962-)，男，浙江建德人，工程師，從事水文水資源保護管理工作.

TV131.1

A

1008-536X(2015)09-0028-05