(19)国家知识产权局 (12)发明 专利 (10)授权公告 号 (45)授权公告日 (21)申请 号 202210647514.0 (22)申请日 2022.06.09 (65)同一申请的已公布的文献号 申请公布号 CN 114722645 A (43)申请公布日 2022.07.08 (73)专利权人 广东省科 学院广州地理研究所 地址 510075 广东省广州市越秀区先烈中 路100号大院之一 (72)发明人 徐飞　赵玲玲　陈彩霞　王钧　 宫清华　 (74)专利代理 机构 广州骏思知识产权代理有限 公司 44425 专利代理师 潘桂生 (51)Int.Cl. G06F 30/20(2020.01)G06F 113/08(2020.01) (56)对比文件 CN 1079089 22 A,2018.04.13 CN 1094720 04 A,2019.0 3.15 CN 10746 3730 A,2017.12.12 CN 102567635 A,2012.07.1 1 郭军庭等.气候和土地利用变化对潮河流 域 产流产沙的影响. 《农业工程学报》 .2012,(第14 期),正文第237-243页. 刘士余等.降雨与植被 变化对彭冲涧小流 域 径流影响的定量研究. 《江西农业大 学学报》 .2018,(第01期),正文第189-197页. 刘卉芳等.黄土区小流 域水沙对降雨及土地 利用变化响应 研究. 《干旱 地区农业研究》 .2010, (第02期),正文第237-242页. 审查员 翟紫伶 (54)发明名称 影响径流变化的因素贡献率获取方法、 系 统、 介质和设备 (57)摘要 本发明提供一种影响径流变化的因素贡献 率获取方法、 系统、 介质和设备， 所获取方法包 括： 获取待研究流域的数字高程模 型和水文气象 等数据； 分析径流演变规律， 并进行突变检验； 构 建待研究流域的WEP ‑L分布式水文模型； 对径流 变化期的各气象要素进行去趋势处理； 设置基准 情景、 去气象要素变化趋势情景、 土地利用变化 情景； 将各情景输入到WEP ‑L模型进行模拟， 得到 各情景径 流， 将去气象要素变化趋势情景和土地 利用变化情景与基准情景进行对比分析， 量化各 气象要素、 土地利用变化和直接人类活动对径流 变化的贡献率。 从而获得多个气象要素 （包括降 水、 气温、 日照时数、 相对湿度、 平均风速） 、 土地 利用变化和直接人类活动对径流变化的影响的 贡献率。 权利要求书4页 说明书15页 附图4页 CN 114722645 B 2022.08.16 CN 114722645 B 1.一种影响径流变化的因素贡献率获取 方法， 其特 征在于， 包括： 获取待研究流域内及预设的周边地区气象站的逐日气象数据， 包括降水、 气温、 日照时 数、 相对湿度、 平均风速， 待研究流域的数字高程模型、 土地利用图、 土壤类型图、 水系图、 水 文站实测逐日径流； 依据待研究流域水文站 的实测径流长系列资料， 分析年径流系列演变规律， 并对径流 系列进行突变检验， 依据 径流突变时间点将突变前后径流系列划分为径流基准期和径流变 化期； 根据所述逐日气象数据、 所述待研究流域的数字高程模型、 所述土地利用图、 所述土壤 类型图和所述水系图构建WEP ‑L分布式水文模型， 并采用径流基准期的实测径流对所述 WEP‑L分布式水文模型进行率定和验证； 分析所述径流变化期的各个气象要素的变化趋势， 对所述径流变化期的各个所述气象 要素的实际逐日数据进行去趋势处 理， 得到径流变化期各个气象要素的去趋势逐日数据； 将各个所述气象要素的实 际逐日数据和所述径流基准期的土地利用图确定为基准情 景数据； 根据各个气象要素 的去趋势逐日数据、 各个所述气象要素 的实际逐日数据和所述 径流基准期的土地利用图生成多组去气象要 素变化趋势情景； 根据各个所述气象要素的实 际逐日数据和所述径流变化期的土地利用图， 生成土地利用变化情景； 将所述基准情景数据、 各组所述去气象要素变化趋势情景和所述土地利用变化情景输 入到所述WEP ‑L分布式水文模型进行模拟， 得到所述基准情景以及各组所述去气象要素变 化趋势情景和所述土地利用变化情景的模拟径流； 根据所述基准情景、 去气象要素变化趋势情景、 土地利用变化情景的模拟径流， 以及研 究流域的实际径流变化量， 量化气候变化和人类活动对径流变化的贡献率， 以及各个所述 气象要素、 土地利用变化和直接人类活动对径流变化的贡献率。 2.根据权利要求1所述的影响径流变化的因素贡献率获取方法， 其特征在于， 所述分析 所述径流变化期的各个气象要 素的变化趋势， 对所述径流变化期的各个所述气象要素的实 际逐日数据进行去趋势处理， 得到径流变化期各个气象要素 的去趋势逐日数据的步骤， 包 括： 通过以下公式， 计算出变化期各个所述气象要素的径流变化期 气象序列： ； ； 其中， 为变化期各个所述气象要素的径流变化期气 象序列， 为气象序列年长度时 间节点， 为气象序列年长度， 为气象序列线性回归方程斜率， 为气象序列线性回归方 程截距； 通过以下公式， 计算出变化期各个所述气象要素的年去趋势气象序列： ； ； 其中， 为变化期各个所述气象要素的年去趋势气象序列； 通过以下公式， 计算出变化期各个所述气象要素的去趋势逐日数据：权　利　要　求　书 1/4 页 2 CN 114722645 B 2； 其中， 为变化期各个所述气象要 素的去趋势逐日数据， 为变化期各个所述气 象要素的实际逐日数据。 3.根据权利要求1所述的影响径流变化的因素贡献率获取方法， 其特征在于， 所述依据 待研究流域水文站的实测径流长系列资料， 分析年径流系列演变规律， 并对径流系列进行 突变检验， 依据 径流突变时间点将突变前后径流系列划分为径流基准期和径流变化期的步 骤， 包括： 通过以下公式， 得到所述待研究流 域各个年份的径流变化 量序列统计结果： ； ； ； 其中， 为所述待研究流域各个年份的径流变化量序列统计结果； 和 为欲 进行假设检验的径流随机变量， 为径流序列年长度时间节点， 为径流序列年长度； 通过以下公式， 得到径流 突变时间点 参考值； ； 其中， 为所述径流序列年长度时间节点的径流 突变时间点 参考值； 根据小于预设的径流突变时间点阈值的径流突变时间点参考值， 将对应的所述径流序 列年长度时间节点确定为径流 突变时间点； 将属于所述径流突变时间点之前的时期确定为径流基准期， 将属于所述径流突变时间 点之后的时期确定为径流变化期。 4.根据权利要求1所述的影响径流变化的因素贡献率获取方法， 其特征在于， 计算出气 候变化对径流变化的贡献率， 包括： 根据各个所述气象要素的对径流变化的影响量之和， 得到气候变化对径流变化的影响 量； 其中， 所述气候变化对径流变化的影响量为各个所述气象要素 的对径流变化的影响量 之和； 根据所述实际径流变化量和所述气候变化对径流变化的影响量， 得到所述气候变化对 径流变化的贡献率。 5.根据权利要求4所述的影响径流变化的因素贡献率获取方法， 其特征在于， 计算出人 类活动对径流变化的贡献率， 包括：权　利　要　求　书 2/4 页 3 CN 114722645 B 3