海面溢油油包水乳化液多种探测参数的BRRDF仿真研究

张晓丹，孔德明，袁 丽，孔德瀚，孔令富

1. 燕山大学信息科学与工程学院，河北 秦皇岛 066004 2. 燕山大学电气工程学院，河北 秦皇岛 066004 3. 河北环境工程学院信息工程系，河北 秦皇岛 066000

引 言

随着海上石油运输和石油开采业的快速发展，海面溢油事故频发[1]。及时监测海面溢油信息、鉴别溢油种类、估算溢油量及污染程度是快速有效治理溢油污染的基础与前提。近年来遥感技术已广泛应用于海面溢油监测中，其中激光诱导荧光(laser induced fluorescence, LIF)探测技术不仅能识别溢油种类，且对溢油范围和污染程度的评估同样具有巨大潜力，已被公认是目前海面溢油遥感领域最先进的探测技术[2-3]。

溢油进入海域将迅速扩散成油膜，并受海面风浪、光氧化、生物降解等作用逐渐形成了水包油或油包水乳化液。近几年国内外对海面溢油检测已开展了大量研究[4-5]。但都是针对海面溢油污染中未乳化阶段油膜和乳化阶段水包油进行的建模与研究，而对乳化阶段油包水这一重要内容的研究较少，且目前仅集中在微观结构特征、分散相水滴粒径分布方面[6]，光学特性方面的研究鲜有相关文献报道。双向反射再辐射分布函数(bidirectional reflectance and reradiation distribution functions, BRRDF)可表征目标在激光照射下的荧光特性[7]。张晓丹等[8]利用BRRDF分别对未乳化阶段油膜和乳化阶段水包油进行了荧光特性研究，将BRRDF引入到了海洋遥感领域的研究中。

我们基于蒙特卡罗方法和Mie散射理论，应用双向反射再辐射分布函数，针对油包水乳化液建立光子传输模型，开创性地模拟油包水乳化液的fBRRDFcosθrcosθi，分析fBRRDFcosθrcosθi与油包水的含油率、厚度、探测时入射接收角度等相关探测参数的关系。

1 油包水乳化液的仿真模型描述

采用蒙特卡罗方法建立光子传输模型。此模型的设定条件为： 气水界面水平、海水无限深、水体内无其他光源； 为能接收到各个方向的荧光信号，探测器位置尽可能远。在给定角度处入射1亿光子，在海面上半球1 296个立体角内分别固定一个虚拟探测器以接收逃离海面的荧光光子，通过统计所有荧光光子在各立体角内的权重得到BRRDF数值。图1(a)为光子在介质中的传输示意图。其中油包水乳化液的显微结构放大图如图1(b)所示。

图1 光子在介质中传播的切面示意图(a)及油包水乳化液的显微结构图(b)Fig.1 Cross-section of photon propagation in medium (a) and microstructure of water-in-oil emulsion (b)

如图1(a)所示，激光入射到介质，可能发生的折射、散射、吸收和辐射荧光等由介质的吸收系数、散射系数和体散射相函数来确定。由于油包水乳化液水滴粒径尺度集中分布在0.01～15 μm范围内，这与可见光波段相当，因此利用Mie散射理论可计算得到介质的光学参数。本工作的油包水乳化液和文献[8]中水包油乳化液的荧光物质分别以连续相油和分散相油滴形式存在，故在整个传输模型中涉及的吸收系数、散射系数等光学参数存在差异性。油包水乳化液的吸收系数、散射系数表达式为

a(λ)=ρao(λ)+aw(λ)

(1)

b(λ)=bw(λ)

(2)

式(1)和式(2)中：λ为波长，ρ为油包水乳化液的含油率，ao(λ)为纯油的吸收系数。aw(λ)和bw(λ)分别为相应浓度乳化液水滴粒子系的吸收系数、散射系数。表达式分别为

(3)

(4)

式(3)和式(4)中：Qabs和Qsca分别为吸收、散射系数效率因子，α为水滴粒子的尺寸参数，m为水滴粒子的复折射率，r为水滴粒子半径，f(r)为水滴粒径分布概率密度函数。

fBRRDF是描述目标表面出射微分辐亮度与入射微分辐照度比值关系的函数。基于LIF技术探测海面溢油时，相同入射接收条件，LIF系统在某一接收波长处的荧光功率Pr正比于fBRRDF与入射接收天顶角余弦的乘积即fBRRDFcosθrcosθi[8]。故建立海面溢油油包水乳化液的BRRDF仿真模型，获得fBRRDFcosθrcosθi与基于LIF技术油包水乳化液探测参数的关系，为实际LIF探测及通过LIF系统接收的荧光信号判断海面溢油油包水污染情况提供理论指导和重要依据。

2 油包水乳化液的仿真研究

2.1 仿真参数处理

计算不同油品不同含油率的油包水乳化液在405和500 nm波段的吸收系数、散射系数如表1和表2所示。可看出，油包水乳化液随含油率的升高其吸收系数逐渐增大，散射系数逐渐减小。这是由于原油主要表现为吸收作用且吸收系数远大于海水吸收系数的原故。

表1 1号油包水乳化液的吸收系数、散射系数Table 1 Absorption coefficient and scattering coefficient of No.1 water-in-oil emulsion

表2 2号油包水乳化液的吸收系数、散射系数Table 2 Absorption coefficient and scattering coefficient of No.2 water-in-oil emulsion

乳化液发射的荧光强度受荧光量子产率影响，荧光量子产率与乳化液含油率的关系如图2所示。可看出荧光量子产率随含油率的升高呈下降趋势。

图2 荧光量子产率与含油率的关系曲线Fig.2 Relationship between fluorescence quantum yield and oil content

2.2 仿真与分析

目标的辐射特性是其相关参数对LIF探测及系统设计等有重要影响。研究油包水乳化液含油率、厚度及探测时入射接收角等参数下的fBRRDFcosθrcosθi对基于LIF的海面溢油监测、溢油量评估、溢油污染处理等有重要意义。

2.2.1 不同浓度的仿真与分析

仿真中构建三维圆柱坐标系表示fBRRDFcosθrcosθi的分布，纵坐标为fBRRDFcosθrcosθi，上方圆盘表示极坐标和角度坐标，分别代表光子的天顶角和方位角。

在激光波长405 nm，荧光波长500 nm时，对厚度30 μm，不同含油率的油包水乳化液进行fBRRDFcosθrcosθi仿真，其中含油率为60%时fBRRDFcosθrcosθi的三维分布如图3所示。并在文献[8]的基础上得到水包油乳化液高浓度下的fBRRDFcosθrcosθi，与油包水的对比结果如图4所示。根据油包水和水包油乳化液各自物理结构与分布等因素，对含油率分别为60%～90%、10%～40%的乳化液进行研究。

图3 油包水乳化液在含油率为60%时fBRRDFcosθrcosθi的三维分布图Fig.3 Three dimensional distribution map of fBRRDFcosθrcosθi in water-in-oil emulsion with oil content of 60%

图4 两种乳化液不同含油率的对比结果(a)： 油包水乳化液； (b)： 水包油乳化液； (c)： 45°接收角处各含油率的fBRRDFcosθrcosθi值Fig.4 Comparison of two emulsions with different oil contents(a)： Water-in-oil emulsion; (b)： Oil-in-water emulsion;(c)： Value of fBRRDFcosθrcosθi of each oil content at the acceptance angle of 45°

由图3可看出，油包水乳化液fBRRDFcosθrcosθi的形状大致关于中心对称，与光子出射的方位角无关，这是由于辐射荧光是各向同性的。由图4可知，两种乳化液的fBRRDFcosθrcosθi值均随含油率的升高呈下降趋势，且油包水的下降速率及fBRRDFcosθrcosθi的整体值较水包油的低。这主要是随乳化液含油率的升高，荧光量子产率下降导致的(见图2)。

为进一步验证仿真的正确性，利用实验室搭建的便捷式LIF系统对不同含油率的原油油包水进行荧光光谱测量，并与fBRRDFcosθrcosθi进行趋势对比。实验中激发光源由NDV4542激光二极管和电源组成，激发波长为405 nm，光谱仪采用AvaSpec-ULS2048型号光纤光谱仪。将Span80和Tween80两种乳化剂以10∶1比例分别与一定量的原油和水混合，再将两者混合倒入烧杯，利用高速搅拌机搅拌30分钟以制备出一定含油率的油包水乳化液。制备的样本含油率分别为60%，70%，80%，90%。在采集光谱过程中，固定激光二极管、光纤探头与样本间的角度、距离。选用Avasoft8软件采集光谱信号，波长范围为420～750 nm，积分时间为400 ms，每组样本进行多点多次采集光谱，将采集到的光谱扣除背景光后求得平均值作为最终的光谱数据如图5(a)所示，取各光谱波长500 nm处荧光与仿真进行对比如图5(b)所示。油包水乳化液的荧光强度随含油率的升高基本呈下降趋势，因此利用fBRRDFcosθrcosθi函数得到的油包水乳化液随含油率增长的变化趋势是正确的。

图5 实验数据与仿真对比(a)： 不同含油率下的荧光光谱； (b)： 趋势对比Fig.5 Comparison of experimental and simulation data(a)： Fluorescence spectra under different oil content;(b)： Trend comparison

2.2.2 不同入射接收角度的仿真与分析

为找到基于LIF系统探测海面油包水乳化液适宜的入射接收角度，必须探讨入射接收角度与fBRRDFcosθrcosθi的关系。由图3可知LIF系统接收的荧光功率与方位角无关，故对入射角分别为0°，25°，45°，65°和85°各接受角的fBRRDFcosθrcosθi进行仿真，结果如图6(a)所示。同时利用2.2.1所述的实验装置和材料采集不同入射角度处油包水的荧光光谱以验证仿真正确性。其他条件固定，采集激光入射角度分别为0°，20°，40°，60°，70°，80°时的荧光光谱，处理后的光谱数据如图6(b)所示，取波长为500 nm处荧光的实测值与仿真进行对比结果如图6(c)所示。

由图6(a)可知接收角θr一定时，fBRRDFcosθrcosθi随入射角θi的增大先平稳后减小，尤其当θi>65°时下降迅速。当入射角θi一定，fBRRDFcosθrcosθi值随接收角θr的增大逐渐减小，在θr=0°时获得最大值。由图6(b)和(c)发现实验中获得的不同入射角处的荧光光谱与仿真结果具有基本一致性的变化趋势。以上说明基于蒙特卡罗光子传输模型进行油包水乳化液的fBRRDFcosθrcosθi仿真是正确的，且可利用此仿真研究其他参数与fBRRDFcosθrcosθi的关系。

图6 仿真与实验数据及对比结果(a)： 不同接受角的fBRRDFcosθrcosθi； (b)： 不同入射角的荧光光谱； (c)： 趋势对比Fig.6 Comparison of simulation and experimental data(a)： fBRRDFcosθrcosθi at different acceptance angles; (b)： Fluorescence spectra at different incident angles; (c)： Trend comparison

2.2.3 不同厚度的仿真与分析

在以上实验基础上探讨油包水乳化液的厚度与fBRRDFcosθrcosθi的关系，选择1号重质、2号轻质原油油品为对象，对不同厚度的油包水乳化液进行仿真，结果如图7所示。

图7 两种油的油包水乳化液不同厚度的fBRRDFcosθrcosθi(a)： 1号油品； (b)： 2号油品； (c)： 45°接收角处不同厚度的fBRRDFcosθrcosθiFig.7 fBRRDFcosθrcosθi of different thicknesses of water-in-oil emulsions of two oil products(a)： No.1 oil product; (b)： No.2 oil product; (c)： fBRRDFcosθrcosθi with different thicknesses at a acceptance angle of 45°

由图7看出1号、2号的fBRRDFcosθrcosθi随乳化液厚度的增加均表现为先上升后饱和的趋势。这主要是随乳化液厚度的增加，虽会产生更多荧光光子，但由于吸收作用使其不足以通过很厚的介质被接收，造成fBRRDFcosθrcosθi逐渐趋于饱和。同时，由表1、表2及图2发现，1号乳化液的吸收系数较大，荧光量子产率较低，故导致1号的fBRRDFcosθrcosθi值更快饱和，且整体值低于2号。由此可通过fBRRDFcosθrcosθi得到海面油包水乳化液的最小厚度，并在一定程度上区分溢油类型。而在乳化液厚度未知的情况下，由于fBRRDFcosθrcosθi存在重叠区，故需结合新信号源进行溢油类型判别，这将作为下一步的工作重点。

3 结 论

结合蒙特卡罗方法和Mie散射理论，建立了海面溢油污染中油包水乳化液的BRRDF仿真模型，利用该模型得到基于LIF接收的荧光信号与探测各参数的关系。分析可知，海面油包水受激辐射的荧光信号随含油率的上升呈下降趋势。这可作为基于LIF探测评估海面溢油油包水乳化液含油率(浓度)的依据之一。在乳化溢油油包水的实际LIF探测中，入射接收角均不宜过大，入射角应小于65°，接收角为0°时可获得最佳荧光信号。由fBRRDFcosθrcosθi与油包水乳化液厚度的关系得知，油包水受激辐射的荧光信号随厚度的增加呈先上升后饱和的趋势，且重质油达到饱和点的最小厚度远小于轻质油。这可作为基于LIF探测海面溢油油包水乳化液最小厚度的判定参考，为评估溢油量提供数据源，同时为溢油油种的识别也提供了特征依据。

溢油进入海域，由于自身性质和海面风浪等环境的不断作用与影响，会经历复杂的风化过程，通常被划分为未乳化、乳化等阶段。本文工作仅探讨在气水界面水平条件下，基于LIF技术海面溢油油包水乳化液探测参数的双向反射再辐射分布函数。对未乳化油膜阶段以及海面风浪、波纹等复杂环境下海面溢油的BRRDF研究将是今后的工作重点。