CN109344999A - 一种径流概率预报方法 - Google Patents
一种径流概率预报方法 Download PDFInfo
- Publication number
- CN109344999A CN109344999A CN201811045131.6A CN201811045131A CN109344999A CN 109344999 A CN109344999 A CN 109344999A CN 201811045131 A CN201811045131 A CN 201811045131A CN 109344999 A CN109344999 A CN 109344999A
- Authority
- CN
- China
- Prior art keywords
- probability
- indicate
- model
- runoff
- observation
- 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 44
- 239000011159 matrix material Substances 0.000 claims abstract description 25
- 238000012549 training Methods 0.000 claims abstract description 14
- 239000000203 mixture Substances 0.000 claims abstract description 12
- 230000007704 transition Effects 0.000 claims abstract description 10
- 238000005315 distribution function Methods 0.000 claims abstract description 9
- 230000000977 initiatory effect Effects 0.000 claims abstract description 4
- 238000012360 testing method Methods 0.000 claims description 10
- 238000004364 calculation method Methods 0.000 claims description 9
- 238000007476 Maximum Likelihood Methods 0.000 claims description 4
- 241000208340 Araliaceae Species 0.000 claims description 3
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 claims description 3
- 235000003140 Panax quinquefolius Nutrition 0.000 claims description 3
- 235000008434 ginseng Nutrition 0.000 claims description 3
- 230000005484 gravity Effects 0.000 claims description 3
- 230000015572 biosynthetic process Effects 0.000 claims 1
- 239000004744 fabric Substances 0.000 claims 1
- 238000003786 synthesis reaction Methods 0.000 claims 1
- 230000008569 process Effects 0.000 description 11
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 5
- 238000012952 Resampling Methods 0.000 description 4
- 238000011156 evaluation Methods 0.000 description 3
- 235000013399 edible fruits Nutrition 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 239000002689 soil Substances 0.000 description 2
- PEDCQBHIVMGVHV-UHFFFAOYSA-N Glycerine Chemical compound OCC(O)CO PEDCQBHIVMGVHV-UHFFFAOYSA-N 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000005611 electricity Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000010355 oscillation Effects 0.000 description 1
- 238000010248 power generation Methods 0.000 description 1
- 238000001556 precipitation Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000005514 two-phase flow Effects 0.000 description 1
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
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/23—Clustering techniques
- G06F18/232—Non-hierarchical techniques
- G06F18/2321—Non-hierarchical techniques using statistics or function optimisation, e.g. modelling of probability density functions
- G06F18/23213—Non-hierarchical techniques using statistics or function optimisation, e.g. modelling of probability density functions with fixed number of clusters, e.g. K-means clustering
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/29—Graphical models, e.g. Bayesian networks
- G06F18/295—Markov models or related models, e.g. semi-Markov models; Markov random fields; Networks embedding Markov models
-
- 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
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Landscapes
- Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- Theoretical Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Tourism & Hospitality (AREA)
- Strategic Management (AREA)
- Human Resources & Organizations (AREA)
- Economics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Marketing (AREA)
- Evolutionary Computation (AREA)
- Evolutionary Biology (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Development Economics (AREA)
- Bioinformatics & Cheminformatics (AREA)
- General Business, Economics & Management (AREA)
- Artificial Intelligence (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Engineering & Computer Science (AREA)
- Entrepreneurship & Innovation (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Probability & Statistics with Applications (AREA)
- Game Theory and Decision Science (AREA)
- Educational Administration (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开一种径流概率预报方法,其中,该方法主要包括:采用基于K‑medoids的聚类方法对训练集进行聚类,得到隐含马尔科夫模型HMM的初始化参数;利用Baum–Welch算法对HMM进行学习,得到HMM的状态转移概率矩阵以及观测模型的概率分布;根据贝叶斯信息准则BIC进行模型选择,选择适宜该训练集的HMM状态个数;最终根据给出的预报因子,结合高斯混合回归GMR推理得到条件概率分布函数作为径流概率预报。本发明的概率预报方法引入了径流隐含状态的概念,可利用水文、地形、气象等诸多因素训练得到隐含状态转移概率矩阵,得到有效可靠的未来径流概率预报分布,为水库优化调度决策提供科学依据。
Description
技术领域
本发明涉及水文水资源领域,更具体地,涉及一种径流概率预报方法。
背景技术
径流预报能够在防洪、供水、发电等水库调度运行管理的各个方面发挥重要作用。高精度并且可靠的预报能为流域梯级水库优化调度决策提供科学依据,对流域水资源合理利用具有重大意义。然而,降雨径流的形成过程受到水文、地形、气象等诸多自然因素的影响,呈现高度的非线性、随机性和不确定性特征,给径流预测带来了困难,难以采用传统的方法开展径流预测。因此,如何构建考虑诸多影响径流过程的隐含因素的预测模型是亟需解决的理论和实际工程问题。
隐含马尔科夫模型(hidden Markov model,HMM)由隐含状态序列以及其对应的观测模型构成,其中隐含状态序列是一个典型的马尔科夫链,具有无后效性,对应的观测模型满足独立输出假设:观测模型跟当前隐含状态相关且仅跟当前隐含状态相关。可将隐含马尔科夫模型运用于水文的径流预报当中,视前期径流、土壤、气象等因素为隐含状态序列,对应的径流过程为观测概率模型,能够充分模拟具有非线性、随机性和不确定性的历史径流过程。如何通过训练后的隐含马尔科夫模型来预测未来的径流序列是目前所需解决的问题。
发明内容
针对现有技术的缺陷,本发明的目的在于解决现有降雨径流的形成过程受到水文、地形、气象等诸多自然因素的影响,呈现高度的非线性、随机性和不确定性特征,给径流预测带来了困难,难以采用传统的方法开展径流预测,以及如何通过训练后的隐含马尔科夫模型来预测未来的径流序列的技术问题。
为实现上述目的,本发明提供一种径流概率预报方法,包括如下步骤:
(1)构建由预报径流和其预报因子组成的数据集,对数据进行归一化处理得到训练集;
(2)基于K-medoids的聚类方法对训练集进行聚类,根据聚类后得到的各类的均值向量、协方差矩阵以及径流状态转移概率矩阵作为隐含马尔科夫模型的初始化参数,所述隐含马尔科夫模型包括隐含状态序列及其对应的观测模型;
(3)利用Baum–Welch算法对隐含马尔科夫模型进行学习,得到训练后的隐含马尔科夫模型径流状态转移概率矩阵以及观测模型的概率分布;
(4)设置不同的隐含状态个数K=1~15;重复步骤(2)-(3),并计算相应的贝叶斯信息准则值;
(5)选择贝叶斯信息准则值最小值所对应的K值,作为隐含马尔科夫模型隐含个数;
(6)根据训练得到的隐含马尔科夫模型,结合高斯混合回归(Gaussian MixtureRegression,GMR)方法,利用观测模型的概率分布以及测试集所给的预报因子得到各径流隐含状态下径流条件概率分布,利用径流状态转移概率矩阵求解测试集中每个时间节点落在K个径流隐含状态下的概率,综合推导得到测试集的条件概率分布函数,作为对测试集的径流概率预报结果。
可选地,步骤(2)中K-medoids算法使用核函数代替欧几里得距离来度量两个数据点之间的相似性;
从训练集包括的数据点中随机选择K个中心点,当更新中心点时,每个数据点被分类为与中心点最相似的类别,并且计算每个数据点与其相同类别所有数据点的相似度,选择相似度最高的数据点作为此类别的新的中心点。
可选地,步骤(3)包括:
E步骤:根据现有的隐含马尔科夫模型,计算各个观测数据输入到模型中的计算结果,确定完全数据的对数似然函数:
其中,Q(θ,θold)表示对数似然函数,θ表示隐含马尔科夫模型的所有参数集θ={π,A,φ},θold表示更新前的参数值,t表示数据集的索引,1≤t≤T,T表示数据集个数,πk表示隐含马尔科夫模型在k状态的初始概率,1≤k≤K,K表示隐含状态个数,Aij表示从状态i转移到状态j的概率,φk表示第k个观测概率模型的参数,p(xt|φk)表示第t个观测数据点在第k个观测概率模型下的条件概率值;
γt(k)表示观测数据在时刻t处于第k个状态变量的概率,其概率计算公式如下:
其中,X表示所有观测数据,zt表示t时刻的状态;
ξt(i,j)表示观测数据在时刻t-1处于第i个状态且在时刻t处于第j个状态的概率,其概率计算公式如下:
M步骤:根据以上E步骤得到的参状态参数γt(k)和ξt(i,j),对隐含马尔科夫模型的参数θ={π,A,φ}进行重估计,假设隐含马尔科夫模型中观测概率模型为高斯分布,则φk={μk,∑k},新的参数估计值由下式给出:
可选地,步骤(4)中贝叶斯信息准则值通过以下公式计算:
其中,表示模型的极大似然值,np表示模型参数个数,隐含状态个数不同则对应的模型极大似然数和模型参数个数均不同。
可选地,步骤(6)包括:
在预报的过程中,对于每个隐含状态,观测概率模型的均值向量μk将按照预报因子、预报变量进行拆分:
其中,为表示预报因子均值的子向量,为表示预报变量均值的子向量,与此同时,协方差矩阵也被拆分:
其中,表示预报因子之间的协方差矩阵,表示预报因子与预报变量的协方差矩阵,表示预报变量与预报因子的协方差矩阵,表示预报变量之间的协方差矩阵;
根据联合正态分布的性质推理,给定预报因子yt后,预报值的条件概率服从联合正态分布,其均值向量和协方差矩阵如下:
在原始的高斯混合回归框架中,不同高斯观测模型的权重hk表示为每个观测值属于第k个高斯观测模型的比重;
通过递归计算来估计此权重值,从而不仅考虑预报因子,而且还考虑封装在隐含马尔科夫模型中的时间序列信息:
其中,hk(yt)表示隐含马尔科夫模型前向变量,表示高斯分布的条件概率值;
对于隐含马尔科夫模型,具有K个隐含状态对应的高斯观测模型,因此,给定预报因子yt的情况下,预报变量zt的条件概率分布函数p(zt|yt)为:
其中,表示第k个隐含状态下zt的条件概率分布函数。
总体而言,通过本发明所构思的以上技术方案与现有技术相比,具有以下有益效果:
本发明提供的一种基于隐含马尔科夫模型结合高斯混合回归的径流预报方法,将隐含马尔科夫模型运用与水文的径流预报当中,视前期径流、土壤、气象等因素为隐含状态序列,对应的径流过程为观测概率模型,能够充分的模拟具有非线性、随机性和不确定性的历史径流过程。混合高斯回归是一种概率密度估计的方法,其主要作用是通过联合概率密度函数衍生出该联合概率密度函数的条件概率密度函数作为预报概率。研究将此权重值推广到隐含马尔科夫模型中,通过递归计算来估计权重值,不仅考虑预报因子,而且还考虑封装在HMM中的时间序列信息。
附图说明
图1为本发明提供的隐含马尔科夫模型结合高斯回归预报方法整体流程图;
图2为本发明提供的小湾水库来水预报值、90%预报区间与实际径流对比示意图;
图3为本发明提供的小湾水库来水预报概率积分转换(PIT)统计图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。此外,下面所描述的本发明各个实施方式中所涉及到的技术特征只要彼此之间未构成冲突就可以相互组合。
混合高斯回归是一种概率密度估计的方法,其主要作用是通过联合概率密度函数衍生出该联合概率密度函数的条件概率密度函数作为预报概率,一般运用于混合模型的预测中,而隐含马尔科夫模型可以视为一种条件混合模型,因此可将混合高斯回归拓展到隐含马尔科夫模型作为预测方法。
本发明提供一种基于隐含马尔科夫模型结合高斯混合回归的径流预报方法,可以有效的考虑影响径流过程的各种因素,具有较高的预报精度以及预报可靠性。
图1所示为本发明基于隐含马尔科夫模型结合高斯混合回归的径流预报方法总体流程图,具体包括以下步骤:
(1)构建由预报径流y和其预报因子z组成的数据集x=[y,z],对数据进行归一化处理。
(2)采用基于K-medoids的聚类方法对训练集进行聚类,与传统的K-means方法不同,K-medoids算法使用核函数k(xi,xi′)代替欧几里得距离||xi-x′i||2用于度量两个数据点之间的相似性,本发明采用的是径向基(RBF)核函数:
其中,xi表示第i个数据,x′i表示第i'个数据,d(i,i′)表示第i个数据与第i'个数据的相似性也即核函数值,σ表示RBF的参数。
首先,从原始数据点中随机选择K个中心点m1:K from{1,…,N},当更新中心点时,每个点被分类为与中心点最相似的类别,并且计算它与相同类别所有点的相似度,选择相似度最高(与其他所有点的核函数值之和最大)的点作为此类别的新的中心点mk:
其中,z′i表示第i'数据的类别。
根据聚类后的各类的均值向量μk、协方差矩阵∑k以及状态转移概率矩阵A得到隐含马尔科夫模型的初始化参数。
(3)利用Baum–Welch算法对隐含马尔科夫模型进行学习,得到隐含马尔科夫模型的径流状态转移概率矩阵A以及观测模型的概率分布p(xt|st=k,φk);
E过程:根据现有的模型,计算各个观测数据输入到模型中的计算结果,首先确定完全数据的对数似然函数:
γt(k)表示给定模型和观测数据,在时刻t处于第k个状态变量的概率,其概率计算公式如下:
ξt(i,j)表示给定模型和观测数据,在时刻t-1处于第i个状态且在时刻t处于第j个状态的概率,其概率计算公式如下:
M过程:M过程是重新估计参数。
对于观测概率模型为高斯分布的隐含马尔科夫模型来说,为了得到每个状态的均值μk和Σk项的M步,新的参数估计值由下式给出:
(4)设置不同的隐含状态个数K=1~15;重复步骤(2)-(3),并计算相应的贝叶斯信息准则值BIC:
(5)选择贝叶斯信息准则值最小值所对应的K值,作为径流预测模型的隐含个数;
(6)根据训练得到的隐含马尔科夫模型,结合高斯混合回归方法,推导得到测试集的条件概率分布函数,作为概率预报结果。
隐含马尔科夫模型也称为非独立混合模型,它可以解释为一个混合模型的扩展,其中每个观测的混合分布概率不是独立的,而是取决于先前观察的潜在状态变量。在预报的过程中,隐含马尔科夫模型的所有变量被区分为预报因子(前期径流和气象指标)和预报变量(未来径流)xt=[yt,zt]。相应的对于每个隐含状态,观测概率模型的均值向量将按照预报因子、预报变量进行拆分:
上式中,为表示预报因子均值的子向量,为表示预报变量均值的子向量,与此同时,协方差矩阵也被拆分:
根据联合正态分布的性质推理,给定预报因子yt后,预报值的条件概率服从联合正态分布,其均值向量和协方差矩阵如下:
在原始的高斯混合回归框架中,不同高斯观测模型的权重hk表示,为每个观测值属于第k个高斯观测模型的比重,此权重值并未考虑观测序列的顺序信息。研究将此权重值推广到隐含马尔科夫模型中,通过递归计算来估计此权重值,从而不仅考虑预报因子,而且还考虑封装在HMM中的时间序列信息:
对于隐含马尔科夫模型,具有K个隐含状态对应的高斯观测模型。因此,给定预报因子yt的情况下,预报变量zt的条件概率分布函数p(zt|yt)为:
本发明以澜沧江下游小湾电的月平均径流量时间序列为对象,数据采用1953年到2010年月径流数据。以前1月径流、前12月径流以及南方涛动指数和亚洲经向环流指数为预报因子,按照图1所示的基于隐含马尔科夫模型结合高斯混合回归的径流预报方法进行预报。
表1给出小湾电站入库径流概率预报的评分表,其中,表1中SSCRPS表示连续分级概率评分,SSMSE表示均方误差评分,SSRMSE表示均方根误差评分,SSMAE表示平均绝对误差评分。
表1概率预报评分表
所有评分SS均与历史重采样预报对比,若预报精度高于历史重采样预报则为正值,若预报精度低于历史重采样预报则为负值,评分范围在-1~1之间,评分越高代表预报精度越好。其中SSCRPS为连续分级概率评分,其同属考虑了预报的偏差以及概率预报的范围,是一种概率预报常用的评价指标;SSMSE为均方误差评分,其评分公式与水文预报常用的评价指标“纳什效率系数”相同;SSRMSE为均方根误差评分即Brier评分,再概率预报中广泛使用;SSMAE为平均绝对误差评分。从表中可以看出,隐含马尔科夫模型的预报评分均大于0,代表研究所提出的概率预报精度均高于历史重采样预报。从表中看出:其中12月~4月的评分普遍比5月~11月的评分要高。这一结果表明,月预报精度与降雨的季节性密切相关:枯水期径流主要取决于流域“干湿”情况,以前期径流作为预报因子,能够较好的进行预报;相比而言,汛期径流主要取决于降水情况,由于降水可预报性很低,汛期径流预报精度也不高。
图2展示了本预报方法在测试期间的预报值、90%预报区间与实际月径流的对比图,图中可以看出本发明的预报方法具有较高的精度。
在预报精度评价的基础上,进一步对预报可靠性进行分析。分析观测值相对于预报概率积分转换(PIT)值,如图3所示。可以看到,概率积分转换整体上呈现均匀分布的特征。这一结果说明:观测值基本上可以看作来自对应概率预报的随机采样。这一结果表明,此预报具有较高的可靠性。
本领域的技术人员容易理解,以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。
Claims (5)
1.一种径流概率预报方法,其特征在于,包括如下步骤:
(1)构建由预报径流和其预报因子组成的数据集,对数据进行归一化处理得到训练集;
(2)基于K-medoids的聚类方法对训练集进行聚类,根据聚类后得到的各类的均值向量、协方差矩阵以及径流状态转移概率矩阵作为隐含马尔科夫模型的初始化参数,所述隐含马尔科夫模型包括隐含状态序列及其对应的观测模型;
(3)利用Baum–Welch算法对隐含马尔科夫模型进行学习,得到训练后的隐含马尔科夫模型径流状态转移概率矩阵以及观测模型的概率分布;
(4)设置不同的隐含状态个数K=1~15;重复步骤(2)-(3),并计算相应的贝叶斯信息准则值;
(5)选择贝叶斯信息准则值最小值所对应的K值,作为隐含马尔科夫模型隐含个数;
(6)根据训练得到的隐含马尔科夫模型,结合高斯混合回归方法,利用观测模型的概率分布以及测试集所给的预报因子得到各径流隐含状态下径流条件概率分布,利用径流状态转移概率矩阵求解测试集中每个时间节点落在K个径流隐含状态下的概率,综合推导得到测试集的条件概率分布函数,作为对测试集的径流概率预报结果。
2.根据权利要求1所述的径流概率预报方法,其特征在于,所述步骤(2)中K-medoids算法使用核函数代替欧几里得距离来度量两个数据点之间的相似性;
从训练集包括的数据点中随机选择K个中心点,当更新中心点时,每个数据点被分类为与中心点最相似的类别,并且计算每个数据点与其相同类别所有数据点的相似度,选择相似度最高的数据点作为此类别的新的中心点。
3.根据权利要求1所述的径流概率预报方法,其特征在于,所述步骤(3)包括:
E步骤:根据现有的隐含马尔科夫模型,计算各个观测数据输入到模型中的计算结果,确定完全数据的对数似然函数:
其中,Q(θ,θold)表示对数似然函数,θ表示隐含马尔科夫模型的所有参数集θ={π,A,φ},θold表示更新前的参数值,t表示数据集的索引,1≤t≤T,T表示数据集个数,πk表示隐含马尔科夫模型在k状态的初始概率,1≤k≤K,K表示隐含状态个数,Aij表示从状态i转移到状态j的概率,φk表示第k个观测概率模型的参数,p(xt|φk)表示第t个观测数据点在第k个观测概率模型下的条件概率值;
γt(k)表示观测数据在时刻t处于第k个状态变量的概率,其概率计算公式如下:
其中,X表示所有观测数据,zt表示t时刻的状态;
ξt(i,j)表示观测数据在时刻t-1处于第i个状态且在时刻t处于第j个状态的概率,其概率计算公式如下:
M步骤:根据以上E步骤得到的参状态参数γt(k)和ξt(i,j),对隐含马尔科夫模型的参数θ={π,A,φ}进行重估计,假设隐含马尔科夫模型中观测概率模型为高斯分布,则φk={μk,∑k},新的参数估计值由下式给出:
4.根据权利要求3所述的径流概率预报方法,其特征在于,所述步骤(4)中贝叶斯信息准则值通过以下公式计算:
其中,表示模型的极大似然值,np表示模型参数个数,隐含状态个数不同则对应的模型极大似然数和模型参数个数均不同。
5.根据权利要求1所述的径流概率预报方法,其特征在于,所述步骤(6)包括:
在预报的过程中,对于每个隐含状态,观测概率模型的均值向量μk将按照预报因子、预报变量进行拆分:
其中,为表示预报因子均值的子向量,为表示预报变量均值的子向量,与此同时,协方差矩阵也被拆分:
其中,表示预报因子之间的协方差矩阵,表示预报因子与预报变量的协方差矩阵,表示预报变量与预报因子的协方差矩阵,表示预报变量之间的协方差矩阵;
根据联合正态分布的性质推理,给定预报因子yt后,预报值的条件概率服从联合正态分布,其均值向量和协方差矩阵如下:
在原始的高斯混合回归框架中,不同高斯观测模型的权重hk表示为每个观测值属于第k个高斯观测模型的比重;
通过递归计算来估计此权重值,从而不仅考虑预报因子,而且还考虑封装在隐含马尔科夫模型中的时间序列信息:
其中,hk(yt)表示隐含马尔科夫模型前向变量,表示高斯分布的条件概率值;
对于隐含马尔科夫模型,具有K个隐含状态对应的高斯观测模型,因此,给定预报因子yt的情况下,预报变量zt的条件概率分布函数p(zt|yt)为:
其中,表示第k个隐含状态下zt的条件概率分布函数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811045131.6A CN109344999A (zh) | 2018-09-07 | 2018-09-07 | 一种径流概率预报方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811045131.6A CN109344999A (zh) | 2018-09-07 | 2018-09-07 | 一种径流概率预报方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN109344999A true CN109344999A (zh) | 2019-02-15 |
Family
ID=65305003
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811045131.6A Pending CN109344999A (zh) | 2018-09-07 | 2018-09-07 | 一种径流概率预报方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109344999A (zh) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110910004A (zh) * | 2019-11-18 | 2020-03-24 | 国电湖南巫水水电开发有限公司 | 一种多重不确定性的水库调度规则提取方法及*** |
CN110909943A (zh) * | 2019-11-27 | 2020-03-24 | 淮阴工学院 | 一种多尺度多因子联合驱动的月径流概率预报方法 |
CN110969290A (zh) * | 2019-11-18 | 2020-04-07 | 国电湖南巫水水电开发有限公司 | 一种基于深度学习的径流概率预测方法及*** |
CN113191561A (zh) * | 2021-05-11 | 2021-07-30 | 华中科技大学 | 一种基于高斯混合模型的径流随机模拟方法及*** |
CN113837475A (zh) * | 2021-09-27 | 2021-12-24 | 中水珠江规划勘测设计有限公司 | 有向图深度神经网络径流概率预报方法、***、设备及终端 |
CN115952577A (zh) * | 2022-12-06 | 2023-04-11 | 中国水利水电科学研究院 | 一种梯级水库群溃决风险分析方法 |
CN116720753A (zh) * | 2023-08-07 | 2023-09-08 | 浙江鹏信信息科技股份有限公司 | 水文数据的处理方法、***及可读存储介质 |
CN116960988A (zh) * | 2023-09-12 | 2023-10-27 | 浙江贵仁信息科技股份有限公司 | 一种水库优化调度方法、***、设备以及介质 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106782510A (zh) * | 2016-12-19 | 2017-05-31 | 苏州金峰物联网技术有限公司 | 基于连续混合高斯hmm模型的地名语音信号识别方法 |
CN107705558A (zh) * | 2017-09-15 | 2018-02-16 | 东南大学 | 一种多数据类型的贝叶斯od矩阵估计方法 |
CN107817106A (zh) * | 2017-10-11 | 2018-03-20 | 温州大学 | 基于贝叶斯残余变换‑奇异值分解‑高斯混合隐马尔科夫模型框架的滚动轴承故障诊断方法 |
-
2018
- 2018-09-07 CN CN201811045131.6A patent/CN109344999A/zh active Pending
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106782510A (zh) * | 2016-12-19 | 2017-05-31 | 苏州金峰物联网技术有限公司 | 基于连续混合高斯hmm模型的地名语音信号识别方法 |
CN107705558A (zh) * | 2017-09-15 | 2018-02-16 | 东南大学 | 一种多数据类型的贝叶斯od矩阵估计方法 |
CN107817106A (zh) * | 2017-10-11 | 2018-03-20 | 温州大学 | 基于贝叶斯残余变换‑奇异值分解‑高斯混合隐马尔科夫模型框架的滚动轴承故障诊断方法 |
Non-Patent Citations (1)
Title |
---|
YONGQI LIU等: "Monthly streamflow forecasting based on hidden Markov model and Gaussian Mixture Regression", 《JOURNAL OF HYDROLOGY》 * |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110910004A (zh) * | 2019-11-18 | 2020-03-24 | 国电湖南巫水水电开发有限公司 | 一种多重不确定性的水库调度规则提取方法及*** |
CN110969290A (zh) * | 2019-11-18 | 2020-04-07 | 国电湖南巫水水电开发有限公司 | 一种基于深度学习的径流概率预测方法及*** |
CN110969290B (zh) * | 2019-11-18 | 2023-09-26 | 国电湖南巫水水电开发有限公司 | 一种基于深度学习的径流概率预测方法及*** |
CN110909943A (zh) * | 2019-11-27 | 2020-03-24 | 淮阴工学院 | 一种多尺度多因子联合驱动的月径流概率预报方法 |
CN113191561A (zh) * | 2021-05-11 | 2021-07-30 | 华中科技大学 | 一种基于高斯混合模型的径流随机模拟方法及*** |
CN113837475A (zh) * | 2021-09-27 | 2021-12-24 | 中水珠江规划勘测设计有限公司 | 有向图深度神经网络径流概率预报方法、***、设备及终端 |
CN113837475B (zh) * | 2021-09-27 | 2024-04-05 | 中水珠江规划勘测设计有限公司 | 有向图深度神经网络径流概率预报方法、***、设备及终端 |
CN115952577A (zh) * | 2022-12-06 | 2023-04-11 | 中国水利水电科学研究院 | 一种梯级水库群溃决风险分析方法 |
CN116720753A (zh) * | 2023-08-07 | 2023-09-08 | 浙江鹏信信息科技股份有限公司 | 水文数据的处理方法、***及可读存储介质 |
CN116720753B (zh) * | 2023-08-07 | 2023-10-31 | 浙江鹏信信息科技股份有限公司 | 水文数据的处理方法、***及可读存储介质 |
CN116960988A (zh) * | 2023-09-12 | 2023-10-27 | 浙江贵仁信息科技股份有限公司 | 一种水库优化调度方法、***、设备以及介质 |
CN116960988B (zh) * | 2023-09-12 | 2023-12-15 | 浙江贵仁信息科技股份有限公司 | 一种水库优化调度方法、***、设备以及介质 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109344999A (zh) | 一种径流概率预报方法 | |
Cannon | A flexible nonlinear modelling framework for nonstationary generalized extreme value analysis in hydroclimatology | |
CN106650767B (zh) | 基于聚类分析和实时校正的洪水预报方法 | |
CN107766990B (zh) | 一种光伏电站发电功率的预测方法 | |
CN105354646B (zh) | 一种混合粒子群优化与极限学习机的电力负荷预测方法 | |
CN109840587B (zh) | 基于深度学习的水库入库流量预测方法 | |
CN105701572B (zh) | 一种基于改进高斯过程回归的光伏短期出力预测方法 | |
CN110363349B (zh) | 一种基于ascs的lstm神经网络水文预测方法及*** | |
CN109583565A (zh) | 基于注意力模型长短时记忆网络的洪水预测方法 | |
CN113554466A (zh) | 一种短期用电量预测模型构建方法、预测方法和装置 | |
CN112733996A (zh) | 基于GA-PSO优化XGBoost的水文时间序列预测方法 | |
CN109299812A (zh) | 一种基于深度学习模型和knn实时校正的洪水预测方法 | |
CN111260126A (zh) | 计及天气与气象因子关联度的短期光伏发电预测方法 | |
CN115130741A (zh) | 基于多模型融合的多因素电力需求中短期预测方法 | |
Malakouti et al. | Evaluation of the application of computational model machine learning methods to simulate wind speed in predicting the production capacity of the Swiss basel wind farm | |
Akpinar et al. | Forecasting natural gas consumption with hybrid neural networks—Artificial bee colony | |
CN116435998A (zh) | 一种光伏发电功率的预测方法 | |
CN111860974A (zh) | 一种基于状态空间和联合分布的干旱多级预测方法 | |
CN114219126B (zh) | 基于残差修正的小水电渗透地区网供负荷预测方法 | |
CN114372631A (zh) | 一种基于小样本学习和lstm的缺资料地区径流预测方法 | |
CN116245259B (zh) | 基于深度特征选择的光伏发电预测方法、装置和电子设备 | |
Guo et al. | Wind speed forecasting of genetic neural model based on rough set theory | |
CN116663393A (zh) | 一种基于随机森林的配电网持续高温下故障风险等级预测方法 | |
Ma et al. | Research and application of second-hand housing price prediction model based on LSTM | |
CN115345245A (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 | ||
WD01 | Invention patent application deemed withdrawn after publication | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20190215 |