宋文超,萬 俊,符 強,侯曉輝,邢寶龍
(1.河南省水利勘測設(shè)計研究有限公司,河南 鄭州450016;2.中牟縣黃河河務(wù)局,河南中牟451450)
生態(tài)湖具有雨洪調(diào)蓄、景觀補水、旅游觀光等多種功能。生態(tài)湖調(diào)度運行即在一定時間內(nèi)對湖內(nèi)水體進行更新與替換,使其維持一定水位、水質(zhì)達到某一標(biāo)準(zhǔn)。目前,關(guān)于水庫調(diào)度運行方式的研究較多,但傳統(tǒng)的生態(tài)湖調(diào)度方式多考慮水量要求,而忽視了水質(zhì)目標(biāo)。近年來,隨著計算機技術(shù)的進步,數(shù)值模擬成為了研究生態(tài)湖調(diào)度運行方式的有效方法,魏珂等[1]根據(jù)洪湖多年引水水質(zhì)構(gòu)建了洪湖生態(tài)環(huán)境用水量模型;羅佳翠等[2]通過二維水環(huán)境模型進行了牛欄江引水效果預(yù)測,得到了滇池水質(zhì)基本達標(biāo)所需的調(diào)水量;哈佳等[3]利用MIKE21水動力模型和對流擴散模型系統(tǒng)模擬了人工湖不同點位水體交換完成時間,確定了引水流量、引水周期、湖區(qū)重污染潛在區(qū)域和富營養(yǎng)化高發(fā)區(qū)。
合理的調(diào)度運行方式對維持水體水動力條件以及換水周期,增強水體自凈能力,加速水體中污染物降解,提高水體水環(huán)境質(zhì)量,維持水體健康,滿足生態(tài)和景觀需求等具有重要意義。為了模擬人工生態(tài)湖建成后的水質(zhì)狀況,以鄭州象湖為例,以城市水體水質(zhì)保護為目標(biāo),對人工生態(tài)湖水體運行方式進行研究,旨在尋找符合水質(zhì)要求的大面積水體運行方式,以滿足城市水生態(tài)建設(shè)要求。
白沙園區(qū)位于河南省鄭州市鄭汴新區(qū)中西部,是鄭州都市區(qū)重要拓寬空間和中原經(jīng)濟區(qū)核心增長重要組成部分。象湖生態(tài)工程是白沙園區(qū)生態(tài)水系工程的一個重要組成部分,位于鄭開大道與賈魯河交叉處,分布于鄭開大道兩側(cè),東至前程路,西至楊橋路,南至高莊、白墳兩村,北至升平路、祭城路之間。象湖一期工程湖體水面面積37萬m2,水體117萬m3,包括水源工程、一期湖體開挖工程、防滲工程和退水工程。一期水源工程指引水工程,設(shè)計通過沿賈魯河左岸新開挖渠道從楊橋干渠引黃河水入象湖,渠道長約390 m,設(shè)計引水流量為1 m3/s。一期工程位于賈魯河、鄭開大道以北區(qū)域,工程總開挖面積約為41.4萬m2,開挖方量約為296.25萬m3,一期工程建設(shè)完成蓄水后,總水面面積可達37萬m2,其中淺水區(qū)、過渡區(qū)及深水區(qū)面積分別為4萬、24萬、9萬m2。湖周護岸采用垂直擋墻式護岸和自然生態(tài)式護岸,水質(zhì)控制目標(biāo)為Ⅳ類。
象湖水體不存在明顯分層現(xiàn)象,采用二維模型即可滿足研究需要,經(jīng)過模型比選,水動力模型采用MIKE21模型中的水動力模塊。MIKE21模型中的水動力模塊(HD模塊)是其核心模塊[4],可以模擬因各種作用力而產(chǎn)生的水位及水流變化,并為水質(zhì)變化、泥沙輸移和溢油污染等其他模塊提供必要的水動力學(xué)基礎(chǔ),可用于忽略分層的二維自由表面流的模擬[5-6]。基本控制方程:
式中:t為時間;η為水面相對于未擾動水面的高度;h為靜止水深;u、v分別為流速在x、y方向上的分量;τsx、τsy分別為自由表面風(fēng)在x、y方向上對水面的剪切應(yīng)力;τbx、τby分別為在x、y方向上的底床摩擦應(yīng)力項;Pa為當(dāng)?shù)卮髿鈮?;ρ為水密度;?為參考水密度;f=2Ωsinφ為科氏力參數(shù)(其中Ω為地球自轉(zhuǎn)角速率,Ω=0.729×10-4s-1,φ 為地理緯度);g 為重力加速度;fˉv、fˉu為地球自轉(zhuǎn)引起的加速度;sxx、sxy、syx、syy為輻射應(yīng)力分量;Txx、Txy、Tyx、Tyy為水平黏滯應(yīng)力項;S 為源匯項;us、vs為源匯項水流流速[6]。
水質(zhì)模型采用ECOLab,ECOLab是MIKE系列模擬軟件中的一個功能模塊,可以模擬水生態(tài)系統(tǒng)中變量之間的轉(zhuǎn)化過程和相互作用,同時可用于描述各種狀態(tài)變量的沉降、懸浮、吸附和解吸等物理過程,可用于湖泊、河流、濕地水質(zhì)模擬[7-8]?;究刂品匠探M:
式中:PC、PN、PP分別為浮游植物含碳量、浮游植物含氮量、浮游植物含磷量;CH為葉綠素含量;IN、IP分別為無機氮和無機磷含量;prpc為浮游植物的產(chǎn)碳量;prch為浮游植物光合作用生產(chǎn)的葉綠素量;unpn、uppp分別為浮游植物攝取氮和磷的量;grpc、grpn、grpp分別為浮游植物被浮游動物牧食后損失的碳、氮、磷的量;depc、depn、depp、dech分別為因浮游植物死亡而損失的碳、氮、磷和葉綠素的量;sepc、sepn、sepp、sech 分別為隨浮游植物下沉損失的碳、氮、磷和葉綠素的量;redn、redp分別為水體中的碎屑礦化分解形成的無機氮和無機磷的量;rezn、rezp分別為水體中的浮游動物在呼吸作用的同時釋放出的無機氮和無機磷量;depn2in和depp2ip分別為浮游植物死亡后釋放的無機氮和無機磷的量;unbn和upbp分別為底棲植物攝取無機氮和無機磷的量;rebn、rebp分別為隨底棲植物呼吸作用釋放的無機氮和無機磷的量。
結(jié)合象湖一期工程規(guī)劃設(shè)計地形資料,構(gòu)建象湖二維水動力水質(zhì)模型,模型范圍東西向為1.2 km,南北向為0.6 km,湖區(qū)網(wǎng)格單元格邊長約10 m,局部地形變化較大區(qū)域網(wǎng)格加密至5 m,模型網(wǎng)格節(jié)點總數(shù)1 087個,網(wǎng)格劃分見圖1。
圖1 計算區(qū)域及網(wǎng)格劃分
水動力模型邊界條件:湖區(qū)入口采用給定流量邊界條件,出口采用給定水位邊界條件;象湖水質(zhì)模型考慮了自然風(fēng)、降雨蒸發(fā)及大氣沉降的影響,并根據(jù)不同的水質(zhì)條件設(shè)置入流水質(zhì)邊界,其中引黃水補水水質(zhì)條件采用經(jīng)過沉淀池沉淀后的引黃水實測數(shù)據(jù)資料,見表1。
表1 引黃水經(jīng)沉淀池沉淀后的水質(zhì)指標(biāo)
(1)流場分析。分析了不同入流方案對湖區(qū)流場的影響。方案一入流方向與賈魯河縱向大致平行,往北略微偏向湖區(qū);方案二入流方向與賈魯河相交約30°,往南偏向賈魯河。湖區(qū)流場見圖2、圖3。從圖中可以看出,由于湖體進出口分布在湖區(qū)西部,因此湖區(qū)西部水體流動較為明顯,而湖區(qū)東部水體流動較弱。在湖區(qū)西部,水流從入口進入湖體并擴散到周圍區(qū)域,并逐步流向下游,最后從湖區(qū)出口流出湖體。方案一和方案二湖區(qū)水體流動的主要區(qū)別在湖區(qū)入口區(qū)域,方案一主流在北側(cè)湖區(qū),在南側(cè)靠近賈魯河處形成了回旋區(qū),方案二主流在南側(cè)靠近賈魯河處,在湖區(qū)北側(cè)形成了回旋區(qū)。
圖2 方案一湖區(qū)流場
圖3 方案二湖區(qū)流場
(2)換水周期分析。不同方案下湖區(qū)換水周期分布見圖4、圖5。從圖中可以看出,湖區(qū)西部換水效果較好,而湖區(qū)東部水體流動較弱,換水周期較長。在湖區(qū)西部,主流區(qū)換水時間短,效果較好,回旋區(qū)換水時間相對較長。方案一西部湖區(qū)北側(cè)換水效果較好,南側(cè)換水周期相對較長;方案二西部湖區(qū)南側(cè)換水效果較好,北側(cè)換水周期較長。由不同進口水流入湖方向湖區(qū)換水率(見表2)可知,方案一湖區(qū)平均實際換水率較方案二的大,說明湖區(qū)進口水流入湖方向與賈魯河縱向大致平行相較于與賈魯河大角度相交更有利于湖區(qū)整體換水。
圖4 方案一湖區(qū)換水周期分布
圖5 方案二湖區(qū)換水周期分布
(3)水系補源方案分析。水系補源在方案一(入流方向與賈魯河縱向大致平行)基礎(chǔ)上采用如下方案:方案A,補水、退水同時進行,按1 m3/s流量進水,13.5 d換水一次;方案 B,換水時從正常蓄水位80.36 m先退水0.50 m,水位下降至79.86 m時,開始補水,補水流量為 1 m3/s,13.5 d 換水完畢。
表2 不同方案湖區(qū)換水率
自然風(fēng)條件下,兩種補源方案湖區(qū)換水周期分布見圖6。從圖中可以看出,兩種補源方案湖區(qū)換水周期分布規(guī)律相近。經(jīng)模型計算,自然風(fēng)條件下,在1 m3/s的換水流量下,湖區(qū)西部換水周期小于22 d,湖區(qū)東部換水周期22~27 d,補水方案A湖區(qū)平均換水周期26 d,補水方案B湖區(qū)平均換水周期23 d。
圖6 自然風(fēng)條件下湖區(qū)換水周期分布
經(jīng)計算,自然風(fēng)條件下,采用補水退水同時進行的補源方案實際換水率為65.9%,采用先退水再補水的補源方案實際換水率為66.6%。自然風(fēng)條件下,采用先退水再補水方式的湖區(qū)換水效果略優(yōu)于補水退水同時進行方式的。
(1)水華風(fēng)險分析。MIKE ECO Lab軟件的優(yōu)勢不僅僅在于可以模擬湖區(qū)營養(yǎng)鹽濃度在水體中的分布情況,而且能夠模擬湖區(qū)藻類在合適溫度、光照等外界條件下的生長情況。因此通過ECO Lab水質(zhì)模型模擬在高溫、自然光照等條件及營養(yǎng)鹽較充足情況下,湖區(qū)藻類生長與時間的響應(yīng)關(guān)系,以此得到防止湖區(qū)出現(xiàn)水華風(fēng)險的初步換水周期和換水量,為后續(xù)換水工況設(shè)置提供參考。
在 1 m3/s流量(即理論換水周期 13.5 d)連續(xù)入流條件下,氣溫、光照等自然條件良好及營養(yǎng)鹽較充足時,評估藻類的生長變化。通過ECO Lab水質(zhì)模型模擬得到湖區(qū)葉綠素濃度逐漸增大(見圖7),經(jīng)過18 d左右葉綠素濃度已達到0.030 mg/L,接近《地表水資源質(zhì)量評價技術(shù)規(guī)程》(SL 395—2007)湖泊營養(yǎng)狀態(tài)評價標(biāo)準(zhǔn)及分級方法中葉綠素濃度大于0.026 mg/L的富營養(yǎng)狀態(tài),可以認(rèn)為此時湖區(qū)有一定的水華風(fēng)險。
圖7 葉綠素濃度變化情況
通過以上模擬結(jié)果分析可知,在高溫、營養(yǎng)鹽負(fù)荷充足情況下,湖區(qū)在連續(xù)入流條件下保持時間不宜超過18 d,否則需要加大換水量或者減少入流污染負(fù)荷。
(2)生態(tài)湖水質(zhì)年變化分析。模擬生態(tài)湖正式運行后不同補換水次數(shù)對湖區(qū)水質(zhì)的影響。在正常運行時,根據(jù)黃河調(diào)水調(diào)沙時間,6月下旬至7月中旬停止引水,8月停止10 d引水的基本運行方式來設(shè)置工況,確保在合理換水周期下水質(zhì)達標(biāo)。換水方式采用在一個換水周期保持連續(xù)換水,模擬工況為在入流方向與賈魯河縱向大致平行,先退水再補水的水系補源方案基礎(chǔ)上換水3、4、6、7次4種。工況一換水3次,換水時間為6月上旬、7月下旬、8月中下旬;工況二換水4次,換水時間為6月上旬、7月下旬、8月中下旬、9月上旬;工況三換水6次,換水時間為6月上旬、7月下旬(2次,每次 5 d)、8月中下旬(2 次,每次 5 d)、9 月上旬;工況四換水7次,換水時間為6月上旬、7月下旬(2次,每次 5 d)、8月中旬(2次,每次 5 d)、8 月下旬、9月上旬。
在以上工況條件下模擬計算得到TN、TP、葉綠素濃度年變化情況,見圖8~圖10。春季氣溫較低,藻類生長基本停滯,TN、TP濃度維持在較低水平。5月后,藻類生長加快,湖區(qū)營養(yǎng)鹽濃度較充足,6月和7月開始富集,導(dǎo)致營養(yǎng)鹽濃度較快升高。在湖區(qū)增加換水補水次數(shù)的情況下,營養(yǎng)鹽濃度開始下降,10月氣溫下降,藻類死亡沉積后開始恢復(fù)至正常水平。從全湖TN平均濃度來看,3次換水和4次換水時TN濃度在8月超過地表水Ⅳ類水質(zhì)標(biāo)準(zhǔn);6次換水和7次換水時TN濃度除7月滿足地表水Ⅳ類水質(zhì)標(biāo)準(zhǔn)外,其他時段均在Ⅲ類水質(zhì)標(biāo)準(zhǔn)范圍內(nèi)。從全湖TP平均濃度來看,3次換水時在7月下旬后難以達到地表水Ⅳ類水質(zhì)標(biāo)準(zhǔn),4次換水時TP濃度在7月下旬至8月底超過地表水Ⅳ類水質(zhì)標(biāo)準(zhǔn),6次換水和7次換水能確保TP濃度在Ⅳ類水質(zhì)標(biāo)準(zhǔn)范圍內(nèi)。
圖8 TN濃度年變化曲線
圖9 TP濃度年變化曲線
圖10 葉綠素濃度年變化曲線
根據(jù)《地表水資源質(zhì)量評價技術(shù)規(guī)程》(SL 395—2007)中湖泊營養(yǎng)狀態(tài)評價標(biāo)準(zhǔn)及分級方法,TN濃度大于 1 mg/L、TP 濃度大于 0.1 mg/L、葉綠素濃度大于0.026 mg/L時,可以認(rèn)為水域處于中度富營養(yǎng)化狀態(tài),湖區(qū)水質(zhì)處于風(fēng)險狀態(tài)。3次換水和4次換水時7月和8月為水華風(fēng)險時段;由于氣溫較高,因此在營養(yǎng)鹽充足情況下,6次換水和7次換水雖然TN、TP濃度減小,但是葉綠素濃度在7月底至8月初仍然較高,在高溫天氣條件下也具有一定的水華風(fēng)險。
綜上所述,從指標(biāo)濃度的年變化曲線看,3次換水和4次換水TN和TP濃度在夏季會超過地表水Ⅳ類水質(zhì)標(biāo)準(zhǔn);6次換水和7次換水夏季TN、TP濃度在地表水Ⅳ類水質(zhì)標(biāo)準(zhǔn)內(nèi);7次換水較6次換水對9月以后生態(tài)湖水質(zhì)條件有所改善,但對水華風(fēng)險的防控影響不大。
(3)生態(tài)湖水質(zhì)優(yōu)化措施模擬分析。根據(jù)上述模擬結(jié)果并考慮實際運行管理方式,采用夏季4次換水基礎(chǔ)上削減來流污染物及人工濕地處理方式來保證湖區(qū)水質(zhì)達標(biāo),并對此進行分析。夏季4次換水時間為7月中下旬2次,8月上中旬1次,8月下旬到9月中旬1次,春季及冬季視水質(zhì)情況決定是否換水。水質(zhì)改善措施主要采用地下水補水方式進行換水,或者在環(huán)湖濱帶布置濕地對湖水進行處理。優(yōu)化方案見表3。
表3 優(yōu)化方案
TN、TP及葉綠素濃度年變化曲線見圖11~圖13。優(yōu)化方案在整體上與原方案變化趨勢相似。春季氣溫較低,藻類生長基本停滯,TN、TP濃度維持在較低水平;5月后,藻類生長加快,湖區(qū)營養(yǎng)鹽濃度較高,6—7月營養(yǎng)鹽開始富集,導(dǎo)致營養(yǎng)鹽濃度較快升高;在湖區(qū)增加換水的情況下,營養(yǎng)鹽濃度開始下降,至10月后開始恢復(fù)至正常水平。從7月下旬換水開始,由于來流濃度不同及生態(tài)優(yōu)化措施的實施,各指標(biāo)濃度開始分化。原方案各指標(biāo)濃度最高,優(yōu)化方案一水質(zhì)最好,優(yōu)化方案三次之。
圖11 優(yōu)化方案TN濃度年變化曲線
圖12 優(yōu)化方案TP濃度年變化曲線
通過對4種方案TN、TP及葉綠素全湖濃度的對比,如果采用4次換水方案,原方案只采用黃河水時TP指標(biāo)超出了地表水Ⅳ類標(biāo)準(zhǔn),在夏季無法保證生態(tài)湖湖區(qū)水質(zhì),需要另外補充水質(zhì)較好的地下水或者對引黃水進行處理。從圖11~圖13可以看出,優(yōu)化方案一、優(yōu)化方案二及優(yōu)化方案三均可以滿足生態(tài)湖湖區(qū)水質(zhì)要求。
圖13 優(yōu)化方案葉綠素濃度年變化曲線
為減少地下水使用量、降低運行成本并考慮生態(tài)措施,建議采用優(yōu)化方案三,即引黃水加環(huán)湖濕地的生態(tài)運行方式。濕地處理水量9 400 m3/d,處理出水水質(zhì)達到Ⅲ類水質(zhì)標(biāo)準(zhǔn)下,經(jīng)模型模擬可以保持湖區(qū)地表水Ⅳ類水質(zhì)。根據(jù)《人工濕地污水處理工程技術(shù)規(guī)范》(HJ 2005—2010),采用垂直潛流式人工濕地來保持湖區(qū)水質(zhì),水力負(fù)荷 1 m3/(m2?d),停留時間 1 d,規(guī)劃濕地面積需約0.080 km2。
根據(jù)有機負(fù)荷來核算濕地面積:
A = Q(ln So-ln Sc) /(KDn) (10)
式中:A 為濕地表面積,m2;Q 為進水流量,m3/d;So 為進水指標(biāo)濃度,mg/L;Sc為出水指標(biāo)濃度,mg/L;K 為反應(yīng)速率系數(shù),d-1,取 1.36;D 為濕地床深度,m,取 1.2 m;n為孔隙度,取30%。
計算得到濕地面積為0.048 km2。綜合上述,采用環(huán)湖濕地生態(tài)運行方式,濕地面積為 0.048~0.080 km2。
湖區(qū)進口水流入湖方向與賈魯河縱向大致平行時,湖區(qū)東部水體的換水周期較西部水體的換水周期長。 在理論換水周期13.5 d的換水流量(1 m3/s)下,湖區(qū)西部換水周期小于22 d,湖區(qū)東部換水周期為22~27 d。
湖區(qū)水質(zhì)各指標(biāo)濃度年際變化規(guī)律基本相同。春季氣溫較低,藻類生長基本停滯,TN、TP濃度維持在較低水平。5月后,藻類生長加快,湖區(qū)營養(yǎng)鹽濃度增大,6月和7月開始富集,導(dǎo)致營養(yǎng)鹽濃度較快升高。10月氣溫開始下降,藻類死亡并沉積,湖區(qū)營養(yǎng)鹽濃度下降,并恢復(fù)至正常水平。
為保持湖區(qū)水質(zhì)在夏季惡劣情況下不超過地表水Ⅳ類、平時保持在Ⅲ類水質(zhì)水平,建議采用6月?lián)Q水1次、7月?lián)Q水2次、8月?lián)Q水2次、9月?lián)Q水1次的方式來運行。如遭遇持續(xù)惡劣高溫天氣,可考慮及時換水以保證湖區(qū)良好水質(zhì)環(huán)境??紤]實際運行管理方式,建議采用4次換水方式,通過建設(shè)垂直潛流式人工濕地來保持湖區(qū)水質(zhì)持續(xù)穩(wěn)定,濕地面積為0.048~0.080 km2。