養殖污水對水環(huán)境的影響預測

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