一个简化的潮汐预报准调和分析方法

1、 一个简化的潮汐预报准调和分析方法1，21，22 ，蒋风芝,李慧娟王如云（1 河海大学水文水资源与水利工程科学国家重点实验室，南京210098； 2 河海大学海洋学院，南京210098） 摘要：在用现有的浅水港日潮汐准调和分析预报方法进行潮汐分析预报时，发现最小二乘法的法8，算法不稳)，因此矩阵是坏条件的(或为病态的方程组的系数矩阵条件数很大，数量级在10定。根据潮汐动力学寻找高频潮族与低频潮族之间可能的相互作用关系，在只考虑相角的变化率情况下，建立了一个简化的浅水准调和分析模型。利用连云港的多年实测数据检验，简化的准调和分析模型相对于原准调和分析模型来讲，最小二乘法的法方程组系数矩阵条件数小

2、很多，因此简化后的模型计算更为稳定。在实测数据时间较长的情况下，简化前后的模型预报精度相当。但当实测数据较短时，简化前的原模型却没有传统的调和分析模型的预报结果精度高，而简化后的模型却能保持比传统的调和分析模型的预报结果有一定的改善。特别是简化后比简化前的模型计 。算时间减少了68% ：浅水潮汐；准调和分析；潮汐预报关键词 1 引言 在潮汐预报方面，一般采用调和分析方法，在深水区域此方法可以获得很好的预报效果，但在浅水区域尤其是河口区域，由于浅水潮汐的复杂性，采用此方法往往不能获得满意的效1，其结果不能令人满意。为此，杜德森后来又提出了一个直60个分潮果。例如杜德森提出的2，该方法虽然使高低潮

3、的预报精度有了提高，但把它应用到接对高低潮进行浅水改正的方法逐时潮位预报上则有许多困难和不便之处。在浅水区域由于非线性效应的加大，潮波往往产生畸变。此时，高频振动的作用必须予以充分考虑。为了提高浅水区域潮汐预报的精度，从调和分析方法来讲就必须增加高频的浅水分潮。在水深不太浅的区域，浅水分潮的振幅会随着阶数的增高而迅速减小，所以在一般港口采用较少数目的主要浅水分潮即可满足潮汐预报的要求。但在浅水区，常常需要考虑到六阶甚至更高阶的相互作用，才能满足潮汐预报的要求。上个世纪六十年代，一些潮汐学者试图通过扩充高频分潮的数目以使预报结果获得改进，如Zelter and 34曾将分潮的数目扩充到110Ro

4、ssiter and LennonCumimngs多个，但效果并不理想。方国以及洪等人认为，不理想的原因在于随着频率的增加，高频分潮的数目极速的增加，不可能从中挑选出少数分潮近似代替所有分潮，难以用有限数目的浅水分潮来体现总的浅水效应。可以认为通过增加浅水分潮以改进潮汐预报，其效果可能是比较有限的。 5基于以上分析，方国洪等提出了一个浅水潮汐预报的准调和分析方法，可以用来推算任意时刻的潮高，也可以用来推算高、低潮，效果比传统的调和分析法有了显著的改进。但我们使用此方法对连云港的多年潮位实测数据进行分析预报时，发现最小二乘法的法方程组系数矩阵条件数很大，算法不稳定。为此，我们对浅水准调和分析模型

5、进行了简化，简化后的模型计算更为稳定，计算时间大为减少。 2 准调和分析方法介绍 基金工程：水文水资源与水利工程科学国家重点实验室开放研究基金(2005407411)；中国教育部科学技术研究重点工程（104104）；中国江苏省普通高等学校高新技术产业发展工程(JH03-010) 作者简介：王如云(1963-)，教授，男，安徽芜湖人，从事计算物理学研究，E-mail:wangry 1 / 7 5提出的浅水预报准调和方法思路是把潮高分做两部分，一部分为低频部分，方国洪等人组成，主要是天文源潮波，另一部分为高频部分，属于浅水分潮。分析分两步21，由潮族0，中包含了11列出的分潮计算。表1，2，用如下

6、表进行，首先对于低频部分，即潮族0，?5)(t进行调和分析，求出各分潮的2的分潮。对实测潮汐Doodson所用的所有属于潮族0，1，调和常数。然后用实测水位减去平均水位和0，1，2族的潮位，依据剩余值再作进一步分析。 表1 调和分潮 序号 分潮符号 分析段数据时Doodson数 传统的调和分析 序号 模型分潮符号 1 模型2 Doodson数 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 j1 2 3 4 5 6 7 8 Sa SSa MmSM f Mf Q21? 1* QA1 Q1*

7、 QA1? 1* BO1* AO1 O1* OA1 PM1M 1? 1* KP21? 1P 1S 1K 1? 1? 1? 1J 1* Q2P1 QS1a 3rr?12r2r42r1r2r31?r1间长度（日）365 334 304 273 243 211 205 204 203 200 196 186 185 184 177 168 165 163 注：以上以1967 001000 002000 010100 022000 020000 130200 200031 211100 120100 121100 002121 211000 111000 000011 100011 200011 10

8、0000 001102 000114 001113 112000 100011 110000 111100 112000 211100 120001000413 213000 j j 10 1r 11 2?r12 113 1r? 14 215 2?r 16 22r17 2(cm) 均方差24.4 24.1 24.2 26.6 31.7 31.3 43.0 45.1 46.1 50.0 52.8 51.2 60.0 59.2 97.9 169.3 268.9 3.49.1 年1月1日起的潮位数据为分析资料，预报30 31 32 33 34 35 36 37 38 39 40 41 42 43 4

9、4 45 46 47 48 49 50 51 52 53 54 55 56 57 58 均方差(cm) 21.3 20.8 20.9 23.7 29.4 29.0 41.7 43.9 44.9 49.2 52.6 50.4 61.3 59.7 131.8 1433.8 6751.2 17626.9 OO 1* QS1* QK21 OQ2 SMN2 N22? 2* AN2 N2 NA*2? 2 OP2* AM2 M2* MA2 SMk2? 2L 2* BS2? 2S 2R 2* k2* kA2 NMS2KJ 2 M2S2 MSk2* N2S2a j2r?2r 213r 2r?3r 214r?r

10、124r 2r?4r 215r 26r 2均方差(cm) 21.8 21.3 21.5 24.2 29.6 29.3 41.7 43.8 44.8 48.9 51.7 50.1 59.2 58.3 97.4 169.1 268.8 349.0 1982年潮位 130000 001142 001140 010023 210023 020022 200022 100112 010021 110021 001122 000220 000120 200000 201000 202000 100221 001210 000422 001322 000222 110022 220000 221000 00

11、1232 001230 000424242000 254100 ?(t)，对其进行调和分析，计算各分潮的潮汐调和常数，再依据下假如利用一年潮位资料式将长周期、全日、半日潮族分别加以组合，得出三个基本准调和分潮，然后把三分日及以上的高频潮表示为这些基本准调和分潮的函数。随着频率的增加，只增加少量的准调和项，每一R(t)和2族的振幅，个准调和项可以近似看作一群频率相近的分潮之和。计算每小时01r(t)。各族的振幅和位相的变化已不再是常数，而具有缓慢的变化。称这些量是准调和分量。 2 / 7 5?gu)?t?(Rcosr?vfHcos? ii0ii0i1?i5?g)?(vsinRr?fHsinut

12、ii0ii0i1?i32?t?(v?u)cosRr?fHcosg （1） i1ii1iii?632?gv?ufHcos)t?(Rsinr? ii1i1ii6i?58?g?uHcos)t?(Rcosr?vf i2ii2ii33i?58?t?(v?cosu)?gRsinr?fH i2iii2ii?33从动力学的原因来看，浅水分潮由两种非线性产生。一种是运动方程中的平流项，如?u)?(?uu.，能够产生高级摄动项。另一种是由运动，和连续方程中的非线性项，如， ?x?x?uuu 是摩擦系数）。将上式定义是流速，（这里方程中的摩擦项所产生，如是水位，?)a?kAcos(表示由非线性产生的属于高频部分的各

13、阶摄动的准调和项视为单一潮波，并用项，则潮高可表达为： ?)a?kAcos(?h? （2） jjjj式中 58?g)?(t?v?h?A?Hfcosu ）（3 i0iiii1?i?,)cos(a?KKA为准调为浅水分潮部分，用(2)式中34个准调和项表示浅水效应，式中A,a为已知参数，它们的表达式见表2。预报时，首先由表1中第1到第58个分潮的和常数，rrRRhA,a。和计算按表2 式计算各族的式计算调和常数按(3)。再由(1)和，根据利用文1中的潮汐预报准调和分析方法对连云港多年的实测潮位资料进行了分析预报，得出的结论是改善效果明显，但浅水系数矩阵的条件数极大（见表4），这意味着方程组病态很严

14、重，该算法不稳定。在进行浅水准调和分析时，选择的浅水分潮并非越多越好，选择的分潮达到一定个数后，浅水系数矩阵奇异严重，预报精度不仅没有改善，反而会由于矩阵奇异给预报可靠度带来负面作用。实际上在观测记录数据资料的过程中，会有恶劣天气或仪器磨损等意外情况导致潮位数据缺测或者具有重大误差等情况发生及观测数据长度等的不同，算法若不稳3 / 7 定，这些原始误差就会在预报过程中被扩大，导致预报结果的可靠性下降，这对预报是极为不利的。鉴于这些，我们对原准调和分析方法进行了简化。根据潮汐动力学寻找高频潮族与低频潮族之间可能的相互作用关系，在只考虑相角的变化率情况下，建立了一个简化的浅水准调和分析模型。 A,

15、a表达式表2 浅水准调和项中 AaAajj jjjj223r?r3rRRR 1 18 11212133?rrRRr?rRR 2 19 11221221222r?2rr?rRRRR 20 3 2211221123rr2rRR?R 21 4 222112323?2rrrRRR 22 5 221221433r2r?rRRR 23 6 21222123r?3rrR4RR 24 7 2111213rR?3r2r?RrRR 25 8 22211121324r?rr2?rRRRR 9 26 12122121434r?2rrrRRRR 27 10 11222121324rr2RR 28 11 2222434r

16、2rRR 29 12 222242rr2?4rRRRR 30 13 2122120425rr?r3RRR 31 14 22121255rr?2RrRR 32 15 2112225262?rrrRRR 33 16 212221636?r2rrRRR 34 17 212221 3简化的准调和分析模型 在忽略了摩檫力作用的前提下，研究潮波在一维等深等宽半无限长沟渠中传播的情况，简?u?u?u?g?0? ?t?x?x化后的方程组中，第一个方程的第二项，第二个方程的第三和第四?u?u?u?0?h? ?t?x?x?x?项与相应方程的主要项相比，其量级即使在浅水地区也小于1，作为零级近似，把它们略去了。因此

17、倍潮波的振幅与相应源潮波振幅的整数次方成正比例，且幂次与频率的倍数相等；复合潮波的振幅则正比例于源潮波振幅的乘积，这个结果对由无摩檫引起的线性浅水分潮大致上是成立的，但是受摩檫力非线性作用引起的浅水潮波，它们的振幅与源潮波的振幅则不一定遵6从以上关系，实际的振幅变化幅度要小于理论的变幅。因为摩檫力比例于潮汐振幅的平方，1，不能认为引潮力增大若干倍，潮汐振幅也相应增大相同的倍数。由此我们而不是线性关系a?)?cos(Ka的含义同原来的模型），同时由于简假设浅水分潮的一般表达式为（表达式中4 / 7 A的计算，可节省计算时间。 2中浅水准调和项表达式化后的计算模型避免了对表简化模型求解浅水调和常数

18、时方法同原来的相似。最后进行预报并与原模型进行分析比a相同的项，我们只须保留一项即可，否则浅水系数矩阵的行列式由于存在完中较。对于表2全相同的两行或多行，其行列式为零，从而系数矩阵奇异，方程组无解。简化后的准调和项为a的表达式见表3。 17项，表4为简化前后模型的预报结果比较，仍以1967年实测潮位作为分析资料，来预报其它年份水位。为书写方便，我们称简化前模型为模型1，简化后模型为模型2。 a的表达式准调和项中 表3 rr3? 9 12 表4 模型1和模型2均方差比较表 调和分析均方模型1均方差模型2均方 年份差（） （） 差（）29.5 30.5 28.6 1962 22.6 24.2 21

19、.3 1963 23 21.7 24.5 1964 23.3 22.1 1965 24.8 23 24.6 21.7 1966 18.9 20.8 17.3 1967 22.9 24.2 21.6 1968 27.1 28.3 26.1 1969 24.7 26 1970 23.5 23.5 24.8 1971 22.1 26.4 1972 27.6 25.3 25.9 1973 27.1 24.8 27.3 28.5 1974 26.3 25.5 26.8 1975 24.6 22.6 1982 21.2 24.4 浅水系数矩阵条件数 模型1：35912244 模型2：14 分析表4，预报精

20、度虽没有所提高，但两者均方差相差仅为0.9cm到1.6cm之间，同时模型2的浅水系数矩阵条件数相比原来的算法显然小很多，说明简化后的算法稳定。算法若不稳定，原始误差就会在预报过程中被扩大，导致数据的可靠性下降，这对预报是极为不利的。如果算法稳定，即使用于分析的资料存在误差，在预报时，对这种误差也不会有很明显的放大，这对预报结果的可靠性是有利的。 为了进一步验证模型2稳定性好于模型1，我们改变用于分析的潮位资料的时间长度，比较两种模型的预报结果。 5 / 7 表5显示，当用于准调和分析的数据时间段长度缩减到204日时，模型2的均方差开始小于传统的调和分析方法及模型1。继续缩减到168日（约半年）

21、之后，模型1的预报均方差开始远远大于以传统的调和分析方法预报的均方差，以致大到难以接受。而模型2较传统方法的均方差有所改进，显然简化后的模型稳定性要好得多。 对简化前后模型进行浅水准调和分析计算（包括将各潮族的分潮进行迭加到求出浅水准调和常数的全部过程）所花的时间做了统计（表6）： 表5 不同分析时间长度两个模型预报精度比较 表6 简化前后模型计算时间比较 原模型计算时间改进后模型计算时间 年份（秒） （秒） 0.78 1962 2.48 0.77 2.41 1963 0.77 2.47 1964 0.80 1965 2.41 0.77 2.44 1966 0.78 1967 2.41 0.7

22、7 1968 2.39 0.78 1969 2.41 0.75 2.41 1970 0.75 1971 2.41 0.75 2.39 1972 0.77 2.39 1973 0.77 2.44 1974 0.77 1975 2.45 0.75 2.41 1982 0.77 平均计算时间 2.42 6 / 7 相对原来模型节省计算时间：68% 由表6可以看出，简化后的模型对浅水准调和分析部分计算时间仅有原来的68%左右，当分析数据量很大时，模型二的优势体现的就更为明显了。 5 结论 利用连云港潮位资料，对简化后的模型进行了分析检验，首先其法方程系数矩阵条件数为814，而简化前的模型高达10，说明简化后的模型其算法更为稳定，而原模型计算不稳定。其次，当分析时段潮位资料时间长度短于204日时，模型1的预报误差较传统调和方法来讲，不但没有减小，反而更大了，而模型2较传统方法来讲，仍然有所改进。由于简化后的模型避免A的计算，且经简化后准调和项只剩下17了对浅水准调和项表达式中项，从而在进行浅