海面原油泄漏数学建模.docx

《海面原油泄漏数学建模.docx》由会员分享，可在线阅读，更多相关《海面原油泄漏数学建模.docx（11页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、海面原油泄漏数学建模海面原油泄漏的数学建模 论文作者1：（程序设计） 论文作者2：（构建模型） 论文作者3：（论文写作） 摘要 本文主要探讨的是海上原油泄漏所形成油膜的面积和形态的估计问题。以墨西哥湾溢油事务为例，建立了在海面流场、季风风场的作用下不安静的海面上溢油扩散、离散和迁移的组合模型。针对问题一，建立了在安静湖面上志向的状况下油膜在重力、表面张力、粘性力、惯性力这几个主要确定因素的作用下的自身扩展模型。利用fay公式得出油墨在各个阶段随时间改变的扩展尺度，计算出油膜面积为A=17.78t3/4，并用matlab模拟出了海上油膜各向同性扩展的圆形模型。针对问题二，在问题一志向化模型的基础

2、上进行改进。在考虑海面流场和季风风场的作用下建立了由油膜扩展、离散和迁移这三种过程所产生的油膜组合运动数学模型。以二维移流扩散方程为基本限制方程，导出了油膜厚度的表达式，从而建立了油膜的虚拟椭圆的长轴和短轴随时间改变的理论公式和图线。首先在问题一扩展尺度的基础上利用菲克定律增加来了油膜的离散尺度而推导出油膜的扩延迟度。然后又引出“虚拟椭圆”的概念得到实际油膜可观测的尺度和面积公式。并依据资料所述美国墨西哥湾泄漏到海面的原油为400万桶，漏油口第84天被堵住,换算可得墨西哥湾原油的平均泄漏速率也许为Q=9100m3/d，最终得出其漏油面积为A=1.171010m2。最终考虑油膜的迁移运动使油膜的

3、离散产生随机性（本模型对离散量加入随机数），从而使油膜形态发生变更，呈现不规则性，此模型认为漂移并不干脆变更油膜的面积。用欧拉一拉格朗日追踪质点法求得油膜的漂移扩散。针对问题三，了解溢油对经济、环境的各方面的负面影响，针对不同的溢油环境实行合理的治理方式，使溢油对经济和环境的影响降到最小。关键字：溢油 fay公式 菲克定律 扩展 离散 迁移 拉格朗日速度 一、 问题重述 近年来环境问题备受关注，2010年美国墨西哥湾“深水地平线”钻井平台发生的爆炸造成原油泄漏事务引起了国际社会的高度关注。题目要求建立数学模型正确估计原油扩散在海面上的形态和面积，并在此基础上评价美国墨西哥湾海上原油泄漏对环境和

4、经济的影响。问题一，在安静的湖面上利用汽油模拟海上原油泄漏事务，假设漏油点漏油的速度是0.01mL/s，湖面足够大，只考虑汽油的扩散，不考虑风等其它因素建立志向化模型计算水面油膜的形态和面积。问题二，建立在海面上考虑季风、洋流等因素的原油的运动模型并再查找资料，利用所建立模型估算美国墨西哥湾原油泄漏中油膜的形态和面积。问题三，依据问题二的的计算结果，评价一下美国墨西哥湾原油泄漏对环境和经济的影响。二、 模型假设 l 假设溢油模型为为连续匀称溢油； l 被原油污染的区域水质等状况相同，原油的扩散速度相同，即不考虑水体在不同水域的差别； l 海水总体积保持不变，不考虑雨雪渗漏以及自然状况下的自净过

5、程（原油在一般状况下，无法自然消退沉积），海面无任何设施（包括船舶、围油设施）； l 油在大规模扩散仅由重力、表面张力、粘滞力确定，扩散过程中性质不变，如密度、张力等，且在竖直方向处于平衡状态。l 问题一中，在安静的湖面上利用汽油模拟海面原油泄漏，基于湖面足够大，不考虑风等其他因素的假设，忽视油膜的扩散及扩延，即扩展起主要作用，假设漏油点的速度为0.01ml/s，假设。依据fay公式，在不同的扩展阶段忽视次要作用项，保留主要作用项，得到油膜扩展直径的理论公式。l 问题二，原油在海面上的扩散是相当困难的，包括溢油的蒸发、乳化、溶解、生物降解、吸附沉降和氧化等诸多难以量化的困难过程，不是任何一种物

6、理化学现象所能说明的。在此我们忽视这些过程，只考虑油膜的扩展，离散，迁移这三个主要过程。并假设季风、洋流等因素的作用（使油膜产生漂移）体现在两方面：一是使油膜的离散产生随机，只变更油膜形态和位置并不干脆变更油膜的面积。二是假设油膜的质心在油膜上的相对位置变更很小以至可以忽视不计，则油膜的轨迹或运动方式可以等效于质心的运动。三、 符号说明 g 重力加速度 水的运动粘性系数 V 溢油总体积 油的密度 空气和油的表面张力系数 海水的密度 空气和水的表面张力系数 A 油膜面积 油和水的表面张力系数 油膜可视厚度 油膜各向同性的扩展直径 Q 平均漏油速率 油膜迁移速度矢量 表面流速 风速 风速引起的漂流

7、速度矢量 柯氏系数 摩擦影响深度 四、 模型的建立与求解 问题一 在安静的湖面上利用汽油模拟海面原油泄漏事务，溢油的主要过程是在惯性力、重力、粘性力和表面张力的作用下油膜的扩展模型。最终得到的湖面油膜形态应当是圆形。依据fay公式，油膜各阶段的扩展直径计算公式为： 惯性扩展阶段： df=2K1(gV)1/4t1/2 （1） 粘性扩展阶段： df=2K2(gVvw)1/6t1/4 （2） 表面张力扩展阶段： df=2K3(wvw1/2)1/2t3/4 （3） 志向条件下当净表面张力系数小于等于零时，原油停止扩散，此时扩展结束所形成的最大直径为： df=356.8V3/8 （4） 其中，t为从溢油

8、起先计算所经验的时间，=1-o/w，=wa-ao-ow, V=10-5t, K1,K2,K3为阅历系数。据此可推出油膜面积A=105V3/4 在安静的湖面上，以原油泄漏点为中心，向东为x轴，向北为y轴建立直角坐标系， （5） 由matlab作图可得如下图型 图1在安静湖面上的油膜形态模拟图 t/s A/m2 图2在安静湖面上油膜面积时间图像 问题二 在不安静的海面上的原油扩散一共包括三个过程： 1. 第一问模型下的自身扩展。2. 由于受到海面流场和季风作用的不匀称作用、海面破裂波及大尺度漩涡等因素作用而引起的油膜分散或破裂分散的离散现象。3. 由于海面流、风力、波浪等因素作用而引起的油膜漂移的

9、油膜迁移运动。基于上诉三种过程建立不安静海面上油膜组合运动的数学模型，并估算美国墨西哥湾原油泄漏中油膜的形态和面积。油墨在不安静海面上的扩延迟度应当等于扩展尺度与离散尺度的和。(一) 油膜扩展尺度 (二) 离散尺度 油膜在不安静海面上的离散具有随机的性质。有质量守恒原理导出油膜厚度的公式为 （6） 式中油膜中所含第i组分的浓度为Ci；h为油膜的厚度；u、v分别为沿X、Y方向的剪切离散速度；、分别为第i种组分中通过油膜表面和底面的油份通量。求解（6）的关键是确定离散通量与其它时均特性的关系，一般是将其与分子扩散相比拟，采纳Fick定律可得 （7） 式中，Kx、Ky，分别为油膜沿x、Y方向的离散系

10、数。假定Kx、Ky不随空间改变，将（7）式代入（6）式中可得 （8） 上式的解为 （9） 上式又可表示为 （10） 式中，、分别为x、y方向油膜厚度的标准差。油膜沿x、y方向的离散尺度为 （11） 其中，一般取， (三) 油膜的扩延 在溢油的扩展离散过程中，溢油的扩延尺度可表示为扩展尺度和离散尺度的叠加。则油膜沿长轴方向(油膜漂移方向)的扩延尺度 （12） 油膜沿短轴方向(垂直与油膜漂移方向)的扩延尺度 （13） 式中k=A/A（A为客观油膜面积，A为油膜面积，包括观测不到的部分） 为了避开在扩展离散后期，计算的油膜扩延尺度包括油膜边缘观测不到的部分，为了得到实际油膜的尺度，引入“虚拟椭圆”、

11、“虚拟长轴”和“虚拟短轴”。椭圆的实际面积为 A=4dxdy （14） 其中， （15） （16） 据此，以原油泄漏点为中心，以原油的主要扩散方向额为为x轴，以垂直于x轴的方向为y轴建立直角坐标系， （17） 其中，Q=0.01ml/s；=10-810-7；=122.23610-6t1.17； =122.23610-6/10t1.17，用matlab模拟图形得： 图3 海面上原油泄漏形态模拟图 (四) 油膜的迁移 原油泄漏后油膜在水面的迁移运动主要受表面流和风力限制。迁移使油膜的离散产生随机性，（本模型对离散量加入随机数，）从而使油膜形态发生变更，呈现不规则性，此模型认为漂移并不干脆变更油膜的

12、面积。故油膜面积仍按式（15）求得。并依据资料所述美国墨西哥湾泄漏到海面的原油为400万桶，漏油口第84天被堵住,换算可得墨西哥湾原油的平均泄漏速率也许为Q=9100m3/d，最终得出其漏油面积为A=1.171010m2。由于漂移的存在，使得油膜会向一个方向移动，记为方向，假设油膜的质心在油膜上的相对位置变更很小以至可以忽视不计，则油膜的轨迹或运动方式可以等效于质心的运动，其中迁移速度为 （18） （19） （20） 其中，是表面流漂移系数，它包含了油膜的内应力、厚度等因素的影响，可由试验确定，一般取1.0；是风漂流系数，试验数据给出值为35；35；，是风速引起的漂流速度矢量，可由下式给出：

13、（21） 其方向为右偏角度，为柯氏偏转角，是由于地球旋转引起的风矢量偏转角，一般采纳Ekman的解析解确定： （22） 式中，D称为摩擦影响深度，可取为 （23） 是海洋上层的垂直涡粘性系数(近似为常数)；是柯氏系数， （24） 是地球自转角速度，是纬度。 图4 迁移速度矢量合成示意图 (五) 油污污染海疆的面积 设油膜某质点在r，时刻的坐标为R(ti)，该质点在扩展、离散、风力、海流、蒸发和挥发等因素的共同作用下，ti+1时刻的坐标为R(ti+1),则有 （25） 为扩展、离散和漂移在t时段内位移的合成，即 （26） 其中为随机数，值为-1到1，为在t时段内的漂移矢量，采纳积分拉格朗日速度得

14、， （27） 简化模型得， （28） 下图为油膜漂移示意图，假设漂移过程中无面积改变，则时间t内：受污染的海疆面积S=dl。其中， 图5 油膜组合漂移模拟图 问题三 墨西哥湾海上漏油由于受到环境影响因素较大，在扩展，离散，漂移的作用下，影响的海疆面积约为S=dl（2.3万平方公里）。此次大面积的漏油对波及地区的环境、经济造成了很大的破坏。 (一)、对环境的影响 这次漏油事务从起先漏油到基本堵漏胜利持续了84天，总溢油量约700000-800000立方米。溢油量之大，加上油膜扩展离散以及油膜迁移，导致了被污染水域面积远大于油膜面积，这大大增加了海洋环境被污染甚至被破坏的程度。海洋生态面临着严峻的

15、考验，在受污染海疆的656类物种中，造成了大约28万只海鸟，数千只海獭、斑海豹、白头海雕等动物死亡，将有10种动物面临生存威逼，3种珍稀动物面临灭顶之灾（2010年6月数据）。溢油覆盖在海面上，阻碍了浮游植物的光合作用，这从食物链最低端切断了能量的来源，使得处在食物链上层的生物的生存受到了威逼。海面上的溢油对鸟类的危害最大，尤其是潜水摄食的鸟类。这些鸟类以海洋浮游生物及鱼类为食，当接触到油膜后，羽毛能浸吸油类，导致羽毛失去防水、保温实力，另一方面它们因不能觅食而用嘴整理自己的羽毛，摄取溢油，损伤内脏。最终它们会因饥饿、寒冷、中毒而死亡。成年鱼有着特别敏感的器官，因此，它们一旦嗅到油味，会很快地

16、游离溢油水域，而幼鱼生活在近岸浅水域简单受到溢油污染，这种情形会导致鱼群的衰减，从而间接影响捕鱼业的发展。海洋哺乳动物受溢油污染危害的状况是不同的，如鲸鱼、海豚和成年海豹对油特别敏感，它们能刚好地逃离溢油水域，可以避开遭遇污染。但成年海豹和小海狗栖息海滩时，会被油的污染所困，以至于死亡。浅水域通常是海洋生物活动最集中的场所，如贝类、幼鱼、珊瑚等活动在该区域，也包括海草层。溢油对该类水域的污染异样敏感，造成的危害在社会上反应剧烈。假如在这类水域运用溢油分散剂，造成危害会更大。因此，当溢油污染会波及到该类水域时，决策者的首选对策应是如何避开污染，而不是待污染后再实行清除措施，更不适于运用分散剂。溢

17、油对环境的危害程度与溢油类型和环境因素有关。溢油种类不同，对环境造成的危害也不同。例如，重燃料油含芳香烃化合物的量比其他油高，而芳香烃化合物与油中的其他烃类成分相比，具有易溶于水、易扩散、生物难以降解、对生物毒性大等特点。假如溢油为重燃料油，不能刚好得到处理，对海洋环境的污染危害要比其他类型的溢油大。溢油对环境的危害程度还与环境自身的特征有关。溢油发生地点是否是敏感区，溢油发生的季节是否是鱼类产卵期、收获期，不同的海况，都影响溢油的危害程度。相同规模的溢油事故，发生在开阔水域要比发生在封闭水域的危害程度要低；发生在海洋生物生长期要比发生在其产卵繁殖期的危害低；发生在荒芜人烟的偏僻地带，靠其自身

18、的净化复原，不用实行任何措施，那么，这种危害可以不计。(二)、对经济的影响 英国BP公司作为肇事者令自己承受了巨额的损失，市场估计，BP担当的赔偿金、清理油污费、建立减压井等相关堵漏设备以及受影响各州的相关支出高达370亿美元。事故还将影响BP公司的油气开采业务，严峻损害了公司信誉。美国墨西哥湾沿岸各州的经济也受到了很大影响。这次漏油事故波及墨西哥湾及周边地区的旅游、农业、航运、能源等多个领域，影响到高达2320亿美元的经济体系的运转。事故发生后，海面污染严峻，波及区域的旅游业受到严峻冲击，同时超过8.8万平方千米水域的渔业被迫停止。美国政府为此增加的各类费用必定影响到上述各州和联邦政府的财政

19、预算，有估计称这次漏油事故将会使美国联邦政府预算增加数百亿甚至千亿美元，这对处在脆弱的经济发展中的美国将会是一个危机。同时和高达数千亿美元的各方面总损失相比，BP公司设立的200亿美元赔偿基金明显是杯水车薪，那么这次事故对经济所产生的负面影响将持续很长一段时间。墨西哥湾海产品占美国市场总量的20，因为漏油事故的影响，消费者将不会购买来自污染区的海产品；同时大面积水域污染必定导致旅游业的亏损。在航海方面，由于运输距离的增加，产品成本必定抬高，这无论是对消费者还是美国的经济都是不利的。石油是工业的血液，对全球经济的发展起着至关重要的作用，但这也确定了它对人类的牵制力。所以，人类在利用石油进步的同时，要加强对开采石油的监督与调整，务必在平安、环保的环境下合理的开采石油，让人类发展与自然界和谐共存。五、 参考文献 海面溢油、离散和迁移的组合模型 武周虎 赵文谦 海上溢油运动数值模拟方法的探讨与改进 刘伟峰 孙英兰 2的行为和归宿数学模型基本理论与建立方法的探讨 李冰绯 附：部分程序见program.txt文件 TXT文件附件无