海上原油泄漏问题

海上原油泄漏问题问题

摘要

本文分析了渤海水体中的行为表现，根据基本控制方程和海洋石油的物性，借助费克定律和物质守恒定律可以导出油滴在水中的扩散微分方程，因此建立起的原油在水中扩散的数学模型。采用CFD商业软件PHOENICS3.6.1对在水流恒定情况下渤海石油在海面和海下扩散的稳态情况进行数值模拟，研究表明了渤海原油扩散扩展性状的特点。可得当在静水海面上，即水流为0，水面油层形状是圆形。当水流速不为0时，水面油层形状视为椭圆形。

在漏油点被控制之前，由于油层很厚，扩展是溢油的最主要行为，可视全部为扩展过程。基于原油在海水中必受到水流、季风等影响，其扩散的形状是椭圆，可建立油膜椭圆扩散模型。在漏油点控制后，原油受到重力、季风、水流、潮汐、海浪等因素影响，其过程是原油的漂移、粉碎、蒸发和乳化过程。可建立油粒子模型，把石油模拟成油粒子云团，采用随机走动法来模拟粒子云团的扩散过程。由于每个粒子的随机运动而导致整个云团在水体中的扩散过程。通过两个阶段模型的建立，可求得渤海水面油层面积，进而以此对渤海漏油对我国环境和经济的影响作出具体评价。并对我国海洋石油管理、海洋污染防治、海洋经济等方面提出合理化建议。

关键词: 扩散; 海洋; 数值模拟，油膜椭圆扩展，油粒子模型

一．问题重述

2010年4月，英国石油公司(BP)位于渤海区域的钻井平台发生爆炸并引发大火，造成美国历史上“最严重的一次”漏油事故；壳牌公司日前承认，其位于英国北海海面上一座钻井平台也发生泄漏事故，已有200吨石油漏入北海；康菲石油公司在始于今年6月发生的两起渤海湾溢油事故中，约115立方米(约700桶)的原

油渗漏到海面。海上原油泄漏事故频发，给周边甚至世界环境和经济带来了巨大的损失。因此正确估计原油扩散在海面上的形状和面积对油污的清理将有着直接意义。

请你们用数学建模方法研究解决：

⑴ 假设漏油点漏油的速度是0.01mL/s ，请你们建立数学模型计算水面油膜的形状和面积(假设湖面足够大，只考虑汽油的扩散，不考虑风等其它因素)； ⑵ 在海面上，原油的扩散相当复杂，还要考虑季风、洋流等因素(当然还有其它因素影响)，请你们在⑴的基础上建立相应模型，再查找资料，利用所建立模型估算我国渤海湾原油泄漏中油膜的形状和面积；

⑶ 根据⑵的计算结果，评价一下我国渤海湾原油泄漏对环境和经济的影响。

二．问题分析

近年来，随着经济的发展，对石油等需求的依赖。渤海湾港口的大型油轮将会越来越多，而发生船舶油污事故的几率将增加。海上大量溢油危害海洋环境资源和地区经济发展。从工程观点来看，一旦溢油事故发生，首先迫切要了解的是油层扩展有多快，将被油膜覆盖的海域有多大，以及在某一时间范围内油膜会扩展到何处等等一系列实际问题。这些信息对寻求控制和处理这类突发性

环境危机的最佳对策是至关重要的，因此建立溢油扩散的数学模型，研究渤海油海上溢油扩散的分布规律具有重大意义。

原油至漏油点刚进入水体后，由于油层很厚，会迅速向四周扩展。在溢油最初至漏油点被控制住，扩展是溢油动态行为最主要的过程。该过程的长短与油的种类、品质、粘性、温度等自身性质密切相关。同时溢油量越大，持续时间也越长。在流体问题求解中所需求解的主要变量的控制方程都可以表示成以下通用形式[3]：

()()()div U div grad S t

φφρφρφφ∂+=Γ+∂ 式中：φ为通用求解变量；φΓ为广义扩散系数；S φ为广义源项。

根据基本控制方程和海洋石油的物性，借助费克定律和物质守恒定律可以导

出油滴在水中的扩散微分方程，因此建立起的原油油在水中扩散的数学模型。采用CFD商业软件PHOENICS3.6.1对在水流恒定情况下渤海石油在海面和海下扩散的稳态情况进行数值模拟，研究表明了渤海原油扩散的特点。模拟得到海洋面原油泄漏在水面形成的形状。

图1 油的风化过程

由于渤海区域原油，其粒径在10~20μm，密度为大于淡水而稍低于海水，常温下粘度大、流动性极差，因此原油在水体中的扩散行为、其油膜及浓度分布情况有其特点。影响这种混合特性的主要因素[2]有燃料浓度，海水的盐度和外界环境温度等。当海水的盐度较高时，沥青微粒数量会增多，但是结合动态环境因素、水流及波动，沥青微粒会连续不断混合到海水中，以阻止沥青颗粒表面形成油膜，表层下水流和流势是影响沥青微粒在海水分散的主要因素。对于进入水体的原油，只有充分了解其在水体中诸多的动力学和非动力学过程，才能合理地模拟其扩散过程。

图2原油在水面受到力的作用分析

Fay首先提出了油膜在平静水面的自身扩展理论。该理论认为溢油进入水体后在重力、惯性力和表面张力作用下，迅速扩散。并根据扩展期间主导力的不同，而将扩展划分为3个过程，认为油膜呈圆形扩展。

W.T.lehr.R.T.Fraga和M.S.Belen对Fay理论进行了修正，得出油膜在海洋里德扩散不可能不受到风的作用。因而其扩散不是呈圆形而是椭圆形的。其主轴同风向一致，这就是油膜椭圆扩展模型。此过程适用于漏油点控制前，由于两次漏出的石油总量为115立方米，漏油点漏油速率为0.01ml/s。可计算出每天漏油量为0.86m3/d，漏油源被控制所需要的时间为115/0.86=133.2d。即在此时间段建立油膜椭圆扩散模型。

在漏油点被控制后，扩散不是其主要过程。渤海石油泄漏进入水体环境，涉及溢油发生时的扩展、风化、在风和海流作用下的漂移、在波浪作用下油膜的破碎和结块、沉降和上浮、在海岸的搁浅附着、在海底的附着、生物降解等等一系列物理、化学和生物过程。随着时间的推移，会逐步消失，这就是水体的自净作用。这个过程，即物理和化学变化过程，统称为风化，它包括了蒸发、溶解、乳化、分散、沉降、光氧化和生物降解等过程，如图1所示[1]。从短期来看，蒸发和乳化过程是重要的风化过程，对石油的残留量及其组成、性质、状态起决定性作用。而从长期来看，光氧化和生物降解作用决定了石油的最终归宿。

在此阶段，可以把浓度场模拟为由大量的粒子组成的“云团”，其中每一个粒子携带一定数量的示踪物质。每一个模型粒子在特定的流场条件下发生平移，在重力和浮力的作用下发生垂直方向上的位移(涌升或者沉降)，即所谓的平流(或称对流)，适宜用拉格朗日法模拟。模型粒子的扩散过程是由于剪流和湍流而引起，适宜采用随机走动法来模拟粒子云团的扩散过程。湍流可视为随机流场，而每个模型粒子在湍流场中的运动则类似于流体分子的布朗运动。由于每个粒子的随机运动而导致整个云团在水体中的扩散过程。这种方法实际上是确定性方法和随机方法的结合，即采用确定性方法模拟平流过程，采用随机性方法模拟扩散过程。

油粒子定义为一些很小的圆球、直径在10-1000∥m之间。考虑到油滴直径的变化范围，要精确地表示一个溢油膜所需要的实际粒子数应是相当大的。因