海上石油平台泄漏、运输船舶事故等造成的海洋溢油污染往往波及海域范围大且影响时间长。大量原油进入海洋, 与海水混合极易形成油包水型乳液^[1-3]。乳化将改变油膜物理性质, 其体积能增大至原始溢油的2至5倍, 黏度至少增大3个数量级, 呈现出介于液体和固体间的状态^[4-5], 在海洋环境中很难自然消失, 加大了溢油回收清理工作实施的难度^[6-7]。因此, 及时监测溢油范围, 准确识别原油及乳化油膜, 能有效地制定海上溢油应急策略^[8-9], 对于保护海洋生态环境、减少溢油造成的经济损失具有重要意义。

大量原油乳化机制研究表明, 原油中的沥青含量是形成乳化的重要因素, Fingas等^[10]以沥青含量定义乳化物为稳定(沥青含量 > 7%), 半稳定(3% < 沥青含量 < 7%)以及不稳定(沥青含量 < 3%)。溢油形成乳化所需的动力能量不高, 在较平静的海面上由波浪所产生的的能量足矣。Fingas^[11]通过实验发现旋转法制备乳化物时所需最小能量约为3×10^–5~5×10^–5 J, 搅拌作用越剧烈乳化程度越大。温度对原油乳化有复杂的影响^[12-13]: 高温有利于油水混合, 可加速乳化物的形成, 而低温抑制水滴的凝聚, 利于乳化物的稳定; 温度影响沥青质的溶解与沉积状态, 因此有时较低温度更利于乳化的形成; 同时, 温度也通过影响蒸发快慢来影响溢油乳化的趋势。

目前SAR溢油监测研究主要集中于溢油与疑似溢油的识别, 涉及乳化进程对油膜散射特性影响方面的研究较少。SAR传感器具有高效性、实时性、不受云雨干扰的优势, 在海面溢油监测中发挥着主力军的作用。单极化SAR基于多尺度影像分割, 利用模糊逻辑、神经网络等方法实现溢油范围的半自动或自动识别^[14-18]。随着全极化SAR卫星的发展, 极化分解参数, 平均散射角、散射熵等参数被用于溢油与海水的区分^[19-22]。相关研究表明SAR数据具有识别较厚乳化油的潜力^[23-24]。阻尼比定义为无油海面与油膜散射强度之比, 能够反映油膜对海水强烈的平滑作用^{[8, 25-27]}, 在一些研究中用来对溢油区域内部进行分类^{[8, 28]}, 识别较厚油层以及混合海水的溢油^[28-29]。Minchew等^[29]认为SAR能够观测较厚浮油的体积浓度, 溢油与表层海水的混合物通过平滑海表毛细重力波、降低海水介电常数来改变海表后向散射, 并以此提出了一种基于布拉格散射理论的油水混合指数来估算溢油乳化物中油的体积分数, 并采用“深水地平线”钻井平台溢油事故中获取的L波段UAVSAR数据进行了测试评估。Skrunes等^[20]利用同步溢油实验获取含有原油、乳化油和生物油膜的RADARSAT2数据, 通过极化矩阵提取的特征向量获得的几何强度和双极化的后向散射截面实部的对数值, 成功识别了原油、乳化油和生物油膜。Garcia-Pineda等^[23]基于纹理分类器神经网络算法确定深水地平线溢油事故中溢油的泄漏范围, 通过开发的乳化油膜识别算法检测泄漏范围内雷达后向散射增强的区域, 从而识别乳化油膜, 识别结果在很大程度上与船只以及光学遥感影像观测结果一致。然而由于海上溢油事故期间对真实溢油进行采样测量乳化程度的难度较大, 上述研究对于乳化油膜乳化程度的定量观测缺乏验证。

SAR是监测海上溢油的主力军, 而乳化是贯穿海上原油溢油的重要过程。因此开展原油乳化外场实验、探索不同乳化程度对微波散射特性的影响, 对于乳化原油的识别具有重要意义。Guo等^[30]认为含水率是描述原油乳化进程的最佳参数, 并发现一定程度的乳化能够增大油膜后向散射系数(normalized radar cross section, N_RCS), 而原油乳化进程对N_RCS的影响未被定量描述。本文利用实验室制备的不同含水率的乳化油品模拟原油乳化进程, 开展C波段全极化散射计原油乳化观测外场实验, 研究原油乳化进程中N_RCS的变化, 为SAR监测乳化溢油提供数据支持。

1 材料和方法 1.1 实验材料

实验使用的油品有半稳定乳化原油A, 稳定乳化原油H、B以及不稳定乳化原油C, 其中A、B两种油采自石油平台油井(油、水及杂质混合), C、H是脱水原油。为定量分析乳化进程对原油N_RCS的影响, 实验制备了不同含水率的乳化油样品, 均以1 kg原油进行配制。根据预实验中不同原油最大稳定含水率(A: 20%; H: 70%; B: 50%; C: 22%)设置乳化原油的含水率间隔如表 1所示。

1.2 实验仪器

乳化原油散射特性实验设置如图 1所示。玻璃钢水槽(长6 m, 宽2.2 m, 高0.7 m)内盛4 700 kg密度为1.02 g/mL的过滤海水至池内0.35 m深度处; C波段全极化散射计置放于水槽短边高台处, 在入射角25°至60°范围内每隔5°扫描一次, 扫描足印直径0.8 m, 扫描中心距短边池壁距离如图 1所示。散射计具体参数如表 2所示。

1.3 实验方法及实验条件

本次研究共开展两次原油乳化实验Ⅰ、Ⅱ, 地点位于烟台牟平海岸带环境综合试验站。实验Ⅰ(2017年8月23—27日)观测C油的自然乳化进程。分9次向池中累加泼入原油, 累计加油量达6 996.7 g。每次加油后待油膜扩散10 min后、散射计进行扫描, 之后根据不同油量人工搅拌2~10 min不等, 观测油膜N_RCS。实验Ⅱ(2018年9月15—18日)利用人工制备的溢油乳化样品模拟原油的乳化进程, 观测不同含水率的乳化油A、H、B。散射计分别观测洁净海水与乳化原油。每个油样观测结束后对水面进行除油处理, 观测除油水面及下一油样N_RCS。两实验具体过程及实验油量如图 2所示, “初始海水”和“除油水面”分别代表初始洁净海水和除油后水面; “搅拌”表示该次观测的油膜经过人工搅拌; “放置一夜”表示该次观测的是放置一晚的乳化油膜。“A-5%”代表含水率为5%的乳化原油A。实验过程中气温与风速由距离水槽200 m的小型气象站记录。为探究玻璃钢水槽池壁对水面风速的影响, 同时利用手持风速计围绕水池两侧的长边池壁在池壁上方与接近水面位置测量了风速。图 3中红色及绿色实线分别表示气象站观测的气温与风速, 深蓝色实线表示池壁上方位置平均风速, 浅蓝色虚线表示近水面位置平均风速, 曲线下方的百分比代表了近水面风速相比池壁上方风速衰减的比率, 风速衰减率平均为12%。这表明池壁对外部的风有一定的遮挡作用, 但不严重。水面波高由测波杆利用海水的导电性, 通过电测方法测量测波杆浸泡于海水中的高度, 测量精度为0.3 mm。在散射计入射角40°扫描区域的左右两侧选择了4个位置进行测量, 观测结果显示实验水面波高在3 mm以内。

2 结果与分析

实验Ⅰ观测了原油的自然乳化过程, 含水率是能够界定其乳化状态的重要参数, 预实验测得C油在28 ℃下最大稳定含水率为22%。不同入射角观测的油膜N_RCS随乳化进程的变化如图 4所示, 图中序列号6显示油膜在温度28 ℃时达到最大稳定含水率(21.9%), 根据现场照片油膜主要分布于入射角35°至45°扫描范围内, 观测到35°至45°油膜N_RCS明显增加。对比风速结果显示风速对油膜的N_RCS变化并没有起主导作用, 不稳定原油的乳化则一定程度增大油膜N_RCS。因此, 乳化能够一定程度影响油膜后向散射机制, Guo等^[30]在相同实验条件下也发现一盯程度的乳化增大N_RCS。

进一步研究乳化程度对油膜N_RCS的影响, 实验Ⅱ观测了不同含水率的原油样品。油膜现场照片如图 5所示, 含水率越高油膜乳化程度越高, 油膜越接近棕褐色且油层中小液滴更密集。

内场实验结果表明温度较低的油膜有更大的密度、黏度与表面张力, 且密度、黏度与表面张力随含水率增大更快, 可见较低的温度更利于油膜乳化的稳定^[31]。同时, 温度还通过影响油膜的挥发量影响着溢油乳化的趋势。由于不同含水率油膜的观测时间相隔较大, 观测结果受环境差异影响。而实验中油膜与除油水面的观测间隔仅在0.5 h以内, 因此利用油膜与除油水面N_RCS的差值(△σ⁰)与阻尼比(D_R)分析乳化程度对N_RCS产生的影响。上述散射特征参数定义如下:

$ \Delta {\sigma^{0}}={\sigma^{0}}_{\mathrm{Oil}}-{\sigma^{0}}_{\text {Water }}, $

(1)

$ D_{\mathrm{R}}={\sigma^{0}}_{\text {Water }} / {\sigma^{0}}_{\text {Oil }}, $

(2)

其中, σ⁰_Oil为除油后水面的N_RCS, σ⁰_Oil为油膜N_RCS。由于布油点位于水面中心扫描入射角45°范围处, 以下分析均采用45°观测数据进行分析。

VV、HH、VH/HV极化方式下△σ⁰及D_R随含水率变化分别如图 6、7所示。对于H及B油, 结果显示△σ⁰及D_R随原油乳化程度增加, 含水率越高对应△σ⁰, D_R越大, 且这种趋势在VV极化下最为显著。对于A油, △σ⁰及D_R随含水率的变化趋势并不显著, 现场照片图 5a显示油膜A-0%呈大小不一的油块堆叠在水面, 增加了水面粗糙度, 使得油水NRCS差异较大, 这种差异并非乳化引起。而乳化油膜A-10%主要分布于扫描范围四周, 扫描范围内分布着平滑的甚薄油膜(图 5b), 造成较低的∆σ⁰_A-10%与D_RA-10%。VV、HH、HV/VH三种极化下∆σ⁰_A-15%均高于∆σ⁰_A-5%, VV极化下D_RA-15%大于D_RA-5%, 即乳化程度较大的油膜△σ⁰及D_R值更高。

因此, 参数△σ⁰、D_R能够反映油膜乳化程度的变化, 且在VV极化下效果最佳。进一步探究上述参数对乳化原油与未乳化原油的识别能力, 不同含水率的乳化原油相对未乳化原油△σ⁰与D_R的变化率分别为:

$ d_{\Delta \sigma 0}(\%)=\left(\Delta {\sigma^{0}}_{\text {emul }}-\Delta {\sigma^{0}}_{\text {crude }}\right) / \Delta {\sigma^{0}}_{\text {crude }} \times 100, $

(3)

$ d_{\mathrm{DR}}(\%)=\left(D_{\text {Remul }}-D_{\text {Rcrude }}\right) / D_{\text {Rcrude }} \times 100, $

(4)

式中, △σ⁰_emul及D_Remul为不同含水率乳化原油计算的△σ⁰与D_R值; △σ⁰_crude及D_Rcrude为未乳化原油计算的△σ⁰与D_R值。如表 3所示, 乳化油膜较未乳化原油△σ⁰与D_R值在VV极化下的变化率随含水率增大。对于沥青含量7.5%的原油B, 乳化油膜含水率达10%便可以通过油水散射差值与原油进行明显区分, 其△σ⁰较原油增加近一倍。含水率达最大稳定值40%时, △σ⁰较初始原油增加约2.2倍; 对于沥青含量 > 7%的脱水原油H, 乳化油膜含水率达40%时能够通过△σ⁰与原油明显区别, △σ⁰较原油增加1.5倍, 达最大稳定含水率50%时较初始原油增加2.2倍; 对于沥青含量介于3%至7%的A油, 乳化油膜含水率达最大稳定值15%后△σ⁰相对含水率5%的乳化油膜增加约1倍, 可以与未乳化原油进行区别。D_R对乳化油与未乳化原油的区分度相对较低, B油含水率高达40%时其D_R较原油增长12.3%, 含水率达50%的H油D_R较原油增长4.6%, 而含水率达15%的A油相对A-5%仅增加1.4%。

油膜后向散射主要取决于表面粗糙度及介电常数的影响^[32]。乳化导致油膜含水率与介电常数的变化, 乳化油层中的小液滴一定程度增加了油膜表面粗糙度(图 5a), 从而影响油膜后向散射机制, 使乳化原油N_RCS高于未乳化原油。散射特征参数△σ⁰与D_R能够增强乳化原油及未乳化原油的这种散射差异, 帮助实现原油及乳化油的区分识别。

3 结论

本文开展内场实验探索了原油乳化进程中物理特性的变化, 含水率能够定量描述原油的乳化过程。利用C波段全极化散射计分别观测了溢油的自然乳化进程以及人工控制下的定量原油乳化进程。在风速小于7 m/s, 波高 < 3 mm条件下, 研究乳化进程对油膜N_RCS的影响。实验观测结果表明乳化通过改变油膜的表面粗糙度与介电常数增大其N_RCS, 使得乳化原油N_RCS高于未乳化原油。油水差值△σ⁰与阻尼比D_R能够反映油膜乳化程度的变化, 油膜乳化程度越大, 其对应△σ⁰与D_R值越高。VV极化下, 油水散射差值△σ⁰能够较好地应用于乳化原油与未乳化原油识别, 沥青含量高的原油在乳化程度较低时便能与原油显著区别。在VV极化下, 阻尼比D_R能够对原油与乳化程度较高的乳化原油进行一定程度区别, 沥青含量越高区别度越显著。

目前实验还存有一些不足, 缺少实测的油膜表面粗糙度参数, 乳化进程增大N_RCS的理论机制需要进一步探索。实验需要进一步优化条件, 在更接近自然海况的条件下, 应用不同波段的微波传感器对比不同波段对原油乳化程度响应的敏感性, 不断实验分析, 为星载SAR识别乳化油膜提供数据参考。

致谢: 感谢中国科学院牟平海岸带环境综合试验站为本实验提供场地支持及气象站观测数据。