长江三江营短期潮位统计预报的研究.pdfVIP

第四篇应用实践与行动措施 长江三江营短期潮位统计预报研究 梁忠民戴荣 (河海夫学水资源环境学院南京210098) 摘要呆用门限回归模型建立长江三江营特征潮位与上游大通站径流和下游吴凇蔫j位的坑计方程，根据前一目 大通流量和相应的矍凇特征潮位预报三江营来束一十完整潮位过程的特征潮位和潮时，再利用不等距三次样紊函 鼓计算整点嘲位，从而实现三江营短期潮位过程的预报。其中，大适站径浇栗用蜜涮日平均流量资料，吴耪潮位 采用无文潮调和分析进行预撮。模拟及jji报验证结果表明．所建立的模型具有较高曲精度。 关键词 门限回归模型特征潮位样条函数调和舟析 1 5l言 长江下游在涨潮时，海水从长江人海El沿长江干流上溯倒流，潮流界一般在江苏省江阴市附近。潮流界 以上受上游洪水影响为主，潮流界以下受长江人海口潮位影响为主。长江三江营位于江苏省江都市，在潮流 界以上，属感潮河段，为半日期形，即受潮汐周期变化影响，又受上游洪水影响。长江中下游的洪水过程变 化缓慢，而海洋潮汐运动对三江营潮位变化的影响比较固定，所以为采用统计分析预报三江营的潮位提供了 物理基础。上游来水的影响可以采用大通站为代表站，下游潮汐影响以吴凇口为代表站。分析表明．相同上 游来水情况下，三江营潮位与吴淞口的潮位具有良好的对应关系，但吴淞口不同量缴潮位下其规律不同，所 以，本丈采用门限回归分析建立以吴凇口潮位为仃限变量的三江营潮位预报模型。使用模型进行三江营特征 潮位预报时(高高潮、高低潮、低高潮、低低潮)，大通站径流量采用每日8时的流量代表，下游潮汐变量 采用吴淞口相应特征潮位。采用经验统计分析途径预报特征潮位的潮时，即根据历年吴淞口和三江营潮位资 料，分别统计确定吴淞13的特征潮位传播到三江营的时间，再根据吴凇口的相应潮时进行三江营潮时的预 报。根据预报的四个特征潮位和潮时，采用样条函数插值计算各整点潮位，得到三江营一个完整的潮位过程 检验，提供了精度评定结果。 2吴凇口特征潮位的推求 吴淞口潮傅过程属于天丈瑚，可以采用天文潮计算中的潮汐调和分析方法进行预报。潮汐调和分析中的 潮高计算公式为： J ^(f)=So+∑Z吩cos[ojt+(％+u)，一毋](1) 式中：h为l时刻的潮位；So为平均潮位；H，g为调和常数；，，“为交点因子和订正角；盯为分潮的角速 率；vo为t=0时刻的天文相角。 上式也可表示为： 、 ， ， ^(‘)=岛+王斗。∞咋l+王斗8in听‘ (2) 式中：吩《弓c∞【岛一(‰+u)，]，yj《Hjsin[毋一(％+n)，]。 式(2)表明．实际潮位可以看作是若干调和分潮叠加的结果。由于在实际分析中只能选取有限个较为 主要的分潮，所以，正确合理地选择分潮数是提高预报精度的关键。一般情况下，高频分潮能较好地描述地 形变化对潮波所含有的短周期波变形的影响．而低频分潮能有救地反映气象状况年过程对径流长波季节变化 的影响。调和分析就是确定一组分潮，使得各分潮叠加后能较好拟台浅海瑚汐过程。在本文中，采用的是固 定分潮数为128的调和分析计算。调和分析的步骤是取连续365天或369天的正点实测潮位资料，运用最小 第十九章长江流域 ·887· 二乘法确定128个分潮的未知系数耳及竹，从这些系数，再确定出各分潮的平均振幅q及位相迟角爵，包 括平均海面S，这样便完成了调和分析。 在调和分析的基础上，利用式(2)进行正点潮位预报，给定正点时『日lt值，即可计算相应的潮高。而 0，可得 』 ∑蛳只cos[4￡+(Vo+n)，一曲]=0 (3) 求解式(3)得到一系列￡值，代表高、低潮出现的时刻。将这些I值代人式(1)。即可计算出相应的高、 低潮潮位值。事实