戴荣-江三江营短期潮位统计预报研究.pptVIP

一、概述 长江下游在涨潮时,海水从长江入海口沿长江干流上溯倒流,潮流界一般在江苏省江阴市附近。潮流界以上受上游洪水影响为主,潮流界以下受长江入海口潮位影响为主。长江三江营位于江苏省江都市,在潮流界以上,属感潮河段,为半日潮形,即受潮汐周期变化影响，又受上游洪水影响。因此，采用门限回归模型建立长江三江营特征潮位与上游大通站径流和下游吴淞特征潮位的统计方程，根据前一日大通流量和相应的吴淞特征潮位预报三江营未来一个完整潮位过程的特征潮位和潮时，再利用不等距三次样条函数计算整点潮位，从而实现三江营短期潮位过程的预报。其中，大通站径流采用实测日平均流量资料，吴淞潮位采用天文潮调和分析进行预报 三江营地理位置概图 二、模型方法－总体设计框架 2.1 模型方法－吴淞特征潮位预报 吴淞口潮位过程属于天文潮,可以采用天文潮计算中的潮汐调和分析方法进行预报。潮汐调和分析中的潮高计算公式(2.1)为： 模型方法－吴淞特征潮位预报 (1)调和分析 式（2.2）表明，实际潮位可以看作是若干调和分潮叠加的结果。在此次研究中,采用的是固定分潮数为128的调和分析计算。取连续365天或369天的正点实测潮位资料,运用最小二乘法确定128个分潮的未知系数Xj及Yj,从这些系数,再确定出各分潮的平均振幅Hj及位相迟角gj,包括平均海面S0,这样便完成了调和分析。 (2)高、低潮位计算 在高、低潮计算中, t值须经过计算确定。由式(2.1),对时间取一次导数并设定为零,d(h(t))/dt=0,可得式（2.3) 模型方法－吴淞特征潮位预报 ………………………………………式(2.3) 求解式（2.3）得到一系列t值,代表高、低潮出现的时刻。将这些t值代入式(2.1),即可计算出相应的高、低潮潮位值。事实上,这个零导数方程不易求解,一般采用近似方法逐步逼近高、低潮出现的时刻t值。这些近似算法包括〝二分法〞、〝插值法〞、和〝二分–插值法〞等。本模型采用的是〝二分法〞算法。 2.2 模型方法－三江营特征潮位潮时预报 一、三江营短期特征潮位预报 门限回归模型基本思路就是依据某变量不同取值的范围，采用若干个线性回归模型来描述非线性关系。 门限回归模型一般形式 模型方法－三江营特征潮位潮时预报 设定门限 将吴淞特征潮位按0.5m划分，然后在每一个小区间内建立三江营特征潮位预报方程。 模型方法－三江营特征潮位潮时预报 二、三江营短期特征潮时预报 三江营特征潮位发生时间（潮时）预报采用经验方法，先统计吴淞特征潮位发生后传播到三江营的时间，预报时，以吴淞发生的特征潮位的时间加上吴淞至三江营的潮水传播时间即为预报的三江营特征潮位的发生时间（特征潮时） 高潮传播时间: 6小时40分 低潮传播时间: 7小时50分 2.3 模型方法－三江营短期潮位过程预报 三、 模型率定 根据2000－2003年的大通日平均流量、三江营特征潮位和吴淞特征潮位资料，建立以相应的吴淞特征潮位为门限变量的门限回归模型 采用最小二乘拟合方法进行模型参数率 模型率定 特征潮位预报方案曲线 四、精度评定 4.1 特征潮位及潮时 模拟精度评定 特征潮位及潮时 模拟精度评定 高高潮及低低潮率定过程线 4.2 特征潮位及潮时 预报精度评定 预报检验期：2004～2005年（共17个月） 特征潮位及潮时 预报精度评定 2004年三江营高高潮和低低潮的预报精度统计表 特征潮位及潮时 预报精度评定 高高潮及低低潮预报实测对比图 五、主要结论 本文根据长江三江营潮位主要受上游径流和下游天文潮 双重影响，以及不同量级天文潮汐对三江营潮位影响显著 的特点，建立了以下游吴淞特征潮位为门限变量、上游大通日平均流量为自变量的门限回归三江营特征潮位预报方案。 同时，考虑到目前实时预报时的资料限制，采用天文潮调和分析方法对吴淞特征潮位进行预报；并利用不等距三次样条函数插值整点潮位，得到三江营短期潮位过程预报。 根据200～2003年的潮位资料对模型进行了率定，并采用2004～2005年潮位资料进行了预报检验。结果表明，建立的潮位预报方案具有较高的精度。 本文提出的感潮河段潮位短期统计预报方法，可供类似研究和应用借鉴。 * * 长江三江营短期潮位统计预报研究 概述 模型方法 模型率定 精度评定 主要结论 大通站 三江营 吴淞 吴