河网地区洪水预报模型开发及在长江下游地区的应用.pdf

Mar.2024 NO.3 VOL.342024年3月 第3期 第34卷CHINA FLOOD&DROUGHT MANAGEMENT中国防汛抗旱特邀主编专栏GUEST EDITOR-IN-CHIEFS COLUMN0 引言长江下游各地区人类活动频繁，下垫面与水网条件复杂1，防洪除涝工程面广、量大，蓄滞洪区、圩区、泵站、调水等工程措施和防洪非工程措施相结合，运行状态多变2。此外，洪涝过程还受本地暴雨、长江过境洪水与沿海风暴潮影响3，诸多因素导致产汇流过程和洪涝灾情复杂，工程群与地区间防洪除涝联合调度难度大。以往研究中，该区域模拟主要采用水动力模型，涉及控制方程组简化、方程组离散和求解、初边条件确定等一系列问题4-5。此类方法需要详细的河道地形资料及特征参数，且计算收敛稳定性受输入边界影响较大，规模庞大、时间长。传统建模及求解方式难以同时满足精准度与时效性的要求，实时预报调度难度极大6-7。全球各地的水文学家开发了许多分布式或半分布式流域水文模型，Freeze和Harlan8首先提出了分布式水文模型设想，最具有代表性的是 SHE（Systme HydrologiqueEuropen）模型9。这类模型从水文循环过程的物理机制入手，将蒸发、下渗、地表水运动、土壤水运动、地下水运动、产汇流等过程联系在一起研究考虑水文变量的空间变异性，通常也称为“白箱”模型。英国Wallingford 水力研究所研发的 InfoWorks10和丹麦水力学研究所（DHI）研发的MIKE 系列软件11可以模拟流域产流、河流水力过程。中国水利水电科学研究院自主研发的城市洪水分析系统IFMSURBAN耦合集成了一维河网、管网模型和二维水动力学模型12，具备模拟城市排水和地表洪水演进的能力。但以上模型缺少耦合流域水文、水动力、工程调度的一体化模型，无法实现人类活动影响下的全流域水循环实时预报调度。基于河海大学太湖流域模型理论及应用实践13，本文河网地区洪水预报模型开发及在长江下游地区的应用陈钢1赵艳红2王晓书2陈秋实1王船海1（1.河海大学水文水资源学院，南京210098；2.沂沭泗水利管理局水文局（信息中心），常州221018）摘要：平原水网地区的洪水预报一直是水文学科的难点问题，尤其受人类活动影响剧烈。为了更好地模拟流域水文循环过程，解决工程实际应用中的挑战，提出了一种新颖的分布式架构水文模型。该模型具有分布式特征，能够根据流域特征和需求灵活选择最合适的水文特征单元进行组合，从而支持对不同水文特征单元精确化、理论化的研究。为了实现实时预报功能，还研发了数据库耦合技术、大型河网实时校正技术及水利工程调度技术。在此基础上，应用分布式架构思想建立了太湖流域模型，对其进行了概化处理，并选取2016年实况数据对模型参数进行了率定。通过对太湖水位1 d、2 d、3 d预见期的实时预报，计算成果基本能够反映平原水网地区太湖流域的水流运动实际状况。综合而言，该分布式架构水文模型方法不仅在科学性上得到验证，而且在实现技术上具有可行性，成功地解决了从理念到现实的应用实践问题，为平原水网地区洪水预报提供了新的思路和方法。关键词：河网地区；洪水预报；水文模型；分布式架构中图分类号：P338；TV87文献标识码：A文章编号：1673-9264（2024）03-08-08DOI:10.16867/j.issn.1673-9264.2024063陈钢，赵艳红，王晓书，等.河网地区洪水预报模型开发及在长江下游地区的应用J.中国防汛抗旱，2024，34（3）：8-15.CHEN Gang,ZHAO Yanhong,WANG Xiaoshu,et al.Development of flood forecasting models in river network areas and the application in the lower reaches of the Yangtze RiverJ.China Flood&Drought Management,2024,34（3）：8-15.（in Chinese）收稿日期：2024-02-28第一作者信息：陈钢，男，副教授，E-mail：。基金项目：国家重点研发计划资助项目（2021YFC3000104）；国家自然科学基金重点项目（U2240209）。8Mar.2024 NO.3 VOL.342024年3月 第3期 第34卷中国防汛抗旱CHINA FLOOD&DROUGHT MANAGEMENT特邀主编专栏GUEST EDITOR-IN-CHIEFS COLUMN提出了一个通用的河网地区洪水预报模型框架，并进行了实例应用。该模型能够根据需要反映平原水网地区复杂的机理变化，对分布式水文模型与实时预报的发展提供了有价值的探索。1 模型架构设计1.1 分布式架构水循环在不同的阶段和时间表现出不同的形式，其迁移和转化规律不同。即使在同一地区，不同阶段的水循环规律也不同。水循环可以分为垂直与水平循环2个过程。在垂直水循环中，水在不同介质中以不同状态进行循环迁移，主要包括2个阶段（图1）：水在空气中迁移和转化的阶段，主要由气象学家研究；陆地水循环，已得到水文学家的深入研究。在空中，水分的迁移规律相对简单，介质一致；而在陆地上，受地形、地貌和下垫面等影响，水循环可分为产流和汇流2个阶段。水汽从海洋和陆地蒸发升至天空，经大气层迁移后以降水形式分布在不同地区。在陆地上，水流遵循由高到低的原则，从山区到平原，最终流入海洋。不同地形地貌影响水流规律：山区流速快，河系单一；平原地区河道纵横交错，流向不定，湖泊和人工建筑物多，洪水运动复杂14；河口感潮区受潮汐和海水影响，水流更为复杂。这种区域性差异导致了陆地水分迁移转化规律的多样性，无法用统一理论描述。因此本文设计了分布式架构模型来解决复杂下垫面的水循环精细模拟，其具备以下特点：模型中具有吸收多种结构（集总式、分布式）、多种类型（黑箱、概念、物理基础、地形因素）已有模型的能力；模型可以根据需要采用多种类型求解方法；模型具有可以采用多种离散尺度的能力；模型参数要求能够尽量系数化，能够直接从下垫面信息提取；具有能够充分模拟全球水循环过程中垂向与横向循环的能力；模型还具有强大的可扩展能力，能够吸收其他优秀模型的能力，以及能够扩展研究区域的能力，如扩展到大气迁移转化与陆气耦合模型等。1.2 水文特征单元由于流域下垫面的复杂性，不同地形、地貌、植被等会影响产汇流规律，因此在分布式架构模型中引入水文特征单元的概念，即具有相同产汇流机理的水文地理区域。水文特征单元可分为产流型、汇流型和混合型，与传统水文单元不同，水文特征单元用于反映产汇流机制的不同。在进行水文单元划分时，需考虑最佳模型离散尺度的影响。根据水文特征单元的概念，结合分布式架构，可以构建描述流域水循环模拟的分布式架构水文模型（图2）。流域被分解为不同类型的水文特征单元，图3展示了流域主要水文特征单元组成。通过引入水文特征单元，将流域水循环分解为不同阶段、不同区域的耦合体，并针对每种水文特征单元采用最适合的模拟模型算法进行求解。在求解每个水文特征单元时，可选择最适合的模型离散时空尺度，进一步将水文特大气水分迁移高空低空地表土壤地下大气水分迁移中游平原河网区下游感潮河网区陆地水流运动海洋上游山丘区陆地水流运动图1 垂直水循环示意图图2 分布式架构水文模型结构示意图降雨、融雪冠层截留分层融雪地表产汇流壤中流运动饱和地下水运动大孔隙流上游山丘区中游平原河网区下游感潮河网区冰雪融水蒸散发损失耕地植被闸泵城镇湖泊河流图3 水文特征单元逻辑结构降雨融雪区特征单元冠层截留区特征单元蒸散发山丘区坡面特征单元喀斯特区特征单元壤中流区特征单元平原区坡面特征单元饱和地下水运动区特征单元山丘区河道特征单元城市管网特征单元平原区河道特征单元城市产汇流特征单元闸坝特征单元湖泊水库等特征单元9Mar.2024 NO.3 VOL.342024年3月 第3期 第34卷CHINA FLOOD&DROUGHT MANAGEMENT中国防汛抗旱特邀主编专栏GUEST EDITOR-IN-CHIEFS COLUMN征单元划分为水文计算单元进行模型计算。这种分布式架构模型解决了由水文特征单元组成的流域水循环模型的问题。例如，在山丘区可以采用新安江模型，将整个山丘区作为一个水文特征单元，对其进行子流域划分，各子流域即为水文计算单元；在平原区，可有河道一维特征单元、河道二维特征单元、湖泊零维特征单元、湖泊二维特征单元、堰闸特征单元和坡面特征单元等。在水文特征单元内部可进一步细分区域，这些区域即为水文计算单元，如河道一维特征单元中的断面、河道二维特征单元中的网格等。水文特征单元的分类不是唯一的，随着研究的深入和对流域水循环机理的认识不断深化，可进一步细分。各种水文特征单元的详细模拟原理详见参考文献15至18。1.3 水网地区多要素多过程耦合1.3.1 产流与汇流特征单元的耦合长江下游地区地势平坦，水系发达，河流纵横交错，多呈网状分布。其汇流过程不同于山丘区，没有统一的汇流出口。为了研究平原区产水量的分配过程，将产水汇流到周边河道里，本文提出河网多边形的概念，实现了产流过程与汇流过程的耦合。河网多边形是指由概化河网、各类分界线所围成的封闭区域（图4）。图4（a）为概化河道围成的多边形，图4（b）为概化河道与分界线围成的多边形，图4（c）为山丘区分水线与概化湖泊边界线围成的多边形。图4（a）多边形面积上的产水量或需水量必须经周边的概化河道排泄或补给；图4（b）多边形产水量或需水量只可能通过概化河道周边排出或补给，其他性质的边，如沿江、沿海、湖泊岸边、山丘分水线或平原区的分界线等，是不能排泄多边形水量的；图 4（c）是一些特殊的多边形，其周边无概化河道，如滨湖的山丘、滨湖一侧的山丘区降雨产流只可能汇入湖泊。为满足数值模拟的需要，将计算区域划分为网格（如1 km1 km）。图5是由概化河道构成的多边形，包括概化河道1、2、3、4、5，多边形面积为A，多边形的面积及各类下垫面面积，可以用其覆盖的网格计算得到。多边形所包含面积上的数值只能分配到多边形周围概化河道1、2、3、4、5。多边形中任一网格，需要明确它的产流流向哪一条河道，其灌溉需水量又是取自于哪一条河道。在没有详细地形的情况下，假定该网格与其距离最近的河道相联系，即取图中距离s最小的概化河道作为与该网格相联系的概化河道，如图5中黑色网格的产水量流入与其相联系的概化河道5，该网格的灌溉需水量亦只能从概化河道5引取，该网格与概化河道1、2、3、4无关。仅把距离远近作为网格与周围概化河道相联系的唯一依据，有时会产生不合理的结果，往往会出现一条小河与很多网格相联系的情况。如在洪水时，有太多的雨洪汇集到该条小河，来不及排泄而形成高水位，或需要太多的灌溉水量导致河道干涸。从模型计算的稳定性方面考虑，在分配中考虑河道过水能力，取综合系数：=sAR0.67（1）式中：为综合系数；A为过水面积，m2；R为水力半径，m；s为计算网格到概化河道的最小距离，m。将网格分配到多边形周边综合系数最小的概化河道，网格所属的下垫面亦随着网格的分配而归属于相应的概化河道。1.3.2 汇流特征单元间耦合水流运动模拟可由零维、一维或二维（准三维）模拟所组成，各部分模拟必须耦合联立才能求解。对模拟区域中控制水流运动的堰、闸、泵等进行耦合，控制建筑物的过流产流量可以用水力学方法来模拟，对于不同的控制建筑物采用相应的水力学公式，采用局部线性化离散出产流量与上下游水位的线性关系。分界线湖泊边界线分水线概化河道山丘边界线（b）概化河道与分界线（c）分界线（a）概化河道概化河道图4 河网多边形图5 河网多边形栅格化处理10Mar.2024 NO.3 VOL.342024年3月 第3期 第34卷中国防汛抗旱CHINA FLOOD&DROUGHT MANAGEMENT特邀主编专栏GUEST EDITOR-IN-CHIEFS COLUMN在河网水流模拟计算中河道的交汇点称为节点，对于河网节点，按节点处的蓄水面积可分为2类：节点处有较大的蓄水面积，节点水位变化导致的蓄水量变化不可忽略，这类节点称之为有调蓄节点；节点处的水量调蓄面积较小，水位变化导致的节点蓄水量变化可以忽略不计，这类节点称之为无调蓄节点。在河网一维水流计算中，节点实际上有1个基本假定：与河节点相通的河道断面水位相等，即节点水面是平的，不考虑其中的水头差。河网节点水量平衡式为：Q=A(z)Zt（2）式中：A(z)为节点调蓄面积，m2；Q为包括降雨产汇流、河道出入流在内的所有出入节点的产流量，m3/s。在河网模拟计算中，节点水位的求解是整个河网水流模型计算的关键，节点水位获得后可以计算得到河道任一断面的水位产流量。节点水位也是全流域耦合的关键。对于每一个节点，根据节点类型分别建立相应的水量平衡方程，结合边界节点的边界条件，可构建关于节点水位的完备线性方程组，从而实现水动力过程之间的耦合求解。2 预报调度功能模型的规划模拟功能不同于实时决策功能，不是简单地改变边界，需要融合多源数据，如风暴潮、台风、数值天气预报、下垫面信息等数据，并通过技术手段进行实时水雨情、工情校正。2.1 数据库耦合技术数据库管理系统主要解决模型与外部数据库的接口问题，包括数值天气预报、风暴潮、水文测站站网、水利工程工情信息的实时引入，同时负责处理模型计算成果的入库。通过设计通用数据库接口，可解决系统与数据库的接口问题，系统所有与数据相关的操作（读、写数据库）均由此接口负责完成。同时该接口在设计时考虑到数据库的开放性、数据库库表结构的可扩充性要求，即当数据库表结构发生变化时，系统只需对数据接口部分做相应调整，系统就能正常运行，而无需修改系统源代码。模型中采用ODBC标准接口与数据库相连接，可以连接任意网络数据库，如 Ms SQL Server 2000、Oracle 8 等。ODBC定义了系统与SQL数据库之间的程序设计接口，通过设置ODBC数据源，即可实现系统用户对外部数据库的读写。2.2 实时校正技术对预报模型进行实时校正，是提高复杂水网地区预报精度的有效途径之一。卡尔曼滤波理论能够同时对多维的时变状态量进行校正，因此比较适合应用到河道洪水预报实时校正中，并且卡尔曼滤波的计算特点（如在线估计时间短、存储量小等）适合实时处理和计算机运算。本研究采用卡尔曼滤波方法作为河网地区预报模型实时校正功能开发的核心算法，实现步骤为：根据河网模型各河段有无实测站点及与实测站点所在河段间的拓扑关系，将所有河段自动分级；依据基于单一河道的滤波方程求解结果，将单一河道各断面实时校正结果简化为校正单位线，在此基础上构建河道多模式单位线生成及管理模块；通过有实测站点的水位预报误差（又称为新息），对实测站点所在河道进行校正，并根据河段分级逐级传播至整个河网。实时校正流程见图6。水量预报模型根据潮位预报模型的预报潮位结果及气象部门的预报降雨资料进行预报模拟后，从实时数据库自动读取实测站点水位、报汛水位等资料，与模型相应位置计算值对比获得新息源，按照河段分级顺序，根据各河段校正单位线，从各实测站点逐级沿河段将新息传播到全河网，对所有节点、断面水位预报水位进行校正。2.3 水利工程调度技术调度模型设计了不同触发规则，支持水量条件、水质条件、时间条件等，运行情况模拟支撑闸泵工程的分孔控制、逐步启闭、不同开启高度等。水利工程的具体控制调度方法在模型中是以控制条件的形式实现，创建联系控制条件即实现联系的具体控制调度，相同调度规则的闸泵工程图6 实时校正流程图河段拓扑关系分级卡尔曼滤波所有河道的实时校正单位线河段分级实时校正断面水位误差逐级传播河网实时校正河网水量预警预报11Mar.2024 NO.3 VOL.342024年3月 第3期 第34卷CHINA FLOOD&DROUGHT MANAGEMENT中国防汛抗旱特邀主编专栏GUEST EDITOR-IN-CHIEFS COLUMN北水体环太湖闸泵圩区河道图例太湖可设置在同一个控制条件中。联系控制有2种方法，GATE格式控制和具体控制。以GATE格式控制为例，控制条件设置步骤如下。（1）设置水位产流量统计。设置统计河道或湖泊的水位产流量状态量，明确控制代表站的断面或节点位置、间隔时段等，设定状态量名称，作为控制条件的判别依据。若判别条件不涉及水位产流量，可不建立水位产流量统计。（2）设定水位相关线和调度线。通常与代表站点调度线有关的闸泵调度需设置此项，若其他闸泵涉及代表站以外的水位相关线，也可再进行分级设置。（3）编写控制条件。闸门运行的规则调度标准，闸门依据控制条件列表中的1条或多条的任意组合调度，返回值为闸门的启闭。（4）确定控制步长。闸门开启后，对开启的方式进行具体控制，以满足实际运行中的特殊要求。（5）在GATE格式控制步长标签下，按照调度规则编写工程具体控制逻辑，并在联系要素列表中选定按此调度规则进行调度的闸泵工程，增加到具体控制信息中。3 太湖流域应用3.1 研究区域太湖流域地处我国东部长江河口段南侧及钱塘江与杭州湾之间，三面滨江临海，一面环山，其北滨长江，南濒钱塘江，东临东海，西以天目山、茅山等山区为界。太湖流域行政区划分属江苏省、浙江省、上海市和安徽省19。流域总面积36 895 km2，其中江苏省19 399 km2，占比为52.6；浙江省12 095 km2，占比为32.8；上海市5 176 km2，占比为14.0；安徽省225 km2，占比为0.6。3.2 模型概化根据上述分布式架构水文模型提出的水文特征单元概念，结合太湖流域的具体情况，太湖流域可以概化为由下面几类特征单元组合构成（表1）。根据太湖流域平原水网的特点，将流域内影响水流运动的因素分别概化为零维模型（湖、荡、圩等零维调蓄节点）、一维模型（一维河道）和联系要素3类。联系要素包括水闸、船闸、涵洞（立交）及泵站等工程，联系要素的上下游都设有节点，节点之间的水位差与产流量之间的关系取决于堰流公式选取及运行方式。太湖流域河网、湖泊概化情况见图7。共包括概化河道（段）1 483条；计算断面4 274个；节点116个，其中边界节点105个，调蓄节点76个；控制建筑物（闸、泵等）168个。3.3 模型率定选取2016年实况数据对模型参数进行率定，率定水位站选取太湖和均匀分布于湖西区、武澄锡虞区、阳澄淀泖区域内的代表站点，率定结果见图 8和表 2。率定结果显示，计算与实测水位过程拟合程度非常高，太湖计算日均最高水位与实测水位基本一致。地区代表站日均最高水位误差基本在0.1 m以内。研究范围内所有站点Nash系数fNSE均值为0.97，确定性系数R2均值为0.96，其中太湖fNSE均值为 0.97，R2均值为 0.99，湖西区 fNSE均值为 0.99，R2均值为0.99，武澄锡虞区fNSE均值为0.95，R2均值为0.98，阳澄淀泖区fNSE均值为0.97，R2均值为0.91。各站水位预报合格率均高于85%，预报项目精度达到甲级等级。其中6个代表站的预报精度在90%以上，太湖水位预报精度最高，达98.9%。特征单元类型山丘区子流域及坡面单元平原区坡面单元平原区河道单元平原区湖泊、水库、行蓄洪区单元闸坝工程单元模型算法新安江模型平原区降雨径流模型一维水动力学模型零维模型准三维模型闸坝工程模拟算法概化详情浙西山丘区，共分为10个子流域湖西山丘区，共分为10个子区域主要为平原区部分，共分为16个子区域共概化了1 483条内河道，总河长7 879.4 km，断面数4 274个共概化了104个湖泊，总面积为871 km2太湖湖区共概化了168个大闸表1 太湖流域主要概化基本信息图7 太湖流域概化12Mar.2024 NO.3 VOL.342024年3月 第3期 第34卷中国防汛抗旱CHINA FLOOD&DROUGHT MANAGEMENT特邀主编专栏GUEST EDITOR-IN-CHIEFS COLUMN5.53.0（c）无锡（大）水位水位/m日期（月-日）4.04.55.03.506-0111-0102-0109-0105-0101-0108-0103-0104-0107-0110-0112-01计算值实测值12-31图8 2016年太湖水位及部分代表站率定成果表（d）枫桥水位水位/m日期（月-日）06-0111-0102-0109-0105-0101-0108-0103-0104-0107-0110-0112-01计算值实测值12-315.53.0（a）太湖水位水位/m日期（月-日）6.04.04.55.03.506-0111-0102-0109-0105-0101-0108-0103-0104-0107-0110-0112-01401801401201002008060计算值太湖流域降雨量实测值5.53.0（b）王母观水位水位/m日期（月-日）6.04.04.55.03.506-0111-0102-0109-0105-0101-0108-0103-0104-0107-0110-0112-01计算值实测值6.57.07.512-3112-31平均降雨量/mm水利分区骨干河湖湖西区武澄锡虞区阳澄淀泖区站点名称太湖平望王母观坊前常州（三）无锡（大）青阳枫桥湘城陈墓最高日均水位/m计算值4.864.106.605.906.265.195.224.714.003.82实测值4.864.106.535.796.175.175.204.673.983.82差值0.000.000.070.110.090.020.020.040.020.00fNSE0.970.980.980.990.960.950.940.970.960.99R20.990.980.990.990.990.980.970.980.800.96合格率（%）98.994.589.988.586.686.993.492.997.094.3表2 2016年水位率定成果表3.4 实时预报为评估实时校正在实时预报中的效果，同步对比设置校正点和不设置校正点的预报成果，在太湖流域模型中对太湖代表站水位设置实时校正点、计算时段选取2016年58月，预报期降雨、边界条件选用实测资料，降低外部输入误差，预报期工程调度选用规则调度，以计算时段内每日8：00作为预报依据时间，预报未来3 d太湖及各区域代表站水位，分别选取 1 d、2 d、3 d预报成果与实测水位对比，计算确定性系数、最大误差、洪峰峰值误差等。如2016年 5 月1日8:00作为预报依据时间，预报2016年5月2日8:00、2016年5月3日8:00、2016年5月4日8:00这3个时刻的水位成果，分别作为1 d、2 d、3 d预见期成果，与对应时刻实测水位做对比。对不同预报期的预报成果使用确定性系数R2和合格率QR来进行评估，其中合格率QR的计算公式为：QR=n/m100%（3）式中：m为预报总次数；n为合格预报次数。对合格预报次数的判定标准为：实测值减去计算值的绝对值小于0.1。太湖作为太湖流域主要调蓄水体，其水位预报精度对太湖流域防洪等至关重要，太湖水位不同预见期成果误差分析如表3所示。当不设置实时校正点时，太湖水位1 d、2 d、3 d预见期预报成果与对应时刻太湖实测水位相比，确定性系数只能达到0.607、0.570、0.524，当设置校正点之后，确定性系数提高到0.996、0.988、0.974；最大误差由 0.52 m、0.56 m、0.59 m 降低到 0.08 m、0.12 m、0.17 m；峰值误差由0.08 m、0.09 m、0.10 m降低到0.03 m、0.06 m、0.04 m，预报合3.04.04.55.03.513Mar.2024 NO.3 VOL.342024年3月 第3期 第34卷CHINA FLOOD&DROUGHT MANAGEMENT中国防汛抗旱特邀主编专栏GUEST EDITOR-IN-CHIEFS COLUMN格率由35.5%、35.5%、35.5%升至99.2%、95.2%、80.6%。从太湖水位1 d、2 d、3 d预见期预报成果与实测数据对比图（图9）可以看出，不管是水位变化趋势和峰值水位，在设置校正点之后，预报精度都有明显改善。在没有设置实时校正点的情况下，太湖水位随着预报时间的推进，预报误差逐步增大，而从设置了实时校正点的预报成果来看，预报误差得到了有效控制。同时在设置校正点和不设置校正点2种计算工况下，预报期条件一致，因此误差来源主要是预报依据时间所在时刻模型状态和流域实际状态的不同，由此可见，该实时校正方法可以大幅提升太湖水位预报精度。预报方案1 d预见期、校正1 d预见期、无校正2 d预见期、校正2 d预见期、无校正3 d预见期、校正3 d预见期、无校正确定性系数0.9960.6070.9880.5700.9740.524最大误差/m0.080.520.120.560.170.59峰值误差/m-0.030.08-0.060.09-0.040.10QR（%）99.235.595.235.580.635.5表3 太湖水位不同预见期预报成果误差分析表50（c）预见期3 d水位/m日期（月-日）234106-2008-0905-1107-2006-1005-0107-1005-2105-3106-3007-3008-19太湖3 d预见期校正太湖实测太湖3 d预见期无校正08-2950（b）预见期2 d水位/m日期（月-日）234106-2008-0905-1107-2006-1005-0107-1005-2105-3106-3007-3008-19太湖2 d预见期校正太湖实测太湖2 d预见期无校正08-2950（a）预见期1 d水位/m日期（月-日）234106-2008-0905-1107-2006-1005-0107-1005-2105-3106-3007-3008-19太湖1 d预见期校正太湖实测太湖1 d预见期无校正08-29图9 太湖水位1 d、2 d、3 d预见期预报成果与实测数据对比4 结论综合考虑平原河网区的汇流特征和汇流影响因素，本文提出了基于水文特征单元的分布式架构模型。以太湖流域为原型，构建了一个完整的多要素、多尺度、多过程的水循环精细模拟模型。（1）提出了分布式架构水文模型，能够反映某一区域机理变化的耦合模型。该模型模拟的对象是考虑人类活动影响的流域水循环系统，结构具有分布式特点，考虑了降雨和下垫面条件空间分布不均匀对流域降雨径流的影响。（2）引入水文特征单元，使得模型在计算时可以根据流域特征和实际需要采用最合适的模块组合，如集总式、概念性或分布式模型。（3）实时预报调度功能，可耦合数值天气预报、风暴潮、水文测站站网、水利工程工情信息，并进行数据同化。这些改进提高了模型的精度和适用性，为平原水网地区流域水资源管理和防洪减灾提供了有力支撑。参考文献1 钱真，孔令婷，陈长太，等.城市平原感潮河网地区涉河工程防洪影响计算方法研究J.水利规划与设计，2022（9）：41-45.2 陈娟，毛雨，刘一帆，等.张建云：2020年长江洪涝灾害应对工程与非工程措施J.中国防汛抗旱，2021，31（11）：35-36.3 周焕，闪丽洁，揭梦璇，等.基于混合回归模型的沿海平原河网洪水位预报研究J.中国水利水电科学研究院学报（中英文），2022，20（2）：145-152.4 杨松彬，董志勇.河网概化密度对平原河网水动力模型的影响研究J.浙江工业大学学报，2007（5）：567-570.5 梁犁丽，冶运涛，龚家国，等.分布式水文模型在短期水文预报中应用的可行性探讨J.中国水利水电科学研究院学报，2013，11（3）：210-215.6 李菁，王船海，张萍霞，等.基于节点传播的河网水动力学模型的实时校正研究J.水动力学研究与进展A辑，2015，30（1）：75-82.7 曹大岭，李建州，喻海军.山丘区城市实时洪涝分析预警系统研发及应用J.中国防汛抗旱，2022，32（1）：84-89.8 Simmons T C，Brunner P，Therrien R，et al.Commemorating the 50thanniversary of the Freeze and Harlan（1969）Blueprint for a physically-based，digitally-simulated hydrologic response modelJ.Journal of Hydrology，2020，584.9 Salim G，David M，Joseph H.A Generalized Multistep Dynamic（GMD）14Mar.2024 NO.3 VOL.342024年3月 第3期 第34卷中国防汛抗旱CHINA FLOOD&DROUGHT MANAGEMENT特邀主编专栏GUEST EDITOR-IN-CHIEFS COLUMNtopmodelJ.Water Resources Research，2023，59（1）.10Koudelak P，West S.Sewerage network modelling in Latvia，use of InfoWorks CS and Storm Water Management Model 5 in Liepaja cityJ.Water and Environment Journal，2007，22（2）：81-87.11Kumar P，Lohani A K，Nema A K.Rainfall runoff modeling using MIKE11 nam modelJ.Current World Environment，2019，14（1）：27-36.12喻海军，马建明，张大伟，等.IFMS Urban软件在城市洪水风险图编制中的应用J.中国防汛抗旱，2018，28（7）：13-17.13王船海，郑世威，李小宁，等.基于太湖流域模型的城市内涝过程高效模拟J.水科学进展，2022，33（3）：462-473.14关铜垒，刘佳嘉，周祖昊，等.考虑山区平原区的分布式水文模型计算单元划分方法研究J.水文，2022，42（5）：7-13.15WANG C H，HUA W J，CHEN G，et al.Distributed-framework basin modeling system：.Overview and model couplingJ.Water，2021，13（5）：678.16CHEN G，HUA W J，FANG X，et al.Distributed-framework basin modeling system：.Hydrologic modeling systemJ.Water，2021，13（5）：744.17LI X N，WANG C H，CHEN G，et al.Distributed-framework basin modeling system：.Hydraulic modeling systemJ.Water，2021，13（5）：649.18CHEN G，WANG C H，FANG X，et al.Distributed-framework basin modeling System：.Application in Taihu BasinJ.Water，2021，13（5）：611.19程文辉，王船海，朱琰.太湖流域模型M.南京：河海大学出版社，2006.Development of flood forecasting models in river network areas and theapplication in the lower reaches of the Yangtze RiverCHEN Gang1,ZHAO Yanhong2,WANG Xiaoshu2,CHEN Qiushi1,WANG Chuanhai1(1.College of Hydrology and Water Resources,Hohai University,Nanjing 210098;2.Yishusi River Basin Administration Hydrological Bureau(Information Center),Changzhou 221018)Abstract：Flood forecasting in plain river networks has always been a challenging issue in hydrology,especially with thesignificant impact of human activities.In order to better simulate the hydrological cycle in the basin and address the challengesin engineering applications,this paper proposes a novel distributed architecture hydrological model.The model has distributedcharacteristics,allowing for flexible selection of the most suitable hydrological feature units based on basin characteristics andrequirements,thereby supporting the precise and theoretical study of different hydrological feature units.To achieve real-timeforecasting functionality,this paper also developed database coupling technology,large-scale river network real-time correctiontechnology,and hydraulic engineering scheduling technology.Based on this,a distributed architecture model of the Taihu Lakebasin was established,which was generalized,and the model parameters were calibrated using 2016 actual data.By conductingreal-time forecasts of Taihu Lake water levels with lead times of 1 d,2 d,and 3 d,the results of this study can basically reflectthe actual water flow conditions in the Taihu Lake basin of the plain river network area.In summary,the distributed architecturehydrological model method proposed in this paper has not only been verified in terms of scientificity but also has feasibility inimplementation technology.It successfully solves the application problems from concept to reality and provides new ideas andmethods for flood forecasting in plain river network areas.Keywords：river networ