基于SARIMA模型的月径流量预测

1、基于SARIMA模型的月径流量预测基于SARIMA模型的月径流量预报3791202102c-0042-03数据驱动模型不以水文过程作为模型建立的基础，而是着重于数据关系的分析，建模方便高效，预报精度高。常见的数据驱动模型有多元回来分析、时间序列分析、神经网络等。其中多元回来分析预报因子选择困难，对预报精度影响较大，神经网络需要数据量大，预报结果不稳定。时间序列模型通过少量数据就可做出精确的短期预报，因此被广泛应用，其中SARIMA模型是一种考虑到序列周期改变的时间序列建模方法。某一观测或统计数值按其发生的时间先后顺序形成的数列称为时间序列1。水文观测值按期观测时间形成的数列是一种典型的时间序列

2、2。该文在水文时间序列分析的基础上建立SARIMA模型对某水文站历史月径流量进行拟合，选择最优模型对将来一段时间的月径流量进行预报，以期为该地区水文工作提供参考。1讨论方法对于平稳时间序列4，一般可对其建立自回来滑动平均ARMAp，q模型，它是自回来ARp模型和移动平均MAq模型的组合，其模型的基本形式如式1。yt=c+1yt-1+2yt-2+.+pyt-p+et+1et-1+2et-2+qet-p1式1中，yt为时间序列第t时刻的观看值;yt-1，yt-2，yt-p为时序yt的滞后序列;et，et-1，et-2，et-q为模型在第t期，第t-1期，第t-q期的误差;1，2，p，1，2，q为待

3、估量参数;c为常数项。ARIMAp，d，q模型是ARMA模型的变换形式，其中d表示原始序列经过d次差分后转化为平稳时间序列。某些时间序列中，存在明显的周期性改变，这些周期是由于季节性改变包括季度、月度等改变或一些外部因素引起的，这类序列称为季节性序列。对时间序列经过有限次差分和季节差分，将随机误差的长期影响变成临时影响，将其转化为平稳时间序列，在此基础上建立季节性差分自回来滑动平均SARIMAp，d，q×P，D，Qs模型。考虑到水文站月径流量的影响因素众多，流域下垫面条件冗杂，建立过程驱动模型所需数据较多，因此在其径流特性分析的基础上考虑对其建立时间序列模型并对其径流量进行预报。2实

4、例分析该文收集了某水文站1987年1月到2021年6月共计378个月的径流量资料，使用1987年1月到2021年12月共30年360个月的流量数据作为训练数据，使用2021年1月到2021年6月共18个月的流量数据作为校验数据，通过对训练数据进行分析处理并建立SARIMA模型。最终把预报值和期望值进行比较并评定预报精度。2.1试验数据预处理原始数据存在明显的以12個月为周期的波动性并且不同时刻波动浮动差异较明显，不满足时间序列建模的必要条件。所以，考虑对原始数据进行对数一阶12次差分改变，以减小原始数据的波动性和趋势性。绘制变换后序列的自相关和偏相关图，从自相关图和偏自相关图观看该序列近似为一

5、个平稳过程且为非白噪声序列，满足时间序列建模的要求。2.2模型识别和参数估量依据自相关图和偏相关图的拖尾性和截尾性进行初步的模式识别。在初步识别的基础上，建立可能的模型，通过最小二乘法估量其模型参数并计算统计量P值如表1所示，P值越小表示模型参数越显著，其中红色标出的参数的P值不显著为零一般当P值小于0.05则认为该参数显著，其中模型6，8的全部参数都显著。在全部参数都显著的模型中，依据AIC值和SBC值进行模型优选，AIC和SBC值越小代表模型的效果越优良。在模型6、8中，模型8的AIC值和SBC值均最小，所以模型8为所求的最优模型，其模型表达式为：yt=0.5208yt-1+0.4625y

6、t-2-0.8253yt-12-0.9994et-2-0.9155et-2422.3模型检验对模型拟合值残差进行白噪声检验，残差序列的自相关系数和偏自相关系数基本都在二倍置信区间内，通过计算残差序列的P值全部大于0.05，没有拒绝假设性试验，相关系数与基本与零没有显著差异，该序列可被看作是白噪声序列，说明模型拟合效果较好，信息提取较充分。2.4模型预报模型建立完毕后，将训练数据输入模型。除个别峰值处整体拟合效果比较理想，对2021年1月到2021年6月的径流量进行动态预报，预报结果见图1。2.5预报精度评定对水文预报方案的有效性评定接受以下确定性系数dy进行。dy越大，方案的有效性越高。345式中，Se为预报的均方差;y为预报要素值的均方差;yi为实测值;y为预报值;为实测系列的均值;n为实测系列的点据数。评定方案的有效性时按表4中标准进行。预报方案有效性到达甲dy0.90、乙dy处于0.700.90时两个等级时，方案可被用于作业预报;方案等级为丙等dy处于0.50.69时时，只可用作参考性预报;方案有效性低于丙等时不能用于作业预报，只能做参考性评估。计算此次预报确实定性系数，得dy=0.8594，预报方案的有效性为乙等，说明该次预报效果较好，可用于作业预报。3结语该文在对某水文站30年历年逐月径流量序列分析的基础上，建立了SARIMA模型对历史径流量进行拟合。利用建