CN114707689A - 一种参数化的水文气象极值序列的趋势检测方法 - Google Patents
一种参数化的水文气象极值序列的趋势检测方法 Download PDFInfo
- Publication number
- CN114707689A CN114707689A CN202210047362.0A CN202210047362A CN114707689A CN 114707689 A CN114707689 A CN 114707689A CN 202210047362 A CN202210047362 A CN 202210047362A CN 114707689 A CN114707689 A CN 114707689A
- Authority
- CN
- China
- Prior art keywords
- trend
- sequence
- extreme value
- distribution
- parameterized
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 25
- 238000001514 detection method Methods 0.000 claims abstract description 48
- 238000005070 sampling Methods 0.000 claims abstract description 8
- 238000004422 calculation algorithm Methods 0.000 abstract description 6
- 230000000694 effects Effects 0.000 description 6
- 238000012360 testing method Methods 0.000 description 6
- 238000010586 diagram Methods 0.000 description 5
- 230000035945 sensitivity Effects 0.000 description 3
- 238000010998 test method Methods 0.000 description 3
- 238000010792 warming Methods 0.000 description 3
- 230000001174 ascending effect Effects 0.000 description 2
- 230000002146 bilateral effect Effects 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 238000010206 sensitivity analysis Methods 0.000 description 2
- 238000013459 approach Methods 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 238000005315 distribution function Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 239000005431 greenhouse gas Substances 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 230000003121 nonmonotonic effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0635—Risk analysis of enterprise or organisation activities
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/10—Services
- G06Q50/26—Government or public services
Landscapes
- Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- Engineering & Computer Science (AREA)
- Strategic Management (AREA)
- Economics (AREA)
- Tourism & Hospitality (AREA)
- General Business, Economics & Management (AREA)
- Entrepreneurship & Innovation (AREA)
- Marketing (AREA)
- Theoretical Computer Science (AREA)
- Development Economics (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Educational Administration (AREA)
- Game Theory and Decision Science (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- General Health & Medical Sciences (AREA)
- Complex Calculations (AREA)
Abstract
本发明公开了一种参数化的水文气象极值序列的趋势检测方法,包括以下步骤:(1)采用广义极值分布GEV函数建立水文气象极值序列的趋势模型;(2)基于贝叶斯推断框架估计水文气象极值序列的趋势模型的参数;(3)根据趋势模型的参数估计值获得贝叶斯对数似然比,然后得出水文气象极值序列的趋势。本发明从参数化角度有效保证水文气象极值序列的趋势检测精度;基于贝叶斯统计推断框架的使用一方面充分保证复杂趋势模式情形下对于水文气象极值序列的分布参数估计精度不会因为引入过多的参数而受到影响,另一方面采用Gibbs采样算法简化了水文气象极值序列分布参数的估计过程。
Description
技术领域
本发明涉及一种趋势检测方法,尤其涉及一种水文气象极值序列的趋势检测方法。
背景技术
全球变暖是对地球的主要威胁,气候变暖是人为温室气体排放增加的结果,正在逐步改变地球气候模式。根据政府间气候变化专门委员会IPCC的说法,温度升高导致的水文水循环特征急剧变化。因此,全球变暖将影响水文过程,并导致世界不同地区出现水文气象极端事件的风险增加,为了有效预测和防范水文气象极端事件发生风险,最为关键一步就是对水文气象极端事件进行有效的趋势性检测和识别。目前国内外使用最多的是非参数法的趋势性检测方法:
Mann-Kendall趋势性检测方法目前应用于水文气象极值序列趋势性检测最为广泛,由于其仅仅适用于单调性趋势的检测,对于非线性,如先增加后降低的非单调性趋势趋势的检测效果存在明显不足,且容易出现假设检验中的两类错误:“弃真”和“纳伪”错误。
发明内容
发明目的:本发明旨在提供一种水文气象极值序列的趋势检测方法,解决水文气象极值序列趋势性检测只能单调性趋势监测的问题。。
技术方案:本发明的一种参数化的水文气象极值序列的趋势检测方法,包括以下步骤:
(1)采用广义极值分布GEV函数建立水文气象极值序列的趋势模型;
(2)基于贝叶斯推断框架估计水文气象极值序列的趋势模型的参数;
(3)根据趋势模型的参数估计值获得贝叶斯对数似然比,然后得出水文气象极值序列的趋势。
步骤(1)的趋势模型:
式中,Xt为水文气象数据,GEV分布的位置μt,尺度σt和形状参数kt是随着时间t逐年变化的。
水文气象极值序列的趋势模型GEV分布的位置μt的趋势结构包括:
式中,μ0、μ1、μ2、μ3为未知待定参数。
水文气象极值序列的趋势模型GEV分布的尺度σt的趋势结构包括:
式中,σ0、σ1、σ2、σ3为未知待定参数。
水文气象极值序列的趋势模型GEV分布的形状参数κt的趋势结构包括:
式中,κ0、κ1、κ2、κ3为未知待定参数。
步骤(2)包括以下步骤:
(21)趋势模型的参数先验分布的假定;
(22)基于Gibbs的蒙特卡洛随机采样法获得潜在的后验分布参数;
(23)获得水文气象极值序列的趋势模型参数估计值。
步骤(2)中假设先验分布参数集为β,后验分布按照公式(5)求得:
式中,f(Xt|β,t)各时间步时变GEV参数的后验分布概率密度值。
步骤(3)中获得贝叶斯对数似然比
Bay_LR=2(LogL1-LogL0) 6
式中,LogL0为平稳模型的对数似然值,而LogL1代表了其余趋势性模型的对数似然值。
有益效果:与现有技术相比,本发明具有如下显著优点:
(1)有效克服了非参数趋势检测法对于水文气象极值序列非单调性趋势检测识别效果不佳的问题,从参数化角度有效保证水文气象极值序列的趋势检测精度。
(2)基于贝叶斯统计推断框架的使用一方面充分保证复杂趋势模式情形下对于水文气象极值序列的分布参数估计精度不会因为引入过多的参数而受到影响,另一方面采用Gibbs采样算法简化了水文气象极值序列分布参数的估计过程。
附图说明
图1为本发明中流程示意图;
图2为本发明各站点极值序列的自相关结果;
图3为本发明非参数化Mann-Kendall趋势检测结果;
图4为本发明非参数化Mann-Kendall趋势检测统计变量Z值得空间分布经验密度图;
图5为本发明特定站点的后验分布参数计算过程;
图6为本发明特定站点的后验分布的趋势模型拟合效果图;
图7为本发明基于Bay_LR趋势检测法进行参数化的极值序列检测结果;
图8为本发明参数非单调性结果;
图9为本发明参数化和非参数化对于样本长度的敏感性分析结果。
具体实施方式
下面结合附图对本发明的技术方案作进一步说明。
为便于理解本发明,首先对水文气象极值序列非参数趋势检验法—Mann-Kendall趋势检验法做以下说明:
式中,xj和xk分别为极值序列中第j和第k个观测值;sgn(.)为依据xj与xk相对大小的符号判断函数;同时统计量S的方差可以按照如下公式求得:
式中,p是数据集中观测值相同的组数,tj是第j个数值相同组。那么标准化后的Z值统计变量可以按如下公式计算:
在双边趋势检验中,对于给定的置信水平(显著性水平)α,若|Z|≥Z1-α/,则原假设是不可接受的,即在置信水平α(显著性检验水平)上,时间序列数据存在明显的上升或下降趋势。Z为正值表示上升趋势,负值表示减少趋势,Z的绝对值在大于等于1.645,1.96,2.576时表示分别通过了置信度90%,95%,99%的显著性检验。
如图1所示,本发明的一种参数化的水文气象极值序列的趋势检测方法,包括以下步骤:
步骤1.采用广义极值分布GEV函数建立水文气象极值序列的趋势模型
首先采用极值分析领域常用的广义极值分布函数GEV进行水文气象极值序列的拟合,
式中,Xt为水文气象数据,GEV分布的位置μt,尺度σt和形状参数κt是随着时间t逐年变化的;
位置μt的趋势结构包括:
式中,μ0、μ1、μ2、μ3为未知待定参数。
尺度σt的趋势结构包括:
式中,σ0、σ1、σ2、σ3为未知待定参数。
形状参数κt的趋势结构包括:
式中,k0、k1、k2、κ3为未知待定参数。
水文气象极值序列的趋势模型一共包括63种趋势模式结构和1种平稳即无趋势模式。
步骤2:基于贝叶斯推断框架估计水文气象极值序列的趋势模型的参数;
(1)首先对水文气象极值序列平稳模式的参数进行估计。通过参数先验分布模式假定,基于Gibbs采样算法的后验分布推求进行。
(2)接着对水文气象极值序列63种趋势模式下的参数进行估计。包括如下步骤:
(21)参数先验分布的结构如下:
式中,位置和形状参数的先验分布模式服从均值为0方差为1000的正太分布,而尺度参数服从参数为25的柯西分布。
(22)基于Gibbs采样算法进行后验分布的求解:
(221)为参数集Φt设定初始化值Φ0。
(222)将初始化值代入到先验分布模型中求得先验分布的似然值LRprio,进而基于观测值样本序列获得GEV分布似然值LRGEV,进而得到后验密度的对数似然值为前两者之和。
(223)不断地更新迭代,重复M步,并结合Gibbs采样算法当后验密度对数似然值逐渐趋近定值时,采用算法终止。
(23)基于M个后验分布参数集求得Φt的中值和参数上下界。获得水文气象极值序列的趋势模型参数估计值。
步骤3:基于Bay_LR趋势检测法进行参数化的水文气象极值序列检测
获得了27种情形下的参数估计结果,接下来Bay_LR趋势检测法可以按照如下公式求解贝叶斯对数似然比Bay_LR
Bay_LR=2(LogL1-LogL0)
其中,LogL0为平稳模型的对数似然值而LogL1代表了其余26种趋势性模型的对数似然值;由于以上的Bay_LR统计量收敛于(卡方)分布,其中卡方分布参数Δ为平稳模型和非平稳模型间的差异,为此可以用于检验零趋势假设(分布参数没有趋势)是否存在的重要依据:Bay_LR的值大于的α分位数,上述零假设可在α显著性水平上被拒绝,也就代表所研究的水文气象极值序列检测序列存在趋势性。当α=0.95时,就可以获得95%显著水平的水文气象极值序列趋势性,Bay_LE会返回固定的p-value值,p-vlaue值小于5%表明具有95%显著水平的趋势性。
为了有效验证本发明提出的参数化趋势检测方法的有效性,给出以下实施例。以中国范围内839个气象站点的年极大日降雨序列作为研究对象。以中国范围内的一般气象站点观测到的日降雨量为基础数据,按照年极大值提取法获得年极大日降雨序列,样本长度为1951-2018年一共68年的数据长度。
(1)非参数的Mann-Kendall趋势检测法的识别结果分析
为了对比分析本发明的参数化的水文气象极值序列趋势检测法的效果,计算出了非参数的Mann-Kendall趋势法的检测结果。但由于Mann-Kendall的检验结果很大程度受到序列的自相关的影响,为此需要先分析降雨极值序列的自相关情况,计算结果如图2所示。839个站点中ACF的p-value值都是大于0.05的表明大部分降雨序列自相关性较弱可以采用Mann-Kendall趋势检验法对原始极值序列进行计算分析。
由图3可知,839个气象站点的年极大日雨量序列只有69个站点具有95%显著水平的上升或者下降趋势性,只占到了69/839=8.2%。图4中Z值的空间分布密度函数图也很好地体现了这一点。
(2)基于贝叶斯的参数估计结果
由于站点众多,且趋势模型种类众多,本节只是显示了站点编号为52203(42.5°N,93.3°E)的最优趋势模型的分析过程,如图5和图6所示。图5说明了由于有效的先验分布假定使得后验分布在特定的采样步骤之后能够收敛于定值,图6也很好地体现了最优趋势模型估计出的理论分布模型和经验分布最为接近的。
(3)基于参数化的水文气象极值序列趋势性检验结果
图7显示了最终的水文气象极值序列Bay_LR趋势检测法的优选结果。由图7中结果可知,一共有330个站点体现出单调性或者非单调性趋势,总体来看参数化的趋势性检测方法可以有效检测出330/839=39.5%。且Bay_LR趋势检测法检测出的站点编号大部分都包括了Mann-Kendall趋势检测的站点编号,换句话说,本发明提出Bay_LR趋势检测法不但能够检测出Mann-Kendall检测出的单调性结果,同时也可以有效识别出多余非单调趋势和参数间存在反向趋势的情况,相比Mann-Kendall检测法效果提升了39.5%-8.2%=31.3%。同时由于参数非单调性趋势和参数间反向趋势的存在图8中所示使得Mann-Kendall检测效果不佳,无法识别此类的极值序列趋势性,也是参数化的Bay_LR趋势检测法能够比Mann-Kendall检测法多出31.3%的原因所在。
(4)非参数法和参数法对于样本长度的敏感性分析
图9中显示出了两种水文气象极值序列趋势检测方法对于样本长度的敏感性响应关系。图9中上面4幅图是参数化的趋势检测方法分别对选定的4个站点的降雨极值序列按照样本长度由41变化到60时,步长为1,Bay_LR趋势检测法检测出的趋势性模型,图中LR等于0代表了序列没有趋势性,而1-26分别代表了26种的趋势性模式结构,可见总的来看非线性趋势1比较常见,而三次型的非线性趋势2模式基本上没有符合的站点。而图9下面4幅图分别是Mann-Kendall趋势检测法对于样本长度的敏感性分析结果,可以发现随着样本长度增加,Mann-Kendall的关键统计量Z值是逐步增大的,并且由没有趋势逐步过渡到趋势性的情形,可见其对于样本长度过于敏感,而参数化的Bay_LR趋势检测法则表现出较弱的敏感性,不管是样本长度增大还是减小,序列的趋势性都能够有效得检测出来。
Claims (9)
1.一种参数化的水文气象极值序列的趋势检测方法,其特征在于:包括以下步骤:
(1)采用广义极值分布GEV函数建立水文气象极值序列的趋势模型;
(2)基于贝叶斯推断框架估计水文气象极值序列的趋势模型的参数;
(3)根据趋势模型的参数估计值获得贝叶斯对数似然比,然后得出水文气象极值序列的趋势。
6.根据权利要求1所述的参数化的水文气象极值序列的趋势检测方法,其特征在于:步骤(2)包括以下步骤:
(21)趋势模型的参数先验分布的假定;
(22)基于Gibbs的蒙特卡洛随机采样法获得潜在的后验分布参数;
(23)获得水文气象极值序列的趋势模型参数估计值。
8.根据权利要求1所述的参数化的水文气象极值序列的趋势检测方法,其特征在于:步骤(3)中获得贝叶斯对数似然比
Bay_LR=2(LogL1-LogL0) (6)
其中,LogL0为平稳模型的对数似然值,而LogL1代表了其余趋势性模型的对数似然值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210047362.0A CN114707689A (zh) | 2022-01-17 | 2022-01-17 | 一种参数化的水文气象极值序列的趋势检测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210047362.0A CN114707689A (zh) | 2022-01-17 | 2022-01-17 | 一种参数化的水文气象极值序列的趋势检测方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN114707689A true CN114707689A (zh) | 2022-07-05 |
Family
ID=82167550
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210047362.0A Pending CN114707689A (zh) | 2022-01-17 | 2022-01-17 | 一种参数化的水文气象极值序列的趋势检测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114707689A (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116663392A (zh) * | 2023-04-23 | 2023-08-29 | 中国水利水电科学研究院 | 水文气象灾害二维复合事件模拟方法、装置、设备及介质 |
-
2022
- 2022-01-17 CN CN202210047362.0A patent/CN114707689A/zh active Pending
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116663392A (zh) * | 2023-04-23 | 2023-08-29 | 中国水利水电科学研究院 | 水文气象灾害二维复合事件模拟方法、装置、设备及介质 |
CN116663392B (zh) * | 2023-04-23 | 2024-02-20 | 中国水利水电科学研究院 | 水文气象灾害二维复合事件模拟方法、装置、设备及介质 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107092744B (zh) | 基于emd-svr的地表沉降量预测方法 | |
CN110197020B (zh) | 一种环境变化对水文干旱影响的分析方法 | |
Caires | Extreme value analysis: wave data | |
CN109815627B (zh) | 一种径流变异影响因子确定方法及系统 | |
CN110081923B (zh) | 野外基线环境参数自动采集系统故障检测方法与装置 | |
Chikobvu et al. | Modelling of extreme minimum rainfall using generalised extreme value distribution for Zimbabwe | |
CN111367959B (zh) | 一种零时滞非线性扩展Granger因果分析方法 | |
Lu et al. | Application of profile likelihood function to the uncertainty analysis of hydrometeorological extreme inference | |
CN113484813B (zh) | 多环境应力下的智能电表故障率预估方法及系统 | |
CN114707689A (zh) | 一种参数化的水文气象极值序列的趋势检测方法 | |
CN105205217A (zh) | 一种判别水文时间序列非平稳性的方法 | |
CN111832176B (zh) | 降雨条件下全极化微波辐射计海面风场反演方法及系统 | |
CN113985496B (zh) | 一种基于lstm-gm神经网络模型的风暴潮智能预报方法 | |
CN117310118B (zh) | 一种地下水污染可视化监测方法 | |
CN113642255A (zh) | 基于多尺度卷积循环神经网络的光伏发电功率预测方法 | |
CN102542284A (zh) | 辨识光谱的方法 | |
Whalen et al. | An evaluation of the self-determined probability-weighted moment method for estimating extreme wind speeds | |
CN117194872A (zh) | 一种风电场风速校正方法、系统、计算机设备及存储介质 | |
CN111651427A (zh) | 一种基于gamlss模型的非一致性水文频率计算方法 | |
Petrov et al. | Maximum entropy estimates of extreme significant wave heights from satellite altimeter data | |
CN114691661B (zh) | 一种基于同化的云导风、温湿廓线预处理分析方法及系统 | |
CN114819260A (zh) | 一种水文时间序列预测模型动态生成方法 | |
Junhong et al. | Software reliability nonlinear modeling and its fuzzy evaluation | |
Yang et al. | Modeling stream flow extremes under non-time-stationary conditions | |
CN108052727B (zh) | 一种金属栅功函数变化导致栅电容统计分布的估计方法 |
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 |