王文才，李一平 ，杜 薇，曾偉峰，許益新

(1.河海大學(xué)環(huán)境學(xué)院，江蘇 南京 210098； 2.河海大學(xué)淺水湖泊綜合治理與資源開發(fā)教育部重點(diǎn)實(shí)驗(yàn)室，江蘇 南京 210098)

長(zhǎng)江感潮河段潮汐變化特征

王文才1,2，李一平1,2，杜 薇1,2，曾偉峰1,2，許益新1,2

(1.河海大學(xué)環(huán)境學(xué)院，江蘇 南京 210098； 2.河海大學(xué)淺水湖泊綜合治理與資源開發(fā)教育部重點(diǎn)實(shí)驗(yàn)室，江蘇 南京 210098)

基于2011年長(zhǎng)江下游大通水文站的逐日流量以及南京、鎮(zhèn)江、江陰、天生港、徐六涇、共青圩6個(gè)潮位站的逐日高潮、低潮資料，研究海洋潮汐從河口向河段上游傳播過(guò)程中的潮汐特征變化規(guī)律。結(jié)果表明，各潮位站25 h滑動(dòng)平均的潮位變化與流量呈現(xiàn)相同的變化趨勢(shì)，南京和鎮(zhèn)江兩個(gè)潮位站的水位變化最為明顯，江陰、天生港、徐六涇的水位變化相對(duì)較小，共青圩潮位站的潮位基本不隨徑流發(fā)生相應(yīng)的變化趨勢(shì)；各潮位站之間均有一定的相關(guān)性，且相鄰兩個(gè)潮位站之間的水位相關(guān)性較高；隨著潮位站之間距離的增大，水位的相關(guān)性呈現(xiàn)出明顯的下降趨勢(shì)；低水位的傳播延遲時(shí)間均大于高水位的傳播延遲時(shí)間；隨著站點(diǎn)與入?？诰嚯x的增大，漲潮歷時(shí)出現(xiàn)逐漸減小、落潮歷時(shí)逐漸增大的變化趨勢(shì)。

徑流；潮汐；感潮河段；預(yù)報(bào)

地球上數(shù)量眾多的河流與海洋直接相連，這些河流在靠近河口的區(qū)段受到下游海洋天文潮汐和上游徑流量的共同作用，水動(dòng)力變化過(guò)程比一般的內(nèi)陸河流更加復(fù)雜[1-3]。潮波在感潮河段向上游傳播時(shí)會(huì)受到底部摩擦，因河道地形寬深發(fā)生變化以及上游徑流作用而發(fā)生改變[4-7]。揭示和理解外海潮汐通過(guò)河口向河流上游傳播過(guò)程中的水動(dòng)力變化規(guī)律，對(duì)感潮河段區(qū)間內(nèi)的潮汐預(yù)報(bào)及河道、河口生態(tài)系統(tǒng)的管理具有重要意義[8-9]，也可以為區(qū)域內(nèi)的防洪、灌溉、航運(yùn)、引排水等提供決策支持。

圖1 研究區(qū)域

長(zhǎng)江下游是典型的感潮河段，且該河段主流上并沒(méi)有阻隔水流的大型水利工程設(shè)施，故可以很好地讓人們研究潮汐沿河口向上游傳播過(guò)程中的變化特征。國(guó)內(nèi)有眾多的學(xué)者通過(guò)各種方法分析長(zhǎng)江下游各潮位觀測(cè)站的潮位變化規(guī)律，并給出觀測(cè)站所在位置的潮位預(yù)測(cè)方法。林勛勵(lì)[10]對(duì)1948—1980年大通水文站的流量以及長(zhǎng)江下游江陰、天生、吳淞各潮位站的水位資料進(jìn)行分析，發(fā)現(xiàn)江陰和天生港受徑流影響而存在高水位期和低水位期，并通過(guò)建立大通水文站月平均流量與各潮位站月平均水位的相關(guān)關(guān)系，對(duì)各潮位站逐月、逐旬的潮汐預(yù)報(bào)進(jìn)行了修正。李國(guó)芳等[11]基于長(zhǎng)江口天生港、徐六涇、楊林3個(gè)潮位站的潮位以及上游大通站的日均流量資料，建立了各潮位站日均增水與大通流量增量之間的關(guān)系，并利用所得關(guān)系式，根據(jù)上游流量的大小，對(duì)潮汐預(yù)報(bào)模型計(jì)算的過(guò)去若干年的潮位進(jìn)行校正，使日均潮位的預(yù)報(bào)精度得到了明顯的提高。宋正逵等[12]采用傳統(tǒng)的海洋潮汐調(diào)和分析方法對(duì)長(zhǎng)江下游天生港潮位站的潮位進(jìn)行預(yù)報(bào)，發(fā)現(xiàn)精度不高，經(jīng)分析后建立了潮位預(yù)報(bào)誤差與同期流量差值的關(guān)系，為天生港站潮位實(shí)時(shí)預(yù)報(bào)修正提供了依據(jù)。但這些研究大多注重對(duì)相關(guān)站點(diǎn)的潮位預(yù)測(cè)結(jié)果，而缺乏對(duì)潮汐沿程及潮汐本身特征的細(xì)致分析。本文以2011年長(zhǎng)江下游大通水文站的逐日流量以及南京、鎮(zhèn)江、江陰、天生港、徐六涇、共青圩6個(gè)潮位站的逐日高潮、低潮資料為基礎(chǔ)，細(xì)致研究海洋潮汐從河口向河段上游傳播過(guò)程中的潮汐特征變化規(guī)律，以期為下游潮汐預(yù)報(bào)提供參考，為水利工程的優(yōu)化調(diào)度提供決策支持。

1 研究區(qū)域概況

長(zhǎng)江的感潮段通常指大通至入?？谶@一河段，其受上游徑流和下游潮汐的共同作用(圖1)。一般認(rèn)為潮區(qū)界(水位隨潮汐漲落的上界)在大通附近，并隨著徑流量、潮汐的變化而有所變化[13]，潮流界(水流所能達(dá)到的上界)位置隨徑流的變化而出現(xiàn)較大的擺動(dòng)[14-16]。大通水文站作為下游流量控制站，流量隨季節(jié)發(fā)生變化，多年平均流量為28 700 m3/s，歷年最大徑流量為92 600 m3/s，最小徑流量4 260 m3/s[17]。下游南京、鎮(zhèn)江、江陰、天生港、徐六涇、共青圩6個(gè)潮位站與大通水文站的距離分別為220、320、405、458、490、586 km[18]。本文通過(guò)這幾個(gè)站的潮位資料來(lái)研究潮汐變化特征。

2 材料與方法

基于2011年大通水文站逐日的日均徑流量以及下游南京、鎮(zhèn)江、江陰、天生港、徐六涇、共青圩6個(gè)潮位站的高潮、低潮潮位及發(fā)生時(shí)刻的資料，采用cubic插值方法得到逐時(shí)潮位，另外考慮到一個(gè)太陰日為24 h 50 min，計(jì)算25 h滑動(dòng)平均潮位，這可以過(guò)濾M2、K2等高頻潮波，從而保留一些相對(duì)低頻、周期較長(zhǎng)的潮波[19]。這樣處理有助于反映徑流影響下的水位變化特征。計(jì)算各潮位站滑動(dòng)平均水位之間的相關(guān)系數(shù)，就可得到相應(yīng)的相關(guān)系數(shù)矩陣，以便分析各潮位站之間的關(guān)聯(lián)程度。筆者以主要受潮汐影響的共青圩潮位站為參考，統(tǒng)計(jì)分析潮差、高低潮延時(shí)等的沿程變化，并給出相應(yīng)的擬合公式。

3 結(jié)果與討論

3.1 徑流及水位特征

2011年全年呈現(xiàn)出偏枯情況，大通年初徑流量為15 300 m3/s。前5個(gè)月徑流量基本穩(wěn)定在15 000 m3/s左右，從六月份開始徑流量呈現(xiàn)逐漸增大的趨勢(shì)，并于6月22日達(dá)到峰值46 000 m3/s，之后徑流量逐漸下降至年末的12 400 m3/s，當(dāng)年平均徑流量為21 143 m3/s(圖2)。各潮位站25 h滑動(dòng)平均的潮位變化與流量呈現(xiàn)相同的變化趨勢(shì)，南京和鎮(zhèn)江兩個(gè)潮位站的水位變化最為明顯，江陰、天生港、徐六涇的水位變化相對(duì)較小，共青圩潮位站的潮位基本不隨徑流發(fā)生相應(yīng)的變化趨勢(shì)(圖3)。

圖2 2011年逐日的日均流量變化過(guò)程

圖3 各潮位站滑動(dòng)平均水位變化過(guò)程

分析各潮位站之間25 h滑動(dòng)平均水位之間的相關(guān)性，分析結(jié)果見表1。從表1可以看出，各潮位站之間均有一定的相關(guān)性，滿足顯著性檢驗(yàn)plt;0.05，且相鄰兩個(gè)潮位站之間的水位相關(guān)性較高；隨著潮位站之間距離的增大，水位的相關(guān)性呈現(xiàn)出明顯的下降趨勢(shì)；共青圩潮位站與各潮位站水位之間的相關(guān)性均小于0.9。

表1 各潮位站滑動(dòng)平均水位相關(guān)系數(shù)矩陣

3.2 潮位延時(shí)特征

以共青圩為參考潮位站，統(tǒng)計(jì)各潮位站高潮位、低潮位到達(dá)時(shí)間，并計(jì)算高潮位、低潮位從共青圩站傳播到相應(yīng)站點(diǎn)所需要的時(shí)間，以此作為高低潮位的延遲時(shí)間。從結(jié)果(圖4)可以看出，低水位的傳播延遲時(shí)間均大于高水位的傳播延遲時(shí)間；除南京和鎮(zhèn)江潮位站外，其他潮位站的水位延遲并不隨徑流量的變化而呈現(xiàn)顯著的變化。

計(jì)算全年的高水位、低水位傳播的平均延時(shí)，計(jì)算結(jié)果見表2。從表2可以發(fā)現(xiàn)，與共青圩相比，徐六涇、天生港、江陰、鎮(zhèn)江、南京的低水位延時(shí)分別為3.46、4.71、5.89、9.42和11.09h，高水位的延時(shí)分別為2.87、4.05、4.68、7.91和9.65h。采用二次多項(xiàng)式進(jìn)行擬合，得到高、低水位傳播延時(shí)與相應(yīng)站點(diǎn)距共青圩距離的擬合公式，其具有較好的效果，可以反映高、低水位沿程傳播的過(guò)程(圖5)。

(a) 高水位傳播延時(shí)

(b)低水位傳播延時(shí)

站點(diǎn)名稱距離/km低水位平均延時(shí)/h高水位平均延時(shí)/h共青圩000徐六涇963.462.87天生港1284.714.05江 陰1815.894.68鎮(zhèn) 江2669.427.91南 京36611.099.65

圖5 高低水位傳播的平均延時(shí)隨傳播距離的變化

3.3 漲落潮歷時(shí)及潮差特征

統(tǒng)計(jì)潮位站的漲落潮平均歷時(shí)，可以發(fā)現(xiàn)隨著站點(diǎn)與入?？诰嚯x的增大，漲潮歷時(shí)出現(xiàn)逐漸減小，落潮歷時(shí)逐漸增大的變化趨勢(shì)(表3)，這也符合河口地區(qū)潮汐的變化特征[20]。其中，到最上游南京站與到鎮(zhèn)江站相比，漲潮歷時(shí)并沒(méi)有減小，反而出現(xiàn)微弱上升。這與俞慕耕[21]對(duì)長(zhǎng)江下游潮汐的統(tǒng)計(jì)結(jié)果基本一致，但目前并沒(méi)有相關(guān)研究來(lái)揭示其原因。與前人的統(tǒng)計(jì)結(jié)果對(duì)比，本研究發(fā)現(xiàn)各潮位站的平均潮差要偏大，這是因?yàn)?011年為偏枯水年，徑流量較小，潮汐動(dòng)力作用較強(qiáng)，造成平均潮差相對(duì)較大。

表3 各潮位站漲落潮平均歷時(shí)與平均潮差的情況

對(duì)各站潮位進(jìn)行歸一化處理后建立的相對(duì)潮差與站點(diǎn)距共青圩潮位站距離的二次擬合關(guān)系，可以較好地反映潮差的變化，可以在一定程度上根據(jù)共青圩的潮位預(yù)測(cè)其上游給定位置的潮差情況(圖6)。

圖6 相對(duì)潮差隨距離的變化關(guān)系

4 結(jié) 論

a. 潮汐從共青圩向上游傳播的過(guò)程中，高水位的延遲小于低水位延遲，高、低水位各自的平均延遲時(shí)間可以通過(guò)二次多項(xiàng)式建立其與相應(yīng)位置距共青圩距離的函數(shù)關(guān)系，用于預(yù)測(cè)河段不同位置的高、低潮到達(dá)時(shí)間。

b. 南京、鎮(zhèn)江的高、低水位的延遲時(shí)間受到徑流增大影響而相應(yīng)增長(zhǎng)，江陰及下游潮位站的高、低水位延遲受徑流影響較小。

c. 鎮(zhèn)江以下，長(zhǎng)江下游河段沿入海口向上游的漲潮歷時(shí)逐漸縮短、落潮歷時(shí)逐漸增加；鎮(zhèn)江以上，漲落潮歷時(shí)與下游相比并不出現(xiàn)明顯的縮短和增加，具體原因還需進(jìn)一步研究。

d. 上游各站的相對(duì)潮差可以和其與共青圩的距離建立二次擬合關(guān)系，用于預(yù)測(cè)沿程的相應(yīng)潮差變化。

本文主要是基于2011年水文資料對(duì)長(zhǎng)江感潮段的潮位變化特征進(jìn)行分析研究，結(jié)果可能存在一定的局限性，后續(xù)還需要進(jìn)行更深入完整的研究。

[1]DRONKERS J J.Tidal computations in rivers and coastal waters[M].Amsterdam：North-Holland Pub.Co.,1964.

[2]LEBLOND P H.On tidal propagation in shallow rivers[J].Journal of Geophysical Research:Oceans(1978-2012),1978,83(C9):4717-4721.

[3]JAY D A.Green’s law revisited:tidal long-wave propagation in channels with strong topography[J].Journal of Geophysical Research:Atmospheres,1991,96(C11):20585-20598.

[4]FRIEDRICHS C T,AUBREY D G.Tidal propagation in strongly convergent channels[J].Journal of Geophysical Research:Oceans(1978-2012),1994,99(C2):3321-3336.

[5]HORREVOETS A C,SAVENIJE H,SCHUURMAN J N,et al.The influence of river discharge on tidal damping in alluvial estuaries[J].Journal of Hydrology,2004,294(4):213-228.

[6]王遠(yuǎn)坤,王棟.基于樣本熵理論的長(zhǎng)江干流徑流序列復(fù)雜性分析[J].河海大學(xué)學(xué)報(bào)(自然科學(xué)版),2015,43(3):203-207.(WANG Yuankun,WANG Dong. Analysis of complexity of runoff series based on sample entropy in Yangtze River[J]. Journal of Hohai University (Natural Sciences),2015,43(3):203-207. (in Chinese))

[7]田威,邱利,李一平.基于EFDC模型的感潮江段溢油事故風(fēng)險(xiǎn)預(yù)測(cè)[J].水資源保護(hù),2015,31(6):98-102.(TIAN Wei,QIU Li,LI Yiping. Prediction of oil spill risk on tide sect of Yangtze River based on EFDC model[J].Water Resources Protection,2015,31(6):98-102.(in Chinese))

[8]KUKULKA T,JAY D A.Impacts of Columbia River discharge on salmonid habitat:1.a nonstationary fluvial tide model[J].Journal of Geophysical Research:Atmospheres,2003,108(C9):351-378.

[9]KUKULKA T,JAY D A.Impacts of Columbia River discharge on salmonid habitat:2.Changes in shallow-water habitat[J].Journal of Geophysical Research:Oceans,2003,108:351-378.

[10]林勛勵(lì).長(zhǎng)江口天文潮預(yù)報(bào)修正的一種方法[J].海洋預(yù)報(bào),1985(1):38-46.(LIN Xunli.A method of updating the astronomical tide forecast in the Yangtze River Estuary[J].Marine Forecast Service,1985(1):38-46.(in Chinese))

[11]李國(guó)芳,譚亞,張秀菊.感潮河段上游流量對(duì)潮位預(yù)報(bào)的影響[J].河海大學(xué)學(xué)報(bào)(自然科學(xué)版),2006,34(2):144-147.(LI Guofang,TAN Ya,ZHANG Xiuju.Influence of upstream discharge in tidal level prediction for tidal reaches[J].Journal of Hohai University(Natural Sciences),2006,34(2):144-147.(in Chinese))

[12]宋正逵,黃利亞,仇娟娟.長(zhǎng)江上游徑流對(duì)天生港站潮位影響淺析[J].水文,2010,30(4):56-58.(SONG Zhengkui,HUANG Liya,QIU Juanjuan.Influence of runoff from upper reaches of Yangtze River on tidal level at Tianshenggang Station[J].Journal of Chinese Hydrology,2010,30(4):56-58.(in Chinese))

[13]徐漢興,樊連法,顧明杰.對(duì)長(zhǎng)江潮區(qū)界與潮流界的研究[J].水運(yùn)工程,2012(6):15-20.(XU Hanxing,Fan Lianfa,GU Mingjie.On tidal mark and tidal current mark in the Yangtze River[J].Port amp; Waterway Engineering,2012(6):15-20.(in Chinese))

[14]侯成程,朱建榮.長(zhǎng)江河口潮流界與徑流量定量關(guān)系研究[J].華東師范大學(xué)學(xué)報(bào)(自然科學(xué)版),2013(5):18-26.(HOU Chengcheng,ZHU Jianrong.Study on the quantitative relationship of the location of the tidal current limit and the river discharge in the Changjiang Estuary[J].Journal of East China Normal University(Natural Science),2013(5):18-26.(in Chinese))

[15]楊云平,李義天,韓劍橋,等.長(zhǎng)江口潮區(qū)和潮流界面變化及對(duì)工程響應(yīng)[J].泥沙研究,2012(6):46-51.(YANG Yunping,LI Yitian,HAN Jianqiao,et al.Variation of tide limit and tidal current limit in Yangtze Estuary and its impact on projects[J].Journal of Sediment Research,2012(6):46-51.(in Chinese))

[16]宋蘭蘭.長(zhǎng)江潮流界位置探討[J].水文,2002,22(5):25-26.(SONG Lanlan.Research on the location of tidal current limit of Yangtze River[J].Journal of Chinese Hydrology,2002,22(5):25-26.(in Chinese))

[17]夏云峰,聞云呈,梁婷.長(zhǎng)江下游典型代表性潮型選擇研究初探[C].//水道與港口工程國(guó)際學(xué)術(shù)研討會(huì)論文集.天津：交通運(yùn)輸部天津水運(yùn)工程科學(xué)研究院,2011:501-505.

[18]GUO L,WEGEN M V D,JAY D A,et al.River-tide dynamics:exploration of nonstationary and nonlinear tidal behavior in the Yangtze River estuary:river tidal dynamics[J].Journal of Geophysical Research:Oceans,2015,120(5):3499-3521.

[19]GODIN G.The analysis of tides[M].Toronto,Buffalo:University of Toronto Press,1972:1005-1012.

[20]陳宗鏞.潮汐學(xué)[M].北京:科學(xué)出版社,1980.

[21]俞慕耕.我國(guó)沿海潮汐的特點(diǎn)[C]//第一屆潮汐與海平面學(xué)術(shù)討論會(huì)論文集.天津:海洋局海洋科技情報(bào)研究所,1986.

TidalvariationfeaturesoftidalreachofChangjiangRiver

WANGWencai1,2,LIYiping1,2,DUWei1,2,ZENGWeifeng1,2,XUYixin1,2

(1.CollegeofEnvironment,HohaiUniversity,Nanjing210098,China; 2.KeyLaboratoryofIntegratedRegulationandResourceDevelopmentonShallowLakesofMinistryofEducation,HohaiUniversity,Nanjing210098,China)

Based on the daily discharge of Datong station in the downstream of the Yangtze river and the high and low tidal water level of Nanjing, Zhenjiang, Tianshenggang, Xuliujing, Gongqingwei stations in 2011, this paper studied the tidal characteristics of ocean tides from the estuary to the upper reaches of the river. The results show that the tidal mean change of 25h sliding average of each tidal station shows the same trend as the discharge, the water level changes of two tide stations in Nanjing and Zhenjiang are the most obvious. Water level change of Jiangyin, Tianshenggang, Xuliujing stations is relatively small. The tidal water level of Gongqingwei station does not change with the runoff. There is a certain correlation between the tidal stations, and the correlation between the water level of the two adjacent tidal stations is higher. With the increase of the distance between tidal stations, the correlation of water levels shows a significant downward trend. Propagation delay time of low tidal level is greater than high tidal water level. Furthermore, with the increase of distance between the tidal stations and the estuary, flood tide duration displays the trend of gradual decrease while ebb tide duration presents the trend of gradual increase.

runoff; tide; tidal reach; forecast

10.3880/j.issn.1004-6933.2017.06.19

國(guó)家自然科學(xué)基金(51579071)；創(chuàng)新研究群體科學(xué)基金(51421006)；2016年度江蘇省普通高校專業(yè)學(xué)位研究生實(shí)踐創(chuàng)新計(jì)劃項(xiàng)目(SJZZ16_0089)；江蘇省“創(chuàng)新團(tuán)隊(duì)計(jì)劃”；江蘇高校優(yōu)勢(shì)學(xué)科建設(shè)工程資助項(xiàng)目

王文才(1992—)，男，碩士研究生，研究方向?yàn)榄h(huán)境科學(xué)。E-mail：591069468@qq.com

李一平，教授。E-mail:liyiping@hhu.edu.cn.

P333

A

1004-6933(2017)06-0121-04

2017-01-09 編輯：彭桃英)