CN116882572B - 一种降水年代际变化的预测方法、装置及存储介质 - Google Patents
一种降水年代际变化的预测方法、装置及存储介质 Download PDFInfo
- Publication number
- CN116882572B CN116882572B CN202310857764.1A CN202310857764A CN116882572B CN 116882572 B CN116882572 B CN 116882572B CN 202310857764 A CN202310857764 A CN 202310857764A CN 116882572 B CN116882572 B CN 116882572B
- Authority
- CN
- China
- Prior art keywords
- precipitation
- chronology
- component
- preset
- order
- 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.)
- Active
Links
- 238000001556 precipitation Methods 0.000 title claims abstract description 142
- 238000000034 method Methods 0.000 title claims abstract description 41
- 230000000737 periodic effect Effects 0.000 claims abstract description 30
- 238000004364 calculation method Methods 0.000 claims abstract description 21
- 238000007781 pre-processing Methods 0.000 claims abstract description 15
- 238000012937 correction Methods 0.000 claims description 32
- 238000001228 spectrum Methods 0.000 claims description 13
- 238000001914 filtration Methods 0.000 claims description 11
- 230000006870 function Effects 0.000 claims description 11
- 238000004590 computer program Methods 0.000 claims description 7
- 230000003247 decreasing effect Effects 0.000 claims description 5
- 238000013519 translation Methods 0.000 claims description 5
- 238000011144 upstream manufacturing Methods 0.000 description 20
- 230000000694 effects Effects 0.000 description 10
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 4
- 238000012360 testing method Methods 0.000 description 3
- 230000007774 longterm Effects 0.000 description 2
- 238000005192 partition Methods 0.000 description 2
- 239000013074 reference sample Substances 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 230000007547 defect Effects 0.000 description 1
- 230000002950 deficient Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 239000000284 extract Substances 0.000 description 1
- 230000001788 irregular Effects 0.000 description 1
- 238000010248 power generation Methods 0.000 description 1
- 230000001105 regulatory effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000010183 spectrum analysis Methods 0.000 description 1
- 238000011282 treatment Methods 0.000 description 1
- 239000013598 vector Substances 0.000 description 1
- 238000010792 warming Methods 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
- G01—MEASURING; TESTING
- G01W—METEOROLOGY
- G01W1/00—Meteorology
- G01W1/10—Devices for predicting weather conditions
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/14—Fourier, Walsh or analogous domain transformations, e.g. Laplace, Hilbert, Karhunen-Loeve, transforms
-
- 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)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Theoretical Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Business, Economics & Management (AREA)
- Environmental & Geological Engineering (AREA)
- Strategic Management (AREA)
- General Engineering & Computer Science (AREA)
- Software Systems (AREA)
- Databases & Information Systems (AREA)
- Algebra (AREA)
- Economics (AREA)
- Human Resources & Organizations (AREA)
- Life Sciences & Earth Sciences (AREA)
- Tourism & Hospitality (AREA)
- Environmental Sciences (AREA)
- Entrepreneurship & Innovation (AREA)
- Biodiversity & Conservation Biology (AREA)
- Atmospheric Sciences (AREA)
- Game Theory and Decision Science (AREA)
- Ecology (AREA)
- Marketing (AREA)
- Development Economics (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- General Business, Economics & Management (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本申请涉及一种降水年代际变化的预测方法、装置及储存介质,其包括:通过对目标时间内的历史降水资料进行预处理,得到降水时间序列对应的年代际分量;根据所述年代际分量遍历预设的阶数集,得到傅里叶级数的最优阶数,所述阶数集包括若干个预设的傅里叶级数的阶数;通过所述最优阶数下的傅里叶级数对所述年代际分量进行拟合,得到降水年代际分量的周期变化特征,根据所述周期变化特征预测降水的年代际变化。本申请具有满足降水年代际变化预测简便易行、计算成本低的要求。
Description
技术领域
本申请涉及降水年代际变化预测技术领域,尤其是涉及一种降水年代际变化的预测方法、装置及存储介质。
背景技术
水资源的多寡以及水电站的发电量通常都与径流量的大小相关,而径流量的大小在很大程度上受到降水的调节。降水的未来变化特征一直是众多学者关注的重点。近年来,全球变暖日益加剧,导致大气环流变化加剧,意味着降水变率显著增强,加大了降水预测的难度。
目前,国内对降水的预测多采用数值模拟的方法,但是该方法通常需要巨型计算机进行解算,计算费用高,耗时冗长。
发明内容
针对上述问题,本发明的目的是提供一种降水年代际变化的预测方法、装置及存储介质,弥补了降水年代际变化预测计算费用高、耗时冗长的不足,更好的满足了降水年代际变化预测简便易行、计算成本低的要求。
为实现上述目的,本发明采取以下技术方案:
一种降水年代际变化的预测方法,包括以下步骤:
通过对目标时间内的历史降水资料进行预处理,得到降水时间序列对应的年代际分量;
根据所述年代际分量遍历预设的阶数集,得到傅里叶级数的最优阶数,所述阶数集包括若干个预设的傅里叶级数的阶数;
通过所述最优阶数下的傅里叶级数对所述年代际分量进行拟合,得到降水年代际分量的周期变化特征,根据所述周期变化特征预测降水的年代际变化。
进一步地,所述根据所述年代际分量遍历预设的阶数集,得到傅里叶级数的最优阶数包括将预设的阶数集分别输入傅里叶级数中,得到当前阶数对应的校正拟合度,判断相邻阶数对应的所述校正拟合度之间的差值是否大于预设的第一拟合度阈值,选择大于预设第一拟合度阈值的阶数中最大的阶数为最优阶数。
进一步地,所述根据所述年代际分量遍历预设的阶数集,得到傅里叶级数的最优阶数包括将预设的阶数集分别代入傅里叶级数,得到当前阶数对应的校正拟合度,判断阶数对应的所述校正拟合度是否大于预设的第二拟合度阈值,选择大于预设第二拟合度阈值的阶数中最小的阶数为最优阶数。
进一步地,所述对目标时间内的历史降水资料进行预处理包括:
根据目标时间内的历史降水资料,采集降水时间序列;
通过滑动平均对所述降水时间序列进行低通滤波,得到年代际分量,所述滑动平均的时间窗口为预设的固定值。
进一步地,所述通过滑动平均对所述降水时间序列进行低通滤波包括通过顺序逐期增减新旧降水时间序列求算移动平均值Yn,计算公式为:
其中,X(n-k)为降水时间序列;k,M为预设的滑动平均的参数。
进一步地,所述滑动平均的时间窗口为9,9年滑动平均Yn的表达式如下:
上述公式中,X(n-k)为降水时间序列;k为预设的滑动平均的参数。
进一步地,所述得到年代际分量之后,所述对目标时间内的历史降水资料进行预处理还包括通过小波分析对所述年代际分量进行周期分析,确定所述年代际分量在预设的阈值范围内。
为实现上述目的,本申请另一方面还提供一种降水年代际变化预测装置,包括:采集模块,用于通过对目标时间内的历史降水资料进行预处理,得到降水时间序列对应的年代际分量;
遍历模块,用于根据所述年代际分量遍历预设的阶数集,得到傅里叶级数的最优阶数,所述阶数集包括若干个预设的傅里叶级数的阶数;
拟合模块,用于通过所述最优阶数下的傅里叶级数对所述年代际分量进行拟合,得到降水年代际分量的周期变化特征,根据所述周期变化特征预测降水的年代际变化。
为实现上述目的,本申请另一方面还提供一种计算机可读存储介质,所述计算机可读存储介质包括计算机程序,所述计算机程序用于实现上述的方法。
为实现上述目的,本申请另一方面还提供一种计算设备,所述计算设备包括存储器和处理器,所述存储器用于存储计算机程序,当所述计算机程序被所述处理器执行时,实现上述的方法。
本发明由于采取以上技术方案,其具有以下优点:
通过对历史降水资料进行预处理,提取降水时间序列的主要变化部分,滤除降水时间序列中的随机部分,得到该地区年代际分量的序列,从而抓住了降水的主要变化特征;通过最优阶数下的傅里叶级数对年代际分量进行拟合,计算量小、计算成本低、耗时短;傅里叶级数具有明显的周期性,能很好地复现降水年代际的周期变化特征,避免了线性拟合方法导致拟合结果发散的问题,能够较好预测降水的年代际变化,简便易行。
附图说明
图1是本申请实施例提供的降水年代际变化预测方法的流程示意图。
图2是本申请实施例1提供的长江上游降水年代际变化预测方法的流程示意图。
图3是本申请实施例1中长江上游夏季降水区域平均的年代际分量(单位:mm/d)曲线示意图。
图4是本申请实施例1中长江上游夏季降水区域年代际分量的小波功率谱示意图。
图5是本申请实施例1中长江上游降水指数的傅里叶级数拟合的校正拟合度和拟合度示意图。
图6是本申请实施例1中长江上游降水指数的年代际分量及其5阶傅里叶级数拟合曲线示意图。
具体实施方式
下面将参照附图更详细地描述本发明的示例性实施方式。虽然附图中显示了本发明的示例性实施方式,然而应当理解,可以以各种形式实现本发明而不应被这里阐述的实施方式所限制。相反,提供这些实施方式是为了能够更透彻地理解本发明,并且能够将本发明的范围完整的传达给本领域的技术人员。
降水的未来变化特征一直是众多学者关注的重点,降水变化除了具有年际变化外,年代际变化也是降水变化的一个非常重要的时间尺度,年际变化的周期通常为1年,而年代际变化的周期通常为10年。因此,对降水年代际变化进行预测也是十分有必要的。
本申请针对降水预测难度大,通过数值模拟的方法预测降水,计算费用高,耗时冗长的特点,提出了一种通过傅里叶级数拟合复现降水周期变化特征的降水年代际变化的预测方法。
参照图1,本申请实施例公开一种降水年代际变化的预测方法,包括以下步骤:
S101:通过对目标时间内的历史降水资料进行预处理,得到降水时间序列对应的年代际分量。
具体地,通过对目标时间内的历史降水资料进行预处理包括;
S201:根据目标时间内的历史降水资料,采集目标时间内的降水时间序列。
S202:通过滑动平均对降水时间序列进行低通滤波,得到年代际分量,滑动平均的时间窗口为预设的固定值。
其中,滑动平均的时间窗口为奇数,通过滑动平均对降水时间序列进行低通滤波包括通过顺序逐期增减新旧降水时间序列求算移动平均值Yn,计算公式为:
上述公式中,X(n-k)为降水时间序列;k,M为预设的滑动平均的参数。
在其它实施方式中,滑动平均的时间窗口可以为偶数,其滑动平均结果Yn的计算公式如下:
上述公式中,X(n-k)为降水时间序列;k,M为预设的滑动平均的参数。
滑动平均法是指在简单平均数法基础上,通过顺序逐期增减新旧数据求算移动平均值,借以消除高频变动因素,找出事物的长期变化趋势,并据此进行预测的方法。
在本发明的其中一种实施方式中,滑动平均的时间窗口为9年,采用9年滑动平均提取降水时间序列的年代际分量,其9年滑动平均Yn的表达式如下:
上述公式中,X(n-k)为降水时间序列;k为预设的滑动平均的参数。
通过滑动平均的方法提取降水时间序列的年代际变化分量,滤除降水时间序列中的随机变化分量,能够获得降水的主要变化特征。其中,滑动平均法可用卡尔曼滤波法替换。
S203:通过小波分析对年代际分量进行周期分析,确定所得年代际分量在预设的阈值范围内。
其中,通过小波分析对年代际分量进行周期分析包括通过设置小波母函数为Morlet小波将年代际分量分解到时间频率域内,得到时间系列的显著波动模式,即周期变化特征。预设的阈值范围是本领域技术人员通过现有技术可以确定的。在其它实施方式中,小波分析可替换为功率谱分析。
Morlet小波既是非正交的,又是Gaussian调节的复数小波,能够很好地平衡时间与频率的局部变化。Morlet小波的表达式如下:
在实数域中,信号f(t)的小波变化离散表达式如下:
小波功率谱的计算公式如下:
Pf(s)=∫|Wf(s,τ)|2dτ
上述公式中,t是时间,ω是无量纲频率,s是伸缩因子,Δt是样本的间隔,i是虚数,τ是平移量。
具体地,通过信号f(t)的小波变化离散表达式,将年代际分量用小波分析转换成频域形式,再对频域形式的年代际分量计算功率谱,对年代际分量进行滤波,以分析年代际分量的周期变化特征。
小波分析将信号分解成一系列小波函数的叠加,这些小波函数都由小波母函数经过平移与尺度伸缩得来的。用这种不规则的小波母函数可以逼近那些非稳定信号中尖锐变化的部分,也可以去逼近离散不连续具有局部特性的信号,从而更为真实的反映信号在某一时间尺度上的变化。
在本发明的其中一种实施方式中,目标时间内的历史降水资料可以取自于英国东英吉利大学气候研究中心(Climatic Research Unit,CRU)的CRU v4.05逐月全球陆地降水资料,资料的空间分辨率为0.5°×0.5°。
S102:根据年代际分量遍历预设的阶数集,得到傅里叶级数的最优阶数,阶数集包括若干个预设的傅里叶级数的阶数;通过最优阶数下的傅里叶级数对年代际分量进行拟合,得到降水年代际分量的周期变化特征,根据周期变化特征预测降水的年代际变化。
通过最优阶数下的傅里叶级数对年代际分量进行拟合,具有计算量小、计算成本低、耗时短的优点;而傅里叶级数具有明显的周期性,能很好地复现降水年代际的周期变化特征,避免了线性拟合方法导致拟合结果发散的问题,能够较好地预测降水的年代际变化,简便易行。
在本发明的其中一种实施方式中,所述步骤S102包括:
S301:构建傅里叶级数的阶数集,阶数集包括若干个预设的傅里叶级数的阶数。
具体地,阶数集中的阶数为自然数,即阶数m依次可以取1、2、3、4……n。
S302:将年代际分量输入MATLAB软件中进行傅里叶级数拟合,得到傅里叶级数拟合曲线。
其中,傅里叶级数的表达式如下:
上述表达式中,m是傅里叶级数的阶数,an是常数项,bn和w都是系数。
傅里叶级数是一类特殊的函数项级数,是以三角函数为基对周期信号的无穷级数展开,即对周期性现象进行数学上的分析。
傅里叶级数拟合是一种基于函数形态和函数参数拟合的方法,可以将实际观察到的数据进行拟合,以及对离散数据进行插值、滤波等处理,能够精准表达复杂曲线。
S303:通过拟合度和校正拟合度度量傅里叶级数拟合的效果。
拟合度是通过比较预测结果与实际发生情况的吻合程度来度量拟合效果。拟合度和校正拟合度都是无量纲的,取值范围为0到1,数值越接近1说明拟合效果越好。其中,拟合度R2等于拟合值和拟合对象相关系数的平方。
由于傅里叶级数变量的个数大于1,拟合度的大小不会随变量的增加而减少,因而拟合度不能反映新增变量是否有意义;而当校正拟合度不再因变量的增加而增加时,证明新增变量是没有意义的,此时校正拟合度和拟合度的差距增大;反之,新增变量是有意义的,此时校正拟合度和拟合度的差距减小。因此,采用校正拟合度来度量拟合效果更为合理。
由上可知,当校正拟合度不再增加时,变量个数的增加已不能有效提高拟合效果,说明此时变量个数是恰当的。因此,利用校正拟合度确定合适的变量个数,合理减少计算量。校正拟合度的计算公式如下:
其中,p为变量个数,n为样本数,R2是拟合度。
S304:根据年代际分量遍历预设的阶数集,得到傅里叶级数的最优阶数。
其中,根据年代际分量遍历预设的阶数集包括将预设的阶数集分别代入傅里叶级数,得到当前阶数对应的校正拟合度,判断相邻阶数对应的校正拟合度之间的差值是否大于预设的第一拟合度阈值,选择大于预设第一拟合度阈值的阶数中最大的阶数为最优阶数。
在本发明的其中一种实施方式中,预设的第一拟合度阈值可以为15%,20%等。相邻阶数对应的校正拟合度之间的差值大于预设的第一拟合度阈值时,表明相邻阶数对应的校正拟合度增长大,而其余相邻阶数对应的年代际分量的拟合效果较为稳定,即可选择大于第一拟合度阈值的最大的阶数为最优阶数。
在本发明的其它实施方式中,根据年代际分量遍历预设的阶数集还可以包括将预设的阶数集分别代入傅里叶级数,得到当前阶数对应的校正拟合度,判断阶数对应的校正拟合度是否大于预设的第二拟合度阈值,选择大于预设第二拟合度阈值的阶数中最小的阶数为最优阶数。
S305:将年代际分量输入MATLAB软件中进行最优阶数下的傅里叶级数拟合,得到傅里叶级数拟合曲线,以预测降水的年代际变化。
傅里叶级数拟合曲线可以较好地体现降水年代际的周期变化特征,从而可以很好地预测降水的年代际变化。并且,拟合值和拟合对象的相关系数(0.9),通过了99%的信度检验。
其中,周期变化特征是指对目标时间内,若干个时间分区对应的降水年代际变化规律进行数字化,形成格式化数据或向量,即降水年代际变化的周期性特征。假设目标时间100年内,降水年代际每隔30年周期变化,将30年一个周期的变化规律进行数字化,形成符合规定格式的数据。
实施例1:
长江上游地区水资源极度匮乏,部分地区的人均水资源量仅为国际缺水警戒线(人均1000立方米)的三分之一左右。此外,长江上游地区是中国最大的水电基地,其发电量约占全国的三分之一。该地区夏季降水存在明显的周期变化特征,其中,该地区1995–2004年的平均降水量显著大于2005–2014年;在1960–1990年呈增加趋势。
因此,本实施方式中,选取长江上游经纬度范围为25°–36°N、100°–112°E的地区作为参考样本,参照图2,长江上游地区降水年代际变化的预测方法,包括以下步骤:S1:根据长江上游地区目标时间1910年-2020年内的夏季历史降水资料,采集长江上游地区的降水时间序列。
S2:根据长江上游地区的降水时间序列通过9年滑动平均,得到长江上游夏季降水的年代际分量;
其中,9年滑动平均结果Yn,计算公式为:
上述公式中,X(n-k)为降水时间序列;k为预设的滑动平均的参数。
参照图3,具体地,将得到的长江上游地区的降水时间序列输入9年滑动平均的计算公式中,以9为滑动平均的周期通过顺序依次计算降水时间序列的滑动平均结果,得到年代际分量及其曲线图。长江上游地区1910年-2020年的夏季降水具有显著的年代际变化特征。
S3:将得到的年代际分量通过小波分析,得到年代际分量的频域形式,根据年代际分量的频域形式通过小波功率谱,得到年代际分量的功率谱,以确定所得年代际分量在预设的阈值范围内。
其中,小波母函数为Morlet小波,Morlet小波的表达式如下:
在实数域中,信号f(t)的小波变化离散表达式如下:
小波功率谱的计算公式如下:
Pf(s)=∫|Wf(s,τ)|2dτ
上述公式中,t是时间,ω是无量纲频率,s是伸缩因子,Δt是样本的间隔,i是虚数,τ是平移量。
参照图4,具体地,将得到的年代际分量通过信号f(t)的小波变化离散表达式,得到年代际分量的频域形式,将年代际分量的频域形式通过小波功率谱,得到年代际分量的功率谱。从小波功率谱的结果可知,长江上游地区夏季降水存在显著的30–50年周期变化,且通过了95%的信度检验。
S4:构建傅里叶级数的阶数集,阶数集中的阶数为1-7的自然数。
S5:根据年代际分量将不同的阶数1-7分别代入傅里叶级数,得到当前阶数对应的校正拟合度,判断相邻阶数对应的校正拟合度之间的差值是否大于预设的第一拟合度阈值15%,选择大于预设第一拟合度阈值的阶数中最大的阶数为最优阶数。
具体地,傅里叶级数的表达式如下:
上述表达式中,m是傅里叶级数的阶数,an是常数项,bn和w都是系数。
校正拟合度的计算公式如下:
上述公式中,p为变量个数,n为样本数,R2是拟合度。
相邻阶数对应的校正拟合度之间的差值大于预设的第一拟合度阈值,证明年代际分量的拟合效果明显增加;而相邻阶数对应的校正拟合度之间的差值小于预设的第一拟合度阈值,证明年代际分量的拟合效果趋于稳定。
参照图5,长江上游地区年代际分量的傅里叶级数拟合的校正拟合度随着阶数的增加,校正拟合度增加,长江上游地区降水的拟合效果越好。而当阶数从4到5时,年代际分量的傅里叶级数拟合的校正拟合度陡增,阶数大于5以后,校正拟合度不再明显增加,降水指数的拟合效果趋于稳定,此时校正拟合度约等于80%,得到阶数5为最优阶数。
S6:根据年代际分量通过5阶傅里叶级数对年代际分量进行拟合,得到降水年代际分量的周期变化特征,根据周期变化特征预测降水的年代际变化。
参照图6,具体地,将得到的1910年-2020年的年代际分量通过5阶傅里叶级数拟合,得到长江上游参考样本地区夏季降水的周期变化特征,拟合值和拟合对象的相关系数(0.9),通过了99%的信度检验。
在本实施方式中,采样时间为1910年-2020年,降水年代际每隔30-50年周期变化,将一个周期的变化规律进行数字化,形成曲线,根据拟合效果,可以得到1910年-2020年降水年代际的变化特征,从而可以预测未来的降水年代际变化特征。长江上游地区夏季的降水量在2020年之后开始增加,而在2030年达到峰值之后开始减少,也就是说该地区的降水量将在2025-2035年高于降水的多年平均值。
本申请还提供一种降水年代际变化预测装置,包括:
采集模块,用于通过对目标时间内的历史降水资料进行预处理,得到降水时间序列对应的年代际分量;
遍历模块,用于根据年代际分量遍历预设的阶数集,得到傅里叶级数的最优阶数,阶数集包括若干个预设的傅里叶级数的阶数;
拟合模块,用于通过最优阶数下的傅里叶级数对年代际分量进行拟合,得到降水年代际分量的周期变化特征,根据周期变化特征预测降水的年代际变化。
本实施方式中提供的装置是用于执行上述方法实施例的,具体流程和详细内容参考上述实施例,此处不再赘述。
本申请还提供一种计算机可读存储介质,所述计算机可读存储介质包括计算机程序,所述计算机程序用于实现上述的方法。
本申请还提供一种计算设备,所述计算设备包括存储器和处理器,所述存储器用于存储计算机程序,当所述计算机程序被所述处理器执行时,实现上述的方法。
由上可见,本申请提供的技术方案通过滑动平均的方法提取降水时间序列的主要变化部分,滤除降水时间序列中的随机部分,得到该地区年代际分量的序列,从而抓住了降水的主要变化特征;通过小波分析将年代际分量分解到时间频率域内,通过小波功率谱得到年代际分量的功率谱,得到时间系列的显著波动模式,从而确定所得年代际分量的合理性;根据年代际分量遍历预设的阶数集,得到傅里叶级数的最优阶数;通过最优阶数下的傅里叶级数对年代际分量进行拟合,得到降水年代际的周期变化特征,傅里叶级数具有明显的周期性,能很好地复现降水年代际的周期变化特征,避免了线性拟合方法导致拟合结果发散的问题,从而能够较好的预测降水年代际变化,简便易行;同时,最优阶数下的傅里叶级数拟合具有计算量小、计算成本低、耗时短的优点。
以上均为本申请的较佳实施例,并非依此限制本申请的保护范围,故:凡依本申请的结构、形状、原理所做的等效变化,均应涵盖于本申请的保护范围之内。
Claims (7)
1.一种降水年代际变化的预测方法,其特征在于,包括以下步骤:
通过对目标时间内的历史降水资料进行预处理,得到降水时间序列对应的年代际分量,所述对目标时间内的历史降水资料进行预处理包括,根据目标时间内的历史降水资料,采集降水时间序列,通过滑动平均对所述降水时间序列进行低通滤波,得到年代际分量,所述滑动平均的时间窗口为预设的固定值,所述通过滑动平均对所述降水时间序列进行低通滤波包括通过顺序逐期增减新旧降水时间序列求算移动平均值Yn,计算公式为:
其中,X(n-k)为降水时间序列;k,M为预设的滑动平均的参数;
所述得到年代际分量之后,所述对目标时间内的历史降水资料进行预处理还包括通过小波分析对所述年代际分量进行周期分析,确定所述年代际分量在预设的阈值范围内,其中,小波母函数为Morlet小波,Morlet小波的表达式如下:
在实数域中,信号f(t)的小波变化离散表达式如下:
小波功率谱的计算公式如下:
Pf(s)=∫|Wf(s,τ)|2dτ
上述公式中,t是时间,ω是无量纲频率,s是伸缩因子,Δt是样本的间隔,i是虚数,τ是平移量;
根据所述年代际分量遍历预设的阶数集,得到傅里叶级数的最优阶数,所述阶数集包括若干个预设的傅里叶级数的阶数;
通过所述最优阶数下的傅里叶级数对所述年代际分量进行拟合,得到降水年代际分量的周期变化特征,根据所述周期变化特征预测降水的年代际变化。
2.根据权利要求1所述的一种降水年代际变化的预测方法,其特征在于:所述根据所述年代际分量遍历预设的阶数集,得到傅里叶级数的最优阶数包括将预设的阶数集分别输入傅里叶级数中,得到当前阶数对应的校正拟合度,判断相邻阶数对应的所述校正拟合度之间的差值是否大于预设的第一拟合度阈值,选择大于预设第一拟合度阈值的阶数中最大的阶数为最优阶数。
3.根据权利要求1所述的一种降水年代际变化的预测方法,其特征在于:所述根据所述年代际分量遍历预设的阶数集,得到傅里叶级数的最优阶数包括将预设的阶数集分别代入傅里叶级数,得到当前阶数对应的校正拟合度,判断阶数对应的所述校正拟合度是否大于预设的第二拟合度阈值,选择大于预设第二拟合度阈值的阶数中最小的阶数为最优阶数。
4.根据权利要求1所述的一种降水年代际变化的预测方法,其特征在于,所述滑动平均的时间窗口为9,9年滑动平均Yn的表达式如下:
上述公式中,X(n-k)为降水时间序列;k为预设的滑动平均的参数。
5.一种降水年代际变化预测装置,其特征在于,包括:
采集模块,用于通过对目标时间内的历史降水资料进行预处理,得到降水时间序列对应的年代际分量,所述对目标时间内的历史降水资料进行预处理包括,根据目标时间内的历史降水资料,采集降水时间序列,通过滑动平均对所述降水时间序列进行低通滤波,得到年代际分量,所述滑动平均的时间窗口为预设的固定值,所述通过滑动平均对所述降水时间序列进行低通滤波包括通过顺序逐期增减新旧降水时间序列求算移动平均值Yn,计算公式为:
其中,X(n-k)为降水时间序列;k,M为预设的滑动平均的参数;
所述得到年代际分量之后,所述对目标时间内的历史降水资料进行预处理还包括通过小波分析对所述年代际分量进行周期分析,确定所述年代际分量在预设的阈值范围内,其中,小波母函数为Morlet小波,Morlet小波的表达式如下:
在实数域中,信号f(t)的小波变化离散表达式如下:
小波功率谱的计算公式如下:
Pf(s)=∫|Wf(s,τ)|2dτ
上述公式中,t是时间,ω是无量纲频率,s是伸缩因子,Δt是样本的间隔,i是虚数,τ是平移量;
遍历模块,用于根据所述年代际分量遍历预设的阶数集,得到傅里叶级数的最优阶数,所述阶数集包括若干个预设的傅里叶级数的阶数;
拟合模块,用于通过所述最优阶数下的傅里叶级数对所述年代际分量进行拟合,得到降水年代际分量的周期变化特征,根据所述周期变化特征预测降水的年代际变化。
6.一种计算机可读存储介质,其特征在于,所述计算机可读存储介质包括计算机程序,所述计算机程序用于实现如权利要求1至4中任一权利要求所述的方法。
7.一种计算设备,其特征在于,所述计算设备包括存储器和处理器,所述存储器用于存储计算机程序,当所述计算机程序被所述处理器执行时,实现如权利要求1至4中任一权利要求所述的方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310857764.1A CN116882572B (zh) | 2023-07-12 | 2023-07-12 | 一种降水年代际变化的预测方法、装置及存储介质 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310857764.1A CN116882572B (zh) | 2023-07-12 | 2023-07-12 | 一种降水年代际变化的预测方法、装置及存储介质 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN116882572A CN116882572A (zh) | 2023-10-13 |
CN116882572B true CN116882572B (zh) | 2024-01-30 |
Family
ID=88260031
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310857764.1A Active CN116882572B (zh) | 2023-07-12 | 2023-07-12 | 一种降水年代际变化的预测方法、装置及存储介质 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116882572B (zh) |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110348040A (zh) * | 2019-05-21 | 2019-10-18 | 合肥通用机械研究院有限公司 | 一种基于均方收敛条件的不均匀沉降函数拟合方法 |
KR102409155B1 (ko) * | 2021-10-08 | 2022-06-16 | 주식회사 다올 | Lstm 기반 지하수 수위 예측 시스템 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20230123573A1 (en) * | 2021-10-18 | 2023-04-20 | Oracle International Corporation | Automatic detection of seasonal pattern instances and corresponding parameters in multi-seasonal time series |
-
2023
- 2023-07-12 CN CN202310857764.1A patent/CN116882572B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110348040A (zh) * | 2019-05-21 | 2019-10-18 | 合肥通用机械研究院有限公司 | 一种基于均方收敛条件的不均匀沉降函数拟合方法 |
KR102409155B1 (ko) * | 2021-10-08 | 2022-06-16 | 주식회사 다올 | Lstm 기반 지하수 수위 예측 시스템 |
Non-Patent Citations (5)
Title |
---|
基于傅里叶拟合的多种季节划分法及相似度量;戴婷;;湖南工程学院学报(自然科学版)(03);第34-40页 * |
基于小波组合模型的短期城市用水量预测;刘志壮;吕谋;周国升;;给水排水(10);第111-115,132页 * |
相对湿度及其变化的年循环研究进展;念达;邓琪敏;付遵涛;;地球科学进展(07);第102-114页 * |
雷州半岛年降水量序列分析及预测模型;邹仁爱;陈俊鸿;;水文(01);第55-59, 33页 * |
频谱分析法在挠力河流域年降水量预报中的应用;李平;卢文喜;王福林;;水文(04);第25-27,30页 * |
Also Published As
Publication number | Publication date |
---|---|
CN116882572A (zh) | 2023-10-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Paternina et al. | Identification of electromechanical oscillatory modes based on variational mode decomposition | |
O’dea et al. | An operational ocean forecast system incorporating NEMO and SST data assimilation for the tidally driven European North-West shelf | |
Yang et al. | Probabilistic short-term wind power forecast using componential sparse Bayesian learning | |
Donelan et al. | On the growth rate of wind‐generated waves | |
Li et al. | A model-based climatology analysis of wind power resources at 100-m height over the Bohai Sea and the Yellow Sea | |
Xu et al. | Wavelet analysis and nonparametric test for climate change in Tarim River Basin of Xinjiang during 1959–2006 | |
CN107423544B (zh) | 一种流域水文极值的非平稳性影响指标的确定方法及装置 | |
Zhang et al. | Multivariate multifractal detrended fluctuation analysis of 3D wind field signals | |
Pandhiani et al. | Time series forecasting using wavelet-least squares support vector machines and wavelet regression models for monthly stream flow data | |
Yan et al. | Changes in extreme wind speeds in NW Europe simulated by generalized linear models | |
CN104619009A (zh) | 一种定位数据采样周期的调整方法、装置及移动终端 | |
CN111651896B (zh) | 一种基于实际风速及实际复杂地形的风电场流场计算方法 | |
Deng et al. | The impact of inter-annual variability of annual cycle on long-term persistence of surface air temperature in long historical records | |
Reeve et al. | Statistical analysis and forecasts of long-term sandbank evolution at Great Yarmouth, UK | |
CN104268408A (zh) | 一种基于小波系数arma模型的能耗数据宏观预测方法 | |
Frison et al. | Chaos and predictability in ocean water levels | |
CN115239105A (zh) | 一种在役风电场的风资源的评估方法以及装置 | |
Liu et al. | On selection of the optimal data time interval for real-time hydrological forecasting | |
Zhu et al. | Validation of rainfall erosivity estimators for mainland China | |
CN116882572B (zh) | 一种降水年代际变化的预测方法、装置及存储介质 | |
CN113808015A (zh) | 输电线路区域气象参数的空间插值方法、装置及设备 | |
Hisaki | Inter-comparison of wave data obtained from single high-frequency radar, in situ observation, and model prediction | |
Zhang et al. | Multifractal analysis of measure representation of flood/drought grade series in the Yangtze Delta, China, during the past millennium and their fractal model simulation | |
Li et al. | Periodic identification of runoff in Hei river based on predictive extension method of eliminating the boundary effect of the wavelet transform | |
Xiao et al. | An improved model to estimate annual sand transport rate by sand-driving winds |
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 |