(19)中华 人民共和国 国家知识产权局 (12)发明 专利申请 (10)申请公布号 (43)申请公布日 (21)申请 号 202111129808.6 (22)申请日 2021.09.26 (71)申请人 中水北方勘测设计 研究有限责任公 司 地址 300222 天津市河西区洞庭路6 0号 (72)发明人 张育德　回晓莹　冯宇鹏　陆建宇　 王杰　韩璞璞　郑永路　石宝　 丛娜　程增辉　 (74)专利代理 机构 天津市鼎和专利商标代理有 限公司 12101 代理人 李凤 (51)Int.Cl. G06F 30/27(2020.01) G06N 3/00(2006.01) G06N 3/12(2006.01)G06F 113/08(2020.01) (54)发明名称 一种无资料地区的逐月径流系列推求方法 (57)摘要 本发明公开了一种无资料地区的逐月径流 系列推求方法， 涉及水文学径流设计领域。 通过 选取参证站流域， 整理参证站流域的降雨径流系 列， 建立参证站流域的产流模型， 拟合参证站流 域的产流模 型参数， 使用参证站流域拟合的产流 模型参数以及设计流域的逐月降水量， 计算设计 流域的逐月径流系列。 本方法资料需求较少， 可 以在仅有逐月降雨资料地区推求逐月径流系 列， 对于资料缺乏地区有一定实际意 义。 权利要求书1页 说明书6页 附图2页 CN 113919211 A 2022.01.11 CN 113919211 A 1.一种无资料地区的逐月径流系列推求方法， 其特 征在于， 包括 步骤如下： S1， 选取与设计流 域产流特点相似的参证站流 域； S2， 将参证站流域降水资料整理为逐月降水系列， 将参证站径流资料整理为逐月径流 系列； S3， 使用S2中整理的逐月降水系列和逐月径流系列， 拟合模型参数， 确定参证站流域降 雨径流模型； S4， 整理设计流 域的降雨资料， 统计为逐月降水系列； S5， 采用S3中确定的 降雨径流模型、 拟合的模型参数以及S4中的逐月降水系列， 推求设 计流域的逐月径流系列。 2.根据权利要求1所述无资料地区的逐月径流系列推求方法， 其特征在于， S1中选取的 参证站流 域与设计流 域在下垫面条件、 产汇流类型、 气候特点特 征上存在相似性。 3.根据权利要求1所述无资料地区的逐月径流系列推求方法， 其特征在于， S2中降水逐 月系列为流域内的面平均降水系列， 根据流域内雨量站使用算数平均法或泰森多边形法计 算， 降雨、 径流资料的同步系列长度在10年以上。 4.根据权利要求1所述无资料地区的逐月径流系列推求方法， 其特征在于， S3参证站流 域的降水径流模型为： W蓄i＝[W蓄i‑1+(Pi参‑P损)*α月]*(1‑k)        (1) Ri月＝[W蓄i‑1+(Pi参‑P损)*α月]*k          (2) 式中： W蓄i、 W蓄i‑1为第i月、 第i ‑1月的月末流域存蓄水量； Pi参为参证站流域第i月的降水 量； P损为月损失水量； α月为月径流系数； k 为月退水系数； Ri月为月径流 量。 5.根据权利要求4所述无资料地区的逐月径流系列推求方法， 其特征在于， S3 中拟合模 型参数的方法为： S301， 根据参数的物理意 义及经验确定一个参数初始值， 包括 W蓄i、 W蓄i‑1、 α月、 k； S302， 计算Ri月， 确定定量描述模型， 建立Ri月计算结果与实测流 量拟合度的目标函数； S303， 通过优化方法确定目标函数的最小值， 得到模型参数α月、 k， 将其作 为设计流域降 雨径流模型参数。 6.根据权利要求5所述无资料地区的逐月径流系列推求方法， 其特征在于， S302中， 所 述目标函数选择枯水期加权均方根误差函数、 残差平方和 函数、 残差绝对值函数或径流百 分比误差函数， S3 03中， 优化方法选择遗传算法或粒子群算法。 7.根据权利要求1所述无资料地区的逐月径流系列推求方法， 其特征在于， S4中逐月降 雨系列为设计流域面平均降水系列， 根据流域内雨量站使用算数平均法或泰森多边形法计 算。权　利　要　求　书 1/1 页 2 CN 113919211 A 2一种无资料地区的逐月径流系列推求方 法 技术领域 [0001]本发明涉及水文学径流过程模拟技术领域， 尤其涉及一种无资料地区的逐月径流 系列 推求方法。 背景技术 [0002]径流系列是水利工程确定规模和工程效益的重要基础。 目前径流设计的主要方法 有水 文比拟法、 水文模型法。 水文比拟法是以参证站的降雨径流为基础， 采用降雨量将参 证流 域的径流比拟至设计流域， 本方法一般均以年降雨量或多年平均降雨量为统计参数 进行比 拟， 无法反应参证流域与设计流域的降雨月分配差异。 现行水文模型种类较多， 基 本都是 基于水量平衡的产汇流模型， 现行水文模型率定和径流计算一般需要降雨、 蒸发、 径流、 植被情况等多种资料， 对资料要求较高， 对于资料缺乏地区的工程设计适用性较差 。 [0003]随着通信技术和自动化设备的发展， 精确的水文数据获取更加便捷， 基于此数据 的水 文学产汇流理论研究更趋 完善， 但是对于资料缺乏地区的径流模拟研究仍然较少。 发明内容 [0004]本发明的目的是提供一种无资料地区的逐月径流系列推求方法， 从而解决现有技 术中 存在的前述问题。 [0005]为实现上述目的， 本发明采用的技术方案如下： 一种无资料地区的逐月径流系列 推求 方法， 包括 步骤如下： [0006]S1， 选取与设计流 域产流特点相似的参证站流 域； [0007]S2， 参证站流域降水资料整理为逐月降水系列， 将参证站径流资料整理为逐月径 流系 列； [0008]S3， 使用S2中整理的逐月降水系列和逐月径流系列， 拟 合模型参数， 确定参证站流   域降雨径流模型； [0009]S4， 整理设计流 域的降雨资料， 统计为逐月降水系列； [0010]S5， 采用S3中确定的降雨径流模型、 拟合的模型参数以及S 4中的逐月降水系列， 推   求设计流 域的逐月径流系列。 [0011]优选， S1中选取的参证站流域与设计流域在下垫面条件、 产汇流类型、 气 候特点特   征上存在相似性。 [0012]优选， S2中降水逐月系列为流域内的面平均降水系列， 根据流域内雨量站使用算 数 平均法或泰 森多边形法计算， 降雨、 径流资料的同步系列长度在10年以上。 [0013]优选， S3参证站流 域的降雨径流模型为： [0014]W蓄i＝[W蓄i‑1+(Pi参‑P损)*α月]*(1‑k)               (1) [0015]Ri月 ＝[W蓄i‑1+(Pi参‑P损)*α月]*k                        (2) [0016]式中： W蓄i、 W蓄i‑1为第i月、 第i ‑1月的月末流域存蓄水量； Pi参为参证站流域第i月的   降水量； P损为月损失水量； α月为月径流系数； k 为月退水系数； Ri月为月径流 量。说　明　书 1/6 页 3 CN 113919211 A 3