养殖污水对水环境的影响预测

第24卷第2期海洋学研究Vol 24 No. 22006年6月JOURNAL OF MARINE SCIENCESune文章编号:1001-909X(2006)02-0039-10养殖污水对水环境的影响预测骆晓明(浙江省钱塘江管理局勘测设计院,浙江杭州310016)摘要:浙江三门湾宁海下洋涂围垦工程投入养殖后排放的养殖污水将对周边水环境产生的影响预测是该项目环境评价的重要內容。在海湾潮流动力模拟的基础上,建立了污染物对流扩散的数学模型,用以预测养殖污水排放后的扩散范围及强度大小,为环境评价提供科学依椐关键词:海涂围垦;养殖污水;污染物质量浓度;水质;对流-扩散方程;数值计箅中图分类号:X830.3文献标识码:A0引言浙江沿海在理论深度基准面以上的海涂资源总面积约为28.9万hm2。海涂围垦已经给浙江省带来了可观的经济效益,但也对环境产生了一些负面的影响,围垦使局部水域的水动力条件改变,出现海床的冲淤变化外,围垦区养殖污水对周边水环境也会产生一定的影响门。为此,本文从水动力数学模型出发,预测围垦工程投入养殖后排放的污水对周边水环境的影响情况。污染物扩散的模拟计算,在河流中的应用较早也较多,建立江河中污染物的动态数学模型,能较好地模拟污染物的散扩过程-η对海洋中污染物质迁移模式的研究也已有报道,一般都是用有限水体划分方法,建立具有统计学特征的污染物混合迁移模式·。水污染的计算机模拟适用于源或汇,对河道和海湾中污染物扩散的情况进行数值模拟,为环境的预测和评价提供科学依据本文以浙江三门湾宁海下洋涂围垦工程投入养殖后排放的养殖污水对周边水环境的影响为例,进行数值模拟,在海湾潮流动力模拟旳基础上,建立污染物扩散旳数学模型,以便预测养殖污水排放后的扩散范围及强度大小。1数学模型浙江三门湾是一个支港多、潮差大、垂直混合良好的海湾,在潮汐、径流作用下,其40海洋学研究24卷2期潮流表现为非恒定的流动。控制方程采用经垂线积分含水平涡动粘滞项的浅水潮波方1.1模型控制方程潮流的计算方程如下:az a(Hu)a(Hv)(1)Or tu ay tg arf乙)+(e,(2)+a如+by++8CH=a()×、3。(3)式中:ξ为水位;h为海底高程;H为总水位,H=h+;a,v分别为x,y方向上的垂线平均流速分量,m/s;g为重力加速度;∫为柯氏力参量(∫=2 CUsing,ω为地球自转角速率);C为谢才系数,取C=1H,n为糙率系数;,,分别为x,y方向的水平涡动扩散系数;W:,Wy为x,y方向的风应力分量,不考虑风应力作用时,取W=0,W,=0;t为时间式(1)为潮流连续方程,式(2)和(3)分别为x,y方向的动量守恒方程。考虑到三门湾边界及周边地形形状较为复杂,为了较好地模拟地形,对上述方程组的求解采用正交曲线坐标。1.2微分方程的离散格式计算模型采用平面二维曲线正交网格的有限差分模式。正交曲线坐标系下的控制方程与原方程相比,除了增加了一些系数之外,形式上是类似的1。对于上述方程,利用传统的ADⅠ法求解,其离散格式与矩形网格下的格式基本一致。2计算条件2.1计算范围为了便于给定边界条件,本模型南边界取壳塘山,东边界取石浦作为研究区域(图1)全水域面积约775km2。采用曲线网格对计算域进行剖分,与一般的矩形网格剖分相比,曲线网格可以更好地贴近边界,从而可以较好地模拟边界处的流态,减小边界所造成的计算影响。计箅域内剖分成300×30、总共有9000个网格,最大的网格边长取150m左右,围垦工程区附近的网格尺度控制在30m之內,计算时间步长为4min。模型计算网格见图1骆晓明:养殖污水对水环境的影响预测41121°2342.4122742.|"E29°2151.3石浦南H金28°4319.1图1数模计算网格图Fig. 1 Computational grid of the numerical model2.2参数的选取在数值模型中选取的参数分别为柯氏参量∫、粘滞系数A、水容重ρ和糙率系数n。其中: Coriolis参量∫=2 asing,φ为29.1°、a为7.29×10s-;A为50~60m2/s;P为1020kg/m3;n为0.028。2.3计算条件2.3.1初始条件初始条件取为(4)lu(x. y, t)lr=o=v(x, y, t)l=o=02.3.2边界条件开边界采用水位控制,即用潮位预报的方法得到开边界条件。计算区开边界采用潮位预报边界条件(图1)S=A+2H F cos[, t-(vo+u)+gI(5)式中:A为平均海面;F和(τ+a)为天文要素;H和g为调和常数。调和常数选取11个分潮计算,其中日分潮为4个(Q1、O1、P1和K1),半日分潮为4个(N2、M2、S2和K2),浅水分潮为3个(M1、MS4和M)。海洋学研究24卷2期3模型验证据本区域的潮位和潮流的大、小潮现场观测资料,对模型进行验证,从而评估模型的可靠性3.1潮位验证选择2003年4~5月水文测验期间,用获取的核电厂址、南田岛、牛山、花岙岛、巡检司、白玉湾岛临时潮位站同步的潮位观测资料进行验证,实测潮位与模拟计算的潮位之间拟合得较好,最高、最低潮位的模拟误差一般在10cm以内3.2潮流验证对工程区流速流向的验证采用2003年国家海洋局第二海洋研究所收集的有关资料,选取其中6个测点的潮流資料加以验证比较。单站实测结果与摸拟结果相比:大潮涨急、落急流向一般相差在10以内;大潮涨急、落急和涨潮平均、落潮平均流速相差在10%以内。总体而言,单站流向和流速的模拟结果令人满意。该模型可以用来预测工程投入养殖后排放污染物的扩散分布等。3.3潮流场分析大量研究表明,三门湾是一个强潮海湾,潮流及外海泥沙运动是海湾地貌发育的主导121°317.8121°54′6.3"E2915′19]纳潮下洋涂1000m2990′35.1图2工程实施后大潮涨急流场骆晓明:养殖污水对水环境的影响预测动力因素及环境糸件。计算域內单站潮流模拟验证的计算结果较好,基本反映了工程区海域潮流的实际变化。为进一步了解计算域內总体流场的分布,绘岀了计算域内涨急、落急的流矢分布示意图(图2和图3)。由图2和图3可见:(1)涨潮时,石浦港和珠门港两水道的潮流,主要影响三门湾內下涂洋一带的滩地和白礁水道。进入三门湾口的涨潮流,右侧要强于左侧,换言之,右侧满山水道的涨潮流要强于左侧的猫头水道。此外,满山水道的部分涨潮水通过青山门北侧深沟进入猫头水道。但通过猫头水道的落潮流泄出的潮量要比满山水道多些。(2)落潮时,各水道流速较大的特点都得到了很好地模拟。工程区临近水域的涨落潮流主要有猫头水道、满山水道、蛇盘水道和白礁水道4股。工程区被满山水道和白礁水道包围,以往复流为主。(3)三门湾内落潮流流速普遍要大于涨潮流流速,最大流速出现在口门,并往里有逐渐减小的趋势总之,计算域內流场模拟计算结果基本反映了该海堿潮流和潮波的实际变化,基本反映了三门湾潮流埸的总体特征。121037821°5463"E29°15′19.1m下洋涂子二石图3工程实施后大潮落急流场Fig 3 Flow field at ebb torrent time for spring tide under Scheme4养殖污水对水环境的影响预测4.1污染物质量浓度的对流-扩散计算模式在潮流计算的基础上,用二维非定态对流-扩散方程进行数值计算,预测污染物分布及海洋学研究24卷2期式中:c为污染物的质量浓度;K,K3分别为x,ν方向的湍流扩散系数;∫。为污染源的污染强度,∫=QC/△x△yH,Q为排水闸流量,C为质量浓度增量。(6)式的初始条件为c(x,y)|=0=c0(x,y)。该方程的边界条件:闭边界上,由于没有a c物质通量,取其质量浓度值为0;开边界上,当流向向外时,要求满足0…,当流向内流时,取边界上的质量浓度值为0。4.2污染物的源强4.2.1进、排水方式及排放口位置下洋涂围涂养殖区在涨潮时纳潮进水,海水通过西堤上的纳潮闸进入养殖区(图2);落潮时养殖污水通过东堤南端的排水闸排到围堤外的海域(图3)根据下洋涂围垦工程围区规划,围区养殖的估算面积为2082hm2,养殖净面积按估算面积的70%计算,约为1461.6hm2,养殖时日换水量为811万m3。工程投入养殖后养殖污水的排放口位于东堤排水闸处(图3)4.2.2污染物排放量氮、磷排放量围区内海水养殖基地是一个半人工控制的生态系,人工投饵是养殖动物的主要能量来源。投放的饵料一部分被养殖动物食用;另一部分沉入池底,池底残饵和动物排泄物等有机物经微生物分解后可产生氮、磷等营养物质,其中小部分被其它低等生物摄取,而大部分则在养殖池塘换水时流失进入周边海域。氮、磷排放量可由下式计算排放量一养殖净面积×饵料投放量×饵料中氮、磷的质量分数ⅹ(1-鱼虾贝对饵料的利用率)×流失率(7)式中:养殖净面积按估箅面积的70%计算,投饵量为33.3~53.3kg/hm2·a),按保守量计,取为53.3kg/hm2·a)。不同饵料中氮、磷的质量分数见表1表1不同饵料中氮、磷的质量分数Tab. 1 Mass ratio of n and P in different baits元素人工配合饵料活饵料0.015根据调査统计,养殖饲料中人工配合饵料一般占13%,鲜活饵料占87%,本评价以保守估算,全部按鲜活饵料计算。鱼虾贝对饵料的利用率一般在20%~35%1。本评价按30%计算,氮、磷的流失率分别以30%和5%计算按(7)式计算,可得本工程养殖区污水中氮的估算排放量为21924×800×2.8%×(1-30%)×30%kg/a,即10.31万kg/a;磷的估箅排放量为21924×800×1.2%1-30%)×5%kg/a,即0.74万kg/a本项目建成后围区内养殖塘日换水量为811万m3,每年养殖天数按300d计,由此可以估算岀养殖区因过剩饵料和动物粪便的分解而产生的养殖污水中氮的质量浓度增量为骆晓明:养殖污水对水环境的影响预测10.31×104×103×103/(81×104×103×300)mg/L,即0.042mg/L;磷的质量浓度增量为0.74×104×103×103/(811×104×103×300)mg/L,即0.003mg/L。化学需氧量(CODM)的排放量本工程按养殖污水最大排放量计,养殖区CODM的排放量为55.96万kg/a。根据以往对浙江省不同地域养殖池塘水质CODM的监测结果,养殖污水中CODλn的质量浓度平均增量约为0.23mg/L4.3水质模型的计算条件养殖污水经排水干渠汇集后由东堤排水闸排放入海,水闸最大排水量为188m3/s。养殖污水排放为非连续排放,只是当潮位处于毎个潮周期的平均潮位以下时才排放,时间约6h,其余时间水闸被关闭养殖污水中氮、磷的存在形态比较复杂,本评价近似地将所有流失的氮、磷分别视为无机氮和活性磷酸盐。无杋机氮、活性磷酸盐和O等因子的排放量、质量浓度、污染源强以及水质现状如表2所示。表2养殖污水中污染物的排放特征Tab 2 Pollutant discharging characteristic in cultivation wastewater项自水质预测评价因子无机活性磷酸盐COD排放量/万kg·a-1)质量浓度增量(mg·.-1)0.003水质质量浓度/大潮0.5741.840.310控制质量浓度标准(二类海水)/0.030充分考虑污染物在水体中长期排放的累积效应,在水质预测计算中,水流模型的设计潮型选取大、小潮连续进行计算。根据该工程污染源的特点和水环境污染的现状,水质影响的评价因子确定为:无机氮、活性磷酸盐和COD。由于围垦区附近海域无机氮和活性磷酸盐的质量浓度均超过二类海水的水质标准,水质已受到一定程度的污染(表2)。因此对水质影响的预测只计算养殖污水在围垦区附近水体中污染物的质量浓度增量,并以此来分析养殖污水对围垦区附近海域水环境的影响。4.4养殖污水对水环境影响的预测结果及分析应用水质模型模拟得到该围垦工程投入养殖后养殖污水排放的各污染因子质量浓度増量的分布。受文章篇幅所限,本文仅绘出COD质量浓度增量分布图(图4、图5)。从图4和图5可以清楚看见,在养殖污水的排放过程中,污染物主要集中在排水口附4.4.1COD质量浓度增量的分布大潮期间,养殖污水中COD的最大质量浓度增量为θ.10~0.15mg/,小潮期间为海洋学研究24卷2期120°35′20.8H2.3日:88排水闸纳潮闸纳潮闸纳潮闸图4大潮时养殖污水中COD质量浓度增量(mg/)分布Fig 4 COD distribution during spring tide120°3520.8121°508.6"E29°1319.3”纳潮闸排水纲纳潮闸纳潮闸纳潮闸图5小潮时养殖污水中COD质量浓度增量(mg/L)分布骆晓明:养殖污水对水环境的影响预测47污水中COD质量浓度的最大增量仅为本底值的8.2%,小潮期间为本底值的9.2%(表2、图4和图5)4.4.2无机氮、活性磷酸盐质量浓度增量的分布大潮期间,养殖污水中无杋氮的最大质量浓度增量为θ.020~0.030mg/L-,小潮期间为.020~0.026mg/L,分布范围均局限于东堤排水口附近。与环境本底值相比,大潮期间养殖污水中无机氮量浓度的最大增量仅为本底值的5.2%,小潮期间仅为本底值的大潮期间,养殖污水中活性磷酸盐的最大质量浓度增量为O.0012~0.0023mg/,小潮期间为0.o010~0.0015mg/Lλ,分布范围均局限于东堤排水口附近。与环境本底值相比大潮期间养殖污水中活性磷酸盐质量浓度的最大增量仅为本底值的6.4%,小潮期间仅为本底值的6.3%。4.4.3预测结果的分析以上计算结果可知,由于养殖污水的污染物质量浓度较低,与海域水质的本底值相比对水环境的“贡献”较小。因此,对海域水环境质量的影响也较小,养殖污水的排放不会改变邻近海域的水质类别;在养殖污水停止排放后,随着水体的流动及污染物质的扩散,污染物的质量浓度值将有所减小。5结语(1)在潮差大、垂直混合良好的海湾对潮流的模拟,采用含水平涡动粘滞项的经垂线积分的浅水潮波方程求解是合适的。(2)本文以三门湾下洋涂围涂工程投入养殖后排放的养殖污水对周边水环境的影响进行数值模拟。污染物的对流-扩散模型是建立在潮流模拟基础上的,并用二维非定态垂直平均扩散方程进行数值计算,计算污染物的质量浓度分布、预测养殖污水排放后的扩散范围及强度大小,可为环境保护研究提供依据。(3)本文以养殖污水中的无机氮、活性磷酸盐、COD污染物的源强和养殖进、排水方式的分折结果,作为污染物扩散数学模型的输入条件,计算结果表明,该工程投入养殖后排放的养殖污水的质量浓度较低,与海域水质的本底值相比,对水环境的“贡献”不大,因此,对海域水环境质量的影响也较小,养殖污水的排放不会改变海域的水质类别。参考文献Ⅰ]高爱根,杨俊毅,曾江宁,等.玉环坎门排污口邻近岩相潮间带生物分布特征[J.东海海洋,2004,22(4):24~30.2]胡礁星.黄浦江表层沉积物中有机氯农药的分布特征及风险评价「冂].环境科学学报,2005,26(3):44~48.3]Kinκ elbach w,候然杰,李患眀.河流中BOD-DO动态的数学模拟方法[J.环境科学学报,1981,1(2):166海洋学研究24卷2期[6]郭震运,王华东,刘培桐.铅山河金属污染物(Cu、Fe)迁移规律及污染预测研究[冂].环境科学学报,1983⑦]林玉环,汞污染河流底质迁移模式研究匚J.环境科学学报,1985,5(3):2768〗曾健业,吴瑜端.纳污港湾降解污染物质的有限水体混合迁移模式—厦门港污染物质的单纯混合模式匚J].环境科学学报,1985,5(4):395~403.]曾健业,吴瑜端,纳污港湾降解污染物质旳有限水体混合迁移模式∏厦门港Cn、φ迁移模式探讨匚冂.环境科学学报,1987,7(1):49~59,10]宋林松.水污染旳计算机模拟「J.环境科学学报,1983,3(2):141~-147[11]王益鸣,雪晓华,胡颢琰,等.浙江沿海产品中有机农的残留水平[J.东海海洋,2005,23(1):54~-64.12]杨晓兰,张键,叶新荣,等.南麂列岛自然保护区潮间带环境质量现状评价J].东海海洋,1994.12(2):70[13]陈增奇,陈飞星,李占玲,等.滨海湿地生态经济的综合评价模型[J.海洋学研究,2005,23(3):47~-514]冯士笮,孙文心.物理海洋数值计箅[M].郑州:河南科学技术出版社,1990[15]许卫忆.实际海域的赤潮生消过程数值模拟[J].海洋与湖沼,2001,32(6):598~60416]胡方西,曹沛錱.三门湾潮波运动特征及其与地貌发育的关系匚J.海洋与湖沼,1981,12(3):225~234[17]中国海湾志编纂委员会.中国海湾志第五分册[M].北京:海洋出版社,1992[18]谢忠明,鱼虾贝高产养殖技术[M].北京:农业出版社,199Prediction of impact of breeding wastewater dischargingon water environmenLUO Xiao-ming(Qiantang River Administration of Zhejiang Province, Hangzhou 310016, China)Abstract: When Xiayangtu Reclamation in Sanmenwan Bay, Ninghai County, Zhejianggoes into operation, the wastewater from animal breeding in the Reclamation Area will affect the water environment around. The prediction of the impact is very important part forenvironmental impact assessment. Based on the simulation of current in the bay, this paper establishes advection-diffusion equations for the pollutant, and then the model is applied to predict the areal distribution of pollutant due to the discharging of wastewaterfrom animal breeding, which provides scientific evident for environmental impact assessmentKey words reclamation; breeding wastewater; wastewater mass concentration; water quality; advection-diffusion equations numerical calculation