海洋污染模型.doc

土地资源网-专业的土地流转平台 PAGE  土地资源网--/专业的土地交易平台、从事土地事务代理、土地咨询、土地买卖等业务 溢油模型的现状及研究方法 溢油模型的概念溢油模型是根据环境科学原理建立的以水为载体，用于模拟溢油轨迹、归宿及其环境影响的综合模型。溢油模型涉及到溢油在环境中经历的风化(风化包括蒸发、溶解、乳化、沉降、化学氧化和生物降解)扩散、漂移，最终将如何归宿等整个过程，可以在计算机上快速模拟并动态显示不同时刻、不同厚度的油膜轨迹，并能评估溢油对海洋环境生态系统的影响。 溢油模型的基本原理根据溢油模型的数学计算理论不同，可以将溢油模型分为两种类型。一种为利用欧拉－拉格朗日理论体系建立的溢油模型，用欧拉理论建立的流场和温盐场计算模型，利用拉格朗日理论体系建立的溢油扩散漂移轨迹的计算模型；另一种为利用蒙特卡罗方法建立的溢油扩散漂移轨迹计算方法，它是在欧拉—拉格朗日追踪的方法求得溢油扩散漂移的基础上，增加了溢油扩散漂移随机数的计算。 欧拉－拉格朗日理论模型用欧拉体系的数值模拟方法计算潮流预报模型，用拉格朗日体系的数值方法计算全动力轨迹模型和溢油归宿模型。全动力轨迹模型是将溢油看作由很多小的运动质点组成，每一质点都在一定的动力作用下在自然界中运动。动力是由潮流、湍流、风海流、重力沉降组成。计算溢油浓度和油膜厚度。溢油归宿模型包括油膜重力扩展模型、溢油蒸发模型、岸线吸附溢油模型、溢油水体溶出模型和溢油沉降再悬浮模型。 蒙特卡罗方法模型是通过运算溢油扩展预测程序显示溢油归宿轨迹。溢油扩展预测程序包括风场预报子程序、蒸发子程序、乳化子程序、潮流场预报子程序、油膜扩展输移子程序和图象处理子程序等6个子程序。风场预报子程序运算风场模型，蒸发子程序运算蒸发模型，乳化子程序运算乳化油的入水通量模型，潮流场预报子程序运算二维或三维潮流场模型，油膜扩展输移子程序运算油膜扩展输移模型，图象处理子程序将预测结果、油膜扩展图形与输移过程用彩色图象动画显示，并同时显示有关文字数据，油膜厚度用彩色表示。 、三维海洋溢油预测模型的建立两篇文章对海洋油污染模型的建立进行了说明。溢油物理化学性质灰色模式的研究，也是油污染的研究方法。有限时间排放的污染扩散模式及污染特征 潮汐流动中污染物排放的一种远区计算模型（可以认为是欧拉计算方法的潮流预报模型） 该模型是适合于潮汐流动中污染物排放的远区动态平面二维数值计算模型。 模型的计算特点是将流场和浓度场分开计算，即所谓“非耦合解”，运用有限节点法计算了流场，用改进的有限节点法计算了浓度场。该模型应用于长江大厂镇江段，其计算出的流场和该江段典型污染物 F-的浓度场与实测一致。模型用于计算潮汐对污染物浓度的影响。 潮汐流动中污染物排放的数学模型的基本方程 Navier-Stokes方程和质量守恒方程。在紊流情况下，采用雷诺的时间平均法，并对方程沿垂向积分，引入 Boussinesq假定、涡粘性假定和静压假定，并且忽略不计垂向加速度、粘性切应力在水平方向上的变化率和水表面切应力，可得到平面二维浓度场计算的基本方程： 式中： u, v分别为 x, y向的垂线平均流速 t：时间 ξ：水位即基准面至水面的垂直距离 h：为基面至河底的垂直距离 f：柯氏力系数，f =2Ωsinθ,Ω为地转角速度，θ为纬度 C：为谢才系数 εx, εy ：分别为 x, y向水流紊动扩散系数 c：污染物浓度 Dx, Dy ：分别为 x, y向污染物紊动扩散系数， S：为污染物的源项。 上述基本方程为初边值问题，要有适当的定解条件才能求解。闭边界上，取法向流速为零，污染物浓度的法向梯度为零。与其它水域相通的开边界上，用水位或流速过程和污染物浓度过程来控制。初始时刻的水位、流速、污染物浓度均可取常值。 2.2　计算方法 为了求解基本方程组，选择适当而有效的计算方法是非常必要的。目前数学模型中应用的计算方法很多，如有限单元法、有限差分法、控制体积法、边界拟合法、破开算子法等等，各种方法各有其优点缺点。我们选用一种介于有限差分法和有限单元法之间的方法——有限节点法。将平面计算区域划分为有限个三角形网格。在每个三角形单元上，未知函数用线性插值函数来近似。任意一个函数 f( x, y) C(水位 ,流速 ,污染物浓度 )可用下式来表示。 f ( x, y) =a + bx + cy ( 5) f ( x, y)：所指内容 将三角形单元三个节点的坐标值及函数值代入方程 ( 5)，可求得 ( 5)式中三系数 a, b, c,则三角形单元内任一点的函数可表示为： 其中 ai, aj, ak, bi, bj, bk, ci, cj, ck 是由三节点的坐标值完全确定，fi, fj, fk 为三节点的函数值，△e为三角形单元面积。则一阶偏导数为： 图1　三角形单