关于海上溢油世件的数学模型.docVIP

论文题目： 墨西哥湾原油泄露的图形模拟及环境经济评估方法 论文作者1： 刘铄 在现阶段建模中你善长：□写作 □程序设计 □√数学思维，好突发奇想 □√构建模型的应用能力强 论文作者2： 刘喆 在现阶段建模中你善长：□√写作 □程序设计 □数学思维，好突发奇想 □构建模型的应用能力强 论文作者3： 刘易 在现阶段建模中你善长：□写作 □√程序设计 □数学思维，好突发奇想 □构建模型的应用能力强 摘要： 本论文首先在平静的湖面上利用汽油模拟海上原油泄露事件，分析静止点连续漏油的条件下的漏油过程，可分为重力——惯性力扩展阶段，重力——粘性力扩展阶段，表面张力——粘性力阶段。利用Fay公式在静点持续漏油模型下的优化，推导出各阶段的扩散公式，然后由实际数据发现三个阶段一般只持续5至6个小时，之后油膜的厚度不变，从而得出扩散阶段结束后面积随时间是正比例增大的，这样就为第二问在海面上原油的扩散做好了铺垫。对海面原油的扩散中，我们研究了风场和流场对原油面积和形状的影响，建立了半圆形——半椭圆形模型。然后我们根据墨西哥湾漏油的实际情况下的风场和流场，对每个时间段内的油膜形状和面积进行分析，结果与实际数据吻合很好。在评价原油泄漏对环境和经济的影响的问题上我们进行了全面精细的分析，做出了毒性参数，波及系数，环境生态损失系数的研究模型，运用矩阵的方式列出数据，最后通过计算本次原油泄漏的数据矩阵与其他漏油事故的数据矩阵的欧氏距离，从而对原油泄漏问题做出了较科学的评价。 关键词： 静点连续漏油Fay公式；风场；流场；半椭圆；平面几何分析；综合评价矩阵；欧氏距离 前言 2010年4月21日，美国墨西哥湾“深水地平线”钻井平台发生爆炸造成原油泄漏，引起了国际社会的高度关注．海上原油泄漏给环境和经济带来巨大的损失．正确估计原油扩散在海面上的形状和面积对油污的清理将有着直接意义。 当油被集中泄入水域某一特定地点后，油在重力作用下将向四周扩展，先在水面形成较厚的油膜，而后油膜面积不断扩大，这时重力是油膜扩展的动力，惯性力则阻止油膜扩展，这一阶段为重力——惯性力扩展阶段；随着油膜范围的进一步扩大，油膜变薄，油层的惯性也随之减小，油与水面界层逐渐起控制作用，这时油层的位能和粘性阻力成为主要作用力，这一阶段称为重力——粘性力扩展阶段；当油膜进一步变薄后，表面张力影响增大，超过粘性阻力，油层继续扩展这一阶段称为表面张力——粘性力阶段。 受重力、惯性力和粘性力控制而扩展的阶段历时较短，一般数十分钟至数十小时，视溢油规模而异。 在流动水域中，油膜将随流漂移，其漂移速度和方向大体上与水流一致。如果水面有风时，油膜大体上以风速的3.5％随风漂移。 第一问 本题的条件是： 在平静的湖面上，所以可以认为水为纯水，即密度为1000千克/立方米，水面是没有波纹的，且水面足够大，不考虑风等其他因素，即油在扩散过程中除了重力和粘滞阻力外，不受其他方面的力； 要求利用汽油模拟，93#汽油的密度为725千克/立方米； 漏油点漏油的速度是0.01mL/s，所以该模型为静点连续漏油模型； 在上述假设的条件下，油膜分子团所受力分别为重力、表面张力、粘滞力这几个力所决定。假设是在湖面的漏油是静止点漏油，以漏油点为原点以南北为y轴，以东西为x轴建立坐标系，沿原点竖直向上方向为z轴。 因为汽油扩散时周围各点所受力相同，所以各点的扩散速度相同，所以这种情况下，汽油扩散的形状是以原点为圆心的圆形。 由于随着时间的变化各个阶段油扩散被驱动的主要力不同，所以我们将有的扩散过程划分成三个阶段。 重力惯性力扩展阶段 溢油进入水体后，在开始的一段时间重力—惯性力占主导，由于油水密度差引起的油膜加速塌落，形成有的初始运动，运动方程为： 式中：r为油膜中心径向坐标； u为t时刻r处的油膜运动速度； h为t时刻r处的油膜厚度； 分别为水和油的密度； 有膜前缘运动为： 所以这一阶段的油膜的形状在z，y平面上的投影可以看做正态分布，而在x，y上的投影为圆。 根据上式求解可得： ············Fay公式第一阶段 D为油膜扩展直径； g为重力加速度； ； v为溢油体积； 为待测系数 粘性扩展阶段 油膜扩展一定时间后，其尺度已经足够大，油膜底部粘性切应力与重力逐渐平衡，并占据主导作用，此时运动方程为： 油膜前缘