大港驗(yàn)潮站潮汐分析與國家高程基準(zhǔn)面變化

吳富梅，魏子卿，李迎春

1.地理信息工程國家重點(diǎn)實(shí)驗(yàn)室，陜西 西安 710054；2.西安測(cè)繪研究所，陜西 西安 710054

1 引 言

高程基準(zhǔn)是高程測(cè)量的起算依據(jù)，是大地測(cè)量基準(zhǔn)的重要組成部分［1－2］。青島大港驗(yàn)潮站處于東部黃海沿岸，在我國海岸線的中部，1957年確定為我國基本驗(yàn)潮站。世界各地區(qū)的高程基準(zhǔn)一般通過沿岸驗(yàn)潮站的平均海面高實(shí)現(xiàn)，我國將大港驗(yàn)潮站的平均海水面定義為我國高程基準(zhǔn)的起算面。我國目前采用的“1985國家高程基準(zhǔn)”，是由大港驗(yàn)潮站1952—1979年驗(yàn)潮數(shù)據(jù)確定的。計(jì)算時(shí)以19年為周期，滑動(dòng)步長(zhǎng)為1年，得到10組平均海面值，取其平均值242.89cm作為當(dāng)?shù)攸S海平均海面高，然后將這個(gè)平均海面作為我國的高程基準(zhǔn)面，從而求出國家水準(zhǔn)原點(diǎn)的高程值為72.260m［2－5］。

海水受到月球和太陽等天體的引潮力作用產(chǎn)生規(guī)律性的海洋潮汐運(yùn)動(dòng)，研究驗(yàn)潮站處的潮汐變化規(guī)律對(duì)于確定國家高程基準(zhǔn)具有理論和實(shí)際意義［6－7］。18世紀(jì)80年代達(dá)爾文對(duì)引潮位進(jìn)行了近似調(diào)和展開，得到了60多個(gè)主要分潮的頻率和平均振幅［4］。19世紀(jì)20年代杜德森引用月球運(yùn)動(dòng)的Brown系數(shù)和Newcomb表對(duì)引潮位進(jìn)行了純調(diào)和展開，獲得400多個(gè)分潮的頻率和振幅［8－10］。文獻(xiàn)［11—12］采用更新的天文變量按照杜德森方法將引潮位展開至3000多個(gè)分潮。這些分潮對(duì)解釋潮汐現(xiàn)象起到了很大作用，但是這些展開的分潮振幅和相位并不能表征實(shí)際海域的分潮特征，因此需要采用大量潮汐觀測(cè)數(shù)據(jù)進(jìn)行潮汐調(diào)和分析來解釋更符合實(shí)際的潮汐現(xiàn)象、分析平均海水面變化以及進(jìn)行潮汐預(yù)報(bào)。大港驗(yàn)潮站自1952年至今已有60多年的觀測(cè)數(shù)據(jù)，這對(duì)進(jìn)一步分析我國高程基準(zhǔn)面的變化具有重要作用。

另外，溫室效應(yīng)造成的冰川消融、局部地區(qū)地下水的開采以及海底地形的變遷等，都會(huì)引起海水面的升降變化［13］。實(shí)際計(jì)算1985國家高程基準(zhǔn)時(shí)，受當(dāng)時(shí)數(shù)據(jù)限制，并未考慮這種變化。隨著大地測(cè)量精度的提高，國家高程基準(zhǔn)如同坐標(biāo)基準(zhǔn)一樣，應(yīng)該考慮各種影響因素，以便保持現(xiàn)勢(shì)性［14－16］，因此對(duì)國家高程基準(zhǔn)面最近60年的變化進(jìn)行研究也有積極的意義。

為了分析和研究大港驗(yàn)潮站處潮汐的特征以及國家高程基準(zhǔn)面在60年中的變化，本文首先采用1980—2011年每小時(shí)潮汐數(shù)據(jù)和1952—2011年月平均數(shù)據(jù)進(jìn)行頻譜分析獲得影響大港驗(yàn)潮站處平均海水面的主要分潮；然后采用最小二乘調(diào)和分析對(duì)每小時(shí)潮汐數(shù)據(jù)進(jìn)行逐年分析獲得主要分潮的振幅和變化；并采用調(diào)和分析和18.61年移動(dòng)平均對(duì)60年間國家高程基準(zhǔn)面的變化進(jìn)行計(jì)算和分析；最后給出關(guān)于我國國家高程基準(zhǔn)定義的建議。

2 采用的數(shù)據(jù)

數(shù)據(jù)1：大港驗(yàn)潮站1980—2011年共計(jì)11 688d的每小時(shí)潮位數(shù)據(jù)，對(duì)其中缺少的72h數(shù)據(jù)（1985－06－01、2001－04－02、2002－03－02各24h），首先采用間隔25h和50h的數(shù)據(jù)平滑內(nèi)插進(jìn)行填補(bǔ)，之后采用3節(jié)調(diào)和分析的結(jié)果進(jìn)行填補(bǔ)。

數(shù)據(jù)2：大港驗(yàn)潮站（1952—2007年）月均潮位數(shù)據(jù)，每月一個(gè)數(shù)據(jù)，共計(jì)672個(gè)數(shù)據(jù)。

對(duì)這兩組數(shù)據(jù)共有年份（1980—2007年）的年平均數(shù)據(jù)進(jìn)行比較（如圖1，藍(lán)色為第1組數(shù)據(jù)，紅色為第2組數(shù)據(jù)），發(fā)現(xiàn)這兩組數(shù)據(jù)在1995年之前相差較小，最大差距為0.34cm；1995年之后相差較大，最大差距在2006年，為3.45cm?？紤]到每小時(shí)潮位資料的原始性和真實(shí)性，認(rèn)為第1組數(shù)據(jù)是可靠的，因此用第1組數(shù)據(jù)代替第2組數(shù)據(jù)中從1980年之后部分的月平均和年平均數(shù)據(jù)，從而獲得1952—2011年的月平均和年平均潮位數(shù)據(jù)（圖2）。

圖1 兩組數(shù)據(jù)的比較Fig.1 Comparison of two sets of data

圖2 1952—2011年年平均潮位數(shù)據(jù)Fig.2 Mean tidal data from 1952to 2011

3 驗(yàn)潮數(shù)據(jù)頻譜分析

3.1 小時(shí)潮汐數(shù)據(jù)的頻譜分析

根據(jù)奈魁斯特定律，小時(shí)潮汐數(shù)據(jù)頻譜分析可以獲得最高頻率為2h的高頻分潮。頻譜分析的結(jié)果如圖3所示。

圖3 每小時(shí)潮汐數(shù)據(jù)頻譜圖Fig.3 Energy－frequency for hourly tidal data

圖3中可見6處分潮集中頻段。第1處是長(zhǎng)周期分潮，主要表征周期大于半個(gè)月的分潮，其中可清晰判別年周期分潮Sa，半年周期分潮Ssa，月周期分潮Msm、Mm等；第2處是日周期分潮，主要表征周期在24h左右的分潮，比較有代表性的有K1、O1、P1、Q1、S1、ψ1等；第3處是半日分潮，如M2、S2、N2、K2、T2、R2等；第4處是8h左右分潮，如MO3、MK3等；第5處是6h左右分潮，如M4、MS4、S4等；第6處是4h分潮，如2MN6、M6、2MS6等。比較明顯的分潮遠(yuǎn)不止這些，將其中振幅較大（大于數(shù)值0.45cm）、周期在一年內(nèi)的分潮選出，共計(jì)有180個(gè)，基本包含第3節(jié)調(diào)和分析的170個(gè)分潮。受篇幅限制，各分潮不詳細(xì)列出。

3.2 月平均和年平均潮汐數(shù)據(jù)頻譜分析

理論上，海洋潮汐除了包括上述年周期以下的分潮外，還包含更長(zhǎng)周期的分潮，如交點(diǎn)潮（周期18.61年）、太陽黑子潮（周期11.13年）、近點(diǎn)潮（周期8.85年）、3.57年周期潮汐和極潮（周期1.19年）等。2.1節(jié)中的小時(shí)潮汐數(shù)據(jù)長(zhǎng)度較短，只有32年，對(duì)于交點(diǎn)潮不到兩個(gè)周期，對(duì)于太陽黑子潮不到3個(gè)周期，不適合長(zhǎng)周期潮汐分析，因此采用合并之后的60年月平均潮汐數(shù)據(jù)進(jìn)行頻譜分析，結(jié)果如圖4所示。

從圖4可以清晰地看出兩處峰值，第1處是年周期分潮，第2處是半年周期分潮。由于離散化原因?qū)е骂l譜分析獲得的周期與相應(yīng)的天文分潮周期不一致，如交點(diǎn)潮周期實(shí)際為18.61年，但分析出來為20年，見圖5。

實(shí)際上，以30d或31d為周期平均是不合理的，因?yàn)橛绊懽畲蟮脑路殖敝芷谑?7.554個(gè)太陽日，因此又利用年平均潮汐數(shù)據(jù)進(jìn)行頻譜分析，如圖5所示。表1給出主要分潮的周期和振幅。

圖4 月平均潮汐數(shù)據(jù)頻譜圖Fig.4 Energy－frequency for monthly tidal data

圖5 年平均潮汐數(shù)據(jù)頻譜圖Fig.5 Energy－frequency for yearly tidal data

表1 年平均潮汐數(shù)據(jù)主要分潮Tab.1 Main tides of yearly tidal data

如果認(rèn)為分析得到的潮波與天文分潮接近，就認(rèn)為海面高受這種天文分潮影響，那么從圖4和表2可看出，大港驗(yàn)潮站處海水面受到交點(diǎn)潮的顯著影響。綜合表1和表2，大港驗(yàn)潮站海水面受到交點(diǎn)潮、太陽黑子潮、6.667年、4.615年等的影響，近點(diǎn)潮和3.57年分潮影響不明顯，其中6.667年和4.615年分潮沒有找到與之對(duì)應(yīng)的分潮。

對(duì)月平均數(shù)據(jù)進(jìn)行自相關(guān)分析，如圖6所示，受年周期分潮影響，其自相關(guān)值表現(xiàn)出很強(qiáng)的年周期特性，其他分潮影響被掩蓋。如果認(rèn)為在整年處自相關(guān)值受周年分潮影響相同，那么自相關(guān)值表現(xiàn)出來的不一致性可以認(rèn)為是由附近的分潮引起的，為此給出整年處的自相關(guān)值，如圖7所示。從圖7可以清晰地看出，在19年處自相關(guān)值最大，可以認(rèn)為是由交點(diǎn)潮引起的［17－19］，另外在5年、7年、12年、15年和23年處自相關(guān)值也較大。5年、7年、12年與前面分析獲得的4.615年、6.667年、12年一致，同樣這里也沒有分析得到近點(diǎn)潮。另外，尚不清楚15年和23年的影響原因，需要累積更多的數(shù)據(jù)或者尋找更合適的方法來分析。

圖6 月平均潮汐數(shù)據(jù)頻自相關(guān)圖Fig.6 Autocorrelation for monthly tidal data

圖7 月平均潮汐數(shù)據(jù)頻自相關(guān)圖（整年部分）Fig.7Autocorrelation for monthly tidal data（full year）

4 驗(yàn)潮數(shù)據(jù)調(diào)和分析

潮汐的引潮位理論只能給出海洋潮汐變化的基本規(guī)律和特點(diǎn)，并不能完全解釋潮汐的各種現(xiàn)象。通過引潮位推導(dǎo)獲得的潮汐振幅與某一海域潮汐的真實(shí)振幅也不一樣，若想準(zhǔn)確了解其大小及變化規(guī)律，必須通過實(shí)際觀測(cè)進(jìn)行潮汐分析和計(jì)算。潮汐調(diào)和分析可以根據(jù)潮汐觀測(cè)數(shù)據(jù)計(jì)算各個(gè)分潮的調(diào)和常數(shù)，以了解該海區(qū)的潮汐特點(diǎn)、分析平均海水面變化機(jī)制，同時(shí)還可進(jìn)行潮汐預(yù)報(bào)。目前，潮汐調(diào)和分析最常用的方法是最小二乘估計(jì)法［8，20］。潮汐調(diào)和分析既適用于短時(shí)段的潮汐數(shù)據(jù)分析也適用于長(zhǎng)時(shí)段數(shù)據(jù)分析。短時(shí)段分析是一種準(zhǔn)調(diào)和分析，估計(jì)的分潮種類少、精度低，不需要考慮海面升降速率（每年1～2mm的升降速率在短時(shí)間內(nèi)可以忽略）；長(zhǎng)時(shí)段調(diào)和分析，估計(jì)的分潮種類多，待估參數(shù)也多，例如19年潮汐數(shù)據(jù)參數(shù)多達(dá)800多個(gè)至幾千個(gè)，還需要考慮海面升降速率。這里對(duì)1980—2011年的潮汐數(shù)據(jù)逐年進(jìn)行年潮汐調(diào)和分析，不考慮海面升降速率，可計(jì)算170個(gè)分潮［8］，基本包含2.2節(jié)分析得到的大部分分潮。

設(shè)大港驗(yàn)潮站t時(shí)刻的潮高為

式中，A0為年平均海水面高；Rj、σj、θj分別為第j個(gè)分潮的振幅、頻率和初相角；x（t）為隨機(jī)誤差。

式中，A0、aj、bj作為未知參數(shù)，這樣每年大約有8760（365d×24h）個(gè)觀測(cè)量，341個(gè)待估參數(shù)，在計(jì)算機(jī)還不成熟的時(shí)代，一般將a、b系數(shù)分開解算，并且進(jìn)行多次迭代，以節(jié)約計(jì)算時(shí)間［5］，現(xiàn)在直接通過矩陣計(jì)算這些數(shù)據(jù)已不成問題。對(duì)缺少的72h數(shù)據(jù)，通過5次迭代分析計(jì)算，進(jìn)行了更科學(xué)的修補(bǔ)。

為顯示驗(yàn)潮站處海水面的變化規(guī)律，作出了主要分潮在這32年間的振幅變化圖，圖8是Q1、O1、K1、K2的振幅變化，圖9是M2振幅的變化，圖10是Sa、N2、S2振幅的變化，圖11是P1、M4、MS4、M6振幅的變化?？梢钥闯觯琈2是影響大港海面高變化的主要分潮，其振幅最大達(dá)136cm，其次是S2、N2、K1、Sa等，其振幅分別達(dá)到40 cm、25cm、25cm、20cm等。其他沒有列出的分潮振幅很小，在1cm量級(jí)甚至更小。

從圖8可看出，Q1、O1、K1、K2的振幅具有非常明顯的約19年周期變化；M2盡管具有19年周期變化，但高點(diǎn)與低點(diǎn)的振幅不一致；其他分潮的周期變化不顯著。應(yīng)該指出，一年潮汐數(shù)據(jù)的最小二乘估計(jì)并不是完全的調(diào)和分析，估計(jì)獲得的分潮的振幅和相位并不僅代表本身的振幅，還受一部分長(zhǎng)周期分潮的影響，因此具有長(zhǎng)周期變化。在潮汐學(xué)中，還需要通過交點(diǎn)因子對(duì)各個(gè)分潮的振幅進(jìn)行訂正獲得分潮的平均振幅，來分析海域的潮汐變化，但在這里意義并不大。

圖8 Q1、O1、K1、K2 振幅的變化Fig.8 Changes of amplitudes of Q1、O1、K1、K2

圖9 M2振幅的變化Fig.9 Changes of amplitudes of M2

圖10 Sa、N2、S2 振幅的變化Fig.10 Changes of amplitudes of Sa、N2、S2

式（2）中的A0代表大港從1980—2011年的年平均海面。為了檢驗(yàn)這一結(jié)果的正確性，將該算法獲得的年平均海面高與直接平均獲得的年平均海面高進(jìn)行比較，如圖12所示。最大差值在1.5mm，平均值是0.071mm，標(biāo)準(zhǔn)差是0.66 mm，可知通過最小二乘調(diào)和分析獲得的年平均海面高是可信的。

圖11 P1、M4、MS4、M6 振幅的變化Fig.11 Changes of amplitudes of P1、M4、MS4、M6

圖12 A0與直接平均獲得的年平均海面高的差值Fig.12 Difference of A0and mean sea level

5 國家高程基準(zhǔn)面的變化

國家高程基準(zhǔn)由青島大港多年平均海水面定義。溫室效應(yīng)造成的冰川融化、局部海底地形變遷等使得國家高程基準(zhǔn)面產(chǎn)生變化。由圖2可看出，大體說來，1952—1980年平均海水面處于下降期，1980—2011年處于上升期。

國家高程基準(zhǔn)面的變化包含兩部分，一部分是長(zhǎng)期性變化，另一部分是周期性變化。研究高程基準(zhǔn)主要關(guān)心平均海面高和長(zhǎng)期變化，應(yīng)盡可能減小周期性變化的影響。下面通過兩種方法研究高程基準(zhǔn)面的變化。

5.1 最小二乘調(diào)和分析

以年平均海面高為基礎(chǔ)數(shù)據(jù)，在估計(jì)升降變化率k時(shí)，同時(shí)考慮交點(diǎn)潮（周期18.61年）、太陽黑子潮（周期11.13年）、近點(diǎn)潮（周期8.85年）、極潮（周期1.19年）以及具有3.57年周期的分潮。這里沒有考慮2.2節(jié)中分析出的6.67年、4.62年、15年和23年周期分潮的影響，因?yàn)檫@些分潮尚沒有找到與之對(duì)應(yīng)的天文分潮，它們?cè)陬l譜分析中的形成機(jī)制并不清楚，并且分析出的周期也不一定正確。

這樣，年平均海面高可表達(dá)成［21］

式中，A′0為t0時(shí)的平均海水面高。

圖13給出1952—1980年擬合的曲線與實(shí)際海水面圖，圖14給出1980—2007年擬合的曲線與實(shí)際海水面圖，表2給出這兩個(gè)時(shí)段分別計(jì)算的平均海面高、海水面上升速率以及各分潮的振幅。

從表2中可知，在1952—1980年期間太陽黑子潮、3.57年分潮和極潮對(duì)大港驗(yàn)潮站海面影響較大，在1980—2011年期間交點(diǎn)潮、太陽黑子潮、極潮對(duì)其影響較大。大港驗(yàn)潮站海水面在1952—1980年下降速率為1.07mm／a，在1980—2007年期間上升速率為1.59mm／a。

圖13 1952—1980年擬合曲線與實(shí)際海水面Fig.13 Fitting curving line and actual sea level for 1952to 1980

圖14 1980—2011年擬合曲線與實(shí)際海水面Fig.14 Fitting curving line and actual sea level for 1980to 2011

表2 1952—1980年間和1980—2011年間的平均海面、海面上升速率以及各分潮振幅Tab.2 Mean sea level，ascending rates and tides amplitudes for 1952to 1980and 1980to 2011

5.2 18.61年移動(dòng)平均

在4.1節(jié)中，僅考慮了5種分潮的影響，這是不符合實(shí)際情況的，剩余的分潮勢(shì)必會(huì)影響海面升降速率和t0時(shí)的平均海水面高的計(jì)算。在計(jì)算1985國家高程基準(zhǔn)時(shí)采用的是19年移動(dòng)平均法，這里考慮采用交點(diǎn)潮周期18.61年（224個(gè)月）移動(dòng)平均估計(jì)國家高程基準(zhǔn)面的變化，其主要目的是盡最大可能來削弱交點(diǎn)潮的影響，同時(shí)減小周期小于18.61年的分潮影響。

這樣，在19年中（228個(gè)月）必然有4個(gè)月的數(shù)據(jù)不需要使用，為了充分利用數(shù)據(jù)，在19年中以224個(gè)月為周期進(jìn)行移動(dòng)平均，獲得5個(gè)均值，再取其均值作為這19年的平均海面。將1952—1980年和1980—2011年的數(shù)據(jù)分別進(jìn)行處理，在1952—1980年期間獲得11個(gè)移動(dòng)平均值，1980—2011年期間獲得14個(gè)移動(dòng)平均值，然后以這些移動(dòng)平均值為觀測(cè)值，進(jìn)行一元線性回歸，獲得這兩個(gè)時(shí)段的海面升降速率以及對(duì)應(yīng)歷元的平均海面高，見圖15和圖16及表3。

圖15 1952—1980年擬合直線與平均海面Fig.15 Fitting line and actual sea level for 1952to 1980

圖16 1980—2011年擬合直線與平均海面Fig.16 Fitting line and actual sea level for 1980to 2011

表3 1952—1980年間和1980—2011年間的平均海面、海面上升速率Tab.3 Mean sea level，ascending rates for 1952to 1980 and 1980to 2011

由18.61年移動(dòng)平均獲得大港驗(yàn)潮站海水面在1952—1980年期間下降速率為0.76mm／a，在1980—2007年期間上升速率為1.62mm／a。與最小二乘調(diào)和分析相比，在1952—1980年差異較大，在1980—2007年差異不大。

5.3 比較與分析

從以上分析可以看出：

（1）最小二乘調(diào)和分析和18.61年移動(dòng)平均都能夠削弱一些主要分潮的影響，最小二乘調(diào)和分析能夠削除交點(diǎn)潮、太陽黑子潮、近點(diǎn)潮、極潮和3.57年周期分潮的影響，也能削弱部分其他分潮的影響，18.61年移動(dòng)平均能夠消除交點(diǎn)潮和以18.61年為公倍數(shù)的分潮影響（如半日潮、日潮等），也能夠削弱其他分潮的影響，但這兩種方法都不能消除所有分潮的影響。

（2）最小二乘調(diào)和分析法可以同時(shí)估計(jì)出海面升降以及某些分潮的振幅，但在同時(shí)估算時(shí)，分潮的周期性變化會(huì)影響升降速率的計(jì)算。

（3）18.61年移動(dòng)平均采用先平均削弱分潮影響后估計(jì)升降速率，分潮對(duì)速率計(jì)算影響較小，但這種方法會(huì)重復(fù)使用數(shù)據(jù)，一旦某一時(shí)段數(shù)據(jù)出現(xiàn)問題會(huì)放大影響面，如1952—1956年之間的數(shù)據(jù)觀測(cè)精度較低且是經(jīng)過修補(bǔ)的，是否存在問題是值得懷疑的，在最小二乘調(diào)和分析中這段數(shù)據(jù)只使用了1次，對(duì)整個(gè)計(jì)算影響較小，但在18.61年移動(dòng)平均中，這段數(shù)據(jù)影響了半數(shù)左右的平均海面的計(jì)算，這也是在1952—1980年兩種方法估計(jì)速率差異較大的部分原因。

（4）相比而言，在潮位數(shù)據(jù)不存在問題的前提下，例如1980—2011年，18.61年移動(dòng)平均比最小二乘調(diào)和分析法要穩(wěn)定、可靠。

6 1985高程基準(zhǔn)的檢核

利用18.61年移動(dòng)平均對(duì)1985國家高程基準(zhǔn)加以檢核是有意義的。定義1985國家高程基準(zhǔn)時(shí)，沒有考慮海面升降的因素，采用1952—1979年10組19年周期的移動(dòng)平均值的均值242.89cm作為基準(zhǔn)。這里筆者利用18.61年移動(dòng)平均對(duì)1952—1979年潮位數(shù)據(jù)進(jìn)行處理，得到平均海面高為243.03cm，與1985國家高程基準(zhǔn)相差0.14cm。這說明18.61年移動(dòng)平均相比于19年移動(dòng)平均是有差別的，但差別不大，也說明18.61年移動(dòng)平均用于計(jì)算國家高程基準(zhǔn)是可行的。

7 結(jié) 論

本文利用青島大港1980—2011年每小時(shí)潮汐數(shù)據(jù)和1952—2007年月平均數(shù)據(jù)對(duì)大港驗(yàn)潮站海洋潮汐進(jìn)行分析、對(duì)國家高程基準(zhǔn)進(jìn)行研究，得出如下結(jié)論：

（1）通過每小時(shí)潮汐數(shù)據(jù)、月平均以及年平均潮汐數(shù)據(jù)的頻譜分析，獲知大港海水面主要受到6個(gè)長(zhǎng)周期和180個(gè)較短周期分潮的影響，其中包括交點(diǎn)潮、太陽黑子潮、近點(diǎn)潮、極潮、年周期分潮Sa，半年周期分潮Ssa，月周期分潮 Msm、Mm 等，日周期分潮K1、O1、P1、Q1、S1、ψ1等，半日分潮M2、S2、N2、K2、T2、R2，8h 分潮 MO3、MK3等，6h 分 潮M4、MS4、S4等，4h 分 潮2MN6、M6、2MS6等。

（2）通過對(duì)每小時(shí)潮汐數(shù)據(jù)的調(diào)和分析，獲知M2是影響大港海面高變化的主要分潮，其振幅達(dá)到130cm，其次是S2、N2、K1、Sa等，其振幅分別達(dá)到40 cm、25cm、25cm、20cm 等，并且M2、Q1、O1、K1、K2的振幅具有非常明顯的19年周期變化。

（3）通過數(shù)據(jù)分析發(fā)現(xiàn)，國家高程基準(zhǔn)面在1952—1980年呈下降趨勢(shì)，在1980—2011年呈上升趨勢(shì)；利用最小二乘調(diào)和分析和18.61年移動(dòng)平均分別計(jì)算出相應(yīng)歷元的平均海面高和升降速率，并對(duì)這兩種方法進(jìn)行了比較。

（4）利用18.61年移動(dòng)平均對(duì)1952—1979年潮位數(shù)據(jù)進(jìn)行處理，得到平均海面高為243.03cm，與1985國家高程基準(zhǔn)相差0.14cm。

［1］ GUO Hairong，JIAO Wenhai，YANG Yuanxi，et al.Systematic Error of the 1985National Height Datum［J］.Geomatics and Information Science of Wuhan University，2004，29（8）：715－719.（郭海榮，焦文海，楊元喜，等.1985國家高程基準(zhǔn)的系統(tǒng)差［J］.武漢大學(xué)學(xué)報(bào)：信息科學(xué)版，2004，29（8）：715－719.）

［2］ YANG Y X，WEN Y N，XIONG J，et al.Robust Estimation for a Dynamical Model of Sea Surface［J］.Survey Review，1999，35（271）：2－10.

［3］ CHEN Zongyong.The Method of Calculating the Mean Sea Level［J］.Journal of Ocean Institute of Shandong，1960（1）：65－73.（陳宗鏞.計(jì)算平均海面的方法［J］.山東海洋學(xué)院學(xué)報(bào)，1960（1）：65－73.）

［4］ CHEN Zongyong.Tide［M］.Beijing：Science Press，1980.（陳宗鏞.潮汐學(xué)［M］.北京：科學(xué)出版社，1980.）

［5］ CHEN Zongyong，ZHOU Tianhua，YU Yifa，et al.A Study on the 1985National Datum Level of Altitude［J］.Journal of Ocean University of Qingdao，1988，18（1）：9－14.（陳宗鏞，周天華，于宜法，等.1985國家高程基準(zhǔn)的研究［J］.青島海洋大學(xué)學(xué)報(bào)，1988，18（1）：9－14.）

［6］ WEI Ziqing.On the Permanent Tide in Geodesy［J］.Acta Geodaetica et Cartographica Sinica，2012，41（4）：475－482.（魏子卿.大地測(cè)量中的永久潮汐問題［J］.測(cè)繪學(xué)報(bào)，2012，41（4）：475－482.）

［7］ CHEN Junyong.Permanent Tides and Geodetic Datum［J］.Acta Geodaetica et Cartographica Sinica，2000，29（1）：12－16.（陳俊勇.永久性潮汐與大地測(cè)量基準(zhǔn)［J］.測(cè)繪學(xué)報(bào)，2000，29（1）：12－16.）

［8］ HUANG Zuke，HUANG Lei.Tidal Principles and Calculation［M］.Beijing：Chinese Ocean University Press，2005.（黃祖柯，黃磊.潮汐原理與計(jì)算［M］.北京：中國海洋大學(xué)出版社，2005.）

［9］ DOODSON A T.The Harmonic Development of the Tide－Generating Potential［J］.Proceedings of the Royal Society of London，1921，100（704）：305－329.

［10］ DOODSON A T.The Analysis of Tidal Observations［J］.Philosophical Transactions of the Royal Society of London，1928，227：223－279.

［11］ XI Qinwen，HOU Tianhang.A New Complete Development of the Tide－generating Potential for the Epoch J2000.0［J］.Acta Geophysica Sinica，1987，30（4）：349－362.（郗欽文，侯天航.新的引潮位完全展開［J］.地球物理學(xué)報(bào)，1987，30（4）：349－362.）

［12］ XI Qinwen.The Evaluation of the Precision of the Development of the Tidal Generation Potential［J］.Acta Geophysica Sinica，1992，35（2）：150－153.（郗欽文.精密引潮位展開的精度評(píng)定［J］.地球物理學(xué)報(bào)，1992，35（2）：150－153.）

［13］ HUANG Zuke，CHEN Zongyong，SI Hongye，et al.Analysis of 19Year Tidal Data［J］.Science in China（Series D），1997，27（2）：174－179.（黃祖柯，陳宗鏞，司鴻業(yè)，等.我國沿海若干驗(yàn)潮站的19年潮汐分析［J］.中國科學(xué)：D輯，1997，27（2）：174－179.）

［14］ CHEN Junyong.Thinking on the Establishment of Modern Geodetic Datum of China［J］.Geomatics and Information Science of Wuhan University，2003，28（S）：1－6.（陳俊勇.建設(shè)我國現(xiàn)代大地測(cè)量基準(zhǔn)的思考［J］.武漢大學(xué)學(xué)報(bào)：信息科學(xué)版，2003，28（S）：1－6.）

［15］ CHEN Junyong.On the Chinese Modern Geodetic Coordinate System and Height System ［J］.Bulletin of Surveying and Mapping，2002（8）：1－5.（陳俊勇.對(duì)我國建立現(xiàn)代大地坐標(biāo)系統(tǒng)和高程系統(tǒng)的建議［J］.測(cè)繪通報(bào)，2002（8）：1－5.）

［16］ WEN Yuanlan，YANG Yuanxi.Research on Sea Level and Rising Trend of Coastal Waters in China［J］.Geomatics and Information Science of Wuhan University，2001，26（2）：127－131.（文援蘭，楊元喜.我國近海平均海面及其變化的研究［J］.武漢大學(xué)學(xué)報(bào)：信息科學(xué)版，2001，26（2）：127－131.

［17］ QIAO Jianrong.Mean Sea Level Changes on the Yellow Sea for 60Years［J］.Coastal Engineering，1985，4（2）：83－88.（喬建榮.六十年來黃海平均海平面的變化［J］.海岸工程，1985，4（2）：83－88.）

［18］ YU Ling，F(xiàn)ANG Guohong.Tidal Harmonic Analysis for Long Term Sea Level Data of Tidal Gauge Stations in Seas Adjacent to China［J］.Advances in Marine Science，2010，28（2）：133－141.（于玲，方國洪.中國近海驗(yàn)潮站長(zhǎng)期觀測(cè)資料的潮汐調(diào)和分析［J］.海洋科學(xué)進(jìn)展，2010，28（2）：133－141.）

［19］ ZHENG Wenzhen.Distribution of Annual Rates of Sea Level and Variation of Long－period Constituents in China［J］.Marine Science Bulletin，1999，18（4）：1－10.（鄭文振.我國海平面年速率的分布和長(zhǎng)周期分潮的變化［J］.海洋通報(bào)，1999，18（4）：1－10.）

［20］ VAN ETTE A C M，SCHOEMAKER H J.Harmonic Analyses of Tides Essential Features and Disturbing Influences［M］.［S.l.］：Netherland Hydrographer，1966：79－107.

［21］ TIAN Hui，CHEN Zongyong.Stochastic Dynamical Analyses for Monthly Mean Sea Level along the Coast of China［J］.Acta Oceanologica Sinica，1988，28（4）：9－15.（田暉，陳宗鏞.中國沿岸近期多年月平均海面隨機(jī)動(dòng)態(tài)分析［J］.海洋學(xué)報(bào)，1988，28（4）：9－15.）