石油是全球重要能源之一, 在“一带一路”合作倡议的大背景下, 中国石油迎来了良好的发展契机, 与各国的能源合作渠道、合作模式呈多元化、多样化发展(余晓钟等, 2021), 而能源的进口主要通过海上运输与管道运输(于伟娜, 2012; 孙宁宁, 2019)。海上贸易的发展以及海上石油的开发、运输导致海上溢油污染风险增大, 及时准确地监测溢油对保护海洋环境和生态资源具有重要意义(刘朋, 2012)。卫星遥感技术以其大范围、高频次、同步、快速、长时序等观测优势, 为溢油的监测、识别及清理工作提供了极大的辅助(郑本昌, 2018; 蒋兴伟等, 2019; 沈亚峰等, 2020)。

海上油膜可分为天然来源和人为来源(吴传雯, 2014), 天然来源如浮游植物或鱼类产生的生物油膜等(Gade et al, 1998)¹⁸⁸⁵¹, 人为来源主要有(Yang et al, 2020)⁷¹⁶³: (1) 石油平台溢油, 如2010年4月美国墨西哥湾“深水地平线”钻井平台发生爆炸, 造成约2 300万加仑原油泄漏(张明亮, 2015); (2) 输油管道的损坏事故造成的溢油, 如2010年7月大连新港发生输油管线爆炸事故, 原油泄露量1 500 t, 污染海域430 km² (温艳萍等, 2013); (3) 大型油轮的碰撞事故等造成的溢油, 如2018年1月巴拿马籍油船“桑吉轮”与中国香港籍散货船在长江口以东发生碰撞, “桑吉轮”发生燃烧及爆炸, 造成136 000 t凝析油泄露(金戈等, 2019); (4) 游轮或货船非法排放燃油, 溢油类型主要为重柴油和轻柴油。溢油进入海洋环境, 会发生蒸发、光氧化、溶解、乳化、颗粒物质的吸附沉降以及微生物降解等风化过程(赵云英等, 1997)。乳化是重要的风化过程, 石油产品中的沥青质、胶质和蜡对乳化物的形成起作用, 乳化使油膜形成不同浓度油包水乳化物, 显著改变油的物理化学性质(严志宇等, 2002²; Guo et al, 2019)。

为了最大程度降低溢油灾害带来的环境影响, 需要采取快速应急反应处置溢油, 正确识别溢油种类对于应急处理策略的制定具有重要意义(陆应诚等, 2019)。在光学遥感方面, 油膜的光谱吸收特性是区分并识别不同油膜的重要依据(Lu et al, 2019), 高光谱因具有丰富的光谱信息而在油种识别方面发挥着重要作用(Yang et al, 2020)⁷¹⁶⁴, 许多学者基于油膜光谱特性提取进行了油膜类型的识别研究(Wettle et al, 2009; 谈爱玲等, 2011; Shi et al, 2018)。虽然光学传感器尤其是高光谱可以提供丰富的光谱信息, 但溢油事故通常发生在天气恶劣的环境, 限制了光学遥感在溢油事故中监测作用的发挥。

合成孔径雷达(synthetic aperture radar, SAR)以其全天时全天候的优势在海面溢油遥感监测过程中发挥着主力军的作用(许晨琪等, 2021)¹³。SAR依赖于海面粗糙度所形成的后向散射对海面进行成像, 海面油膜抑制海面的毛细重力波, 使海面粗糙度降低, 从而减弱后向散射, 在SAR图像上形成暗斑(李煜等, 2019; 陈韩, 2020)。油膜的抑制效应与油膜厚度、表面弹性及风速风向等因素相关(Zhang et al, 2015⁷¹; Zheng et al, 2021¹²)。基于此理论基础, Wismann等(1998)在北海的受控溢油实验中, 分别释放了重质燃料油和轻质燃料油, 发现在S、C、X和Ku波段重质燃料油的抑制作用大于轻质燃料油, L波段两者抑制比差异不大, 且重质燃料油抑制比随油膜厚度增加而增加。Gade等(1998)利用星载SIR-C/X-SAR图像对生物油膜和矿物油膜进行区分研究, 结果发现同一油膜的抑制比强烈依赖于风速, 相对于C波段和X波段, L波段下生物油膜的抑制作用更大。Jones等(2016)基于2015年挪威海岸溢油实验获取了L波段无人机载合成孔径雷达(uninhabited aerial vehicle synthetic aperture radar, UAVSAR)图像, 分析可知在释放后的较短时间内, 矿物油膜在不同极化方式下的抑制比均表现出明显的分带, 而植物油在释放后分区并不明显。此外, 基于极化分解和特征参数的油种区分也有研究, Tian等(2010)基于中国海南海上溢油实验获取了RADARSAT-2图像, 利用非相干目标分解结果对矿物油、生物油和似然物进行区分与识别研究。段冰等(2013)在分析生物油和矿物油极化比的基础上, 提出基于交叉极化比的生物油和矿物油的区分方法。Skrunes等(2014)获取了海上矿物油和生物油观测实验的全极化C波段合成孔径雷达(synthetic aperture radar, SAR)数据, 发现几何强度和同极化向量积的实部可以用于区别生物油膜和矿物油。

综上所述, 随着微波技术的快速发展, SAR在油种识别方面已存在极大的潜能。本文基于海面溢油来源的主要类型, 选取柴油(石油衍生品代表)、原油、油水混合物(乳化油代表)和棕榈油(生物油代表)四种油品, 开展C波段全极化散射计溢油观测实验, 探究微波识别油膜的敏感特征参数, 并将其应用于海上油膜实验获取的SAR图像进行油种识别。

1 实验与方法 1.1 外场实验 1.1.1 实验油品及实验仪器

基于海上溢油类型, 选择了四种油品进行实验, 分别为: (1) 柴油; (2) 原油, 其沥青含量为0.73%, 密度为0.92 g/mL; (3) 油水混合物; (4) 棕榈油, 其密度为0.85 g/mL。本次实验位于青岛市崂山区南姜码头陆基海水池(45 m× 40 m×2 m)设置的围油栏内(为防止油膜溢出, 实验设置双层围油栏), 围油栏区域大小为6.8 m×3.2 m, 美国ProSensing公司生产的C波段全极化散射计设置在池边钢板平台上, 极化方式为垂直(vertical transmission vertical reception, VV)极化、水平(horizontal transmission horizontal reception, HH)极化和交叉(vertical transmission horizontal reception/ horizontal transmission vertical reception, VH/HV)极化, 其在25°~60°入射角范围内以5°为间隔进行扫描, 扫描足印直径0.8 m, C波段全极化散射计参数如表 1所示, 实验设置如图 1所示, 不同入射角下的扫描中心与池边距离如表 2所示。

1.1.2 实验条件及过程

本次外场实验时间为2020年9月26~29日, 4天时间内分别进行柴油、原油、油水混合物、棕榈油的后向散射系数(normalized radar cross section, N_RCS)观测, 实验台附近设置小型风场、气温观测仪, 每次实验伴随测量海水、油膜温度。9月26日测量海水温度后首先使用C波段全极化散射计扫描海水的N_RCS, 分别得到四种极化方式下的N_RCS, 即σ_VV⁰, σ_HH⁰, σ_VH⁰和σ_HV⁰, 其中σ_VH⁰ = σ_HV⁰。为了减小测量误差, 散射计每次测量时扫描3次, 将3次扫描结果求均值得到最终的N_RCS, 记录相关数据; 然后分12次向围油栏内倒入柴油伴随油温测量, 油膜的N_RCS处理与海水的相同; 柴油实验完毕后清理油膜。其他油品的N_RCS测量与柴油相同, 为了使油膜尽可能均匀分布, 在倒入原油和油水混合物后进行了搅拌操作, 而后进行温度测量和散射计观测, 具体过程及实验油量如图 2所示。

1.2 油种识别特征参数选择 1.2.1 极化差

海面N_RCS表示为镜面散射、布拉格散射以及非布拉格散射三者之和, 中等入射角下, 镜面散射在总散射中占比小, 可以忽略不计, 因此, 可以将中等入射角下海面N_RCS表示为布拉格散射与非布拉格散射之和(Zheng et al, 2021)²。极化差(polarization difference, P_D)定义为VV极化的N_RCS与HH极化之差, 计算公式如式(1)所示, 其去除了非布拉格散射的影响, 主要包含接近布拉格波数的短波分量的贡献。P_D在清洁海面时值最大, 随着溢油影响增大而慢慢降低(Angelliaume et al, 2017)。

1.2.2 抑制比

油膜在SAR图像上形成暗斑, 原因之一为油膜抑制海面的毛细重力波, 使海面粗糙度降低, 从而减弱海面的后向散射。抑制比(damping ratio, D_R)可以定量描述油膜对海浪谱的抑制效应, 定义为无油海水与油膜覆盖海表N_RCS的比值(Wismann et al, 1998³⁶¹⁵; 许晨琪, 2021³⁴), 其计算公式如式(2)所示。

1.2.3 极化比

海水的介电常数大于油膜, 油膜的存在使得有效介电常数变小, 同样也会导致后向散射减小。极化比(polarization ratio, P_R)定义为HH极化与VV极化的N_RCS之比, 计算公式如式(3)所示, 在布拉格散射理论下, 极化比与海面粗糙度无关, 仅取决于入射角和海面的介电特性。当油膜厚度较小时(如生物油膜和薄矿物油膜), 微波会穿透油膜而返回海水的特征, 因此介电常数的变化可以忽略, 其极化比与无油海水的极化比基本相同; 当海面覆盖厚油(对C波段而言油膜厚度大于1 mm)或乳化油膜, 介电常数将影响极化比, 在同一入射角下, 随着含油率的增加(即含水率的减小), 极化比呈现增加趋势(Angelliaume et al, 2018); 同时通过观测极化比值与理论值的对比, 发现非布拉格散射对N_RCS产生重要影响(Skrunes et al, 2015)。

σ_VV⁰为VV极化的N_RCS, σ_HH⁰为HH极化的N_RCS, σ^{0, Oil}为油膜的N_RCS, σ^{0, Water}为海水的N_RCS。

2 结果与分析 2.1 南姜码头微波实验数据特征分析 2.1.1 海水和油膜N_RCS分析

在布拉格理论中, 随着入射角增大N_RCS呈现减小的趋势(Garcia-Pineda et al, 2020)³, 且在20°~70°入射角下, σ_VV⁰>σ_HH⁰> σ_VH/HV⁰ (Jones et al, 2018)。本次实验海水和柴油N_RCS随入射角的变化如图 3所示, σ_VV⁰、σ_HH⁰和σ_VH/HV⁰随着入射角的增大呈现减小的趋势; 入射角在25°、30°、55°和60°下的N_RCS呈现σ_VV⁰>σ_HH⁰>σ_VH/HV⁰, 35°~50°入射角下由于仪器或环境影响出现σ_HH⁰>σ_VV⁰在本文研究中不考虑。由图 1及表 2可知, 25°入射角与池壁距离较近, 60°入射角与池壁最远, 且油膜分布少, N_RCS易受环境因素影响, 因此本文选择30°和55°入射角下的数据进行研究。图 4a和4c为30°和55°入射角在不同实验序列下油膜的N_RCS与风速的对应变化图, 分析可知, 风速的变化对部分实验序列的N_RCS产生了明显的影响, 风速的增大或减小改变了水面油膜厚度的分布, 进而使N_RCS增大或减小, 如序列2~4对应油膜N_RCS随风速增大而增大, 序列16~18和25~27对应油膜N_RCS随风速增大而增大, 序列21~24对应油膜N_RCS与风速的变化一致性较高, 其他序列下的油膜N_RCS与风速的变化一致性并不显著; 当风速大于3 m/s时, 油膜N_RCS并未随风速的增大而增大到最大值, 当风速小于1 m/s时油膜N_RCS同样未出现最低值, 且在风速小于1 m/s时, 油膜N_RCS在VV极化方式下的降低更为明显。根据油膜N_RCS随温度变化图 4b和4d可知, 在序列18~22和序列26~30即原油和油水混合物实验中, 水温与对应N_RCS的变化趋势较为一致(风速整体呈下降趋势), 此时实验时间均处于每日温度较高的中午, 温度的升高加快乳化进程(严志宇等, 2002)³, 因此, N_RCS呈现增加趋势(Guo et al, 2020)⁴⁸。

油水差(N_RCS differences between oil-water, ∆σ⁰)可以反映油膜相对于周围海水的变化(许晨琪, 2021)³³, 其计算公式为

柴油、原油、油水混合物和棕榈油VV极化方式下的油水差如图 5所示, 在30°和55°入射角下对VV极化的数据进行分析发现: 油水混合物、原油和柴油均出现了油膜N_RCS大于海水即∆σ⁰>0的现象, 这与油膜抑制作用增强使得∆σ⁰<0的结论并不一致。在不同极化方式下对30°和55°入射角的∆σ⁰>0的数据进行统计如表 3所示, 可以发现油水混合物作为油膜乳化物在不同极化方式下均出现了油膜N_RCS增大的现象, 且其比例均大于50%, 说明一定程度乳化能够增大油膜的N_RCS, 从而使得∆σ⁰>0, 印证了Guo等(2020)⁴²的结论; 原油在三种极化方式下出现油膜N_RCS增大, 这是由于油膜出现部分乳化的现象(原油和油水混合物由于黏度大, 为保证扫描区域有油膜分布, 试验过程进行了人工搅拌); 柴油和棕榈油(实验过程没有进行人工搅拌)不含沥青, 不易发生油包水的乳化(严志宇等, 2002)², 但柴油在不同极化方式下存在大于50%的∆σ⁰>0且不同序列下∆σ⁰的值比较接近, 这是由于柴油进水后扩散较快, 黏度较小的柴油在风场的作用下向围油栏左长边聚集(面向水池), 散射计扫描轨迹中油膜分布不均匀且油量少, 致使∆σ⁰>0且差值大小比较接近。棕榈油在实验过程中风速在逐渐减小(图 4), 油膜在扫描轨迹上分布较均匀, 随着油量的增加N_RCS不断减小, 因此棕榈油∆σ⁰<0。由实验数据分析可知, 棕榈油∆σ⁰均小于0, 柴油的∆σ⁰变化小, 因此, 基于∆σ⁰可以将二者与原油和油水混合物区分, 原油和油水混合物由于乳化使得∆σ⁰均呈现增大趋势, 在VV极化方式下比较容易识别(表 3)(许晨琪等, 2021)¹⁶, 而在HH和VH/HV极化方式下较难将两者区分。如果在相同条件下比较, 油水混合物∆σ⁰>0的比例应远远高于原油。

2.1.2 油膜散射特征参数变化分析

由于棕榈油实验序列为6个, 因此, 基于特征参数的分析选取4种油品的前6个序列, 序列1为无油海水, 序列2~6为油膜特征参数变化, 该实验序列下的风速、水温如图 6所示。

图 7a, 7b为30°和55°入射角下柴油、原油、油水混合物和棕榈油的极化差P_D, 图 7c~7e和图 7f~7h分别为30°和55°入射角下四类油膜不同极化方式下以dB表达的N_RCS计算的抑制比D_R, 图 7i, 7j分别为30°和55°入射角下四类油膜的极化比P_R。随着油膜影响增大(油膜厚度相对增大), N_RCS减小, P_D和D_R均应呈现减小趋势, 而P_R则应呈现增加趋势。柴油在序列3~4下风速的增加使得油膜分散, 油膜对N_RCS的影响降低, 因此P_D和D_R异常增加, P_R减小; 而在序列4~6下随着油膜影响的增大, P_D和D_R呈现减小趋势, P_R呈现增加趋势。原油的P_D、D_R和P_R变化较小, 说明在此过程中原油的油膜厚度变化小且该入射角下为薄油分布, 序列3~5下的P_D和D_R逐渐增大, P_R逐渐减小, 这是由于风速的增加导致扫描区域油膜较薄, 其对N_RCS的影响减小, 序列6不同入射角的油膜分布不同, 因此, P_D和D_R出现不同的变化趋势。油水混合物的P_D和D_R呈现先减小后增大的趋势, P_R的变化趋势在两个入射角下并不相同; 在30°入射角下序列2~4的P_D随着油膜影响增大而减小, P_R随之增大, D_R则在序列3~4出现了异常增大, 这与乳化增大N_RCS有关; 序列4~6下的风速呈现增大趋势, 油膜分散使得其影响减小, N_RCS增加迅速使得P_D、D_R增加, P_R减小, 55°入射角下序列6的P_D存在较小的增加趋势, 而P_R均呈现增加趋势, 这可能与风速影响油膜分布而导致的异常有关。棕榈油随着油膜影响的增加, P_D和D_R呈现减小趋势, P_R呈现增加趋势, 30°入射角下的序列4和6风速减小, 油膜聚集使得油膜影响增大, 因而P_D和D_R在VV极化和VH/HV极化方式下异常减小, P_R异常增加; 序列5风速增加, 油膜分散使得P_D增加, VV极化和VH/HV极化方式下的D_R增加, P_R减小, 55°入射角下的变化正常; 实验过程中D_R在HH极化方式下出现了异常增加趋势, 且在30°入射角下的增加趋势更大, 这可能与风速影响油膜分布导致的数据异常有关。

综上所述, 柴油作为原油的衍生品, 其P_D、D_R两类散射特征参数值与原油接近; 棕榈油在四类油膜中的P_D和D_R特征参数最小并与其他油种差距较大, 最易识别; 油水混合物在P_D和D_R特征参数下的变化趋势为增大, 且变化趋势大于原油和柴油; 棕榈油的P_R呈现增加趋势, 其他油种则在不同入射角下的变化趋势不一致。对散射特征参数的均值进行统计, 其均值大小比较如图 8所示, 分析可知: (1) 不同极化方式的原油、油水混合物、柴油三类油膜的D_R均值均大于棕榈油, 可以有效识别棕榈油; VV极化和HH极化下柴油、原油和油水混合物的D_R均值接近, 而在VH/HV极化方式下三者D_R均值尤其是原油和油水混合物的差异较大, 因此VH/HV极化方式下对于原油和油水混合物的识别有效; (2) 柴油作为原油的衍生品, 其P_D与原油接近, 但原油、油水混合物和棕榈油的P_D均值依次减小, 因此P_D能识别原油、油水混合物和棕榈油; (3) 矿物油在不同入射角下的P_R变化趋势不相同, 且不同入射角下的均值大小关系相反, 因此在本实验中不适于四类油种的识别。

2.2 SAR图像油种识别应用 2.2.1 SAR验证图像获取

2011年6月6~9日, 挪威海洋清洁运营协会公司(Norwegian Clean Seas Association for Operating Companies, NOFO)在北海(59°59′ N, 2°27′ E)开展海上溢油实验, 实验区域如图 9所示。实验选取原油、乳化油和植物油三种油品, 原油是Balder油, 密度为0.914 g/mL; 乳化油由Oseberg blend原油混合5%的IFO380制成, 实验室测得含水率为69%; 植物油为Radiagreen ebo油, 具有与天然油膜类似的双亲结构, 用于模拟天然单分子生物成因油膜; 其他油膜特性可参考Skrunes等(2014)⁵³⁰⁴和舒思京(2020)。三种油膜的释放过程如表 4所示, 释放后对乳化油进行机械回收, 剩余1 m³以成像, 对原油进行分散处理。

本文获取了一景6月8日的RADARSAT-2图像进行油种识别应用, 其具体参数如表 5所示, 风速数据为参加实验的船只或附近的石油平台提供。经过定标、滤波处理后的VV极化N_RCS图像如图 10a所示, 基于VV极化图像的N_RCS利用阈值法协同目视解译的方法提取出油膜范围如图 10b所示。图中左上为植物油, 中间为乳化油, 右下为原油, 三类油膜清晰可见, 图像获取时间分别为植物油释放后13 h, 乳化油释放后29 h, 原油释放后9 h。

2.2.2 基于敏感特征参数的识别

南姜码头实验结果分析可知, 油水差∆σ⁰、抑制比D_R和极化差P_D可以应用于图 10油种的识别。基于ENVI 5.3软件对经过定标、滤波后的图像选取溢油邻近海水, 利用MATLAB R2018b按列计算海水样本的N_RCS均值, 而后计算油膜与对应列海水均值的差值形成∆σ⁰图像, 如图 11所示; VV和HH极化方式的图像显示植物油的∆σ⁰比原油和乳化油小, 且原油、乳化油油膜边缘∆σ⁰值大于油膜中间, 植物油不显示此差异, 这是由于原油、乳化油边缘油膜厚度小, 易于发生乳化, 使得边缘∆σ⁰值大, 这一结论与Guo等(2022)¹³一致, 植物油不发生乳化, 因此油膜边缘与中间∆σ⁰值无明显差异; VH/HV极化方式下原油、乳化油和植物油的∆σ⁰值分布并无明显差异。

与∆σ⁰计算方法类似, 基于dB表达的N_RCS计算不同极化方式下的D_R (N_RCS[dB])如图 12a~12c所示, 基于linear units表达的N_RCS计算的D_R (N_{RCS[linear units]})如图 12d~12f所示; VV极化和HH极化方式的D_R (N_RCS[dB])显示原油、乳化油边缘值大于油膜中间, 植物油不显示此差异, 这是由于原油、乳化油边缘油膜厚度小, 易于发生乳化, 使得边缘D_R (N_RCS[dB])值大, 这一结论与Guo等(2022)¹³一致, VH/HV极化方式下原油、乳化油和植物油的D_R (N_RCS[dB])值分布并无明显差异。VV极化和HH极化方式的D_R (N_{RCS[linear units]})显示植物油小于原油、乳化油, 矿物油膜边缘值明显小于油膜内部, 植物油的图像较为均匀; VH/HV极化方式下的图像可以看出乳化油的D_R均值大于原油, 两者存在差别。D_R (N_RCS[dB])和D_R (N_{RCS[linear units]})均可以在VV和HH极化方式下较好地识别植物油, 但D_R (N_{RCS[linear units]})在VV和HH极化方式下植物油和矿物油的差异更明显, 且VH/HV极化方式下可以识别原油和乳化油, 因此, D_R (N_{RCS[linear units]})更有利于油种识别。

基于ENVI 5.3计算的P_D图像如图 13a所示, 可以看出原油和乳化油的边缘P_D值大于油膜中间, 植物油的边缘和中间的差异不如矿物油明显, 但在油膜左侧也出现了边缘值大的情况, 这是由于边缘油膜厚度小, 油膜影响小, 进而导致边缘P_D大; 为了进一步增大油种之间的差异, 尝试计算油水差的极化差(polarization difference of oil-water difference, P_DOW), 即VV极化与HH极化油水差的差值, 其计算公式如式(5)所示, 计算结果如图 13b所示, 由图 13b可知, 基于P_DOW的特征图像可以看出原油、乳化油和植物油三类油膜的差异。

3 讨论

南姜码头散射计实验发现VH/HV极化在区分识别矿物油膜方面有优势。星载SAR图像的VH/HV极化方式下受噪声影响较大(Garcia-Pineda et al, 2020)⁴, 油膜抑制毛细重力波使得N_RCS更小, 因此, VH/HV极化在油种识别方面不敏感, 但在对NOFO溢油实验获取的图像计算特征参数时发现, VH/HV极化方式下的D_R (N_{RCS[linear units]})有利于原油和乳化油的识别。

油膜的抑制效应与油膜厚度、表面张力、黏度、风速风向等有关(Zhang et al, 2015)⁷¹, 此外, 一定程度乳化增大N_RCS, 这对油膜的抑制效应也会产生影响。Guo等(2022)¹³提出D_R (N_RCS[dB])能够一定程度反映原油与乳化原油混合条件下的相对油膜厚度大小, 其值越小厚度越大; 同时D_R (N_RCS[dB])能够反映SAR图像中同一油膜的相对乳化程度, 乳化程度越高值越大。Zhang等(2015)提出油膜厚度小于0.01 mm时, D_R (N_{RCS[linear units]})变化不大, 厚度大于0.01 mm时, D_R (N_{RCS[linear units]})总体上随厚度增加而增加。随着油膜影响增大(油膜厚度的增大), P_D呈现减小趋势, 而P_R呈现增加趋势, 因此, 结合各参数的变化特征可以有效判断南姜码头实验中的相对油膜厚度大小。由NOFO获取的SAR图像计算原油和乳化油的D_R图像可以看出油膜边缘油膜厚度小, 内部厚度大, 存在明显分区, 同时油膜边缘乳化作用较强, 对D_R的变化也会产生影响; 植物油厚度小, 且不发生乳化, 因此D_R的变化较小。乳化油释放距成像时间29 h, 且进行了回收操作, 油膜厚度较小, 原油释放9 h, 在海水的作用下发生乳化; 因此, 原油和乳化油的特征参数图显示相似(Skrunes et al, 2014)⁵³¹⁷。

4 结论

本文开展基于C波段全极化散射计的外场实验, 获取了柴油、原油、油水混合物和棕榈油的不同实验序列下的N_RCS, 探索可以有效识别不同油种的敏感特征参数, 并将外场实验获取的敏感特征参数运用于SAR图像中进行不同油种的识别。结果发现: (1) VV和HH极化方式下的油水差∆σ⁰有利于识别植物油和矿物油; (2) VV和HH极化方式下的抑制比D_R (N_RCS[dB])和D_R (N_{RCS[linear units]})均有利于识别植物油和矿物油, 但D_R (N_{RCS[linear units]})可以更清晰地显示植物油和矿物油的差异, 且在VH/HV极化方式下的D_R (N_{RCS[linear units]})对于原油和乳化油的识别更有效, 因此, 基于linear units的D_R更有利于识别油种; (3) 极化差P_D在南姜码头实验中有效识别了原油、油水混合物和棕榈油, 在SAR图像中发现油水差的极化差P_DOW可以有效识别三类油种。

未来还需要进一步控制实验条件, 获取其他油膜的数据, 增加其他特征参数, 探索具有普适性的油种区分识别方法, 为星载SAR油种区分与识别提供数据参考和方法支撑。

致谢 感谢加拿大国家航天局RADARSAT-2数据支持。