大连湾易溶保守液体化学品污染扩散初探研究.pdfVIP

大连湾易溶保守液体化学品污染扩散初探 雷孝平 朱国林 (大连海事大学 中国大连) 摘要 本文根据船运液体化学品的性质，将其溢漏后运动形式分成五类 以大连湾为依托 对其中易溶保守液体化学品的扩散形式进厅初步研究。根据浅水潮波理论，采用AD[法，计 算大连湾m:分潮流场。利用即得潮流的潮位资料，采用深度平均二维扩散模型，对易溶保 守液体化学品进行扩散计算，取得直观上较为准确的污染物浓度时空分布。为环境应急反应 决策，污染损失计算提供了依据 ，引言 近年来大连湾海上散装液体化学品运输业发展迅猛，除大连港该类货品吞吐量稳步增长 外，以多种融资形式兴建的液体储罐码头规模宏大。如大连港湾液体储罐码头公司 D〔MTTC) 和大连化学工业公司等。鉴于该类货品的危害特性，为保护大连湾环境及资源，9-Tk待制订溢 漏事故应急反应计划，解决损害赔偿处理中污染损失计算等实际问题。这就必须以获取溢漏 物在海水中的运动状况为前提。 2大连湾易溶保守液体化学品污染扩散数学模型的建立 2.，船运散装液体化学品溢漏扩散的主要运动形式 一般而言，其溢漏后的某一阶段以某一运动形式为主，且辅以多种形式在环境中输移、 扩散。可将其主要运动形式分成五类:(I)强挥发性类蒸汽在空气中的输移和扩散:(2)不溶 于水，在水面以二维形式进行输移和扩散:(3)溶于水，在水白以三维形式输移、扩散;(4)比 重大且不溶于水，沉降于水底:(5)与空气和水发生反应 对某一类物质而言，应根据其理化性质，分析其在环境中主要的运动形式，给出该类物质在 不同时期、介质中的输移、扩散的数学模型。对其进行模拟，以获取其运动的时空分布和归 宿。 2.2模型的建立和流场的准备 因水域环境的不同，易溶保守物质扩散形式不同。适当的输移、扩散和水动力模型取决 于液体化学品的特性以及关心海域水动力特胜。大连湾内水流土要是潮汐流，正规的半日潮 流，水体混合较为充分，对易溶保守液体化学品可采用深度平均的二维潮流数学模型和输移、 扩散模型。 为做好流场的准备，根据浅水潮波理论 建立了大连湾潮流数学模型，采用ADI法、开 边界强迫水位和二维信息控制，得出大连湾M:分潮流场、潮位值 计算结果的实侧资料吻 合较好且缩短了计算时间口 3基本方程及定解条件 根据大连湾的水深条件及扩散范围，本文选用深度平均浓度扩散方程: a(HP) a(HUP)a(HVP) +~一下一一一+ ax ay 二acxH(K,豁氦、Laapy)·、 (1) 其中，H=h+4;h为平均海平面以下的水深:J为平均海平面以上的海面高度;P为 液体化学品在水中的浓度:U.V为X,Y方向上潮流的速度分量:凡、凡为X,Y方向上的 扩散系数;:i1为液体化学品单位时间内溢漏的数量。 在大连湾潮流和潮汐的数值计算基础上，由方程 (1)求解P的定解条件为 在r-0的初始时刻有:尸=口 (2) 在计算区域的开边界P:上，有: -aaPt+·u((n“`”,·1))-oanPe一。4(m) ) 3 尸=0 (涨潮) 叼 在计算区域的闭边界r，上，有: 几 护 n 一 ︸ 曰 抓 ︸ (5) 其中:nn2分别为开边界及闭边界的单位法向矢量;u=L十V:D,为液体化学品. 浓度扩散系数在n:上的投影· 在建立方程 (功 的AD[有限差分格式时，使用 同潮流计算一样的交错式网格 (图 1)aP的差分点 与弓的相同，K,.K,的差分点分别与潮流分量 U和 V 的差分点相同。由于方程 (1)中的平流输运项 汉HUP) 以HVP) 的空间非对称性而不可使用空 ay falT嵘 分，否