一种油包水乳化液溢油量等效估算方法研究

袁 丽， 孔德明， 陈基亮， 仲美玉， 张晓丹， 谢贝贝， 孔令富

1. 成都信息工程大学计算机学院， 四川 成都 610225 2. 燕山大学电气工程学院， 河北 秦皇岛 066000 3. 燕山大学信息科学与工程学院， 河北 秦皇岛 066000

溢油进入海域由于自身性质和海面环境的不断影响， 会经历一系列复杂的物理化学变化[1]。

其乳化阶段， 溢油与海水发生作用， 会形成不同油水比的溢油乳化物。

根据溢油的乳化程度可将溢油污染分为“油包水”、 “多重”、 “水包油”三种。

乳化状态下的部分溢油， 会呈“巧克力慕斯状”或“黑褐色泡沫状”[2]， 比重和粘度变大， 且不容易消散， 从而妨碍大多数机械回收设备的有效工作， 给海洋环境造成更大损害。

所以对海面乳化溢油进行及时与准确地识别、 快速估算其溢油范围和溢油量， 对于海上溢油污染治理和减少海洋环境的损害具有十分重要的意义。

近年来多种遥感技术可应用于溢油监测， 主要包括光学遥感技术和微波遥感技术。

其中激光诱导荧光(laser induced fluorescence, LIF)探测技术属于主动光学遥感， 指依靠激光源发射激光到海面, 海水及各种溢油成分受激发射含各种成分信息的荧光[3]。

通过它能够探测到独特的溢油荧光光谱特征， 这使得在未乳化溢油油种鉴别以及油膜厚度评估方面有了可依据的信息， 但对于乳化溢油的监测尚未形成一套可以有效使用的理论研究成果和成熟的探测方法。

因此， 利用LIF探测技术评估海面乳化溢油范围和污染程度是目前该领域的一个研究热点。

对于油包水型乳化溢油， 随溢油中水滴的逐渐增多， 其尺寸和分布也会不断发生变化， 这些变化影响着油包水乳化液的固有光学性质[4]， 从而使LIF系统发射激光到海面获取的荧光信号不断变化。

最直接的方法是建立两者之间的量值关系模型进而通过获取的荧光值来反演对应的油水比， 实现乳化溢油量的评估。

但关于乳化液中水滴的吸收、 散射作用的计算还没有相关的研究成果可供使用, 尤其是散射系数定量化难度很大， 所以建立精确的量值模型非常困难。

鉴于此， 设计了油包水型乳化溢油的等效估算模型， 提出了一种基于LIF探测机理的等效估算方法， 并通过实例对该方法的有效性进行了验证， 给出了方法对应的适用范围。

1.1 等效模型

油包水型乳化溢油是由于海浪、 风对海面溢油的不断作用， 一部分水滴进入溢油， 并悬浮在溢油中而形成。

当发射激光到海面， LIF系统获取的油包水乳化液荧光信号主要为乳化液中油组分受激产生， 乳化液中的水会对荧光产生影响作用。

如果能将水的影响作用进行恰当的描述和处理， 那么乳化溢油量的估算问题即可得到解决。

等效法[5]是常用的一种科学研究方法， 当面对一个较为复杂的问题时， 通常可以尝试提出一个简单的方案， 并使它们的效果几近相同， 从而将问题化难为易， 求得解决。

因此考虑采用等效思想来简化乳化溢油的溢油量估算问题。

根据油包水型乳化溢油， 油为连续相， 可将其等效看成一定厚度的油膜；

而水为分散相， 将油中分散的水滴进行聚集， 会形成一片片的薄水， 将成片的薄水也看成连续， 这样所有连在一起的水滴可看作覆盖在油膜上面的均匀水层；

另外考虑油包水乳化液与大气接触的主要为溢油， 所以在等效的水层上面覆盖一层薄油， 并且此处的薄油只看作一个油面， 而不考虑厚度的影响， 以达到实际乳化液与等效模型在激光入射和荧光出射时具有相同的外部环境。

最终建立如图1所示的模型， 即将油包水乳化液等效看成上面覆盖油面和薄水层的油膜。

图1 油包水乳化溢油与等效模型对比图Fig.1 Comparison between water in oil emulsion spillage and equivalent model

1.2 等效方法

对于等效模型中水层产生的影响， 不考虑其厚度， 即将水层的吸收散射作用忽略；

和覆盖的油面一起， 均只考虑由其界面产生的反射折射作用。

这样可将油包水乳化液中水的影响处理为等效模型中的界面作用。

基于等效模型， 根据图2所示光的辐射传输过程， 假设探测系统激光的发射功率为Ρ0， 经过大气的传输吸收到达油面， 考虑大气为均匀分布， 激光强度衰减为Ρ0exp(-βLH)， 其中，βL为激光在大气中的衰减系数(也称消光系数， 是吸收系数和散射系数之和)，H为探测系统距离海面的高度。

图2 等效模型的光线传输与界面作用示意图Fig.2 Schematic diagram of light transmission and interface action of equivalent model

在气油界面， 入射光发生反射， 除部分反射进入大气， 大部分透射进入油中， 根据光学原理[6]， 将入射光除去反射光即为折射光能量。

设ρL1为气油界面的反射系数， 则在气油界面下， 激光强度转变为Ρ0exp(-βLH)(1-ρL1)。

油面下方的薄水层， 消光作用不考虑， 可视为水面， 所以认为入射光又陆续经过油水界面、 水油界面， 然后到达油膜层。

设激光在油水界面的反射系数为ρL2， 在水油界面的反射系数为ρL3， 将水油界面下激光强度记为PL， 则

ΡL=Ρ0exp(-βLH)(1-ρL1)(1-ρL2)(1-ρL3)

(1)

(2)

式(2)中，βF为荧光在大气中的衰减系数，ρF1为荧光在油气界面的反射系数，ρF2为荧光在水油界面的反射系数，ρF3为荧光在油水界面的反射系数。

油膜层的荧光转换效率为η， 是受激发射的荧光功率与激发激光功率的比值， 则

ηΡL{1-exp[-(ke+ki)d]}=KF

(3)

式(3)中，ke和ki分别为激发波长、 发射波长处油层的消光系数。

发射的荧光只有一部分被探测系统接收， 由于油层产生的荧光均匀分布于4π空间， 假设接收器的接收面积为S， 溢油的折射率为mF， 则对应的立体角为

(4)

将探测系统接收的荧光强度记为PF， 则根据式(4)可求得

(5)

将式(1)、 式(2)、 式(3)和式(5)进行联立， 可得

(1-ρF2)(1-ρF3)

(6)

式(6)中，S和H与探测系统相关， 看作系统常量；

当溢油油种确定， 各个界面的反射系数也为常量值， 所以设定

exp(-βFH)(1-ρF1)(1-ρF2)(1-ρF3)=R

(7)

式(7)中，R为一个固定值。

将式(7)代入式(6)， 经过运算整理得出

(8)

油种已知的前提下，ke和ki为已知量；

当探测系统确定，R值也可对应求出， 因此只要将系统测得的荧光值代入该计算公式， 即可求出油包水乳化液等效的油膜厚度。

2.1 误差表示

假设油包水乳化液的含油率(体积浓度)为yo、 厚度为d1(激光能够到达)， 我们将D作为乳化液的实际溢油厚度， 如图3所示， 根据油量相等， 那么含油率yo、 厚度d1与实际溢油厚度D的关系可表示为

图3 实际溢油厚度示意图Fig.3 Schematic diagram of actual oil spill thickness

d1Syo=SD

(9)

式(9)中，S为海面溢油乳化液的面积。

将式(9)进行约简， 可推出

d1yo=D

(10)

实际溢油厚度D与等效油膜厚度d相差的值即为误差， 设为m， 则

(11)

可见， 当油包水乳化液和等效模型获取的荧光值相等时， 误差的变化取决于油包水乳化液的含油率yo、 厚度d1。

而具体值与油种， 系统入射激光等因素相关。

2.2 实例设计

根据式(11)， 只能得到影响误差的相关因素。

所以需要设计实例将误差的值具体化， 以分析和验证等效方法的适用性和有效性。

以油品Romashkino原油的油包水乳化液为例， 分别测定相同含油率不同厚度和相同厚度不同含油率的油包水乳化液一系列荧光值， 并建立实际溢油厚度与荧光值的关系曲线；

在等效油膜厚度公式里代入相关参数值， 可得到等效油膜厚度与荧光值的关系曲线。

通过两种关系曲线的比较， 进一步获得和分析误差的具体值。

设油包水乳化液的厚度d1为4 μm， 配置相同厚度的不同含油率的油包水乳化液， 含油率范围为0.6～0.9， 并获取其荧光强度值；

设油包水乳化液的含油率yo为95%， 配置相同含油率的不同厚度的油包水乳化液， 厚度范围为5～15 μm， 并获取其荧光强度值。

图4(a)为厚度4 μm的不同D的油包水乳化液的荧光值结果， 用非线性最小二乘法[7]进行拟合。

通过拟合曲线可以看出， 该原油的油包水乳化液荧光值随实际溢油厚度的增加而逐渐降低；

图4(b)为含油率95%的不同D的油包水乳化液的荧光值结果， 也用非线性最小二乘法进行拟合。

通过拟合曲线可以看出， 该原油的油包水乳化液荧光值随实际溢油厚度的增加先增加而后基本不变。

图4 实际溢油厚度与荧光值关系曲线示意图(a)：
等效厚度为4 μm不同浓度油包水乳化液；

(b)：
含油率为0.95不同厚度油包水乳化液Fig.4 Fluorescence intensity vs.actual oil spill thickness(a)：
Water in oil emulsions with different concentrations；

(b)：
Water in oil emulsion with different thickness

根据表1和表2中所列出的原油和海水的相关光学参数：

表1 原油的光学参数[8]Table 1 Optical parameters of crude oil

表2 海水的光学参数[9]Table 2 Optical parameters of seawater

可计算得出(1-ρL1)(1-ρL2)(1-ρL3)(1-ρF1)(1-ρF2)(1 -ρF3)=0.911 86；

另外对应的ke=393 960，ki=149 540，η=0.038[4, 10]。

将这些参数值代入等效油膜厚度计算式(8)， 经过整理， 可将荧光值与等效厚度的关系近似表示为

PF=0.034 65-0.034 65exp(-0.543 5d)

(12)

该关系用曲线表示如图5， 可看出荧光值随等效油膜厚度的增加而增加， 当油膜厚度达到10 μm左右时， 荧光值趋于饱和状态。

图5 等效油膜厚度与荧光值关系曲线示意图Fig.5 Fluorescence intensity vs.equivalent thickness

2.3 误差比较

将荧光值与实际溢油厚度、 等效油膜厚度的关系曲线放在同一坐标系， 则当纵轴值相同， 横轴值的差即为两者的误差值。

两种关系曲线的交叉点， 代表等效油膜厚度和实际溢油厚度值相同， 此交叉点的上下， 分别表示等效油膜厚度值比实际溢油厚度值偏大与偏小。

交叉点以上或以下， 离交叉点越远， 误差值越大。

交叉点及其上下一定区域， 就是误差允许的范围。

在不同浓度油包水对应实际溢油厚度时， 含油率取极限值1时， 等效油膜厚度比实际溢油厚度偏小的值最大；

在不同厚度油包水对应实际溢油厚度时， 荧光值达到饱和时， 等效油膜厚度比实际溢油厚度偏大的值最大。

具体分析如图6所示。

图6(a)的实际溢油厚度曲线反映油包水乳化液不同含油率yo与荧光值的变化趋势， 所以两种关系曲线的对比， 也反映油包水乳化液不同含油率与误差的关系：

当等效厚度值d与实际溢油厚度值D相等， 即两者厚度均为3.52 μm时， 对应乳化液含油率yo大约为0.88；

当含油率yo大于0.88时， 等效厚度值d小于实际溢油厚度D， 且误差值随含油率的增加会逐渐增大。

当含油率趋近极限值1时， 即含水为零， 两者的误差值最大。

此时的实际溢油厚度可看成4 μm， 计算的等效油膜厚度大约为3.5 μm， 误差为10%左右。

当含油率yo小于0.88时， 等效厚度值d大于实际溢油厚度D， 且误差值随着含油率的降低而逐渐增大， 此时的最大误差值所对应的含油率由该油包水最低含油或者最低含水量决定。

假设当误差也为10%时， 此时对应的含油率为0.8。

以上说明， 当油包水乳化液含油或含水在一定范围时， 采用等效方法可有效处理乳化液中水对荧光的影响作用， 实际溢油厚度与等效油膜厚度之间的误差较小。

图6(b)的实际溢油厚度曲线反映油包水乳化液不同厚度d1与荧光值的变化趋势， 所以两种关系曲线的对比， 也反映油包水乳化液不同厚度与误差的关系：

图6 实际溢油厚度与等效油膜厚度误差比较图(a)：
不同浓度时的实际与等效关系图；

(b)：
不同厚度时的实际与等效关系图Fig.6 Error comparison between actual oil spill thickness and equivalent oil film thickness(a)：
Water in oil emulsions with different concentrations；

(b)：
Water in oil emulsion with different thickness

当等效厚度值d与实际溢油厚度值D相等， 即两者厚度均为6 μm时， 对应的油包水乳化液厚度d1大约为6.3 μm。

当油包水乳化液厚度d1小于6.3 μm时， 等效厚度值d小于实际溢油厚度D， 且误差值随厚度的增加会逐渐减小。

当油包水乳化液厚度d1大于6.3 μm时， 等效厚度值d大于实际溢油厚度D， 且误差值随厚度的增加会逐渐增大， 根据等效厚度曲线的变化趋势， 当等效厚度在大于10 μm时， 荧光值变化趋于平缓， 此时对应最大误差值， 且误差约为40%。

假设当误差也为10%时, 此时对应的油包水乳化液厚度大约为8 μm。

以上说明， 当油包水乳化液厚度在一定范围时， 可使得实际溢油厚度与等效油膜厚度的误差在很小范围， 即等效处理方法是有效的。

综合以上， 当油包水型乳化溢油的含油率和厚度均满足相应范围时， 如此例中油包水乳化液含油率在0.8以上， 并且厚度在6.3μm左右， 等效处理的方法产生的误差会较小， 可有效评估乳化溢油量。

海面乳化溢油量的估算是一个难点问题。

目前关于乳化液的光谱定量分析理论和方法十分欠缺， 因为油水比的变化与荧光值之间具有复杂的关系， 使得直接量化极其复杂。

对于该问题的解决， 本文提出的等效模型， 将油包水型乳化溢油中水的消光作用转换成界面的折射影响， 从而将乳化溢油溢油量的计算转变为油膜厚度的计算， 并推导出计算公式， 最终实现对海面乳化溢油的估算分析。

提出的等效方法简化了对海面油包水型乳化溢油的溢油估算， 是该问题的一种创新性解决途径。

下一步将研究采用相似思想对水包油型乳化溢油建立等效方法。