生态动力学模型及其在长江口及毗邻海域水污染防治规划中的应用邵一平－环境科学年会.pptxVIP

主要研究目的 建立长江口及毗邻海域水质模型，模拟控制污染物水质浓度场的时空变化规律，建立陆域污染源与水质浓度场的动态响应关系，为容量计算和总量控制方案的制定提供关键技术依据。 建立长江口三维生态动力学模型，模拟并校核二维线性水质模型计算形成的总量分配方案的合理性。 研究范围陆域包括江苏、浙江两省所属沿海沿江14个地市及上海市全市，面积共计10.36万km2。 海域包括三省市近岸海域，即长江口、杭州湾及毗邻海域，范围为32°00′～29°30′N，123°E以西海域，面积约3.8万km2。 水动力模型1. 长江口及其邻近海域三维流场计算；2. 长江口及其邻近海域三维盐度场计算；3. 长江口及其邻近海域三维温度场计算；4. 长江口及其邻近海域三维悬沙场计算；－水动力基本方程－－三维温度场和盐度场基本方程盐度为保守量，在动力场的作用下求解对流扩散方程 温度考虑表层太阳辐射的作用 －－三维悬沙场控制方程泥沙考虑底层的沉降和再悬浮以及盐度对沉降速度的影响 --模型范围网格数147×132，平均格距2km，最小格距在长江口内约300m，最大格距在外海约6km 进行三维计算时，模型垂向采用 坐标，共分14层。 －－初边条件选择水位和流速对外界动力响应较快，初值均取为零；盐度取洪季或枯季的典型分布值；温度初值取当月平均实测值；泥沙的初值取0。外海潮位采用M2、S2、K1、O1 、K2、N2和P1 等7个分潮调和常数。风场和温度场采用NOAA提供的NCEP/NCAR再分析风场和气象资料，每6小时一次，本区域精度约2度。对上边界长江径流取大通站的当月平均值，盐度、温度和泥沙的边界采用流入定常，流出无梯度的辐射边界条件。上游边界盐度取0，悬沙浓度取0.5 kg/m3；外海边界盐度取观测的典型值，悬沙浓度较小，取0.02 kg/m3；温度的开边界输入取实测月平均值。－－水动力模型率定和验证1996年9月长江口区域测站分布 碧海专项监测长江口二维线性水质模型基本方程－对流扩散方程－模型研究范围长江上游边界取在徐六泾；海域东边界为东经123°；南边界在北纬29°30’；北边界在江苏吕四港以北北纬32°15’长江口－杭州湾平面二维数学模型范围－网格剖分、模型边界网格剖分： 198×170。网格大小： 1500m*1500m。模型基面： 统一采用城建吴淞基面。水质边界：上游流入边界：（长江和钱塘江）取同步实测浓度过程线。外海边界：1）不同水期的实测浓度边界－（率定验证阶段） 2）零浓度梯度边界－（方案研究阶段）污染源负荷输入 33个概化点污染源 干湿沉降 海水养殖排污 富营养化控制指标（DIN，DIP） 有机污染物指标（CODMn，NH3-N）其他输入条件长江口风速、风向、气温、太阳辐射等外部输入资料同水动力模型保持一致，采用NOAA提供的资料。长江口毗邻海域同步水文水质监测站位图 模拟值与实测值的比较－－COD误差统计分析--模型系数取值污染物扩散系数 K2：0.03～0.17。本研究取0.1。 长江口污染物扩散系数取值范围：0～450m2/s。 －－污染物降解系数－－各排污源水质浓度响应场的模拟结果COD－－各排污源水质浓度响应场的模拟结果COD干湿沉降引起的背景浓度场COD氨氮 海水养殖产生的背景浓度场DIN DIP长江口及毗邻海域三维生态动力学模型氮、磷营养盐的迁移转化过程模型变量浮游植物（以C计）－Phytoplankton carbon (PC)浮游植物（以N计）－ Phytoplankton nitrogen (PN)浮游植物（以P计）－ Phytoplankton phosphorus (PP)叶绿素-a——Chlorophyll-a (CH)浮游动物——Zooplankton (ZC) 碎屑碳——Detritus carbon (DC)碎屑氮——Detritus nitrogen (DN)碎屑磷——Detritus phosphorus (DP) 无机氮——Inorganic nitrogen (IN－－NO3＋NH3)无机磷——Inorganic phosphorus (IP)溶解氧——Dissolved oxygen (DO)底栖植物（以C计）——Benthic vegetation carbon (BC) 北至北纬32°42‘，南至北纬28°，东至东经124°，西起江阴。 共1641个计算单元 垂向分六层，从上到下各层相对厚度分别为0.1H，0.15H，0.15H，0.15H，0.2H，0.25H。DO模拟结果和实测结果比较DIN模拟结果和实测结果比较DIP模拟结果和实测结果比较Chl-a—05-11Chl-a—05-7－－关键模型参数率定结果－－计算误差分析－－