长乐市防洪排涝规划中M IKE11模型参数率定分析

沈明星

（河海大学设计研究院有限公司,江苏 南京 210098）

0 前言

长乐市隶属福建省福州市，东濒东海，北临闽江，城区位于闽江口感潮区、江海交界处，是准半岛型沿海城市。市域内河网密布，平原港道纵横，主要河网有洞江、北洋、南洋三片水系，涝水通过河道、水闸自排至闽江、东海。由于其特殊的地理位置及气候条件，常年受洪、涝、潮、台风威胁，且随着城镇化步伐的加快，洪涝灾害造成的损失日益凸显，近年来多次遭遇严重洪涝灾害，进行全市系统的防洪排涝规划十分必要。本规划中防洪排涝计算采用DHI的MIKE11河网水力分析模型，本文就模型参数率定进行初步分析。

1 模型参数率定

因规划区域内无水文站，缺乏水位监测资料，为率定模型，选取2005年“龙王”台风作为典型降雨，实地调查获得河道洪水淹没水深和淹没时间。考虑到现状长乐市各小流域河道壁面种类及状况具有一定的相似性，故选取上洞江、南洋和北洋三个片区作为模型参数率定实地调查区。

1.1 糙率选取

水动力模型构建需要对主要参数河床阻力（糙率）进行率定，率定的参数用于规划方案继续保留的原河道，对于规划新开的河道或规划重新整治的河道，粗糙系数取0.025~0.030。

1.2 “龙王”台风模拟

1.2.1 计算分区

根据实地调查区现状河道走向、道路分布和地块高程等特点，将洞江片区、南洋片区和北洋片区分别划分为17个、22个和14个计算单元。

1.2.2 降雨量

根据长乐市雨量站的观测资料，2005年10月2号9：00至2005年10月3号8:00“龙王”台风24小时降水量为550 mm，降雨主要集中在18时至21时降水达405mm，占24小时降水量的73.6%，19时出现峰值，时段降水深度190mm。“龙王”台风24小时降水过程见图1。

图1“龙王”台风24小时降水过程图

1.2.3“龙王”台风洪水过程

以洞江片为代表，采用无因次单位线法计算出“龙王”台风洪水过程，“龙王”台风在洞江片区产生的洪水总量约为3458万m3。洞江片“龙王”台风洪水相关参数见表1。

表1 洞江片“龙王”台风洪水相关参数分析

1.2.4 洪水淹没情况调查

长乐市2005年10月2日“龙王”洪水灾害，是长乐市有史以来最严重的一次，受灾面积广、受灾人口达16万人，受淹房屋达2.87万间，倒塌房屋200多间。通过实地访问、调查，获取具可靠性的洪痕共20个。

1.2.5 现状地块蓄量分析

现状河道多为自然河道，岸边不同程度地分布着低洼地和农田，低洼地和农田在洪水期间蓄积一定的水量，降低河道的洪水位，是影响模型模拟水位的重要因素之一。根据万分之一地形图上的地块高程标识，分析计算出不同水位对应的低洼地和农田的蓄水量。洞江片现状地块蓄量分析结果见表2。

表2 洞江片现状地块蓄量分析结果

2 参数率定结果

将“龙王”台风期间洞江片各计算单元洪水过程、边界条件，输入模型中进行分析计算，在可行域内不断调整模型水动力学糙率等参数，直至调查点模拟最高水位、淹没时间与实地调查基本吻合为止。选取洞江片福北线公路桥、塘屿村、丰山村、珠湖村以及大溪村处共5个洪痕调查点进行模拟，经过多次调试最后选定糙率为0.045，各调查点洪水过程模拟结果见图，调查点模拟最高洪水位、周边地块淹没时长与调查对比情况见表3。洞江片区各调查点“龙王”台风洪水过程模拟结果见图2。

图2 洞江片区各调查点“龙王”台风洪水过程模拟结果

表3 调查点最高洪水水位、周边地块淹没时长调查与模拟对比表

对比结果显示，调查点模拟最高洪水位、退水时间基本与调查结果相吻合，表明模型参数已经与实际基本相符，可以用该模型及其参数对其他规划分区进行不同情景条件下防洪排涝系统运行情况的演算与分析。

3 结语

MIKE11模型计算过程中，参数的选定十分关键，率定验证的数据越多，模型就越可靠。在缺乏水文监测资料区域，可选取典型降雨，并通过实地调查获得河道洪水淹没水深和淹没时间，为模型率定提供可靠的验证数据。模型通过不断试算，直至调查点模拟最高水位、淹没时间与实地调查基本吻合，则参数率定基本完成。