CN108920427B - 一种月尺度水文模型时变参数的动态估计方法 - Google Patents
一种月尺度水文模型时变参数的动态估计方法 Download PDFInfo
- Publication number
- CN108920427B CN108920427B CN201810557800.1A CN201810557800A CN108920427B CN 108920427 B CN108920427 B CN 108920427B CN 201810557800 A CN201810557800 A CN 201810557800A CN 108920427 B CN108920427 B CN 108920427B
- Authority
- CN
- China
- Prior art keywords
- month
- scale
- data
- monthly
- hydrological model
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Data Mining & Analysis (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Physics (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Life Sciences & Earth Sciences (AREA)
- Operations Research (AREA)
- Probability & Statistics with Applications (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Algebra (AREA)
- Evolutionary Biology (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Bioinformatics & Computational Biology (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种月尺度水文模型时变参数的动态估计方法,步骤为:根据实测水文数据,采用数据同化方法估计月尺度水文模型参数θ的连续序列;通过相关分析从候选的流域特征表征因子筛选出月尺度水文模型参数θ的流域特征表征因子;拟合月尺度水文模型参数θ与其流域特征表征因子的经验公式;根据经验公式以及流域特征因子的月尺度数据,计算流域月尺度水文模型的参数序列。本发明能够实现月尺度水文模型时变参数的动态估计,可提高变化环境下“非稳态”流域的径流模拟精度。
Description
技术领域
本发明属于流域水文模型领域,特别涉及了一种月尺度水文模型时变参数的动态估计方法。
背景技术
流域水文模型是研究流域水文自然规律和解决水文实际问题的重要工具。它基于数学物理方程和参数对流域水文循环过程进行概化,构建流域降水——径流的转换关系,实现流域径流的模拟和预报,从而为流域洪水预报、水资源科学利用与管理等提供技术支撑。
现有的流域水文模型参数估计方法一般根据水文实测资料采用优化算法进行率定,其假定条件为模型参数在流域内为常数,即参数不随时间变化。现行方法的主要实施步骤为:(1)选取完整的实测水文资料序列,如降雨、潜在蒸发、径流资料等;(2)选取优化参数的目标函数,一般采用模拟径流与实测径流的误差平方和最小为目标函数;(3)采用优化算法对水文模型的待估计参数进行优选,其中较为常用的优化算法有SCE-UA算法、基因算法等。
因此,现有的水文模型参数估计方法存在问题:(1)水文模型的参数被当作常数,无法准确地反映变化环境背景下流域特征条件所呈现的动态变化;(2)模型参数依赖于足够的实测水文资料序列,通过优化算法反演的方式获得,且估计值为固定值,无法实现参数时间序列的连续估计。
发明内容
为了解决上述背景技术提出的技术问题,本发明旨在提供一种月尺度水文模型时变参数的动态估计方法,在水文模型中更好地反映流域特征条件变化,提高变化环境下“非稳态”流域的径流模拟精度。
为了实现上述技术目的,本发明的技术方案为:
一种月尺度水文模型时变参数的动态估计方法,包括以下步骤:
(1)根据实测水文数据,采用数据同化方法估计月尺度水文模型参数θ的连续序列;
(2)通过相关分析从候选的流域特征表征因子筛选出月尺度水文模型参数θ的流域特征表征因子;
(3)拟合月尺度水文模型参数θ与其流域特征表征因子的经验公式;
(4)根据步骤(3)得到的经验公式以及流域特征因子的月尺度数据,计算流域月尺度水文模型的参数序列。
进一步地,在步骤(1)中,选取有足够水文资料的流域,将降水、潜在蒸发和径流处理为月尺度数据,确定月尺度水文模型,建立基于月尺度水文模型的数据同化状态转移方程和观测方程:
yt+1=h(xt+1,θt+1)+ξ
上式中,θt+1、θt分别为第t+1、t月的模型参数值,ε为均值为0的正态分布误差;xt+1、xt分别为第t+1、t月的模型状态变量,η为均值为0的正态分布误差;ut+1为第t+1月的模型驱动数据,包括降雨和潜在蒸散量;yt+1为第t+1月的模型径流模拟值,ξ为均值为0的正态分布误差;f、h均表示水文模型。
进一步地,在步骤(2)中,候选的流域特征表征因子包括1、3、6个月前期影响雨量P1、P3、P6,归一化植被指数NDVI以及经济社会指标GDP、人口Pop。
进一步地,为了消除数据同化估计的参数序列的波动以及综合考虑流域特征数据的可利用性,在流域特征表征因子筛选的过程中:
(a)由于估计的月尺度水文模型参数θ的序列以及1、3、6个月前期影响雨量P1、P3、P6,归一化植被指数NDVI数据均为月尺度数据序列,即每月有一个数值,将它们处理为多年月平均值,即对多年的每个月份进行平均:
(b)经济社会指标GDP以及人口Pop为年尺度数据,将估计的月尺度水文模型参数θ的序列处理为年尺度的月平均值,即分别对每年进行月平均计算:
上式中,i表示月份,i=1,2,…,12,j表示第j年,j=1,2,…,m,m表示资料序列的总年数;再分别对θj与GDPj和Popj进行相关分析。
进一步地,在步骤(2)中,计算月尺度水文模型参数θ与各候选的流域特征表征因子之间的Spearman秩相关系数rs,筛选出|rs|>0.6的候选流域特征表征因子作为最终的流域特征表征因子。
进一步地,在步骤(3)中,采用多元线性回归方法拟合月尺度水文模型参数θ与其流域特征表征因子的经验公式:
当步骤(2)筛选出的最终的流域特征表征因子为1、3、6个月前期影响雨量或/和归一化植被指数NDVI时,经验公式如下:
θi=aWi+b
或者
θi=aWi 2+bWi+c
当步骤(2)筛选出的最终的流域特征表征因子为1、3、6个月前期影响雨量或/和归一化植被指数NDVI以及经济社会指标GDP或/和人口Pop时,经验公式如下:
θi=aWi+bZj+c
或者
θi=aWi 2+bWi+cZj+d
上式中,θi表示θ的多年月平均值;Wi表示1、3、6个月前期影响雨量或/和归一化植被指数NDVI的多年月平均值;Zj表示经济社会指标GDP或/和人口Pop的年尺度的月平均值;a、b、c、d为待定系数。
采用上述技术方案带来的有益效果:
(1)现有技术一般仅能考虑流域气候条件和下垫面条件不发生变化的情况,认为模型的参数为常数,而本发明方法考虑了水文模型参数随时间变化的特性,更能准确地反映变化环境下流域特征条件的变化;
(2)现有技术需要较为详细和足够长度的实测水文资料序列进行参数估计,本发明方法根据月尺度水文模型时变参数与流域特征因子之间的经验关系,构建以流域特征表征因子为自变量的时变参数函数型式,无需通过参数率定来进行月尺度模型参数的估计,可为无(缺)资料地区的月径流模拟提供参考依据。
附图说明
图1是本发明的基本流程图。
具体实施方式
以下将结合附图,对本发明的技术方案进行详细说明。
本发明根据月尺度水文模型参数估计值序列,采用统计方法建立时变参数与流域特征表征因子的经验公式,然后基于构建的经验公式和流域特征数据计算月尺度水文模型参数。如图1所示,本发明的技术方案具体包括以下步骤:
步骤1,根据实测水文数据,采用数据同化方法估计月尺度水文模型参数θ的连续序列。
选取有足够水文资料的流域,将降水、潜在蒸发和径流处理为月尺度数据,确定月尺度水文模型。建立基于水文模型的数据同化状态转移方程和观测方程,如下:
yt+1=h(xt+1,θt+1)+ξ (2)
式中:θt+1、θt分别为第t+1、t月的模型参数值,ε为均值为0的正态分布误差;xt+1、xt分别为第t+1、t月的模型状态变量,η为均值为0的正态分布误差;ut+1为第t+1月的模型驱动数据,一般包括降雨P、潜在蒸散发PET;yt+1为第t+1月的模型径流模拟值,ξ为均值为0的正态分布误差;f、h在此处均表示水文模型。
步骤2,采用统计方法分析参数与流域特征因子的相关关系,筛选出月尺度水文模型参数θ的流域特征表征因子W或/和Z。
选取与月尺度水文模型参数θ具有潜在物理关联的流域特征因子,包括1、3、6个月前期影响雨量P1、P3、P6,归一化植被指数NDVI以及经济社会指标GDP、人口Pop;采用统计指标Spearman秩相关系数rs分析模型时变参数序列与流域特征因子的相关关系,筛选出与月尺度水文模型时变参数有较强相关关系的流域特征因子。当|rs|>0.6时认为两种要素之间具有较强相关关系。
为了消除数据同化估计的参数值序列的波动以及综合考虑流域特征数据的可利用性,在流域特征因子筛选的过程中:
(1)由于获取的月尺度水文模型参数估计值序列以及1、3、6个月前期影响雨量、归一化植被指数NDVI数据均为月尺度数据序列,即每月有一个数值。因此,将它们处理为多年月平均值,即对多年的每个月份进行平均, i表示月份,i=1,2,…,12,j表示第j年,j=1,2,…,m,m表示资料序列的总年数。随后分别对θi与以及NDVIi进行相关分析。
(2)经济社会指标GDP以及人口Pop一般为年尺度数据。因此,将月尺度水文模型参数估计值序列处理为年尺度的月平均值,即分别对每年进行月平均计算,i表示月份,i=1,2,…,12,j表示第j年,j=1,2,…,m,m表示资料序列的总年数。随后分别对θj与GDPj和Popj进行相关分析。
步骤3,采用多元线性回归方法拟合月尺度水文模型参数θ分别与挑选出的流域特征表征因子之间的经验公式。
当θi仅与前期影响雨量或/和归一化植被指数呈现较强相关关系时,经验公式如下:
θi=aWi+b (3)
或者
θi=aWi 2+bWi+c (4)
当θi与前期影响雨量或/和归一化植被指数以及经济社会指标均呈现较强相关关系时,经验公式如下:
θi=aWi+bZj+c (5)
或者
θi=aWi 2+bWi+cZj+d (6)
式中:θi表示θ的多年月平均值;Wi表示与参数有较强相关关系的前期影响雨量或(和)归一化植被指数的多年月平均值;Zj表示与参数有较强相关关系的经济社会指标GDP或/和人口Pop的年尺度的月平均值;a、b、c、d表示函数表达式的待定系数。
步骤4,根据步骤3中的经验公式以及流域特征因子的月尺度数据Wt′或/和Z′,计算流域月尺度水文模型的参数序列θt′。
当构建的时变参数函数型式的自变量为前期影响雨量、归一化植被指数时,月尺度水文模型时变参数的计算公式为:
θt′=a′Wt′+b′ (7)
或者
θt′=a′Wt′2+b′Wt′+c′ (8)
当构建的时变参数函数型式的自变量为前期影响雨量、归一化植被指数以及经济社会指标时,月尺度水文模型时变参数的计算公式为:
θt′=a′Wt′+b′Z′k+c′ (9)
或者
θt′=a′Wt′2+b′Wt′+c′Z′k+d′ (10)
式中:θt′表示待估计的月尺度水文模型时变参数序列;Wt′表示研究流域中时变参数流域特征表征因子数据序列,即前期影响雨量或/和归一化植被指数;Zt′表示研究流域中时变参数流域特征表征因子数据序列,即经济社会指标GDP或/和人口Pop;a′、b′、c′、d′表示步骤3中已估计出的经验公式的待定系数。
实施例仅为说明本发明的技术思想,不能以此限定本发明的保护范围,凡是按照本发明提出的技术思想,在技术方案基础上所做的任何改动,均落入本发明保护范围之内。
Claims (2)
1.一种月尺度水文模型时变参数的动态估计方法,其特征在于,包括以下步骤:
(1)根据实测水文数据,采用数据同化方法估计月尺度水文模型参数θ的连续序列;
选取有足够水文资料的流域,将降水、潜在蒸发和径流处理为月尺度数据,确定月尺度水文模型,建立基于月尺度水文模型的数据同化状态转移方程和观测方程:
yt+1=h(xt+1,θt+1)+ξ
上式中,θt+1、θt分别为第t+1、t月的模型参数值,ε为均值为0的正态分布误差;xt+1、xt分别为第t+1、t月的模型状态变量,η为均值为0的正态分布误差;ut+1为第t+1月的模型驱动数据,包括降雨和潜在蒸散量;yt+1为第t+1月的模型径流模拟值,ξ为均值为0的正态分布误差;f、h均表示水文模型;
(2)通过相关分析从候选的流域特征表征因子筛选出月尺度水文模型参数θ的流域特征表征因子;
候选的流域特征表征因子包括1、3、6个月前期影响雨量P1、P3、P6,归一化植被指数NDVI以及经济社会指标GDP、人口Pop;
为了消除数据同化估计的参数序列的波动以及综合考虑流域特征数据的可利用性,在流域特征表征因子筛选的过程中:
(a)由于估计的月尺度水文模型参数θ的序列以及1、3、6个月前期影响雨量P1、P3、P6,归一化植被指数NDVI数据均为月尺度数据序列,即每月有一个数值,将它们处理为多年月平均值,即对多年的每个月份进行平均:
(b)经济社会指标GDP以及人口Pop为年尺度数据,将估计的月尺度水文模型参数θ的序列处理为年尺度的月平均值,即分别对每年进行月平均计算:
上式中,i表示月份,i=1,2,…,12,j表示第j年,j=1,2,…,m,m表示资料序列的总年数;再分别对θj与GDPj和Popj进行相关分析;
(3)拟合月尺度水文模型参数θ与其流域特征表征因子的经验公式;
采用多元线性回归方法拟合月尺度水文模型参数θ与其流域特征表征因子的经验公式:
当步骤(2)筛选出的最终的流域特征表征因子为1、3、6个月前期影响雨量或/和归一化植被指数NDVI时,经验公式如下:
θi=aWi+b
或者
θi=aWi 2+bWi+c
当步骤(2)筛选出的最终的流域特征表征因子为1、3、6个月前期影响雨量或/和归一化植被指数NDVI以及经济社会指标GDP或/和人口Pop时,经验公式如下:
θi=aWi+bZj+c
或者
θi=aWi 2+bWi+cZj+d
上式中,Wi表示1、3、6个月前期影响雨量或/和归一化植被指数NDVI的多年月平均值;Zj表示经济社会指标GDP或/和人口Pop的年尺度的月平均值;a、b、c、d为待定系数;
(4)根据步骤(3)得到的经验公式以及流域特征因子的月尺度数据,计算流域月尺度水文模型的参数序列。
2.根据权利要求1所述月尺度水文模型时变参数的动态估计方法,其特征在于,在步骤(2)中,计算月尺度水文模型参数θ与各候选的流域特征表征因子之间的Spearman秩相关系数rs,筛选出|rs|>0.6的候选流域特征表征因子作为最终的流域特征表征因子。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810557800.1A CN108920427B (zh) | 2018-06-01 | 2018-06-01 | 一种月尺度水文模型时变参数的动态估计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810557800.1A CN108920427B (zh) | 2018-06-01 | 2018-06-01 | 一种月尺度水文模型时变参数的动态估计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108920427A CN108920427A (zh) | 2018-11-30 |
CN108920427B true CN108920427B (zh) | 2020-07-24 |
Family
ID=64418046
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810557800.1A Active CN108920427B (zh) | 2018-06-01 | 2018-06-01 | 一种月尺度水文模型时变参数的动态估计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108920427B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109636030B (zh) * | 2018-12-10 | 2023-05-26 | 浙江大学 | 一种基于分形理论的水文模型多目标参数优化方法 |
CN109492342B (zh) * | 2018-12-25 | 2022-08-19 | 河海大学 | 水文模型参数空间尺度效应分析方法 |
CN112883558B (zh) * | 2021-01-27 | 2022-04-26 | 长江水利委员会水文局 | 一种水文模型参数时变形式构造方法 |
CN112800636B (zh) * | 2021-04-08 | 2021-07-02 | 水利部交通运输部国家能源局南京水利科学研究院 | 一种估算无资料地区流域地表水资源量的方法及系统 |
CN114357737B (zh) * | 2021-12-21 | 2024-04-05 | 武汉大学 | 针对大尺度水文模型时变参数的代理优化率定方法 |
CN116932990B (zh) * | 2023-09-14 | 2024-01-09 | 南方科技大学 | 一种基于转换函数方法的水文地质参数估计方法及系统 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104765981A (zh) * | 2015-05-07 | 2015-07-08 | 中国科学院地理科学与资源研究所 | 一种基于分布式时变增益水文模型的蒸散发数据同化方法 |
CN107169204A (zh) * | 2017-05-17 | 2017-09-15 | 武汉大学 | 不同时间尺度上水文模型参数转换的计算方法 |
WO2018003621A1 (ja) * | 2016-06-29 | 2018-01-04 | 学校法人東京理科大学 | 河川水位流量計算装置及びそのプログラム |
CN107818238A (zh) * | 2017-09-28 | 2018-03-20 | 河海大学 | 一种确定蒸散发变化主因及判别因素间耦合作用的方法 |
-
2018
- 2018-06-01 CN CN201810557800.1A patent/CN108920427B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104765981A (zh) * | 2015-05-07 | 2015-07-08 | 中国科学院地理科学与资源研究所 | 一种基于分布式时变增益水文模型的蒸散发数据同化方法 |
WO2018003621A1 (ja) * | 2016-06-29 | 2018-01-04 | 学校法人東京理科大学 | 河川水位流量計算装置及びそのプログラム |
CN107169204A (zh) * | 2017-05-17 | 2017-09-15 | 武汉大学 | 不同时间尺度上水文模型参数转换的计算方法 |
CN107818238A (zh) * | 2017-09-28 | 2018-03-20 | 河海大学 | 一种确定蒸散发变化主因及判别因素间耦合作用的方法 |
Non-Patent Citations (2)
Title |
---|
关于气候变化与人类活动对径流影响研究的评述;刘春蓁 等;《水利学报》;20140430;第45卷(第4期);第379-393页 * |
双参数月水量平衡模型在汉江流域上游的应用;乐通潮 等;《资源科学》;20041130;第26卷(第6期);第97-103页 * |
Also Published As
Publication number | Publication date |
---|---|
CN108920427A (zh) | 2018-11-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108920427B (zh) | 一种月尺度水文模型时变参数的动态估计方法 | |
CN101480143B (zh) | 一种预测灌区作物单产量的方法 | |
CN108304668B (zh) | 一种结合水文过程数据和历史先验数据的洪水预测方法 | |
Liu et al. | Coupling the k-nearest neighbor procedure with the Kalman filter for real-time updating of the hydraulic model in flood forecasting | |
Nourani et al. | Estimation of prediction interval in ANN-based multi-GCMs downscaling of hydro-climatologic parameters | |
Matsumura et al. | Maize yield forecasting by linear regression and artificial neural networks in Jilin, China | |
CN108876021B (zh) | 一种中长期径流预报方法及系统 | |
Zhang et al. | Explicitly integrating parameter, input, and structure uncertainties into Bayesian Neural Networks for probabilistic hydrologic forecasting | |
CN107169204B (zh) | 不同时间尺度上水文模型参数转换的计算方法 | |
Ye et al. | Hydrologic post-processing of MOPEX streamflow simulations | |
CN112288164B (zh) | 一种计及空间相关性和修正数值天气预报的风功率组合预测方法 | |
CN103235096A (zh) | 污水水质检测方法与装置 | |
CN112883558B (zh) | 一种水文模型参数时变形式构造方法 | |
CN114254802B (zh) | 气候变化驱动下植被覆盖时空变化的预测方法 | |
CN113887972A (zh) | 一种基于水文过程的综合旱情监测与评估方法 | |
Zhang et al. | Bayesian neural networks for uncertainty analysis of hydrologic modeling: a comparison of two schemes | |
Samantaray et al. | Sediment assessment for a watershed in arid region via neural networks | |
CN114357737B (zh) | 针对大尺度水文模型时变参数的代理优化率定方法 | |
Rocha et al. | Early prediction of durum wheat yield in Spain using radial basis functions interpolation models based on agroclimatic data | |
CN109190810B (zh) | 基于tdnn的中国北方草原地区ndvi的预测方法 | |
CN109272144B (zh) | 基于bpnn的中国北方草原地区ndvi的预测方法 | |
CN111311026A (zh) | 一种顾及数据特征、模型和校正的径流非线性预测方法 | |
CN110851897A (zh) | 一种多因素关联下的渡槽应力应变预测方法 | |
Maqsood et al. | Statistical downscaling and projection of climatic extremes using machine learning algorithms | |
Liu et al. | Unravelling and improving the potential of global discharge reanalysis dataset in streamflow estimation in ungauged basins |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |