船运铁矿粉污染海冰反射光谱特征研究

刘丙新，郭 刚，吴东来，刘成玉，谢 锋

1.大连海事大学航海学院，辽宁 大连 116026 2.中国科学院上海技术物理研究所，上海 200083

引 言

环渤海地区是我国最具经济发展活力的地区之一，有秦皇岛港、黄骅港、京唐港、曹妃甸港和天津港等大型港口，每年都有大量铁矿运输船舶进出[1]。根据MARPOL关于防止船舶垃圾污染的附件V，国际海事组织引入了新的分类标准，以便能够识别对海洋环境有害的物质(HME)[2]，对于含有来自HME残留物的干燥残留物和/或洗涤水必须在适当的港口接收设施中排放。一些研究数据表明铁矿在港口船舶转运过程中，产生的铁矿粉、铁渣等会对海洋环境带了影响[3-5]。

目前，对海冰表面船运铁矿粉的污染监测研究较少。澳大利亚国家研究机构联邦科学与工业研究组织利用光学方法评估港口附近粉尘的成分及比例以及他们的大小分布，但对于铁尘沉降的范围和数量没有相关的研究[6]。利用遥感技术对海冰表面的铁矿粉污染情况进行监测，能够获取其面积、范围等信息，是最具潜力的方法。应用光学遥感进行海冰信息提取[7-8]、海洋溢油监测[9]、海冰中塑料污染[10]以及土壤中物质反演都有较多研究[11]，但海面(海冰、海水)的矿物粉污染研究仍然较少。

本文旨在提取铁矿粉与海冰的光谱反射特征，找出有效识别海冰表面铁矿粉污染的波段范围，并探讨光谱特征与铁矿粉颗粒几何截面积比例的相关性，建立海冰表面铁矿粉面积比例反演模型，为海冰受铁矿粉污染的范围以及沉降量的估算提供参考和依据。

1 实验部分

1.1 实验环境

选取位于中国渤海湾海域进行实地测量，研究区内的海冰是由一年生海冰组成，平均厚度为28 cm。为了减少反射阳光对实验结果的干扰，实验者穿着深色衣服。该实验是在晴天，无云的天气条件下进行的，以确保充足和稳定的照明，现场温度为-11 ℃。

1.2 材料

实验所使用的铁矿粉为来自鞍钢鲅鱼圈分公司生产使用的铁矿石(Fe2O3)，经过粉碎后，用细筛筛选，所使用的筛子为100目。将铁矿粉放到漏网里洒到海冰表面，尽量接近实际冰区的铁矿粉洒落过程，共有109个样区，每个样区洒落的铁矿粉面积比例见表1。

表1 样本面积分布情况Table 1 Area fraction distribution of the observed samples

1.3 仪器及光谱采集

光谱测量设备为ASD FieldSpec®3地物光谱仪，光谱范围为350～2 500 nm，重采样后光谱分辨率1 nm。测量时FieldSpec 3光谱仪的视场角为25°，观测天顶角为0°。探头距离冰面约为45 cm，那么，探头实际覆盖的冰面范围是一个半径约为10 cm的圆。圆内由裸冰和铁矿粉覆盖，对于不同样本来说，洒在冰面的铁矿粉面积比例不同。光谱测量时，先测量参考白板的辐射，然后顺次测量11个样本的光谱，每个样本大约需要40 s左右，最后测量参考白板的辐射。重复测量三次，整个实验进行了约28 min，在此过程，太阳高度角变化小于2°。为了避免铁矿粉对海洋环境的污染，整个实验结束之后，将铁矿粉收集并带回实验室。

1.4 数据处理

在野外反射率光谱测量中，部分近红外和短波红外波段会受到水汽强烈的吸收作用，所获取数据的噪声很大，所以需要去除这些噪声较大的波段，最终选取350～1 300，1 420～1 780和2 050～2 300 nm三个区域，每一条反射率光谱曲线共计1 563个数据。将所选取的反射率光谱曲线用进行均值滤波，去除噪声，窗口大小设置为5。

1.5 端元提取

为了尽可能克服不同厚度和不同内部结构的海冰对海冰-铁矿粉像元的提取的干扰，考虑使用光谱曲线的形状来提取铁矿粉-海冰像元。

(1)光谱矢量角余弦

选取若干个铁矿粉、海冰(包括不同厚度海冰)反射率曲线的特征波长(包括吸收、反射、跃升等)。将选取的波长反射率作为矢量，每个波段为矢量的一个维度，计算待提取像元与铁矿粉光谱的光谱矢量角。根据阈值，判定是否为铁矿粉-海冰像元。矢量角的计算式为

(1)

其中，Acos为光谱矢量角余弦值；λi为第i波段的波长；N为波段数；ρc(λi)为铁矿粉在λi的反射率；ρ(λi)为待定像元在λi的反射率；|·|为对矢量求模。

(2)光谱吸收指数

海冰在短波红外波段反射率很小，也可以考虑单独使用短波红外来提取铁矿粉-海冰像元，海冰在1 460 nm处有一个吸收谷，选取1 425，1 550和1 460 nm三个波段计算，1 460 nm波段的吸收指数SAI(spectral absorption index)，计算式为

(2)

1.6 混合像元分解

线性混合光谱模型认为混合像元的光谱吸收、反射强度等于内部各端元的光谱吸收、反射强度的光谱矢量与各端元面积系数矢量的乘积。

(3)

式(3)中，x为观测的混合光谱向量；S为端元光谱向量；a为丰度(即面积比例系数)；w为噪声光谱向量(即误差项)；M为端元数量；k为端元序号。

求解各端元面积比例系数即：立足于混合光谱和端元光谱，通过方程组解算得到结果。但在实际自然环境中，线型混合模型受到两个约束条件的限制

(4)

2 结果与讨论

2.1 反射光谱特征

图1为根据覆盖不同面积比例铁矿粉的海冰反射率光谱曲线图，纯净海冰的反射率有两个大的反射峰(590，1 090 nm左右)，一个小的反射峰(802 nm左右)，两个吸收相对较强的谷(805，1 015 nm左右)。在1 220 nm之后，海冰的反射率基本呈现一条缓慢下降的直线，在1 500～1 780和2 050～2 300 nm两个波段，其反射率不足0.01。铁矿粉的反射率波变化也较大，有三个吸收峰(598，720和1 110 nm)，在885 nm附近有一个很强的吸收谷，在1 460 nm附近有一个较小的吸收谷，铁矿粉的最大反射率为0.193，最小值为0.048。9个样本的反射率曲线则基本介于海冰和铁矿粉的反射率之间。

图1 不同面积比例铁矿粉覆盖海冰反射率光谱曲线Fig.1 Reflectance spectra of iron-powder covered sea ice

2.2 端元提取

表2、表3分别是观测样本和不同厚度海冰的Acos和SAI。根据式(1)计算样本的光谱适量角余弦得到海冰Acos最小，为0.588，样本1至样本9的Acos值位于海冰和铁矿粉之间，随着铁矿粉面积的增大而增大。铁矿粉的光谱曲线形状近似呈现倒立三角形，海冰的呈直线或者顶角极大的正立三角形。铁矿粉的SAI都大于1，海冰的SAI小于铁矿粉。

表2 不同样本的Acos和SAITable 2 Acos and SAI values of samples

表3 不同厚度海冰Acos和SAITable 3 Acos and SAI values of sea ice with different thicknesses

经对比不同厚度海冰的Acos和SAI发现，利用SAI的M+2SD作为阈值，可以将所有观测厚度的海冰与铁矿粉进行区分，而Acos的M+2SD无法将所有海冰都与铁矿粉区分开，因此，在进行端元提取时，采用的是SAI的M+2SD作为阈值。

2.3 光谱特征与面积相关性

图2分别为面积比例与反射率的相关系数，在350～698 nm，反射率与铁矿粉面积的相关系数绝对值都大于0.90，在918～1 400，1 500～1 780和2 250～2 300 nm波段，相关系数都大于0.90。86%以上波段范围的反射率与铁矿粉面积的相关系数在0.90以上，91.75%以上的波段范围相关系数在0.80以上。铁矿粉面积与反射率在可见光近红外波段的下降和在短波红外波段内的上升呈现出强线性相关。随着铁矿粉的面积的增大，样本在可见光波段的吸收逐渐增强，而在短波红外波段的反射逐渐增强，由于铁矿粉吸收系数较大，太阳辐射无法透过铁矿粉颗粒层，冰底海水反射和冰内粒子散射的能量也无法透过铁矿粉颗粒层。铁矿粉在传感器视场中的面积对反射率的变化起主要作用，海冰-铁矿粉混合像元的反射特征可以用线性混合像元来描述[12-13]。

图2 铁矿粉覆盖面积比例与光谱特征相关性(a)：面积比例与波长相关性；(b)：面积比例与反射率相关性Fig.2 Correlation of iron-powder covered area fraction and spectral features(a)：Correlation of area fraction and wavelength； (b)：Correlation of area fraction and reflectance

2.4 面积比例反演

海冰-铁矿粉像元是线性混合像元，在提取出了海冰-铁矿粉像元之后，可通过选取一个或多个波段的组合参量估算铁矿粉面积比例。经过多次试验，最终选取1 610～1 630 nm波段的平均反射率来估算铁矿粉比例。根据线性混合像元理论可以得到如方程组(5)

(5)

aF=aρ1 610-1 630-β

(6)

如果已知铁矿粉(可以根据先验知识得知)厚度hF、孔隙率P、密度DF，就可以估算铁矿粉的体积，计算式为

mF=AashFPDF

(7)

图3为样本的铁矿粉面积比例反演及其误差，面积较小的误差一般大于面积较大的。可以看出，在克服海冰背景反射率的影响下，可比较精确地估算海冰-铁矿粉像元中铁矿粉面积比例，平均反演精度可达94.23%。

图3 样本的铁矿粉面积比例反演及其误差Fig.3 Predicted area fraction and the error

3 结 论

开展了真实海冰条件下散落铁矿粉的反射光谱观测实验，提取分析了不同面积比例覆盖海冰的光谱曲线特征。利用反射光谱的SAI能够将测试所用不同厚度海冰和被铁矿粉污染的海冰进行区分，实现海冰端元提取。海冰表面的铁矿粉截面面积比例与918～1 400，1 500～1 780和2 250～2 300 nm波段的反射率有强相关性，相关系数都大于0.90。其中，1 610～1 630 nm波段的平均反射率对铁矿粉面积比例的估算精度较高，平均反演精度达到94.23%。

本工作仅考虑了海冰表面受到铁矿粉污染的情况，未对铁矿粉进入海冰内部的情况进行考虑。未来工作中，将逐步完善海冰受污染场景的研究，为海冰受铁矿粉污染的监测提供技术参考。