1. 引言
平原区暴雨洪水特征与所在地理区位密切相关，因水利工程调度运行，导致洪水受自然流域、外调洪水等多种因素控制，是我国多数平原区洪水形成的原因。在暴雨洪水特性研究方面，国内外学者常对典型暴雨洪水开展深入分析研究，洪水特性研究是当前研究热点，如2018年新沐河特大洪水 [1] 、沂沭泗流域2019年台风洪水 [2] 、淮河流域颍河2021洪水 [3] 、桃江上游2022年洪水 [4] 、四川泸州市2014年暴雨洪水 [5] 等。平原区暴雨洪水特点及调控历来是国内外学者研究重点，V. S. LEBEDEVA、Godwin E.等分别在伏尔加–阿赫图、乌干达东部低地平原地区开展了洪水的水工程时空调度及风险评估研究 [6] [7] 。我国平原区承载了大部分人口及经济产值，其暴雨洪水特征已有较长的研究历史，如20世纪50年代龚召熊以5次大洪水为实例对长江中游平原区洪水特性进行了总结分析 [8] ，近年来，杨劲松等对平原区暴雨洪水沉积重建研究进行了归纳分析 [9] ，刘洋采用灾害场次频率分析法和实测水面线反推法等探索了平原区洪涝灾害调查评价方法 [10] 。然而，对平原区因工程调度引起的多源洪水特征及应对策略有待进一步深入研究。
另外，近年大数据技术发展迅速，水利科研人员在洪涝灾害预警、洪水预报等方面积极探索大数据技术应用 [11] [12] 。应用粗糙集(Rough Sets)、决策树(Decision Tree)、神经网络(Neural Network)、相关规则、K-nearest邻居、联机分析处理、数据可视化等方法 [13] ，对数据进行总结、分类、关联及聚类是数据挖掘有效途径，其中相关规则法描述了一个事物中某些属性同时出现的规律和模式，为本研究的基本方法。
位于沂沭泗流域的新沭河下段平原区是多水源洪水典型代表区域。新沭河是在原沙河基础上开挖的漫滩行洪入海人工河道，西起大官庄枢纽，向东流经大兴镇入江苏省石梁河水库，出库后东流至临洪河口入海。新沭河河道全长80 km，其中山东境内长20 km，入江苏境内至石梁河水库坝址长15 km，石梁河水库泄洪闸至入海口长45 km，本文称该段为新沭河下游，区间流域面积2850 km2。入海口三洋港闸2008年10月动工兴建，2013年12月建成投运。
据2021年1月鉴定评估及相关资料分析 [14] ，三洋港闸下游回淤严重，影响工程安全，过流能力安全评级为B级，安全评定为二类闸。针对新沭河下游河段凸显的问题，利用洪水资源进行水闸工程调度是防淤减淤重要的非工程措施。开展新沭河暴雨洪水特征分析有利于掌握洪水资源利用途径与方式，为工程调度取得先机条件。三洋港闸调度控制的蓄泄过程、上游降雨等均属于随机事件的非平稳时间序列，开展数据集制定与处理方法研究，其成果可为多水源平原区开展暴雨洪水关系大数据方法研究奠定良好基础。
2. 分析方法
2.1. 水源分区
根据新沭河下游区域地理区位特点，以地形特征及水利工程调度为分区依据，将新沭河下游洪水来源分为沂河、沭河、石梁河水库、蔷薇河、石安河、范河、连云港城区等7个区域，图1为分区概化示意图。
图1. 暴雨分析区域与新沭河下游关系示意图
2.2. 降雨特征属性
本文以降雨历时、总降雨量、站点24 h最大降雨强度、覆盖面积、最强点距新沭河距离等特征反映区域暴雨状况，分别定义如下。
1) 降雨历时：以区域内各代表站进行分析，区域内本场降雨起始时间至结束时间。
2) 总降雨量：降雨历时内，区域内各站降雨量代数和的平均值，即平均降雨量。
3) 站点24 h最大降雨强度：降雨历时内，区域内各代表站中单站24 h最大降雨量。
4) 覆盖面积：若区域内各代表站均有降雨，则该覆盖面积为区域总面积，否则为降雨测站占比数对应的面积。
5) 暴雨中心距新沭河距离：最大降雨量站点与新沭河下游临洪断面的距离，或称为中心与新沭河距离，如图1所示的距离概化线。
本文对7个区域各场次暴雨分别计算出以上特征，以新沭河下游洪水最大水位及出现时间作为洪水特征，用于分析各区域暴雨对新沭河下游洪水影响。
2.3. 分析方法
洪水影响采用相关性分析方法，本研究采用皮尔逊相关系数法，计算总降雨量、最强点距新沭河距离的洪水相关系数，并进行对比分析，得出各区域对新沭河下游洪水影响程度。该方法步骤如下：
1) 整理暴雨洪水资料，得出场次洪水的暴雨分布特征；
2) 计算各区域暴雨的总雨量(P)、暴雨中心距新沭河距离(D)，计算各场次洪水新沭河下游最大水位(Z)、出现时间与最大雨量发生时间差值(T)；
3) 建立(P,Z)、(P,T)、(D,Z)、(D,T)等散点系列；
4) 数据合理性分析与整理；
5) 计算各区域(P,Z)、(P,T)、(D,Z)、(D,T)相关系数；
6) 分析各区域的相关系数，确定各区域降雨与新沭河下游洪水关系。
3. 数据资料
本研究以2013年至2022年期间实测降雨、洪涝水情等实时雨水情数据 [15] 及水文年鉴为资料来源，洪涝水情有临洪枢纽闸坝水情及三洋港闸实测闸上水位、闸下水位、过闸流量等。
根据2013至2022年的雨洪资料，三洋港闸出现行洪过程年份有2018年、2019年、2020年、2021年、2022年，其中2018年新沭河出现2次、2020年3次洪水过程，其他年份各1次，依次编号为2018-1、2018-2、2019-1、2020-1、2020-2、2020-3、2021-1、2022-1，研究区各洪水过程降雨时段及24 h最大雨量如表1所示，各场次洪水暴雨特征如图2所示。
表1. 研究区各洪水过程降雨时段及24 h最大雨量
图2. 研究区场次洪水暴雨特征值
三洋港闸过闸流量观测资料不足，本研究以水位为主要要素反映新沭河下游洪水水情，代表站为三洋港闸站，采用其上游水位。以各场次洪水三洋港闸上水位及最大值发生时间反映洪水特征。各场洪水闸上水位特征值分别如图3所示。
图3. 各场次洪水三洋港闸上洪水特征信息图
4. 数据处理分析
4.1. 相关性计算
以各区域暴雨特征、新沭河下游洪水特征等为数据源，构建各区域暴雨与新沭河下游洪水特征值散点系列，即(P,Z)、(P,T)、(D,Z)、(D,T)。分别计算相关系数。计算成果如图3所示。从相关系数计算成果分析可知，沂河区、沭河区的总雨量与新沭河下游水位具有较好的正相关，其他区域各散点相关性较差。
图4. 各区域暴雨洪水散点相关系数分布图
4.2. 合理性分析
从表1、图2各暴雨特征可看出，沂沭区各暴雨总量及覆盖面积均大于其他区域，占明显优势，对新沭河下游洪水具有控制性作用，其相关性计算成果符合实际。
5. 主要结论
新沭河下游地处沂沭泗下游平原区，洪水来源可划分为7个区域，属典型多水源平原区，经对2013年三洋港闸建成运行以来的近10年的暴雨洪水关系研究，形成以下主要结论：
1) 新沭河下游，需接纳蔷薇河区域、石安河区域、范河区域、城区等周边洪水，同时是沂沭河区域、石梁河水库区域等上游由水利工程调度泄洪通道，洪水水源组成复杂。
2) 沂沭区暴雨量与三洋港闸上游水位具有较好的正相关，受水利工程调度分洪驱动，对于新沭河下游，沂河区、沭河区的外调洪水具有控制性作用，三洋港闸调度运行应密切关注沂沭区暴雨状况。
3) 本文给出了多水源洪水与暴雨之间的特征值相关性分析方法，可作为洪水与暴雨数据挖掘模拟计算研究的数据处理基础。
基金项目
江苏省水利科技项目(2022020)资助。