伶仃洋河口橫向環(huán)流

易侃，龔文平*

（1.中山大學(xué)海洋學(xué)院近岸海洋科技研究中心，廣東廣州 510275）

伶仃洋河口橫向環(huán)流

易侃1，龔文平1*

（1.中山大學(xué)海洋學(xué)院近岸海洋科技研究中心，廣東廣州 510275）

橫向環(huán)流存在于大多數(shù)河口中，在河口的動量平衡和物質(zhì)輸運過程中起著重要的作用，其形成的機(jī)制眾多。本文利用三維水動力模型EFDC（Environmental Fluid Dynamics Code），結(jié)合1989年和2002年的實測數(shù)據(jù)對伶仃洋河口進(jìn)行研究，得到了伶仃洋河口在洪季和枯季時期不同的水動力和鹽度分布結(jié)構(gòu)。在文中選取了具有代表性的兩個斷面，針對洪枯季、大小潮以及漲落潮等不同時期的橫向流結(jié)構(gòu)和鹽度剖面分布特點進(jìn)行比較和分析。結(jié)果顯示，受局地的徑、潮流條件和地形約束影響，橫向環(huán)流結(jié)構(gòu)并不能充分形成，在不同的區(qū)域表現(xiàn)出截然不同的結(jié)構(gòu)特點。伶仃洋河口橫向環(huán)流主要是在徑、潮流與地形三者之間的相互作用下形成的，不同的區(qū)域的形成機(jī)制存在差異。綜合分析后，本文認(rèn)為伶仃洋河口的橫向環(huán)流在該區(qū)域的動力過程中相對較為次要。

伶仃洋河口；EFDC模型；河口橫向環(huán)流；水動力

1 引言

河口橫向環(huán)流作為河口地區(qū)特有的動力過程，能夠影響河口縱向擴(kuò)散速率［1—3］，引起強(qiáng)烈的垂向混合［4］，造成泥沙的橫向輸運［5］，以及調(diào)整動量收支平衡［6］，對于河口地區(qū)的動量和物質(zhì)輸運具有十分重要的意義。已有的研究表明，橫向環(huán)流的形成機(jī)理大致可以歸納為以下幾種：（1）彎道效應(yīng)［4—6］；（2）地球偏轉(zhuǎn)的科氏力效應(yīng)［7］；（3）在高度分層且具有陡峭橫剖面的河口的邊界混合［7—8］；（4）縱向密度梯度的差異平流輸運［6］；（5）風(fēng)驅(qū)動作用［9］；（6）潮汐與河道橫剖面水深變化的相互作用［10］。在實際的河口中，橫向環(huán)流的形成常常是各種機(jī)理共同作用的結(jié)果，而在不同的河口系統(tǒng)中由不同的機(jī)理產(chǎn)生的橫向環(huán)流的動力特征以及各機(jī)理相對重要性都存在很大的差異。相較于國外的研究成果而言，國內(nèi)針對河口橫向環(huán)流的研究十分有限。

珠江口伶仃洋為珠江三角洲內(nèi)最大的河口灣，其動力條件與地形均十分復(fù)雜。在不同的季節(jié)（洪季與枯季）與不同的潮汐條件（大潮與小潮）下，河口的分層與混合發(fā)生著變化，同時河口的縱向環(huán)流也發(fā)生相應(yīng)的變化［11］。但橫向環(huán)流的變化如何，目前還鮮有研究。本文即針對這一研究空白，采用EFDC模型進(jìn)行數(shù)值模擬研究。在對模型進(jìn)行較充分的驗證后，選取伶仃洋的不同斷面，研究不同季節(jié)、不同潮汐條件下斷面的橫向環(huán)流的變化，探討橫向環(huán)流的形成機(jī)制。

珠江口伶仃洋是位于珠江三角洲東南側(cè)，呈NNW－SSE向的喇叭形河口灣，也是華南最大的河口灣。伶仃洋地勢由灣頂向灣口傾斜，水深由西向東遞增，整個伶仃洋河口形似漏斗。灣口外島嶼棋布。東江、北江以及西江帶來的大量的淡水和泥沙由虎門、蕉門、洪奇瀝和橫門4個口門匯入伶仃洋，其中除虎門為潮流優(yōu)勢型河口外，其他3個口門都為徑流優(yōu)勢型。伶仃洋水下地形表現(xiàn)出“三灘兩槽”的基本格局，由西向東依次為西灘、西槽、中灘、東槽以及東灘（圖1）。

圖1 伶仃洋河口平面示意圖Fig.1 The map of the Lingding Estuary

伶仃洋河口灣受潮流和徑流的相互作用，徑流年際變化不大，從北到南4個口門的年平均徑流量為1 910、1 790、663和1 158 m3／s［12］，總徑流量約1 741.1×108m3／a。洪、枯季流量變化較大，洪季的徑流量約占年徑流量的80%以上［13］。伶仃洋潮汐類型屬于不正規(guī)半日潮類型，平均潮差在1.1 m左右（橫門），潮差較小，潮差沿程逐漸增大，在灣頂最大可達(dá)3.7 m。由于伶仃洋的平均深度較淺，橫斷面較寬（7～40 km），且大于內(nèi)羅斯貝半徑（小于15 km），科氏力較為重要。在科氏力的作用下漲落潮流路分離，漲潮流主要從內(nèi)伶仃島東面的深槽入侵，加上珠江淡水入海的4個口門位于伶仃洋的西北，故東南部的鹽度比西北部高，東西橫向鹽度差最大可達(dá)20［14］。珠江徑流從西側(cè)注入和水流橫向密度差異等因素使得伶仃洋存在著由西向東傾斜的橫向水位比降，其水位差與潮汐漲落、潮流大小有關(guān)，最大值出現(xiàn)在潮流落急的前后［15］。

2 研究方法

2.1 資料

本文利用交通部廣州海上監(jiān)督局出版的1989年海圖以及中華人民共和國海事局出版的1998年海圖進(jìn)行數(shù)字化作為地形條件分別用于枯季和洪季的研究。在伶仃洋區(qū)域進(jìn)行了大量的現(xiàn)場觀測，由于1989年枯季與2002年洪季時觀測站位較多，且資料較為全面，在這里作為數(shù)據(jù)源，對模型加以驗證，測站位置見圖1。

2.2 模型及計算條件

本文利用EFDC模型模擬伶仃洋河口的水動力和鹽度輸運過程。EFDC模型是由美國弗吉利亞海洋研究所（VIMS）的Hamrick開發(fā)［16－17］。模型中水平方向采用正交曲線網(wǎng)格，垂向上采用sigma坐標(biāo)。對時間積分采用二階精度的三步蛙跳格式，并在一定時間間隔插以二步的歐拉格式，對平流項采取二階的迎風(fēng)格式。擴(kuò)散系數(shù)采用Mellor Yamada 2.5階的湍閉合模型進(jìn)行計算。

本文主要研究伶仃洋河口的水動力過程，建立EFDC模型，為了使它能夠用于伶仃洋河口的多項問題的研究，并且充分考慮鹽水入侵對河流上游的影響，因此計算區(qū)域的徑流邊界在伶仃洋河口的四大口門分別向上游延伸5～10 km不等，各支流岔口向上游延伸1～5 km不等。在水平方向上共81 081個網(wǎng)格單元，空間分辨率在橫向上平均為270 m，最小為20 m，縱向上平均為360 m，最小為16 m；垂向網(wǎng)格采用σ坐標(biāo)分為15層；時間步長為10 s。模型網(wǎng)格如圖2所示，由于伶仃洋內(nèi)地形變化較大，我們分別采用1989年與1998年的地形進(jìn)行枯季與洪季的計算。邊界條件的給定（包括水位和鹽度）采用模型嵌套的方法，將包括整個珠三角的大模型的計算結(jié)果內(nèi)插到伶仃洋的開邊界處，大模型在外海邊界的網(wǎng)格分辨率為10 km。在洪、枯季，上游的開邊界條件利用西江的梧州、北江的石角、東江的博羅的實時流量序列，采用大模型計算所得。在本次模擬中，沒有考慮風(fēng)的影響?？菁镜挠嬎銜r間為1988年2月1日至3月31日，洪季的計算時間為2002年6月1日至7月31日。

圖2 研究區(qū)域枯季和洪季網(wǎng)格圖Fig.2 Model grid of the study area for the dry and flood season

3 模型驗證

采用實際觀測數(shù)據(jù)對模型進(jìn)行驗證，枯季時為1988年3月初在伶仃洋進(jìn)行的定點觀測數(shù)據(jù)，共18個測站；洪季時為2002年6月在伶仃洋東西槽進(jìn)行的定點觀測，共12個站，觀測分別在大小潮期進(jìn)行，各站點的地理位置如圖1所示。本文中挑選具有代表性的站點將模型計算結(jié)果與實測值進(jìn)行比較，如圖3、4，圖中點線表示實測值，實線表示模擬值。

圖3 枯季大潮期間比對情況Fig.3 Model calibration during the spring tide in the dry season

圖4 洪季大潮期間對比情況Fig.4 Model calibration during the spring tide in the flood season

從對比的結(jié)果來看，枯季時模擬的水位、流速流向及鹽度與實測情況基本吻合。洪季時大潮期模擬的水位、流速流向與實測值比對結(jié)果良好，模擬的鹽度的量值和變化趨勢與實測情況較為一致。模型計算中還存在一定誤差，如枯季時流速的模擬值與實測值相比，還存在振幅與相位上的差別。由于伶仃洋的地形與邊界條件異常復(fù)雜，總體而言，本次模型模擬的結(jié)果與實測情況吻合良好，各變量在量值上和變化趨勢上都與實際情況基本一致，可以認(rèn)為，模擬的結(jié)果具有較高的可信度。

4 結(jié)果與分析

4.1 伶仃洋水動力和鹽度分布特征

為更加直觀地了解伶仃洋的水動力結(jié)構(gòu)及鹽分輸運的特點，本文中選取了幾個有代表性的時刻，描述其垂向平均的流速和鹽度在整個伶仃洋河口的分布（圖5，圖6）。

模擬的結(jié)果顯示，伶仃洋河口的鹽度分布存在橫向差異，由于4個大口門的位置偏西，再加上科氏力的作用，上游來的徑流主要從西邊流入，由外海輸入的高鹽水在科氏力作用下向東偏移，形成鹽度東高西低的分布特征。水動力在橫向上存在差異，深槽的流速大于淺灘，高鹽水主要沿深槽上溯，深槽的鹽度大于淺灘。等鹽度線基本上成東北—西南向排列，并且在東部沿著深槽處具有雙峰結(jié)構(gòu)。在枯季（圖5），鹽水入侵的距離較遠(yuǎn)，在4個口門處，鹽度在15～20之間，其中，虎門的鹽水入侵情況最為顯著，鹽度可達(dá)25。枯季大小潮期間的鹽度分布也存在差異，在小潮期間，由于鹽度的分層加強(qiáng)，重力環(huán)流更為顯著，有利于底部的高鹽水入侵（未用圖顯示）。洪季時（見圖6），大量的徑流從4個口門注入伶仃洋，在河口的上段，內(nèi)伶仃島以北區(qū)域主要受淡水控制，鹽度整體偏低。在河口的下段，鹽度的分布同樣表現(xiàn)出東高西低，深槽高淺灘低的特征。

圖5 枯季大潮時期漲急（a）和落急（b）垂向平均的流速和鹽度分布圖Fig.5 The distribution of the vertically averaged current and salinity for the maximum flood tide（a）and the maximum ebb tide（b）during the spring tide in the dry season

伶仃洋地區(qū)的水動力受到徑流和海洋相互作用的影響，同時受到地形的約束，因此十分復(fù)雜。漲潮時，由外海傳入的潮流在進(jìn)入伶仃洋后，由于地形的展寬以及科氏力的作用，在內(nèi)伶仃島處分成東西兩支，西分支沿著伶仃水道向北流動，東分支與從香港暗士頓水道流入的潮流匯合后沿礬石水道上溯，最后兩股水流匯于川鼻水道，流入虎門內(nèi)。落潮時的流路相反。

4.2 伶仃洋橫向水流結(jié)構(gòu)

圖6 洪季大潮時期漲急（a）和落急（b）垂向平均流速和鹽度分布圖Fig.6 The distribution of the vertically averaged current and salinity for the maximum flood tide（a）and the maximum ebb tide（b）during the spring tide in the flood season

為了研究伶仃洋的橫向水流結(jié)構(gòu)特征，探討橫向環(huán)流的形成特點，在伶仃洋河口中選取了兩個具有代表性的橫剖面進(jìn)行研究，剖面的位置如圖1所示，C1表示剖面1，C2表示剖面2。橫剖面1位于虎門口門附近，該地區(qū)受虎門流入的徑流和外海傳入的潮流兩方面的影響較為顯著，且只有單個深槽。橫剖面2位于伶仃洋河口的中段，包括伶仃水道和礬石水道兩個深槽，并且橫向鹽度差異明顯。由于縱向水流主要沿深槽流動，而伶仃洋中兩條深槽向北不斷輻聚，成“八”字形結(jié)構(gòu)，單一的橫剖面無法同時與兩條水槽垂直，由此得到的橫向流速分量實際上包含縱向流速投影到橫剖面所產(chǎn)生的誤差。為了減少這種誤差，剖面2并不是一條直線，而是由分別垂直于兩個深槽的折線組成。根據(jù)選取的橫剖面，得到不同時期的橫向水流結(jié)構(gòu)和鹽度分布的結(jié)構(gòu)（見圖7～10，圖中垂向流速放大100倍）。

在枯季（見圖7），由于鹽水入侵的距離較遠(yuǎn)，剖面1的鹽度較高，并且在漲落潮都表現(xiàn)出深槽鹽度高，淺灘鹽度低，底層鹽度高，表層鹽度低的特征。在小潮漲急時（見圖7a），鹽度垂向混合較好，在深槽的上部存在一個逆時針的橫向環(huán)流，其他區(qū)域的橫向流基本表現(xiàn)為從左右兩邊指向深槽，并且表底層的流速流向一致；在小潮落急時（見圖7b），鹽度在深槽處出現(xiàn)分層，在深槽的左上方存在一個順時針的橫向環(huán)流，在深槽的右邊，橫向水流指向右側(cè)淺灘，表底層流向一致。枯季大潮漲急時期（見圖7c），鹽度的垂向混合增強(qiáng)，在深槽的右上方存在一個小的順時針的環(huán)流圈，其他區(qū)域橫向流速在表底層都由淺灘指向深槽；在大潮落急時（見圖7d），深槽的左上方存在順時針的橫向環(huán)流結(jié)構(gòu)，其他區(qū)域流向與漲潮時相反。在洪季小潮時期（見圖8a、8b），整個剖面1的鹽度基本為零，只在右邊界附近鹽度大于0.5，無論漲落潮時期，橫向流速在表底層基本一致，漲落潮流向相反，沒有發(fā)現(xiàn)橫向環(huán)流。在洪季大潮期（見圖8c、8d），鹽水入侵有所增強(qiáng)，在剖面1存在橫向的鹽度差異，在漲急時期在深槽的左上方存在逆時針的環(huán)流，而在落急時期，環(huán)流結(jié)構(gòu)消失，整個剖面的橫向水流流向基本一致，表現(xiàn)為從深槽指向淺灘。

橫剖面2的橫向動力結(jié)構(gòu)與剖面1不同，在枯季（見圖9），鹽度的垂向混合強(qiáng)烈，表底層鹽度基本一致，表底層橫向流速流向也基本一致，在漲急時都指向西，落急時都指向東，在東槽左側(cè)的斜坡上一段小區(qū)域出現(xiàn)表底層相反的水流。洪季小潮時期（見圖10a、10b），高鹽度出現(xiàn)在深槽的底部，在淺灘區(qū)域鹽度混合充分。洪季小潮漲急時，在東槽存在一個順時針的橫向環(huán)流，在西槽的底層，圍繞著高鹽度區(qū)域也存在一個小的順時針環(huán)流；在落急時，環(huán)流結(jié)構(gòu)消失，表底層流向一致。在洪季大潮時期（見如圖10c、10d），鹽度的垂向混合增強(qiáng)，在漲急時，在深槽處同樣出現(xiàn)表底層流向相反的橫向水流結(jié)構(gòu)，而在落急時，表底層流向又重新恢復(fù)一致。

圖7 橫剖面1枯季小潮潮期間漲急（a）、落急（b）以及大潮期間漲急（c）和落急（d）橫向流速（矢量）和鹽度（等值線）分布圖Fig.7 The distribution of current（vector arrows）and salinity（contours）for the maximum flood（a），the maximum ebb（b）during the neap tide and the maximum flood（c），the maximum ebb（d）during the spring tide in the dry season in cross section 1

參照Lerczak和Rockwell［7］的研究方法，根據(jù)公式（1）將橫向流速的大小沿橫剖面平均得到的橫向流速振幅用來度量橫向流的大?。?/p>

式中，｜v｜是橫向流速v的絕對值，A是橫剖面的面積。從計算得到的不同時期的橫向流速振幅在一個潮周期內(nèi)的變化過程中（見圖11）可以看到，在枯季，剖面1的橫向水流的振幅在一個潮周期內(nèi)的變化沒有明顯的規(guī)律，并且大小潮時都在0.06 m／s左右，最大值達(dá)0.08 m／s。在洪季小潮期，剖面1的橫向水流的振幅變化呈周期性變化，在落潮期間出現(xiàn)最大值，可達(dá)0.075 m／s，最小值在漲潮期間（0.02 m／s）；在洪季大潮期，這種周期性的變化特征明顯減弱，橫向水流的振幅較小潮期更大，最大可達(dá)0.08 m／s，最小為0.03 m／s左右。剖面2的橫向流振幅相比于剖面1表現(xiàn)出完全不同的變化特征，無論是洪季大小潮還是枯季大小潮都具有明顯的周期性?？菁敬蟪睍r期的最大振幅為0.08 m／s，大于小潮時期（0.07 m／s）；洪季的情況則相反，小潮時期的振幅大于大潮時期，落潮時期的振幅大于漲潮期，然而變化的相位基本一致。

圖8 橫剖面1洪季小潮期間漲急（a）、落急（b）以及大潮期間漲急（c）和落急（d）橫向流速（矢量）和鹽度（等值線）分布圖Fig.8 The distribution of current（vector arrows）and salinity（contours）for the maximum flood（a），the maximum ebb（b）during the neap tide and the maximum flood（c），the maximum ebb（d）during the spring tide in the flood season in cross section 1

從剖面1和剖面2的模擬結(jié)果可以看到，由于伶仃洋復(fù)雜的地形和水動力條件的影響，模擬得到的橫向環(huán)流的形成并不充分，只在剖面的局部區(qū)域出現(xiàn)環(huán)流結(jié)構(gòu)。在剖面的大部分區(qū)域，橫向水流在表底層流向一致，并且隨漲落潮變化。顯然，這些表底層流向一致的橫向流其實是縱向流沿剖面投影的分量，主要是由于潮流受地形的影響，流路發(fā)生變化。在剖面1，漲潮流隨著深槽向上游收縮變窄而向中間匯聚，流速投影在剖面1上，則產(chǎn)生表底層都指向深槽中間的橫向分量，落潮時結(jié)果相反。剖面2處，由于深槽的地形變化，潮流也會在剖面上投影出表底層流向一致的流速分量。橫向流速的振幅隨漲落潮變化的特征與縱向流速的變化特征相似，也可反映出其中的相關(guān)聯(lián)性。接下來，我們將通過對不同剖面中具有明顯橫向環(huán)流結(jié)構(gòu)的時段進(jìn)行研究分析，探究其形成的機(jī)制。

圖9 橫剖面2枯季小潮期間漲急（a）和落急（b）、大潮期間漲急（c）和落急（d）橫向流速（矢量）和鹽度（等值線）分布圖Fig.9 The distribution of current（vector arrows）and salinity（contours）for the maximum flood（a），the maximum ebb（b）during the neap tide and the maximum flood（c），the maximum ebb（d）during the spring tide in the dry season in cross section 2

圖10 橫剖面2洪季小潮期間漲急（a）和落急（b）、大潮期間漲急（c）和落急（d）橫向流速（矢量）和鹽度（等值線）分布圖Fig.10 The distribution of current（vector arrows）and salinity（contours）for the maximum flood（a），the maximum ebb（b）during the neap tide and the maximum flood（c），the maximum ebb（d）during the spring tide in the flood season in cross section 2

圖11 橫向流速振幅在潮周期內(nèi)從漲憩開始隨時間變化過程圖Fig.11 The intertidal amplitude variation of the transverse velocity beginning from the high slack water

5 討論

如前所述，河口橫向環(huán)流主要的形成機(jī)制包括：彎道效應(yīng)，地轉(zhuǎn)偏向力，邊界層混合，差異對流，風(fēng)以及潮流與橫剖面地形相互作用等。其中，彎道效應(yīng)，地轉(zhuǎn)偏向力以及差異對流在充分混合的水流中能夠驅(qū)動橫向環(huán)流［7］，而邊界層混合只有在鹽度分層較強(qiáng)的系統(tǒng)中才顯得重要［8］。本文研究的伶仃洋河口成喇叭形，且在模擬的過程中沒有考慮風(fēng)的因素，可排除彎道效應(yīng)和風(fēng)的影響。模擬的結(jié)果顯示，在枯季鹽度的垂向混合較充分，在洪季則為部分混合，不符合邊界混合機(jī)制形成的條件，且得到的環(huán)流結(jié)構(gòu)與該機(jī)制的結(jié)果不符，也可排除。

模擬結(jié)果中，剖面1的鹽度結(jié)構(gòu)在枯季的漲落潮都表現(xiàn)為深槽更大，并且在漲潮期間在深槽處形成的橫向環(huán)流的流向與差異對流產(chǎn)生的斜壓梯度力和地轉(zhuǎn)偏向力驅(qū)動的橫向流方向相反，由此可以判斷斜壓梯度力和地轉(zhuǎn)偏向力也不是其最主要的驅(qū)動力。已有的研究表明［7］，橫向動量平衡方程中主要是橫向壓力梯度與垂向剪切擴(kuò)散之間的平衡，因此在垂向混合較強(qiáng)的河口，由差異對流產(chǎn)生的橫向流的大小滿足：

式中，g表示重力加速度，H表示河口的平均深度，Av表示垂向渦動黏性系數(shù)，B表示河口寬度，Δρ（y）表示由差異對流引起的橫向密度差異，ρo表示平均密度。潮流的橫向剪切大小uy可表示成2Δu（y）／B作用在一個潮周期（σ－1）內(nèi)，σ表示半日潮頻率。由公式（2）可知，隨著河口的寬度增大，由差異對流產(chǎn)生的橫向流相應(yīng)減小。由此可見，在伶仃洋河口較寬闊的特征不利于橫向環(huán)流的形成。

在此，本文主要針對剖面1枯季小潮和剖面2洪季小潮期間這兩個橫向環(huán)流結(jié)構(gòu)較為顯著的時期進(jìn)行分析。研究區(qū)域的水流結(jié)構(gòu)在整個潮周期內(nèi)的變化過程（見圖12）顯示，枯季小潮時期剖面1的水流在淺灘和深槽出現(xiàn)漲落潮不對稱現(xiàn)象，深槽的漲落潮相對滯后；表底層也存在漲落潮不對稱的現(xiàn)象，表層由于受到徑流的影響，漲潮歷時短而落潮歷時長。在落潮時期，深槽的水流向兩邊擴(kuò)散，淺灘的水流則向深槽匯聚；在漲潮時期則相反。考慮到在表底層橫向流表現(xiàn)出流向基本一致的特點，可判定淺灘的橫向流主要是由潮流的漫灘和歸槽引起，深槽的橫向流則主要是由于槽道向陸收縮所致。

如前面求取橫向流振幅的方法一樣，求得橫向動量方程中各項絕對值的面積平均值（〈｜φ｜〉，其中φ表示動量方程中的任意項），由此來表示各動量項在動量平衡中的相對重要性，得到的各橫向動量分量在潮周期內(nèi)的變化過程如圖13［7］，圖中Fu表示科氏力項，UVx表示縱向?qū)α黜?，VVy＋WVz表示橫向?qū)α黜?，Pay／ρo表示橫向正壓梯度，Pcy／ρo表示橫向斜壓梯度，AvVzz表示摩擦項。

圖12 剖面1枯季小潮（a）和剖面2洪季小潮（b）期間從漲憩開始表層（淺色）和底層（深色）平面流在潮周期內(nèi)的變化過程Fig.12 The intertidal variation of plane flows for the surface layer（light color）and the bottom layer（dark）beginning from the high slack water during the neap tide in the dry season（a）and in the flood season（b）in cross section 1 and section 2

從結(jié)果（見圖13a）中可以看出，剖面1在枯季小潮時期整個潮周期內(nèi)基本上都是正壓梯度和斜壓梯度在動量平衡中占主導(dǎo)作用，在漲急時，正壓梯度項〈｜Pay／ρo｜〉最大，斜壓梯度項〈｜Pcy／ρo｜〉和科氏力項〈｜Fu｜〉都相對較小。剖面1的水位在枯季小潮時期漲急時，表現(xiàn)出兩邊高，中間低，而在落急時則為西高東低（見圖14），由于從虎門河口流出的淡水主要位于剖面的西側(cè)，而外海傳入的潮流受到地轉(zhuǎn)偏向力的影響出現(xiàn)漲潮時偏東，因而導(dǎo)致漲潮時橫剖面的兩側(cè)水位較高，而落潮時西高東低。由此可見，由于徑流、潮流與地形的相互作用，在漲潮時，由于槽道向北收縮兩側(cè)的水流在正壓梯度力作用下向中間流動，在深槽處兩股水流匯聚，由于東側(cè)的水位梯度較大，驅(qū)動?xùn)|側(cè)的表層橫向流向西越過深槽，故在深槽的上端偏西的位置產(chǎn)生環(huán)流圈。落潮時，由于槽道向海方向變寬，水流沿槽道輻散，在深槽產(chǎn)生向兩側(cè)的橫向流分量，在表層，由于正壓梯度的作用，產(chǎn)生自西向東的橫向流，故在深槽的上層產(chǎn)生一個小的橫向環(huán)流。從圖13a可以看到，由于深槽的表底層潮汐相位的不一致，在轉(zhuǎn)流階段，同樣可以形成表底層流向不一致的橫向環(huán)流。

圖13 各動量項的絕對值（10－4）沿面積平均從漲憩開始在潮周期內(nèi)的變化過程Fig.13 The intertidal variation of the absolute value of momentum terms beginning from the high slack water

圖14 剖面1在枯季小潮漲急和落急的相對平均水位分布圖，圖中橫坐標(biāo)表示距剖面西邊界的距離Fig.14 The distribution of the relative average water level for the maximum flood and maximum ebb during the neap tide in the dry season in cross section 1，x axis denotes the distance from the western boundary of the cross section

剖面2在洪季小潮時期的結(jié)果（圖12b）與剖面1一樣，受徑流影響同樣存在淺灘與深槽漲落潮不一致以及深槽的表底層漲落潮不一致。在洪季小潮的漲急時期，縱向流速的分布（見圖15a）顯示中底層大，表層小的特點，而在落急時期（見圖15b）表層流速更大。通過與枯季時期漲急時的橫向水流結(jié)構(gòu)（見圖9a、9c）進(jìn)行比較后發(fā)現(xiàn)，在枯季時期，剖面2在漲急時橫向流整體指向西側(cè)，而在洪季，剖面2的橫向流分布（見圖10a、10c）表現(xiàn)出在剖面西側(cè)指向西，東側(cè)指向東，同時深槽底層指向西的特征。進(jìn)行動量分析后發(fā)現(xiàn)，在整個潮周期內(nèi)，正壓梯度項〈｜Pay／ρo｜〉在整個潮周期內(nèi)在動量平衡中占據(jù)主導(dǎo)地位，其變化的規(guī)律與橫向流振幅的變化規(guī)律基本一致，摩擦力項〈｜AvVzz｜〉變化趨勢與正壓梯度項基本一致，但在量值上較正壓梯度項小，主要在落潮時才起一定的作用，在漲潮期明顯減小，斜壓梯度項〈｜Pcy／ρo｜〉在整個潮周期內(nèi)都很?。▓D13b）。由于剖面2的位置靠近口門，在漲潮時期，在剖面2的西北側(cè)（見圖7），從口門流出的大量徑流受到潮汐的頂托作用出現(xiàn)壅水現(xiàn)象，水位抬高，流速減小，由此形成的正壓梯度力促使?jié)q潮流向東偏轉(zhuǎn)，并且在之后提前進(jìn)入落潮階段，而在剖面2深槽底層的潮汐受影響較小，繼續(xù)沿河道上溯，形成指向西側(cè)的橫向流分量，由于表底層水流流向的不一致從而形成橫向環(huán)流。在洪季大潮期，隨著潮汐動力的增強(qiáng)，在剖面2的東側(cè)深槽區(qū)的中底層，更多的橫向流開始指向西側(cè)。

通過上述分析發(fā)現(xiàn)，無論是剖面1的枯季小潮時期還是剖面2的洪季小潮時期形成的橫向環(huán)流，都是徑流、潮汐以及地形三者共同作用的結(jié)果。橫向流在河口動力的重要性的判定可以參照公式（3）［7］，當(dāng)橫向流能夠驅(qū)動水質(zhì)點在一個潮周期內(nèi)運動的距離超過半個河道寬度時，便可認(rèn)為橫向流在河口動力中是重要的，

式中，〈｜v｜〉表示橫向流速絕對值在橫剖面上的平均，σ表示半日潮頻率，B表示河道寬度。根據(jù)計算顯示，剖面1和剖面2的任何時刻，橫向輸送的規(guī)模都小于河口的半個寬度，可認(rèn)為橫向流在伶仃洋河口動力中相對較為次要。

圖15 剖面2洪季小潮漲急（a）和落急時期（b）縱向流速分布圖，圖中等值線表示縱向流速，單位為m／sFig.15 The distribution of the vertical velocity for the maximum flood（a）and maximum ebb（b）during the neap tide in the flood season in cross section 2，and contours denote velocity with the unit of m／s

6 結(jié)論

（1）伶仃洋河口的地形變化，包括河口的向陸收縮以及槽道的彎曲等使得水流的流路復(fù)雜，影響橫向流的結(jié)構(gòu)，使其具有明顯的局地特征。

（2）在不同的徑流和潮汐條件下，伶仃洋河口的橫向流隨著漲落潮變化而變化，環(huán)流結(jié)構(gòu)只在部分時間存在，如剖面1的枯季時期以及剖面2的洪季小潮期的部分區(qū)域，總體上橫向流表現(xiàn)出表底層流向一致的特征。

（3）伶仃洋河口的橫向環(huán)流結(jié)構(gòu)與經(jīng)典的差異對流產(chǎn)生的斜壓驅(qū)動以及地轉(zhuǎn)偏向力驅(qū)動的環(huán)流結(jié)構(gòu)存在明顯的差異。經(jīng)研究分析，其形成機(jī)理主要是徑流、潮流以及地形三者之間的相互作用形成的，具有明顯的局地特征。

（4）伶仃洋內(nèi)的橫向環(huán)流相對較弱，在伶仃洋的河口動力中較為次要。

［1］ Smith R.Longitudinal dispersion of a buoyant contaminant in a shallow channel［J］.Journal of Fluid Mechanics，1976，78（04）：677－688.

［2］ West J R，Mangat J S.The determination and prediction of longitudinal dispersion coefficients in a narrow，shallow estuary［J］.Estuarine，Coastal and Shelf Science，1986，22（2）：161－181.

［3］ Guymer I，West J R.Longitudinal dispersion coefficients in estuary［J］.Journal of Hydraulic Engineering，1992，118（5）：718－734.

［4］ Seim H E，Gregg M C.The importance of aspiration and channel curvature in producing strong vertical mixing over a sill［J］.Journal of Geophysical Research：Oceans（1978－2012），1997，102（C2）：3451－3472.

［5］ Geyer W R，Signell R P，Kineke G C.Lateral trapping of sediment in a partially mixed estuary［C］／／Physics of Estuariesand Coastal Seas：Proceedings of the 8th International Biennial Conference on Physics of Estuaries and Coastal Seas.1998：115－124.

［6］ Nidzieko N J，Hench J L，Monismith SG.Lateral circulation in well-mixed and stratified estuarine flows with curvature［J］.Journal of Physical O-ceanography，2009，39（4）：831－851.

［7］ Lerczak J A，Rockwell G W.Modeling the lateral circulation in straight，stratified estuaries［J］.Journal of Physical Oceanography，2004，34（6）：1410－1428.

［8］ Chen S N，Sanford L P.Lateral circulation driven by boundary mixing and the associated transport of sediments in idealized partially mixed estuaries［J］.Continental Shelf Research，2009，29（1）：101－118.

［9］ Chen S N，Sanford L P，Ralston D K.Lateral circulation and sediment transport driven by axial winds in an idealized，partially mixed estuary［J］.Journal of Geophysical Research：Oceans（1978－2012），2009，114（C12）：C12006.

［10］ Valle-Levinson A，Li C，Wong K C，et al.Convergence of lateral flow along a coastal plain estuary［J］.Journal of Geophysical Research：Oceans（1978－2012），2000，105（C7）：17045－17061.

［11］ Wong L A，Chen J C，Xue H，et al.A model study of the circulation in the Pearl River Estuary（PRE）and its adjacent coastal waters：1.Simulations and comparison with observations［J］.Journal of Geophysical Research：Oceans（1978－2012），2003，108（C5）：3156.

［12］ 袁麗蓉，楊清書，謝莉莉，等.人類活動干擾下伶仃洋河口灣地形演變趨勢變化［C］／／第十六屆中國海洋（岸）工程學(xué)術(shù)討論會（下冊）.2013.

Yuan Lirong，Yang Qingshu，Xie Lili，et al.Trend of the Lingding Estuary topographic evolution under human interference［C］／／The 16th China Ocean（Shore）Engineering Symposium（the 2nd volume）.2013.

［13］ 應(yīng)秩甫，陳世光.珠江口伶仃洋咸淡水混合特征［J］.海洋學(xué)報，1983，5（1）：1－10.

Ying Zhifu，Chen Shiguang.The characteristics of saline and fresh water mixing in the Lingding Estuary［J］.Haiyang Xuebao，1983，5（1）：1－10.

［14］ 趙煥庭.珠江河口灣伶仃洋的地形［J］.海洋學(xué)報，1981，3（2）：255－274.

Zhao Huanting.The topography of the Lingding，Pearl River Estuary［J］.Haiyang Xuebao，1981，3（2）：255－274.

［15］ 應(yīng)秩甫，伶仃洋橫向動力平衡［M］／／珠江口海岸帶和灘涂資源綜合調(diào)查研究文集（四）.廣州：廣東科技出版社，1986.

Ying Zhifu.The lateral dynamic balance of the Lingding Estuary［M］／／The Comprehensive Investigation and Research Anthology of the Coastal Zone and Salt Marsh Resources in the Pearl River Estuary（4th）.Guangzhou：Guangdong Science and Technology Press，1986.

［16］ Hamrick J M.User's manual for the environmental fluid dynamics computer code［R］.Special Report NO.331 in Applied Marine Science and O-cean Engineering，Virginia Institute of Marine Science，College of William and Mary，1996.

［17］ Hamrick J M.A three-dimensional environmental fluid dynamics computer code：the oretical and computational aspects［R］.Special Report 337，Virginia Institute of Marine Science，College of William and Mary，1992.

Lateral circulation in the Lingding Estuary

Yi Kan1，Gong Wenping1
（1.Center for Coastal Ocean Science and Technology Research，School of Marine Science，Sun Yat-sen University，Guangzhou 510275，China）

Lateral circulation，which can be generated in estuarine systems by a variety of mechanisms，plays an important role in the momentum balance and mass transport in estuaries.A three-dimensional hydrodynamic model（EFDC），combined with field measurement data in 1989 and 2002，is used to investigate the lateral circulation in the Lingding Estuary.The model have successfully reproduced the hydrodynamics and salinity distribution during the flood and dry seasons.Two cross-sectional profiles in two different regions are chosen to examine the lateral circulation and salinity dynamics in different periods.The results show that lateral circulation is poorly developed in the Lingding Estuary and features more local characteristics because of its special bathymetry and hydrodynamic conditions.The lateral circulationin the Lingding Estuary is generated by theinteraction among river flow，tidal currents and bathymetry，and the mechanisms in two sections are different.In this study，we conclude that lateral circulation is not significant for the dynamics of Lingding Estuary.

Lingding Estuary；EFDC model；lateral circulation；estuarine hydrodynamics

P731.2

A

0253-4193（2015）03-0001-14

，龔文平.伶仃洋河口橫向環(huán)流［J］.海洋學(xué)報，2015，37（3）：1—14，

10.3969／j.issn.0253-4193.2015.03.001

Yi Kan，Gong Wenping.Lateral circulation in the Lingding Estuary［J］.Haiyang Xuebao，2015，37（3）：1—14，doi：10.3969／j.issn.0253-4193.2015.03.001

2014-05-12；

2014-10-03。

國家自然科學(xué)基金（41061130542）。

易侃（1992—），男，江西省宜春市人，從事河口動力學(xué)研究。E-mail：yikan302＠foxmail.com

*通信作者：龔文平。E-mail：gongwp＠m(xù)ail.sysu.edu.cn