卓玉生，王 巖

（中交第一航務(wù)工程勘察設(shè)計(jì)院有限公司，天津 300220）

引言

我國(guó)對(duì)理論最低潮面的計(jì)算方法的規(guī)定有一個(gè)發(fā)展過(guò)程，而通常情況下，理論最低潮面一旦確定，不宜輕易變動(dòng)，這就導(dǎo)致國(guó)內(nèi)理論最低潮面采用值復(fù)雜多樣。同時(shí)，相同計(jì)算方法下，計(jì)算時(shí)采用不同的數(shù)據(jù)長(zhǎng)度、年份，也可能得到不同的結(jié)果，給港口與航道工程及航運(yùn)帶來(lái)較多困擾。在此，對(duì)確定理論最低潮面的計(jì)算方法進(jìn)行總結(jié)，并對(duì)其影響因素進(jìn)行匯總分析，為確定理論最低潮面給出方法和建議。

1 理論最低潮面及相關(guān)概念

理論最低潮面，是根據(jù)前蘇聯(lián)弗拉基米爾斯基提出的方法[1]計(jì)算出的平均海面以下理論上能達(dá)到的最低潮位，通常記為平均海面到理論最低潮面的距離。理論最低潮面是國(guó)際上常用的深度基準(zhǔn)面之一，也是常用的潮高基準(zhǔn)面之一。深度基準(zhǔn)面是水深測(cè)量工作中水深的起算面，現(xiàn)行國(guó)家標(biāo)準(zhǔn)《海道測(cè)量規(guī)范》（GB12327-1998）明確規(guī)定“以理論最低潮面作為深度基準(zhǔn)面”[2]1，其算法為弗拉基米爾斯基方法的改進(jìn)算法?，F(xiàn)行國(guó)家標(biāo)準(zhǔn)《中國(guó)海圖圖式》（GB12319-1998）也做出規(guī)定：“中國(guó)沿海地區(qū)一般采用‘理論最低潮面’作為深度基準(zhǔn)”[3]。根據(jù)國(guó)家海洋信息中心編制的潮汐表[4]說(shuō)明中描述，“潮高基準(zhǔn)面與海圖深度基準(zhǔn)面一致”。

港口與航道工程及航運(yùn)中，確定適宜的理論最低潮面數(shù)值是一項(xiàng)重要的基礎(chǔ)性工作。確定深度基準(zhǔn)面時(shí)，既要考慮船舶的航行安全，又要考慮航道的利用率，這是確定海圖深度基準(zhǔn)面時(shí)的基本原則。深度基準(zhǔn)面定得過(guò)高，船舶容易發(fā)生事故；深度基準(zhǔn)面定得過(guò)低，則降低了航道的利用率。

2 理論最低潮面計(jì)算原理

理論最低潮面概念是前蘇聯(lián)弗拉基米爾斯基提出的，其計(jì)算原理概括為：利用潮位調(diào)和常數(shù)，選取振幅較大的13 個(gè)分潮，通過(guò)一些規(guī)定和算法，求取平均海面下能達(dá)到的最低潮位。其算法中的13 個(gè)分潮包括三部分：M2、S2、N2、K2、K1、O1、P1、Q1等8 個(gè)主要分潮，M4、MS4、M6等3 個(gè)淺水分潮和Sa、Ssa等2 個(gè)長(zhǎng)周期分潮。為方便進(jìn)行算法對(duì)比分析，以下內(nèi)容對(duì)算法中間過(guò)程進(jìn)行了提煉和概化，詳細(xì)過(guò)程可參見(jiàn)文獻(xiàn)[1]或文獻(xiàn)[5]。

1）8 個(gè)主要分潮部分

8個(gè)主要分潮在平均海面上t時(shí)的潮高h(yuǎn)8(t)可表示為：

式中f、H、g、σ分別為各分潮的交點(diǎn)因子、振幅、遲角、角速度，v0+u為各分潮在格林威治1 月1 日0 時(shí)的相角，t為平太陽(yáng)時(shí)。

通過(guò)形式變換(1)式可變換為以下形式：

式中φ為各分潮的相角，均可用以φK1為變量的函數(shù)表示，R、ε為輔助變量，其中R始終為正值。當(dāng)φM2-ε1=180°、φS2-ε2=180°φN2-ε3=180°時(shí)，cos180°=-1，h8(φK1)可能為極小值，用L8(φK1)表示，即

φK1從0°到360°變化，選取L8(φK1)的最小值，即為8 個(gè)主要分潮組合出的最低潮位，記為L(zhǎng)8，對(duì)應(yīng)的φK1相角記為(φK1)8min。

2）3 個(gè)淺水分潮改正部分

3 個(gè)淺水分潮在平均海面上t時(shí)的潮高改正可表示為以φK1為變量的函數(shù)：

式中φ為各分潮的相角，均可用以φK1為變量的函數(shù)表示。

根據(jù)L8對(duì)應(yīng)的相角(φK1)8min,求出3 個(gè)淺水分潮的相角φ，計(jì)算各分潮的fHcosφ值，最后求其總和，記為ΔL3。

3）2 個(gè)長(zhǎng)周期分潮改正部分

對(duì)于平均海面季節(jié)變化較大的海區(qū)，可根據(jù)由一年以上潮位資料分析所得潮汐調(diào)和常數(shù)的長(zhǎng)周期分潮Sa和Ssa來(lái)進(jìn)行改正。

2 個(gè)長(zhǎng)周期分潮在平均海面上t 時(shí)的潮高改正可表示為以φK1為變量的函數(shù)，記為

式中φ為各分潮的相角，均可用以φK1為變量的函數(shù)表示。

根據(jù)L8對(duì)應(yīng)的相角(φK1)8min在其±60°范圍內(nèi)以15°間隔為樣本，分別通過(guò)查表（或通過(guò)公式）求取ε2等相關(guān)輔助角[1]，進(jìn)而分別求出2 個(gè)長(zhǎng)周期分潮的相角，計(jì)算2 個(gè)長(zhǎng)周期分潮的貢獻(xiàn)之和，選取最小值，記為ΔL2。

該項(xiàng)改正也可采用海平面季節(jié)改正數(shù)。

以上三部分，L8部分算法較穩(wěn)定；ΔL3部分在3 個(gè)淺水分潮振幅之和在20 cm 附近時(shí)，由于其算法的不連續(xù)性，淺水分潮調(diào)和常數(shù)的較小變化，可能出現(xiàn)理論最低潮面的較大變化；ΔL2部分對(duì)“平均海面季節(jié)變化較大”沒(méi)有量化的判別，對(duì)其理解容易產(chǎn)生分歧，從而導(dǎo)致不同人使用相同潮位數(shù)據(jù)計(jì)算理論最低潮面結(jié)果不唯一。

3 確定理論最低潮面的影響因素

3.1 理論最低潮面算法演變

通常深度基準(zhǔn)面一旦被正式采用，不宜輕易改變?，F(xiàn)行國(guó)家標(biāo)準(zhǔn)《海道測(cè)量規(guī)范》（GB 12327-1998）明確規(guī)定[2]：“深度基準(zhǔn)面一經(jīng)確定且在正規(guī)水深測(cè)量中已被采用者，一般不得變動(dòng)”。該規(guī)定弱化了現(xiàn)行規(guī)范算法的強(qiáng)制執(zhí)行要求，基于此點(diǎn)，很多港口或海域采用的理論最低潮面數(shù)值并未因規(guī)范規(guī)定算法的變化而改變，從而加劇了我國(guó)沿海采用的深度基準(zhǔn)面的復(fù)雜多樣。在此，列出我國(guó)理論最低潮面算法的主要演變過(guò)程，并對(duì)比各算法的特點(diǎn)和相互之間的差異。由于對(duì)理論最低潮面的理解有分歧，諸多文獻(xiàn)發(fā)布了不同的改進(jìn)算法[6-8]，在此不作描述。

1956 年以前，我國(guó)采用的深度基準(zhǔn)面不統(tǒng)一，主要采用略最低低潮面作為深度基準(zhǔn)面，也有采用實(shí)測(cè)最低潮位、實(shí)測(cè)中等水位、平均海面等其他潮面的。

1956 年以后，海軍司令部海道測(cè)量部在全國(guó)海洋測(cè)繪中，統(tǒng)一采用理論最低潮面（初期稱“理論深度基準(zhǔn)面”），其算法主要依據(jù)前蘇聯(lián)弗拉基米爾斯基方法（以下簡(jiǎn)稱“弗拉基方法”）。1959 年海軍司令部海道測(cè)量部出版的《實(shí)用潮汐學(xué)》是早期重要的參考文獻(xiàn)。

1975 年，《海道測(cè)量規(guī)范》-水深測(cè)量（以下簡(jiǎn)稱“75 規(guī)范”）規(guī)定:“長(zhǎng)、短期驗(yàn)潮站的深度基準(zhǔn)面，根據(jù)八個(gè)分潮（M2、S2、N2、K2、K1、O1、P1、Q1）的潮汐調(diào)和常數(shù)H、g值，按理論深度基準(zhǔn)面的公式計(jì)算。當(dāng)M4、MS4、M6三個(gè)淺海分潮振幅之和大于20 cm 時(shí)，還要進(jìn)行淺水分潮改正”[9]。75 規(guī)范中考慮的分潮共計(jì)11 個(gè)，規(guī)范中未給出理論最低潮面具體算法。不考慮長(zhǎng)周期分潮改正，保證了采用相同潮位數(shù)據(jù)計(jì)算結(jié)果的唯一性，但該算法降低了深度基準(zhǔn)面的保證率。

1990 年，《海道測(cè)量規(guī)范》（GB 12327-90）（以下簡(jiǎn)稱“90 規(guī)范”）規(guī)定:“以理論最低潮面作為深度基準(zhǔn)面”[10]。同時(shí)，90 規(guī)范給出了理論最低潮面計(jì)算公式[10]，除計(jì)算公式（明顯不合理，應(yīng)屬排版疏漏）外，與《海道測(cè)量規(guī)范》（GB 12327-98）基本一致，在此不做詳述。

1998 年，《海道測(cè)量規(guī)范》（GB 12327-1998）（以下簡(jiǎn)稱“98 規(guī)范”）修正了90 規(guī)范理論最低潮面算法中的明顯疏漏，是現(xiàn)行國(guó)家強(qiáng)制標(biāo)準(zhǔn)。98規(guī)范明確規(guī)定了用13 分潮計(jì)算理論最低潮面的計(jì)算公式[2]，可概化如下：

98 規(guī)范算法是弗拉基方法的改進(jìn)算法，改進(jìn)后的計(jì)算結(jié)果有時(shí)會(huì)與原算法有較大差異，在做國(guó)外項(xiàng)目時(shí)應(yīng)注意與當(dāng)?shù)氐睦碚撟畹统泵嫠惴ǖ牟顒e。

對(duì)比不同時(shí)期深度基準(zhǔn)面算法之間的差異，將不同時(shí)期算法列于下表1：

表1 不同深度基準(zhǔn)面算法對(duì)比

以下實(shí)例中，用2018～2019 年兩年潮汐表[4]所載逐時(shí)潮位預(yù)報(bào)值進(jìn)行調(diào)和分析，取得黃、渤海29 個(gè)站位的調(diào)和常數(shù)，分別按以上三種算法計(jì)算理論最低潮面，繪制與潮汐表采用值對(duì)比圖（見(jiàn)圖1）。其中的弗拉基方法，對(duì)比計(jì)算過(guò)程中均認(rèn)為“平均海面季節(jié)變化較大”，“平均海面季節(jié)變化較小”情況與75 規(guī)范算法一致，不再單獨(dú)列出。

圖1 黃、渤海各站潮汐表與三種理論最低潮面算法結(jié)果對(duì)比

為方便顯示黃、渤海各站三種規(guī)范算法計(jì)算結(jié)果與潮汐表所載理論最低潮面差異，繪制差值圖見(jiàn)圖2。

圖2 黃、渤海各站三種理論最低潮面算法結(jié)果與潮汐表采用值差值

從以上計(jì)算及對(duì)比結(jié)果可以看出，大部分潮汐表采用值與98 規(guī)范規(guī)定算法計(jì)算結(jié)果存在明顯差異，應(yīng)非計(jì)算誤差所致，而是采用了不同算法。

3.2 采用不同數(shù)據(jù)資料引起的差異

1）采用不同長(zhǎng)度的潮位數(shù)據(jù)

采用不同長(zhǎng)度的潮位數(shù)據(jù)，調(diào)和分析結(jié)果會(huì)有差異。一般情況下，采用的潮位數(shù)據(jù)資料越長(zhǎng)，調(diào)和分析可靠性越高，理論最低潮面計(jì)算結(jié)果穩(wěn)定性越好。

由于月球升交點(diǎn)作18.61 年周期變化，較多文獻(xiàn)提出引入潮汐歷元和深度基準(zhǔn)面歷元，長(zhǎng)度為19 年。通常情況下，采用歷元周期長(zhǎng)度的潮位資料，既可以提高計(jì)算理論最低潮面的可靠性，也能修正非周期性變化影響。國(guó)際上，美國(guó)、澳大利亞等國(guó)均有潮汐歷元和深度基準(zhǔn)面歷元長(zhǎng)度的規(guī)定，如美國(guó)曾將潮汐和深度基準(zhǔn)歷元確定為1960～1978 年，2003 年起重定為1983～2001 年,長(zhǎng)度為19 年；而澳大利亞潮汐和平均海面委員會(huì)建議采用的潮汐基準(zhǔn)面歷元為1992～2011 年,長(zhǎng)度為20 年[11]。我國(guó)尚未有歷元規(guī)定，為提高計(jì)算理論最低潮面的可靠性，可盡量增加采用的潮位資料長(zhǎng)度。鑒于我國(guó)海平面上升引起的潮波改變、岸線及水下地形變化、局部地區(qū)潮差非周期性變化等因素，有條件時(shí)，宜分析理論最低潮面的非周期性變化趨勢(shì)，確定適宜的理論最低潮面數(shù)值。

在此呼吁掌握長(zhǎng)期潮位資料的權(quán)威部門，確定潮汐和深度基準(zhǔn)面歷元時(shí)段，并給出不同海域主港的歷元時(shí)段平均海面、理論最低潮面和深度基準(zhǔn)面保證率等數(shù)值，將會(huì)對(duì)規(guī)范國(guó)內(nèi)深度基準(zhǔn)面體系起到關(guān)鍵性作用。不宜變動(dòng)深度基準(zhǔn)面時(shí)，也會(huì)對(duì)評(píng)價(jià)采用的深度基準(zhǔn)面提供依據(jù)。

2）采用不同時(shí)段的潮位數(shù)據(jù)

搜集夏威夷大學(xué)海平面中心（UHSLC）發(fā)布的連云港等8 站1975 年～1997 年連續(xù)23 年逐時(shí)實(shí)測(cè)潮位資料，分別以每年潮位資料為樣本進(jìn)行調(diào)和分析，然后分別以98 規(guī)范算法計(jì)算理論最低潮面，繪制各站理論最低潮面逐年計(jì)算結(jié)果見(jiàn)圖3。

圖3 各站1975 年～1997 年逐年潮位計(jì)算理論最低潮面結(jié)果

根據(jù)計(jì)算結(jié)果，統(tǒng)計(jì)各站23 年逐年計(jì)算理論最低潮面時(shí)，不同年份之間最大值、最小值和最大互差，并給出用23 年資料計(jì)算出的理論最低潮面，結(jié)果見(jiàn)表2。

表2 各站23 年逐年計(jì)算理論最低潮面統(tǒng)計(jì)結(jié)果

通過(guò)對(duì)比分析可以看出，采用不同年份資料推算理論最低潮面可能出現(xiàn)明顯差異，本例中按現(xiàn)行98 規(guī)范算法最大互差為25cm。

3.3 潮位觀測(cè)位置變化引起的差異

長(zhǎng)期潮位觀測(cè)過(guò)程中，驗(yàn)潮站位置發(fā)生變化的情況時(shí)有發(fā)生，也有驗(yàn)潮站所在位置因工程建設(shè)導(dǎo)致潮汐特征改變情況出現(xiàn)。在變動(dòng)前、后位置潮汐特征相似的情況下，為延長(zhǎng)潮位資料序列，應(yīng)進(jìn)行變動(dòng)前、后位置同步驗(yàn)潮，對(duì)原位置潮位資料進(jìn)行必要訂正。變動(dòng)前、后位置潮汐性質(zhì)變化較大情況下，應(yīng)作為兩個(gè)測(cè)站使用。

隨著我國(guó)港口建設(shè)的發(fā)展，一些長(zhǎng)期驗(yàn)潮站所在位置的潮汐性質(zhì)因地形地貌變化（如在原本開(kāi)敞環(huán)境下變?yōu)楠M長(zhǎng)通道內(nèi)部情況）也會(huì)發(fā)生變化，進(jìn)行潮汐分析時(shí)也要考慮。

3.4 調(diào)和分析的算法影響

調(diào)和分析的算法間接影響計(jì)算理論最低潮面結(jié)果。調(diào)和分析的算法中，能夠分離出的調(diào)和分潮個(gè)數(shù)和方程組求解方法不同，都會(huì)造成計(jì)算結(jié)果差異。由于現(xiàn)在的潮汐調(diào)和分析理論已十分成熟，這類差異通常較小。

4 理論最低潮面計(jì)算結(jié)果檢驗(yàn)

4.1 深度基準(zhǔn)面保證率驗(yàn)證

深度基準(zhǔn)面保證率是常用的評(píng)價(jià)海圖深度基準(zhǔn)面是否適宜的標(biāo)準(zhǔn)之一。

深度基準(zhǔn)面保證率的計(jì)算公式為：

深度基準(zhǔn)面保證率=高于深度基準(zhǔn)面的低潮次數(shù)/低潮總次數(shù)×100 %。

通常認(rèn)為深度基準(zhǔn)面保證率的適宜數(shù)值為 95 %。也有深度基準(zhǔn)面的保證率應(yīng)大于98 %，特殊情況下不得小于95 %的說(shuō)法[12]，一般要求統(tǒng)計(jì)深度基準(zhǔn)面保證率應(yīng)有一年以上的水位觀測(cè)資料。實(shí)際上，統(tǒng)計(jì)采用的資料越長(zhǎng)，保證率的可靠性越高，遇到異常年份，僅采用一年水位資料統(tǒng)計(jì)深度基準(zhǔn)面保證率可能產(chǎn)生較大偏差,這種情況在保證率較低時(shí)較明顯。

4.2 附近港口或海域長(zhǎng)期站采用值對(duì)照

計(jì)算理論最低潮面時(shí)，如附近港口或海域有長(zhǎng)期潮位站，可與長(zhǎng)期站采用值進(jìn)行對(duì)照，必要時(shí)進(jìn)行同步觀測(cè)，與長(zhǎng)期站潮位資料做相關(guān)分析，修正觀測(cè)時(shí)段潮位代表性影響，以提高計(jì)算理論最低潮面的可靠性。

5 結(jié)語(yǔ)

通過(guò)以上分析，在當(dāng)前形勢(shì)下，建議按以下思路確定理論最低潮面：

1）確定理論最低潮面算法；

首先確定理論最低潮面的算法，以保證與附近港口或海域理論最低潮面算法協(xié)調(diào)一致。

2）計(jì)算理論最低潮面；

計(jì)算時(shí)應(yīng)采用可靠的方法或軟件。

3）采用的潮位資料的代表性分析和修正；

理論最低潮面計(jì)算結(jié)果應(yīng)與附近港口或長(zhǎng)期驗(yàn)潮站采用值進(jìn)行對(duì)比分析，有條件應(yīng)進(jìn)行短期同步觀測(cè)，通過(guò)相關(guān)分析等方法，確定采用的潮位資料的代表性，必要時(shí)應(yīng)進(jìn)行修正。

4）安全性評(píng)價(jià)。

對(duì)計(jì)算的理論最低潮面結(jié)果應(yīng)進(jìn)行安全性評(píng)價(jià)。安全性和與附近港口或海域采用值的協(xié)調(diào)性沖突時(shí)，應(yīng)以安全性為重，重新確定算法。

總之，理論最低潮面的確定不應(yīng)該是簡(jiǎn)單的計(jì)算，而是要在保證港口與航道安全的前提下，兼顧利用率，同時(shí)盡可能與附近港口或海域的采用值保持協(xié)調(diào)一致。