陳涇，朱建榮*

（1.華東師范大學(xué)河口海岸學(xué)國家重點(diǎn)實(shí)驗(yàn)室，上海 200062）

長江河口青草沙水庫鹽水入侵來源

陳涇1，朱建榮1*

（1.華東師范大學(xué)河口海岸學(xué)國家重點(diǎn)實(shí)驗(yàn)室，上海 200062）

應(yīng)用改進(jìn)的三維數(shù)值模式ECOM-si，從模式計(jì)算的鹽度和流向的變化過程、漲憩和落憩時(shí)刻鹽度等值線和淡水區(qū)域的變化，分析在一般動力條件下青草沙水庫取水口鹽水入侵來源。計(jì)算結(jié)果表明，小潮后中潮、大潮、大潮后中潮和小潮期間北支倒灌占青草沙水庫取水口表層鹽水入侵比例分別為69.5%、89.3%、98.5%和99.5%，占底層鹽水入侵比例分別為34.9%、88.9%、98.5%和99.5%。除了小潮后中潮期間底層鹽水入侵來源主要來自下游外海（占65.1%），青草沙水庫取水口表層和底層鹽水入侵來源主要來自北支鹽水倒灌，尤其是大潮后中潮和小潮期間幾乎全部來自北支鹽水倒灌。

青草沙水庫；鹽水入侵；來源；數(shù)值計(jì)算

1 引言

長江河口地區(qū)，經(jīng)濟(jì)發(fā)達(dá)，人口密集，需要大量優(yōu)質(zhì)原水。2010年前上海的用水主要取自黃浦江，水量不足，水質(zhì)較差，為典型的水質(zhì)性缺水城市。長江河口水量充沛，水質(zhì)優(yōu)良，要從根本上解決上海的用水難題，必須從長江河口取水。為了解決用水難題，上海市在長江口建成了大型河口江心水庫——青草沙水庫。青草沙水庫位于南北港分汊口附近、長興島西北側(cè)（見圖1），水域面積達(dá)到66.15 km2，于2010年建成并開始向上海供水，承擔(dān)了上海市約50%的原水供應(yīng)，規(guī)劃供水規(guī)模為7.19×106m3／d，受益人口超過1 000萬人［1］。青草沙水庫的興建，大大緩解了上海供水緊張的局面，為上海城市發(fā)展和工農(nóng)業(yè)生產(chǎn)提供有利的保障。但從長江河口取水面臨的主要問題是枯季鹽水入侵，如果取水口鹽度大于飲用水標(biāo)準(zhǔn)0.45（實(shí)用鹽度單位，按國際慣例單位一般不標(biāo)注）就不能取水。

長江河口枯季經(jīng)常面臨鹽水入侵，最大的特點(diǎn)是南支除受外海鹽水入侵外，還受上游北支鹽水倒灌的入侵［2］。已有的大量觀測和研究表明，潮汐和徑流量是影響鹽水入侵的主要原因，另外還受風(fēng)應(yīng)力［3—4］、口外陸架環(huán)流［5］和河勢變化［2］等影響。以往對長江河口鹽水入侵來源的研究，沈煥庭等［2］指出潮汐河口外海鹽水入侵導(dǎo)致的鹽度周日變化一般規(guī)律是鹽度最高、最低值分別出現(xiàn)在漲憩、落憩附近，且周日變幅大。長江口南支、南北港受北支鹽水倒灌后，改變了鹽度周日變化規(guī)律，具體反映為鹽度周日變幅小，周日的鹽度峰值和谷值出現(xiàn)在落憩和漲憩附近。茅志昌等［6］通過分析現(xiàn)場觀測資料，指出青草沙水源地的鹽水源自外海鹽水入侵和北支鹽水倒灌，其中以受北支倒灌鹽水團(tuán)過境作用為主。顧玉亮等［7］根據(jù)觀測資料分析，指出影響南支水源地的鹽水來源有兩個(gè)，即北支鹽水倒灌和南北港外海鹽水入侵。樂勤等［1］根據(jù)多年監(jiān)測數(shù)據(jù)指出，影響青草沙水庫水域氯化物變化的咸潮入侵源有3個(gè)，即北支咸潮倒灌、南港咸潮入侵和北港咸潮入侵，其中最主要的入侵因素是北支咸潮倒灌。以往對長江河口鹽水入侵的研究，主要針對陳行水庫鹽水入侵來源。陳行水庫于1992年建成啟用，運(yùn)營時(shí)間較早，對其研究較多。青草沙水庫于2010年建成啟用，運(yùn)營時(shí)間較晚，對鹽水入侵來源研究相對較少。開展青草沙水庫鹽水入侵來源的研究，定量給出不同潮型下北港外海鹽水入侵和上游北支鹽水倒灌在水庫取水口鹽水來源的百分比，可為水庫避咸蓄淡提供科技依據(jù)。

本文應(yīng)用長江河口鹽水入侵三維數(shù)值模式，在模式驗(yàn)證的基礎(chǔ)上，開展青草沙水庫鹽水入侵來源的研究。

圖1 長江河口形勢圖和測站位置Fig.1 Map of the Changjiang Estuary and the distribution of observation stations

2 三維數(shù)值模式的設(shè)置和驗(yàn)證

2.1 模式設(shè)置

本文應(yīng)用改進(jìn)的三維數(shù)值模式ECOM-si［8—10］，該模式長期應(yīng)用于長江河口地區(qū)水動力過程和鹽水入侵等方面的研究，并取得諸多成果［3，11—14］。

模式采用水平曲線非正交網(wǎng)格，范圍包括整個(gè)長江河口、杭州灣和鄰近海區(qū)，上游邊界設(shè)在長江枯季潮區(qū)界大通，外海開邊界東邊到124.5°E附近，北邊到33°N附近，南邊到28°N附近（見圖2）。對長江河口區(qū)域，包括南北支分汊口和深水航道工程區(qū)域的網(wǎng)格進(jìn)行局部加密，并且較好地?cái)M合了岸線和導(dǎo)堤?？趦?nèi)網(wǎng)格分辨率為100～500 m不等，口外網(wǎng)格較疏，分辨率最大為10 km左右。垂向采用σ坐標(biāo)均勻分為10層，時(shí)間步長取60 s。長江河口區(qū)域淺灘較多，模式運(yùn)用干濕判別法實(shí)現(xiàn)潮灘移動邊界的模擬，臨界水深取0.2 m。模式地形采用2010年岸線及水深資料。外海開邊界由潮位驅(qū)動，考慮16個(gè)分潮（M2，S2，N2，K2，K1，O1，P1，Q1，MU2，NU2，T2，L2，2N2，J1，M1，OO1），由各分潮調(diào)和常數(shù)合成得到。分潮調(diào)和常數(shù)由全球潮汐數(shù)值模式NAOTIDE計(jì)算的結(jié)果得到（http：//www.miz.nao.ac.jp／）。初始水位和流速取零。關(guān)于初始鹽度場，從大通到徐六涇620 km的河道，均為淡水，鹽度為0，從徐六涇到口門由枯季各月多次實(shí)測資料插值得到，在長江口外由《渤海黃海東海海洋圖集（水文）》各月分布圖數(shù)字化得到［15］，考慮海表面風(fēng)應(yīng)力的作用。

圖2 模式計(jì)算區(qū)域和網(wǎng)格（a），南北支分汊口及其附近水域網(wǎng)格（b），口門深水航道及其鄰近水域網(wǎng)格（c）Fig.2 Domain and mesh for numerical model（a），enlarged views of the model mesh around the bifurcation of the North Branch and the South Branch（b），and the Deep Waterway Project（c）

2.2 模式驗(yàn)證

改進(jìn)的ECOM-si模式已在長江河口進(jìn)行了大量的驗(yàn)證，模式計(jì)算的水位、流速、流向和鹽度與實(shí)測資料吻合良好［3，10，14］。本文采用2011年12月24日至2012年1月13日長江口觀測資料對模式作進(jìn)一步驗(yàn)證，測站的分布圖見圖1。模式從2011年11月1日起算，徑流量采用大通水文站每日的實(shí)測值，風(fēng)場采用QSCAT／NCEP提供的時(shí)間分辨率為6 h、空間分辨率為0.5°×0.5°的數(shù)值產(chǎn)品。

圖3和圖4分別為船測站點(diǎn)B和C表底層流速、流向和鹽度觀測值和模擬計(jì)算值隨時(shí)間變化，測點(diǎn)B位于南北支分汊口附近，測點(diǎn)C位于北港北汊口。測點(diǎn)B實(shí)測最大流速表層約為1.4 m／s，底層約為1 m／s，因底摩擦作用表層流速大于底層流速，模式計(jì)算流速大于實(shí)測值。落潮流歷時(shí)大于漲潮流歷時(shí)，流向同實(shí)測流向一致。觀測表層鹽度小于底層鹽度，計(jì)算鹽度與觀測鹽度較為接近。測點(diǎn)C模式計(jì)算流速、流向和鹽度與實(shí)測值吻合良好，漲憩時(shí)刻鹽度達(dá)到最大，落憩時(shí)刻鹽度達(dá)到最小，表明鹽水入侵源自下游外海。同樣，落潮流歷時(shí)大于漲潮流歷時(shí)，底層鹽度大于表層鹽度。

圖5為浮筒和水文站觀測和模式計(jì)算表層鹽度隨時(shí)間變化。浮筒測點(diǎn)A位于北支口，鹽度在25～30之間波動，模式計(jì)算值略高于實(shí)測值。浮筒測點(diǎn)D位于北港攔門沙區(qū)域，該處存在鹽度鋒面，鹽度隨時(shí)間變化大，模式計(jì)算值略大于實(shí)測值。浮筒測點(diǎn)青草沙水庫取水口在大部分時(shí)間鹽度低于飲用水標(biāo)準(zhǔn)0.45，在2012年1月6－10日出現(xiàn)大的波動和峰值，模式很好地模擬出了這個(gè)過程。在堡鎮(zhèn)水文站，在2011年12月24至2012年1月4日期間，鹽度在0.4上下變化，波動很小。在1月5－13日之間鹽度出現(xiàn)大的波動和峰值，實(shí)測鹽度最大值達(dá)到9.0，同樣模式很好地再現(xiàn)了這個(gè)過程。

從上面的驗(yàn)證結(jié)果可以看出模式計(jì)算的流速、流向和鹽度同實(shí)測資料吻合良好，可以較好地模擬長江河口水動力和鹽水入侵過程。

2.3 數(shù)值試驗(yàn)

本文研究青草沙水庫鹽水入侵來源，影響水庫取水口鹽水入侵的因素有徑流量、潮汐和風(fēng)應(yīng)力等。先設(shè)計(jì)一個(gè)徑流量和風(fēng)況為冬季1－2月一般狀況的控制數(shù)值試驗(yàn)，模擬長江河口和水庫取水口鹽水入侵狀況。數(shù)值模式計(jì)算時(shí)段為鹽水入侵嚴(yán)重的1－2月，模式從1月1日開始運(yùn)行，至2月底結(jié)束，輸出2月的計(jì)算結(jié)果作分析和比較。徑流量取大通1950年以來1月和2月平均值，分別為1.11×104m3／s和1.2 ×104m3／s。風(fēng)場取NCEP多年半月平均風(fēng)場，外海開邊界考慮16個(gè)主要分潮。

為了解在不同動力條件下青草沙水庫鹽水入侵狀況和來源，設(shè)計(jì)在北支上段封堵的數(shù)值試驗(yàn)，將沒有北支鹽水倒灌情況下計(jì)算結(jié)果與控制試驗(yàn)結(jié)果比較，定量分析青草沙水庫鹽水入侵來源。

圖3 船只測點(diǎn)B處流速（a、b）、流向（c、d）和鹽度（e、f）隨時(shí)間變化Fig.3 Variations of current speed（a、b），direction（c、d）and salinity（e、f）from17：30 on January 11，2012 to 7：00 on January 13，2012 at the station B

圖4 船只測點(diǎn)C處流速（a、b）、流向（c、d）和鹽度（e、f）隨時(shí)間變化Fig.4 Variations of current speed（a、b），direction（c、d）and salinity（e、f）from 12：30 on January 10，2012 to 7：30 on January 13，2012 at the station C

圖5 浮筒和水文站測點(diǎn)2011年12月24日0∶00至2012年1月13日24∶00表層鹽度隨時(shí)間變化Fig.5 Variations of surface salinity of the observed data and model results at the buoy stations and hydrometric stations from 00∶00 on December 24，2011 to 24∶00 on January 13，2012.

3 結(jié)果與分析

3.1 控制試驗(yàn)

圖6為青草沙水庫取水口水位、流速、流向和表層、底層鹽度隨時(shí)間變化，潮汐為半日潮，日不等現(xiàn)象明顯，半月大小潮變化顯著。大潮最高潮位約2.1 m，最低潮位約－1.1 m（黃海85基面）。水流為往復(fù)流、表層最大流速約130 cm／s。取水口大部分時(shí)間鹽度低于0.45，2月24日至3月2日鹽度大于0.45，水庫不宜取水。在2月4－8日小潮后中潮期間底層鹽度出現(xiàn)大的峰值，最大值達(dá)到4.8，量值遠(yuǎn)大于表層鹽度。峰值出現(xiàn)在漲憩時(shí)刻，表明取水口底層鹽水入侵來自北港外海。除了小潮后中潮期間表底層鹽度差異顯著，其他時(shí)段差異微小，表明垂向混合均勻。

為清楚地表示不同潮型期間取水口鹽水來源，將2月11－26日表層流向和鹽度放大（見圖7）。在大潮后中潮期間11－14日（結(jié)合圖6中水位過程），前2天落憩時(shí)刻鹽度達(dá)到峰值，漲憩時(shí)刻達(dá)到谷值，表明鹽水入侵來自上游，即北支鹽水倒灌。后2天鹽度大于0.45，但波動很小，表明倒灌鹽水團(tuán)正處于取水口水域，水平混合均勻。在小潮期間15－18日，絕大部分鹽度小于0.45，水庫能取淡水，鹽度波動小，但峰值出現(xiàn)在漲憩時(shí)刻，表明鹽水入侵來自下游外海。但該低鹽水源自北支倒灌、小潮期間移至取水口下游，故本質(zhì)上還是源自上游北支倒灌。在小潮后中潮期間19－22日，鹽度小于0.45，鹽度峰值出現(xiàn)在漲憩，鹽水入侵來自下游外海。在大潮期間23－26日，鹽度隨時(shí)間逐漸升高，鹽度峰值出現(xiàn)在落憩時(shí)刻，鹽度入侵來自上游北支倒灌。因此，大潮和大潮后中潮期，水庫取水口鹽水入侵來自上游北支鹽水倒灌，小潮和小潮后中潮期來自下游外海鹽水入侵。

圖6 控制數(shù)值試驗(yàn)計(jì)算的青草沙水庫取水口2月1日至3月2日水位、表層流速和流向、鹽度隨時(shí)間變化Fig.6 Temporal variation of water level，surface current speed and direction，and salinity from February 1 to March 2 at the water intake of Qingcaosha Reservoir，calculated by the control experiment

2月4－7日小潮后中潮期間漲憩時(shí)刻（t1）、落憩時(shí)刻（t2）表層和底層鹽度平面圖可見圖8，在青草沙水庫附近北港水域鹽度從口外向口內(nèi)遞減，北港鹽度鋒面明顯。在漲憩時(shí)刻，表層等鹽度線0.45正在經(jīng)過取水口，上游為鹽度低于0.45的淡水，下游為鹽水，底層北港北汊高鹽水入侵北港十分明顯，鹽度顯著升高，取水口下游附近出現(xiàn)等鹽度線3，表明時(shí)刻底層鹽水入侵強(qiáng)烈，這與圖5給出的小潮后中期期間底層高鹽度峰值是一致的。在落憩時(shí)刻，等鹽度線0.45向下游移動，取水口附近淡水面積增大，上游倒灌的鹽水距取水口尚遠(yuǎn)。故從取水口漲憩時(shí)刻鹽度增大、淡水面積減少，落憩時(shí)刻鹽度減小、淡水面積增加的鹽度平面分布，可得出小潮后中潮期間取水口鹽水入侵來自下游外海。

圖7 局部時(shí)段放大的青草沙取水口表層流向和鹽度隨時(shí)間變化（2月11－26日）Fig.7 A zoomed view of temporal variation of surface salinity and current direction at the water intake of Qingcaosha Reservoir from February 11 to 26

圖8 小潮后中潮期間漲憩時(shí)刻（t1）表層（a）、底層（b），落憩時(shí)刻（t2）表層（c）、底層（d）鹽度平面分布Fig.8 Salinity distributions at surface（a）and bottom（b）layer at flood slack（t1），and at surface（c）and bottom（d）at ebb slack（t2）in moderate tide following neap tide

在大潮期間2月25日漲憩和落憩時(shí)刻（圖9），北支倒灌進(jìn)入南支的鹽度高于1.0鹽水位于南支上段北側(cè)，鹽度低于0.45的淡水位于橫沙小港北側(cè)附近水域，在漲憩和落憩時(shí)刻取水口鹽度均大于0.45，不宜取水。從取水口鹽度變化看，漲憩時(shí)刻下游淡水區(qū)域趨近，落憩時(shí)刻下游淡水區(qū)域遠(yuǎn)離，而上游高鹽水趨近，所以大潮期間取水口鹽水入侵源自上游北支倒灌。這與上面從鹽度變化過線分析結(jié)果一致。

圖9 大潮期間漲憩（t3）表層（a）、底層（b），落憩（t4）表層（c）、底層（d）鹽度平面分布Fig.9 Salinity distributions at surface（a）and bottom（b）layer at flood slack（t3），and at surface（c）and bottom（d）at ebb slack（t4）during spring tide

3.2 北支鹽水倒灌影響

在控制實(shí)驗(yàn)的基礎(chǔ)上，封堵北支上段，完全消除北支倒灌的影響，模擬和對比分析青草沙水庫鹽水入侵。圖10為封堵北支上段后青草沙水庫取水口表層和底層鹽度隨時(shí)間變化過程，對比未封堵時(shí)情況（見圖6），除了2月4－7日底層鹽度仍出現(xiàn)峰值、量值基本一致外，其他時(shí)段鹽度大幅下降。這表明取水口表層鹽水入侵來自北支鹽水倒灌，底層小潮后中潮來自下游外海，其他時(shí)段主要來自上游北支鹽水倒灌。

為定量分析北支倒灌在青草沙水庫取水口附近水域鹽水入侵中的比例，我們計(jì)算了沿北港橫斷面不同潮型期間漲潮和落潮總的平均鹽通量（為絕對值之和平均，見表1），在控制數(shù)值試驗(yàn)中小潮后中潮、大潮、大潮后中潮和小潮期間平均鹽通量分別為20.66、20.76、34.11和15.28 t／s，封堵北支上段后分別變?yōu)?2.08、2.80、0.11和0.00 t／s。通過北港橫斷面的鹽通量由來自上游北支鹽水倒灌和下游外海鹽水入侵共同引起，那么封堵北支上段后鹽通量的減少部分可認(rèn)為是由北支倒灌產(chǎn)生，從整個(gè)斷面角度看北支倒灌占鹽水入侵的百分比在小潮后中潮、大潮、大潮后中潮和小潮期間分別為41.5%、86.5%、99.7%和100%，也即除了在小潮后中潮期間青草沙水庫取水口鹽水入侵58.5%來自下游外海，其他潮型主要來自北支鹽水倒灌，尤其是大潮后中潮和小潮幾乎全部來自北支鹽水倒灌。

圖10 封堵北支上段情況下青草沙水庫取水口2月1日至3月2日表層鹽度（黑線）、底層鹽度（紅線）隨時(shí)間變化Fig.10 Variation of surface（black line）and bottom（red dashed line）salinity at the water intake of Qingcaosha Reservoir from February 1 to March 2，after the North Branch was blocked

從青草沙水庫取水口表層和底層鹽度變化看，小潮后中潮、大潮、大潮后中潮和小潮期間北支倒灌占表層鹽水入侵比例分別為69.5%、89.3%、98.5%和99.5%，占底層鹽水入侵比例分別為34.9%、88.9%、98.5%和99.5%（表2），也就是說青草沙水庫取水口表層鹽水入侵來源主要來自北支鹽水倒灌，尤其是大潮后中潮和小潮期間幾乎全部來自北支鹽水倒灌。底層鹽水入侵來源小潮后中潮期間主要來自下游外海，占65.1%，其他潮型主要來自北支鹽水倒灌，尤其是大潮后中潮和小潮期間幾乎全部來自北支鹽水倒灌。長江河口鹽水入侵最大特點(diǎn)是北支鹽水倒灌，嚴(yán)重影響南支鹽水入侵和水源地。因?yàn)榍嗖萆乘畮烊∷诰嚯x北港口門較遠(yuǎn)，北支鹽水倒灌成為了除小潮后中潮期間外鹽水入侵外的主要來源。

表1 北港斷面不同潮型期間平均鹽通量（t／s）Tab.1 Averaged salinity fluxes in different tide patterns at the cross section of the North Channel

表2 青草沙水庫取水口不同潮型期間平均鹽度Tab.2 Averaged salinity in different tide patterns at the water intake of the Qingcaosha Reservoir

4 結(jié)論

應(yīng)用改進(jìn)的三維數(shù)值模式ECOM-si，考慮徑流量、潮汐和風(fēng)應(yīng)力等的作用，研究青草沙水庫取水口鹽水入侵的來源。采用2011年12月至2012年1月長江口觀測資料對模式作流速、流向和鹽度的驗(yàn)證，模式計(jì)算結(jié)果與實(shí)測值吻合良好。

設(shè)計(jì)控制數(shù)值試驗(yàn)，根據(jù)鹽度和流向的變化過程、漲憩和落憩時(shí)刻鹽度等值線和淡水區(qū)域的變化，分析在一般動力條件下青草沙水庫取水口鹽水入侵來源。設(shè)計(jì)封堵北支上段數(shù)值試驗(yàn)，計(jì)算不同潮型下通過北港橫斷面的漲潮和落潮平均鹽通量和取水口表層和底層平均鹽度，并與控制數(shù)值試驗(yàn)結(jié)果比較，定量得出鹽水入侵來源的百分比。

從通過北港斷面鹽通量計(jì)算結(jié)果看，青草沙水庫取水口附近水域鹽水入侵中北支倒灌所占百分比在小潮后中潮、大潮、大潮后中潮和小潮期間分別為41.5%、86.5%、99.7%和100%。從青草沙水庫取水口表層和底層鹽度變化看，小潮后中潮、大潮、大潮后中潮和小潮期間北支倒灌占表層鹽水入侵比例分別為69.5%、89.3%、98.5%和99.5%，占底層鹽水入侵比例分別為34.9%、88.9%、98.5%和99.5%。表層鹽水入侵來源主要來自北支鹽水倒灌，底層鹽水入侵來源小潮后中潮期間主要來自下游外海，占65.1%。除了小潮后中潮，其他潮型表層和底層鹽水入侵來源主要來自北支鹽水倒灌，尤其是大潮后中潮和小潮期間幾乎全部來自北支鹽水倒灌。

需要指出的是本文結(jié)果是在一般動力因子情況下模式結(jié)算的結(jié)果，若冬季出現(xiàn)偏北大風(fēng)，北港的鹽水入侵會大幅加強(qiáng)，期間鹽水入侵來源可能主要來自下游外海。

［1］樂勤，關(guān)許為，劉小梅，等.青草沙水庫取水口選址與取水方式研究［J］.給水排水，2009（2）：46－51.

［2］沈煥庭，茅志昌，朱建榮.長江河口鹽水入侵［M］.北京：海洋出版社，2003：16.

［3］Li L，Zhu J R，Wu H.Impacts of wind stress on saltwater intrusion in the Yangtze Estuary［J］.Science China Earth Sciences，2012，55（7）：1178－1192.

［4］朱建榮，傅利輝，吳輝.風(fēng)應(yīng)力和科氏力對長江河口沒冒沙淡水帶的影響［J］.華東師范大學(xué)學(xué)報(bào)：自然科學(xué)版，2008（6）：1－8.

［5］項(xiàng)印玉，朱建榮，吳輝.冬季陸架環(huán)流對長江河口鹽水入侵的影響［J］.自然科學(xué)進(jìn)展，2009，19（2）：192－202.

［6］茅志昌，沈煥庭，肖成獻(xiàn).長江口北支鹽水倒灌南支對青草沙水源地的影響［J］.海洋與湖沼，2001，32（1）：58－66.

［7］顧玉亮，吳守培，樂勤.北支鹽水入侵對長江口水源地影響研究［J］.人民長江，2003，34（4）：1－4.

［8］朱建榮.海洋數(shù)值計(jì)算方法和數(shù)值模式［M］.北京：海洋出版社，2003：72－126.

［9］朱建榮，朱首賢.ECOM模式的改進(jìn)及在長江河口、杭州灣及鄰近海區(qū)的應(yīng)用［J］.海洋與湖沼，2003，34（4）：364－374.

［10］Wu H，Zhu J.Advection scheme with 3rd high-order spatial interpolation at the middle temporal level and its application to saltwater intrusion in the Changjiang Estuary［J］.Ocean Modelling，2010，33（1）：33－51.

［11］Wu H，Zhu J，Chen B，et al.Quantitative relationship of runoff and tide to saltwater spilling over from the North Branch in the Changjiang Estuary：A numerical study［J］.Estuarine，Coastal and Shelf Science，2006，69（1）：125－132.

［12］Chen B，Zhu J，F(xiàn)u L.Formation mechanism of freshwater zone around the Meimao Sandbank in the Changjiang estuary［J］.Chinese Journal of O-ceanology and Limnology，2010，28：1329－1339.

［13］Li L，Zhu J，Wu H，et al.A numerical study on water diversion ratio of the Changjiang（Yangtze）estuary in dry season［J］.Chinese Journal of O-ceanology and Limnology，2010，28：700－712.

［14］Wu H，Zhu J，Ho Choi B.Links between saltwater intrusion and subtidal circulation in the Changjiang Estuary：A model-guided study［J］.Continental Shelf Research，2010，30（17）：1891－1905.

［15］海洋圖集編委會.渤海黃海東海海洋圖集（水文）［M］.北京：海洋出版社，1992：13－168.

Sources for saltwater intrusion at the water intake of Qingcaosha Reservoir in the Changjiang Estuary

Chen Jing1，Zhu Jianrong1
（1.State Key Laboratory of Estuarine and Coastal Research，East China Normal University，Shanghai 200062，China）

Based on temporal variation of the simulated salinity and current direction，change of the isohaline and freshwater zone at flood slack and ebb slack，this study used the modified numerical 3-D model ECOM-si to analyze sources of saltwater intrusion at the water intake of Qingcaosha Reservoir under normal dynamic condition.The modeled results showed that saltwater-spill-over from the North Branch into the South Branch（SSO）accounts for 69.5%，89.3%，98.5%and 99.5%of the surface saltwater intrusion at the water intake of the Qingcaosha Reservoir during moderate tide following neap tide，spring tide，moderate tide following spring tide and neap tide，respectively，and 34.9%，88.9%，98.5%，99.5%of the bottom saltwater intrusion at the water intake of the reservoir during the same four-type tides mentioned above.Except for the bottom saltwater intrusion mainly from downstream open sea（accounts for 65.1%）during moderate tide following neap tide，the source of the saltwater intrusion at surface and bottom layer at the water intake of the reservoir is mainly from the SSO，especially during moderate tide following spring and neap tide.

Qingcaosha Reservoir；saltwater intrusion；source；numerical simulation

P731.12

A

0253-4193（2014）11-0131-11

2013-08-30；

2014-06-05。

國家自然科學(xué)基金項(xiàng)目（41176071）；上海市科學(xué)技術(shù)委員會重大項(xiàng)目（12231203101）；水利部公益性項(xiàng)目（201201068）。

陳涇（1988－），男，福建省龍海市人，從事河口海岸動力學(xué)研究。E-mail：chenjingfjxm＠sina.com

*通信作者：朱建榮，男，教授，從事河口海洋動力學(xué)研究。E-mail：jrzhu＠sklec.ecnu.edu.cn

陳涇，朱建榮.長江河口青草沙水庫鹽水入侵來源［J］.海洋學(xué)報(bào)，2014，36（11）：131－141，

10.3969／j.issn.0253-4193.2014.11.015

Chen Jing，Zhu Jianrong.Sources for saltwater intrusion at the water intake of Qingcaosha Reservoir in the Changjiang Estuary［J］.Acta Oceanologica Sinica（in Chinese），2014，36（11）：131－141，doi：10.3969／j.issn.0253-4193.2014.11.015