王志一，等: 高分辨率光学和 SAＲ 遥感影像在地震地质灾害调查中的应用 2018 年seismo-geological disasters． The investigation was carried out from different perspectives including road damageanalysis，extraction of secondary disasters ( e． g． landslides) ，and coseismic deformation field calculation． Theresults show that combining the characteristics and advantages of high resolution optical and SAＲ remotesensing images，can play an important role in the comprehensive investigation of seismic and geologicaldisasters．Keywords: high resolution remote sensing image; In SAＲ; earthquakes geological disasters; emergency survey遥感技术因具有覆盖范围广、图像获取方便等特点，能够客观、全面地反映地震后灾区的景观，能为震害调查、损失快速评估提供科学依据［1 － 2］。我国曾对1966 年邢台地震、1975 年海城地震、1976 年唐山地震、1988 年澜沧 － 耿马地震和 1989 年大同地震等大震进行了震后灾区航空摄影，并积累了丰富的震害遥感影像判读经验［3 － 4］。在 2003 年的伽师 6. 8 级地震中，王晓青等［3］根据以往震害影像统计并结合本次地震震害遥感特征，提出了遥感震害分级分类标准和地震烈度划分标准，进而得到了基于震害遥感影像的伽师地震等震线图，然而，光学影像由于受多云多雨天气状况的限制，如果震后出现云雨天气，则光学影像在地震灾害调查就会受到限制。雷达影像具有全天候、全天时的特点，就可以克服上述不足。在 2008 年的汶川8. 0 级大地震中，刘云华等［2］利用多时相 SAＲ 图像对映秀及周边地区做了地震灾害识别研究，利用幅度图像做比值变化检测以及利用相位信息做干涉处理可以在地震灾害评估中取得更好效果。2017 年 8 月 8 日 21 时 19 分 46 秒在四川阿坝州九寨沟县( 103. 82E，33. 2N) 发生 7. 0 级地震，震源深度约为 20km，震中距九寨沟县 39 km( 图 1) 。紧接而来的降雨又造成灾区发生大面积的崩塌、滑坡等次生地质灾害。在灾区余震与次生地质灾害不断发生的情况下，遥感成为获得震区灾情的主要手段。地震发生后第一时间，中国地质环境监测院根据预报震中位置坐标，启动应急响应机制，迅速搜集查询到震区卫星数据存档影像资料，编制地震灾区震前遥感影像图( 时相: 2013 年 8 月 3 日，资源三号卫星，分辨率 2. 0 m) 。同时，与资源卫星应用中心、二十一世纪空间技术、环球星云科技、航天世景等数据公司收集地震灾区震后卫星数据，于 8 月 10 日晚成功获取了灾后北京二号、高分一号、高分二号等国产卫星数据，及时开展了数据处理及地质灾害解译工作。但受到震后灾区多云、阵雨等天气影响，拍摄区多被云层覆盖，对地质灾害解译有一定的影响。因此，在综合利用不同高分辨率光学卫星数据的同时，对高分辨率 SAＲ 遥感影像也进行了地质灾害调查应用。1 地震地质背景九寨沟位于四川省阿坝藏族羌族自治州九寨沟县漳扎镇境内，地处青藏高原向四川盆地过渡地带，距离成都市 400 多千米，是一条纵深 50 余千米的山沟谷地，总面积 64 297 hm2，森林覆盖率超过 80% 。隶属于四川省阿坝藏族羌族自治州，位于青藏高原东部边缘，阿坝州东北部。东、北与甘肃省文县、舟曲县、迭部县交界，西、南与四川省若尔盖县、平武县、松潘县接壤。九寨沟县地势西北、西南高，东南低，属高原湿润气候。该区域地质背景复杂，碳酸盐岩分布广泛，褶皱断裂发育，新构造运动强烈，地壳抬升幅度大，多种营力交错复杂［5］( 图 1) 。此次地震震中位于岷江断裂、塔藏断裂和虎牙断裂附近( 图 2)［6］。根据中国地震局消息，九寨沟县M7. 0 级地震震源机解结果为: 矩震级 Mw 约 6. 5。断层节面 1 走向 = 326° /倾角 = 62° /滑动角 = － 15°; 节面 2 走向 = 64° /倾角 77° /滑动角 = － 151°。波形拟合的质心深度为 11km。发震构造推测为塔藏断裂南侧分支和虎牙断裂北段，初步推断此地震为一次走滑型为主事件。2 地震地质灾害高分辨率光学遥感解译地震发生后，从北京二十一世纪空间技术股份有限公司获取到北京二号、高分一号、高分二号等国产高分卫星遥感数据，经过辐射定标、正射校正、图像融合、图像裁剪等预处理，为后续地震地质灾害遥感解译提供基础数据。根据震后高分卫星遥感解译，此次地震造成漳扎镇周边及景区多处山体滑坡。从图 3 和图 4( a 为震前2017 年 8 月 3 日影像，b 为震后 2017 年 8 月 9 日影像) 九寨沟黄龙机场至渣扎镇 G544 公路震前震后遥感影像对比可以发现，公路沿线山体发生多处滑坡，震前图像上的道路明显，具有连续性; 震后道路由于受到山体滑坡的影响，造成道路中断。 中国地质灾害与防治学报 ·83·图 1 九寨沟 M7. 0 级地震区地质构造图Fig. 1 The geological structure diagram of Jiuzhaigou M7. 0 earthquake area图 2 主震震中位置及周边历史地震分布图Fig. 2 The location of the main earthquake and the map of the surrounding history王志一，等: 高分辨率光学和 SAＲ 遥感影像在地震地质灾害调查中的应用 2018 年图 3 九寨沟 G544 公路上四寨村附近滑坡震前震后遥感影像对比图Fig. 3 The comparison of remote sensing imagesfor landslides near Sizhai before and after earthquake3 地震地质灾害高分辨率雷达遥感解译3. 1 In SAＲ 地震形变场观测3. 1. 1 雷达干涉测量( In SAＲ) 基本原理合成孔径雷达干涉测量( In SAＲ) 是利用对同一地区观测的两幅 SAＲ 复影像数据进行相干处理，通过相位信息获取地表高程信息及形变信息的技术。合成孔径差分干涉测量( D-In SAＲ) 技术是 In SAＲ 技术的拓展，利用同一地区的两幅干涉图像，其中一幅是形变前的干涉图像，另一幅是形变后获取的干涉图像，然后通过差分处理来获取地表形变的测量技术。D-In SAＲ 几何示意图如图 5所示，在重复轨道条件下，S1 为卫星第一次过境点，卫星轨道高度为 H，S2 为卫星第二次过境点，两次过境之间的距离 B 成为基线，其与水平方向的夹角为 α。一般情况下，垂直基线即 B⊥= Bcosα 应小于一个临界值才能获取质量较好的干涉结果。卫星两次过境与地面目标点的距离分别为 Ｒ1和 Ｒ2，如果在两次过境时间段内，对应的目标点发生了位移 Δd 时，会在卫星视线方向产生相应的位移 Δr，进而反应为电磁波的相位差，通过对相位差的计图 4 九寨沟 G544 公路如意坝附近滑坡震前震后遥感影像对比图Fig. 4 The comparison of remote sensing images forlandslides near Ｒuyiba before and after earthquake算，可以获取两景 SAＲ 影像期间地表发生的形变。图 5 D-In SAＲ 几何示意图Fig. 5 The geometrics ketch of the D-In SAＲ两幅影像进行干涉处理，得到干涉条纹图。干涉图相位 m包含着多种因素的相位［7］，即平地相位 e、地形相位 t、地表形变相位 d及噪声相位 n，公式表示为:m= e+ t+ d+ n高分辨率光学和 SAＲ 遥感影像在地震地质灾害调查中的应用———以九寨沟 M7. 0 级地震为例王志一1，徐素宁1，王 娜2，马秀强3( 1. 中国地质环境监测院，北京 100081; 2. 北京环球星云遥感科技有限公司，北京 100101;3. 中国地质大学( 北京) ，北京 100083)摘要: 快速全面的获取地震灾情信息对于减少灾害损失、预防次生灾害发生具有重要的意义。遥感技术凭借其快速、大范围、高精度地获取地面信息的特点在震后信息提取中发挥了重要作用。传统光学遥感数据由于受到震后多云雨天气的影响，往往难以获取有效数据。合成孔径雷达( SAＲ) 能够全天候获取数据，已成为震害评估的重要数据源。随着遥感图像向着高分辨率方向发展，数据所包含的信息也越来越丰富，更有利于地震地质灾害调查。本文以 2017 年九寨沟地震为例，采用高分辨率光学和 SAＲ 遥感图像为数据源，从道路震害分析、滑坡次生灾害提取、同震形变场计算等多个角度开展应用研究，对此次地震地质灾害进行调查。结果表明，结合高分辨率光学和 SAＲ 遥感影像各自的特点和优势，能够在地震地质灾害信息全面调查中发挥重要作用。关键词: 高分辨率遥感; In SAＲ; 地震地质灾害; 应急调查中图分类号: P642; P237 文献标识码: A 文章编号: 1003-8035( 2018) 05-0081-08收稿日期: 2018-01-03; 修订日期: 2018-03-19基金项目: 国家重点研发计划课题: 城市群地质环境演化多源遥感监测与预警( 2017YFB0503803)第一作者: 王志一( 1984-) ，男，高级工程师，主要从事资源与环境遥感研究工作。E-mail: wangzy@ mail． cigem． gov． cnApplication of the high resolution optical and SAＲ remotesensing data images induced by the Jiuzhaigou M7. 0earthquake geological hazards surveyWANG Zhiyi1，XU Suning1，WANG Na2，MA Xiuqiang3( 1． China Institute for Geo-Environmental Monitoring，Beijing 100081，China; 2． Beijing Earth Star Ｒemote SensingTechnology Co． LTD． ，Beijing 100101，China; 3． China University of Geosciences ( Beijing) ，Beijing 100083，China)Abstract: Ｒapid and comprehensive acquisition of earthquake geological disaster information is of greatsignificance to reduce the loss of disaster and prevent secondary disasters． Ｒemote sensing technology hasplayed an important role in the extraction of earthquake information due to its ability of fast，wide range andhigh precision acquisition of ground information． The traditional optical remote sensing data are easily affectedby the cloudy and rainy weather after earthquake，so it is normally difficult for traditional optical remotesensing to obtain effective data． Synthetic aperture radar ( SAＲ) ，which is capable of obtaining data in all-weather condition，has already become an important data source for earthquake damage assessment． As theresolution of remote sensing images are becoming higher，the information contained in the data are alsoincreasingly richer，which is more favorable to post-earthquake investigation． Jiuzhaigou earthquake is taken asan example，high-resolution optical and SAＲ remote sensing images are used as data source to investigate the 中国地质灾害与防治学报 ·85·干涉图中的噪声相位，可以通过多视、滤波去噪等步骤处理。干涉图中的平地相位，也称参考趋势面相位，在干涉图中占主导地位。利用轨道基线求出影像重叠区的同名像点的多个参考斜距差，选择 N 介多项式模型描述可以去除平地相位［8］，生成新的干涉图。干涉图中的地形相位，是由于地形起伏引起的相位差，为了获取地表的形变信息，必须将地形相位去除。通常采用外部的 DEM 数据结合轨道参数去除地形相位，即“二轨法”，最终得到同震形变干涉图［9］。根据 D-In SAＲ 测量的基本原理，“二轨法”的主要步骤包括: 主辅影像预处理、影像配准、重采样、干涉图生成、滤波、形变相位分离、相位解缠、地理编码等，具体流程如图 6 所示。图 6 D-In SAＲ“二轨法”处理流程Fig. 6 The process flow of D-In SAＲ “two-track”3. 1. 2 九寨沟地震 In SAＲ 形变场观测结果地震发生后的 8 月 10 日晚在环球星云德清遥感卫星地面接收站成功获取了灾后第一组同震重轨干涉SAＲ 影像对，利用 ＲADAＲSAT-2 震前( 2017 年 5 月 30日) 、震后( 2017 年 8 月 10 日) 的数据，外部 DEM 采用SＲTM( 30 m) 数据，基于 D-In SAＲ“二轨法”处理方法，最终得到此次地震的同震形变场，结果如图 7 左图所示。同震形变场的局部细节如图 7 右图所示，可以看出断层东侧失相干比较严重，西侧约有 4 个干涉条纹，最外侧条纹距离中国地震台网中心确定的震中心约为36. 64km，从而可以大致推断出此次地震形变场的范围。对 In SAＲ 获取的同震形变场进一步反演，可以为理解此次地震的发震机制以及未来地震灾害演化趋势提供重要的依据。3. 2 雷达影像地震地质灾害解译3. 2. 1 雷达影像地震地质灾害解译基本原理相对于光学遥感图像受天气和光照条件的限制来说，雷达遥感影像具有全天候、全天时、不受云层覆盖影响等优点。地震发生后，需要对大范围的灾情信息进行宏观掌握，圈定极震区，以及对各类震害快速定位，如地质灾害、次生灾害、建筑物震害等信息，由于SAＲ 图像是侧视成像，当地形起伏有变化时，就会产生透视收缩、叠掩和雷达阴影等像点位移，因此利用SAＲ 图像进行震害信息提取的首要工作是影像图像几何校正与配准。利用 SAＲ 图像进行地质灾害识别主要是利用SAＲ 图像的幅度信息，对比分析、解译灾害发生前后的影像，从而提取灾害发生的位置、分布、面积大小等信息。其关键技术包括: SAＲ 图像滤波辐射校正等预处理、灾害前后影像精配准、构造差异图像、地理编码和变化区域提取等。对 SAＲ 图像进行几何纠正和配准主要是为 SAＲ 图像赋予准确的空间坐标信息，消除由 SAＲ 斜距成像和地表地型起伏等因素引起的几何变形，将 SAＲ 图像重投影至指定的坐标系统中，便于 SAＲ影像中灾情信息的定位以及与其它数据进行分析对图 7 In SAＲ 提取九寨沟地震形变场及局部图Fig. 7 Jiuzhaigou seismic deformation field and local enlarged diagram extracted by In SAＲ王志一，等: 高分辨率光学和 SAＲ 遥感影像在地震地质灾害调查中的应用 2018 年比［10］。基于雷达影像地震地质灾害解译流程如图8 所示。图 8 基于雷达影像地质灾害解译流程图Fig. 8 Flow chart of geological disasterinterpretation based on radar image图 9 九寨沟地震滑坡 SAＲ 监测图Fig. 9 SAＲ monitoring diagram for landslidescoused by the Jiuzhaigou earthquake3. 2. 2 九寨沟地震雷达影像地质灾害解译结果此次地震造成了多处山体滑坡，造成了震区道路多处损毁，影响地震应急救援的效率。图 9 显示震区一处滑坡造成了道路损毁的图像变化，从地震前( 图 9a) 和地震后( 图 9b) 的图像对比可以很明显的看出，震前图像上的道路明显，具有连续性; 震后道路由于受到山体滑坡的影响，造成道路中断，暗色调原有的道路被碎石覆盖，散射特性发生了变化。地震发生后，九寨沟景区内部分景点损毁，其中的火花海、诺日朗瀑布等景点损坏严重，几乎“消失”。火花海水位整体下降，已经见底，湖边发生小面积滑坡，也裸露出了黄土。如图 10 所示，从地震前( 图 10a) 和地震后( 图 10b) 的 ＲADAＲSAT-2 图像对比可以看出，火花海湖面原有的暗色调水体在震后消失，湖水消失见底。图 10 九寨沟地震湖泊损毁 SAＲ 监测图Fig. 10 SAＲ monitoring diagram for lakedamaged by the Jiuzhaigou earthquake从地震前( 图 11a) 和地震后( 图 11b) 的光学卫星影像对比也可以明显看出，红框内火花海湖面原有水体在震后消失，湖水消失见底。4 讨论与结论近年来，国产高分辨率遥感数据随着数据源的不断 中国地质灾害与防治学报 ·87·图 11 九寨沟景区火花海震前震后光学遥感影像对比图Fig． 11 The comparison of remote sensing imagesfor landslides before and after the Jiuzhaigou earthquake丰富，数据质量的不断提升，在国民经济社会发展的各行各业均得到了大力的应用，尤其在一些特大自然灾害应急工作中，更是发挥了极其重要的作用，为抗震救灾及灾害评估提供了可靠的数据支撑。同时，在地震发生往往受到阴雨天气的影响，造成光学遥感影像在救灾决策中的作用受到了很大的限制的情况下，SAＲ 图像由于其全天候的特点也成为了地震发生后重要的遥感数据源。利用光学和 SAＲ 遥感数据，特别是高分辨率遥感影像，对全方面掌握和理解地震灾情具有十分重要的意义。本次研究采用高分辨率光学遥感数据与 SAＲ 数据结合的分析解译技术，利用多时相光学和雷达数据对受地震影响较大的渣扎镇周边地区做了地震灾害识别，同时在利用雷达干涉处理识别地震破裂带方面取得了较好的效果。研究表明，利用震前震后的光学影像和 SAＲ 幅度图像进行变化监测，能够在山体滑坡、堰塞湖堵塞等震害信息提取取发挥重要的作用。通过结合光学影像和SAＲ 影像各自的特点和优势，，快速进行地震地质灾害的应急调查，将成为未来地面应急调查以及灾害评估的主要数据和技术支撑手段。参考文献:［1］ 王志一，甘甫平，李贤庆，等． CBEＲS － 02B 星数据在“5·12”地震灾害监测中的应用———以唐家山堰塞湖动态监测为例［J］． 国土资源遥感，2009，21( 1) : 89 － 92．WANG Zhiyi，GAN Fuping，LI Xianqing，et al． Theapplication of CBEＲS-02B data to the disastermonitoring of the “5. 12 ”Wenchuan Earthquake: Thedynamic monitoring of the Tangjiashan barrier lake［J］．Ｒemote Sensing for Land ＆ Ｒesources，2009，21 ( 1 ) :89 － 92.［2］ 刘云华，屈春燕，单新建，等． SAＲ 遥感图像在汶川地震灾害识别中的应用［J］． 地震学报，2010，32( 2) :214 － 223.LIU Yunhua， QU Chunyan， SHAN Xinjian， et al．Application of SAＲ data to damage identification of theWenchuan earthquake［J］． Acta Seismologica Sinica，2010，32( 2) : 214 － 223.［3］ 王晓青，魏成阶，苗崇刚，等． 震害遥感快速提取研究: 以 2003 年 2 月 24 日巴楚 － 伽师 6. 8 级地震为例［J］． 地学前沿，2003，10 ( S1) : 285 － 291.WANG Xiaoqin，WEI Chengjie，MIAO Conggang，etal， The Extraction of seismic damage from remotesensing images—A case study of Bachu-Jiashiearthquake with Ms = 6. 8 occurred on Feb． 24［J］．Earth Science Frontiers，2003，10( S1) : 285 － 291.［4］ 魏成阶，先博勤，张宗科，等． 地震灾害航空遥感快速调查技术研究———以唐山地震区作为模拟试验场［M］/ /何建邦，田国良，王劲峰． 重大自然灾害遥感监测与评估研究进展． 北京: 科学出版社，1993.WEI Chengjie，XIAN Boqin，ZHANG Zongke，et al．Study on fast investigation of earthquake disasters fromairborne remote sensing ［M］/ / HE Jianbang， TIANGuoliang，WANG Jinfeng． The Progress of monitoringand evaluation of natural disasters by remote sensing．Beijing: Science Press，1993．［5］ 中国地震局． 区域构造图［EB/OL］． ［2017 － 08 －10］ http: / / www． cea． gov． cn / publish / dizhenj /468 /553 /101710 /101720 /201708082344275 23689626 /index． html．China Earthquake Administration． Ｒegion geologicalstructure map［EB / OL］． ［2017 － 08 － 10］． http: / /www． cea． gov． cn / publish / dizhenj /4 68 /553 /101710 /·88· 王志一，等: 高分辨率光学和 SAＲ 遥感影像在地震地质灾害调查中的应用 2018 年101720 /201708082344275 23689626 / index． html.［6］ 中国地震局． 主震震中位置及基本参数［EB/OL］．［2017 － 08 － 08］．China Earthquake Administration． Location and basicparameters of main earthquake epicenter ［EB / OL］．［2017 － 08 － 08］．［7］ 张景发，龚丽霞，姜文亮． PS In SAＲ 技术在地壳长期缓慢形变监测中的应用［J］． 国际地震动态，2006( 6) : 1 － 6.ZHANG Jingfa， GONG Lixia， JIANG Wenliang．Application of PS In SAＲ technique to measurement oflong-term crustal deformation［J］． Ｒecent Developmentin World Seismology，2006( 6) : 1 － 6.［8］ 刘国祥． In SAＲ 数据处理及关键算法［J］． 四川测绘，2005，28( 1) : 45 － 47.LIU Guoxiang． Data reduction and key algorithms inIn SAＲ［J］． Surveying and Mapping of Sichuan，2005，28( 1) : 45 － 47.［9］ 王微，薛腾飞，张景发，等． 基于 In SAＲ 同震形变场的烈度评估研究［C］/ /地壳构造与地壳应力文集．2017: 117 － 130.WANG Wei， XUE Tengfei， ZHANG Jingfa， et al．Ｒesearch on seismic intensity evaluation based onIn SAＲ coseismic deformation field［C］/ / Bulletin of theInstitute of Crustal Dynamics． 2017: 117 － 130.［10］ 袁文龙． 星载 SAＲ 图像几何校正原理及软件设计与实现［D］． 北京: 中国科学院研究生院( 电子学研究所) ，2005.YUAN Wenlong． Theory and software realization ofsatellite SAＲ image geo-rectification ［D］． Beijing:Graduate University of Chinese Academy of Sciences( Institute of Electrics) ，2005.( 上接第 80 页)［4］ 王鹏． 地基 SAＲ 干涉测量原理及其形变监测应用研究［J］． 测绘信息与工程，2012，37( 4) : 22 － 25 + 28.WANG Peng． Ground-based SAＲ interferometryprinciples and its applications to displacementmonitoring［J］． Journal of Geomatics，2012，37 ( 4 ) : 22－ 25 + 28.［5］ 吕刚． 贵州重大地质灾害及影响因素分析［J］． 贵州地质，2016，33( 2) : 108 － 112.LYU Gang． Analyses of important geo-disastersdistribution rules and influence factors of Guizhou［J］．Guizhou Geology，2016，33( 2) : 108 － 112.［6］ 李强． PSIn SAＲ 技术在地表沉降监测中的应用研究［D］． 太原: 太原理工大学，2012.LI Qiang． Application of PSIn SAＲ in surface subsidencemonitoring ［D ］． Taiyuan: Taiyuan University ofTechnology，2012.［7］ 陈怡曲． 基于 In SAＲ 的形变监测技术研究［D］． 成都:电子科技大学，2013.CHEN Yiqu． The study of deformation monitoring basedon SAＲ interferometry ［D］． Chengdu: University ofElectronic Science and Technology of China，2013.