第24卷第2期海洋學研Vol.24 No. 22006年6月JOURNAL OF MARINE SCIENCESJune，2006文章編號: 1001-909X(2006)02 - 0039-10養(yǎng)殖污水對水環(huán)境的影響預測駱曉明(浙江省錢塘江管理局勘測設計院，浙江杭州310016)滴要:浙江三門灣寧海下洋涂圍墾工程投入養(yǎng)殖后排放的養(yǎng)殖污水將對周邊水環(huán)境產(chǎn)生的影響預測是該項目環(huán)境評價的重要內(nèi)容。在海灣潮流動力模擬的基礎(chǔ)上,建立了污染物對流擴散的數(shù)學模型，用以預測養(yǎng)殖污水排放后的擴散范圍及強度大小，為環(huán)境評價提供科學依椐。關(guān)鍵詞:海涂圍墾;養(yǎng)殖污水;污染物質(zhì)量濃度;水質(zhì);對流擴散方程;數(shù)值計算中圖分類號: X830. 3文獻標識碼: A0引言浙江沿海在理論深度基準面以上的海涂資源總面積約為28.9萬hm2D。海涂圍墾已經(jīng)給浙江省帶來了可觀的經(jīng)濟效益，但也對環(huán)境產(chǎn)生了一些負面的影響，圍墾使局部水域的水動力條件改變，出現(xiàn)海床的沖淤變化外，圍墾區(qū)養(yǎng)殖污水對周邊水環(huán)境也會產(chǎn)生一定的影響。為此，本文從水動力數(shù)學模型出發(fā)，預測圍墾工程投入養(yǎng)殖后排放的污水對周邊水環(huán)境的影響情況。污染物擴散的模擬計算，在河流中的應用較早也較多，建立江河中污染物的動態(tài)數(shù)學模型,能較好地模擬污染物的散擴過程2~7]。對海洋中污染物質(zhì)遷移模式的研究也已有報道，-般都是用有限水體劃分方法，建立具有統(tǒng)計學特征的污染物混合遷移模式(8.9]。水污染的計算機模擬適用于源或匯[10],對河道和海灣中污染物擴散的情況進行數(shù)值模擬，為環(huán)境的預測和評價提供科學依據(jù)11~13]。本文以浙江三門灣寧海下洋涂圍墾工程投入養(yǎng)殖后排放的養(yǎng)殖污水對周邊水環(huán)境的影響為例，進行數(shù)值模擬，在海灣潮流動力模擬的基礎(chǔ)上，建立污染物擴散的數(shù)學模型，以便預測養(yǎng)殖污水排放后的擴散范圍及強度大小。中國煤化工1數(shù)學模型MYHCNM HG浙江三門灣是一個支港多、潮差大、垂直混合良好的海灣，在潮汐、徑流作用下，其收稿日期: 2005-10-13作者簡介:駱曉明(1966-)， 男，浙江臨安市人，高級工程師，主要從事水利水電和海洋工程研究與設計。①浙江省水利水電勘測設計院.浙江省灘涂圍墾總體規(guī)劃報告，2005.●40●海洋學研究24卷2期潮流表現(xiàn)為非恒定的流動。控制方程采用經(jīng)垂線積分含水平渦動粘滯項的淺水潮波方程[14.15].1. 1模型控制方程潮流的計算方程如下:aZ」a(Hu)_ a(Hv)=0(1)atTdxdy)uduauxW,√u+v__.d。d、」a。u、Xt:+u十v+gx- fo-H+guC:H==Exax(2)Xxbxdy'h√u2+v2Jua+fu-(3)at「0 dyHC2Hdyey式中: ζ為水位; h為海底高程; H為總水位，H=h+ζ ; u, v分別為x，y方向上的垂線平均流速分量，m/s; g為重力加速度; f為柯氏力參量(f= 2wsing，w為地球自轉(zhuǎn)角速率); C:為謝才系數(shù)，取C.=-Ht，n為糙率系數(shù); e,， ey分別為x, y方向的水平渦動擴散系數(shù); W，W,為x，y方向的風應力分量，不考慮風應力作用時，取W.=0, W,=0; t為時間。式(1)為潮流連續(xù)方程，式(2)和(3)分別為x,y方向的動量守恒方程?？紤]到三門灣邊界及周邊地形形狀較為復雜，為了較好地模擬地形，對上述方程組的求解采用正交曲線坐標。1.2微分 方程的離散格式計算模型采用平面二維曲線正交網(wǎng)格的有限差分模式。正交曲線坐標系下的控制方程與原方程相比，除了增加了一些系數(shù)之外，形式上是類似的14.15]。對于上述方程，利用傳統(tǒng)的ADI法求解，其離散格式與矩形網(wǎng)格下的格式基本一致。2計算條件2.1計算范圍為了便于給定邊界條件，本模型南邊界取殼塘山，東邊界取石浦作為研究區(qū)域(圖1)。全水域面積約775 km'。采用曲線網(wǎng)格對計算域進行剖分，與中國煤化工，曲線網(wǎng)格可以更好地.貼近邊界，從而可以較好地模擬邊界處的MYHCNMHG計算影響。計算域內(nèi)剖分成300X300、總共有90000個網(wǎng)格，最大的網(wǎng)格邊長取150m左右，圍墾工程區(qū)附近的網(wǎng)格尺度控制在30 m之內(nèi)，計算時間步長為4 min。模型計算網(wǎng)格見圖1。駱曉明:養(yǎng)殖污水對水環(huán)境的影響預測●41●121"23'42. 4”12297'42.1"E29921'51.3”N石浦災洋島南田健跳28° 43'19. 1 "圖1數(shù)模計算網(wǎng)格圖Fig. 1 Computational grid of the numerical model2.2參數(shù)的選取在數(shù)值模型中選取的參數(shù)分別為柯氏參量f、粘滯系數(shù)A、水容重ρ和糙率系數(shù)n。其中: Coriolis參量f= 2wsinp,φ為29.1°、w為7. 29X10~*s-';A為50~ 60 m2/s;p為1 020kg/m'; n為0.028。2.3計算條件2.3.1 初始條件(5(x，y，t)|2=o=ζ。初始條件取為:lu(x， y, t)|r=o=v(x， y，t)|r=o=0(4)2.3.2 邊界條件開邊界采用水 位控制，即用潮位預報的方法得到開邊界條件。計算區(qū)開邊界采用潮位預報邊界條件(圖1):ζ=Ao2rHF;MYH中國煤化工(5)式中: A。為平均海面; F; 和(vo+u);為天C NMH G常數(shù)。調(diào)和常數(shù)選取11個分潮計算，其中日分潮為4個(Q、 O、P1和K)，半日分潮為4個(N2、 M2、Sz和K2)，淺水分潮為3個(M、MS和M;)。閉邊界在閉邊界， 取法向?qū)?shù)為0，即Vn=0;在潮灘區(qū)采用漫灘邊界處理?；鶞拭鎸⑺钯Y料的深度基準面統(tǒng)一到平均海平面上。●42●海洋學研究24卷2期3模型驗證據(jù)本區(qū)域的潮位和潮流的大、小潮現(xiàn)場觀測資料，對模型進行驗證，從而評估模型的可靠性。3.1潮位驗證選擇2003年4~5月水文測驗期間，用獲取的核電廠址、南田島、牛山、花岙島、巡檢司、白玉灣島臨時潮位站同步的潮位觀測資料進行驗證，實測潮位與模擬計算的潮位之間擬合得較好，最高、最低潮位的模擬誤差一般在10 cm以內(nèi)。3.2潮流驗證對工程區(qū)流速流向的驗證采用2003年國家海洋局第二海洋研究所收集的有關(guān)資料,選取其中6個測點的潮流資料加以驗證比較。單站實測結(jié)果與摸擬結(jié)果相比:大潮漲急、落急流向一般相差在10°以內(nèi);大潮漲急、落急和漲潮平均、落潮平均流速相差在10%以內(nèi)?？傮w而言，單站流向和流速的模擬結(jié)果令人滿意。該模型可以用來預測工程投入養(yǎng)殖后排放污染物的擴散分布等②。3.3潮流場分析大量研究表明，三門灣是一個強潮海灣，潮流及外海泥沙運動是海灣地貌發(fā)育的主導121*31'7.8"12154'6.3"E729*15'19.]"N青山港、納潮閘納下洋涂潮閘。以滿山樂遭，1.00m/m中國煤化工CHCNM H GiI 290'35.1”圖2工程實 施后大潮漲急流場Fig.2 Flow field at flood torrent time for spring tide under Scheme②國家海洋局第二海洋研究所，寧?？h下洋涂圍墾工程對周邊水流及泥沙沖淤影響的研究報告，2005.駱曉明:養(yǎng)殖污水對水環(huán)境的影響預測●43●動力因素及環(huán)境條件[16.17]。計算域內(nèi)單站潮流模擬驗證的計算結(jié)果較好，基本反映了工程區(qū)海域潮流的實際變化。為進一步了解計算域內(nèi)總體流場的分布，繪出了計算域內(nèi)漲急、落急的流矢分布示意圖(圖2和圖3)。由圖2和圖3可見: (1)漲潮時，石浦港和珠門港兩水道的潮流，主要影響三門灣內(nèi)下涂洋-帶的灘地和白礁水道。進入三門灣口的漲潮流，右側(cè)要強于左側(cè)，換言之，右側(cè)滿山水道的漲潮流要強于左側(cè)的貓頭水道。此外，滿山水道的部分漲潮水通過青山門北側(cè)深溝進入貓頭水道。但通過貓頭水道的落潮流泄出的潮量要比滿山水道多些。(2)落潮時，各水道流速較大的特點都得到了很好地模擬。工程區(qū)臨近水域的漲落潮流主要有貓頭水道、滿山水道、蛇盤水道和白礁水道4股。工程區(qū)被滿山水道和白礁水道包圍，以往復流為主。(3) 三門灣內(nèi)落潮流流速普遍要大于漲潮流流速，最大流速出現(xiàn)在口門，并往里有逐漸減小的趨勢。總之，計算域內(nèi)流場模擬計算結(jié)果基本反映了該海域潮流和潮波的實際變化，基本反映了三門灣潮流埸的總體特征。121931'7.8"121954'6.3"E29915' 19.1”N黃山港下洋涂水縮永道遇山水道行強1.000 ws290 0'35.1"圖3工程實施后大潮落急流場Fig.3 Flow field at ebb torrent time for spring tide under Scheme4養(yǎng)殖污水對水環(huán)境的影響預測中國煤化工,MYHCNMHG4.1污染物質(zhì)量濃度的對流~擴散計算模式在潮流計算的基礎(chǔ)上，用二維非定態(tài)對流-擴散方程進行數(shù)值計算，預測污染物分布及影響范圍。污染物質(zhì)量濃度的二維非定態(tài)的垂直平均對流-擴散方程為an+u)caxdc. _ 1rd(HK.Taadc、二) +:(HK,5)]+f.a(6)dy H-axay"dy●44●海洋學研究24卷2期式中:c為污染物的質(zhì)量濃度;K;,K,分別為x,y方向的湍流擴散系數(shù);f.為污染源的污染強度，f.=QgCo/Ox△yH, Q。為排水閘流量，C, 為質(zhì)量濃度增量。(6)式的初始條件為c(x，y)|=o=co(x, y)。該方程的邊界條件:閉邊界上，由于沒有物質(zhì)通量，取其質(zhì)量濃度值為0;開邊界上，當流向向外時，要求滿足Xc+u o+vaac=0，當Xtaa流向內(nèi)流時，取邊界上的質(zhì)量濃度值為0。4.2污染物的源強4.2.1 進、排水方式及排放口位置.下洋涂圍涂養(yǎng)殖區(qū)在漲潮時納潮進水,海水通過西堤上的納潮閘進入養(yǎng)殖區(qū)(圖2);落潮時養(yǎng)殖污水通過東堤南端的排水閘排到圍堤外的海域(圖3)。根據(jù)下洋涂圍墾工程圍區(qū)規(guī)劃，圍區(qū)養(yǎng)殖的估算面積為2082 hm',養(yǎng)殖凈面積按估算面積的70%計算，約為1461.6hm’，養(yǎng)殖時日換水量為811萬m°。工程投入養(yǎng)殖后養(yǎng)殖污水的排放口位于東堤排水閘處(圖3)。4.2.2 污染物排放量氮、磷排放量圍區(qū) 內(nèi)海水養(yǎng)殖基地是一個半人工控制的生態(tài)系，人工投餌是養(yǎng)殖動物的主要能量來源。投放的餌料一部分被養(yǎng)殖動物食用，另-部分沉入池底，池底殘餌和動物排泄物等有機物經(jīng)微生物分解后可產(chǎn)生氮、磷等營養(yǎng)物質(zhì)，其中小部分被其它低等生物攝取，而大部分則在養(yǎng)殖池塘換水時流失進入周邊海域[18]。氮、磷排放量可由下式計算:排放量=養(yǎng)殖凈面積X餌料投放量X餌料中氮、磷的質(zhì)量分數(shù)X(1一魚蝦貝對餌料的利用率) X流失率(7)式中:養(yǎng)殖凈面積按估算面積的70%計算，投餌量為33.3~53.3 kg/(hm2●a),按保守量計，取為53.3 kg/(hm?●a)圍。不同餌料中氮、磷的質(zhì)量分數(shù)見表1。表1不同餌料中氮、 磷的質(zhì)量分數(shù)Tab.1Mass ratio of N and P in different baits元素人工配合餌料鮮活餌料0.0150. 028P0.004 .0. 012根據(jù)調(diào)查統(tǒng)計，養(yǎng)殖飼料中人工配合餌料一般占13 %，鮮活餌料占87 %，本評價以保守估算，全部按鮮活餌料計算。魚蝦貝對餌料的利用率-般在20%~35%[16]。本評價按30%計算，氮、磷的流失率分別以30%和5%計算@。安(7) 式計算，可得本工程養(yǎng)殖區(qū)污水中氮的估算排放量為21 924X 800X2.8%X(1-30%) X30% kg/a，即10.31萬kg中國煤化工21 924X 800X1.2%x(1-30%) X5% kg/a ，即0.74萬kg/al.YHCNMHG本項目建成后圍區(qū)內(nèi)養(yǎng)殖塘日換水量為811萬m'，每年養(yǎng)殖天數(shù)按300d計，由此可以估算出養(yǎng)殖區(qū)因過剩餌料和動物糞便的分解而產(chǎn)生的養(yǎng)殖污水中氮的質(zhì)量濃度增量為③岱山縣科協(xié)，海水池塘養(yǎng)殖實用技術(shù)，2001.④國家海洋局第二海洋研究所.溫州龍灣海濱圍墾工程環(huán)境影響報告書，2004.駱曉明:養(yǎng)殖污水對水環(huán)境的影響預測●45.10.31X10*X10*X10*/ (811X10*X10*X300) mg/L，即0. 042 mg/L;磷的質(zhì)量濃度增量為0.74X10*X10*X10*/ (811X10*X10*X300) mg/L，即0.003 mg/L?；瘜W需氧量(CODmn) 的排放量本工程 按養(yǎng)殖污水最大排放量計，養(yǎng)殖區(qū)CODMn的排放量為55.96萬kg/a。根據(jù)以往對浙江省不同地域養(yǎng)殖池塘水質(zhì)CODwn的監(jiān)測結(jié)果，養(yǎng)殖污水中CODMn的質(zhì)量濃度平均增量約為0.23 mg/L⑧。4.3水質(zhì)模型的計算條件養(yǎng)殖污水經(jīng)排水干渠匯集后由東堤排水閘排放入海，水閘最大排水量為188m2/s。養(yǎng)殖污水排放為非連續(xù)排放，只是當潮位處于每個潮周期的平均潮位以下時才排放，時間約6h，其余時間水閘被關(guān)閉。養(yǎng)殖污水中氮、磷的存在形態(tài)比較復雜，本評價近似地將所有流失的氮、磷分別視為無機氮和活性磷酸鹽。無機氮、活性磷酸鹽和COD等因子的排放量、質(zhì)量濃度、污染源強以及水質(zhì)現(xiàn)狀如表2所示，表2養(yǎng)殖污水中污染物的排放特征Tab.2 Pollutant discharging characteristic in cultivation wastewater水質(zhì)預測評價因子項目無機氮活性磷酸鹽COD排放量/(萬kg. a-1)10. 310.7455.96質(zhì)量濃度增量/(mg. L-1)0.0420.0030. 23水質(zhì)質(zhì)量濃度/大潮0. 5740.0361.84.(mg●L-1)小潮0.3100. 0241.42控制質(zhì)量濃度標準(二類海水)/0. 300.0303.0充分考慮污染物在水體中長期排放的累積效應，在水質(zhì)預測計算中，水流模型的設計潮型:選取大、小潮連續(xù)進行計算。根據(jù)該工程污染源的特點和水環(huán)境污染的現(xiàn)狀，水質(zhì)影響的評價因子確定為:無機氮、活性磷酸鹽和COD。由于圍墾區(qū)附近海域無機氮和活性磷酸鹽的質(zhì)量濃度均超過二類海水的水質(zhì)標準，水質(zhì)已受到一定程度的污染(表2)。因此對水質(zhì)影響的預測只計算養(yǎng)殖污水在圍墾區(qū)附近水體中污染物的質(zhì)量濃度增量，并以此來分析養(yǎng)殖污水對圍墾區(qū)附近海域水環(huán)境的影響。4.4養(yǎng)殖污水對水環(huán)境影響的預測結(jié)果及分析應用水質(zhì)模型模擬得到該圍墾工程投入養(yǎng)殖后養(yǎng)殖污水排放的各污染因子質(zhì)量濃度增量的分布。受文章篇幅所限，本文僅繪出COD質(zhì)量濃度增量分布圖(圖4、圖5)。從圖4和圖5可以清楚看見，在養(yǎng)殖中國煤化工:物主要集中在排水口附近。YHCNMHG4.4.1 COD質(zhì)量濃度增量的分布大潮期間，養(yǎng)殖污水中COD的最大質(zhì)量濃度增量為0. 10~0. 15 mg/L，小潮期間為0.10~0.13mg/L,分布范圍均局限于東堤排水口附近。與環(huán)境本底值相比，大潮期間養(yǎng)殖●46●海洋學研究24卷2期120*35'20. 8"121*50'8.6" E29*13'19.3"N納潮閘e是持水閘29*3'45. 5"圖4大潮時 養(yǎng)殖污水中COD質(zhì)量濃度增量(mg/L) 分布Fig. 4 COD distribution during spring tide120"35'20. 8"T 1 29*13'19.3"0.1納湖閘排水閘中國煤化工)MHCNMHG29"3'45. 5*圖5小潮時 養(yǎng)殖污水中COD質(zhì)量濃度增量(mg/L) 分布Fig.5 COD distribution during neap tide駱曉明:養(yǎng)殖污水對水環(huán)境的影響預測●47.污水中COD質(zhì)量濃度的最大增量僅為本底值的8.2%，小潮期間為本底值的9.2%(表2、圖4和圖5)。4.4.2無機氮、活性磷酸鹽質(zhì)量濃度增量的分布大潮期間，養(yǎng)殖污水中無機氮的最大質(zhì)量濃度增量為0.020~0.030mg/L，小潮期間為0.020~0.026 mg/L，分布范圍均局限于東堤排水口附近。與環(huán)境本底值相比，大潮期間養(yǎng)殖污水中無機氮量濃度的最大增量僅為本底值的5.2%，小潮期間僅為本底值的8.3%。大潮期間，養(yǎng)殖污水中活性磷酸鹽的最大質(zhì)量濃度增量為0.001 2~0. 002 3 mg/L,小潮期間為0.0010~0.0015mg/L,分布范圍均局限于東堤排水口附近。與環(huán)境本底值相比，大潮期間養(yǎng)殖污水中活性磷酸鹽質(zhì)量濃度的最大增量僅為本底值的6.4%，小潮期間僅為本底值的6.3%。4.4.3預測結(jié)果的分析以上計算結(jié)果可知，由于養(yǎng)殖污水的污染物質(zhì)量濃度較低，與海域水質(zhì)的本底值相比，對水環(huán)境的“貢獻”較小。因此，對海域水環(huán)境質(zhì)量的影響也較小，養(yǎng)殖污水的排放不會改變鄰近海域的水質(zhì)類別;在養(yǎng)殖污水停止排放后,隨著水體的流動及污染物質(zhì)的擴散,污染物的質(zhì)量濃度值將有所減小。5結(jié)語(1)在潮差大、垂直混合良好的海灣對潮流的模擬，采用含水平渦動粘滯項的經(jīng)垂線積分的淺水潮波方程求解是合適的。(2)本文以三門灣下洋涂圍涂工程投入養(yǎng)殖后排放的養(yǎng)殖污水對周邊水環(huán)境的影響進行數(shù)值模擬。污染物的對流-擴散模型是建立在潮流模擬基礎(chǔ)上的，并用二維非定態(tài)垂直平均擴散方程進行數(shù)值計算，計算污染物的質(zhì)量濃度分布、預測養(yǎng)殖污水排放后的擴散范圍及強度大小，可為環(huán)境保護研究提供依據(jù)。(3)本文以養(yǎng)殖污水中的無機氮、活性磷酸鹽、COD污染物的源強和養(yǎng)殖進、排水方式的分折結(jié)果，作為污染物擴散數(shù)學模型的輸入條件，計算結(jié)果表明，該工程投入養(yǎng)殖后排放的養(yǎng)殖污水的質(zhì)量濃度較低，與海域水質(zhì)的本底值相比，對水環(huán)境的“貢獻”不大，因此，對海域水環(huán)境質(zhì)量的影響也較小，養(yǎng)殖污水的排放不會改變海域的水質(zhì)類別。參考文獻:中國煤化工[1] 高愛根，楊俊毅，曾江寧，等.玉環(huán)坎門排污C.MHCNMHG東海海洋，2004. 22 (4):24~ 30.[2]胡礁星、黃浦江表層沉積物中有機氯農(nóng)藥的分布特征及風險評價J]. 環(huán)境科學學報，2005， 26 (3): 44~48.[3]Kinzelbach W，候然杰，李惠明.河流中BOD DO動態(tài)的數(shù)學模擬方法[J]. 環(huán)境科學學報，1981，1 (2): 166~ 170.[4]候然杰，李惠明.描述北方河流封凍期BOD-DO動態(tài)的數(shù)學模型[J].環(huán)境科學學報，1982, 2 (2): 113~119.[5]雷孝思，鄧玉珍，一個適用于中距離范圍的湍流擴散模型[J]. 環(huán)境科學學報，1984， 4 (4): 376~ 385. .●48●海洋學研究24卷2期[6]郭震運，王華東，劉培桐.鉛山河金屬污染物(Cu、 Fe)遷移規(guī)律及污染預測研究[J]. 環(huán)境科學學報，1983.3 (4): 298~ 309.[7]林玉環(huán).汞污染河流底質(zhì)遷移模式研究[J]. 環(huán)境科學學報，1985，5 (3): 276~285. .[8]曾健業(yè)，吳瑜端，納污港灣降解污染物質(zhì)的有限水體混合遷移模式-廈門港污染物質(zhì)的單純混合模式 [J].環(huán)境科學學報，1985， 5 (4): 395~ 403.[9]曾健業(yè)，吳瑜端.納污港灣降解污染物質(zhì)的有限水體混合遷移模式-I1廈門港Cn、Pn遷移模式探討[J].環(huán)境科學學報，1987， 7 (1): 49~59.[10]宋林松、水污染的計算機模擬[J].環(huán)境科學學報，1983， 3 (2): 141~147.[11]王益鳴，雪曉華，胡顥琰，等.浙江沿海產(chǎn)品中有機農(nóng)的殘留水平[J]. 東海海洋，2005， 23 (1): 54~64.[12] 楊曉蘭，張鍵，葉新榮，等.南麂列島自然保護區(qū)潮間帶環(huán)境質(zhì)量現(xiàn)狀評價U].東海海洋，1994，12 (2): 70~ 76.[13]陳增奇，陳飛星，李占玲，等.濱海濕地生態(tài)經(jīng)濟的綜合評價模型[J]. 海洋學研究，2005， 23 (3): 47~55.[14]馮士笮，孫文心物理海洋數(shù)值計算[M]. 鄭州:河南科學技術(shù)出版社，1990.[15]許衛(wèi)憶.實際海域的赤潮生消過程數(shù)值模擬J]. 海洋與湖沼，2001. 32 (6): 598~604.[16] 胡方西，曹沛奎.三門灣潮波運動特征及其與地貌發(fā)育的關(guān)系[J]. 海洋與湖沼，1981， 12 (3): 225~234.[17]中國海灣志編纂委員會，中國海灣志第五分冊[M]. 北京:海洋出版社，1992.謝忠明，魚蝦貝高產(chǎn)養(yǎng)殖技術(shù)[M].北京:農(nóng)業(yè)出版社，1995. .Prediction of impact of breeding wastewater dischargingon water environmentLUO Xiao- ming(Qiantang River Administration of Zhejiang Province, Hangzhou 310016，China)Abstract: When Xiayangtu Reclamation in Sanmenwan Bay, Ninghai County, Zhejiang,goes into operation, the wastewater from animal breeding in the Reclamation Area will af-fect 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 pa-per establishes advection -diffusion equations for the pollutant, and then the model is ap-plied to predict the areal distribution of pollutant due to the discharging of wastewaterfrom animal breeding， which provides scientific evident for environmental impact assess-ment.Key words : reclamation; breeding wastewatey- atenno m.. r-ncentration ; water qual-中國煤化工ity; advection diffusion equations ; numeriTYHCNMHG