中國沿岸驗潮站潮汐調(diào)和常數(shù)的置信度確定?

李 榮 暴景陽 許 軍 劉 聚

（海軍大連艦艇學(xué)院海洋測繪系 大連 116018）

中國沿岸驗潮站潮汐調(diào)和常數(shù)的置信度確定?

李 榮 暴景陽 許 軍 劉 聚

（海軍大連艦艇學(xué)院海洋測繪系 大連 116018）

對中國沿岸七個具有代表性的長期驗潮站按照年觀測序列進行了調(diào)和分析結(jié)果的精度統(tǒng)計，并根據(jù)間接平差的條件對分潮振幅進行了t分布檢驗，還對具有缺測數(shù)據(jù)的年觀測資料得來的主要分潮調(diào)和常數(shù)進行了置信度確定。結(jié)果表明，由完整年觀測資料得來的分潮振幅具有毫米級精度。根據(jù)α=0.1進行t檢驗后，由年觀測序列得來的122個分潮有近一半的分潮對實際計算水位高度沒有較大作用。由具有缺測數(shù)據(jù)的年觀測資料得來的長周期分潮振幅的置信度下降明顯。

驗潮站；潮汐；調(diào)和常數(shù)；假設(shè)檢驗

1 引言

潮汐作為海洋里較典型的現(xiàn)象，掌握其規(guī)律性變化對海洋工程建設(shè)、航道交通、海洋防災(zāi)減災(zāi)等具有重要的作用［11］。對潮汐調(diào)和分析的結(jié)果進行精度評定或檢驗對于準確計算調(diào)和常數(shù)具有一定的必要性，能為調(diào)和常數(shù)的有效應(yīng)用提供一定的質(zhì)量參考依據(jù)。本文將利用多個驗潮站長期觀測數(shù)據(jù)的調(diào)和分析結(jié)果對計算出來的調(diào)和常數(shù)進行假設(shè)檢驗，驗證其計算結(jié)果的可信度，為后續(xù)的深度基準計算、潮汐預(yù)報與分析等提供參考。

2 長期數(shù)據(jù)下的調(diào)和分析

2.1 調(diào)和分析原理

調(diào)和分析的類型根據(jù)觀測數(shù)據(jù)的時間長度可分為長期、中期、短期觀測分析［1～2］，進行年觀測分析時基于某時的實測水位高度表達式：

利用某時的實測水位高度表達式，將測得的n個數(shù)據(jù)代入表達式中，列出n個方程求解2m+1個未知數(shù)（n＞2m+1）。 ζ(t)為t時刻的水位；a0為觀測時段的平均海面；i代表某一個分潮，即分潮序號；f、u為交點因子與交點訂正角，是因月球軌道變化對振幅與相角的訂正；σ是某分潮的角速率，通常情況下是已知的；(V0+u)為天文相角；r(t)為水位觀測時的噪聲，也可以理解為觀測時的誤差。剩下的兩項H、g即為需要求解的潮汐調(diào)和常數(shù)。

2.2 調(diào)和常數(shù)的檢驗指標

利用方程組對調(diào)和常數(shù)進行求解時，總是會有誤差的存在，從數(shù)據(jù)的獲取到數(shù)據(jù)的計算處理，并且根據(jù)誤差傳播律進行傳播［3～6］，對誤差的控制與評估可以為計算結(jié)果的有效利用提供一定的參考依據(jù)，從而使運用更加合適準確。

利用式（1）求解出來的參數(shù)為振幅的正余弦值［7］，因此所得出的中誤差為對應(yīng)正余弦值的中誤差，在利用t檢驗對振幅H進行檢驗時所需要用到的協(xié)因數(shù)陣則需要根據(jù)誤差傳播律對其進行計算，以得到相應(yīng)振幅的協(xié)因數(shù)。可根據(jù)下式進行協(xié)因數(shù)的計算［5］：

利用長期觀測數(shù)據(jù)進行調(diào)和常數(shù)解算時，可以得出一百多個分潮的振幅和遲角，但它們并不是都具有比較重要的貢獻［12］，有些數(shù)量級較小的分潮振幅可以認為近似等于零，即在水位高度的計算時意義不大。在運用最小二乘法計算調(diào)和常數(shù)時實際上就是一個間接平差的過程［8～9］，因此可以采用t檢驗法，作t統(tǒng)計量：

作 原 假 設(shè)為 H0：E()=0，備選 假 設(shè)H1：E()≠0 ，式中 t(f)為服從t分布的統(tǒng)計量，f就是多余觀測數(shù)，也稱自由度；為單位權(quán)中誤差，為平差參數(shù)的協(xié)因數(shù)陣。

利用上式就可以作一個概率表達式：

表1 t檢驗臨界值表

3 調(diào)和常數(shù)的假設(shè)檢驗

3.1 有效分潮數(shù)的確定

本次檢驗所采用的數(shù)據(jù)來源于夏威夷大學(xué)海平面中心，有東方、海口、連云港、呂泗、汕尾、坎門、北海七個長期驗潮站的逐時水位觀測資料，部分驗潮站可能存在少量缺測數(shù)據(jù)，但不影響其總體分析。一年觀測數(shù)據(jù)分析可得出一百多個分潮的調(diào)和常數(shù)，由于篇幅限制，在此僅列出部分驗潮站以及部分主要分潮的調(diào)和分析結(jié)果。

表2 各驗潮站年調(diào)和分析結(jié)果統(tǒng)計表 單位：cm

從上表可以看出，由年觀測資料獲得的四個主要分潮振幅中誤差普遍在毫米級，維持在±0.2cm～±0.4cm之間。個別年份的振幅中誤差達到了厘米級，分別出現(xiàn)在北海站的76、77年，連云港站的76年，呂泗站的77年，四個主要分潮振幅中誤差介于±0.9cm～±3.3cm之間。細究其原因發(fā)現(xiàn)，四份振幅中誤差偏大的年觀測資料里均出現(xiàn)了缺測數(shù)據(jù)，誤差最大的出現(xiàn)在北海站的76年，該年實測數(shù)據(jù)8491個，缺測293個；北海站77年實測數(shù)據(jù)8712個，缺測48個；連云港站76年實測數(shù)據(jù)8748個，缺測36個；呂泗站77年實測數(shù)據(jù)8734個，缺測26個。綜上可得，振幅中誤差大小與缺測數(shù)據(jù)個數(shù)呈現(xiàn)出一定的正相關(guān)。

經(jīng)過以上的數(shù)據(jù)準備，將計算出的調(diào)和常數(shù)用式（4）計算t值，再根據(jù)查表可得，當自由度大于200時，t分布已趨于正態(tài)分布，并且在α等于0.1，即置信度為90%時t檢驗值為1.645，而由年觀測資料進行調(diào)和分析時自由度遠大于200，因此統(tǒng)計t值大于1.645的分潮個數(shù)，即拒絕H0：E()=0，認為H1：E()≠0成立，這些即是對水位高度具有較大貢獻的分潮，統(tǒng)計值列于表3。

表3 調(diào)和常數(shù)假設(shè)檢驗結(jié)果統(tǒng)計表

上表結(jié)果統(tǒng)計于七個長期驗潮站各15年觀測數(shù)據(jù)的結(jié)果，在不出現(xiàn)缺測數(shù)據(jù)時，平年逐時觀測數(shù)據(jù)達8760個，閏年逐時觀測數(shù)據(jù)達8784個。通過上表可以發(fā)現(xiàn)，在不出現(xiàn)缺測數(shù)據(jù)時，經(jīng)過置信度為90%的t檢驗后，各驗潮站由年觀測資料獲得的分潮振幅不等于零的分潮個數(shù)介于54～81之間，也就是說由年觀測資料獲得的122個分潮在實際計算水位高度時我們可以認為有近一半的分潮可以不參與計算，只需要利用剩下的有效分潮的調(diào)和常數(shù)，即可對水位高度進行較為準確的計算。呂泗站在1988年出現(xiàn)大量的缺測數(shù)據(jù)，有效分潮數(shù)僅為1，當缺測數(shù)據(jù)較多時不再符合年調(diào)和分析的條件。

通過上表也可以發(fā)現(xiàn)，當年觀測資料出現(xiàn)缺測數(shù)據(jù)時，大于α=0.1的t檢驗值的分潮個數(shù)會出現(xiàn)較大程度的減少。如連云港站在不出現(xiàn)缺測數(shù)據(jù)的年份有效分潮數(shù)為70個左右，當出現(xiàn)缺測數(shù)據(jù)時則減少到29、23、9個，即有效分潮個數(shù)隨著缺測數(shù)據(jù)的增多而減少，這也與表2中振幅中誤差偏大的年份相對應(yīng)。

3.2 主要分潮調(diào)和常數(shù)的置信度確定

對主要分潮振幅進行t檢驗時，自由度大于200，查表得當 α=0.01，t值為2.576；當年觀測資料沒有缺測數(shù)據(jù)時，7個驗潮站的13個主要分潮的t值遠大于2.576，即由完整年觀測資料得來的13個主要分潮的調(diào)和常數(shù)置信度在99%以上。這里主要對出現(xiàn)缺測數(shù)據(jù)的年觀測資料進行主要分潮調(diào)和常數(shù)的置信度評估。出現(xiàn)缺測數(shù)據(jù)的有連云港站的76、81年，呂泗站的77年，北海站的76、77年，?？谡镜?9年。對這6份年觀測資料得來的主要分潮振幅進行置信度確定，統(tǒng)計結(jié)果列于表4。

在年觀測資料出現(xiàn)缺測數(shù)據(jù)時，Ssa、M4、MS4、M6分潮振幅置信度明顯下降，或者說其振幅結(jié)果可能存在較大誤差；Q1、K2分潮振幅置信度從99%分別下降到95%、98%，雖然下降不明顯，但其振幅結(jié)果受到一定影響。其它的主要分潮振幅在出現(xiàn)缺測數(shù)據(jù)的情況下置信度沒有受到明顯影響［10］，這類分潮周期均為一日或半日；而置信度明顯下降的分潮均為長周期分潮，這也驗證了分潮周期與數(shù)據(jù)長度之間的關(guān)系。

表4 主要分潮振幅置信度結(jié)果統(tǒng)計表

4 結(jié)語

由年觀測資料得來的分潮振幅均可達到毫米級精度，若年觀測資料出現(xiàn)缺測數(shù)據(jù)時，分潮振幅中誤差會明顯增大，且中誤差大小與缺測數(shù)據(jù)個數(shù)呈現(xiàn)出正相關(guān)。由年觀測資料可計算得到122個分潮調(diào)和常數(shù)，在實際計算水位高度時具有較大作用的分潮個數(shù)大約為前者的一半，而當年觀測資料出現(xiàn)缺測數(shù)據(jù)，有效分潮數(shù)明顯減少，反映出年觀測資料不完整時對分潮調(diào)和常數(shù)結(jié)果具有較大影響。在有條件的情況下，建議年調(diào)和分析時選取沒有缺測數(shù)據(jù)的年觀測資料進行。

［1］黃祖珂，黃磊.潮汐原理與計算［M］.青島：青島海洋大學(xué)出版社，2005：75-89.

［2］PATRICK F C I，JOSEF Y C I，MICHAEL G G，et al.North Pacific internal tides from the Aleutian Ridge：altim?eter observations and modeling［J］.Journal of Marine Re?search，2001，59（2）：167-191.

［3］A S Franco，J Harari.On the Stability of Long Series Tidal Analyses［J］.International Hydrographic Review，1993（1）：77-89.

［4］暴景陽，劉雁春，肖付民.潮汐分析與預(yù)報的誤差分析［J］. 海洋測繪，1995，15（1）：31-37.

［5］暴景陽，許軍.中國沿岸驗潮站潮汐調(diào)和常數(shù)的精度評估［J］. 海洋測繪，2013，33（1）：1-4.

［6］武漢大學(xué)測繪學(xué)院測量平差學(xué)科組.誤差理論與測量平差基礎(chǔ)［M］.武漢：武漢大學(xué)出版社，2009：216-221.

［7］方國洪，鄭文振，陳宗鏞.等.潮汐和潮流的分析和預(yù)報［M］.北京：海洋出版社，1986：67-75.

［8］于玲，方國洪.中國近海驗潮站長期觀測資料的潮汐調(diào)和分析［J］. 海洋科學(xué)進展，2010，02：133-141.

［9］暴景陽，晁定波，李建成.南中國海TOPEX_POSEID軌跡交叉點測高數(shù)據(jù)的潮汐調(diào)和分析［J］.測繪學(xué)報，2000，29（1）：17-23.

［10］王驥，方國洪.不完整逐時潮汐觀測記錄的分析［J］.海洋學(xué)報，1981，03（2）：193-210.

［11］劉雁春.海洋測繪空間結(jié)構(gòu)及數(shù)據(jù)處理［M］.北京：測繪出版社，2003：55-68.

［12］暴景陽，劉雁春.海道測量水位控制方法研究［J］.測繪科學(xué)，2006，06：49-51.

Hypothesis Test of Harmonic Constants for Long Term Tidal Stations Along the Coast of China

LI Rong BAO Jingyang XU Jun LIU Ju
（Department of Hydrography and cartography，Dalian Naval Academy，Dalian 116018）

For seven typical long term tidal stations along the coast of China，the accuracy evaluations are performed based on statistics of the analyzed results series from yearly observations.The T-distribution of the tidal amplitudes is tested according to the conditions of indirect adjustment method.The confidence level of the harmonic constants from yearly observations with the missing data is determined.The statistics shows that the harmonic constants from the complete year observation data can reach mm accuracy level.According to the T-test，nearly half of the harmonic constants from the yearly observation do not have a significant effect on the actual calculated water level.The confidence level of the amplitude of long period constituents from yearly observations with the missing data is significantly reduced.

tidal station，tide，tidal harmonic constants，hypothesis test

P731

10.3969/j.issn.1672-9730.2017.11.032

Class Number P731

2017年5月20日，

2017年6月29日

李榮，男，碩士研究生，研究方向：海洋大地測量理論與方法。暴景陽，男，博士，教授，博士生導(dǎo)師，研究方向：海洋大地測量、海道測量、測量數(shù)據(jù)處理等。許軍，男，博士，講師，研究方向：水位控制、潮汐模型與海域垂直基準模型構(gòu)建的理論及其應(yīng)用。劉聚，男，博士研究生，研究方向：海洋測繪水位控制技術(shù)。