錢塘江河口抗咸流量預(yù)報(bào)模式及其檢驗(yàn)

李若華，史英標(biāo)，張舒羽

（1.河海大學(xué)環(huán)境學(xué)院，江蘇南京 210098；2.浙江省水利河口研究院，浙江杭州 310020）

錢塘江河口抗咸流量預(yù)報(bào)模式及其檢驗(yàn)

李若華1，2，史英標(biāo)2，張舒羽2

（1.河海大學(xué)環(huán)境學(xué)院，江蘇南京 210098；2.浙江省水利河口研究院，浙江杭州 310020）

錢塘江河口咸水入侵的影響因素有徑流、潮汐、江道地形、取水量等，其中徑流和潮汐的相互作用是取水口鹽度是否超標(biāo)的決定因素?；谙趟肭值闹饕绊懸蛩?，分析抗咸流量預(yù)報(bào)的關(guān)鍵點(diǎn)，建立了取水口鹽度超標(biāo)時(shí)間與徑流、潮差的定量響應(yīng)關(guān)系，構(gòu)建了抗咸流量預(yù)報(bào)模式，給出了模式中各參數(shù)的計(jì)算方法。將該預(yù)報(bào)模式應(yīng)用于錢塘江河口2012—2014年的最小抗咸流量模擬預(yù)報(bào)，經(jīng)實(shí)踐檢驗(yàn)，按照預(yù)報(bào)值調(diào)度下泄流量，可滿足杭州市抗咸要求?？瓜虒?shí)踐表明，不僅大潮期需保證足夠的抗咸流量，中、小潮期泄放流量也不得小于預(yù)報(bào)流量。

抗咸流量；咸水入侵；徑流；預(yù)報(bào)模式；錢塘江河口；潮汐

咸水入侵是潮汐河口普遍存在的一種自然現(xiàn)象，我國(guó)許多潮汐河口如長(zhǎng)江口［1］、珠江口［2］、錢塘江口［3］等的周邊城鎮(zhèn)供水經(jīng)常受咸潮災(zāi)害的威脅，為保證供水，上游水庫(kù)的壓咸調(diào)度是河口地區(qū)抵御咸水入侵以保證供水的主要措施之一。杭州市約85%的供水取自錢塘江河口段，每年的7—11月都要受到咸水入侵的威脅，主要通過(guò)新安江水庫(kù)和富春江水庫(kù)的聯(lián)合調(diào)度以增大下泄流量來(lái)抗咸，抗咸流量由從事多年抗咸工作的老專家根據(jù)以往抗咸經(jīng)驗(yàn)進(jìn)行預(yù)報(bào)，沒(méi)有固定的預(yù)報(bào)模式，這不利于此項(xiàng)預(yù)報(bào)技術(shù)的推廣，因此建立普適的抗咸流量預(yù)報(bào)模式非常必要。

錢塘江河口是強(qiáng)涌潮、高含沙量、河床變形劇烈的游蕩性河口，咸水入侵規(guī)律極其復(fù)雜。對(duì)于錢塘江河口的咸水入侵規(guī)律研究，20世紀(jì)80年代以來(lái)已開(kāi)展了許多工作，取得了大量的研究成果［4-7］。模擬咸水入侵的數(shù)學(xué)模型從短歷時(shí)的定床模型［8］和可考慮涌潮作用的定床模型［9］，發(fā)展到可進(jìn)行長(zhǎng)歷時(shí)鹽度預(yù)報(bào)的動(dòng)床模型［10］，并提出了“大潮多放、小潮少放”（相差±25%）［11］的徑流調(diào)度方式。然而，抗咸預(yù)報(bào)不是預(yù)報(bào)取水口處的鹽度，而是預(yù)報(bào)為保證取水口鹽度達(dá)標(biāo)所需的抗咸流量，是鹽度預(yù)報(bào)問(wèn)題的反問(wèn)題，需要根據(jù)潮汐及咸水入侵規(guī)律，基于當(dāng)年的江道形勢(shì)、咸水入侵程度、抗咸目標(biāo)等約束條件進(jìn)行預(yù)報(bào)。因預(yù)報(bào)涉及因素眾多，且多個(gè)因素是未知數(shù)，目前抗咸流量靠經(jīng)驗(yàn)預(yù)報(bào)，尚未建立固定化的預(yù)報(bào)模式，對(duì)該預(yù)報(bào)技術(shù)的推廣是非常不利的。本文從咸水入侵的影響因素著手，推導(dǎo)抗咸流量預(yù)報(bào)需要的參數(shù)，建立預(yù)報(bào)模式，給出模式中各參數(shù)的計(jì)算方法，并應(yīng)用于杭州抗咸的最小流量預(yù)報(bào)，結(jié)合抗咸實(shí)際效果檢驗(yàn)?zāi)Ｊ降倪m用性。

1　咸水入侵影響因素及抗咸期

1.1錢塘江河口概況

錢塘江河口是一個(gè)典型的強(qiáng)潮河口，浦陽(yáng)江口—澉浦為河口段（圖1），受徑流和潮汐共同作用。錢塘江河口多年平均徑流流量約952 m3/s，洪汛期3—6月的徑流量占全年的60%以上。徑流不僅年內(nèi)變化大，年際變化也很大，年最大徑流量可達(dá)年最小徑流量的4倍以上。潮汐受M2分潮控制，一天兩漲兩落，澉浦站最大潮差可達(dá)9 m以上，澉浦以上潮差逐漸減小。錢塘江河口段水體含沙量大，泥沙以粉砂為主，河床易沖易淤。因潮強(qiáng)流急，河口受咸水入侵影響嚴(yán)重，枯水大潮期咸水可隨潮上溯至浦陽(yáng)江口，其中閘口—倉(cāng)前段咸、淡水摻混最為劇烈，而杭州市取水口就布置在閘口河段。

圖1　錢塘江河口段示意圖

1.2影響因素分析

咸水入侵受諸多因素的影響，包括徑流、潮汐、江道地形和沿程取水量等，其中徑流和潮汐是決定性的因素，取水口鹽度是否超標(biāo)主要取決于短時(shí)間內(nèi)徑流和潮汐的相互作用。錢塘江河口咸水入侵的影響因素分述如下：

a.徑流。徑流對(duì)咸水入侵的影響反映在不同時(shí)間尺度的變化上。從年尺度講，年均徑流量大，則江道容積大，潮汐動(dòng)力強(qiáng)，為保證取水口鹽度達(dá)標(biāo)則需要更大的徑流流量；從季節(jié)尺度講，徑流分兩個(gè)層面，洪汛前及洪汛期的徑流總量決定汛后秋季大潮汛期的江道地形，進(jìn)而決定下半年河口段的潮汐強(qiáng)弱及咸水入侵程度，大潮汛期的徑流流量則決定咸水上溯的距離及取水口鹽度超標(biāo)的情況以及江道淤積量，進(jìn)而影響潮汐強(qiáng)度。同時(shí)河口段是一個(gè)敏感的小系統(tǒng)，初始鹽度的影響可在數(shù)天內(nèi)消除，根據(jù)以往抗咸經(jīng)驗(yàn)，3日徑流流量的調(diào)整即會(huì)引起咸水入侵的明顯差異，這也是保障錢塘江河口御咸調(diào)度及時(shí)性的自然基礎(chǔ)。

b.潮汐。潮汐對(duì)咸水入侵的影響主要反映在季節(jié)間潮汐大小的變化，以及一個(gè)潮汛期（15 d）內(nèi)大小潮的差異等，漲潮前期的強(qiáng)涌潮可促使咸水入侵加劇，鹽分快速升高。此外，即時(shí)徑流和江道地形又決定了取水口河段的潮汐強(qiáng)度及咸水入侵程度，取水口河段的潮汐強(qiáng)度可由下游七堡站的逐月最大潮差代表，該站潮差對(duì)該河段咸水入侵程度起控制作用。

c.江道地形。江道地形對(duì)河口咸水入侵的影響主要反映在江道沖淤及尖山河灣走直、走彎河勢(shì)對(duì)咸水入侵的影響，一般而言，江道容積大，河勢(shì)走直，潮汐上溯能力強(qiáng)，則有利于咸水入侵，反之則不利于咸水入侵。江道地形的變化對(duì)咸水入侵的影響可同步反映在潮差大小上，即江道地形的變化帶來(lái)的影響可由七堡站潮差代表。

d.取水流量。在下半年的枯水大潮期，上游來(lái)水較少，而此時(shí)高溫干旱，兩岸農(nóng)業(yè)、環(huán)境用水等需求量大，占徑流量的比重較高，因此取水流量亦會(huì)對(duì)咸水入侵產(chǎn)生明顯影響。根據(jù)浙江省水利廳制定的《錢塘江河口水資源配置工程調(diào)度辦法（試行）》，抗咸期間富春江流量在450 m3/s以下時(shí)，沿線允許取水總量不超過(guò)50 m3/s，該調(diào)度辦法限制了抗咸期間的無(wú)序取水，從而使得取水量對(duì)咸水入侵的影響減弱。

e.前期鹽度。前期鹽度指大潮之前累積在取水口下游的咸水，若咸水鹽分較高，在大潮期來(lái)臨時(shí)將隨著潮流強(qiáng)勁上溯，造成取水口鹽度超標(biāo)，但若中小潮期保證了足夠的抗咸流量，則可以將咸水輸送至下游，避免鹽度在河口段累積。

1.3抗咸期

通常情況下，錢塘江河口1—3月河床淤積，且外海天文潮小，潮流上溯能力弱；4—6月降水量增多并進(jìn)入梅雨季，徑流量大，從而壓制了咸水入侵，因此上半年河口段沿程鹽度均較小，咸潮一般不會(huì)影響取水口的正常取水。在7—11月這段時(shí)間，一方面雨量驟減，又適逢夏秋之際的高溫期，農(nóng)業(yè)灌溉用水、環(huán)境用水和生活用水量猛增，徑流量較小，因而壓制咸潮的能力減弱；而另一方面又恰遇秋季大潮汛，外海鹽度較高的水體隨潮上溯，因此這段時(shí)間是取淡水最不利的時(shí)期。11—12月，雖然降雨量不多，但灌溉用水、環(huán)境用水和居民生活用水明顯減少，更主要的是徑流的減小導(dǎo)致下游河道淤積，潮汐上溯強(qiáng)度減弱，咸水入侵程度減輕。因此抗咸期主要為7—11月。

枯水年洪汛期水量小，雖然新安江水庫(kù)在洪汛期儲(chǔ)存的用于下半年抗咸的水量少，但汛后江道容積小，下半年咸水入侵力度弱，抗咸形勢(shì)反而比較樂(lè)觀，豐水年洪汛期的大水量沖開(kāi)了江道，導(dǎo)致下半年潮汐上溯能力大幅加強(qiáng)，咸水入侵加劇，抗咸形勢(shì)比較嚴(yán)峻。因此，豐水年的7—11月是抗咸的重點(diǎn)時(shí)期。

2　抗咸流量預(yù)報(bào)模式

2.1抗咸流量預(yù)報(bào)模式示意圖

據(jù)GB5749—2006《生活飲用水衛(wèi)生標(biāo)準(zhǔn)》，當(dāng)江水含氯度超過(guò)250 mg/L時(shí)定為超標(biāo)。由前述可知，取水口鹽度是否超標(biāo)主要取決于徑流和潮汐的相互作用，為達(dá)到抗咸目標(biāo)，可根據(jù)潮汐強(qiáng)度（可由七堡站潮差代表）下泄相應(yīng)的徑流量壓潮拒咸，因此需建立取水口鹽度超標(biāo)時(shí)間（即抗咸目標(biāo)）與徑流、潮汐之間的定量響應(yīng)關(guān)系，同時(shí)為確定抗咸流量需預(yù)測(cè)出下半年的潮汐強(qiáng)度；因外海潮汐年際變化較小，取水河段的潮汐強(qiáng)弱年際變化主要由江道地形決定，梅汛期后7月的江道地形對(duì)下半年咸水入侵程度起到至關(guān)重要的作用，因此為得到下半年的潮汐強(qiáng)度需預(yù)測(cè)出梅汛期后的江道地形；錢塘江河口每年汛前4月進(jìn)行一次大范圍的水下地形測(cè)驗(yàn)，因此汛后江道形勢(shì)可根據(jù)汛前水下地形和汛期徑流量通過(guò)一定的計(jì)算方法估算。根據(jù)上述分析，可設(shè)計(jì)抗咸流量預(yù)報(bào)模式如圖2所示。

圖2　抗咸流量預(yù)報(bào)模式示意圖

2.2抗咸目標(biāo)

受咸水入侵影響的保護(hù)目標(biāo)主要為沿線的取水口。杭州市九溪、南星橋、清泰、赤山埠4個(gè)取水口位于錢塘江閘口河段，與之相應(yīng)的水廠的總制水量約140萬(wàn)m3/d，占杭州市總制水量的85%。珊瑚沙水庫(kù)位于閘口上游，為避咸蓄淡調(diào)節(jié)水庫(kù)，總庫(kù)容近190萬(wàn)m3，當(dāng)外江鹽度超標(biāo)時(shí)，可經(jīng)輸水渠道直接供給相應(yīng)的水廠。綜合珊瑚沙水庫(kù)的庫(kù)容及各水廠的蓄水量，杭州市水業(yè)集團(tuán)根據(jù)供水要求，提出抗咸目標(biāo)為：①南星橋取水口鹽度連續(xù)超標(biāo)時(shí)間不大于2 d；②珊瑚沙水庫(kù)取水口鹽度連續(xù)超標(biāo)不超過(guò)0.5 d。該抗咸目標(biāo)可根據(jù)蓄水能力和供水要求改變。

2.3抗咸目標(biāo)與徑流、潮汐的關(guān)系

借助考慮涌潮作用下的平面二維鹽度數(shù)學(xué)模型［9］，進(jìn)行多組徑流和潮汐的組合計(jì)算，得到南星橋、珊瑚沙取水口鹽度超標(biāo)時(shí)間與凈徑流、潮汐的定量響應(yīng)關(guān)系曲線如圖3所示（圖中潮差為七堡站潮差）。由圖3可見(jiàn)，當(dāng)七堡站潮差為2.8 m時(shí)，大潮期徑流流量需大于320 m3/s，疊加沿線總?cè)∷?0 m3/s后，得到最小抗咸流量為370 m3/s。因此，若已知七堡站最大潮差，則可根據(jù)圖3預(yù)報(bào)大潮期的最小凈抗咸流量，然后再根據(jù)“大潮多放25%流量、小潮少放25%流量”的調(diào)度方式，疊加沿線取水量，預(yù)報(bào)出中、小潮期的最小抗咸流量。

圖3　南星橋、珊瑚沙取水口鹽度超標(biāo)時(shí)間對(duì)徑流流量、潮差的響應(yīng)

為提升關(guān)系曲線的應(yīng)用價(jià)值，對(duì)圖3曲線通過(guò)回歸分析進(jìn)行擬合，構(gòu)建相關(guān)函數(shù)關(guān)系。當(dāng)七堡站潮差為3.9 m時(shí)，南星橋、珊瑚沙取水口鹽度超標(biāo)時(shí)間對(duì)凈徑流、潮汐的響應(yīng)關(guān)系分別為

式中：y為取水口最長(zhǎng)連續(xù)超標(biāo)時(shí)間，d；x為徑流流量，m3/s。

2.4七堡站月最大潮差預(yù)測(cè)

據(jù)多年研究經(jīng)驗(yàn)，七堡站月最大潮差與閘口—鹽官河段吳淞高程7 m下的江道容積相關(guān)關(guān)系較好，將多年的江道容積與同期的七堡站月最大潮差數(shù)據(jù)點(diǎn)繪在圖上，連成曲線后如圖4所示［11］。因此，若已知汛后的江道容積，則可根據(jù)圖4預(yù)測(cè)出下半年七堡站的月最大潮差。

圖4　不同江道容積條件下七堡站月最大潮差預(yù)報(bào)曲線

2.5汛后江道容積預(yù)測(cè)

錢塘江鹽官以上河段沖淤規(guī)律是洪沖潮淤，即豐水期鹽官以上河段沖刷，沖起的泥沙被帶至下游河段淤積，閘口—鹽官河段的江道容積增大，潮汐強(qiáng)度增強(qiáng)；枯水期則相反，潮流將下游的泥沙帶至鹽官以上河段淤積，江道容積減小，潮汐強(qiáng)度減弱。一般情況下，錢塘江河口段每年都要經(jīng)歷淤積—沖刷—淤積的自然演變的過(guò)程。錢塘江河口河床在徑、潮流的作用下，沖淤幅度大、速度快，尤其是汛期的大流量造床作用非常強(qiáng)。決定下半年潮汐強(qiáng)度的汛后江道容積可根據(jù)下式［11］計(jì)算：

式中：Vb為吳淞高程7 m條件下汛后7月的江道容積；ˉQ為汛期4—7月的平均徑流流量；Vi為吳淞高程7 m條件下汛前4月的江道容積。據(jù)實(shí)測(cè)資料驗(yàn)證，該式相關(guān)系數(shù)可達(dá)0.93。

此外，也可以借助動(dòng)床數(shù)學(xué)模型［7］，在4月江道地形的基礎(chǔ)上，輸入汛期的逐日徑流流量，計(jì)算出汛后7月的江道地形，然后統(tǒng)計(jì)出閘口—鹽官河段吳淞高程7 m下的江道容積。因泥沙沖淤的復(fù)雜性，目前動(dòng)床模型計(jì)算結(jié)果的精度并不比經(jīng)驗(yàn)公式算出的結(jié)果高，但工作量卻大幅增加。

3　預(yù)報(bào)模式的檢驗(yàn)

3.12012年抗咸過(guò)程模擬預(yù)報(bào)

錢塘江河口徑流基本以10年為周期進(jìn)行豐、枯水文年的轉(zhuǎn)換［12］。2009年之前的10年，以枯水年為主，汛后江道容積大多小于4億m3；2010年徑流轉(zhuǎn)豐后，江道容積大幅增加，2010年、2011年汛后江道容積分別達(dá)5.47億m3和5.23億m3，已經(jīng)屬于豐水江道。

2012年汛前江道容積為4.78億m3，屬豐水江道，汛期平均徑流流量為1980 m3/s，采用式（3）算出汛后江道容積為5.39億m3。根據(jù)預(yù)測(cè)的汛后江道容積，對(duì)照?qǐng)D4插值可得2012年7—12月的七堡站月最大潮差（表1），其中8月潮差最大，為3.5 m。預(yù)測(cè)的2012年汛后江道容積比2011年（5.23億m3）略大，2011年下半年七堡站實(shí)測(cè)最大潮差為3.28 m，實(shí)測(cè)日徑流流量為886 m3/s。若徑流量較小，則潮差會(huì)更大，因此模擬預(yù)測(cè)的2012年下半年七堡站月最大潮差3.5 m是基本合理的。根據(jù)七堡站月最大潮差預(yù)報(bào)值，對(duì)照?qǐng)D3，結(jié)合不同潮差時(shí)取水口超標(biāo)時(shí)間對(duì)凈徑流的響應(yīng)關(guān)系，得到各月大潮期最小抗咸徑流流量模擬預(yù)報(bào)值（表1）。

表1　2012年7—12月七堡站月最大潮差與最小抗咸徑流流量

根據(jù)2012年杭州灣天文潮預(yù)報(bào)，控制河口段短期潮汐規(guī)律的澉浦站下半年望（農(nóng)歷月半）潮略高于朔（農(nóng)歷月初）潮，如農(nóng)歷七月十六（9月1日）、七月初四（8月20日）澉浦月最高潮位預(yù)報(bào)值分別為6.76 m和6.67m。結(jié)合天文潮的月內(nèi)分布及各月大潮期最小抗咸流量，基于“大潮多放25%流量，小潮少放25%流量”的調(diào)度方式，對(duì)逐日的抗咸流量進(jìn)行分配，并統(tǒng)一疊加沿線總?cè)∷?0 m3/s，模擬預(yù)報(bào)出2012年下半年最小抗咸流量如圖5所示。由圖5可知，下半年（7—12月）抗咸共需要最少51.6億m3徑流量，其中8—9月抗咸徑流流量較大，大潮時(shí)最大為430 m3/s，小潮時(shí)最小為260 m3/s。

圖5　2012年下半年逐日最小抗咸徑流流量

根據(jù)相關(guān)水廠的記錄，2012年南星橋取水口、珊瑚沙取水口最長(zhǎng)連續(xù)超標(biāo)時(shí)間分別為19.25 h、7.50 h，勉強(qiáng)滿足抗咸要求，大潮期間預(yù)報(bào)與實(shí)際下泄流量對(duì)比及抗咸情況列于表2。由表2可見(jiàn)，鹽度超標(biāo)時(shí)間主要出現(xiàn)在10月中旬，大潮期實(shí)際下泄流量略大于模擬預(yù)報(bào)流量，但15 d平均流量比預(yù)報(bào)流量小10%，中、小潮期實(shí)際流量比預(yù)報(bào)流量平均小96 m3/s，平均相差幅度達(dá)30.8%（表3）。大潮期過(guò)后，河口段累積的鹽分在徑流的作用下逐步移到下游，咸水在小潮期被推送到最遠(yuǎn)，小潮期過(guò)后的中潮期鹽分又逐漸上移累積，因此中、小潮期也必須保證足夠的抗咸流量，以免造成鹽分累積而導(dǎo)致下一個(gè)大潮期鹽度超標(biāo)影響取水。由此可見(jiàn)，2012年最小抗咸流量模擬預(yù)報(bào)成果是比較準(zhǔn)確的，只要按照預(yù)報(bào)調(diào)度下泄流量，則可以保證抗咸要求，調(diào)度流量小于預(yù)報(bào)流量，即使是中小潮期，也會(huì)造成鹽度超標(biāo)。

表2　2012年大潮期間下泄流量及頂潮拒咸情況

表3　2012年10月中小潮期徑流實(shí)際流量與預(yù)報(bào)流量差異

3.22013—2014年抗咸模擬預(yù)報(bào)

2013年的抗咸實(shí)踐中，10月上旬也出現(xiàn)了較長(zhǎng)時(shí)間的鹽度超標(biāo)，南星橋取水口、珊瑚沙取水口最長(zhǎng)連續(xù)超標(biāo)時(shí)間分別為30.5 h、4 h，雖然在抗咸目標(biāo)要求范圍內(nèi)，但已經(jīng)造成了杭州水廠的緊張，其原因是在大潮期前的中、小潮期下泄的抗咸流量不足，實(shí)際流量比預(yù)報(bào)流量平均小164 m3/s，相差幅度為53.4%。由此可見(jiàn)，不僅在大潮期需要保證足夠的抗咸流量，中小潮期的抗咸流量也必須保證。

2014年抗咸期間，大、中、小潮期均保證了足夠的抗咸流量，南星橋、珊瑚沙取水口均沒(méi)有出現(xiàn)鹽度超標(biāo)的情況。

4　結(jié) 論

a.影響錢塘江河口段取水口鹽度的主要因素有徑流、潮汐、江道地形、取水流量、前期鹽度等，徑流減小、潮汐增強(qiáng)、江道容積增加、取水流量加大、大潮前期鹽度較高都將加劇咸水入侵，其中影響咸水入侵最重要的因素是徑流和潮汐的相互作用。

b.基于咸水入侵的影響因素，分析了抗咸流量預(yù)報(bào)的關(guān)鍵，構(gòu)建了抗咸預(yù)報(bào)模式，并采用數(shù)學(xué)模型建立了取水口鹽度超標(biāo)時(shí)間與徑流、潮差的定量響應(yīng)關(guān)系，在整合前人研究成果的基礎(chǔ)上，給出了預(yù)報(bào)模式中各參數(shù)的計(jì)算方法。

c.應(yīng)用預(yù)報(bào)模式對(duì)2012—2014年的抗咸流量進(jìn)行了模擬預(yù)報(bào)，結(jié)果表明，按照預(yù)報(bào)流量調(diào)度下泄徑流，可滿足杭州市抗咸要求?？瓜虒?shí)踐表明，不僅需在大潮期保證足夠的抗咸流量，中、小潮期也必須保證足夠的泄放流量。

［1］吳繼偉，劉新成，潘麗紅.長(zhǎng)江口北支咸潮倒灌控制工程水動(dòng)力數(shù)值模擬［J］.水利水電科技進(jìn)展，2006，26（4）：43-45.（WU Jiwei，LIU Xincheng，PAN Lihong. Hydrodynamic simulation of control project on saltwater intrusion from north branch of Yangtze River Estuary［J］. Advances in Science and Technology of Water Resources，2006，26（4）：43-45.（in Chinese））

［2］方神光，崔麗琴.磨刀門水道枯季咸潮入侵特性及規(guī)律［J］.水利水電科技進(jìn)展，2015，35（4）：14-18.（FANG Shenguang，CUILiqin.Characteristicsandrulesof seawaterintrusionduringdryseasoninModaomen Waterway［J］.Advances in Science and Technology of Water Resources，2015，35（4）：14-18.（in Chinese））

［3］尤愛(ài)菊，韓曾萃，史英標(biāo)，等.浙北引水對(duì)錢塘江河口段主要水環(huán)境的影響分析［J］.河海大學(xué)學(xué)報(bào)（自然科學(xué)版），2006，34（2）：152-156.（YOU Aiju，HAN Cengcui，SHI Yingbiao，et al.Impact of water diversion from Xinanjiang Reservoir to North Zhejiang Province on water environment at Qiantang Estuary［J］.Journal of Hohai University（Natural Sciences），2006，34（2）：152-156.（in Chinese））

［4］HAN Zengcui，SHAO Yaqin.Salt water intrusion and counter measures in some coastal cities in China［J］. China Ocean Engineering，1986，3（2）：177-193.

［5］韓曾萃，潘存鴻，史英標(biāo)，等.人類活動(dòng)對(duì)咸水入侵的影響［J］.水科學(xué)進(jìn)展，2002，13（3）：333-336.（HAN Zengcui，PAN Cunhong，SHI Yingbiao，et al.Impact of human activities on salt water intrusion in estuaries［J］. Advances in Water Science，2002，13（3）：333-336.（in Chinese））

［6］趙渭軍，韓繼靜，史英標(biāo)，等.杭州市抗咸二期取水工程氯度研究［J］.中國(guó)農(nóng)村水利水電，2007（5）：26-29.（ZHAO Weijun，HAN Jijing，SHI Yingbiao，et al.Study on the chlorinity of the second-stage salinity-resisting and water-intake project in Hangzhou City［J］.China Rural Water and Hydropower，2007（5）：26-29.（in Chinese））

［7］史英標(biāo)，潘存鴻，程文龍，等.錢塘江河口段鹽度時(shí)空分布與預(yù)測(cè)模型［J］.水科學(xué)進(jìn)展，2012，23（3）：409-418.（SHI Yingbiao，PAN Cunhong，CHENG Wenlong，et al.Temporal-spatial variation and numerical forecast model of salt water intrusion of the Qiantang Estuarine Reach［J］.Advances in Water Science，2012，23（3）：419-428.（in Chinese））

［8］韓曾萃，程杭平.錢塘江江水含鹽度計(jì)算的研究［J］.水利學(xué)報(bào)，1981，12（6）：46-50.（HAN Zengcui，CHENG Hangping.Study on the salinity calculation of the Qiantang Estuary［J］.Journal of Hydraulics Engineering，1981，12（6）：46-50.（in Chinese））

［9］潘存鴻，張舒羽，史英標(biāo)，等.涌潮對(duì)錢塘江河口鹽水入侵影響研究［J］.水利學(xué)報(bào)，2014，45（11）：1301-1309.（PAN Cunhong，ZHANG Shuyu，SHI Yinbiao，et al.Study on saltwater intrusion in Qiantang Estuary affected by tidal bore［J］.Journal of Hydraulic Engineering，2014，45（11）：1301-1309.（in Chinese））

［10］史英標(biāo)，李若華，姚凱華.錢塘江河口一維鹽度動(dòng)床預(yù)報(bào)模型及應(yīng)用［J］.水科學(xué)進(jìn)展，2015，26（2）：212-220.（SHI Yingbiao，LI Ruohua，YAO Kaihua.1-D movable bed umerical forecast model of salinity of the Qiantang River Estuarine Reach and its application［J］. Advances in Water Science，2015，26（2）：212-220.（in Chinese））

［11］韓曾萃，程杭平，史英標(biāo)，等.錢塘江河口咸水入侵長(zhǎng)歷時(shí)預(yù)測(cè)和對(duì)策［J］.水利學(xué)報(bào)，2012，43（2）：232-240.（HAN Zengcui，CHENG Hangping，SHI Yingbiao，et al. Long-term predictions and countermeasures of saltwater intrusion in the Qiantang Estuary［J］.Journal of Hydraulic Engineering，2012，43（2）：232-240.（in Chinese））

［12］曾劍，孫志林，潘存鴻，等.錢塘江河口徑流長(zhǎng)周期特性及其對(duì)河床的影響［J］.浙江大學(xué)學(xué)報(bào)（工學(xué)版），2010，44（8）：1584-1588.（ZENG Jian，SUN Zhilin，PAN Cunhong，et al.Long periodic feature runoff and its effect on riverbed in Qiantang Estuary［J］.Joural of Zhejiang University（Engineering Science），2010，44（8）：1584-1588.（in Chinese））

Prediction model of salinity-resisting discharge in Qiantang Estuary and its validation

LI Ruohua1，2，SHI Yingbiao2，ZHANG Shuyu2
（1.College of Environment，Hohai University，Nanjing 210098，China；2.Zhejiang Institute of Hydraulics&Estuary，Hangzhou 310020，China）

Runoff，tides，riverbed morphology，and water intake influence saltwater intrusion in the Qiantang Estuary.Of these factors，runoff and tides are the decisive factors that determine whether the salinity at the intake exceeds the standard. Based on the main factors in saltwater intrusion，key issues in predicting the salinity-resisting discharge were analyzed，and a quantitative relationship regarding the response of the standard-exceeding time to the runoff and tidal range was established.A model for predicting the salinity-resisting discharge was built，and computational methods for parameters in the model were proposed.The prediction model was applied in the prediction of minimum salinity-resisting discharge in the Qiantang Estuary from 2012 to 2014.The validation results show that reservoir operation based on the predicted discharge can satisfy the salinity-resisting demand of Hangzhou City，and a discharge sufficient to resist saltwater intrusion must be guaranteed not only during spring tides but also during mid-and neap tides.

salinity-resisting discharge；saltwater intrusion；runoff；prediction model；Qiantang Estuary；tide

X522

A

10067647（2016）05007506

10.3880/j.issn.10067647.2016.05.014

浙江省科技計(jì)劃（2014F10036，2015F50010）

李若華（1977—），男，高級(jí)工程師，博士，主要從事河口環(huán)境水力學(xué)研究。E-mail：liruohua2000@163.com

（20150907 編輯：鄭孝宇）