水文气象灾害二维复合事件模拟方法、装置、设备及介质制造方法及图纸

本申请涉及一种水文气象灾害二维复合事件模拟方法、装置、设备及介质，其方法包括：采集水文气象灾害二维复合事件中两个变量时间序列数据；对每个变量时间序列数据随时间变化及其之间的变化趋势进行量化；基于趋势量化结果构建一致性和非一致性计算方法，确定最佳拟合分布函数；对每个变量时间序列数据自身的时变性及其之间的相关性进行检验；基于相关性检验结果描述两个变量时间序列数据之间相关性的最佳Copula函数；基于最佳拟合分布函数和最佳Copula函数，构建两个变量的联合概率分布函数；基于联合概率分布函数的联合超越概率的分位数确定水文气象灾害二维复合事件的重现期。本申请能够准确模拟气候变化下的复杂复合事件。件。件。

【技术实现步骤摘要】
水文气象灾害二维复合事件模拟方法、装置、设备及介质


[0001]本申请涉及水文气象模拟
，尤其是涉及一种水文气象灾害二维复合事件模拟方法、装置、设备及介质。

技术介绍

[0002]水文气象灾害二维复合事件，可定义为具有一定时空关联关系的可造成严重社会环境风险的两个过程或致灾因子的组合，例如暴雨
‑
风暴潮遭遇下的洪水事件、高温
‑
无雨遭遇下的干旱事件等，其发生频率一般采用统计模型方法。
[0003]近年来，联合概率模型发展迅速，成为水文气象灾害复合事件研究的主要统计建模方法。联合概率模型中最常见和典型的是Copula模型，它是一种连接多维联合分布及其边缘分布的函数，也是一种灵活、稳健的非线性相关性分析工具。联合概率Copula模型模拟主要分为两个部分：单变量(例如暴雨、风暴潮)的边缘分布模拟和关联性(例如暴雨
‑
风暴潮)模拟。
[0004]随着气候变化的加剧，形成复合事件的变量都不再是一成不变的，它们都在气候变化或社会经济变化的作用下发生变化，且不仅自身会变化，它们之间的关联性也有可能发生变化，例如暴雨和风暴潮的相关性变得更强或更弱。
[0005]而现有的Copula模型都是基于“一致性”假设的，即假设仅单一变量的边缘分布是时变的，或者仅变量之间的关联性是时变的，缺乏二者均变化的模型框架，只能模拟单一要素变化，或者只能模拟要素间关联性变化，无法模拟要素和关联性都变化的情况，因此，无法准确模拟气候变化下的复杂复合事件，例如暴雨和风暴潮共同作用下产生的洪水事件。

技术实现思路

[0006]为了解决上述技术问题，本申请提供一种水文气象灾害二维复合事件模拟方法、装置、设备及介质。
[0007]第一方面，本申请提供一种水文气象灾害二维复合事件模拟方法，采用如下的技术方案：
[0008]一种水文气象灾害二维复合事件模拟方法，包括：
[0009]采集水文气象灾害二维复合事件中两个变量时间序列数据；
[0010]对每个变量时间序列数据随时间变化以及所述两个变量时间序列数据之间的变化趋势进行量化，生成趋势量化结果；
[0011]基于所述趋势量化结果构建一致性和非一致性计算方法，确定每个变量时间序列数据的最佳拟合分布函数；
[0012]对每个变量时间序列数据自身的时变性以及所述两个变量时间序列数据之间的相关性进行检验，生成相关性检验结果；
[0013]基于所述相关性检验结果描述所述两个变量时间序列数据之间相关性的最佳Copula函数；
[0014]基于所述最佳拟合分布函数以及最佳Copula函数，构建所述两个变量时间序列数据的联合概率分布函数；
[0015]基于所述联合概率分布函数的联合超越概率的分位数确定所述水文气象灾害二维复合事件的重现期。
[0016]可选的，所述对每个变量时间序列数据自身的时变性以及所述两个变量时间序列数据之间的相关性进行检验，生成相关性检验结果包括：
[0017]结合Kendall
’
s tau相关系数和Spearman相关系数以及滚动窗口方法，对每个变量时间序列数据自身时变性以及所述两个变量时间序列数据之间的相关性进行检验，生成所述相关性检验结果；和/或，
[0018]所述基于所述相关性检验结果描述所述两个变量时间序列数据之间相关性的最佳Copula函数包括：
[0019]若所述相关性检验结果为存在显著趋势，则假设一个时变参数，利用赤池信息准则和贝叶斯信息准则选取所述最佳Copula函数；
[0020]若所述相关性检验结果为不存在显著趋势，则假设一个常数参数，利用赤池信息准则和贝叶斯信息准则选取所述最佳Copula函数。
[0021]可选的，所述对每个变量时间序列数据随时间变化以及所述两个变量时间序列数据之间的变化趋势进行量化，生成趋势量化结果包括：
[0022]利用Mann
‑
Kendall趋势检验法对显著性水平为0.05的每个变量时间序列数据的单调趋势进行描述，生成所述趋势量化结果；和/或，
[0023]所述基于所述趋势量化结果构建一致性和非一致性计算方法，确定每个变量时间序列数据的最佳拟合分布函数包括：
[0024]若所述趋势量化结果为不存在显著趋势，则假设概率分布模型的边缘分布参数不变，利用每个概率分布模型对每个变量时间序列数据进行拟合，得到每个变量时间序列数据的边缘分布函数，利用赤池信息准则和贝叶斯信息准则对每个变量时间序列数据的边缘分布函数进行评估，确定每个变量时间序列数据的最佳拟合分布函数；
[0025]若所述趋势量化结果为存在显著趋势，则创建多个不同非一致性假设的非平稳模型，利用通过最大似然估计方法和贝叶斯马尔可夫链蒙特卡罗方法拟合非平稳边缘分布以及估计参数，利用赤池信息准则和贝叶斯信息准则对每个变量时间序列数据的边缘分布函数进行评估，确定每个变量时间序列数据的最佳拟合分布函数。
[0026]可选的，设J表示一致性假设下的所述两个变量时间序列数据的联合概率分布函数，C表示由θ
C
参数化的Copula函数，所述联合概率分布函数通过以下公式计算：
[0027]J(x1,x2|θ
c
)＝C(F1(x1|θ1),F2(x2|θ2)|θ
C
)＝C(u,v|θ
C
)；
[0028]式中，F1和F2表示两个变量时间序列数据x1和x2的最佳拟合边缘分布函数；θ1和θ2分别表示两个变量时间序列数据x1和x2的最佳拟合分布函数的估计参数；θ
C
表示Copula的参数；u和v表示将两个边缘分布均一化处理到(0,1)范围内的所对应的F1和F2的边缘概率。
[0029]可选的，设J
t
表示非一致性假设下的所述两个变量时间序列数据的联合概率分布函数，C表示由具有超参数θ
C0
和θ
C1
的控制的Copula函数，两个最佳拟合边缘分布函数分别由θ
1t
和θ
2t
控制，所述联合概率分布函数通过以下公式计算：
[0030][0031][0032]式中，F
1t
和F
2t
表示两个变量时间序列数据x
1t
和x
2t
的最佳拟合边缘分布函数；θ
1t
和θ
2t
分别表示两个变量时间序列数据x
1t
和x
2t
的最佳拟合分布函数的估计参数；表示Copula的参数；u
t
、v
t
表示转换到均匀分布u[0,1]下的两个非平稳边缘概率分布，t标识参数或变量随时间变化。
[0033]可选的，所述联合超越概率的分位数Q(p)的计算公式为：
[0034]Q(p)＝(x1＝F1‑1(u|θ1),x2本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种水文气象灾害二维复合事件模拟方法，其特征在于，包括：采集水文气象灾害二维复合事件中两个变量时间序列数据；对每个变量时间序列数据随时间变化以及所述两个变量时间序列数据之间的变化趋势进行量化，生成趋势量化结果；基于所述趋势量化结果构建一致性和非一致性计算方法，确定每个变量时间序列数据的最佳拟合分布函数；对每个变量时间序列数据自身的时变性以及所述两个变量时间序列数据之间的相关性进行检验，生成相关性检验结果；基于所述相关性检验结果描述所述两个变量时间序列数据之间相关性的最佳Copula函数；基于所述最佳拟合分布函数以及最佳Copula函数，构建所述两个变量时间序列数据的联合概率分布函数；基于所述联合概率分布函数的联合超越概率的分位数确定所述水文气象灾害二维复合事件的重现期。2.根据权利要求1所述的方法，其特征在于，所述对每个变量时间序列数据自身的时变性以及所述两个变量时间序列数据之间的相关性进行检验，生成相关性检验结果包括：结合Kendall
’
s tau相关系数和Spearman相关系数以及滚动窗口方法，对每个变量时间序列数据自身时变性以及所述两个变量时间序列数据之间的相关性进行检验，生成所述相关性检验结果；和/或，所述基于所述相关性检验结果描述所述两个变量时间序列数据之间相关性的最佳Copula函数包括：若所述相关性检验结果为存在显著趋势，则假设一个时变参数，利用赤池信息准则和贝叶斯信息准则选取所述最佳Copula函数；若所述相关性检验结果为不存在显著趋势，则假设一个常数参数，利用赤池信息准则和贝叶斯信息准则选取所述最佳Copula函数。3.根据权利要求1所述的方法，其特征在于，所述对每个变量时间序列数据随时间变化以及所述两个变量时间序列数据之间的变化趋势进行量化，生成趋势量化结果包括：利用Mann
‑
Kendall趋势检验法对显著性水平为0.05的每个变量时间序列数据的单调趋势进行描述，生成所述趋势量化结果；和/或，所述基于所述趋势量化结果构建一致性和非一致性计算方法，确定每个变量时间序列数据的最佳拟合分布函数包括：若所述趋势量化结果为不存在显著趋势，则假设概率分布模型的边缘分布参数不变，利用每个概率分布模型对每个变量时间序列数据进行拟合，得到每个变量时间序列数据的边缘分布函数，利用赤池信息准则和贝叶斯信息准则对每个变量时间序列数据的边缘分布函数进行评估，确定每个变量时间序列数据的最佳拟合分布函数；若所述趋势量化结果为存在显著趋势，则创建多个不同非一致性假设的非平稳模型，利用通过最大似然估计方法和贝叶斯马尔可夫链蒙特卡罗方法拟合非平稳边缘分布以及估计参数，利用赤池信息准则和贝叶斯信息准则对每个变量时间序列数据的边缘分布函数进行评估，确定每个变量时间序列数据的最佳拟合分布函数。
4.根据权利要求1至3任一项所述的方法，其特征在于，设J表示一致性假设下的所述两个变量时间序列数据的联合概率分布函数，C表示由θ
c
参数化的Copula函数，所述联合概率分布函数通过以下公式计算：J(x1,x2|θ
c
)＝C(F1(x1|θ1),F2(x2|θ2)|θ
C<...

【专利技术属性】
技术研发人员：王晗，李青，卢文静，马强，张顺福，李辉，孙涛，刘荣华，何秉顺，孙东亚，
申请(专利权)人：中国水利水电科学研究院，
类型：发明
国别省市：