一种基于条件概率的水文干旱恢复预估方法技术

本发明专利技术公开了一种基于条件概率的水文干旱恢复预估方法，包括如下步骤：收集统计流域预研究期的水文数据，构建水文干旱异常指数，识别水文干旱事件；确定水文干旱事件恢复时刻，划定水文干旱事件的发展期和恢复期；提取水文干旱事件的发展期时间、缺水量、恢复期时间和恢复水量；选取边缘分布函数计算水文干旱恢复到正常状态时的恢复期时间和恢复水量。本发明专利技术的有益效果：通过确定干旱恢复时刻，划分干旱过程，提出水文干旱恢复预估方案，定量给出水文干旱事件恢复所需的时间和水量，提出的以时间和水量作为水文干旱恢复的阈值指标，方便监测并预警预报，具有更强的实用性和可操作性，为实际应用提供便利。为实际应用提供便利。为实际应用提供便利。

【技术实现步骤摘要】
一种基于条件概率的水文干旱恢复预估方法


[0001]本专利技术涉及一种灾害风险预估领域，特别涉及一种基于条件概率的水文干旱恢复预估方法。

技术介绍

[0002]水文干旱事件发生后，由于气象条件（降水、蒸散发等因素）变化影响造成干旱过程的变化，其过程变化包括为发展期和恢复期。目前，判定干旱发展和恢复期使用最广泛的方法是将干旱过程线简化为“三点一线”，即确定干旱开始、最严重和结束的三个时间点，把三点用光滑曲线连接，寻找干旱最严重时刻，将干旱恢复时刻与最严重时刻等同，根据干旱最严重时刻划分干旱过程为发展期和恢复期，该方法虽然简单，但考虑到实际干旱过程受气象条件影响，其过程存在明显的波动性和一定的滞后性，简化的干旱过程与真实过程存在较大差异，发展期和恢复期判定主观性较强，精度难免受到影响。
[0003]实际上，气象条件波动影响着干旱过程，即干旱的不同过程反映了气象条件的变化，干旱发展期和恢复期也存在不同的气象条件特征，判定干旱发展期和恢复期本质上是确定气象条件与干旱过程之间的关系。此外，目前关于水文干旱的恢复状态判断方法主要通过分析气象水文干旱与水文干旱之间的关系实现，这种方法考虑到不同类型水文干旱间的外在联系，而忽略水文干旱事件自身发展过程的内在联系，且干旱恢复以烈度等无量纲指标判定为主，难以监测，实际操作性不足。

技术实现思路

[0004]针对现有技术存在的不足，本专利技术提供了一种基于条件概率的水文干旱恢复预估方法。
[0005]为了解决上述技术问题，本专利技术采用的技术方案如下：一种基于条件概率的水文干旱恢复预估方法，包括如下步骤：步骤S1，收集统计流域预研究期的水文数据，构建水文干旱异常指数，识别水文干旱事件；步骤S2，确定水文干旱事件恢复时刻，划定水文干旱事件的发展期和恢复期；步骤S3，基于步骤S2划定的水文干旱事件的发展期和恢复期，提取水文干旱事件的发展期时间、缺水量、恢复期时间和恢复水量；步骤S4，选取边缘分布函数拟合水文干旱特征，建立水文干旱恢复条件概率模型，计算水文干旱恢复到正常状态时的恢复期时间和恢复水量；步骤S2具体为：基于时间的一致性得到水文干旱事件和气象要素序列匹配数据，利用相关性分析对水文干旱事件和气象要素的相关性进行诊断，得到相关性统计量序列，获得相关性统计量序列中统计值最小时刻，见公式（1）所示；（1）；
式中：为相关性统计量序列中统计值最小时刻，H
ij
为第j个水文干旱事件第1时刻到第i时刻的数据序列，M
ij
为第j个水文干旱事件第1时刻到第i时刻的数据序列H
ij
对应的气象要素序列，corr表示为求数据序列和气象要素序列的相关性统计量；步骤S3，具体为：基于步骤S2划定的水文干旱事件的发展期和恢复期，统计水文干旱事件的干旱发展期时间和干旱恢复期时间；如公式（2）、（3）所示；
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
（2）；
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
（3）；式中：D1为水文干旱事件的干旱发展期时间，D2为水文干旱事件的干旱恢复期时间，t1为水文干旱事件开始时刻，t
e
为水文干旱事件结束时刻。
[0006]进一步的，步骤S2中，判断并标识每一次水文干旱事件的相关性统计量序列中统计值最小时刻，划定发展期和恢复期的时间点，将相关性统计量序列中统计值最小时刻干旱过程前定义为干旱发展期，将相关性统计量序列中统计值最小时刻干旱过程后定义为干旱恢复期，其中数据序列为水文干旱期间不同时期干旱强度组成的序列，气象要素序列为水文干旱期间不同时期气象要素组成的序列，常见的有蒸散发量序列、气温序列等。
[0007]进一步的，步骤S3中，当时，水文干旱事件的发展期缺水量如公式（4）所示； （4）；式中，WD为水文干旱事件的发展期缺水量，S
i
为i时刻干旱强度，S
i
‑1为i
‑
1时刻干旱强度；当时，水文干旱事件的恢复期恢复水量如公式（5）所示；（5）；式中，RW为水文干旱事件的恢复期恢复水量，为水文干旱事件结束时刻t
e
的干旱强度，为水文干旱事件截取阈值。
[0008]进一步的，步骤S1中收集统计流域预研究期的水文数据，构建水文干旱异常指数，具体为：收集统计流域预研究期的长时间序列n年（年）的水文数据，为保持年水文数据量一致性，剔除闰年2月29日径流量数据，将每年径流量数据从1到365依次排列，计算日时间尺度的水文干旱异常指数如公式（6）所示；（6）；式中：为第n年的水文干旱异常指数，为第1年第1日径流量，为第2年第1日径流量，依次的为第n年第1日径流量，为第1年第2日径流量，依次的为第1年第365日径流量，为流域预研究期n年第1日平均径流量，为流域预研究期n年第2日平均径流量，依次的为流域预研究期n年第365日平均径流量。
[0009]进一步的，步骤S1中识别水文干旱事件，具体为：设置水文干旱事件截取阈值，对流域预研究期的长时间序列n年上的水文干旱指数进行逐日判断，若日水文干旱强度S大于等于0的连续天数大于等于15天则认定发生一次水文干旱事件，依次识别出流域预研究期的长时间序列n年内所有的水文干旱事件；日水文干旱强度的计算如公式（7）所示；（7）；式中，为第n年第i日的干旱强度，为水文干旱事件截取阈值，为第n年第i日的水文干旱异常指数。
[0010]进一步的，步骤S4中选取边缘分布函数拟合水文干旱特征，建立水文干旱恢复条件概率模型，计算水文干旱恢复到正常状态时的恢复期时间和恢复水量；具体为：使用备选的边缘分布函数分别拟合水文干旱发展期时间和缺水量、恢复期时间和恢复水量数据序列，采用最大似然法估计边缘分布函数参数，并根据拟合评估效果选出最佳的边缘分布函数，选用Frank、Gumbel、Clayton copula函数类型分别构建水文干旱发展期时间和恢复期时间、缺水量和恢复水量的联合分布，使用最大似然法估计copula函数参数，采用赤池信息量准则和平方欧式距离评选最优联合分布函数；基于Copula的联合分布和条件概率公式，以水文干旱发展期时间和缺水量为条件，水文干旱恢复期时间和恢复水量为目标，建立水文干旱恢复条件概率模型如公式（8）所示；
[0011] (8)；式中，为条件下事件发生的概率，u为水文干旱发展期时间（缺水量），v为水文干旱恢复期时间（恢复水量），X为水文干旱发展期时间的边缘分布累积概率，Y为水文干旱恢复期时间的边缘分布累积概率；x(u)和y(v)分别为X≤u和Y≤v的累积概率；C(x(u)，y(v))为X≤u和Y≤v的联合累积概率；确定水文干旱恢复的置信水平为95%，即当水文干旱事件恢复概率大于或等于0.95时，水文干旱恢复的置信水平为95%相应的恢复期时间和恢复水量被视为恢复该水文干旱事件的基础，计算水文干旱恢复所需的时间和水量。
[0012]本专利技术通过识别水文干旱事件，提取水文干旱发展期时间和缺水量、恢复期时间和恢复水量，构建水文干旱恢复条件概率模型本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种基于条件概率的水文干旱恢复预估方法，其特征在于：包括如下步骤：步骤S1，收集统计流域预研究期的水文数据，构建水文干旱异常指数，识别水文干旱事件；步骤S2，确定水文干旱事件恢复时刻，划定水文干旱事件的发展期和恢复期；步骤S3，基于步骤S2划定的水文干旱事件的发展期和恢复期，提取水文干旱事件的发展期时间、缺水量、恢复期时间和恢复水量；步骤S4，选取边缘分布函数拟合水文干旱特征，建立水文干旱恢复条件概率模型，计算水文干旱恢复到正常状态时的恢复期时间和恢复水量；步骤S2具体为：基于时间的一致性得到水文干旱事件和气象要素序列匹配数据，利用相关性分析对水文干旱事件和气象要素的相关性进行诊断，得到相关性统计量序列，获得相关性统计量序列中统计值最小时刻，见公式（1）所示；（1）；式中：为相关性统计量序列中统计值最小时刻，H
ij
为第j个水文干旱事件第1时刻到第i时刻的数据序列，M
ij
为第j个水文干旱事件第1时刻到第i时刻的数据序列H
ij
对应的气象要素序列，corr表示为求数据序列和气象要素序列的相关性统计量；步骤S3，具体为：基于步骤S2划定的水文干旱事件的发展期和恢复期，统计水文干旱事件的干旱发展期时间和干旱恢复期时间如公式（2）、（3）所示；（2）；（3）；式中：D1为水文干旱事件的干旱发展期时间，D2为水文干旱事件的干旱恢复期时间，t1为水文干旱事件开始时刻，t
e
为水文干旱事件结束时刻。2.根据权利要求1所述的一种基于条件概率的水文干旱恢复预估方法，其特征在于：步骤S2中，判断并标识每一次水文干旱事件的相关性统计量序列中统计值最小时刻，划定发展期和恢复期的时间点，将相关性统计量序列中统计值最小时刻干旱过程前定义为干旱发展期，将相关性统计量序列中统计值最小时刻干旱过程后定义为干旱恢复期，其中数据序列为水文干旱期间不同时期干旱强度组成的序列，气象要素序列为水文干旱期间不同时期气象要素组成的序列。3.根据权利要求2所述的一种基于条件概率的水文干旱恢复预估方法，其特征在于：步骤S3中，当时，水文干旱事件的发展期缺水量如公式（4）所示；（4）；式中，WD为水文干旱事件的发展期缺水量，S
i
为i时刻干旱强度，S
i
‑1为i
‑
1时刻干旱强度；当时，水文干旱事件的恢复期恢复水量如公式（5）所示； （5）；
式中，RW为水文干旱事件的恢复期恢复水量，为水文干旱事件结束时刻t
e
的干旱强度，为水文干旱事件截取阈值。4.根据权利要求1所述的一种基于条件概率的水文...

【专利技术属性】
技术研发人员：韩会明，周王莹，雷声，孙军红，王农，简鸿福，龙鹏，刘明超，
申请(专利权)人：江西省水利科学院江西省大坝安全管理中心，江西省水资源管理中心，
类型：发明
国别省市：