金德鋼,孫 堯,夏珊珊

(寧波市水利水電規(guī)劃設(shè)計(jì)研究院,浙江 寧波 315192)

1 問題的提出

寧波市地處我國(guó)東南沿海,沿海和島嶼地區(qū)普遍缺水,加之山低源短,無建造大中型水庫(kù)的條件。但海涂資源豐富,根據(jù)當(dāng)?shù)貙?shí)際情況建設(shè)了許多堵港蓄淡水庫(kù)和海涂水庫(kù),用以攔蓄陸域來水或豐水期抽蓄河網(wǎng)淡水,其中規(guī)模較大的有象山縣大塘港、寧?？h胡陳港、毛嶼港、一市港、車岙港下水庫(kù)及慈溪四灶浦等蓄淡水庫(kù)。正在興建鄭徐海涂水庫(kù)和規(guī)劃興建慈西、臨海浦等海涂水庫(kù)。在短時(shí)間解決堵港水庫(kù)和海涂水庫(kù)咸水淡化問題,已提到議事日程。本文結(jié)合二維水動(dòng)力模型對(duì)慈溪鄭徐海涂水庫(kù)排咸設(shè)計(jì)成果分析,提出對(duì)海涂水庫(kù)排咸設(shè)計(jì)理論計(jì)算方法和方案措施探討。

2 工程概況

鄭徐海涂水庫(kù)工程位于慈溪市范市鎮(zhèn)、觀海衛(wèi)鎮(zhèn)、附海、勝山、崇壽及新浦鎮(zhèn)的灘涂上,建于鄭家浦與徐家浦圍涂區(qū)內(nèi),南靠九塘,北至十塘,東靠徐家浦河道西直堤,西至鄭家浦隔提,占地666.7 hm2(1萬畝),總庫(kù)容4 270萬m3[1]。工程由堤壩、翻水泵站、放水閘、排咸站、配套河道、節(jié)制閘組成,工程于2010年11月開工,計(jì)劃2013年完工。

3 海涂水庫(kù)水質(zhì)咸化主要影響因素分析

(1)庫(kù)底鹽堿底泥:庫(kù)底沉積物中鹽分釋放是一個(gè)緩慢過程,其對(duì)海涂水庫(kù)水質(zhì)影響主要在建庫(kù)初期[2-3]。

(2)咸化淺層地下水:地下水含鹽量高于水庫(kù)水體的含鹽量,是導(dǎo)致水庫(kù)水質(zhì)咸化的主要因素[2-3]。地下水與水庫(kù)水的補(bǔ)給關(guān)系與水庫(kù)蓄水水位的變化有很大關(guān)系,當(dāng)水庫(kù)蓄水水位較高時(shí),由水庫(kù)水補(bǔ)給地下水,水庫(kù)水質(zhì)受地下水的影響較小。反之,庫(kù)底地下水將會(huì)通過土壤孔隙與上覆庫(kù)內(nèi)水體相連,地下水將補(bǔ)給水庫(kù)水。

(3)氣候氣象因素:蒸發(fā)是影響水質(zhì)咸化的一個(gè)重要因素,蒸發(fā)濃縮作用占優(yōu)勢(shì)也引起水庫(kù)水質(zhì)咸化。尤其在炎熱的夏季,地下水、水庫(kù)水體中的鹽分?jǐn)U散系數(shù)增大,可加快咸化的速度,經(jīng)常會(huì)導(dǎo)致突發(fā)性泛咸現(xiàn)象。風(fēng)是影響水質(zhì)咸化的另一個(gè)因素,風(fēng)吹過水面產(chǎn)生剪切應(yīng)力,將形成垂直環(huán)流和水平環(huán)流,使鹽分從庫(kù)底鹽堿土壤和咸化地下水向水庫(kù)水體中擴(kuò)散,從而加快水庫(kù)水質(zhì)的咸化速度。

4 研究思路及模型介紹

4.1 研究思路

通過建立庫(kù)區(qū)水動(dòng)力模型,研究水庫(kù)在可能運(yùn)行條件下的水動(dòng)力特性、鹽度分布規(guī)律,在此基礎(chǔ)上,根據(jù)鄭徐水庫(kù)實(shí)際,研究制定相應(yīng)的排咸工程措施和非工程措施,并用底泥鹽度釋放模型對(duì)水庫(kù)鹽度淡化趨勢(shì)進(jìn)行預(yù)測(cè),為水庫(kù)運(yùn)行管理提供支撐,實(shí)現(xiàn)淡化水體、提高水質(zhì)的目的,最大程度地發(fā)揮水庫(kù)效益。

4.2 數(shù)學(xué)模型介紹

4.2.1 Delft 3D水動(dòng)力模型

使用Delft 3D軟件建立鄭徐水庫(kù)庫(kù)區(qū)二維水動(dòng)力模型,按照不同工況分析庫(kù)區(qū)的水動(dòng)力流場(chǎng),在此基礎(chǔ)上進(jìn)一步研究引水調(diào)度對(duì)水庫(kù)鹽度變化的影響,為鄭徐水庫(kù)運(yùn)行調(diào)度提出指導(dǎo)性意見和建議。

模型的計(jì)算區(qū)域?yàn)閹?kù)區(qū)內(nèi)水域,該水域根據(jù)堤線布置進(jìn)行概化。采用Delft 3D軟件自帶Rgfgrid程序生成正交曲線網(wǎng)格。網(wǎng)格的疏密情況根據(jù)實(shí)際需要而不同,進(jìn)出水口及工程敏感區(qū)域網(wǎng)格較密,邊界附近一般區(qū)域網(wǎng)格相對(duì)較疏?？偩W(wǎng)格數(shù)82×82,最小網(wǎng)格間距15m,最大網(wǎng)格間距80m,平均網(wǎng)格間距30 m,計(jì)算面積約為6.4 km2。

4.2.2 底泥鹽度釋放模型

鄭徐水庫(kù)水體鹽分主要來源于底泥,底泥鹽度釋放是一個(gè)不斷衰減的過程。廣東工業(yè)大學(xué)彭進(jìn)平等對(duì)河口蓄淡水庫(kù)底泥鹽分釋放進(jìn)行了研究,建立了底泥鹽度釋放模型,分析了3種不同底質(zhì)情況下的底泥釋放規(guī)律,確定了釋放衰減系數(shù)、初始釋放通量等參數(shù),并應(yīng)用該模型對(duì)鑒江河口水庫(kù)的淡化進(jìn)行了分析預(yù)測(cè)。

底泥中鹽度釋放量一般隨時(shí)間呈指數(shù)衰減規(guī)律:

式中:F0為初始釋放通量(mg/(d.cm2)),k為衰減系數(shù)(1/d),t為時(shí)間(d)。

以天為計(jì)算時(shí)段進(jìn)行計(jì)算,根據(jù)質(zhì)量守恒原理,采用完全混合系統(tǒng)水質(zhì)模型,模型計(jì)算表達(dá)式為:

式中:S為水庫(kù)庫(kù)面面積(m2),F為底泥中的鹽分釋放通量(mg/(d.cm2)),Cout為完全混合后的水質(zhì)濃度(mg/L),Cin為來水水質(zhì)濃度(mg/L),Qout為流出系統(tǒng)的流量(m3/s),Qin為流入系統(tǒng)的流量(m3/s),V為水庫(kù)的庫(kù)容(m3)。

鄭徐水庫(kù)淡化預(yù)測(cè)擬采用鑒江水庫(kù)的研究成果[4],對(duì)鄭徐水庫(kù)的淡化進(jìn)行定量預(yù)測(cè),并確定水庫(kù)正常運(yùn)行情況下的排咸泵規(guī)模。

相關(guān)參數(shù)的確定:底泥風(fēng)干后的鹽度7.07 g/kg,本底間隙水鹽度8.80mg/L,初試釋放通量為8.97 mg/(d.cm2)。

5 計(jì)算成果分析

5.1 水動(dòng)力模型計(jì)算成果

(1)水庫(kù)中水體交換的主要?jiǎng)恿碜杂谝潘?由于水庫(kù)庫(kù)面面積較大,在各類引排水條件下,水庫(kù)整體水動(dòng)力條件較弱,大部分區(qū)域流速小于0.01 m/s,遠(yuǎn)離引排水閘的岸邊水體基本處于靜止?fàn)顟B(tài)。

(2)在典型引水條件下,水庫(kù)中含鹽份水體被來水稀釋,鹽度有明顯下降,下降幅度與引水口距離成反比,與水深成反比?？梢?在前期水庫(kù)鹽度較大時(shí),采用引排水調(diào)度進(jìn)行集中洗鹽,可以較大程度地降低水庫(kù)水體的含鹽度,在正常運(yùn)行期時(shí),水體靜置時(shí)間較長(zhǎng),庫(kù)內(nèi)水體鹽度出現(xiàn)分層,此時(shí),采用小流量排咸泵進(jìn)行排咸,可以達(dá)到較好的排咸效果。鄭徐水庫(kù)引水閘附近觀察點(diǎn)鹽度分層示意見圖1。

圖1 鄭徐水庫(kù)引水閘附近觀察點(diǎn)鹽度分層示意圖

(3)湖心島位置正對(duì)引水口,其南北兩側(cè)水道的過水?dāng)嗝婷娣e相對(duì)大小對(duì)水庫(kù)流場(chǎng)影響較大?？s窄南側(cè)水道可增大東北側(cè)水體的流動(dòng)速度,擴(kuò)大庫(kù)區(qū)水體流動(dòng)范圍,不僅增強(qiáng)了庫(kù)水的流動(dòng)性,而且更為有效地降低水庫(kù)水體鹽度,鑒于此,推薦縮窄南側(cè)水道方案。

(4)將庫(kù)區(qū)東側(cè)徐家浦九塘北側(cè)的放水閘沿堤線向北移動(dòng)300m,拉大了與引水口的距離,可使南部庫(kù)區(qū)的水體流動(dòng)性較放水閘北移前有明顯增強(qiáng),推薦放水閘北移方案。

5.2 鹽度預(yù)測(cè)模型計(jì)算成果

在水庫(kù)運(yùn)行初期,通過引、放水閘及泵(包括排咸泵)的組合調(diào)度,對(duì)水庫(kù)進(jìn)行洗鹽,逐步降低水庫(kù)水體鹽度,減小底泥鹽度釋放強(qiáng)度,最終實(shí)現(xiàn)用較小流量的排咸泵達(dá)到水庫(kù)供水水質(zhì)要求的目的。根據(jù)GB/T 19923—2005《工業(yè)用水水質(zhì)標(biāo)準(zhǔn)》,氯離子的含量標(biāo)準(zhǔn)為250 mg/L。水庫(kù)鹽度隨時(shí)間變化見圖2。

經(jīng)計(jì)算建成后第1年引排水2次,然后水庫(kù)蓄至4.5 m開始供水,保持水位不變,開始用排咸流量為0.34 m3/s的排咸泵排咸,再運(yùn)行1 a,水庫(kù)達(dá)到工業(yè)供水水質(zhì)要求。

圖2 水庫(kù)鹽度隨時(shí)間變化圖

6 工程措施

6.1 排咸點(diǎn)布置

為達(dá)到排咸效果,取水口布置于水體交換少、高程較低處。本工程進(jìn)水口位于西南角,出水口位于東側(cè)、西南角兩處,相對(duì)而言庫(kù)區(qū)北側(cè)區(qū)域水體交換緩慢,鹽水易在此沉積。且北側(cè)距外海最近,排咸距離短,運(yùn)行成本低。擬定沿北堤布置4處排咸點(diǎn)。將取水口設(shè)置于-2.8 m高程,取水口處需設(shè)置集咸取水坑,尺寸為:以取水口群為中心,半徑20 m,取水坑底暫擬定高程不高于-4.0 m,一方面有利于抽水水流流態(tài),另一方面也可使鹽水更集中于取水口。

6.2 取水流量及水泵選型

水庫(kù)投入使用后排咸設(shè)施主要維持水庫(kù)水質(zhì),降低水庫(kù)返鹽影響程度。根據(jù)上述分析泵總流量0.34 m3/s,布設(shè)4個(gè)抽排點(diǎn),按每個(gè)抽排點(diǎn)4個(gè)吸水口計(jì),每個(gè)抽排點(diǎn)流量為0.085 m3/s,考慮運(yùn)行過程中存在停電等因素,適當(dāng)提高單個(gè)抽排點(diǎn)流量為0.10 m3/s。

進(jìn)出水管用抗腐蝕、小阻尼管材,進(jìn)水管考慮新建壩體沉降較大采用適應(yīng)變形能力較強(qiáng)的PE排污管,出水管埋于海堤內(nèi)采用預(yù)應(yīng)力鋼筋混凝土,經(jīng)濟(jì)管徑采用規(guī)范公式:

式中:Q為流量 (m3/s)。計(jì)算結(jié)果為進(jìn)水管徑150 mm,出水管徑195 mm,根據(jù)現(xiàn)有規(guī)格進(jìn)水管采用DN160(長(zhǎng)度小于600 m)和DN200(長(zhǎng)度大于600 m),出水管采用DN250,進(jìn)水管設(shè)計(jì)流速1.32 m/s(長(zhǎng)度小于600 m)和0.85m/s(長(zhǎng)度大于600 m),出水管設(shè)計(jì)流速1.63 m/s。

進(jìn)水口最低水位為2.0 m,出水口高程為5.0 m,水頭差3 m。進(jìn)水管沿程水頭損失共計(jì)5.35m(長(zhǎng)度小于600m)和3.02 m(長(zhǎng)度大于600 m),進(jìn)水管局部水頭損失忽略,出水管沿程水頭損失共計(jì)1.55 m,出水管局部水頭損失共計(jì)0.26 m,總水頭損失約為10.18 m,水泵采用80～280(Ⅰ)-B(電機(jī)功率11 kW)離心泵,揚(yáng)程19 m,允許氣蝕余量2.3 m,吸程按下式計(jì)算:

Hs=10.33(標(biāo)準(zhǔn)大氣壓)-2.3(氣蝕余量)-0.5(安全富余量)=7.53 m

水庫(kù)排咸最低水位2.0 m,水泵安裝高程不高于5.0 m,本次設(shè)計(jì)擬定為4.8 m。

6.3 其它排咸措施

6.3.1 合理調(diào)度,加大庫(kù)內(nèi)水體交換速率

鄭徐水庫(kù)為典型的海涂平原水庫(kù),自身集雨面積只有水庫(kù)水面,來水主要是從西側(cè)八塘橫江的河道及曹娥江引水工程的境外來水。二維水動(dòng)力模型計(jì)算成果表明,庫(kù)內(nèi)水動(dòng)力條件較弱,局部庫(kù)灣水流不暢,水體環(huán)流主要來自于風(fēng)場(chǎng)、泵站引水及排水閘放水。在工程措施排咸的基礎(chǔ)上,充分利用水庫(kù)引排水,采用合理的調(diào)度方式,稀釋庫(kù)內(nèi)鹽度高的水體,再引入河網(wǎng)中的淡水或天然降水進(jìn)行補(bǔ)充,循環(huán)往復(fù),達(dá)到水質(zhì)凈化的目的。

(1)水庫(kù)建成初期,水體鹽分主要來源于水庫(kù)的底質(zhì)土層,此時(shí)水庫(kù)含鹽度較高,鹽分從底部向庫(kù)面水體擴(kuò)散明顯,此時(shí)應(yīng)利用引排水進(jìn)行集中排咸,具體的調(diào)度方式應(yīng)綜合考慮降雨、河網(wǎng)水位及閘泵規(guī)模等因素。

(2)由于慈溪60%以上的降雨集中在4—9月,在曹娥江引水的補(bǔ)充下,汛期水量相對(duì)充足,應(yīng)充分利用汛期富足的棄水和雨水進(jìn)行連續(xù)置換水。置換水方式選擇排水閘和管道排咸相結(jié)合,排水閘換水起排時(shí)間選擇在東河區(qū)河網(wǎng)水位位于低水位時(shí),水庫(kù)放水時(shí)與徐家浦十塘閘聯(lián)合放水,當(dāng)河網(wǎng)水位達(dá)到警戒水位時(shí)應(yīng)停止放水,并開閘引水,待河網(wǎng)水位降至常水位以下時(shí)再開始放水。

(3)在非汛期,水庫(kù)取排水量相對(duì)減少,水體交換不足,底部水體鹽分堆積,且不易向表層水體擴(kuò)散,使底部水體含鹽度提高,在遭遇一定天氣條件下,枯水期水庫(kù)泛咸的可能性增大,此時(shí),可采用排咸泵進(jìn)行小流量排咸,降低水庫(kù)水質(zhì)災(zāi)害的風(fēng)險(xiǎn)。

6.3.2 科學(xué)布置,優(yōu)化水庫(kù)取用水裝置

由于庫(kù)區(qū)上、中層水體鹽度分布較為均勻,且易于淡化,因而可以采用表層浮球式取水裝置取用庫(kù)區(qū)已達(dá)到水質(zhì)標(biāo)準(zhǔn)的局部合格水,提高水庫(kù)的利用率,同時(shí)也可以加快底部水體淡化。

6.3.3 加強(qiáng)監(jiān)測(cè),及時(shí)采取有效措施防治水質(zhì)災(zāi)害

在正常情況下,液體溶解擴(kuò)散和沉淀使水庫(kù)水含氯度與水位高程成梯度分布。水庫(kù)在運(yùn)行期間,由于局部水體置換率低,若不能及時(shí)排出,將會(huì)形成化學(xué)成分復(fù)雜的高濃度咸水,在干旱及雨量減少等氣候環(huán)境影響下,產(chǎn)生庫(kù)水泛咸,嚴(yán)重影響水庫(kù)供水。因此,及時(shí)掌握水庫(kù)水質(zhì)及分布情況對(duì)于水庫(kù)安全供水尤為重要。

7 結(jié) 語

采用Delft 3D水動(dòng)力模型能很好地模擬庫(kù)區(qū)的水動(dòng)力流場(chǎng)和引排水調(diào)度對(duì)水庫(kù)鹽度變化的影響,能提供水庫(kù)進(jìn)水、排水口合理布置、排咸措施選擇和引排水運(yùn)行調(diào)度提出指導(dǎo)性意見和建議。鹽度預(yù)測(cè)模型根據(jù)引、放水閘及泵(包括排咸泵)的組合調(diào)度,對(duì)水庫(kù)進(jìn)行洗鹽,逐步降低水庫(kù)水體鹽度,減小底泥鹽度釋放強(qiáng)度,最終實(shí)現(xiàn)用較小流量的排咸泵達(dá)到水庫(kù)供水水質(zhì)要求的目的,可按要求合理預(yù)測(cè)水庫(kù)咸水淡化時(shí)間,可供其他海涂水庫(kù)排咸設(shè)計(jì)參考。

[1]金德鋼.慈溪市鄭徐水庫(kù)初步設(shè)計(jì)報(bào)告[R].寧波:寧波市水利水電規(guī)劃設(shè)計(jì)研究院,2010.

[2]毛獻(xiàn)忠,陳甫源,余祈文,等.堵港蓄淡水庫(kù)水體淡化預(yù)測(cè)研究 [J].水利學(xué)報(bào),2004(07):80-85.

[3]毛獻(xiàn)忠,朱小敖,陳甫源,等.沿海堵港蓄淡水庫(kù)加快水體淡化措施的研究[J].水科學(xué)進(jìn)展,2005(06):773-776.

[4]向軍,彭進(jìn)平,逢勇,等.鑒江河口底泥鹽度釋放規(guī)律研究[J].水文,2008,28(4):12-15.