申请/专利权人：成都理工大学;中国铁路设计集团有限公司;中国国家铁路集团有限公司

申请日：2023-09-05

公开（公告）日：2023-11-17

公开（公告）号：CN116882204B

主分类号：G06F30/20

分类号：G06F30/20;G06F113/08

优先权：

专利状态码：有效-授权

法律状态：2023.11.17#授权;2023.10.31#实质审查的生效;2023.10.13#公开

摘要：本发明公开了一种无实测径流资料地区暴雨‑融雪洪峰径流强度估算方法，属于水文预报技术领域，包括确定研究区峰值径流时间节点A；利用GIS栅格化数据图层并划分不同的产流单元类型；结合研究区地表覆盖类型数据与土壤质地数据确定研究区不同水文响应单元，并结合不同水文响应单元的CN值分别计算不同类型产流单元的产流量；利用计算得到的净产流体，结合研究区的DEM地形栅格数据，基于MassFlow模拟获得时间节点A下地表暴雨‑融雪产汇流时空演进结果。本方法通过重点描述无资料地区水文的产流和汇流过程，达到了简化模型结构、减少参数需求的目的，可估算地表多年降雨‑融雪峰值径流的稳定状态值。

主权项：1.一种无实测径流资料地区暴雨-融雪洪峰径流强度估算方法，其特征在于，包括如下步骤：步骤一、结合研究区的国家气象站温度降雨监测数据和积雪覆盖数据，获得暴雨期和融雪期，并根据其时间轴分布特征确定峰值径流时间节点A；步骤二、利用GIS栅格化数据图层并统一栅格大小，按照不同海拔高度降水-产流形式将栅格单元划分为不同的产流类型，包括降雨-融雪产流单元、降雨产流单元以及降雪不产流单元；包括：统计每个栅格单元多年在时间节点A时的积雪覆盖率；基于GIS的空间分析得到统计频率和大于该统计频率积雪覆盖率的多项式拟合模型；通过拟合曲线确定多年时间节点A的地区平均积雪覆盖率所对应的统计频率；以超过该统计频率的栅格区域代表研究区在时间节点A的多年积雪覆盖平均范围，并以此为依据划分降雨产流单元和降雨-融雪产流单元；以国家气象检测站多年时间节点A平均气温数据为基点，结合研究区数字高程数据，根据气温垂直递减率获取研究区地表气温分布图，以0℃作为降水形态的分界线，以此为依据划分降雨-融雪产流单元和降雪不产流单元；步骤三、结合研究区地表覆盖类型数据集细化无积雪区域的降雨产流单元类型，修正时间节点A时的地表覆盖类型数据；步骤四、将修正后的地表覆盖类型数据与研究区的土壤质地数据相结合，确定研究区不同水文响应单元，并结合不同水文响应单元的CN值分别计算不同类型产流单元的产流量；步骤五、利用步骤四计算得到的净产流体，结合研究区的DEM地形栅格数据，基于MassFlow模拟获得时间节点A下地表暴雨-融雪产汇流时空演进结果；包括：基于步骤四计算得到的净产流体数据生成流体深度文件h，基于研究区的DEM地形栅格数据生成三维地形文件z；构建z+h三维数值计算模型，基于MassFlow实现降水净产流体的三维运动模拟，获得研究区地表径流的时空演化特征，包括不同模拟过程中不同时间节点的径流速度以及过程中最大径流深度以及径流强度；不同类型产流单元的产流量按照下式分别计算： ；式中：表示径流量，mm；表示降雨量，mm；表示融雪当量，mm；表示初损量，即发生径流前的初损失量，mm；表示流域当时的最大可能保持量，mm；降雨-融雪产流单元满足；降雨产流单元满足；降雪不产流单元为积雪区； 与之间的经验关系应满足：；；式中，表示初损系数；位于不同高程的降雨-融雪产流单元融雪当量按照下式计算： ；式中：表示不同产流单元24小时融雪当量（mm），i表示产流单元编号；表示冰川或雪的度日因子（mm（d*℃））；表示不同产流单元气温；表示流域积雪面积（km²）；表示流域总面积（km²）；度日因子DDF考虑了坡度、海拔、纬度的影响，其计算公式如下： ；式中：表示产流单元海拔高度（m）；表示产流单元纬度（°）；表示产流单元坡度（°）；不同水文响应单元的CN值为校正后的湿润条件下的土壤参数值，校正公式如下： ；式中：为查表所得的中等湿润条件下的土壤参数值；为湿润条件下的土壤参数值。

全文数据：

权利要求：

百度查询： 成都理工大学;中国铁路设计集团有限公司;中国国家铁路集团有限公司 一种无实测径流资料地区暴雨-融雪洪峰径流强度估算方法

免责声明
1、本报告根据公开、合法渠道获得相关数据和信息，力求客观、公正，但并不保证数据的最终完整性和准确性。
2、报告中的分析和结论仅反映本公司于发布本报告当日的职业理解，仅供参考使用，不能作为本公司承担任何法律责任的依据或者凭证。

阅读全文 双屏查看 官方信息 专利公告 收藏专利 下载PDF 下载WORD