臺州市溫黃平原地區河網水文水動力模型研究及應用

陳 雰，黃 成

（1.臺州學院，浙江 臺州 318000；2.浙江昊信工程咨詢有限公司，浙江 臺州 318000）

沿海感潮城市如臺州、溫州等，是經濟發展的先驅、濱海城市的明珠。由于所處的特定地理位置、氣候條件、地形特點，在盡享水資源帶來的水土膏膚和舟楫之利的同時，也深受水災的威脅。與大流域地區、內陸地區或其他城市的水災相比，其復雜性和特殊性主要表現以下兩個方面：（1）多致災因子作用，澇災危險性增大，一般情況下，由單個因子造成的水災并不嚴重，而多個因子相互影響、共同作用，往往造成巨大的損失，如熱帶氣旋發生時常伴有暴雨，如適遇天文大潮即“三碰頭”，就有可能造成嚴重的水災[1-2]。（2）下墊面環境復雜，澇災承災體脆弱性日趨增加，平原河網調蓄能力降低、地面沉降嚴重等因素加劇了澇災的形成。以上這兩種負反饋作用使得沿海感潮地區成為受洪澇災害影響最為嚴重、風險最高的區域之一。頻繁的澇災不僅會直接造成財產的重大損失，嚴重影響其經濟發展，還會對社會、環境和人們心理帶來十分不利的影響。同時沿海城市人口密集、經濟發達，對水災信息采集、分析、預報的準確性和反應速度都提出了更高的要求[3]。

此次研究區域臺州市溫黃平原區域位于浙江省中部沿海，屬于典型的沿海感潮區域。洪水受上游來水、當地降雨、下游頂托(如河口潮汐頂托、臺風增水頂托等)和人類活動等綜合影響[3]，洪水運動情況復雜，洪水下泄不暢。研究成果可以為沿海感潮城市暴雨內澇減災防災以及城市暴雨內澇風險的管理提供科學依據。

1 概況

此次研究區域臺州市溫黃平原區域位于浙江省中部沿海，介于東經121°2′～121°40′，北緯28°15′～28°44′，東瀕東海，西連溫州市，南臨玉環市，北接椒江、臨海市、黃巖區。溫黃平原為臺州市最大的平原，也是浙江省的四大平原之一。域內水系主要為金清水系河網，部分為永寧江水系河網。降雨時空分布不均，年際變化較大，年均降雨量1250～2200mm，多年平均降雨量1639.8mm。年降水日數140～180d。

臺州市溫黃平原水文水動力模型包括椒南、路橋、黃巖東部、溫嶺4個市（區），東西長約61km，南北寬約53km，研究區域土地總面積1976.5km2。溫黃平原為臺州市最大的平原，也是浙江省的四大平原之一。域內水系主要為金清水系河網，部分為永寧江水系河網。

2 建模原理

根據流域的地理特性，可以將溫黃平原河網流域大部分歸為平原區，西部部分區域歸為山丘區。洪水受上游來水、當地降雨、下游頂托(如河口潮汐頂托、臺風增水頂托等)和人類活動等綜合影響，洪水運動情況復雜，洪水下泄不暢，這部分區域的洪水模擬需采用以水文、水動力學相結合的方法解決[4]。

溫黃平原河網區的洪水模擬采用水文學和水力學相結合的方法。根據平原地區流域的特征，從流域洪水模擬的角度，流域由河道、湖泊、水庫、蓄洪區、行洪區、圩區及工程建筑物(如閘、堤、大壩、鐵路等)等組成，在平水期，水流主要在河道、湖泊中運動；洪水期，隨著水位的上升，上游來流超過了河道泄流能力時，通過分洪、行洪、蓄洪、圩區等措施，使水流通過閘、口門、漫堤等水工建筑物進入行蓄洪區，而行洪區內的工程建筑物如鐵路、公路又影響著水流運動，模擬流域洪水的目的在于正確地模擬上述水流的運動情況。

2.1 水文模型

模型中山丘區水文模型采用分布式三水源新安江模型。在三水源新安江模型中流域以雨量站控制面積為單元，把流域劃分為與雨量站個數相同的子流域（或單元），在每個子流域中，把雨量站的點雨量作為子流域的平均雨量，計算相應的產流。模型中溫黃平原區域西部地區包括山丘區地形，按照流域（根據所收集小流域邊界shapefile文件）劃分成小區域，生成山丘區水文模型并進行設置，共計山丘區水文模型227個，總面積673.8km2；模型中東部為平原地區生成平原區預報要素并進行設置，共計平原區預報要素42個，總面積1302.7km2。平原地區水文模型統計單位以各個鄉鎮為統計單元設置平原區預報要素，根據各個要素的雨量站（根據泰森多邊形生成權重信息）、下墊面等信息進行設置[4]。

2.2 產匯流模型

由于不同下墊面具有不同的產流規律，該流域下墊面分成四類：水面、水田、旱地和城鎮道路。流域水流運動模擬由零維、一維、二維模擬所組成，各部分模擬必須耦合聯立才能求解，各部分模擬的耦合是通過“聯系”來實現的?！奥撓怠本褪歉鞣N模擬區域的連接關系，主要是指流域中控制水流運動的堰、閘、泵等，“聯系”的過流流量可以用水力學的方法來模擬。

2.3 水文水動力學模型耦合

此次模型開發平原區坡面匯流單位線模型，通過GIS地圖統計河網密度，計算水文分區河道匯流長度，率定坡面匯流速度，統計不同時段匯入計算分區周邊河道柵格的概率，從而確定平原區分布式坡面匯流單位線。

通過平原區分布式單位線匯流單位線，將平原產流的結果匯入河網水動力模型中的河道斷面，進而實現產水模型與水動力模型的耦合。

3 計算結果分析

由于松門、金清閘、西江閘等水位站監測數據受潮位和水閘控制工況等影響較大，因此，溫黃平原地區模型共選取大溪、溫嶺、路橋、澤國、椒江石柱和麻車橋6個水位代表站進行率定，分片的起調水位參考實測水位進行設定。將模型計算的區域主要代表站的水位與實測水位進行對比分析，包括水位過程、水位峰值及出現時間等，水位代表站基面統一采用1985國家高程基準[5-6]。代表站計算與實測洪峰水位對比見表1，過程對比見圖1。由圖1、表1可知，2010年各計算洪峰水位與實測洪峰水位總體誤差較小，誤差在10cm以下。

圖1 溫黃平原2010年7月溫嶺站計算水位過程對比

4 結論

文章針對溫黃平原地區河網地區的特征，結合地區水系分布特點，構建臺州市溫黃平原河網水文水動力模型，并對模型進行率定和驗證計算，主要結論如下：（1）通過對2010年溫黃平原地區的汛期水情進行模擬演算，計算地區各代表站的洪峰水位及水位過程線，與實際情況相差較小，水位過程線趨勢與實況擬合較好，模型模擬結果合理，基本反映了溫黃平原地區汛期的水流特點和實際狀況。（2）模型率定驗證成果表明，臺州市溫黃平原水文水動力模型能夠真實模擬自然情況下的實際水流運動，并具有較高的模擬精度。由于模型的精度受閘泵工程的調度等人類活動影響，因此必須掌握準確的閘泵工程調度運行情況，這樣模型才能取得更好的模擬精度。

表1 2010年溫黃平原代表站洪峰水位計算值與實測值比較表