CN114564487A - 预报预测相结合的气象栅格数据更新方法 - Google Patents
预报预测相结合的气象栅格数据更新方法 Download PDFInfo
- Publication number
- CN114564487A CN114564487A CN202210464764.0A CN202210464764A CN114564487A CN 114564487 A CN114564487 A CN 114564487A CN 202210464764 A CN202210464764 A CN 202210464764A CN 114564487 A CN114564487 A CN 114564487A
- Authority
- CN
- China
- Prior art keywords
- time
- model
- meteorological
- time sequence
- prediction
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 63
- 238000012545 processing Methods 0.000 claims description 12
- 241001123248 Arma Species 0.000 claims description 7
- 238000004364 calculation method Methods 0.000 claims description 5
- 238000012360 testing method Methods 0.000 claims description 5
- 238000007689 inspection Methods 0.000 claims description 4
- 238000000342 Monte Carlo simulation Methods 0.000 claims description 3
- 239000011159 matrix material Substances 0.000 claims description 3
- 241001611138 Isma Species 0.000 claims description 2
- 238000001514 detection method Methods 0.000 abstract description 4
- 230000006641 stabilisation Effects 0.000 abstract 1
- 238000011105 stabilization Methods 0.000 abstract 1
- 238000001556 precipitation Methods 0.000 description 6
- 238000007619 statistical method Methods 0.000 description 5
- 238000011160 research Methods 0.000 description 4
- 230000000694 effects Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 2
- 241000282414 Homo sapiens Species 0.000 description 1
- 241000695274 Processa Species 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000002238 attenuated effect Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000013499 data model Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000010348 incorporation Methods 0.000 description 1
- 238000012933 kinetic analysis Methods 0.000 description 1
- 238000007726 management method Methods 0.000 description 1
- 238000012067 mathematical method Methods 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- BULVZWIRKLYCBC-UHFFFAOYSA-N phorate Chemical compound CCOP(=S)(OCC)SCSCC BULVZWIRKLYCBC-UHFFFAOYSA-N 0.000 description 1
- 230000001681 protective effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F16/00—Information retrieval; Database structures therefor; File system structures therefor
- G06F16/20—Information retrieval; Database structures therefor; File system structures therefor of structured data, e.g. relational data
- G06F16/23—Updating
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01W—METEOROLOGY
- G01W1/00—Meteorology
- G01W1/10—Devices for predicting weather conditions
-
- 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
- Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
- Y04S—SYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
- Y04S10/00—Systems supporting electrical power generation, transmission or distribution
- Y04S10/50—Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications
Landscapes
- Engineering & Computer Science (AREA)
- Environmental & Geological Engineering (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Databases & Information Systems (AREA)
- Data Mining & Analysis (AREA)
- Life Sciences & Earth Sciences (AREA)
- Atmospheric Sciences (AREA)
- Biodiversity & Conservation Biology (AREA)
- Ecology (AREA)
- Environmental Sciences (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了预报预测相结合的气象栅格数据更新方法,该方法先从服务器上获取气象栅格数据,并定义气象栅格数据的时间序列;再通过平稳化检测、模型识别选择、模型定阶、模型参数估计建立预测模型;最后对未来气象栅格数据预报预测,实现数据更新。本发明以时间序列模型来对气象栅格数据预测,对历史观测数据进行分析来预测未来的气象数据,根据数据之间时间的依赖性,增加了预测预报的精准性。
Description
技术领域
本发明涉及预报预测相结合的气象栅格数据更新方法,属于气象预测预报技术领域。
背景技术
气象的预测预报对人类的工作生活至关重要,从耕作出行到航天航空,从日常生活到经济建设,几乎都与气象有关。随着越来越多现代化气象装备的出现,使得气象站更加专业化、现代化,气象学的研究技术得到了快速的革新,预报精度得到了极大的提升。气象预测主要是根据定点观测的风速、气温、气压、降雨量等气象数据来确定未来一段时间内气象的变化。
当前的气象预测技术主要分为经验统计法、动力学分析以及模式输出统计法。经验统计法以数理统计学作为基础,通过对气象要素的历史变化做全面***的统计分析,再结合相关因素的影响,总结出一定的变化规律,确定局部气象要素的预报方程;动力学分析是基于应用物理学和数学方法,从理论上研究大气运动、气象***演变过程中的动力、热力变化过程,在研究气象***变化过程中引入水、热平衡方程,进而建立相对应的数学模型;模式输出统计法是利用概率统计方法研究大量气象要素历史观测资料,分析气候变化的统计规律建立合理的预报模型,通过统计数学模式来预测未来的气象状况。
气象栅格是国家气象网格项目启动之后提出的概念,即通过大范围的气象数据特征来描述某一个区域内的整体情况。在气象栅格数据中包含了各种不同气象要素的信息,所以气象栅格数据的质量参差不齐,数据分散。当前气象预测的主要方法是数值预报,需要准确的边界条件,而气象栅格数据对边界的划分比较固定或者模糊,这样就没有一个准确的边界条件,对预测出来的结果有一定的偏差,并且预测出来的结果可能会跨越多个气象栅格,不利于栅格数据的管理与生成。
发明内容
本发明所要解决的技术问题是:提供预报预测相结合的气象栅格数据更新方法,以时间序列模型来对气象栅格数据预测,对历史观测数据进行分析来预测未来的气象数据,提高数据之间的依赖性,增加了预测预报的精准性。
本发明为解决上述技术问题采用以下技术方案:
预报预测相结合的气象栅格数据更新方法,包括如下步骤:
步骤1,获取气象栅格数据,对于待更新的气象要素,将气象栅格数据中该气象要素的值按照时间顺序依次排列,得到该气象要素的一组时间序列;
步骤2,采用ADF单位根检验法检验步骤1得到的时间序列是否平稳,若平稳则进入步骤3,若不平稳则对时间序列进行处理,使处理后的时间序列达到平稳性要求;
步骤3,对于达到平稳性要求的时间序列,计算自相关系数和偏相关系数,根据自相关系数和偏相关系数识别待更新的气象要素的时间序列的预测模型,并对预测模型定阶;
步骤4,利用最小二乘法对步骤3识别的预测模型进行参数估计,完成预测模型的建模;
步骤5,利用建模后的预测模型实现时间序列未来时间的气象要素值的预测与更新。
作为本发明的一种优选方案,所述步骤2的具体过程如下:
上述p阶自回归模型对应的特征方程为:
采用假设检验的方式进行检验,原假设为,即时间序列不平稳;备择假设为,即时间序列平稳;令记为,构造ADF统计量,其中,表示的估计,表示估计标准差,当小于临界值时,认为时间序列平稳,反之不平稳;其中临界值通过蒙特卡洛模拟得到;
若判定待更新的气象要素的时间序列不平稳,则对时间序列进行一阶差分法处理,对一阶差分处理后的时间序列用ADF单位根检验法进行检验,若平稳则进入步骤3,若不能达到平稳性要求,则在一阶差分的基础上进行二阶差分处理,并继续用ADF单位根检验法进行检验。
作为本发明的一种优选方案,所述步骤3中,识别待更新的气象要素的时间序列的预测模型,并对预测模型定阶,具体如下:
若时间序列的自相关系数不满足以上自相关系数对应的条件且偏相关系数不满足以上偏相关系数对应的条件,则确定时间序列的预测模型为ARMA(p,q)模型;采用最小信息准则函数定阶法,即AIC法,对ARMA(p,q)模型进行定阶,ARMA(p,q)模型拟合的AIC准则函数为:,为残差方差的估计值,p,q为模型阶数,当AIC函数取值最小时对应的p,q就为模型的阶数。
作为本发明的一种优选方案,所述步骤4的具体过程如下:
当步骤3识别的预测模型为ARMA(p,q)模型时:,参数估计即求解自回归系数和移动平均系数,为残差,根据条件最小二乘法,t ≤0时,x t = 0,得到残差的有限项表达式:,其中为ARMA(p,q)模型逆转形式中的逆函数;则残差平方和,再通过迭代使残差平方和到达最小值,因为,,联立方程即能估计出的值;
同理,当步骤3识别的预测模型为AR(p)模型或MA(q)模型时,同样按照最小二乘法估计出参数。
作为本发明的一种优选方案,所述步骤5的具体过程如下:
当预测模型为AR(p)模型时:
第t +k个时间点的气象要素值x t+k 表示为:
则预测值为:
当预测模型为MA(q)模型时:
第t +k个时间点的气象要素值x t+k 表示为:
则预测值为:
ARMA(p,q)模型是AR(p)模型与MA(q)模型的结合,当预测模型为ARMA(p,q)模型时:
第t +k个时间点的气象要素值x t+k 表示为:
预测值为:
本发明采用以上技术方案与现有技术相比,具有以下技术效果:
1、本发明不同于以往的预测方式,以时间序列模型来对气象栅格数据预测,对历史观测数据进行分析来预测未来的气象数据,提高了气象栅格数据之间的联系,使得预测更新之后的栅格数据精准率更高。
2、本发明采用平稳的时间序列建立时间数据模型,对于非平稳性数据,采用差分法处理离散的时间序列使其平稳化,以满足模型建立的标准。
3、本发明对气象栅格数据的非线性时间序列实现了预测与更新,结构简单,设计合理,易于实现。
附图说明
图1是本发明预报预测相结合的气象栅格数据更新方法的流程图。
具体实施方式
下面详细描述本发明的实施方式,所述实施方式的示例在附图中示出。下面通过参考附图描述的实施方式是示例性的,仅用于解释本发明,而不能解释为对本发明的限制。
本发明采用预测预报相结合的方式,通过对气象栅格数据的时间序列建模,完成数据的预测预报与更新,便于气象栅格数据的更新。如图1所示,为本发明预报预测相结合的气象栅格数据更新方法的流程图,具体过程如下:
第一步,从服务器上获取气象栅格数据,并定义气象栅格数据的时间序列。
从服务器上获取气象栅格数据,若为一个月中的降水值x,设样本空间为S={x},若对每个x∈S,总有一个确定的函数X(t, x),时间t∈T与它相对应,T为一个区间的时间。这样关于每一个时间t n 对应的函数记为X n ,{X 1, X 2,…, X n }可以看作一个随机过程,而每一个函数X n 称为这个随机过程的样本函数。把一个月内的降水值,作为时间序列的一组观测值。
第二步,为气象栅格数据时间序列进行平稳性检验和处理。
在设计预测模型之前,需要考虑气象栅格数据的平稳性,采用经典的ADF单位根检验法检测时间序列是否平稳,若平稳则进行下一步操作,若不平稳就将时间序列通过一阶差分法、二阶差分法处理,使其达到平稳性的要求。
一个气象栅格数据的时间序列满足p阶自回归,如下:
其中记为,记为,记为。可简记为。ADF单位根检测法的原理是基于p阶自回归模型的特征方程,上面假设的p阶自回归模型所对应的特征方程为,如果该方程的所有特征根都在单位圆内,则序列平稳,反之不平稳。而在时,该特征方程必然存在至少一个单位根。所以在检测时间序列是否平稳时,如果,则为平稳;如果则为不平稳。当时间序列包含单位根时,这个回归系数是有偏差的。所以采用假设检验的方式,原假设为,即,此刻时间序列不平稳:备择假设为,即时间序列平稳。进一步构造ADF统计量,其中表示的估计,表示估计标准差。当小于临界值,则认为时间序列平稳,反之不平稳。其中的临界值是在特定气象栅格数据样本通过蒙特卡洛模拟得到。通过上述方法判断时间序列是否平稳。
除此以外,包含截距项和同时包含截距项和趋势项的另外两种形式为:
其中a 0,a 2为常数,这两种形式也采用相同的检测方法。
若气象栅格数据的时间序列不平稳,就对序列进行一阶差分法处理,一阶差分形式为:。如果一阶差分处理后的时间序列还不能达到平稳性要求,则在一阶差分的基础上进行二阶差分处理,二阶差分形式为:。处理之后还是用ADF单位根检测法检测。
选择用来预测预报的气象栅格数据绝大部分情况都是平稳的。当时间序列判断为平稳,就可以进行模型识别与定阶。
第三步,利用自相关系数与偏相关系数实现模型识别与定阶。
根据第二步中,一个月降水值的平稳性时间序列,以这个序列作为样本,求出它的自相关系数和偏相关系数。自相关系数求解如下:
偏相关系数求解如下:
ARMA模型相关性特征如下表:
拖尾是指时间序列以指数速度向零衰减;截尾是指时间序列从某个时点突然变得很小。在实际运用中,这个定阶原则在操作上有一定的困难。由于气象栅格数据样本的随机性,样本的相关系数不会呈现出理论上截尾的情况,它的时间序列的自相关系数与偏相关系数仍会小幅度波动。当n足够大的时候,这个降水样本的自相关系数和偏相关系数服从正态分布:即,,根据正态分布的性质有:
所以可以用2倍标准差范围来判断所得气象栅格数据的时间序列服从什么预测模型。
如果所得降水时间序列的自相关系数,在最初的q项明显超过2倍标准差范围,即;而在之后几乎95%的自相关系数都在2倍标准差范围以内,即,并且自相关系数衰减为0及其附近小值的过程很突然,这就表现出截尾性,称为自相关系数截尾,且截尾阶数为q,所以这个时间序列的预测模型可确定为阶数为q的MA(q)模型。同理时间序列的偏相关系数,在最初的p项明显超过2倍标准差范围,在之后几乎95%的偏相关系数都在2倍标准差范围以内,并且偏相关系数衰减为0及其附近小值的过程很突然,则偏相关系数截尾,且截尾阶数为p,这个时间序列的预测模型可确定为阶数为p的AR(p)模型。如果时间序列的自相关系数和偏相关系数均有5%以上超过2倍标准差范围,或者自相关系数和偏相关系数衰减为0附近小值的过程比较缓慢,这就是自相关系数和偏相关系数拖尾,所以预测模型为ARMA(p,q)。
在识别预测模型之后,需要对预测模型定阶,模型的阶数对预测结果的精确度有较大的影响。其中AR(p)、MA(q)预测模型在识别中就已经完成定阶,而ARMA(p,q)预测模型需要采用最小信息准则函数定阶法,即AIC法,来对预测模型进行模型定阶。ARMA(p,q)预测模型拟合的AIC准则函数为:。其中n为样本容量,为残差方差的估计值,p,q为模型阶数。当AIC函数取值最小时,此时的p,q就为模型的阶数。
第四步,利用最小二乘法实现模型参数估计。
在对气象栅格数据时间序列确定好模型与定阶后,接下来就需要通过最小二乘法利用观测值来确定预测模型中的未知参数值。最小二乘法的主要思想是计算出残差平方和,当计算结果达到最小时,所对应的参数值即为模型的参数估计值。在ARMA(p,q)模型中:,与分别为自回归系数与移动平均系数,也就是需要求解的系数。根据条件最小二乘法,之前没有观测到的序列值为0,也就是,t ≤ 0时,x t =0,可以得到残差的有限项表达式:,其中为ARMA(p,q)模型逆转形式中的逆函数。所以残差平方和,再通过迭代使残差平方和到达最小值。而又因为,,联立方程即可估计出的值。因为ARMA(p,q)包含了AR(p)、MA(q)模型的性质,同理AR(p)、MA(q)模型同样也按照最小二乘法估计出参数。
到此就完成了气象栅格数据序列的预测模型的建模。
第五步,通过时间序列建模完成气象栅格数据的预报预测与更新。
在已知的历史气象栅格数据信息,对时间序列未来某个时间的x t+k (k =1,2,…)进行预测。一般是根据预测误差的平方达到最小时来预测,可以通过求解条件期望的方式来预测。根据第三步中,选择适合的模型实现预测。
在AR(p)模型中:
第t +k个时间点的x t+k 可以表示为:
因此预测值为:
在MA(q)模型中:
第t +k个时间点的x t+k 可以表示为:
预测值为:
ARMA(p,q)模型是AR(p)模型与MA(q)模型的结合,在ARMA(p,q)模型中:
第t +k个时间点的x t+k 可以表示为:
预测值为AR(p)模型与MA(q)模型叠加:
用已获取的过去气象栅格数据,作为初始样本来预测未来几天的气象栅格数据,例如用前30天的数据预测第31天的数据,再用前29天和未来预测第31天的数据预测第32天的数据,以此类推,达到气象栅格数据的预测预报效果,并实现气象栅格数据更新。
以上实施例仅为说明本发明的技术思想,不能以此限定本发明的保护范围,凡是按照本发明提出的技术思想,在技术方案基础上所做的任何改动,均落入本发明保护范围之内。
Claims (6)
1.预报预测相结合的气象栅格数据更新方法,其特征在于,包括如下步骤:
步骤1,获取气象栅格数据,对于待更新的气象要素,将气象栅格数据中该气象要素的值按照时间顺序依次排列,得到该气象要素的一组时间序列;
步骤2,采用ADF单位根检验法检验步骤1得到的时间序列是否平稳,若平稳则进入步骤3,若不平稳则对时间序列进行处理,使处理后的时间序列达到平稳性要求;
步骤3,对于达到平稳性要求的时间序列,计算自相关系数和偏相关系数,根据自相关系数和偏相关系数识别待更新的气象要素的时间序列的预测模型,并对预测模型定阶;
步骤4,利用最小二乘法对步骤3识别的预测模型进行参数估计,完成预测模型的建模;
步骤5,利用建模后的预测模型实现时间序列未来时间的气象要素值的预测与更新。
2.根据权利要求1所述的预报预测相结合的气象栅格数据更新方法,其特征在于,所述步骤2的具体过程如下:
上述p阶自回归模型对应的特征方程为:
采用假设检验的方式进行检验,原假设为,即时间序列不平稳;
备择假设为,即时间序列平稳;令记为,构造ADF
统计量,其中,表示的估计,表示估计标准差,当小于临界值时,认为时间序列
平稳,反之不平稳;其中临界值通过蒙特卡洛模拟得到;
若判定待更新的气象要素的时间序列不平稳,则对时间序列进行一阶差分法处理,对一阶差分处理后的时间序列用ADF单位根检验法进行检验,若平稳则进入步骤3,若不能达到平稳性要求,则在一阶差分的基础上进行二阶差分处理,并继续用ADF单位根检验法进行检验。
4.根据权利要求3所述的预报预测相结合的气象栅格数据更新方法,其特征在于,所述步骤3中,识别待更新的气象要素的时间序列的预测模型,并对预测模型定阶,具体如下:
6.根据权利要求5所述的预报预测相结合的气象栅格数据更新方法,其特征在于,所述步骤5的具体过程如下:
当预测模型为AR(p)模型时:
第t +k个时间点的气象要素值x t +k 表示为:
则预测值为:
当预测模型为MA(q)模型时:
第t +k个时间点的气象要素值x t +k 表示为:
则预测值为:
ARMA(p,q)模型是AR(p)模型与MA(q)模型的结合,当预测模型为ARMA(p,q)模型时:
第t +k个时间点的气象要素值x t +k 表示为:
预测值为:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210464764.0A CN114564487B (zh) | 2022-04-29 | 2022-04-29 | 预报预测相结合的气象栅格数据更新方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210464764.0A CN114564487B (zh) | 2022-04-29 | 2022-04-29 | 预报预测相结合的气象栅格数据更新方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114564487A true CN114564487A (zh) | 2022-05-31 |
CN114564487B CN114564487B (zh) | 2022-08-02 |
Family
ID=81721359
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210464764.0A Active CN114564487B (zh) | 2022-04-29 | 2022-04-29 | 预报预测相结合的气象栅格数据更新方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114564487B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117111180A (zh) * | 2023-08-15 | 2023-11-24 | 北京大学重庆大数据研究院 | 模式后气象要素发生概率预测方法、可读存储介质及计算设备 |
CN118094127A (zh) * | 2024-04-29 | 2024-05-28 | 南京信息工程大学 | 基于高空数据与gm灰色迭代的降水量预测方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20150302313A1 (en) * | 2014-04-22 | 2015-10-22 | State Grid Corporation Of China | Method of predicating ultra-short-term wind power based on self-learning composite data source |
CN110689179A (zh) * | 2019-09-18 | 2020-01-14 | 北京工商大学 | 基于时空序列混合模型的水华预测方法 |
CN114118511A (zh) * | 2021-08-31 | 2022-03-01 | 中国资源卫星应用中心 | 基于云量预测信息的大区域多星联合覆盖有效性评估方法 |
-
2022
- 2022-04-29 CN CN202210464764.0A patent/CN114564487B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20150302313A1 (en) * | 2014-04-22 | 2015-10-22 | State Grid Corporation Of China | Method of predicating ultra-short-term wind power based on self-learning composite data source |
CN110689179A (zh) * | 2019-09-18 | 2020-01-14 | 北京工商大学 | 基于时空序列混合模型的水华预测方法 |
CN114118511A (zh) * | 2021-08-31 | 2022-03-01 | 中国资源卫星应用中心 | 基于云量预测信息的大区域多星联合覆盖有效性评估方法 |
Non-Patent Citations (1)
Title |
---|
毛钰嘉: "基于时间序列的组合预测模型研究——以江西省CPI为例", 《中国优秀博硕士学位论文全文数据库(硕士)基础科学辑》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117111180A (zh) * | 2023-08-15 | 2023-11-24 | 北京大学重庆大数据研究院 | 模式后气象要素发生概率预测方法、可读存储介质及计算设备 |
CN118094127A (zh) * | 2024-04-29 | 2024-05-28 | 南京信息工程大学 | 基于高空数据与gm灰色迭代的降水量预测方法 |
Also Published As
Publication number | Publication date |
---|---|
CN114564487B (zh) | 2022-08-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Loukatou et al. | Stochastic wind speed modelling for estimation of expected wind power output | |
CN107153874B (zh) | 水质预测方法及*** | |
CN114564487B (zh) | 预报预测相结合的气象栅格数据更新方法 | |
Shamshad et al. | First and second order Markov chain models for synthetic generation of wind speed time series | |
Sartini et al. | Comparing different extreme wave analysis models for wave climate assessment along the Italian coast | |
JP2009294969A (ja) | 需要予測方法および需要予測装置 | |
CN111815065B (zh) | 基于长短时记忆神经网络的短期电力负荷预测方法 | |
CN112285807A (zh) | 一种气象信息预测方法及装置 | |
CN112116149B (zh) | 一种考虑预报不确定性关联演化特征的多站中长期径流滚动概率预测方法 | |
Sahu et al. | A comparison of Bayesian models for daily ozone concentration levels | |
US10770898B2 (en) | Methods and systems for energy use normalization and forecasting | |
CN109558968B (zh) | 风电场出力相关性分析方法及装置 | |
Chen et al. | Postprocessing of ensemble weather forecasts using a stochastic weather generator | |
Chen et al. | Combining stochastic weather generation and ensemble weather forecasts for short-term streamflow prediction | |
CN110084403B (zh) | 架空导线载流量的多时段联合概率密度预测方法及*** | |
KR101463492B1 (ko) | 비정상성 분위 사상을 적용한 전지구 기후모델의 오차보정방법 | |
CN114357670A (zh) | 一种基于bls和自编码器的配电网用电数据异常预警方法 | |
CN115238947A (zh) | 气候变化下旱涝急转事件的社会经济暴露度预估方法 | |
Dkengne et al. | A limiting distribution for maxima of discrete stationary triangular arrays with an application to risk due to avalanches | |
KR101463493B1 (ko) | 추계학적 태풍모의를 적용한 전지구 기후모델의 보완방법 | |
KR101409316B1 (ko) | 다중모델 확률예측 시스템 및 방법 | |
CN112580844A (zh) | 气象数据处理方法、装置、设备及计算机可读存储介质 | |
CN104679939B (zh) | 一种飞机设计经济可承受性评估过程的多准则决策方法 | |
Das et al. | Short-term prediction of land surface temperature using multifractal detrended fluctuation analysis | |
CN115577854A (zh) | 一种基于eemd-rbf组合的分位数回归风速区间预测方法 |
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 |