作者简介: 崔永强, 1979年生, 燕山大学计算机科学与技术博士研究生 e-mail: hbyd_yqc@sina.com
油膜厚度是海面溢油污染评估分析的一个重要指标, 激光诱导荧光(LIF)技术是目前最有效的海面溢油探测技术之一, 基于LIF探测技术的油膜厚度反演算法当下仅有适用于薄油膜(≤10~20 μm)的评估方法, 而对于较厚油膜(>20 μm)的评估目前尚无有效的反演算法。 鉴于此, 提出一种基于LIF技术适用于较厚油膜的反演算法, 该算法采用油膜荧光信号反演油膜厚度, 推导了油膜厚度反演公式, 并给出了基于该反演算法的油膜厚度评估方法。 首先采用最大类间方差算法(Otsu)选取合适的荧光光谱波段, 然后根据选取波段内每个波长的光谱数据反演油膜厚度, 最后采用反演油膜厚度的平均值作为油膜厚度评估结果。 研究了该算法的适用范围, 给出了该算法有效评估范围最大值与测量相对误差的关系, 并结合消光系数给出了在多种测量误差条件下不同消光系数油品有效评估范围的最大值。 通过实验对本文方法进行了验证, 选用原油和白油的混合油(1:50)作为实验油品, 以波长为405 nm的激光作为激发光源, 采集波长范围为420~750 nm, 采集了海水背景荧光和拉曼散射光光谱、 实验油品的荧光特征光谱和多种不同厚度的较厚油膜的荧光光谱。 采用Otsu算法选取420~476 nm波段评估油膜厚度, 在实验油品油膜厚度≤800 μm时, 该算法对油膜厚度的评估具有较高的精度, 平均误差为10.5%; 在油膜厚度>800 μm时, 平均误差为28.8%, 评估误差较大且随油膜厚度的增加快速变大, 该实验结果与利用测量相对误差和消光系数的分析结果一致。 实验结果表明, 该方法可以实现对海面较厚油膜厚度的有效评估, 并可以根据测量相对误差和消光系数判断评估结果的有效性。
Oil film thickness is an important indicator for the assessment and analysis of oil spill pollution on the sea surface. Laser-Induced Fluorescence (LIF) is one of the most effective technologies for oil spill detection at present, the oil film thickness inversion algorithm based on LIF is only suitable for thin oil film (≤10~20 μm), there is no effective method for the evaluation of thicker oil film (>20 μm). In view of this, an inversion algorithm based on LIF detection technology for evaluating thicker oil film is proposed, the algorithm uses oil film fluorescence signal to invert oil film thickness, deduces oil film thickness inversion formula, and gives oil film thickness evaluation method based on the inversion algorithm. First, Otsu algorithm is used to select the appropriate fluorescence spectrum band, and then the oil film thickness is retrieved according to the spectrum data of each wavelength in the selected band, finally, the average value of the retrieved oil film thickness is used as the oil film thickness evaluation result. The applicable range of the algorithm is studied, the relationship between the maximum value of the effective evaluation range of the algorithm and the relative measurement error is given, the maximum value of effective evaluation range of oil with different extinction coefficients under various measurement errors is given. The method in this paper is verified by experiments. The mixture of crude oil and mineral oil (1:50) is selected as the experimental oil and the laser with a wavelength of 405 nm is used as the excitation source. The collection wavelength range is 420~750 nm. The background fluorescence and Raman scattering spectra of sea water, the fluorescence spectra of experimental oil and various thick oil films are collected to invert the oil film thickness. Otsu algorithm is used to select the band of 420~476 nm to evaluate the oil film thickness. The evaluation results show that when the oil film thickness is ≤800 μm, the algorithm has high accuracy, with an average error of 10.5%; when the thickness is >800 μm, the average error is 28.8%, with a large evaluation error and rapidly increases with the increase of the oil film thickness. The analysis results of relative error and extinction coefficient are consistent with the experimental results. The results show that the method can effectively evaluate the thickness of the thick oil film on the sea surface, and judge the effectiveness of the evaluation results according to the measurement relative error and extinction coefficient.
随着海上石油运输、 开采的快速发展, 轮船排污、 运油船只及海上钻井平台泄漏等溢油事故的发生, 海上溢油污染问题已经十分突出[1, 2]。 及时获取海水表面溢油信息, 鉴别溢油种类, 估算溢油范围及污染量, 对溢油污染的快速处理和水体环境的保护具有重要意义[3]。
近年来多种遥感技术被用于溢油监测, 其中激光诱导荧光(laser induced fluorescence, LIF)探测技术是目前被认为最有效的海面溢油探测技术之一[4, 5]。 研究人员可以利用LIF技术鉴别海面溢油的污染情况与溢油种类[6, 7, 8, 9], 但对于油膜厚度的监测目前尚未形成一套可以有效适用于各种不同厚度油膜的评估方法。 因此, 利用LIF探测技术评估海面溢油范围及其污染程度已成为近期该领域的一个研究热点。
对于海面油膜厚度的评估, 传统方法是根据《波恩协议》油膜色彩与油膜厚度的对应关系人工确定油膜厚度。 Kung等提出了利用LIF探测技术评估海面溢油的建模方法[10], 1980年, Hoge等基于该方法提出了一种利用拉曼散射光评估油膜厚度的积分反演算法并进行了溢油探测实验[11]。 我国对海面溢油也进行了大量的研究, 中国海洋大学成功研制了用于海洋参数测量的激光雷达系统并于2006年完成了水面油膜的探测实验[12]。 上述研究仅适用于薄油膜厚度评估, 目前尚无适用于较厚油膜的反演算法。 鉴于此, 本文提出了一种基于LIF技术采用油膜荧光信号评估较厚油膜的反演算法, 研究了该算法的适用范围, 并通过实验对该算法的有效性进行了验证。
如前所述, Hoge等提出了一种适用于评估薄油膜厚度的积分反演算法, 其原理是海水受紫外激光束照射激发的拉曼散射光由于海面油膜的吸收呈指数衰减, 据此依据拉曼散射光衰减程度反演出油膜厚度, 油膜厚度d可表示为[11]
式中: ke和k分别为油膜在激发波长和波长λ 的消光系数; R'为油膜覆盖海面时海水的拉曼峰信号强度; R为纯净海水的拉曼峰信号强度。 其原理如图1所示[11]。
积分拉曼法依赖拉曼散射光信号, 当海面油膜较厚, 拉曼散射光被油膜吸收而无法监测到拉曼散射光信号时, 该算法无法反演油膜厚度, 因此该算法仅适用于具有拉曼散射光信号的薄油膜(≤ 10~20 μ m), 而不适用于较厚油膜的评估。
Kung等提出了利用LIF探测技术对海面溢油建模的方法, 采用LIF技术探测溢油覆盖的海面, 探测器接收到波长λ 的信号K包括油膜荧光信号、 海水背景荧光信号和海水拉曼散射光信号三部分, 该信号可以表示为[10]
式中, P0为入射的激光能量, ke和k分别为油膜在激发波长和波长λ 的消光系数, d为油膜厚度, η , ξ 和Ψ 分别为在波长λ 的油膜荧光转换效率、 海水荧光转换系数和海水拉曼转换系数, δ r为δ 函数。
由式(2)可知, 随着油膜厚度的增加, 海水的拉曼散射光信号由于油膜的吸收呈指数衰减, 而油膜荧光信号呈指数增强, 在无法监测到拉曼散射光信号时, 可以采用油膜荧光信号反演较厚油膜的厚度。
采用LIF技术探测无溢油覆盖的海面, 此时探测器接收到的信号为海水的背景荧光和拉曼散射光信号, 波长λ 的信号Kw可表示为
采用激光照射被测油品, 可以获取该油品的荧光特征光谱, 此时探测器接收到波长λ 的荧光信号Ko可表示为
将式(3)和式(4)代入式(2), 化简可得
式(5)化简可得油膜厚度的反演公式
式(6)中, Kw为海水的背景荧光和拉曼散射光在波长λ 的信号强度, Ko为油品在波长λ 的荧光信号强度, 因此该算法需要获取该油品的荧光特征光谱、 海水的背景荧光和拉曼散射光光谱以及在激发波长和波长λ 处的消光系数。 如文献[11]所述, 采用LIF技术反演油膜厚度需要利用油品的一些物理参数, 本文假定油品的消光系数和荧光转换效率已经在实验室获得。
常见油品的荧光特征峰波长一般在430~520 nm范围内, 因此采用波峰光谱信号即可反演油膜厚度, 考虑到现有探测器荧光接收设备具有较高的光谱分辨率, 仅采用单个波长数据反演油膜厚度会导致较大的误差, 因此选择一组波段, 对其中的每个波长分别反演油膜厚度, 采用反演结果的平均值作为油膜厚度评估结果。
由式(6)可知, Ko-K与油膜厚度d为指数关系, 选择Ko-K的值显著的波段反演油膜厚度可以有效降低测量误差的影响, 本文采用最大类间方差算法(Otsu)选取波段, 将Ko-K转换为0~255的整数值, 采用Otsu算法计算阈值, 筛选出Ko-K的值显著的波段评估油膜厚度。
设选取的波段为D={λ 1, λ 2, λ 3, … , λ n}, 油膜厚度评估方法如式(7)所示
式(7)中,
海面油膜厚度的有效评估范围与测量误差和油品的消光系数有关, 因此本文从测量误差和消光系数两个方面分析评估范围。
1.4.1 测量误差对评估范围的影响
由于设备精度、 实验环境等因素的影响, 采集到的数据存在一定的测量误差, 当油膜荧光强度趋于饱和时, 测量误差对反演结果有较大的影响, 因此该算法能够有效评估油膜厚度的范围与测量误差相关。 一般情况下, 海水的背景荧光信号和拉曼散射光信号Kw远小于较厚油膜的探测信号K, Kw的测量误差对反演计算的结果影响很小, 因此本文仅分析K的测量误差对反演计算的影响。
由式(6)化简可得
设
由式(9)可知, 存在值R, 当
由图2可见, 当油膜荧光信号趋于饱和时, 随着测量误差的增加R值不断降低。 在测量误差为1%的情况下, R=0.95, 即
1.4.2 消光系数对评估范围的影响
由式(8)可知, ke+k与油膜厚度d为反比关系。 如前所述Kw远小于K, 因此忽略Kw的影响, 可得在不同测量误差情况下ke+k与有效反演油膜厚度最大值的关系, 如图3所示。
![]() | 图3 不同测量误差下ke+k与有效反演厚度最大值的关系Fig.3 The relationship between ke+k and maximum inversion thickness with different errors |
由图3可见, 在相同测量误差的影响下, ke+k值越大, 能够有效反演油膜厚度的最大值越小。 在测量误差为1%的情况下, ke+k=20的油品, 能够有效反演的油膜厚度在1 500 μ m左右; ke+k=50的油品, 能够有效反演的油膜厚度在600 μ m左右; ke+k=100的油品, 能够有效反演的油膜厚度在300 μ m左右。 一般情况下, 重油消光系数较大, 中质油次之, 轻质油较小, 可以根据消光系数确定不同种类的油品有效反演厚度的最大值。 因此, 轻质油有效反演厚度的最大值要大于1 000 μ m, 重油有效反演厚度的最大值约为300 μ m。
对于荧光效率较低的油品, 不满足Kw远小于K的条件, 此时能够有效反演油膜厚度的最大值要小于该分析结果。
实验装置和材料包括: AvaSpec-ULS2048光谱仪、 NDV4542激光二极管、 手柄支架、 500 mL烧杯和移液器, 实验油品为原油和白油的混合油(1:50), 实验海水采用渤海海水。 实验装置如图4所示。
由式(8)可知, 以海水的背景荧光和拉曼散射光作为背景光信号, 采集到的不同厚度水面油膜和实验油品的荧光信号即为K-Kw和Ko-Kw, 因此设计实验过程如下:
(1) 在烧杯中放置500 mL海水, 采集海水的背景荧光和拉曼散射光光谱, 并以此作为背景采集其他光谱数据。
(2) 在烧杯中放置足够厚油膜(10 000 μ m), 采集实验油品的荧光特征光谱。
(3) 在烧杯中放置500 mL海水, 在水面分别放置25, 50, 100, 200, 400, 600, 800, 1 000, 1 200和1 400 μ m厚度油膜, 静置约30 min, 待油膜扩散均匀后, 采集样品荧光光谱。
光谱采集过程中, 保持探头与样本距离、 角度不变, 采用Avasoft8软件采集光谱信号, 积分时间为200 ms, 采集波长范围为420~750 nm, 每次采集30组数据, 采用平均值作为光谱数据。
图5给出了实验采集的海水背景荧光和拉曼散射光光谱、 实验油品的荧光特征光谱、 厚度为200 μ m油膜荧光光谱的30次采集光谱以及不同厚度油膜的荧光光谱。
由图5(d)可见, 荧光强度随着油膜厚度的增加不断增加, 但增加值逐渐变小, 最终趋于稳定。 实验油品的主波峰在456 nm附近, 在420~500 nm附近波段, 不同厚度油膜的荧光强度区别明显, 而在500~750 nm波段, 不同厚度油膜的荧光光谱有较多的重叠, 荧光强度不能有效区分不同厚度的油膜。
对不同厚度油膜的荧光光谱计算各波段Ko-K的值, 采用Otsu算法选择差异显著的波段, 结果如表1所示。
![]() | 表1 不同厚度油膜的波长选取范围 Table 1 Wavelength selection of oil films with different thickness |
由表1可见, 在420~476 nm波段范围内, 不同厚度的油膜都具有较大的差异, 因此, 采用该波段范围进行油膜厚度评估。
实验室获取该实验油品的消光系数, 在420~480 nm波段ke+k的值如图6所示。
利用实验获取的油品特征光谱、 海水荧光和拉曼光谱以及实验油品在激发波长和波长λ 的消光系数, 采用所选波段评估油膜厚度, 评估结果如表2所示。
![]() | 表2 不同厚度油膜的评估结果 Table 2 Estimate results of oil films with different thickness |
油膜厚度的评估结果与实际厚度如图7所示。
由表2可见, 当油膜厚度≤ 800 μ m时, 该算法具有较好的评估精度, 平均误差为10.5%; 当油膜厚度> 800 μ m时, 评估误差较大, 并随着油膜厚度的增加而快速变大。
利用测量误差和实验油品的消光系数, 可以分析得到能有效反演的数据结果。
![]() | 表3 不同厚度油膜 |
根据不同厚度油膜的多次采集数据, 估算测量相对误差, 计算方法如式(10)所示
式(10)中, Kmax和Kmin分别为测量的最大值和最小值, 由式(10)估算测量相对误差如表4所示。
![]() | 表4 不同厚度油膜测量误差 Table 4 Error of oil films with different thickness |
由表4可见, 本次实验测量误差在2%范围内, 在
由图6可知, 实验油品在评估波段的ke+k的值约为34 cm-1, 采用油品消光系数分析有效反演厚度, 由图3可得本次实验有效反演厚度最大值约为780 μ m。
对比表2的实验结果, 可知采用测量相对误差和油品消光系数的分析结果, 与本次实验结果基本一致: 油膜厚度≤ 800 μ m时该算法具有较好的评估精度、 油膜厚度> 800 μ m时评估误差较大。
随着荧光信号的增加, 油膜厚度分辨率逐渐降低。 采用实验油品的主波峰波段光谱信号可得, 在25~50 μ m厚度范围内, 分辨率约为0.01 μ m; 在600~800 μ m厚度范围内, 分辨率约为0.23 μ m。
提出一种基于LIF探测技术利用油膜荧光强度评估较厚油膜厚度的反演算法, 采用Otsu算法选择光谱波段进行厚度评估, 并对油膜厚度的有效评估范围进行了分析, 分析了有效评估范围的最大值与测量误差和消光系数的关系, 并给出了评估油膜厚度最大值的估算方法。 采用原油和白油的混合油(1:50)作为实验油品对该算法进行了验证, 实验结果表明, 该算法可以有效评估较厚油膜厚度, 实验结果与分析一致。 本文计算方法并未考虑探测高度、 角度以及探测环境如波浪等因素的影响, 当在高空进行探测时, 需要保持探测角度不变, 同时对接收信号根据探测高度进行校正。
溢油污染受风、 浪等环境因素的影响, 会发生扩散、 乳化、 溶解等一系列物理化学变化, 形成水面浮油膜和悬浮水体中的油, 该算法仅适用于海面油膜的厚度评估, 对溢油乳化过程的评估将是下一步工作的重点。
[1] |
|
[2] |
|
[3] |
|
[4] |
|
[5] |
|
[6] |
|
[7] |
|
[8] |
|
[9] |
|
[10] |
|
[11] |
|
[12] |
|