李希彬，張秋豐，姚志剛，葉風(fēng)娟，李玉杰

（1.國(guó)家海洋局天津海洋環(huán)境監(jiān)測(cè)中心站，天津300451；2.天津市海洋環(huán)境監(jiān)測(cè)預(yù)報(bào)中心，天津300451；3.中國(guó)海洋大學(xué)海洋環(huán)境學(xué)院，山東青島 266100）

天津近海潮汐特征分析

李希彬1,2，張秋豐1,2，姚志剛3，葉風(fēng)娟1,2，李玉杰1,2

（1.國(guó)家海洋局天津海洋環(huán)境監(jiān)測(cè)中心站，天津300451；2.天津市海洋環(huán)境監(jiān)測(cè)預(yù)報(bào)中心，天津300451；3.中國(guó)海洋大學(xué)海洋環(huán)境學(xué)院，山東青島 266100）

對(duì)塘沽海洋環(huán)境監(jiān)測(cè)站從1950～2008年，59年的潮汐資料進(jìn)行調(diào)和分析，分析了其分潮調(diào)和常數(shù)的變化曲線，并利用FFT譜分析方法對(duì)其調(diào)和常數(shù)的變化周期和原因進(jìn)行了分析；之后應(yīng)用FFT譜分析方法對(duì)去除天文潮后的余水位進(jìn)行分解，分析了近50多年來年平均余水位的多層次周期分布，進(jìn)而利用最小二乘法進(jìn)行線性分析，分析了天津近海平均海平面的變化趨勢(shì)。結(jié)果表明，天津沿岸潮汐M2分潮振幅變化存在20年、4～6年、2～3年的周期，且明顯受到改革開放后工程建設(shè)的影響；余水位變化包含著20年、5～6年、2～3年、1年等不同周期的時(shí)間尺度變化，天津近海海平面呈逐漸上升趨勢(shì)，上升幅度約為3.4 mm/a，主要受海平面上升和地面沉降兩部分因素影響。

天津；潮汐；調(diào)和分析；譜分析；海平面

1 資料介紹

潮汐是海洋活動(dòng)中最為顯著的動(dòng)力現(xiàn)象，與人們的生產(chǎn)、生活及軍事活動(dòng)都有著密切的關(guān)系，也是影響海上交通運(yùn)輸、海洋工程施工、海港安全作業(yè)以及航道設(shè)計(jì)維護(hù)的重要因素，它很早就引起了人類的重視。隨著對(duì)海洋的開發(fā)和利用不斷深化，人類對(duì)海洋的依賴性也越來越大，潮汐的觀測(cè)、分析和利用，一直是人類關(guān)注的重點(diǎn)領(lǐng)域，更是科學(xué)考察研究的熱門課題。

塘沽作為京、津的海上門戶，位于渤海灣畔，隸屬于中國(guó)四大直轄市之一的天津市，從1949年3月15日開始，天津市港務(wù)局在六米碼頭建井驗(yàn)潮，并從1950年1月1日起正式驗(yàn)潮，1986年7月1日～1994年8月4日，由塘沽海洋站在天津新港20段與21段碼頭交匯處建臨時(shí)驗(yàn)潮站，1994年8月5日至今，塘沽海洋站在天津港東突堤工作船碼頭（38°59'N,117°47'E），建永久性驗(yàn)潮站。從1950年1月1日起開始記錄驗(yàn)潮數(shù)據(jù)，距今已有59年的驗(yàn)潮歷史，有多年的系統(tǒng)、連續(xù)、高精度的潮位資料，中間驗(yàn)潮井雖有搬遷，但搬遷位置始終位于天津港內(nèi)，且相距不遠(yuǎn)，潮汐性質(zhì)變化不大，本文認(rèn)為這59年的潮汐數(shù)據(jù)基本可以代表天津近海59年的潮汐變化過程，對(duì)科研和工程開發(fā)均具有較高的利用價(jià)值。

目前已經(jīng)有許多關(guān)于近海潮汐特征分析和研究的方法，歐素英等[1]曾經(jīng)用小波變換的方法研究了廣東近海的海平面周期變化，進(jìn)而研究了廣東沿海的海平面變化趨勢(shì)；楊清書等[2]用低通數(shù)字濾波的方法研究了海平面的變化趨勢(shì)；王曉東等[3]用單站潮汐預(yù)報(bào)的譜分析方法對(duì)廈門和東山的潮汐資料進(jìn)行分析；范文靜等[4]對(duì)塘沽海洋站水位變化及其極值與海面風(fēng)的相關(guān)關(guān)系進(jìn)行了研究；但關(guān)于天津近海的潮汐資料周期性變化的分析和海平面趨勢(shì)變化方面的成果較少，本文的目的即通過對(duì)天津近海59年左右的潮位觀測(cè)數(shù)據(jù)的分析研究，以期能夠?qū)μ旖蚪５某奔捌渥兓卣饔懈M(jìn)一步的認(rèn)識(shí)，同時(shí)對(duì)海平面的長(zhǎng)期變化趨勢(shì)進(jìn)行研究，來揭示近59年來渤海海域的海平面變化趨勢(shì)。

2 潮位資料分析

首先，本文對(duì)塘沽海洋站1950～2008年，共59年的潮汐資料進(jìn)行統(tǒng)計(jì)分析，統(tǒng)計(jì)結(jié)果顯示：塘沽海洋環(huán)境監(jiān)測(cè)站觀測(cè)的潮汐類型數(shù)大約為0.6，屬于不規(guī)則半日潮海區(qū)，其年最高高潮位為581 cm，出現(xiàn)在1992年9月1日，最低潮位-108 cm，出現(xiàn)在1957年12月18日出現(xiàn)，年平均高潮位為374 cm，年平均低潮位為134 cm，年平均潮差為240 cm，年最大潮差為437 cm，出現(xiàn)在1980年10月。

2.1 調(diào)和分析

已有的研究表明，渤海灣是一個(gè)潮汐能量明顯占優(yōu)的海區(qū)，其潮汐信號(hào)以經(jīng)由渤海海峽傳入的協(xié)振潮為主，而由本海區(qū)引潮力所導(dǎo)致的強(qiáng)迫潮基本可以忽略不計(jì)。因此，對(duì)塘沽海洋站海平面長(zhǎng)期變化的研究必須包括對(duì)該站位潮波信號(hào)的時(shí)間演變過程的研究。

通過對(duì)塘沽海洋站59年的時(shí)間序列潮汐數(shù)據(jù)進(jìn)行逐年的調(diào)和分析[5～6]，可以得到一個(gè)其調(diào)和常數(shù)的時(shí)間序列，本文選取了該站位最主要的六個(gè)分潮來進(jìn)行分析，分別是半日分潮M2和S2，全日分潮K1和O1以及淺水分潮M4和MS4，結(jié)果見圖1～12：

其中：

（1）半日分潮M2和S2分潮的振幅以及兩者的遲角之間均存在較相似的變化趨勢(shì)，其中，M2和S2分潮的振幅從1960～1970年之后至2000年，大體呈逐漸增大的趨勢(shì)，2000年之后，其振幅逐漸減小，中間有小的震蕩，M2分潮振幅最大值出現(xiàn)在1964年，為117 cm，最小值出現(xiàn)在2003年，為97 cm，S2振幅最大值出現(xiàn)在1969年，為35.2 cm，最小值出現(xiàn)在1999年，M2和S2分潮的遲角從50年代開始，大體呈現(xiàn)逐漸增大的趨勢(shì)，至70～80年增至最大，其中，M2分潮遲角最大值出現(xiàn)在1981年，為93.2°，S2分潮遲角最大值出現(xiàn)在1978年，為168.9°。1980年之后，M2和S2分潮的遲角呈逐漸減小的趨勢(shì)，M2和S2分潮遲角最小值分別為81.7°和152.6°，均出現(xiàn)在2000年；

圖5 K1分潮振幅變化圖 圖6 K1分潮遲角變化圖

圖7 O1分潮振幅變化圖 圖8 O1分潮遲角變化圖

圖9 M4分潮振幅變化圖 圖10 M4分潮遲角變化圖

圖11 MS4分潮振幅變化圖 圖12 MS4分潮遲角變化圖

（2）全日分潮K1和O1分潮振幅和遲角，均呈現(xiàn)大約1年期的震蕩趨勢(shì)，其中，K1分潮振幅維持在35 cm左右，最大值36.9 cm，出現(xiàn)在1971年，O1分潮振幅維持在27 cm左右，最大值29.5 cm，同樣出現(xiàn)在1971年；K1和O1分潮遲角在1980年前均維持震蕩趨勢(shì)，K1分潮遲角維持在154°左右，O1分潮遲角維持在100°左右，其中，K1分潮遲角最大值出現(xiàn)在1960年，為156.8°，O1分潮遲角最大值出現(xiàn)在1977年，為102.3°，1980年后兩者均呈震蕩減小的趨勢(shì)，K1分潮最小值出現(xiàn)在2008年，為149.9°，O1分潮最小值出現(xiàn)在2003年，為94.5°；

（3）淺水分潮M4和MS4分潮的振幅以及兩者的遲角均相似變化趨勢(shì)，兩個(gè)淺水分潮的振幅從1950年開始，至1980年間，均呈現(xiàn)震蕩增加的趨勢(shì)，1980年后，呈現(xiàn)震蕩減小的趨勢(shì)，其中，M4分潮振幅最大值出現(xiàn)在1978年，為10.3 cm，最小值出現(xiàn)在2005年，為6.8 cm，MS4分潮振幅最大值出現(xiàn)在1968年，為5.95 cm，最小值出現(xiàn)在2005年，為3.9 cm，M4和MS4分潮的遲角從1950年開始，呈震蕩減小的趨勢(shì)，M4和MS4分潮遲角的最大值均出現(xiàn)在1961年，分別為82.6°和158.9°，最小值均出現(xiàn)在2008年，分別為51.5°和126.2°。

對(duì)潮汐能量的計(jì)算表明，M2分潮的勢(shì)能（以調(diào)和常數(shù)中振幅的平方表示勢(shì)能）占全部分潮能量總和的70%左右，并且其振幅隨時(shí)間的變化幅度也接近20 cm，遠(yuǎn)大于其他分潮。因此，本文選取M2分潮作為代表，通過譜分析方法對(duì)其年振幅隨時(shí)間變化的變化特征進(jìn)行進(jìn)一步的分析，結(jié)果見圖13。

圖13 M2分潮振幅變化譜

如圖所示，M2分潮的振幅主要存在20年、4～6年，2～3年左右的周期。因?yàn)閿?shù)據(jù)長(zhǎng)度只有59年，而譜分析一般要對(duì)原始fft的結(jié)果進(jìn)行頻域平均以提高可信度，這導(dǎo)致20年周期的譜不僅包含了該周期附近信號(hào)的能量，還反映了其他周期大于20年的低頻信號(hào)的能量，因此在20年左右的周期附近其譜分析的結(jié)果出現(xiàn)了最大值，其中就包括由于升交點(diǎn)向西移動(dòng)而導(dǎo)致的18.61年的長(zhǎng)周期變化，以及海平面的其他低頻和長(zhǎng)期變化。此外，信號(hào)中還包括明顯的4～6年以及2～3年周期的變化，這個(gè)周期與ENSO現(xiàn)象2～7年的變化周期相似，一定程度上反映了全球大氣強(qiáng)迫的周期變化對(duì)潮汐信號(hào)的影響[1～2]。另外，渤海灣的潮波信號(hào)也明顯的受黃河以及海河入海徑流量變化的影響，王海龍等[7]曾指出，黃河徑流量分別存在著2～3年和4～5年左右變化周期，徑流量的這種周期變化必然會(huì)導(dǎo)致天津近海潮波信號(hào)發(fā)生相應(yīng)變化；黃河、海河等沖淤帶來的岸線和地形的變化也會(huì)引起潮波信號(hào)的變異，另外天津近海近年來，大范圍的圍海造地也會(huì)引起潮波信號(hào)的變化，從而導(dǎo)致其潮汐分潮振幅和遲角的變化，運(yùn)用最小二乘法擬合方法計(jì)算得出：1980年前，M2分潮的振幅呈逐漸增加的趨勢(shì)，速度約為0.068 cm/a，而1980年之后，M2分潮的振幅呈迅速減小的趨勢(shì)，減小的速度約為0.345 cm/a，同時(shí)，隨著圍海造地和港口建設(shè)的進(jìn)行，大面積灘涂被填，港口航道疏浚開發(fā)，對(duì)淺水分潮影響也較大，1980年前，淺水分潮M4和MS4的振幅均呈增加的趨勢(shì)，增加的速度分別為0.075 cm/a和0.029 cm/a，而1980年之后，兩個(gè)淺水分潮的振幅均呈現(xiàn)減小的趨勢(shì)，減小的速度分別為0.058 cm/a和0.015 cm/a，以上變化其中固然部分反映了上文譜分析得出的周期性變化的影響，但也與改革開放之后，天津近海的人為建設(shè)有密切關(guān)系，圍海造地使天津近海的地理環(huán)境出現(xiàn)很大變化，造成天津近海的潮汐特征也出現(xiàn)了較大改變。

2.2 余水位分析

前文已經(jīng)指出，天津近海海平面的變化主要由外海傳入的潮波所導(dǎo)致。另外，其他外界強(qiáng)迫如局地的大氣強(qiáng)迫也會(huì)引起余水位的異常變化，其變化幅度有時(shí)甚至接近或者超過0.5 m，因此，有必要對(duì)這些余水位異常增減水及其時(shí)間變化進(jìn)行分析研究。借助調(diào)和分析的結(jié)果，我們可以從原始觀測(cè)中剔除掉潮波所導(dǎo)致的海平面變化，得到其他非潮因素所導(dǎo)致的海平面變化的時(shí)間序列，即余水位，結(jié)果見圖14。

如圖所示，天津近海余水位最大值為26.7 cm，出現(xiàn)在2006年11月，最小值為-44.2 cm，出現(xiàn)在1967年5月。從圖中我們可以看到余水位的時(shí)間序列存在著明顯的變化，有大的周期性的起伏，也存在著小幅的震蕩，為了更好的分析天津近海余水位的周期性變化，我們對(duì)上述余水位時(shí)間序列進(jìn)行能譜分析，結(jié)果見圖15。

與圖13類似，周期為20年左右的譜中反映了其他極低頻信號(hào)的能量以及海平面的長(zhǎng)期變化，我們予以忽略。從圖15上看，天津近海余水位的變化周期主要包括5～6年、2～3年，1年等不同周期的時(shí)間尺度變化，其中5～6年左右周期信號(hào)比較強(qiáng)，次之為1年左右和2～3年的周期。余水位變化主要取決于外界大氣強(qiáng)迫，尤其是風(fēng)場(chǎng)和氣壓場(chǎng)，而余水位中5～6年的周期正是在一定程度上反映了ENSO等外界大氣強(qiáng)迫的變化，大氣強(qiáng)迫的變化會(huì)導(dǎo)致風(fēng)場(chǎng)和氣壓場(chǎng)發(fā)生相應(yīng)變化，從而對(duì)該海區(qū)的余水位的變化產(chǎn)生影響，此外，太陽黑子的周期性變化也可能是影響天津近海余水位變化的因素之一。余水位中還存在著明顯的1年左右的變化周期，該周期反映了季風(fēng)躍遷的影響，夏季盛行的偏南季風(fēng)導(dǎo)致中國(guó)近海北高南低，冬季正好相反。而2～3年的周期反映了與該海區(qū)附近河流徑流量變化的影響，如前文所指出的，黃河等存在著2～3年的變化周期，因而在余水位中可以看到徑流量周期變化對(duì)海平面的影響。另外，極地海冰和陸冰的長(zhǎng)期變化以及氣候的長(zhǎng)期變化也會(huì)導(dǎo)致海平面產(chǎn)生變化。

2.3 海平面變化分析

海平面的變化直接影響制約沿海地區(qū)人民生活和社會(huì)經(jīng)濟(jì)的發(fā)展，海平面的上升將進(jìn)一步影響到這些地區(qū)的安全和經(jīng)濟(jì)發(fā)展。雖然海平面上升是一個(gè)相對(duì)緩慢的過程，但長(zhǎng)期積累的結(jié)果足以對(duì)沿海經(jīng)濟(jì)發(fā)展、城市安危和人民生活環(huán)境帶來多方面的不利影響，這種影響比任何一種自然災(zāi)害都要廣泛和深入。因此對(duì)海平面變化趨勢(shì)的研究和預(yù)測(cè)，具有十分重要的意義[8]。

一般來說，年平均海平面時(shí)間序列會(huì)隨年際變化表現(xiàn)出某種趨勢(shì)性，即認(rèn)為平均海平面 y與時(shí)間 t有如下關(guān)系：y(t)=a+bt。利用年均水位變化序列，通過最小二乘法擬合來確定天津近海的海平面變化率，結(jié)果見圖16。

從圖中可以看出，天津近海的平均海平面是在緩慢上升的，經(jīng)過計(jì)算得到天津近海平均海平面上升的速度約為3.4 mm/a，這種上升的變化，主要是由兩部分因素組成的。一方面是由于海平面的上升；另一方面說明天津港地區(qū)的地面在逐漸下沉，特別是在1986年驗(yàn)潮站搬遷之后，位于吹填區(qū)域的新建碼頭，局部的下沉量大應(yīng)當(dāng)是合理的。1986年以前的海平面變化情況是比較平穩(wěn)的，1986年以后則變化比較劇烈，正是說明了上述變化情況。分別對(duì)1950～1986年和1986～2008年的潮汐數(shù)據(jù)進(jìn)行平均和擬合，發(fā)現(xiàn)：1950～1986年的平均海平面為225.422 cm，平均海平面上升的速度基本為1.4 mm/a，而1986～2008年的平均海平面則達(dá)到了236.923 cm，平均海平面上升的速度基本為4.7 mm/a，相比1950～1986年，平均海平面高出11.501 cm，平均海平面上升的速度高出3.3 mm/a，也說明了天津港地區(qū)處于逐漸沉降中，當(dāng)然，以上結(jié)論是在只有天津港一個(gè)水準(zhǔn)點(diǎn)的水準(zhǔn)資料的情況下做出的，若要論證整個(gè)天津近海地區(qū)的下沉情況，還需要對(duì)周邊大面積地區(qū)的水準(zhǔn)測(cè)量結(jié)果進(jìn)行綜合比較才能得出結(jié)論。

圖14 余水位月均序列表

圖15 余水位變化譜

圖16 天津近海海平面變化趨勢(shì)

3 結(jié)論

綜合以上分析，可以得到以下結(jié)論：

（1）通過對(duì)塘沽海洋站每年的潮汐數(shù)據(jù)進(jìn)行調(diào)和分析，得到關(guān)于分潮調(diào)和常數(shù)的時(shí)間序列，分析了各主要調(diào)和常數(shù)的變化趨勢(shì)，并對(duì)主要影響天津近海的M2分潮的振幅進(jìn)行FFT譜分析，分析其能量變化周期及產(chǎn)生原因；

（2） 用譜分析的方法研究了天津近海余水位在時(shí)域-頻域中的周期分布，很好的揭示了天津余水位的周期變化特征，發(fā)現(xiàn)其變化周期主要包含20年、5～6年、1年等時(shí)間尺度，其中，以5～6年的時(shí)間變化為主，之后分析了產(chǎn)生這些變化周期的可能因素；

（4）針對(duì)近年來較為關(guān)注的海平面問題，本文采用最小二乘法擬合的方法，線性擬合了天津近海的海平面變化率，結(jié)果顯示天津近海海平面呈逐年升高的趨勢(shì)，上升的速度約為3.4 mm/a，并就天津港的平均海平面上升現(xiàn)象從海平面上升和地面沉降兩方面進(jìn)行了分析。

[1]歐素英,陳子燊.小波變換在相對(duì)海平面變換研究中的應(yīng)用[J].地理科學(xué),2004,24(3):358-364.

[2]楊清書,吳超羽.低通數(shù)字濾波在確定海平面變化趨勢(shì)中的應(yīng)用[J].海洋通報(bào),1996,15(1):7-14.

[3]王曉東,蔣國(guó)榮,王英俊.單站潮汐的譜分析預(yù)報(bào)[J].海洋預(yù)報(bào),2008,25(3):53-59.

[4]范文靜,張愛軍.塘沽海洋站水位變化及其極值與海面風(fēng)的相關(guān)關(guān)系研究[J].海洋通報(bào).2006,25(4):29-35.

[5]黃祖珂,黃磊.潮汐原理與計(jì)算[M].青島:中國(guó)海洋大學(xué)出版社,2005:12-88.

[6]王長(zhǎng)海.潮汐調(diào)和分析的一種模式[J].海洋預(yù)報(bào),l995,l2(4):7l-76.

[7]王海龍,李國(guó)勝.近50年來黃河入海水沙通量變化的多尺度分析[P].自然科學(xué)進(jìn)展,2006,16(12):1639-1644.

[8]左軍成,陳宗鏞,周天華.中國(guó)沿岸相對(duì)海面變化的本征分析和預(yù)測(cè)[J].海洋與湖沼,1995,26(3):331-337.

P731

A

1003-0239（2011）02-0041-08

2010-07-15

李希彬(1983-),男,碩士,主要從事海洋預(yù)報(bào)工作。 E-mail:lixb_tj@yahoo.com.cn