摘要：以海口市主城区为例，基于PCSWMM软件构建一维、二维耦合内涝模型，根据多情景模拟结果，分析内涝灾害特征及分类成因，并针对性地提出组合减灾方案。结果表明，基于PCSWMM的内涝模型可提供支撑成因分析和防灾决策的可靠信息。海口市主城区积水成因主要包括降雨量超标、潮位顶托、河道和管道排水能力不足以及地势低洼和地形平坦。通过管道、河道改造结合增设蓄水池和排涝泵站的组合减灾方案可以有效缓解区域内涝。
基于高精度的数值模拟模型进行城市内涝过程模拟和成因分析，是当前的主流思路。笔者基于PCSWMM软件构建海口市主城区一维、二维耦合内涝模拟模型，获得降雨和潮位组合情景下的主城区积水程度及内涝分布特征；在此基础上，揭示内涝形成机制和关键影响因素，并从整体层面提出可行的整治方案，以期为海口市主城区的内涝防治提供科学可靠的决策支持。
01 研究区域概况
海口市地处低纬度热带北缘，属于热带海洋气候，全市多年平均降雨量为1827mm，其中5月—10月为雨季，降雨量占全年降雨量的78.1%；9月为降雨高峰期，平均降雨量为300.7mm，占全年的16.5%。海口市主城区位于海口市北部，濒临琼州海峡，是海南省政治、经济、科技、文化中心和最大的交通枢纽。受强降雨与风暴潮影响，海口市内涝呈突发、频发和高危害的特征。
02 海口市主城区内涝模型构建
根据管网、检查井与河道信息（来源于海口市水务局），以及下垫面利用信息和地形高程数据（遥感解析）等，构建海口市主城区一维、二维耦合内涝模拟模型。所构建模型共包含子汇水区139个、地表网格62920个、管段13122根、检查井13073个、河道10条、排水口47处，其中在龙昆沟出海口等35处设置挡潮闸。挡潮闸运行方式：当闸上游水位高于闸下游潮水位时，开闸排水，反之关闸挡潮。
采用2014年威马逊台风期间海口市的降雨和潮位数据作为输入边界条件，以区域内各内涝点实测深度与模拟深度作为评判对象，对模型参数进行率定。实测与模拟结果的对比情况见图1。实测积水深度范围为20～80cm，模拟积水深度范围为22～75cm，平均值分别为47.3、46.8cm，仅相差0.5cm。线性回归决定系数（R2）、纳什效率系数（NSE）以及相对误差（Re）分别为0.80、0.796、0.9%，一般认为当R2>0.6、NSE>0.5、|Re|<20%时，模型的拟合精度即满足要求。综上可知，构建的内涝模型能较好地模拟海口市主城区的积水与排涝过程。
03 模拟研究及结果分析
利用所构建模型得出不同设计降雨耦合同频率潮位情景下海口市主城区的内涝情况如下：当重现期分别为2、5、10、20、50 年时，积水面积分别为11.13、14.50、16.67、18.29、20.94km2，积水体积分别为227.90×104、348.68×104、419.68×104、491.22×104、620.21×104m3，平均积水深度分别为0.20、0.24、0.25、0.27、0.30m。50年一遇情景下的最大积水深度及其分布如图2所示。
模拟结果表明，在2年一遇降雨情景下，研究区域的积水面积占总面积的12.7%；在50年一遇情景下，积水面积、积水体积、平均积水深度较2年一遇情景下分别增加了88%、172%、50%，积水体积的增加幅度明显大于积水面积，表明积水产生后难以排出，内涝风险随重现期的增加而显著增加。海甸岛及龙昆沟排水子系统内滨海大道、道客沟-龙昆沟、大同沟-大同分洪沟沿线区域为严重内涝区。海甸岛与龙昆沟排水子系统是海口市主城区内涝灾害高风险区域。
04 内涝成因分析
以内涝程度最严重的龙昆沟排水子系统为例，开展内涝成因分析。在50年一遇情景下，龙昆沟排水子系统的内涝分布如图3所示。
4.1 区域积水与区域参数关联性
采取Pearson相关性分析法，以子汇水区为单元，研究龙昆沟区域积水特征与区域参数之间的关联性。区域参数共5个，其中积水面积占比是指最大积水时刻子汇水区内积水网格面积与子汇水区面积之比；积水体积是指最大积水时刻子汇水区内的积水量；平均积水深度是指最大积水时刻子汇水区内积水量与子汇水区面积之比；综合积水程度是通过熵权法对积水面积占比、积水体积及积水时长求取权重后的归一化加权值，其值越小，对应的积水程度越严重。关联性结果如表1所示。可知，坡度与积水面积占比、综合积水程度分别呈显著相关、显著弱相关关系；平均高程与积水面积占比、平均积水深度、综合积水程度分别呈显著强相关、显著弱相关、显著相关关系；积水情况与管道长度的相关关系较弱。因此，区域积水程度受平均高程与坡度的影响较大。
4.2 典型积水区域内涝成因分析
选取典型积水区域，按照雨水传输路径分析积水区域的内涝成因。
4.2.1 降雨量超标和潮位顶托
相同潮位和不同降雨重现期组合情景下龙昆沟下游区域的内涝情况如表2所示。50年一遇重现期时下游积水面积与积水量较2年一遇分别增加了336%、654%。因此，降雨量超标是导致严重内涝的重要原因。
提取潮位顶托下龙昆沟、大同分洪沟出口排涝过程与大同沟水位，如图4所示。可知，在08：45—11：30期间，受潮位顶托作用，河道排涝流量减小，出现间断排涝或者无法排涝情况，导致水位快速上升，大同沟水位超过堤防高程（2.80m），造成漫溢，急剧降低了河道对后续来水的调蓄能力。
为研究潮峰与雨峰相对位置对排涝过程的影响，对比潮峰与雨峰时刻重合以及超前或滞后2h和4h共5种情景下龙昆沟排水子系统的总积水量，如图5所示。可知，潮峰位置对排涝过程存在显著影响，潮峰滞后雨峰时，将导致汇集涝水高峰与高潮位重叠，涝水无法自流排出，加剧内涝。最大影响对应的潮峰滞后时间与区域降雨产汇流过程和排水系统的排水特征均相关。
4.2.2 管道排水能力不足
滨河路某管道流量过程及管道周边地表积水过程、管道与大同沟相交断面水位过程见图6（管道流量为负代表河道涝水倒灌进入管道）。可知，高潮位与强降雨时段，河道水位较高，管道流量较小甚至变为负值，而地表积水深度与河道水位变化趋势一致。因此，当管道排水受河道顶托与倒灌影响时，区域内涝程度显著增加。
龙昆南路沿线涝水漫溢严重，沿程最大流量如图7所示。此段下游排水管道尺寸显著减小且存在逆坡，使得涝水无法有效排出，中间段涝水汇入后反方向流向上游，产生沿程漫溢。
4.2.3 河道排水能力不足
大同沟周边区域为内涝高风险区，其上游和下游断面流量过程如图8所示。可知，在降雨高峰时段，受河道排水能力制约，大同沟入口和出口断面流量峰值相差26.9m3/s，导致涝水沿程漫溢至地表。道客沟下游段也存在相同问题。
4.2.4 地形地势
地形地势决定涝水扩散路径，对内涝演变过程和分布格局具有重要影响。以下游区域子汇水区O2为例，其地形与内涝过程如图9所示。可知，O2子汇水区地势中间高、四周低。西南方向子汇水区地表积水到达图中A点后顺着高程低的龙昆北路与龙华路流入O2子汇水区，汇入量达62.06×104m3（未扣除O2子汇水区进入其他排水分区的涝水量），远超管道漫溢水量。因此，地形平坦、地势过低将加剧区域涝水汇入，加重内涝程度。
除O2子汇水区外，红城湖西侧子汇水区D13_2地势明显低于南部与西部相邻区域，南部与西部涝水向地势较低的北部和东部区域流动，导致大面积、长时间积水（见图10），总积水量达16.6×104m3，中部地势最低点积水深度达到1.7m。因此，区域如存在地势较低情况，将极易形成严重内涝。
05 区域内涝治理方案
依据区域内涝成因，提出龙昆沟排水子系统整治方案，包括管道和河道改造及增设蓄水池和雨水泵站，具体见图11。治理前后积水情况对比见图12。
① 管网系统改造
具体措施：在滨河路、义龙东路、大同路等二级道路，改造管道76根、新建管道36根、设置管道拍门15个；在滨河路、华海路周边各增设1座雨水泵站；在红城湖管道沿线修建3处排水出口。
治理效果：下游管道改造所在子汇水区的积水面积与积水体积削减率分别为19.6%、33.0%；上游红城湖管道沿线排水量较改造前增加111×104m3。
② 河道改造
具体措施：参考海口市蓝线规划，将大同分洪沟暗涵由8m拓宽为20m；大同沟部分断面堤防高程提高至整体堤防高程即2.80m；道客沟下游段暗涵拓宽至10m。
治理效果：大同分洪沟最大排涝流量增加了33m3/s；道客沟排入龙昆沟的水量增加了20m3，上游D13_1子汇水区积水体积减少了30%。
③ 增设蓄水池
具体措施：在滨海新村附近、南大立交桥西南侧、红城湖与龙昆南路之间分别布设蓄水池。
治理效果：蓄水池周边15个子汇水区的积水面积和积水体积分别减少了38.2%、62.2%；红城湖周边子汇水区的汇入水量减少了24.2%。
④ 增设排涝泵站
具体措施：在东西湖周边、上游蓄水池3与道客沟之间分别布设排涝泵站，综合考虑积水程度及成本，上、下游泵站规模分别设置为10、55m3/s。
治理效果：大同沟沿线及东西湖周边7个子汇水区的积水面积和积水体积分别减少了47.6%、59.3%，蓄水池1和蓄水池3的规模分别降低了85.0%、29.1%。
模拟分析结果验证了所提方案的有效性。方案如实施，龙昆沟排水子系统下游面状积水程度将显著降低，红城湖周边、道客沟-龙昆沟、大同沟-大同分洪沟沿线等内涝严重区域的积水深度及积水面积将明显减少。据统计，积水面积和积水体积将分别减少30%和37%，内涝缓解效果显著。
06 结论
① PCSWMM软件满足开展城市内涝精细化模拟分析研究的需求，在完善基础资料的条件下，基于PCSWMM软件的内涝模拟模型可提供支撑成因分析和防灾决策的可靠信息。
② 受多方面因素影响，海口市主城区易发生内涝，且内涝程度随降雨和潮位的提高而显著增加。以下游区域为例，50年一遇情景的积水面积与积水量较2年一遇分别增加了336%、654%。应结合城市定位和发展，确定合理的内涝防治标准。
③ 海口市主城区的内涝成因包括：降雨量超标和潮位顶托、河道和管道排水能力不足、地势低洼及地形平坦等。内涝高风险区域分布在海甸岛及龙昆沟排水子系统内的滨海大道、道客沟-龙昆沟、大同沟-大同分洪沟沿线。区域间因河渠和道路连接而紧密联系，应从整体考虑制定防灾策略和方案才能达到预期效果。
④ 龙昆沟排水子系统分析结果表明，管道和河道改造结合增设排涝泵站和蓄水池的组合减灾方案是沿海城市防灾的优选，但应合理确定各类措施的布设位置和规模，同时还需考虑区域间涝水联系来制定精细化调控方案，以达到最佳防灾效果。
