丁 磊,陳黎明,高祥宇,繳 健,胡 靜
(1. 南京水利科學(xué)研究院港口航道泥沙工程交通行業(yè)重點(diǎn)實(shí)驗(yàn)室,水文水資源與水利工程科學(xué)國家重點(diǎn)實(shí)驗(yàn)室,江蘇南京 210029; 2. 南京市水利規(guī)劃設(shè)計(jì)院股份有限公司,江蘇南京 210006)
圖1 長江口水源地位置Fig.1 Positions of water sources in Yangtze River estuary
長江河口是我國最大的河口,上起徐六涇,在平面上呈三級分汊、四口入海的河勢格局。崇明島將長江口分為北支和南支,長興島和橫沙島將南支分為北港和南港,九段沙將南港分為北槽和南槽(圖1)。長江口水量豐沛,水體自凈能力強(qiáng),為沿岸居民提供了豐富的淡水資源。從1883年開始,上海市取水口位于黃浦江。由于黃浦江處于太湖流域下游,隨著時(shí)間推移,上游污染影響下水質(zhì)較差且不穩(wěn)定性的問題日益凸顯。因此,有必要尋找新水源地。從20世紀(jì)90年代開始,長江口三大水源地——陳行水庫、青草沙水庫和東風(fēng)西沙水庫相繼建成(圖1),目前供水規(guī)模已達(dá)上海市的80%左右。因此長江口水源地安全對上海市經(jīng)濟(jì)發(fā)展和人民生活產(chǎn)生重要影響。
由于三大水源地取水口均位于河口地區(qū),枯季鹽水入侵成為水源地安全最主要的威脅。根據(jù)《地表水環(huán)境質(zhì)量標(biāo)準(zhǔn)》(GB 3838—2002)及相關(guān)規(guī)范要求,鹽度超過0.45‰時(shí)不可作為飲用水。徑流和潮汐是影響河口鹽水入侵的主要因素,20世紀(jì)80年代初的研究即已經(jīng)涉及徑流和潮汐對鹽度分布的影響[1]。安徽大通水文站距離長江口600多千米,是距離長江口最近且不受潮汐影響的水文站,因此大通的徑流量可表示長江口徑流特征[2]。羅小峰等[3]通過不同徑潮組合的數(shù)值模擬表明,徑流直接影響鹽水入侵距離、潮汐影響鹽水回蕩范圍。侯成程等[4]研究了長江口鹽水入侵對大通枯季徑流量變化的響應(yīng)時(shí)間。Qiu等[5]研究了漲落潮、大小潮以及潮汐季節(jié)變化對長江口鹽水入侵影響。長江口三級分汊四口入海的河勢格局又使得鹽水入侵的時(shí)空分布更為復(fù)雜[6],Qiu等[5]認(rèn)為不同的潮汐強(qiáng)度會使各汊道分流比發(fā)生改變,進(jìn)一步影響鹽度的空間分布。北支的鹽水倒灌是長江口鹽水入侵的一個(gè)重要特征。因北支分流比低于5%而進(jìn)潮量占整個(gè)長江口的25%左右,因此北支是長江口鹽水入侵最嚴(yán)重的汊道。北支高濃度鹽水在徑流較小、潮動力較強(qiáng)時(shí)會倒灌進(jìn)入南支。丁磊等[7]對2013,2014年南北支分汊口處鹽度實(shí)測資料進(jìn)行分析,研究了北支鹽水倒灌的影響因素及下泄路徑。Wu等[8]研究了徑潮動力對北支倒灌鹽通量的影響。陳敏建等[9]研究了南支鹽度超標(biāo)面積與徑流、潮差的函數(shù)關(guān)系。孫昭華等[10]將潮差關(guān)系與農(nóng)歷日期建立聯(lián)系,提出了一種僅需知道大通流量就可快速估算北支鹽水倒灌影響下南支特定水域鹽度的方法。
東風(fēng)西沙水庫和陳行水庫取水口位于南支水域,青草沙水庫取水口位于北港水域,均會不同程度受到北支鹽水倒灌的影響。關(guān)于北支鹽水倒灌物理過程的表征,以往是通過基于鹽度場的等值線圖或是縱剖面鹽度圖進(jìn)行表示[11],該方法能夠從場的角度形象體現(xiàn)鹽水倒灌的平面過程,但無法刻畫鹽度隨時(shí)間的連續(xù)變化。而實(shí)際上測點(diǎn)鹽度過程線和潮位過程線的關(guān)系也可反映鹽水來源,此前鮮有這方面的詳細(xì)分析。同時(shí),倒灌鹽水團(tuán)對三大水源地影響時(shí)間的研究也較少。因此,本文通過建立數(shù)學(xué)模型對長江口水-鹽動力特性進(jìn)行模擬,對不同徑潮動力下長江口水源地受鹽水入侵的影響展開研究。研究結(jié)果可作為長江口淡水資源利用及水源地安全研究的依據(jù),為其他河口水源地建設(shè)和運(yùn)行調(diào)度提供參考。
長江口地區(qū)水資源總量為42.33億m3,其中地表水資源量38.00億m3。長江入海水量年內(nèi)分配不均勻,基本表現(xiàn)為洪季(5—10月)流量大,枯季(11月至次年4月)流量小。
根據(jù)大通站1950—2016年資料統(tǒng)計(jì),大通站多年平均流量28 300 m3/s,1954年8月1日出現(xiàn)最大流量92 600 m3/s,1979年1月31日出現(xiàn)最小流量為4 620 m3/s。三峽工程使得長江流量的年內(nèi)分布發(fā)生改變,枯季流量總體增加。2003年三峽水庫蓄水后各年最小流量如圖2所示。最小流量不會低于5 000 m3/s,尤其是到試驗(yàn)性蓄水運(yùn)行期(2008年汛后開始)時(shí),僅2014年出現(xiàn)低于10 000 m3/s的流量。
圖2 大通水文站年最小流量Fig.2 Annual minimum runoff of Datong station
長江口為中等潮差河口,中竣站年平均潮差2.66 m,最大潮差4.62 m。河口口門處總進(jìn)潮量為13億m3(枯季小潮)~53億m3(洪季大潮)。潮汐受外海潮波控制,口外潮汐為正規(guī)半日潮,口內(nèi)潮汐為非正規(guī)淺海半日潮。東海前進(jìn)波系統(tǒng)在本區(qū)域M2分潮為主,起支配作用;其次還受到黃海旋轉(zhuǎn)潮波影響,以K1 和O1分潮較顯著。潮波進(jìn)入長江口區(qū)域后,受邊界條件和上游徑流影響,潮波發(fā)生變形,既非典型的前進(jìn)波,也非典型的駐波。
目前能搜集到的長江口高精度鹽度資料主要來源于同步全潮水文測驗(yàn),優(yōu)點(diǎn)是精度高且為多點(diǎn)同步測量,缺點(diǎn)是成本高因而缺乏連續(xù)性。因此建立數(shù)學(xué)模型并用實(shí)測資料進(jìn)行驗(yàn)證是較為常用的研究手段。長江口枯季大部分區(qū)域鹽淡水混合類型為緩混合型,北支為強(qiáng)混合型,鹽度垂線差異較小,鹽水楔不明顯,因此采用平面二維數(shù)學(xué)模型可以對長江口的鹽水入侵進(jìn)行較好模擬,國內(nèi)學(xué)者有較為成功的經(jīng)驗(yàn)[3,12]。但因長江口人類活動明顯,局部地形變化較快,鹽水入侵會因此受到影響,采取不同地形會使研究結(jié)果產(chǎn)生差異。因此建立采用較新實(shí)測地形的數(shù)學(xué)模型和利用最近的實(shí)測資料對模型進(jìn)行驗(yàn)證較為必要。
利用Delft3D軟件建立長江口大范圍平面二維潮流鹽度數(shù)學(xué)模型對鹽度輸運(yùn)進(jìn)行模擬。De1ft3D是由荷蘭Delft水力研究院開發(fā)的,是目前較為先進(jìn)的水動力、水質(zhì)、泥沙等模型系統(tǒng)之一。Delft3D由6個(gè)模塊組成,各模塊既獨(dú)立又相互聯(lián)系,能較精確地進(jìn)行大尺度水流(Flow)、水動力(Hydro-dynamics)、波浪(Waves)、泥沙(Morphology)、水質(zhì)(Waq)和生態(tài)(Eco)計(jì)算。Delft3D模型的計(jì)算穩(wěn)定性強(qiáng),采用干濕動邊界處理技術(shù),對河口海岸區(qū)域有較好適應(yīng)性,可快速解決網(wǎng)格繪制、水深參數(shù)插值等問題,并具有強(qiáng)大的后處理功能。主要利用其中的Flow模塊對長江口鹽度輸運(yùn)過程進(jìn)行模擬研究。
水動力計(jì)算的淺水方程基于Navier-Stokes方程??刂品匠倘缦拢?/p>
連續(xù)性方程:
(1)
(2)
水平動量方程為:
(3)
(4)
式中:f為科里奧利參數(shù)(1/s);ρ0為水體密度(kg/m3);Pξ,Pη為ξ,η方向的靜水壓力梯度(kg/(m2·s2));Fξ,F(xiàn)η為ξ,η方向的紊動動量通量(m/s2);Mξ,Mη為ξ,η方向的動量源或匯(m/s2)。
在Delft 3D-FLOW模塊中,物質(zhì)輸運(yùn)采用對流擴(kuò)散方程進(jìn)行模擬。輸運(yùn)方程以守恒形式呈現(xiàn):
(5)
式中:c為鹽度;DH為水平擴(kuò)散系數(shù)(m/s2)。
模型采用正交曲線網(wǎng)格,且可以使網(wǎng)格線最大程度地貼合邊界線,避免“階梯”邊界導(dǎo)致發(fā)散。模型基于有限差分?jǐn)?shù)值方法。時(shí)間項(xiàng)的離散采用ADI差分格式,將1個(gè)時(shí)間步長分成2步,每一步為半個(gè)時(shí)間步長,前、后半個(gè)步長分別對不同方向進(jìn)行隱式處理。
計(jì)算范圍及網(wǎng)格如圖3所示,包括長江口、杭州灣及鄰近海域。北邊界位于北緯34.67°,南邊界最遠(yuǎn)位于北緯29.33°,東邊界最遠(yuǎn)位于東經(jīng)124.24°,東、南、北開邊界采用水位邊界,由主要分潮調(diào)和分析所得。上游邊界取在大通,為流量控制,驗(yàn)證時(shí)使用大通站逐時(shí)實(shí)測流量。橫向網(wǎng)格1 431個(gè),縱向163個(gè)。外海處網(wǎng)格尺寸最大達(dá)到2 km×2 km,長江口區(qū)域進(jìn)行了局部加密處理,最小網(wǎng)格尺寸為70 m×60 m。
江陰以下至長江口地形為2011年實(shí)測地形,江陰至大通為2006年實(shí)測地形。外海由海圖數(shù)字化得到。坐標(biāo)為高斯-克呂克坐標(biāo),高程統(tǒng)一為85高程系統(tǒng)。根據(jù)柯朗數(shù)(Courant number)原則,時(shí)間步長取15 s。模型糙率由謝才系數(shù)提供,根據(jù)地形情況采用不同的數(shù)值,范圍為80~200 m1/2/s。模型初始流速和初始水位設(shè)為0。外海邊界鹽度為35‰,長江口內(nèi)根據(jù)實(shí)測資料進(jìn)行插值,大通邊界鹽度為0,大通至徐六涇的鹽度由線性插值得到。擴(kuò)散系數(shù)為250 m2/s。模型計(jì)算運(yùn)行4個(gè)月(8個(gè)完整的半月周期)作為初始場進(jìn)行驗(yàn)證。
圖3 模型整體和局部網(wǎng)格Fig.3 Entire and local model grids
模型采用2016年1月長江口大范圍全潮同步水文測驗(yàn)進(jìn)行水動力與鹽度的驗(yàn)證。本次測驗(yàn)共布置10條垂線,基本覆蓋了長江口水域,同時(shí)還搜集了測驗(yàn)期間連興港、青龍港、崇頭、徐六涇、白茆、六滧(南)、高橋等潮位(圖4)。
圖4 水文測驗(yàn)測點(diǎn)及潮位站Fig.4 Hydrological test vertical lines and tidal stations
對大、中、小潮水位、流速、流向及鹽度進(jìn)行了驗(yàn)證。鑒于篇幅關(guān)系僅給出大潮時(shí)崇頭、六滧潮位驗(yàn)證及Z1,Y8流速、流向和鹽度驗(yàn)證結(jié)果(圖5)??梢钥闯觯撃P湍茌^好地模擬長江口水動力及鹽度輸運(yùn)。
圖5 模型驗(yàn)證結(jié)果Fig.5 Results of model verification
數(shù)學(xué)模型計(jì)算時(shí),上游采用恒定流。潮汐方面,取驗(yàn)證時(shí)的完整半月周期閉合循環(huán)計(jì)算。模型計(jì)算運(yùn)行4個(gè)月(8個(gè)完整的半月周期)作為初始場。分析時(shí)間為16 d,包含1個(gè)完整的半月周期,從最大潮差出現(xiàn)前4 d開始,再到最小潮差出現(xiàn)后4 d 結(jié)束。
長江口水源地取水口鹽度會隨漲落潮變化,若落潮時(shí)鹽度上升,漲潮時(shí)鹽度下降,該現(xiàn)象表明該處鹽水來源于取水口上游,即受北支倒灌鹽水影響;若漲潮時(shí)鹽度上升,落潮時(shí)鹽度下降,則表明該處鹽水來源于取水口下游,為外海鹽水正面影響或已經(jīng)下泄到取水口下游的鹽水團(tuán)隨漲潮流上溯所致。不同徑流下,水源地取水口鹽度大小不同,但鹽度過程變化特征較為相似。以大通徑流為12 500 m3/s(接近1月徑流50%頻率)為例[13-14],分析3個(gè)水源地取水口鹽度過程特征。對不同水源地而言,因受北支倒灌鹽水和正面入侵鹽水程度不同,鹽度過程也有著不同的特征。
3.2.1東風(fēng)西沙水庫 東風(fēng)西沙水庫取水口在半月周期中,鹽度過程與水位過程的關(guān)系根據(jù)潮差大小,表現(xiàn)出兩種不同的特征(圖6(a))。潮差較大時(shí),鹽度過程1 d 內(nèi)出現(xiàn)4個(gè)極大值。漲潮時(shí)鹽度上升,落潮時(shí)鹽度短暫下降后再次上升,且落潮時(shí)出現(xiàn)的鹽度極大值比漲潮時(shí)大得多(圖6(b))。隨著潮差減小,鹽度過程整體呈下降趨勢。潮差小到一定程度時(shí),落潮時(shí)的鹽度極大值消失,除了漲潮時(shí)鹽度有較小的上升外,鹽度總體呈下降趨勢(圖6(c)),即該階段北支不再出現(xiàn)鹽水倒灌。最小潮出現(xiàn)以后,落潮時(shí)鹽度極大值又逐漸出現(xiàn),表現(xiàn)為圖6(b)所示特征。對日平均鹽度進(jìn)行計(jì)算,發(fā)現(xiàn)半月周期中日均鹽度最大值出現(xiàn)在大潮當(dāng)天,較高鹽度持續(xù)到大潮后的中潮,日均鹽度上升時(shí)間略高于下降時(shí)間。
圖6 東風(fēng)西沙水庫取水口鹽度過程與水位過程Fig.6 Salinity and water level hydrograph in Dongfengxisha reservoir intake
3.2.2陳行水庫 陳行水庫取水口鹽度過程沒有出現(xiàn)單一的周期振動規(guī)律。在半月周期中,鹽度過程與水位過程的關(guān)系表現(xiàn)出了4種不同的特征(圖7(a))。大潮前的中潮后期及大潮期,鹽度在1 d 內(nèi)表現(xiàn)為兩漲兩落,漲潮時(shí)鹽度降低,落潮時(shí)鹽度升高。鹽度總體呈上升趨勢,且逐漸變緩(圖7(b))。此后,有半天鹽度特征發(fā)生改變,落潮時(shí)鹽度上升后,漲潮時(shí)并沒有立即下降,該階段鹽度變化較小且穩(wěn)定在一個(gè)較高的值,直到落潮的后期(落急以后)鹽度又明顯下降(圖7(c))。此后鹽度又恢復(fù)1 d 內(nèi)兩漲兩落,但仍然表現(xiàn)為漲潮時(shí)鹽度上升落潮時(shí)下降,鹽度總體呈下降趨勢(圖7(d))。該特征一直持續(xù)到小潮后的中潮,鹽度總體下降到一定程度后不再下降,且出現(xiàn)新的特征。落潮時(shí)再次出現(xiàn)鹽度升高,且越來越明顯,1 d 內(nèi)出現(xiàn)4個(gè)鹽度峰值(圖7(e))。日均鹽度最大值出現(xiàn)在大潮后的中潮,較高鹽度持續(xù)到小潮,且從最小值上升到最大值持續(xù)時(shí)間與最大值下降到最小值基本相同。
圖7 陳行水庫取水口鹽度過程與水位過程Fig.7 Salinity and water level hydrograph in Chenhang reservoir intake
3.2.3青草沙水庫 青草沙水庫取水口鹽度過程表現(xiàn)出4種不同特征(圖8)。如圖8(a)所示,小潮后期至大潮后期,漲潮時(shí)鹽度上升,落潮時(shí)鹽度下降,日鹽度變化幅度在四種特征中最大。大潮后期及大潮后的中潮鹽度日變化不再為兩漲兩落,1 d 內(nèi)出現(xiàn)4個(gè)鹽度峰值,除漲潮時(shí)鹽度上升外,落潮時(shí)鹽度在短暫下降后也會有上升。鹽度總體趨勢為波動中上升(圖8(b)),振幅相對前一階段減小。這種現(xiàn)象持續(xù)2 d后再次發(fā)生改變,鹽度過程再次表現(xiàn)為1 d 內(nèi)兩漲兩落,但與第一階段不同的是,漲潮時(shí)鹽度下降而落潮時(shí)鹽度上升。漲潮時(shí)的鹽度峰值消失而第二階段中出現(xiàn)的落潮峰值繼續(xù)出現(xiàn)。鹽度總體趨勢仍為波動中上升(圖8(c))。該階段持續(xù)3 d左右又發(fā)生改變,漲潮時(shí)鹽度不再下降,而是保持在一個(gè)相對固定的值,該現(xiàn)象持續(xù)1 d左右,在落潮時(shí)鹽度再次下降(圖8(d)),此時(shí)為小潮,鹽度過程再次表現(xiàn)為第一階段的特征。日均鹽度最大值出現(xiàn)在小潮,較高鹽度持續(xù)到小潮后的中潮,表現(xiàn)出了明顯的鹽度下降時(shí)間長于上升時(shí)間。
圖8 青草沙水庫取水口鹽度過程與水位過程Fig.8 Salinity and water level hydrograph in Qingcaosha reservoir intake
大通流量高于20 000 m3/s時(shí)幾乎不會有鹽水倒灌現(xiàn)象。模型上游流量取8 000~20 000 m3/s時(shí),每500 m3/s設(shè)置1個(gè)工況進(jìn)行計(jì)算。表1給出了不同大通流量條件下,半月周期中崇頭、東風(fēng)西沙水庫、陳行水庫以及青草沙水庫最大日平均鹽度出現(xiàn)時(shí)間。受正面入侵鹽水影響時(shí),表現(xiàn)為大潮鹽水入侵強(qiáng),小潮鹽水入侵弱。若南支以下水域半月周期中鹽度最大值出現(xiàn)在大潮后的其他時(shí)段,則為北支鹽水倒灌的影響[6]。水源地取水口日平均鹽度最大值出現(xiàn)時(shí)間與崇頭日平均鹽度最大值出現(xiàn)時(shí)間的差值可以表示北支倒灌最高濃度鹽水團(tuán)運(yùn)動到水源地取水口所需要的時(shí)間。
表1 日平均鹽度最大值出現(xiàn)時(shí)間Tab.1 Occurrence time for maximum value of daily average salinity
流量小于14 500 m3/s時(shí),崇頭日平均鹽度最大出現(xiàn)在分析的第5 d,即最大潮差出現(xiàn)當(dāng)天。流量大于14 500 m3/s時(shí),崇頭最大鹽度出現(xiàn)在分析的第6 d,即最大潮出現(xiàn)的后1天。雖然北支分流比小,但隨著流量增大會使得進(jìn)入北支的徑流量增大,對鹽水上溯起到了頂托作用,使得最大鹽度出現(xiàn)時(shí)間延遲。
東風(fēng)西沙水庫附近水域鹽水來源主要為北支的倒灌,最大鹽度出現(xiàn)時(shí)間均比崇頭滯后1 d。說明倒灌鹽水團(tuán)從崇頭運(yùn)動到東風(fēng)西沙水庫取水口需要1 d時(shí)間,取水口在大潮時(shí)受北支倒灌鹽水團(tuán)影響最嚴(yán)重。
陳行水庫日平均鹽度最大值出現(xiàn)時(shí)間比東風(fēng)西沙水庫滯后,但隨著流量增大,滯后時(shí)間逐漸變短。流量Q<10 000 m3/s時(shí)滯后3 d; 10 000 m3/s≤Q<14 500 m3/s時(shí)滯后2 d;Q≥14 500 m3/s時(shí)滯后1 d。倒灌鹽水團(tuán)從崇頭運(yùn)動到陳行水庫取水口需要2~4 d,取水口在大潮后的中潮時(shí)受北支倒灌鹽水團(tuán)影響最為嚴(yán)重。
青草沙水庫日平均鹽度最大值出現(xiàn)時(shí)間比東風(fēng)西沙水庫滯后,且長于陳行水庫的滯后時(shí)間。隨著流量的增大,滯后時(shí)間也逐漸變短,對徑流的響應(yīng)比陳行水庫更為敏感。當(dāng)大通流量Q<8 500 m3/s時(shí),青草沙水庫取水口最大鹽度出現(xiàn)時(shí)間比東風(fēng)西沙水庫取水口滯后8 d; 8 500 m3/s≤Q<10 000 m3/s時(shí)滯后7 d; 10 000 m3/s≤Q<12 000 m3/s時(shí)滯后6 d; 12 000 m3/s≤Q<14 500 m3/s時(shí)滯后5 d;Q≥14 500 m3/s時(shí)滯后4 d。倒灌鹽水團(tuán)從崇頭運(yùn)動到青草沙水庫取水口需要5~9 d,取水口在小潮時(shí)受北支倒灌鹽水團(tuán)影響最為嚴(yán)重。
徑流量的增大使得北支倒灌水體進(jìn)入南支后下泄速度加快,因此陳行水庫與青草沙水庫附近水域均表現(xiàn)出隨大通流量的增大,受倒灌鹽水團(tuán)影響時(shí)間提前的特點(diǎn)。圖9為水源地取水口在半月周期中日均鹽度的最大值與流量的關(guān)系。如圖所示,徑流越大,鹽水入侵程度越弱,受倒灌鹽水團(tuán)影響的程度也越弱。以往的研究指出流量與鹽度呈指數(shù)關(guān)系,而本文進(jìn)一步研究時(shí)發(fā)現(xiàn),當(dāng)流量在一定范圍內(nèi)時(shí),鹽度與徑流表現(xiàn)出很好的線性關(guān)系:
東風(fēng)西沙水庫:y=-0.000 4x+7.5(8 000≤x≤18 000)
陳行水庫:y=-0.000 4x+6.0(8 000≤x≤14 000)
青草沙水庫:y=-0.000 3x+4.9(8 000≤x≤16 000)
基于Delft3D建立長江口平面二維潮流鹽度數(shù)學(xué)模型。利用2016年長江口大范圍同步全潮水文測驗(yàn)資料對測驗(yàn)期間大、中、小潮的流速、流向、水位、鹽度進(jìn)行驗(yàn)證。驗(yàn)證效果良好,可以對長江口的水-鹽動力特性進(jìn)行模擬,對徑潮動力影響下長江口水源地取水口鹽度進(jìn)行研究,得出如下結(jié)論:
長江口枯季鹽水入侵會威脅水源地安全。北支鹽水倒灌是南支、北港水源地安全威脅的重要來源。三大水源地因位置不同,鹽度過程線特征也不相同。若落潮時(shí)鹽度上升,漲潮時(shí)鹽度下降,表明該處鹽水來源于取水口上游,即受北支倒灌鹽水影響;若漲潮時(shí)鹽度上升,落潮時(shí)鹽度下降,則表明該處鹽水來源于取水口下游,為外海鹽水正面影響或已經(jīng)下泄到取水口下游的鹽水團(tuán)隨漲潮流上溯所致。水源地取水口鹽度過程線與潮位過程線的關(guān)系可作為受北支倒灌鹽水和正面入侵鹽水影響程度的重要依據(jù)。
從南北支分汊口處崇頭開始計(jì)算,北支倒灌鹽水團(tuán)影響到東風(fēng)西沙水庫、陳行水庫、青草沙水庫取水口分別需要1, 2~4和5~9 d。水源地分別在大潮、大潮后的中潮以及小潮時(shí)受倒灌鹽水影響最為嚴(yán)重。隨著徑流增大,水源地受鹽水入侵影響的時(shí)間會提前,但是鹽度隨徑流增大而直線下降。