欢迎访问一起赢论文辅导网
经济管理论文
当前位置:首页 > 经济管理论文
高分辨率光学和SAR遥感影像在地震地质灾害调查中的应用
来源:一起赢论文网     日期:2019-01-29     浏览数:926     【 字体:

 王志一,等: 高分辨率光学和 SAR 遥感影像在地震地质灾害调查中的应用 2018 seismo-geological disastersThe investigation was carried out from different perspectives including road damageanalysisextraction of secondary disasters ( eglandslides) and coseismic deformation field calculationTheresults show that combining the characteristics and advantages of high resolution optical and SAremotesensing imagescan play an important role in the comprehensive investigation of seismic and geologicaldisastersKeywords: 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]利用多时相 SAR 图像对映秀及周边地区做了地震灾害识别研究,利用幅度图像做比值变化检测以及利用相位信息做干涉处理可以在地震灾害评估中取得更好效果。2017 8 8 21 19 46 秒在四川阿坝州九寨沟县( 103. 82E33. 2N) 发生 7. 0 级地震,震源深度约为 20km,震中距九寨沟县 39 km( 1) 。紧接而来的降雨又造成灾区发生大面积的崩塌、滑坡等次生地质灾害。在灾区余震与次生地质灾害不断发生的情况下,遥感成为获得震区灾情的主要手段。地震发生后第一时间,中国地质环境监测院根据预报震中位置坐标,启动应急响应机制,迅速搜集查询到震区卫星数据存档影像资料,编制地震灾区震前遥感影像图( 时相: 2013 8 3 日,资源三号卫星,分辨率 2. 0 m) 。同时,与资源卫星应用中心、二十一世纪空间技术、环球星云科技、航天世景等数据公司收集地震灾区震后卫星数据,于 8 10 日晚成功获取了灾后北京二号、高分一号、高分二号等国产卫星数据,及时开展了数据处理及地质灾害解译工作。但受到震后灾区多云、阵雨等天气影响,拍摄区多被云层覆盖,对地质灾害解译有一定的影响。因此,在综合利用不同高分辨率光学卫星数据的同时,对高分辨率 SAR 遥感影像也进行了地质灾害调查应用。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 area2 主震震中位置及周边历史地震分布图Fig. 2 The location of the main earthquake and the map of the surrounding history王志一,等: 高分辨率光学和 SAR 遥感影像在地震地质灾害调查中的应用 2018 年图 3 九寨沟 G544 公路上四寨村附近滑坡震前震后遥感影像对比图Fig. 3 The comparison of remote sensing imagesfor landslides near Sizhai before and after earthquake3 地震地质灾害高分辨率雷达遥感解译3. 1 In SAR 地震形变场观测3. 1. 1 雷达干涉测量( In SA) 基本原理合成孔径雷达干涉测量( In SA) 是利用对同一地区观测的两幅 SAR 复影像数据进行相干处理,通过相位信息获取地表高程信息及形变信息的技术。合成孔径差分干涉测量( D-In SA) 技术是 In SAR 技术的拓展,利用同一地区的两幅干涉图像,其中一幅是形变前的干涉图像,另一幅是形变后获取的干涉图像,然后通过差分处理来获取地表形变的测量技术。D-In SAR 几何示意图如图 5所示,在重复轨道条件下,S1 为卫星第一次过境点,卫星轨道高度为 HS2 为卫星第二次过境点,两次过境之间的距离 B 成为基线,其与水平方向的夹角为 α。一般情况下,垂直基线即 B= Bcosα 应小于一个临界值才能获取质量较好的干涉结果。卫星两次过境与地面目标点的距离分别为 R1和 R2,如果在两次过境时间段内,对应的目标点发生了位移 Δd 时,会在卫星视线方向产生相应的位移 Δr,进而反应为电磁波的相位差,通过对相位差的计图 4 九寨沟 G544 公路如意坝附近滑坡震前震后遥感影像对比图Fig. 4 The comparison of remote sensing images forlandslides near uyiba before and after earthquake算,可以获取两景 SAR 影像期间地表发生的形变。图 5 D-In SAR 几何示意图Fig. 5 The geometrics ketch of the D-In SAR两幅影像进行干涉处理,得到干涉条纹图。干涉图相位 m包含着多种因素的相位[7],即平地相位 e、地形相位 t、地表形变相位 d及噪声相位 n,公式表示为:m= e+ t+ d+ n高分辨率光学和 SAR 遥感影像在地震地质灾害调查中的应用———以九寨沟 M7. 0 级地震为例王志一1,徐素宁1,王 娜2,马秀强3( 1. 中国地质环境监测院,北京 100081; 2. 北京环球星云遥感科技有限公司,北京 100101;3. 中国地质大学( 北京) ,北京 100083)摘要: 快速全面的获取地震灾情信息对于减少灾害损失、预防次生灾害发生具有重要的意义。遥感技术凭借其快速、大范围、高精度地获取地面信息的特点在震后信息提取中发挥了重要作用。传统光学遥感数据由于受到震后多云雨天气的影响,往往难以获取有效数据。合成孔径雷达( SA) 能够全天候获取数据,已成为震害评估的重要数据源。随着遥感图像向着高分辨率方向发展,数据所包含的信息也越来越丰富,更有利于地震地质灾害调查。本文以 2017 年九寨沟地震为例,采用高分辨率光学和 SAR 遥感图像为数据源,从道路震害分析、滑坡次生灾害提取、同震形变场计算等多个角度开展应用研究,对此次地震地质灾害进行调查。结果表明,结合高分辨率光学和 SAR 遥感影像各自的特点和优势,能够在地震地质灾害信息全面调查中发挥重要作用。关键词: 高分辨率遥感; In SA; 地震地质灾害; 应急调查中图分类号: P642; P237 文献标识码: A 文章编号: 1003-8035( 2018) 05-0081-08收稿日期: 2018-01-03; 修订日期: 2018-03-19基金项目: 国家重点研发计划课题: 城市群地质环境演化多源遥感监测与预警( 2017YFB0503803)第一作者: 王志一( 1984-) ,男,高级工程师,主要从事资源与环境遥感研究工作。E-mail: wangzy@ mailcigemgovcnApplication of the high resolution optical and SAremotesensing data images induced by the Jiuzhaigou M7. 0earthquake geological hazards surveyWANG Zhiyi1XU Suning1WANG Na2MA Xiuqiang3( 1China Institute for Geo-Environmental MonitoringBeijing 100081China; 2Beijing Earth Star emote SensingTechnology CoLTD. ,Beijing 100101China; 3China University of Geosciences ( Beijing) Beijing 100083China)Abstract: apid and comprehensive acquisition of earthquake geological disaster information is of greatsignificance to reduce the loss of disaster and prevent secondary disasters. Remote sensing technology hasplayed an important role in the extraction of earthquake information due to its ability of fastwide range andhigh precision acquisition of ground informationThe traditional optical remote sensing data are easily affectedby the cloudy and rainy weather after earthquakeso it is normally difficult for traditional optical remotesensing to obtain effective dataSynthetic aperture radar ( SA) which is capable of obtaining data in all-weather conditionhas already become an important data source for earthquake damage assessmentAs theresolution of remote sensing images are becoming higherthe information contained in the data are alsoincreasingly richerwhich is more favorable to post-earthquake investigationJiuzhaigou earthquake is taken asan examplehigh-resolution optical and SAremote sensing images are used as data source to investigate the 中国地质灾害与防治学报 ·85·干涉图中的噪声相位,可以通过多视、滤波去噪等步骤处理。干涉图中的平地相位,也称参考趋势面相位,在干涉图中占主导地位。利用轨道基线求出影像重叠区的同名像点的多个参考斜距差,选择 N 介多项式模型描述可以去除平地相位[8],生成新的干涉图。干涉图中的地形相位,是由于地形起伏引起的相位差,为了获取地表的形变信息,必须将地形相位去除。通常采用外部的 DEM 数据结合轨道参数去除地形相位,即“二轨法”,最终得到同震形变干涉图[9]。根据 D-In SAR 测量的基本原理,“二轨法”的主要步骤包括: 主辅影像预处理、影像配准、重采样、干涉图生成、滤波、形变相位分离、相位解缠、地理编码等,具体流程如图 6 所示。图 6 D-In SAR“二轨法”处理流程Fig. 6 The process flow of D-In SAR “two-track3. 1. 2 九寨沟地震 In SAR 形变场观测结果地震发生后的 8 10 日晚在环球星云德清遥感卫星地面接收站成功获取了灾后第一组同震重轨干涉SAR 影像对,利用 RADASAT-2 震前( 2017 5 30) 、震后( 2017 8 10 ) 的数据,外部 DEM 采用STM( 30 m) 数据,基于 D-In SAR“二轨法”处理方法,最终得到此次地震的同震形变场,结果如图 7 左图所示。同震形变场的局部细节如图 7 右图所示,可以看出断层东侧失相干比较严重,西侧约有 4 个干涉条纹,最外侧条纹距离中国地震台网中心确定的震中心约为36. 64km,从而可以大致推断出此次地震形变场的范围。对 In SAR 获取的同震形变场进一步反演,可以为理解此次地震的发震机制以及未来地震灾害演化趋势提供重要的依据。3. 2 雷达影像地震地质灾害解译3. 2. 1 雷达影像地震地质灾害解译基本原理相对于光学遥感图像受天气和光照条件的限制来说,雷达遥感影像具有全天候、全天时、不受云层覆盖影响等优点。地震发生后,需要对大范围的灾情信息进行宏观掌握,圈定极震区,以及对各类震害快速定位,如地质灾害、次生灾害、建筑物震害等信息,由于SAR 图像是侧视成像,当地形起伏有变化时,就会产生透视收缩、叠掩和雷达阴影等像点位移,因此利用SAR 图像进行震害信息提取的首要工作是影像图像几何校正与配准。利用 SAR 图像进行地质灾害识别主要是利用SAR 图像的幅度信息,对比分析、解译灾害发生前后的影像,从而提取灾害发生的位置、分布、面积大小等信息。其关键技术包括: SAR 图像滤波辐射校正等预处理、灾害前后影像精配准、构造差异图像、地理编码和变化区域提取等。对 SAR 图像进行几何纠正和配准主要是为 SAR 图像赋予准确的空间坐标信息,消除由 SAR 斜距成像和地表地型起伏等因素引起的几何变形,将 SAR 图像重投影至指定的坐标系统中,便于 SAR影像中灾情信息的定位以及与其它数据进行分析对图 7 In SAR 提取九寨沟地震形变场及局部图Fig. 7 Jiuzhaigou seismic deformation field and local enlarged diagram extracted by In SAR王志一,等: 高分辨率光学和 SAR 遥感影像在地震地质灾害调查中的应用 2018 年比[10]。基于雷达影像地震地质灾害解译流程如图8 所示。图 8 基于雷达影像地质灾害解译流程图Fig. 8 Flow chart of geological disasterinterpretation based on radar image9 九寨沟地震滑坡 SAR 监测图Fig. 9 SAmonitoring diagram for landslidescoused by the Jiuzhaigou earthquake3. 2. 2 九寨沟地震雷达影像地质灾害解译结果此次地震造成了多处山体滑坡,造成了震区道路多处损毁,影响地震应急救援的效率。图 9 显示震区一处滑坡造成了道路损毁的图像变化,从地震前( 9a) 和地震后( 9b) 的图像对比可以很明显的看出,震前图像上的道路明显,具有连续性; 震后道路由于受到山体滑坡的影响,造成道路中断,暗色调原有的道路被碎石覆盖,散射特性发生了变化。地震发生后,九寨沟景区内部分景点损毁,其中的火花海、诺日朗瀑布等景点损坏严重,几乎“消失”。火花海水位整体下降,已经见底,湖边发生小面积滑坡,也裸露出了黄土。如图 10 所示,从地震前( 10a) 和地震后( 10b) 的 RADASAT-2 图像对比可以看出,火花海湖面原有的暗色调水体在震后消失,湖水消失见底。图 10 九寨沟地震湖泊损毁 SAR 监测图Fig. 10 SAmonitoring diagram for lakedamaged by the Jiuzhaigou earthquake从地震前( 11a) 和地震后( 11b) 的光学卫星影像对比也可以明显看出,红框内火花海湖面原有水体在震后消失,湖水消失见底。4 讨论与结论近年来,国产高分辨率遥感数据随着数据源的不断 中国地质灾害与防治学报 ·87·图 11 九寨沟景区火花海震前震后光学遥感影像对比图Fig11 The comparison of remote sensing imagesfor landslides before and after the Jiuzhaigou earthquake丰富,数据质量的不断提升,在国民经济社会发展的各行各业均得到了大力的应用,尤其在一些特大自然灾害应急工作中,更是发挥了极其重要的作用,为抗震救灾及灾害评估提供了可靠的数据支撑。同时,在地震发生往往受到阴雨天气的影响,造成光学遥感影像在救灾决策中的作用受到了很大的限制的情况下,SAR 图像由于其全天候的特点也成为了地震发生后重要的遥感数据源。利用光学和 SAR 遥感数据,特别是高分辨率遥感影像,对全方面掌握和理解地震灾情具有十分重要的意义。本次研究采用高分辨率光学遥感数据与 SAR 数据结合的分析解译技术,利用多时相光学和雷达数据对受地震影响较大的渣扎镇周边地区做了地震灾害识别,同时在利用雷达干涉处理识别地震破裂带方面取得了较好的效果。研究表明,利用震前震后的光学影像和 SAR 幅度图像进行变化监测,能够在山体滑坡、堰塞湖堵塞等震害信息提取取发挥重要的作用。通过结合光学影像和SAR 影像各自的特点和优势,,快速进行地震地质灾害的应急调查,将成为未来地面应急调查以及灾害评估的主要数据和技术支撑手段。参考文献:1] 王志一,甘甫平,李贤庆,等. CBES 02B 星数据在“5·12”地震灾害监测中的应用———以唐家山堰塞湖动态监测为例[J]. 国土资源遥感,200921( 1) : 89 92WANG ZhiyiGAN FupingLI Xianqinget alTheapplication of CBES-02B data to the disastermonitoring of the 5. 12 Wenchuan Earthquake: Thedynamic monitoring of the Tangjiashan barrier lakeJ].Remote Sensing for Land & Resources200921 ( 1 ) :89 92.2] 刘云华,屈春燕,单新建,等. SAR 遥感图像在汶川地震灾害识别中的应用[J]. 地震学报,201032( 2) :214 223.LIU YunhuaQU ChunyanSHAN Xinjianet alApplication of SAdata to damage identification of theWenchuan earthquakeJ]. Acta Seismologica Sinica201032( 2) : 214 223.3] 王晓青,魏成阶,苗崇刚,等. 震害遥感快速提取研究: 2003 2 24 日巴楚 - 伽师 6. 8 级地震为例[J]. 地学前沿,200310 ( S1) : 285 291.WANG XiaoqinWEI ChengjieMIAO ConggangetalThe Extraction of seismic damage from remotesensing imagesA case study of Bachu-Jiashiearthquake with Ms = 6. 8 occurred on Feb24J].Earth Science Frontiers200310( S1) : 285 291.4] 魏成阶,先博勤,张宗科,等. 地震灾害航空遥感快速调查技术研究———以唐山地震区作为模拟试验场[M/ /何建邦,田国良,王劲峰. 重大自然灾害遥感监测与评估研究进展. 北京: 科学出版社,1993.WEI ChengjieXIAN BoqinZHANG Zongkeet alStudy on fast investigation of earthquake disasters fromairborne remote sensing M/ / HE JianbangTIANGuoliangWANG JinfengThe Progress of monitoringand evaluation of natural disasters by remote sensingBeijing: Science Press1993.[5] 中国地震局. 区域构造图[EB/OL]. [2017 08 10http: / / wwwceagovcn / publish / dizhenj /468 /553 /101710 /101720 /201708082344275 23689626 /indexhtmlChina Earthquake Administration. Region geologicalstructure mapEB / OL]. [2017 08 10]. http: / /wwwceagovcn / publish / dizhenj /4 68 /553 /101710 /·88· 王志一,等: 高分辨率光学和 SAR 遥感影像在地震地质灾害调查中的应用 2018 101720 /201708082344275 23689626 / indexhtml.6] 中国地震局. 主震震中位置及基本参数[EB/OL].[2017 08 08].China Earthquake AdministrationLocation and basicparameters of main earthquake epicenter EB / OL].[2017 08 08].[7] 张景发,龚丽霞,姜文亮. PS In SAR 技术在地壳长期缓慢形变监测中的应用[J]. 国际地震动态,2006( 6) : 1 6.ZHANG JingfaGONG LixiaJIANG WenliangApplication of PS In SAtechnique to measurement oflong-term crustal deformationJ]. Recent Developmentin World Seismology2006( 6) : 1 6.8] 刘国祥. In SAR 数据处理及关键算法[J]. 四川测绘,200528( 1) : 45 47.LIU GuoxiangData reduction and key algorithms inIn SAR[J]. Surveying and Mapping of Sichuan200528( 1) : 45 47.9] 王微,薛腾飞,张景发,等. 基于 In SAR 同震形变场的烈度评估研究[C/ /地壳构造与地壳应力文集.2017: 117 130.WANG WeiXUE TengfeiZHANG Jingfaet al.Research on seismic intensity evaluation based onIn SAcoseismic deformation fieldC/ / Bulletin of theInstitute of Crustal Dynamics2017: 117 130.10] 袁文龙. 星载 SAR 图像几何校正原理及软件设计与实现[D]. 北京: 中国科学院研究生院( 电子学研究所) 2005.YUAN WenlongTheory and software realization ofsatellite SAimage geo-rectification D]. Beijing:Graduate University of Chinese Academy of Sciences( Institute of Electrics) 2005.( 上接第 80 )4] 王鹏. 地基 SAR 干涉测量原理及其形变监测应用研究[J]. 测绘信息与工程,201237( 4) : 22 25 + 28.WANG PengGround-based SAinterferometryprinciples and its applications to displacementmonitoringJ]. Journal of Geomatics201237 ( 4 ) : 2225 + 28.5] 吕刚. 贵州重大地质灾害及影响因素分析[J]. 贵州地质,201633( 2) : 108 112.LYU GangAnalyses of important geo-disastersdistribution rules and influence factors of GuizhouJ].Guizhou Geology201633( 2) : 108 112.6] 李强. PSIn SAR 技术在地表沉降监测中的应用研究[D]. 太原: 太原理工大学,2012.LI QiangApplication of PSIn SAin surface subsidencemonitoring D ]. Taiyuan: Taiyuan University ofTechnology2012.7] 陈怡曲. 基于 In SAR 的形变监测技术研究[D]. 成都:电子科技大学,2013.CHEN YiquThe study of deformation monitoring basedon SAinterferometry D]. Chengdu: University ofElectronic Science and Technology of China2013.

[返回]
上一篇:基于Sentinel-1ASAR的湖南省宁乡市洪水淹没监测
下一篇:中国出口农产品质量测度