河口擋潮閘對(duì)三角洲潮汐不對(duì)稱時(shí)空變化的影響*

沈倩穎, 季小梅, 張蔚, 徐龑文

1. 河海大學(xué), 江蘇省海岸海洋資源開(kāi)發(fā)與環(huán)境安全重點(diǎn)實(shí)驗(yàn)室, 江蘇 南京 210098;

2. 港口航道泥沙工程交通行業(yè)重點(diǎn)實(shí)驗(yàn)室, 江蘇 南京 210029;

3. 河海大學(xué), 疏浚技術(shù)教育部工程研究中心, 江蘇 南京 210098;

4. 上海灘涂海岸工程技術(shù)研究中心, 上海 200061

外海規(guī)則潮波進(jìn)入三角洲內(nèi), 在地形和徑流等因素影響下, 生成淺水分潮, 潮波也發(fā)生變形, 使得漲、落潮在振幅、歷時(shí)、流速上不再對(duì)稱(Dronkers,1986)。潮汐和潮流這種變形在三角洲地區(qū)決定了泥沙的凈輸移趨勢(shì), 影響河床穩(wěn)定性, 從而對(duì)三角洲的地貌演變方向產(chǎn)生重要影響(Speer et al, 1985;Aldridge, 1997; Wang et al, 2002)。三角洲內(nèi)的潮汐特性與海洋和陸地邊界的改變密切相關(guān)(Vellinga et al, 2014), 目前研究集中在氣象條件改變下的海平面抬升(Woodworth et al, 2011)和長(zhǎng)期人類活動(dòng), 如河道采砂疏浚(Zhang et al, 2018)導(dǎo)致的河床下切,上游大壩修建(Yu et al, 2020)對(duì)徑流的季節(jié)性調(diào)節(jié),對(duì)三角洲內(nèi)的潮波特性產(chǎn)生的影響。

擋潮閘作為河口地區(qū)重要的防洪擋潮工程極大地改變了三角洲海洋邊界條件, 對(duì)三角洲內(nèi)部潮汐動(dòng)力和潮汐不對(duì)稱性有著重要影響。荷蘭三角洲大部分陸地低于海平面(低至-6m), 1953年特大洪水后,荷蘭在西南部海岸成功建設(shè)一系列三角洲擋潮閘工程, 抵御北海風(fēng)暴潮。其中, 1970年哈靈河口擋潮閘壩對(duì)萊茵河-默茲河三角洲內(nèi)部潮汐動(dòng)力環(huán)境影響最為顯著。T?nis等(2002)利用多年長(zhǎng)期測(cè)深數(shù)據(jù)研究哈靈河口關(guān)閉后口門(mén)外地形演變, 研究發(fā)現(xiàn)工程建設(shè)后擋潮閘外發(fā)生泥沙淤積, 且地形適應(yīng)時(shí)間尺度約為11年。Vellinga等(2014)研究萊茵河-默茲河三角洲70年歷史水位數(shù)據(jù), 發(fā)現(xiàn)人類活動(dòng)對(duì)河網(wǎng)內(nèi)極端高、低水位的影響超過(guò)了平均海平面上升產(chǎn)生的影響。而工程建設(shè)對(duì)三角洲內(nèi)部河網(wǎng)潮波特性的變化則沒(méi)有深入的研究討論。本文主要利用長(zhǎng)時(shí)間水位數(shù)據(jù)研究河口擋潮閘修建對(duì)潮汐漲落潮歷時(shí)不對(duì)稱性的影響, 借助非平穩(wěn)潮汐調(diào)和分析方法分解萊茵河-默茲河三角洲內(nèi)各站點(diǎn)潮位數(shù)據(jù), 并根據(jù)傳統(tǒng)潮汐不對(duì)稱方法定量分析整個(gè)潮汐河網(wǎng)的潮動(dòng)力和潮汐不對(duì)稱特性在河口擋潮閘壩工程建設(shè)前后的時(shí)空演變規(guī)律, 從而了解局部的河口工程對(duì)潮汐河網(wǎng)整個(gè)動(dòng)力體系產(chǎn)生的影響。

1 數(shù)據(jù)和方法

1.1 研究區(qū)域

萊茵河-默茲河三角洲位于荷蘭西南部, 由萊茵河、默茲河的復(fù)雜支流組成, 萊克河和瓦爾河為萊茵河下游的兩條分汊支流(圖1)。兩大河流最終經(jīng)兩條入海通道流入北海, 分別是北部新水道和南部哈靈水道, 擋潮閘建成后流入北海的平均流量分別為1960m3·s–1、230m3·s–1(Ysebaert et al,2016), 三角洲面積約7500km2。入?？陂T(mén)外平均潮差在2.4m左右, 潮型表現(xiàn)為正規(guī)半日潮特征(T?nis et al, 2002)。

圖1 萊茵河-默茲河三角洲及潮位站分布圖圖b中數(shù)字表示站點(diǎn)序號(hào); 方形塊、圓形塊分別表示站點(diǎn)實(shí)測(cè)潮位數(shù)據(jù)時(shí)間長(zhǎng)度為1961—2018年、1971—2018年。該圖基于國(guó)家測(cè)繪地理信息局標(biāo)準(zhǔn)地圖服務(wù)網(wǎng)站下載的審圖號(hào)為GS(2020)4390的標(biāo)準(zhǔn)地圖制作Fig. 1 Rhine–Meuse Delta and hydrological stations. The numbers in figure b represent the station number. The square blocks and round blocks indicate that the measured water level data of stations are from 1961—2018 and 1971—2018,respectively

1970年在哈靈水道入?？陂T(mén)處建成哈靈水道擋潮閘壩, 長(zhǎng)4.5km, 有17孔泄水閘, 是荷蘭三角洲工程的第六座建筑物。該工程一方面抵御北海的風(fēng)暴潮, 另一方面隨著潮位和上游徑流量變化開(kāi)閉排水閘門(mén), 保證北部新水道徑流穩(wěn)定、航運(yùn)安全。南部哈靈水道和荷蘭深河在工程建設(shè)后由潮汐通道變?yōu)榕c北海隔離的淡水盆地, 通過(guò)連接北部入海通道的中部支流傳播潮動(dòng)力(Hoitink et al, 2017), 潮差由閉合前的2m減小到最大僅20cm (Antheunisse et al,2007)。此外擋潮閘壩改變了三角洲河網(wǎng)上游徑流分配, 北部和中部通道內(nèi)徑流量增大(Vellinga et al,2014), 所以擋潮閘建立對(duì)三角洲內(nèi)河網(wǎng)系統(tǒng)的徑潮動(dòng)力產(chǎn)生了很大影響, 進(jìn)而造成了河口動(dòng)力地形大范圍、長(zhǎng)時(shí)間適應(yīng)變化, 主要表現(xiàn)為在連接南北通道的中部支流內(nèi)發(fā)生沖刷, 在南部通道發(fā)生淤積。但地形的變化主要還是由長(zhǎng)期疏?；顒?dòng)主導(dǎo)(Sloff et al, 2013; Vellinga et al, 2014), 且相較于擋潮閘建立導(dǎo)致三角洲動(dòng)力的突然變化, 地形變化產(chǎn)生的影響較緩慢。

1.2 數(shù)據(jù)

潮位站點(diǎn)實(shí)測(cè)水位數(shù)據(jù)和上游控制站點(diǎn)流量數(shù)據(jù)從荷蘭基礎(chǔ)設(shè)施和環(huán)境部網(wǎng)站(Rijkswaterataat)下載(https://www.rijkswaterstaat.nl)。從20世紀(jì)開(kāi)始荷蘭開(kāi)展高精度水位測(cè)量工作, 本文收集萊茵河-默茲河三角洲河網(wǎng)13個(gè)潮位站點(diǎn)的50～60年長(zhǎng)期水位數(shù)據(jù), 實(shí)測(cè)水位數(shù)據(jù)1970年前為間隔3h數(shù)據(jù), 1970年后為間隔1h數(shù)據(jù), 水位基面為阿姆斯特丹高程基準(zhǔn)面(Normaal Amsterdams Peil)。流量數(shù)據(jù)選取三角洲上游3個(gè)控制站點(diǎn)流量, 圖2表示從北至南萊克河、瓦爾河及默茲河上游1961—2018年月均徑流量。和徑流動(dòng)力, 與海洋中規(guī)則潮汐有很大不同, 呈現(xiàn)顯著的非線性變形和不對(duì)稱性特征。傳統(tǒng)的調(diào)和分析工具T_TIDE (Pawlowicz et al, 2002)基于潮位不包含徑流等非穩(wěn)態(tài)信號(hào)的假設(shè), 在河口上游地區(qū)反演潮動(dòng)力誤差較大。本文選用Matte等(2013)年提出的非平穩(wěn)潮汐調(diào)和分析NS_TIDE方法, 能夠在保存調(diào)和分析優(yōu)勢(shì)的基礎(chǔ)上, 很好地解決徑流影響下非平穩(wěn)潮汐信號(hào)中主要分潮信息的分解問(wèn)題(Matte et al, 2014)。圖3a顯示了13號(hào)站點(diǎn)的T_TIDE和NS_TIDE兩種方法預(yù)測(cè)得到的逐時(shí)潮位數(shù)據(jù)與實(shí)測(cè)數(shù)據(jù)對(duì)比, 圖3b顯示了這兩種方法預(yù)測(cè)得到的潮位時(shí)間序列相對(duì)于實(shí)測(cè)潮位時(shí)間序列的均方根誤差的沿程變化。結(jié)果表明, 對(duì)于三角洲下游站點(diǎn), 兩種方法的計(jì)算誤差相當(dāng), 都在0.16m左右, 但對(duì)于上游站點(diǎn)T_TIDE方法誤差達(dá)到0.3m左右, NS_TIDE方法計(jì)算誤差沿程變化不大, 說(shuō)明NS_TIDE方法能較好反演潮位受徑流影響明顯的站點(diǎn)。因此本文采用非平穩(wěn)調(diào)和分析方法, 以一年數(shù)據(jù)長(zhǎng)度為窗口,計(jì)算13個(gè)水文站點(diǎn)的6個(gè)分潮K1、O1、M2、S2、M4、M6調(diào)和常數(shù)。

圖2 三角洲上游3條通道月均流量圖Fig. 2 Time series of monthly mean river discharges at three upstream channels of the delta

圖3 T_TIDE和NS_TIDE預(yù)測(cè)得到潮位時(shí)間序列(a)和兩種方法預(yù)測(cè)潮位的均方根誤差(b)Fig. 3 Results of T_TIDE and NS_TIDE for the time series of water level hindcasts (a) and root-mean-square errors of water level hindcasts obtained from the two methods (b)

1.3.2 潮汐不對(duì)稱量化方法

1.3 研究方法

1.3.1 潮汐調(diào)和分析方法

三角洲內(nèi)的潮波運(yùn)動(dòng)由于受制于河口特殊地形

式中,A表示振幅比;G表示相位差(單位: °);aM2和aM4分別為M2和M4分潮的振幅(單位: m);gM2和gM4分別為M2和M4分潮的相位(單位: °)。A＞0.01表明潮波發(fā)生較為顯著的變形;G則反映潮汐不對(duì)稱漲潮和落潮之間的優(yōu)勢(shì)關(guān)系, 當(dāng)0＜G＜π時(shí)為“漲潮主導(dǎo)”, 即漲潮歷時(shí)更短; 當(dāng)π＜G＜2π時(shí)為“落潮主導(dǎo)”, 即落潮歷時(shí)更短。

2 結(jié)果

2.1 擋潮閘壩建立前后潮汐分潮變化特征

2.1.1 潮汐分潮空間分布特征

本文使用NS_TIDE調(diào)和分析方法計(jì)算出6個(gè)主要分潮K1、O1、M2、S2、M4、M6振幅和相位, 分析了1970年河口擋潮閘建立前后三角洲內(nèi)部潮汐振幅、相位空間分布特征(圖4, 圖5)。萊茵河-默茲河三角洲河網(wǎng)空間上被分成三部分, 分別為北部通道站點(diǎn)(1, 2, 3, 4, 5, 6)、南部通道站點(diǎn)(10, 11, 12, 13)以及連接南北主通道的中部支流站點(diǎn)(7, 8, 9)。對(duì)比圖4和圖5, 天文分潮M2在擋潮閘建立前后都是最主要的分潮, 其振幅是其他分潮振幅的4～10倍。1970年擋潮閘建立前, M2分潮在進(jìn)入河口后沒(méi)有立即衰減, 站點(diǎn)1到站點(diǎn)8, 站點(diǎn)10到站點(diǎn)11, M2分潮振幅都增加了12cm左右, 且三角洲下游南北通道內(nèi)的潮汐振幅大小相當(dāng), M2分潮振幅都在80cm左右。

圖4 1961—1970年主要分潮(M2、S2、M4、K1、O1、M6)的平均振幅(a、b)和相位(c、d)空間分布圖中數(shù)字為站位號(hào)Fig. 4 Mean tidal amplitudes (a, b) and phases (c, d) for the main constituents (M2, S2, M4, K1, O1, M6) during 1961—1970.The number is the station number

圖5 1971—2018年主要分潮(M2、S2、M4、K1、O1、M6)的平均振幅(a、b)和相位(c、d)空間分布圖中數(shù)字為站位號(hào)Fig. 5 Mean tidal amplitudes (a, b) and phases (c, d) for the main constituents (M2, S2, M4, K1, O1, M6) during 1971—2018.The number is the station number

南部河口擋潮閘建立后, 天文分潮K1、O1、M2、S2進(jìn)入河口后分潮振幅向陸逐漸衰減, 且M2分潮作為主要分潮其沿程衰減幅度大于其他分潮。淺水分潮M4、M6的振幅在北部通道內(nèi)因?yàn)楹涌谑湛s, 在非線性作用下產(chǎn)生并向陸增強(qiáng), 接著在摩擦阻力作用下衰減, 但因?yàn)闇\水分潮振幅較小, 沿程變化不顯著。南北縱向?qū)Ρ? 主要分潮振幅由北向南減小,且分潮振幅沿程衰減幅度也由北向南減小。以M2分潮為例, 北部通道由80cm衰減至49cm, 南部通道由13cm衰減至10cm。擋潮閘建立前后, 分潮相位始終由海向陸沿程逐漸增加, O1分潮上下游相位差最小, 為43°, M6分潮上下游相位差最大, 為232°,M6分潮相位沿程增加最快, 符合頻率越高的分潮在河道內(nèi)衰減得越快的規(guī)律(Godin, 1999)。

2.1.2 潮汐分潮時(shí)間變化特征利用實(shí)測(cè)潮位1961—2018年的6個(gè)站點(diǎn)主要分潮振幅差(Δa=a1971-2018-a1961-1970,其中a為分潮振幅,a1961-1970、a1971-2018分別表示擋潮閘建設(shè)前后的分潮平均振幅)和相位差(Δg=g1971-2018-g1961-1970, 其中g(shù)為分潮相位,g1961-1970、g1971-2018分別表示擋潮閘建設(shè)前后的分潮平均相位)衡量河口擋潮閘對(duì)三角洲內(nèi)分潮振幅和相位的影響。如圖6所示, 1970年河口擋潮閘壩建設(shè)對(duì)中部、南部站點(diǎn)潮動(dòng)力影響最為劇烈, 6種主要分潮振幅明顯減弱, 其中M2分潮作為主要分潮, 其減小幅度大于其他分潮。擋潮閘建設(shè)后, 南部站點(diǎn)M2分潮振幅普遍減小了80%左右, 中部8號(hào)站點(diǎn)M2分潮振幅減小了60%。淺水分潮M4由M2分潮非線性作用產(chǎn)生的, 所以M2分潮振幅的明顯減小也會(huì)導(dǎo)致M4分潮產(chǎn)生相應(yīng)的變化, 在南部和中部通道分別減小了70%和20%左右。因?yàn)?970年哈靈河口擋潮閘壩的建立, 一方面減弱了北海潮動(dòng)力對(duì)南部通道的直接影響, 潮汐僅通過(guò)北部新水道進(jìn)入三角洲流域, 另一方面河網(wǎng)流量的再分配, 北部新水道以及中部支流排放的流量比例大幅度增加, 徑流量增大會(huì)加大天文分潮能量損失, 所以中部站點(diǎn)分潮振幅明顯減小, 但北部站點(diǎn)仍受潮動(dòng)力主導(dǎo), 分潮振幅變化不大。此外圖6b顯示在1970年工程建設(shè)后, 南部站點(diǎn)主要天文分潮相位明顯增大, 其他站點(diǎn)略微減小,反映潮波在南部通道內(nèi)傳播速度下降。

圖6 1970年前后主要分潮的振幅差(a)和相位差(b)1和4是北部站點(diǎn), 8是中部站點(diǎn), 10、11和13是南部站點(diǎn)Fig. 6 Differences in tidal amplitudes (a) and phases (b) for the main constituents before and after 1970. Northern stations (1,4), central station (8), southern stations (10, 11, 13)

2.2 擋潮閘壩建立前后潮汐不對(duì)稱變化特征

2.2.1 潮汐不對(duì)稱空間分布特征

基于分潮振幅比A(M4/M2)和相位差G(2M2-M4)研究三角洲的潮汐不對(duì)稱特性。圖7a為三角洲內(nèi)各站點(diǎn)振幅比空間分布, 潮汐不對(duì)稱在三角洲北部通道內(nèi)向上游沿程先增強(qiáng), 并在5號(hào)站點(diǎn)附近潮波變形最大, 接著向上游略微減弱, 這與M2、M4分潮在上游地區(qū)受到較大徑流作用同時(shí)發(fā)生衰減有關(guān)。潮汐不對(duì)稱在中部和南部總體向上游增強(qiáng),A分別由0.24增大至0.30, 0.11增大至0.13, 除了南部口門(mén)處站點(diǎn)位于擋潮閘附近, 可能受泄水閘開(kāi)閉影響, 潮汐不對(duì)稱較上游站點(diǎn)偏大。中部通道潮汐不對(duì)稱現(xiàn)象最為顯著, 南部通道潮汐不對(duì)稱現(xiàn)象較弱, 因?yàn)楹涌趽醭遍l建立后南部通道潮動(dòng)力較三角洲其他地區(qū)小, 而連接南北主通道的中部支流中通過(guò)的徑流量較大, 對(duì)潮動(dòng)力的非線性作用增大, 潮汐不對(duì)稱現(xiàn)象更加顯著。圖7b為三角洲內(nèi)各站點(diǎn)相位差參數(shù)G, 均在0～180°范圍內(nèi), 表明漲潮時(shí)間要短于落潮時(shí)間, 所以萊茵河-默茲河三角洲是漲潮主導(dǎo)型潮汐不對(duì)稱。Speer等(1985)研究表明, 對(duì)于不同的淺水分潮, 其與產(chǎn)生它的兩個(gè)半日分潮之間的相對(duì)相位基本相同, 說(shuō)明用G判斷三角洲潮汐不對(duì)稱的方向是可靠的。并且研究發(fā)現(xiàn)相位差參數(shù)G向陸沿程增加,站點(diǎn)1到站點(diǎn)6從9°增大到56°, 站點(diǎn)7到站點(diǎn)9從23°增大到77°, 站點(diǎn)10到站點(diǎn)13從32°增大到137°,相位差在三角洲內(nèi)向上游沿程逐漸增加也可能是潮汐不對(duì)稱向上游增強(qiáng)的一個(gè)重要原因。

圖7 1971—2018年M2、M4分潮平均振幅比(a)和相位差(b)空間分布圖中數(shù)字為站位號(hào)Fig. 7 Mean tidal amplitude ratio (a) and relative phase (b) of the constituents M2 and M4 during 1971—2018. The number is the station number

2.2.2 潮汐不對(duì)稱時(shí)間變化特征

根據(jù)6個(gè)站點(diǎn)振幅比A(M4/M2)在1961—2018年期間年平均變化(圖8), 分析得到1970年擋潮閘建成前后三角洲內(nèi)的潮汐不對(duì)稱現(xiàn)象的變化情況。1970年哈靈水道擋潮閘壩建立沒(méi)有改變?nèi)侵逎q潮主導(dǎo)的潮汐不對(duì)稱特征, 但在三角洲河網(wǎng)不同區(qū)域都發(fā)生了不同程度不同性質(zhì)的變化。北部上游4號(hào)和中部8號(hào)站點(diǎn)的振幅比在1970年明顯增大, 4號(hào)站點(diǎn)由0.14增大至0.22, 8號(hào)站點(diǎn)由0.15增大至0.30,中部站點(diǎn)潮汐不對(duì)稱現(xiàn)象在擋潮閘建成后顯著增強(qiáng)。南部站點(diǎn)的潮汐不對(duì)稱變化較復(fù)雜, 下游10號(hào)、11號(hào)站點(diǎn)潮汐不對(duì)稱明顯增強(qiáng), 雖然在后期由于河道淤積的影響, 潮汐不對(duì)稱有個(gè)緩慢減弱并趨于穩(wěn)定的過(guò)程; 上游13號(hào)站點(diǎn)分潮振幅比由0.16減小為0.13, 潮汐不對(duì)稱因?yàn)镸4分潮的顯著減小而減弱。北部入?？?號(hào)站點(diǎn)徑潮動(dòng)力環(huán)境受南部擋潮閘影響較小, 潮汐不對(duì)稱變化不顯著。因?yàn)檎军c(diǎn)分布有限, 對(duì)于三角洲其他區(qū)域潮汐不對(duì)稱特性是否發(fā)生不同變化無(wú)法進(jìn)一步分析, 已有站點(diǎn)潮汐不對(duì)稱的變化情況只能反映局部變化情況。

圖8 1961—2018年期間6個(gè)水文站點(diǎn)年平均振幅比A (M4/M2)逐年變化圖中藍(lán)色站點(diǎn)表示北部站點(diǎn)(a、b), 紅色表示中部站點(diǎn)(c), 綠色表示南部站點(diǎn)(d、e、f)Fig. 8 The mean tidal amplitude ratio (M4/M2) of six hydrological stations during 1961—2018

3 分析與討論

3.1 分潮振幅變化對(duì)潮汐不對(duì)稱的影響

根據(jù)公式(1), 衡量潮汐不對(duì)稱的指標(biāo)振幅比A隨天文潮M2振幅沿程的減小和淺水分潮M4振幅的增大而增大, 這解釋了外海潮波進(jìn)入河口后, 在地形摩擦阻力和徑流動(dòng)力等非線性作用下潮波變形程度隨著M2分潮振幅的向陸衰減和M4分潮的向陸增加增強(qiáng)。在上游地區(qū), 徑流量作用增大, 摩擦阻力進(jìn)一步增加, M2和M4分潮共同衰減導(dǎo)致潮汐不對(duì)稱減弱(童朝鋒 等, 2020)。這一規(guī)律在加拿大圣勞倫斯河口(Godin, 1999)、荷蘭西斯凱爾特河口(Wang et al, 2002)、中國(guó)長(zhǎng)江三角洲(Yu et al, 2020)等潮汐河口都得到驗(yàn)證。而擋潮閘建成導(dǎo)致潮汐振幅都發(fā)生劇烈變化, M2分潮在中部、南部下游站點(diǎn)的減小幅度大于M4分潮的衰減幅度, 最終相應(yīng)位置處潮汐不對(duì)稱增強(qiáng)。

3.2 徑流量和潮動(dòng)力變化對(duì)潮汐不對(duì)稱的影響

擋潮閘建成導(dǎo)致萊茵河-默茲河三角洲內(nèi)的徑潮動(dòng)力環(huán)境產(chǎn)生了短時(shí)間強(qiáng)烈的響應(yīng)。為了研究三角洲內(nèi)擋潮閘建設(shè)前的各個(gè)站點(diǎn)對(duì)于上游徑流量變化和外海潮動(dòng)力變化的響應(yīng), 借助NS_TIDE方法,通過(guò)將輸入的上游徑流量和外海潮差時(shí)間序列分別乘以一定倍數(shù)(0.8,0.9,1.0,1.1,1.2)來(lái)評(píng)估三角洲內(nèi)的潮波對(duì)變化的徑流和潮差的響應(yīng)。總體上, 上游徑流和外海潮差增大后, 三角洲內(nèi)各站點(diǎn)平均水位是逐漸上升的。在變化的上游徑流影響下, 主要分潮M2、M4振幅、相位和潮汐不對(duì)稱指標(biāo)變化情況如圖9所示。上游徑流的增加, 削弱了下游站點(diǎn)M2、M4分潮振幅, 增加了分潮相位, 但對(duì)于上游地區(qū)的8號(hào)和13號(hào)站點(diǎn), M4分潮振幅隨著徑流量增加而增大。三角洲大部分站點(diǎn)潮汐不對(duì)稱隨著徑流量的增大逐漸增強(qiáng), 且上游地區(qū)變化更為顯著, 說(shuō)明徑流量增大促進(jìn)了M2分潮能量向M4分潮轉(zhuǎn)移(圖9c)。擋潮閘建成之后, 如果北部和中部站點(diǎn)處只考慮徑流量增加, 其潮汐不對(duì)稱也相應(yīng)增強(qiáng), 且中部上游站點(diǎn)變化顯著。

圖9 M2 、M4分潮振幅、相位和兩者振幅比(A)、相位差(G)對(duì)不同量級(jí)徑流量的敏感性a. M2分潮振幅; b. M4分潮振幅; c. 振幅比(M4/M2); d. M2分潮相位; e. M4分潮相位; f. 相位差(2M2-M4)。Q為上游徑流量Fig. 9 Sensitivity of tidal amplitudes and phases of the constituents M2 and M4 and tidal amplitude ratio and relative phase of the constituents M2 and M4 to variations in river discharge. (a) M2 tidal amplitude; (b) M4 tidal amplitude; (c) tidal amplitude ratio (M4/M2); (d) M2 tidal phase; (e) M4 tidal phase; (f) relative phase (2M2-M4). Q is river discharge

在變化的外海潮差情況下, 主要分潮M2、M4振幅、相位和潮汐不對(duì)稱指標(biāo)變化情況如圖10所示。外海潮差的緩慢增加會(huì)導(dǎo)致M2分潮振幅減小,相位增大, 從而導(dǎo)致M2分潮的衰減。在三角洲下游地區(qū), 隨著潮差的增加, M4分潮振幅減小、相位增加; 在上游地區(qū), 隨著潮差增加, M4分潮相位卻減小, 這可能與小潮期間較強(qiáng)的底摩阻使其衰減有關(guān)(Matte et al, 2014)。隨著潮動(dòng)力增強(qiáng), 在三角洲下游地區(qū), 潮汐不對(duì)稱減弱, 而上游地區(qū)潮汐不對(duì)稱則增強(qiáng)(圖10c)。實(shí)際上擋潮閘建成后三角洲河網(wǎng)各個(gè)區(qū)域同時(shí)受到徑流量和外海潮差變化的影響, 而且其影響是短時(shí)間且變化劇烈。綜合兩個(gè)因素來(lái)看,南部站點(diǎn)外海潮差的突然減小是其潮汐不對(duì)稱變化的主導(dǎo)因素, 而北部和中部站點(diǎn)潮汐不對(duì)稱的變化則主要受徑流量的變化主導(dǎo)。

圖10 M2 、M4分潮振幅、相位和兩者振幅比(A)、相位差(G)對(duì)不同量級(jí)潮差的敏感性a. M2分潮振幅; b. M4分潮振幅; c. 振幅比(M4/M2); d. M2分潮相位; e. M4分潮相位; f. 相位差(2M2-M4)。R為外海潮差Fig. 10 Sensitivity of tidal amplitudes and phases of the constituents M2 and M4 and tidal amplitude ratio and relative phase of the constituents M2 and M4 to variation in tidal range. (a) M2 tidal amplitude; (b) M4 tidal amplitude; (c) tidal amplitude ratio (M4/M2); (d) M2 tidal phase; (e) M4 tidal phase; (f) relative phase (2M2-M4). R is tidal range

4 結(jié)論

本文以河口擋潮閘修建影響下的荷蘭萊茵河-默茲河三角洲為研究對(duì)象, 通過(guò)分析三角洲內(nèi)多個(gè)站點(diǎn)1960—2018年逐時(shí)潮位數(shù)據(jù), 研究了在河口擋潮閘建立后三角洲內(nèi)的潮汐動(dòng)力和潮汐不對(duì)稱特性的時(shí)空變化規(guī)律, 得出以下結(jié)論:

1) 荷蘭萊茵河-默茲河三角洲主要分潮振幅、相位時(shí)空變化顯著。1970年以前, 潮波進(jìn)入河口后,南北通道內(nèi)振幅大小相當(dāng), 且主要天文分潮M2振幅由于河岸收縮匯聚效應(yīng)先略微增大后衰減。南部河口擋潮閘建設(shè)后, 北部和中部通道徑流量增大,摩擦阻力增加, 主要天文分潮振幅由河口向上游沿程衰減。時(shí)間上, 南部擋潮閘建設(shè)導(dǎo)致三角洲中部、南部站點(diǎn)分潮振幅明顯減小。

2) 荷蘭萊茵河-默茲河三角洲內(nèi)的潮波為漲潮主導(dǎo)型不對(duì)稱, 即漲潮時(shí)間短于落潮時(shí)間。三角洲是正規(guī)半日潮, 主要天文分潮M2與其對(duì)應(yīng)的淺水分潮M4之間相互作用是三角洲潮汐不對(duì)稱的主要來(lái)源, 潮波變形程度總體向陸逐漸增強(qiáng), 到達(dá)某一點(diǎn)后由于受到徑流作用較大, 潮汐能量明顯削弱, 潮汐不對(duì)稱反而減弱。

3) 短時(shí)間高強(qiáng)度的擋潮閘工程對(duì)萊茵河-默茲河三角洲潮動(dòng)力和徑流的傳播產(chǎn)生了較大影響。北部、中部站點(diǎn)潮汐不對(duì)稱在南部擋潮閘建立后增強(qiáng),中部站點(diǎn)主要因?yàn)橥ǖ缽搅髁吭黾? 潮汐能量由M2分潮向M4分潮轉(zhuǎn)變, 潮汐不對(duì)稱性顯著增強(qiáng)。而南部通道主要受潮動(dòng)力變化影響, 下游站點(diǎn)因?yàn)镸2分潮振幅的劇烈下降導(dǎo)致潮汐不對(duì)稱增強(qiáng), 上游站點(diǎn)潮汐不對(duì)稱減弱則是因?yàn)镸4分潮振幅的衰減程度超過(guò)了M2分潮振幅的衰減程度。