[发明专利]一种高精度的河口短期潮位预报方法

权利要求书

1.一种高精度的河口短期潮位预报方法，其特征是，包括以下步骤：

S1，收集河口地区起报时刻前M+m时段内的径流站实测径流数据、分析站和参考站的实测潮位数据，根据收集的参考站潮位数据计算其对应潮差序列，径流站测得的或预报的径流序列用Q(t)表示，分析站潮位数据序列用h(t)表示，参考站测得的或预报的潮位序列计算而得的潮差序列用R(t)表示；

S2，收集河口地区起报时刻前m时段以及后K小时的径流站预报径流数据，其中，M时段为长时段，m时段为短时段；

S3，采用调和分析方法，基于参考站需预报时刻前m时段的前M时段的潮位数据以K小时为预报步长滚动预报其前m时段内的潮位值，并由此计算其对应的预报潮差序列；基于参考站起报时刻前M时段潮位数据预测后K小时内的潮位值并由此计算其对应潮差；

S4，采用非稳态调和分析方法，基于分析站需预报时刻前m时段的前M时段的实测潮位数据及同时段的实测径流站径流序列和参考站潮差序列，求解对应未知系数，结合步骤S2中前m时段的径流站径流预报值以及S3中同时段的参考站预报潮差序列，以K小时作为预报步长滚动预报分析站起报时刻前m时段内的潮位值；基于分析站起报时刻前M时段实测径流站径流序列、分析站潮位序列和参考站潮差序列进行分析，求解对应未知系数，结合步骤S2中后K小时的径流站径流预报值以及S3中后K小时参考站预报潮差，预报分析站起报时刻后K小时的潮位值；

S5，计算S4中前m时段的分析站预报潮位值与S1中分析站对应时段实测潮位数据的差值，作为分析站起报时刻前m时段的预测误差序列；

S6，采用自回归方法对S5算得的误差序列构建自回归模型；

S7，根据S6构建的自回归模型，预测起报时刻后K小时内非稳态调和分析方法的预报误差值；

S8，将S4计算的起报时刻后K小时的潮位值减去S7中预测的预报误差值，作为分析站起报时刻后K小时的最终预报潮位值。

2.根据权利要求1所述的一种高精度的河口短期潮位预报方法，其特征是，S1中M时段为8785小时，m时段取值范围为：360≤m≤1200，单位为小时，基于过去不同m时段的预报误差序列构建自回归模型，以进行K小时滚动预报时达到的最小预报均方误差的目标值为m确定值。

3.根据权利要求1所述的一种高精度的河口短期潮位预报方法，其特征是，S2中K小时取值为K≤48。

4.根据权利要求1所述的一种高精度的河口短期潮位预报方法，其特征是，所述调和分析方法潮位表达式为：

其中h(t)为分析站实测潮位；η₀为平均海平面；t为时间；k为分潮序号；n为调和分析采用的分潮数；σ_k、f_k、u_k、υ_0,k、H_k和g_k分别为第k个分潮的频率、相位角、交点因子、交点订正角、天文初相角、振幅和迟角，其中调和常数H_k和g_k的求解采用最小二乘法。