李 榮 暴景陽(yáng) 許 軍 劉 聚

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

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

李 榮 暴景陽(yáng) 許 軍 劉 聚

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

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

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

1 引言

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

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

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

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

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

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

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

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

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

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

利用上式就可以作一個(gè)概率表達(dá)式：

表1 t檢驗(yàn)臨界值表

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

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

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

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

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

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

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

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

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

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

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

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

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

4 結(jié)語(yǔ)

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

［1］黃祖珂，黃磊.潮汐原理與計(jì)算［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áng)，劉雁春，肖付民.潮汐分析與預(yù)報(bào)的誤差分析［J］. 海洋測(cè)繪，1995，15（1）：31-37.

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

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

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

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

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

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

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

［12］暴景陽(yáng)，劉雁春.海道測(cè)量水位控制方法研究［J］.測(cè)繪科學(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日

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