DB21∕T 2426-2015 污染事故影响数值预测分析技术规程

ICS13.020.40Z01DB21辽宁省地方标准DB21/T2426—2015污染事故影响数值预测分析技术规程ImpactAnalysisandPredictionTechniqueforpollutionaccidentbasedonthenumericalmodel2015-2-17发布2015-4-17实施术监督局发布DB21/T2426-2015I前言本标准按照GB/T1.1-2009标准化工作导则起草。本标准由辽宁省海洋水产科学研究院提出。本标准由辽宁省海洋与渔业厅归口。本标准主要起草单位：辽宁省海洋水产科学研究院。本标准主要起草人：王昆、王年斌、宋伦、吴金浩、马志强、张玉凤、宋永刚、田金、李楠、杨爽。本标准首次发布。DB21/T2426-20151污染事故影响数值预测分析技术规程1范围本标准规定了基于数值模型的污染事故影响预测分析技术规程的术语、操作程序与技术方法。本标准适用于辽宁省管辖海域的污染事故影响预测分析。2规范性引用文件下列文件对于本文件的应用是必不可少的。凡是注日期的引用文件，仅所注日期的版本适用于本文件。凡是不注日期的引用文件，其最新版本（包括所有的修改单）适用于本文件。GB3097海水水质标准GB11607渔业水质标准GB/T12763.2海洋调查规范第2部分：海洋水文观测GB17378.4海洋监测规范第4部分：海水分析3术语下列术语适用于本标准。3.1水动力数值模拟将水动力学的实际问题抽象成满足其物理特征的典型偏微分方程组来描述的数学问题，并采用合适的离散方法对方程组进行数值求解的过程。3.2污染事故影响预测分析对特定海域已经发生的污染事故进行后期影响的预测、评估与分析的过程。3.3污染物输运模型在水动力学模拟的背景流场的基础上，建立以特征污染物浓度为主要未知量的对流扩散模型，模拟预测污染物入海后的浓度场实时分布状况。DB21/T2426-201523.4参数率定利用已有历史数据或补测的海流、潮位以及化学检测分析浓度资料与模拟值进行比较，对数值模型中的参数进行调节确定的过程。3.5源强由事故引起的单位时间内污染物的排放量。4方法与技术针对不同海洋功能区海水水质标准GB3097与渔业水质标准GB11607的要求，计算受事故影响的主要特征污染物浓度的超标范围。4.1站位的布设原则1)站位布设应均匀覆盖事故影响的研究海域，所获信息应代表和反映海域的基本环境特征；2)尊重历史，考虑历史站位的延续性；3)充分考虑海域周边的陆域环境特征和压力，对受污染事故影响和人类活动影响显著的区域需加密布设。4)充分考虑海域本身的环境特征和问题，如水动力环境特征、地形地貌特征等；5)关注沿海地区经济社会发展需求，在开发热点区、敏感区及邻近海域布设监测站位；6)以尽可能少的站位来获取尽可能多的环境信息，但一般不得少于6个。4.2监测指标与方法a)水动力要素：海流与水深,监测方法按照GB/T12763.2规定执行；b)水质要素：营养盐浓度、事故发生处特征污染物指标的浓度,检测分析方法按照GB17378.4规定执行。4.3监测频率和时间a)水动力要素：分四个季节，分别在大、小潮期每隔1h进行至少连续25h的准同步观测在转流时刻前后可每隔0.5h进行一次观测；b)水质要素：与水动力要素观测同步采样。4.4数值模型4.4.1水动力数值方程可采用笛卡尔坐标系下，不可压缩流体的连续性方程和动量方程：0=∂∂+∂∂+∂∂zwyvxuDB21/T2426-201532222()()hVduquuufvgdtxxxyzzhgg∂∂∂∂∂∂=--+++∂∂∂∂∂∂2222()()hVdvqvvvfugdtyyxyzzhgg∂∂∂∂∂∂=---+++∂∂∂∂∂∂2222()()hVdwq∂∂∂∂∂=+++∂∂∂∂∂在水深方向上对连续性方程（1）进行积分得自由面方程：0=⎥⎥⎦⎤⎢⎢⎣⎡∂∂+⎥⎥⎦⎤⎢⎢⎣⎡∂∂+∂∂∫∫--hhhhhvdzyudzxt式中，h(x,y)为静水深；t为时间；u、v、w分别为x、y、z方向的流速；η为水位；q为水体动压值；f为柯氏力系数；γh、γv为水平和垂直向的紊动粘性系数，推荐采用ke-双方程紊流模型来确定。紊流动能输运k方程：esnsnsn-+∂∂∂∂+∂∂∂∂+∂∂∂∂=Gzkzykyxkxdtdkktktkt)()()(紊流特征量输运ε方程：kCGkCzzyyxxdtdttt221)()()(ε-ε+∂ε∂σν∂∂+∂ε∂σν∂∂+∂ε∂σν∂∂=εεεεεε式中，k为紊动动能；e为紊动耗散率；ks、es分别为紊动动能k和紊动耗散率e对应的Prandtl数；e1C、e2C是经验系数；G为紊动产生项。4.4.1.1初始条件u(x,y,t0)=u0(x,y)；v(x,y,t0)=v0(x,y)；η(x,y,t0)=η0(x,y)其中，u0、v0、η0分别为初始流速和潮位，在海域流场计算中，初始条件与计算的最终结果无关。4.4.1.2边界条件固壁边界：给定法向流速为零的有滑移边界。在外海开边界处，给定潮位时间序列文件。可通过以下两种途径来获得：1)利用潮汐表直接查表换算；2)利用潮汐调和分析的方法，即开边界的水位为时间的调和函数：∑++-=kkkkkkkkUVgtHf))(cos(地wh式中，潮汐调和常数振幅与迟角可根据实测数据利用最小二乘法求出。有河流入海的情况，需要考虑感潮河段的流量，以流量边界的形式给出。DB21/T2426-20154在自由表面（z=η）处：)(uuzuaTv-=∂∂gn；)(vvzvaTv-=∂∂gn式中，Tg是风应力系数，au、av分别为风速在x、y方向的分量。在底面（z=-h）处：uzuBvgn=∂∂；vzvBvgn=∂∂式中，Bg是底摩阻系数，根据湍流边界层理论：222nBcvug+=g式中，nc是谢才系数，由曼宁公式1/61ncHn=计算，h+=hH为总水深，n为糙率，其值可取范围0.020~0.028。4.4.2污染物输运模拟方程以水动力学方程模拟出的流场为背景，为了查清事故污染物的影响程度与范围，利用对流扩散输运模型进行数值模拟，三维对流扩散方程为：MzCkxyCkyxCkxzCwwyvCxuCtCvhhs+∂∂∂∂+∂∂∂∂+∂∂∂∂=∂-∂+∂∂+∂∂+∂∂)()()()()()((5)式中，C为浓度；M为源项；kh、kv分别为水平和垂直向的污染物扩散系数；ws为污染物的沉降速度，其值主要与海区及污染物的粒径有关。4.4.2.1数值计算输入条件入流边界按实测给定。出流边界：0=∂∂+∂∂nCutCn式中，nu为边界的法向流速。DB21/T2426-20155在自由表面（z=η）处：)()(CCCwwzCkTTTsv-+=--∂∂ba式中，αT、βT为表面处与C无关的物质交换系数，若自由表面处没有物质交换，则αT=βT=0。在底面（z=-h）处：)()(CCCwwzCkBBBsv-+=-+∂∂-ba式中，αB、βB为底面处与C无关的物质交换系数，若底面处没有物质交换，则αB=βB=0。4.4.2.2源强的确定污染事故引起的主要排污口产生的点源源强根据排污口处特征污染物的浓度和排污量来确定，对于线源可概化成多个点源来进行等效计算。4.4.3网格划分对于计算域比较规则，地形变化比较平缓的区域，宜采用结构化四边形网格进行划分；对于边界岸线变化复杂、地形坡道变化比较大的计算区域，宜采用非结构化三角形网格对计算域进行划分，能够满足尽量真实的刻画计算域、与实际地形情况做到最大程度的吻合的要求；对于受事故影响较大的排污口位于湾顶及边界附近，需要在边界处对网格进行加密，提高精度和分辨细部也是主要靠网格的合理布置与加密来实现。4.4.4离散格式水动力学潮流场的水平和垂向数值模拟，推荐采用在垂向分层的空间交错棱柱形网格的基础上，采用有限差分与有限体积相结合的方法来离散求解。如果是河道、海湾以及浅海区域，可采用只分一层，即简化成二维模型来求解。输运方程的离散，采用半隐的一阶迎风有限体积法离散，同时为了避免施加额外的稳定性条件，这里采用子迭代方法在每个棱柱上对方程进行离散。4.4.5动边界的模拟对于潮滩，水陆交界的位置随着潮位的涨落时而淹没时而干出，动边界内网格节点的干湿变化宜采用限制水深法来进行处理，这里限制水深取0.3m。4.4.6参数率定与模型验证为了证明模型的准确性与稳定性，需利用4.1~4.3的实测资料或已有历史数据对模型进行验证与参数率定。取计算稳定后的典型潮时(如涨急、高潮、落急、低潮时)的潮位、海流或浓度的模拟值与实测值进行比较，获得能满足该海域实际物理特征的数值模型，用于后期的模拟与预报。DB21/T2426-201564.5计算结果与分析通过上述模型的模拟计算，取浓度分布趋于稳定时一个典型潮周期内的的平衡浓度场，为了使结果表现更加地直观化、具体化与形象化，利用面向对象的编程技术或借助软件中提供的后处理程序工具将潮位、潮流、扩散浓度场等模拟结果进行多时间步动画演示，在说明了浓度场分布扩散趋势与流场二者之间相互关系的同时，也进行了污染物输运扩散路径与范围的预测分析，并附以相应的图表。4.6超标范围确定根据辽宁省海洋功能区划，对不同的海洋功能区，按GB3097、GB11607标准规定的限值，利用污染物扩散影响浓度分布场进行事故影响的污染面积计算，即超标准水质的包络线形成的面积，采用相关软件内部嵌套的命令流来求得，并附以相应的图表。注意事项：标准中的数值模型是经过多个实际工程以及有解析解的算例验证得出的。附录A（资料性附录）污染事故影响数值预测分析技术规程计划框架图，见图A.1.图A.1污染事故影响数值预测分析技术规程计划框架图_________________________________