高光谱遥感在断层气监测中的应用

空间对地观测研究是近年来国内外地震科学研究发展的重要方向。由于新的观测技术和理论的快速发展，使得高光谱在断层逸出气方面的应用在强震前兆信息分析中有可能发挥重要作用。本文简单介绍了断层气体地面监测和卫星遥感监测方面的应用现状，概述了高光谱遥感的原理、反演方案和方法等。笔者建议运用地面监测和高光谱遥感技术相结合的方法来监测强震前后活动断层逸出气的光谱变化，根据数学算法反演痕量气体组分浓度，最后，结合震例，探讨强震前后断层逸出气体异常变化量与强震的关系，提取断层气体各组分强震前兆信息。利用空间对地观测技术来研究断层逸出气，这样可以全面提高对强震前兆异常信息的空间分析和可靠性判定能力，对全面发展地震监测预报技术具有重要的意义。

引 言 

断层气是地球内部气体及挥发组分不断沿着活动板块和块体边界以及其他活动性深断裂等地壳薄弱地带向地表迁移和释放的气体，活断层是连接深部和深部气体释放的天然通道。断层气的逸出空间主要集中在断裂带附近，大量资料表明断层活动越强的地区，断层气含量越高。监测断层气释放的浓度变化，已成为活断层探测和地震前兆分析的重要方法。

我国地震监测预报的目的是通过对地壳各种物理量的监测，获取强震前兆异常信息，探索地震预报方法，减轻地震灾害。由于国力和自然条件的限制以及观测环境的干扰和破坏，目前我国地面地震台网密度仍然很低，仅能获取地表离散点的数据，很难对地球物理和地球化学背景场进行完整的描述，要完整捕获地震前兆信息比较困难，急需引进新的观测技术。

高光谱能获取地物目标的精细光谱特征，综合地面目标的空间、时间、光谱等特征，探测各种目标的成分属性及状态属性。其波段多，光谱通道数多且是连续的，光谱分辨率高，信息量大，图谱合一，有利于利用光谱特征研究地物和采用各种光谱匹配模型对地物进行精细分类与识别。因此应用高光谱卫星遥感技术，可以获取更多更准确的断层气体信息。

本文介绍了地面和卫星遥感监测断层气体的研究现状，重点概述了高光谱遥感的基本原理、反演方案和方法等，并归纳总结了基于地面和高光谱遥感卫星监测相结合的断层逸出气强震预测指标异常识别方法和技术途径。

1断层气体研究现状

1.1我国近 40年的近地表断层气的观测研究进展

我国在气体地球化学动态地震异常及其预报探索方面积累了大量的地震断层气前兆异常实例资料，尤其是氡（Rn）、汞（Hg）、Ｈ2、He、CO2等地震敏感气体组分方面，并且在全国强震活动地区各地震台站配置有相应的气体观测设备，为我们今后的工作奠定了实验数据基础。例如：1986年两次门源地震和 1987年张掖地震前长草沟测点的Rn、Hg、Ｈ2、CO2四种气体发生同步异常（孔令昌,1996）；1990年青海茫崖 6.7级地震、共和6.9级地震前断层气Rn、Hg的异常变化特点是临震突变幅度大，异常反应范围广（何跟巧,1990）。1992年 1月 12日在甘肃嘉峪关南5.4级地震，嘉峪关观测点距震中约20ｋｍ,地震前,断层气Rn、Hg、Ｈ2、He和 Ｎe五种微量气体组份出现异常（何跟巧,1994）。近40年的近地表断层气体的观测结果为我们进行空地对比研究奠定了地面监测的数据基础。 

1.2卫星遥感在地震气体监测中的应用现状 

20世纪 90年代以来,由于遥感的理论和技术的发展、新型传感器的出现和进步,地震学者通过卫星遥感数据反演获得地震前后震中及周边地区大气中的某些气体含量的变化,对地震前后排气现象进行了大面积系统的研究。

姚清林等（2005） 对2000年4月30日青藏高原CO大面积逸出的情况进行探究。结论表明CO含量异常升高的区域具形状不规则的圈层结构；整个CO逸出区CO体积分数值为2002年正常值的1.57～4.10倍,且该现象在之前数天内是持续存在的,并认为是2000年6月6日甘肃景泰Ｍ5.9地震和2000年6月8日缅甸北部Ｍ6.9地震的前兆。

郭广猛等（2006）同样研究了2002年3月31日台湾7.5级强震前的CO异常。研究发现2002年3月3日、27日、30日,CO在冲绳海槽和台湾北部出现异常区,其面积随着时间扩大,且与MODIS热红外高温异常区位置吻合,并认为CO、温度异常是地下岩石受力破裂、地下气体逸出所致。

Pulinets等（2007）通过对1985年墨西哥 Miehoaean7.8级地震前后的热红外数据进行分析,发现一些地区的气温和水汽含量出现异常,特别是位于Cocos和Rivera板块边界上的Baja California地区。Pulinets认为水汽的产生和地震排放的Rn气体有关,并证明了Rn排放不仅限于地震震中,在附近板块作用地区也存在。

通过上述观测事实综合分析,可以发现地球是一个统一的系统,地震导致地球岩石圈破裂,使得岩石圈以至地幔的断层气体排放到空气中,由于岩石圈－大气圈耦合,导致地震破裂区上空气体的含量发生变化。随着卫星观测技术和反演技术的完善,地震前后断层气体浓度变化能得到更好的分析研究。

2高光谱遥感在断层气方面的应用研究

高光谱遥感具有非常高的光谱分辨率,它不仅可以探测到比常规遥感更精细的地物信息,而且也能探测到更精细的大气吸收特征（束炯,2006）。大气中的分子和粒子成分在太阳反射光谱中有强烈反应,这些成分包括水汽、二氧化碳、氧气、臭氧、云和气溶胶等。传统宽波段遥感方法无法识别出由于大气成分的变化而引起的光谱差异,而波段很窄的高光谱则能够识别出这种光谱差异,成像光谱仪记录的太阳反射光谱中可反映出大气层中分子和粒子成分的光谱信息（袁迎辉,2007）。

目前,高光谱遥感研究大气的主要目标是水蒸气、云和气溶胶。同样,我们可以推广应用到断层气体监测研究中,利用高光谱遥感的高空间分辨率、高光谱分辨率数据,通过光谱识别技术,将各种气体加以区分,再根据反演方法,计算其浓度。

2.1高光谱遥感在断层气体研究方面的反演方法

光谱是表征物质物理属性的主要手段,利用光谱特征可以实现断层气体的遥感监测。电磁辐射与原子和分子间的相互作用是高光谱遥感探测气体成分及特性的基础,根据断层气体成分在紫外、可见和红外光谱波段的特征吸收性质来反演其浓度。

如图 1所示,高光谱遥感将入射光源通过成像光谱仪变为干涉光通过气体样品,经过光电转换还原出气体的光谱图,光谱图含有丰富的谱元信息,具有高的信噪比和光谱分辨率（T∙Ogawa,1994；Gerald Wetzel,1995；Reinhard Beer,2001）。 

利用高光谱遥感所得到的气体光谱,我们可以根据断层气体成分在紫外、可见和红外光谱波段的特征吸收性质来反演其浓度。

1.实施断层气体反演需要以下几个基本条件

（1） 基于逐线计算辐射传输模型的快速参数化模型，模拟计算气体的辐射传输过程。一个精确的前向模型是成功反演的保障。它必须全面考虑辐射传输的物理过程，在给定大气廓线、光谱数据库、辅助数据情况下，准确快速计算出传感器入瞳处的辐射亮度。Tjemakes等人（2003）做了逐线计算程序的比对试验,共有13个研究小组的7个不同的逐线计算辐射传输模型参加了比对实验,结果表明在大多数红外光谱区这些模型都能在观测噪声水平内复制真实观测值。目前还没有在整个红外波段都表现优异的逐线计算的辐射传输计算模型,大气辐射测量计划的逐线计算辐射传输模式是痕量断层气体反演比较常用的前向模型。通常用于痕量断层气体探测的传感器都会根据自身的科学目的以及达到此目的所使用的波段来选择在所选波段性能优异的前向模型或者在已有前向模型的基础上开发特 定的前向模型。

（2）先验知识：断层气体反演中使用的先验知识主要来源于三个方面：已有数据库的气体廓线数据,如AFGR大气廓线数据库、TIGR初始猜测值数据库等；地面观测数据,如美国国家海洋大气管理局（NOAA）的气候监测与诊断实验室独立观测数据、平流层变化监测网络（NDSC）分布于世界各地的地面观测站的地面观测数据等；一些大气化学模型产生数据,如三维全球化学与传输模型（GEOS-CHEM）、对流层三维化学传输模型（IMAHES）等（CliveD∙Rodgers 1976,2000）。 

最小化代价函数方案。通常,我们可以用下式来描述观测向量和真实状态向量之间的关系（CliveD∙Rodgers 2000）:ｙ＝Ｆ（ｘ，ｂ） ＋ε

公式中，ｙ为观测向量（测量的红外辐亮度）；Ｆ（ｘ，ｂ）为前向模型函数；ｘ为状态向量（待反演的大气参数）；ｂ为前向模型所有其他的参数：压强、温度、其他组分的廓线、地表发射率、仪器的光谱响应函数和光谱数据库；ε为由仪器噪声引起的测量误差。反演就是要从给定ｙ的条件下反求出ｘ。

2.1.1反演方法

目前，遥感地震资料反演的方法很多，大致有如下几种：地震资料粒子群反演方法、最小二乘法、叠前地震反演方法、地震非线性随机反演方法、多尺度地震资料联合反演方法、非线性迭代方法等，下面主要介绍非线性迭代方法。

对于断层气体遥感反演来说，问题是非线性的（Cathy Clerbaux,2002）。目前解决这类问题主要是在利用先验知识的同时,利用非线性迭代方法求得统计意义上的最优解。

对于气体探测领域的研究Rodgers做了大量开创性的工作,尤其在利用优化估计方法反演气体参数方面（CliveD∙Rodgers 1976,2000）。

定义代价函数 φ（ｘ）： 

公式中ｇ（ｘ）为罚函数；λ为正则化因子。当 λ≠ 1时，可通过最小化代价函数得到正则化解。当 λ＝1，可通过最小化代价函数得到最小方差解，将其写成矩阵表达的形式：

由公式可以得到牛顿迭代法的一般形式。忽略二阶导数项，引入先验知识，根据构造罚函数的不同，得到Gauss-Newton迭代法：

和Levenber-Marquardt迭代法

公式中，Sε为观测误差的协方差矩阵；Sａ为先验状态向量 xａ的协方差矩阵；

Ｋｉｊ＝∂Fｉ／∂xi为 权函数矩阵（S∙A∙Clough，1995）。

Gauss-Newton法和Levenberg-Marquardt法都是牛顿方法的变形,它们不需要计算海色（Hessian）矩阵,具有较高的计算效率,虽然它们只具有线性收敛性,但是通过选择合理的初始值,也能达到较快的收敛速度,这对于气体反演问题来说非常有效。我们在处理断层气体时，可以选择合理的地面观测初始值，在运用非线性迭代方法时，使断层气体光谱线性收敛，这样就可以快速进行反演计算。

2.1.2误差分析

由于断层气体遥感反演是非线性的，没有唯一解，所以我们不仅需要从有限个可能的解中选择一个某种意义上的“最优 ”解，而且还需要对这个“最优 ”解进行精度评定，分析各种误差源的传播以及它们对反演精度的影响。

在利用非线性迭代方法线性近似的情况下，真实状态向量和反演的状态向量的差别可以写成以下形式：

其中：

（1）

为平核矩阵。它的迹等于观测所含的独立信息量个数。

（2）

Ｇε称为测量误差，由仪器噪声引起，其协方差为：

（3）

称为模型参数误差，由对模型参数知识的不完全引起。其协方差为：

其中Ｓb为代表模型参数不确定型的协方差矩阵。

总的误差协方差矩阵为：

利用高光谱遥感所得到的气体光谱，我们可以根据断层气体成分在紫外、可见和红外光谱波段的特征吸收性质来反演其浓度。实施断层气体反演时，要先计算入射光在进入传感器时的辐射亮度，即前向模型条件，并要有已知的先验知识数据库，即各气体组分应对应怎样的光谱曲线，最后在利用非线性迭代方法线性近似的情况下，将反演获得的数值减去总的误差，就得到了断层气体各组分的光谱测量结果，提取断层气体各组分的最优光谱数据。

2.2高光谱活动断层逸出气强震预测指标信息提取研究的方法和技术途径 

目前,高光谱遥感在火山气体监测方面取得了一定成果。高光谱图像在红外范围（1900－2100ｎｍ）目前已经投入使用,其中包括二氧化碳吸收曲线。这些航空机载的可见光／红外成像光谱仪图像在夏威夷Kilauea East Rift地区的PuuOo Vent上空获得。利用辐射传输模型和新开发的映射技术模拟测量上升流光谱辐射,获得PuuOo Vent plume二氧化碳浓度图（Claudia Spinetti 2008）。鉴于上述高光谱在火山气体方面的应用，同样，我们可以开展高光谱遥感在断层气方面的强震监测预报。

随着遥感数据处理软件和数据反演算法业已成熟，我们可以采用这些软件和方法来处理高光谱卫星数据，提取断层逸出气信息。同时，进行野外实地断层逸出气光谱测量，对比高光谱遥感和近地表光谱结果，探讨高光谱断层逸出气强震预测指标异常识别方法。笔者归纳总结出具体的技术途径如下：

高光谱卫星遥感数据处理。运用卫星地球化学遥感数据处理软件 （ENVI、PCI、ERDAS） 等，同时借鉴环境科学方面的软件，对卫星数据进行处理，反演计算断层气体各组分的浓度，探究提取H2O、CO、CO2、He、Ｈ2和CH4强震断层气体信息的方法。

收集地震地质资料、各种气体组分的反射特征光谱资料。在地震活动强烈的地区进行野外断层气体测量，确定不同地震地质区域的特征气体组分；实验室标定地震地球化学特征组分的特征波谱，诊断特征组分蚀变和地表异常的特征波谱，在地震断裂带进行特征气体组分通量测定，建立野外光谱数据库。

对比高光谱遥感和近地表光谱结果，结合实验室标准光谱，探讨高光谱断层逸出气强震预测指标异常识别方法。

结合震例，把地震活动与利用高光谱遥感数据获得断层气体地球化学场的变化联系起来，进行构造活动和地震异常判定，归纳总结出高光谱活动断层逸出气强震预测指标，并进行震例检验。

高光谱活动断层逸出气强震预测指标信息提取技术途径如图 2所示。

3结 语

高光谱遥感作为一门新兴的传统遥感发展的前沿技术，具有光谱分辨率高、信息量大、图谱合一的特点‚这是传统的遥感技术无法比拟的。高光谱遥感的这些特点，使得对地物的分类识别，物化信息的提取和预测的广度、深度和精度大大提高，遥感技术从定性为主向更高精度的定量遥感发展。同时其所拥有的巨大的信息量，为其向更广领域的拓展提供了巨大的潜力。可以相信，随着高光谱遥感技术的发展，结合地面观测，我们可以更好地确定与强震断层逸出气体组分有关的变化；结合震例，探讨强震前后断层逸出气体异常变化量与强震的关系；结合地震孕育区其他有关的地震活动信息和前兆异常现象，预判未来地震可能发生的地点，提高地震预报水平。

现在我国地震系统中，四川省地震局已建立接收高光谱遥感数据的平台，当务之急应充分利用现有资源，把高光谱遥感数据与地面台站监测数据结合起来，研究提取断层气异常信息的方法，探索强震前兆信息的判定指标。

利用空间对地观测技术，把高光谱遥感监测和地面监测综合起来，并大力加强地震前兆与震源关系的深入研究。对不同地区、不同类型的地震，其前兆种类、幅度、分布范围、持续时间存在着一定的差异。通过地震前兆场孕育模式及地震前兆场时空分布与孕震动力环境关系的研究，取得地震前兆时空分布与演化规律的正确认识，继续对地震遥感异常的成因机制进行研究，将会使高光谱遥感技术在地震预报中具有更广阔的应用前景。

莱森光学- iSpecHyper-VM系列无人机机载高光谱成像系统

莱森光学-多旋翼无人机高光谱成像系统iSpecHyper-VM系列是一款基于小型旋翼无人机的高性能机载高光谱成像系统。

莱森光学-多旋翼无人机高光谱成像系统光谱范围400-1000nm/900-1700nm,由高性能高光谱相机、稳定云台、高清相机、GNSS模块、机载控制与数据采集模块、机载供电模块、地面站模块等部分组成。核心载荷高光谱相机完全自主研发，采用1英寸大靶面CCD图像传感器具有高光谱分辨率、高灵敏度、大视场及优异的成像性能，配合定制开发的高性能稳定云台，能够有效降低飞行过程中无人机抖动引起的图像扭曲与模糊；同时与GPS同步触发，实现可见光照片匹配同步GPS信息；通过地面实时数据平台实时观测飞行器采样点并可利用地面端设置采集的航线预览及矫正功能；有辐射度校正、反射率校正、区域校正等批处理能力。辅助相机可实时可见、监控拍摄效果，有图像实时回传功能；有辐射度校正、反射率校正、区域校正等批处理能力，光谱相机控制、数据采集、自动曝光、自动扫描速度匹配、辅助摄像头功能、支持远程遥控、支持巡航、惯导采集模式，数据支持Envi等主流遥感软件，实时常用植被指数计算功能、光谱及图像数据预览、建库、建立指标解译库等功能。

莱森光学-多旋翼无人机高光谱成像系统与大疆M600 pro/M300RTK无人机完美适配，支持同类型的多种无人机；系统支持配件升级及定制化开发，为教育科研、智慧农业、林业调查、水环境监测、目标识别、军事反伪装等行业高端应用领域提供高性能解决方案。

详情可见：http://www.lisenoptics.cn/duoxuanyiwurenjigaoguangpuchengxiangxitong/833.html

莱森光学- iSpecHyper-VS系列高光谱成像相机/便携式高光谱成像系统

莱森光学-iSpecHyper-VS系列高光谱成像相机是莱森光学（LiSen Optics）专门用于公安刑侦、物证鉴定、精准农业、矿物地质勘探等领域的最新产品，主要优势具有体积小、帧率高、独有高光通量分光设计、信噪比灵敏度高、大靶面探测器、高像质等特点。iSpecHyper-VS采用了透射光栅外置推扫或内置原理高光谱成像，系统集成高性能数据采集与分析处理系统，高速接口传输，全靶面高成像质量光学设计 ，物镜接口为标准C-Mount，可根据用户需求更换物镜。iSpecHyper-VS系列高光谱成像相机，广泛应用于公安刑侦、物证鉴定、精准农林、遥感遥测、 工业检测、 医学医疗、采矿勘探等各领域。

莱森光学-iSpecHyper-VS系列高光谱成像相机/便携式高光谱成像系统主要技术特点

² 光谱范围400-1000nm或900-1700nm，分辨率优于3nm

² 独有高光通量分光成像设计、信噪比灵敏度高

² 高性能CMOS/CCD/InGaAs(TE Cooled)图像传感器，数据格式支持ENVI,支持多区域ROI

² 可选配电控自动对焦技术、自动曝光、自动成像扫描匹配

² 全靶面高成像质量光学设计，点列斑直径小于0.5像元

² 多种焦距物镜镜头（12.5/25mm/35mm/75mm）可根据用户需求更换物镜

² 多种光学附件可选：户外太阳光源、实验室太阳匀化光源工作台、三脚架（野外专用）、数据线、电源线、电源（外置电源模块）、仪器便携式手提箱（定制三防箱）、专业高光谱采集软件

详情可见：http://www.lisenoptics.cn/gaoguangpuchengxiangxiangji/834.html

莱森光学- iSpecField-系列手持式/便携式地物光谱仪

莱森光学-iSpecField系列地物光谱仪是莱森光学（LiSen Optics）专门用于野外遥感测量、土壤环境、矿物地质勘探等领域的最新明星产品，由于其操作灵活、便携方便、光谱测试速度快、光谱数据准确是一款真正意义上便携式地物光谱仪。iSpecField系列地物光谱仪光谱范围250-1100nm/250-1700nm/250-2500nm,采用了工业级触控显示屏手柄探头，手柄探头同时采用了独有光学设计内置摄像头（相机）、GPS、激光指示器、内置光学快门控制，同时地物光谱仪主机与工业级触控显示屏手柄探头一体化设计，可野外现场直接进行地物光谱操作测量，野外操作更加便捷方便，非常适合复杂的野外地物光谱测量。

莱森光学-iSpecField系列地物光谱仪主要技术特点

² 光谱范围250-1100nm/250-1700nm/250-2500nm，固定全息光栅一次性快速扫描分光

² 2048像素面阵BT-CCD，256/512像元InGaAs，高像素双路探测器同步测量，光谱精度高、分辨率高

² 主机与工业级触控显示手柄探头一体化结构，野外测量无需额外电脑，操作灵活

² 最短积分时间30微秒，测量动态范围大，7吋高清触控显示远程触发一键测量

² 内置＞800万像素自动对焦摄像头（相机）、GPS、激光指示器、内置光学快门控制

² GPS定位、地物图像、角度测量、距离测量、空气温湿度测量与光谱实时对应同步

² SpecAnalysis专用地物分析软件，兼容ENVI、TSG、Arcgis等第三方工具软件，嵌入了USGS数据库和NDVI等13个植被指数

² 丰富测量光学配件：标准白板/灰板、叶片透射夹、矿物土壤专用探头、手枪式光纤探头、室内太阳光源、视场角镜头、BRDF测量附件、透反射实验室支架装置、远程触发器等可满足野外和实验室测量需求，可实现透射反射率、辐照度、辐亮度等测试

² 大容量电池，续航时间4-5小时，供电电池模块可拆卸，可配置备用电池模块，满足长时间野外测量

详情可见：http://www.lisenoptics.cn/bianxieshidiwuguangpuyi/#go