孫 楊,韓龍喜

典型暴雨農(nóng)業(yè)面源總磷入庫水質(zhì)影響過程

孫 楊,韓龍喜

(河海大學(xué)環(huán)境學(xué)院,江蘇南京 210098)

以東北大伙房水庫為典型研究案例,根據(jù)典型暴雨情況下農(nóng)業(yè)面源中總磷入庫的時間變化特征,建立考慮入庫水流影響的庫區(qū)三維水動力數(shù)學(xué)模型及相應(yīng)的三維水質(zhì)模型,模擬分析1次典型暴雨入庫過程中庫區(qū)總磷三維空間分布及隨時間變化的特征。模擬結(jié)果表明,受入庫流量和豐水期主導(dǎo)風(fēng)向東北風(fēng)的影響,總體上總磷隨著水流主要沿西南岸線向壩址方向擴散,總磷濃度超標(biāo)面積逐漸增大,隨著流程增加,其濃度逐漸減小。對面源入庫過程的模擬分析可以為面源污染的預(yù)報及治理提供相關(guān)的技術(shù)依據(jù)。

典型暴雨;農(nóng)業(yè)面源污染;總磷;三維數(shù)值模擬;大伙房水庫

隨著經(jīng)濟發(fā)展和城鎮(zhèn)規(guī)模的擴大,通過點源治理,水體污染已得到明顯控制。但由于城鎮(zhèn)垃圾、農(nóng)田化肥、農(nóng)藥等隨雨水產(chǎn)生的地表徑流進入了水體,其面源污染治理沒有得到重視,使得面源污染在環(huán)境污染中所占的比例越來越大[1]。水體面源污染一般是晴天積累,雨天排放,沒有固定的發(fā)生源。而降水徑流則是其發(fā)生的主要動力,因此當(dāng)一場暴雨發(fā)生時,面源污染隨之發(fā)生。面源入庫源強及庫區(qū)水質(zhì)響應(yīng)是一個動態(tài)過程。對這一過程的進行模擬分析可為面源污染的預(yù)報及治理提供相關(guān)的技術(shù)依據(jù)。

目前,針對農(nóng)業(yè)面源污染的預(yù)報工作,國內(nèi)外取得了一系列研究成果。范麗麗等[2]應(yīng)用SWAT模型計算分析了大寧河流域不同降雨徑流條件下非點源污染的產(chǎn)輸出特性。美國國家環(huán)保局開發(fā)出SWMM暴雨徑流管理模型模擬城市降雨面源污染情況[3-4]。蔣金等[5]通過分析降雨徑流的水文過程、其與非點源污染物的作用機理及其對非點源污染的影響,探討如何有效地減少非點源污染。程紅光等[6]以黑河流域為研究區(qū)對不同降水條件下農(nóng)田、城鎮(zhèn)、草地、灌木及森林等土地利用類型的各亞流域氮的入河系數(shù)進行了研究。呂興娜[7]利用鐵嶺市柴河堡水文站降雨時節(jié)實測的暴雨徑流數(shù)據(jù)、泥沙數(shù)據(jù)、污染物濃度數(shù)據(jù)研究降雨時節(jié)暴雨徑流量與主要污染物濃度及污染物瞬時流失量之間的關(guān)系?？傮w而言,目前研究成果主要集中于分析面源污染物濃度與不同降雨徑流條件的關(guān)系。由于水庫有防洪、蓄水灌溉、發(fā)電、養(yǎng)魚等作用,甚至可作為當(dāng)?shù)氐娘嬘盟?因此研究暴雨情況下面源污染對水庫水質(zhì)的影響有其特殊意義。本文以大伙房水庫為例,分析典型暴雨發(fā)生后,面源污染團進入水庫對水庫水質(zhì)的影響。

1 研究區(qū)域概況

大伙房水庫是一座大型水利樞紐工程,位于遼寧省東部山區(qū),是遼寧省兩個重要工業(yè)城市——沈陽市、撫順市的飲用水水源地。其控制流域面積5437 km2,流域內(nèi)多山地。平均海拔高度175 m,水庫最大庫容21.87億m3,最大水深41.5 m,平均水深20 m,最大水面長30 km,平均水面寬l 000 m,最大水面面積110 km2[8]。其所處流域為暖溫帶大陸性季風(fēng)氣候,冬季漫長寒冷,夏季炎熱多雨,全年日照時數(shù)在2 280～2 670 h。流域多年平均年降水量為650～800 mm,降水量主要集中在7月、8月,占全年降水量的50%左右。流域內(nèi)多年平均年水面蒸發(fā)約為1 100～1 600mm,多年平均相對濕度為60%～70%,多年平均風(fēng)速為1.5～3.8 m/s[9]。

社河是距離大伙房水庫最近的一條入庫河流,流經(jīng)撫順縣多個村鎮(zhèn),農(nóng)業(yè)面源污染較為嚴(yán)重。目前社河周邊地區(qū)農(nóng)田已呈酸化趨勢,由于作物對化肥的吸收一般只有20%～30%,其余大部分營養(yǎng)元素通過地表徑流和淋溶進入水體,導(dǎo)致藻類等生物大量繁殖。水體溶解氧的急劇下降,使水質(zhì)惡化,導(dǎo)致魚蝦大量死亡。近年來,當(dāng)?shù)氐呐！⑸i、羊、雞、鴨等畜禽養(yǎng)殖數(shù)量比前幾年已有了明顯增長,而畜禽糞尿總體流失率在10%左右,社河流域由畜禽糞尿的流失而排入水環(huán)境中的污染物總量已超過該地區(qū)居民生活污染物總量[10]。因此,本文選取社河流入大伙房水庫區(qū)域為典型研究區(qū)域,見圖1所示。

圖1 社河入庫流域

大伙房水庫最大水深為41.5 m,平均水深20 m,屬于深水水庫,其流速、污染物濃度沿垂向分布差異較大。筆者根據(jù)典型暴雨發(fā)生時的入庫河流的水文、水質(zhì)及污染物濃度變化過程的監(jiān)測數(shù)據(jù),建立水庫三維水動力、水質(zhì)模型,模擬分析典型暴雨入庫過程中庫區(qū)總磷三維空間分布及隨時間變化的特征。

2 數(shù)學(xué)模型

2.1 三維水動力數(shù)學(xué)模型

用笛卡爾坐標(biāo)系下描述水庫三維水流運動的控制方程組[11]：

式中：u、v、w分別為坐標(biāo)x、y、z方向上的流速分量;t為時間坐標(biāo);ρ0為水體密度;p為壓強;g為重力加速度;AM為水平渦黏系數(shù);vz為垂向渦黏系數(shù);f為柯氏力系數(shù),f=2ΩsinΦ,其中,Ω=2π/86184為地球自轉(zhuǎn)頻率,Φ為當(dāng)?shù)鼐暥取?/p>

2.2 三維水質(zhì)數(shù)學(xué)模型

在笛卡爾坐標(biāo)系下,污染物三維輸運方程描述如下：

式中：ρ為污染物質(zhì)量濃度;K為污染物降解系數(shù); Dv為垂向擴散系數(shù);Dk為水平擴散系數(shù)。

2.3 模型驗證

選用明渠垂向流速分布的經(jīng)驗公式來驗證數(shù)值模擬的結(jié)果[12]。試驗選取的計算區(qū)域為東西開邊界的矩形明渠水槽,長1000m,寬400m。其他參數(shù)見表1。

表1 明渠恒定流垂向流數(shù)值模擬參數(shù)值

選用明渠垂向濃度分布的指數(shù)律來驗證數(shù)值計算結(jié)果,得出的數(shù)值模擬結(jié)果與垂向流速指數(shù)律分布公式所得的流速分布對比見圖2、圖3?？梢钥闯?模型較好地模擬了明渠水流的流速分布。

圖2 明渠垂向流速分布示意圖

圖3 明渠垂向流速分布解析解與數(shù)值解比較

3 面源入庫庫區(qū)TP質(zhì)量濃度變化特征分析

3.1 庫區(qū)三維水動力特征分析

a.網(wǎng)格布置。水動力模擬采用無結(jié)構(gòu)三角網(wǎng)格劃分計算區(qū)域,在社河入庫區(qū)域適當(dāng)加密。平面共布置2082個節(jié)點,3652個網(wǎng)格單元,垂向上采用等距離分層進行剖分,共分為10層。

b.定解條件及模擬參數(shù)取值。垂向渦黏系數(shù)vz采用Kolmogorov-Prandl方程求解,水平渦黏系數(shù)AM采用Smagorinsky亞網(wǎng)格尺度紊動模型計算。入流邊界中的蘇子河和渾河邊界采用豐水期實測平均入庫流量,分別是138.7 m3/s和100.6 m3/s,面源入庫的社河邊界取的是流量時間序列。模擬參數(shù)取值：粗糙高度ks=0.01m,水平擴散系數(shù)Dk=1m2/s,垂向擴散系數(shù)Dv=1m2/s,TP降解系數(shù)為0.0018 d-1;定解條件：初始流速為0m/s,初始水位為129.85m,出流邊界水位為129.85 m,岸邊界法向流速為0 m/s。

c.模擬時段。由圖4可知,8月4—8日這段時間內(nèi)有1場暴雨發(fā)生,因此選擇此時段作為農(nóng)業(yè)面源污染的研究時段。社河最大入庫流量478 m3/s發(fā)生在8月4日12：00,總磷的入庫質(zhì)量濃度也在8月4日20：00時達到最大值0.72 mg/L。

圖4 社河入庫流量、TP質(zhì)量濃度隨時間變化曲線

以豐水期平均風(fēng)速2.3 m/s和主導(dǎo)風(fēng)向(東北風(fēng))作為自由表面條件,模擬水庫豐水期三維水動力流場。當(dāng)模擬時間為12 h,即社河入庫流量達到最大值時,吞吐流對水庫流場的影響也達到最大。由圖5(a)可知,大伙房水庫社河入庫口水域水流結(jié)構(gòu)呈現(xiàn)受吞吐流支配的水動力特征和顯著的順岸流特征。當(dāng)模擬時間為72 h,此時入庫流量逐漸減小到最小值而主導(dǎo)風(fēng)向東北風(fēng)的大小為2.3 m/s時,風(fēng)應(yīng)力對水庫流場影響最大。由圖5(b)可知,表層流場受風(fēng)場的作用比較大,從社河流入水庫的水流形成順時針的大環(huán)流,入庫口水域水流結(jié)構(gòu)呈現(xiàn)以風(fēng)生流為支配的水動力特征。

圖5 三維水動力模擬表層流速

3.2 庫區(qū)濃度平面分布特征

大伙房水庫富營養(yǎng)化綜合指標(biāo)已接近富營養(yǎng)化水平,總磷等污染物因子呈增加趨勢,因此選取總磷作為污染因子研究[13]。大伙房水庫年均總磷質(zhì)量濃度為0.02 mg/L,湖庫Ⅱ類水TP質(zhì)量水平方向上,取表層的TP質(zhì)量濃度分布狀態(tài)進行特征分析。圖6為TP入庫后不同時間表層水體水平分布濃度場。TP以面源污染形式入庫形成污染團,在東北向水流輸移作用下,污染團主要沿西南岸線向壩址方向擴散。如圖6所示,污染物入庫16 h、20 h、24 h后,社河入庫流量較大,此時吞吐流對水庫流場的影響較大,污染團隨著水流方向向下游擴散。污染物入庫42 h、60 h、78 h后,社河入庫流量逐步減小,風(fēng)生流對水庫流場的影響越來越大,污染團向下游擴散的同時受到東北風(fēng)向的影響,有一部分污染團流向西南岸或入庫口,因此污染團偏西南岸線向下游擴散。隨著流程的增加,由于擴散作用,污染物質(zhì)量濃度逐漸減小。

圖6 表層水體TP質(zhì)量濃度場水平分布

隨著時間推移,污染團質(zhì)量濃度(ρ(TP)＞0.025 mg/L)超標(biāo)面積不斷增大,即對于受納水體,影響范圍不斷增加。面源入庫16h、20h、24h、42h、60h、78h后,TP質(zhì)量濃度最大超標(biāo)面積分別為1.92 km2、3.19 km2、3.48 km2、4.32 km2、4.47 km2、5.92 km2, TP污染團擴散最遠處到社河入庫口的距離依次是1.14km、2.09km、2.21km、3.52km、4.59km和5.19km。

3.3 庫區(qū)TP質(zhì)量濃度縱向分布特征分析

圖7 TP質(zhì)量濃度場縱向分布

垂向方向上,取如圖7所示垂線斷面濃度分布進行分析。由于污染物在社河中垂向擴散達到穩(wěn)定,因此在進入大伙房水庫后,TP質(zhì)量濃度在垂向上的擴散也達到穩(wěn)定。如圖7所示,面源入庫16 h、 20 h、24 h、42 h、60 h、78 h時,TP質(zhì)量濃度污染團中心處到社河入庫口的距離依次是0.25 km、0.73 km、1.09 km、1.40 km、2.96 km和3.93 km,沿著水流方向,污染團中心不斷向下游遷移,TP污染物質(zhì)量濃度沿流程逐漸減小。

4 結(jié) 語

以大伙房水庫為典型研究區(qū)域,建立三維水動力、水質(zhì)模型。該模型有效地模擬了在一次典型暴雨發(fā)生的情況下,農(nóng)業(yè)面源污染因子總磷進入水庫后的遷移特征。模擬結(jié)果表明,受入庫流量和豐水期主導(dǎo)風(fēng)向東北風(fēng)的影響,總體上污染物隨著水流主要沿西南岸線向壩址方向擴散,對于庫區(qū)受納水體,總磷的超標(biāo)面積逐漸增大,但隨著流程增加,污染物濃度逐漸減小。面源入庫過程對庫區(qū)水質(zhì)的影響是一個動態(tài)過程,對這一過程的模擬分析可以為水庫面源污染的預(yù)報及治理提供相關(guān)的技術(shù)依據(jù),在一定程度上提高了面源污染事故處理的水平,能夠有效地控制和減輕事故危害。

[1]李國斌,王焰新,程勝高.基于暴雨徑流過程監(jiān)測的非點源污染負荷定量研究[J].環(huán)境保護,2002(5)：46-48.(LI Guobin,WANG Yanxin,CHENG Shenggao.A quantitative study on non-point source pollutant loading basedonstormrunoffmonitoringresults[J]. EnvironmentalProtection,2002(5)：46-48.(in Chinese))

[2]范麗麗,沈珍瑤,劉瑞民.不同降雨-徑流過程中農(nóng)業(yè)非點源污染研究[J].環(huán)境科學(xué)與技術(shù),2008,31(10)：5-8.(FAN Lili,SHEN Zhenyao,LIU Ruimin.Agricultural non-point source pollution in different rainfall runoff[J]. Environmental Science&Technology,2008,31(10)：5-8. (in Chinese))

[3]CHANG C H,WEN C G,LEE C S.Use of intercepted runoffdepthforstormwaterrunoffmanagementin industrialparksinTaiwan[J].WaterResources Management,2008,22：1609-1623.

[4]PARK S Y,LEE K W,PARK I H,et al.Effect of the aggregation level of surface runoff fields and sewer network for a SWMM simulation[J].Desalination,2008,226(1/ 3)：328-337.

[5]蔣金,安娜,張義,等.水文過程中降雨徑流對非點源污染的影響[J].安徽農(nóng)業(yè)科學(xué),2012,40(6)：3529-3531. (JIANG Jin,AN Na,ZHANG Yi,et al.Influence of rainfall run-off in hydrologic process on non-point pollution [J].Journal of Anhui Agricultural Sciences,2012,40 (6)：3529-3531.(in Chinese))

[6]程紅光,郝芳華,任希巖,等.不同降雨條件下非點源污染氮負荷入河系數(shù)研究[J].環(huán)境科學(xué)學(xué)報,2006,26 (3)：392-397.(CHENG Hongguang,HAO Fanghua,REN Xiyan,et al.The study of the rate loss of nitrogenous nonpoint source pollution loads in different precipitation levels [J].Acta Scientiae Circumstantiae,2006,26(3)：392-397.(in Chinese))

[7]呂興娜.降雨-徑流與非點源污染物濃度及流失量關(guān)系分析[J].環(huán)境保護科學(xué),2010,36(2)：46-49.(LV Xingna.Analysis on relations of rainfall-runoff and nonpoint pollutionconcentrationandlossamount[J]. Environmental Protection Science,2010,36(2)：46-49. (in Chinese))

[8]白玉良,賀國靜,吳玉禾.大伙房水庫環(huán)境水文氣象效應(yīng)[J].水利管理技術(shù),1994(3)：35-37.(BAI Yuliang, HEGuojing,WUYuhe.Theenvironment hydrometeorological effects of Dahuofang reservoir[J]. Water Resources Development&Management,1994(3)：35-37.(in Chinese))

[9]張靜,何俊仕.大伙房水庫汛期氣象水文特征分析[J].人民黃河,2012,34(11)：45-47.(ZHANG Jing,HE Junshi.Study on hydro-meteorology characteristics in flood season of Dahuofang reservoir[J].Yellow River,2012,34 (11)：45-47.(in Chinese))

[10]盧雯.大伙房水庫上游地區(qū)農(nóng)業(yè)面源污染現(xiàn)狀及防治措施[J].環(huán)境保護與循環(huán)經(jīng)濟,2010(6)：72-73.(LU Wen.The status of agricultural nonpoint source pollution and thecontrolmeasuresoftheupstreamareaof Dahuofang reservoir[J].Environmental Protection and Circular Economy,2010(6)：72-73.(in Chinese))

[11]韓龍喜,陸東燕,李洪晶,等.高鹽度湖泊艾比湖風(fēng)生流三維數(shù)值模擬[J].水科學(xué)進展,2011,22(1)：97-103. (HAN Longxi,LU Dongyan,LI Hongjing,et al.3D numericalmodeling of wind-driven flow in the high salinity Aibi Lake[J].Advances in Water Science,2011,22(1)：97-103.(in Chinese))

[12]黃才安.明渠流速分布指數(shù)律與對數(shù)律的統(tǒng)一及轉(zhuǎn)換[J].人民長江,1994,25(1)：42-44.(HUANG Caian. The unification and conversion of open channel flow distribution exponential law and logarithmic law[J]. Yangtze River,1994,25(1)：42-44.(in Chinese))

[13]王利佳,鳳頎,白艷麗.大伙房水庫水質(zhì)污染和富營養(yǎng)化問題[J].遼寧城鄉(xiāng)環(huán)境科技,2001(6)：45-47. (WANG Lijia,FENG Qi,BAI Yanli.The water quality pollutionandeutrophicationproblemofDahuofang reservoir[J].Liaoning Urban and Rural Environmental Science&Technology,2001(6)：45-47.(in Chinese))

Water quality effects of agricultural non-point source total phosphorus input into reservoir in a typical rainstorm

SUN Yang,HAN Longxi
(College of Environment,Hohai University,Nanjing 210098,China)

According to the time-varying characteristics of agricultural non-point source total phosphorus(TP) input into the Dahuofang Reservoir,in the northeastern area of China,for a typical rainstorm case,we established a three-dimensional hydrodynamic mathematical model and the corresponding three-dimensional water quality model,taking into account the influence of the flow into the reservoir,in order to simulate and analyze the threedimensional spatial distribution and time-varying characteristics of TP in the reservoir area during a typical rainstorm storage process.The simulation results show that,under the influence of the inflow and the northeasterly wind direction,which was the dominant wind direction in the wet period,TP diffused towards the dam site along the southwest coastline,the TP-exceeded area increased gradually,and the concentration of TP decreased with the flow process.The simulation analysis of this process can provide a relevant technical basis for the forecast and control of non-point source pollution.

typical rainstorm;agricultural non-point source pollution;total phosphorus;three-dimensional numerical simulation;Dahuofang Reservoir

X524

A

1004-6933(2014)01-0018-05

201304-28 編輯：高渭文)

10.3969/j.issn.1004-6933.2014.01.004

水利部公益性行業(yè)科研專項(201101031)

孫楊(1989—),女,碩士研究生,研究方向為環(huán)境水力學(xué)。E-mail：sunyang817@126.com