提取地震同震形变场特征的最佳雷达波段问题研究

摘要/Abstract

摘要： 利用弹性形变模型, 并以走滑型断层为例, 对地震同震形变场进行了模拟研究。结合雷达成像几何关系, 得到干涉雷达的LOS向形变场以及干涉相位条纹图像。对比分析了不同雷达波段(X波、 C波、 L波)以及不同雷达入射角(30°和50°)的同震形变场和干涉相位图。结果表明, 同一雷达波段不同雷达入射角的LOS向形变场会有不同的效果, 入射角逐渐增大时, LOS向形变量随之增大, 体现了发震断层走滑特征; 采用不同雷达入射角, 其探测到的地表形变垂向分量和水平分量信息量不同, 表明入射角较大或较小时能够分别较好地对地表形变的水平变化和垂向变化进行监测; 同一雷达入射角, 不同雷达波段的干涉相位条纹有明显区别, 波长越短, 干涉条纹越密集, 容易出现去相关现象, 长波长的干涉相位条纹比较清晰。从某种意义上讲, 尽管长波长的干涉相位反映地表形变信息的细节会减少, 但受噪声影响较小, 抗干扰性能较短波长的强, 干涉像对相干性高。地表形变表现为长期缓慢的小幅度渐变性形态, 也表现为短期迅速的大面积突发性形态。不同雷达波长探测到的地表形变有效信息不同, 为了更好地利用干涉雷达对各种尺度和不同变化速率的地表形变进行监测, 需要在满足精度的条件下尽可能地选择波长较长及多样化的雷达数据。

关键词: 同震形变, 线弹性模型, 干涉雷达, 雷达波段, 随机误差

Abstract: In this paper,we modeled earthquake co-seismic deformation field using a linear elastic model in half-space with a strikeslip earthquake fault. Afterward, we got coseismic deformation field in LOS direction, descending and ascending pass InSAR results with two kinds of radar imaging geometry respectively, and made comparative analysis of results with different SAR bands (X, C and L) and different SAR incidence angles (30°and 50°). The results show that, in the same band, the InSAR results contains different effective information with different incidence angles. The deformation field in LOS direction increases with the increase of incidence angle, reflecting the characteristics of the seismic fault. The detected vertical component and horizontal component of the surface deformation by InSAR differ with different incidence angle. Larger or smaller incidence angle can be used to well detect horizontal deformation or vertical deformation respectively. In the same incidence angle, fringes of the InSAR results are obviously different with different wavelengths. With shotee wavelength, the rate of the fringe is higher, resulting in incoherence phenomena. In contrast, the fringe with longer wavelength is relatively clear. In a sense, using longer wavelength may result in more usable interferometric pairs, and the InSAR result are less affected by noise, although it may contain less precise crustalmotion measurements. The surface deformation is not only embodied as long-term, slow and gradual change, but also embodied as shortterm, rapid and large change. Radar with different wave length can detect different effective information about surface deformation. In order to better detect surface deformation with different rate and various scale of changes by using InSAR, it is necessary to use the longest radar wavelengths possible within some limits, and multiple observations.

Key words: Co-seismic deformation, Linear elastic model, SAR interferometry, Radar band, Random noise

中图分类号:

P315

严珍珍, 石耀霖, 范湘涛, 杜小平. 提取地震同震形变场特征的最佳雷达波段问题研究[J]. 地震, 2011, 31(3): 46-55.

YAN Zhen-zhen, SHI Yao-lin, FAN Xiang-tao, DU Xiao-ping. Study on Optimal SAR Band for Mapping Co-seismic Deformation[J]. EARTHQUAKE, 2011, 31(3): 46-55.

参考文献

[1] Massonnet D, Rossi M, Cesar C. The displacement field of the Landers earthquake mapped by radar interferometry[J]. Nature, 1993, 364(6433): 138-141.
[2] 张景发, 刘钊. InSAR技术在西藏玛尼强震区的应用[J]. 清华大学学报(自然科学版). 2002, 42(6): 847-850.
[3] 孙建宝, 石耀霖, 沈正康, 等. 基于线弹性位错模型反演1997年西藏玛尼M_W7.5地震的干涉雷达同震形变场-Ⅱ滑动分布反演[J]. 地球物理学报, 2007, 50(5): 1390-1397.
[4] 孙建宝, 徐锡伟, 沈正康, 等. 基于线弹性位错模型及干涉雷达同震形变场反演1997年玛尼M_W7.5地震参数-Ⅰ均匀滑动反演[J]. 地球物理学报, 2007, 50(4): 1097-1110.
[5] 单新建, 马瑾, 王长林. 利用星载D-INSAR技术获取的地表形变场提取玛尼地震震源断层参数[J]. 中国科学(D辑), 2002, 32(10): 837-847.
[6] 王超, 刘智, 张红. 张北-尚义地震同震形变场雷达差分干涉测量[J]. 科学通报, 2000, 45(23): 2550-2555.
[7] 单新建, 柳稼航, 马超. 2001年昆仑山口西8.1级地震同震形变场特征的初步分析[J]. 地震学报, 2004, 26(5): 474-480.
[8] 孙建宝, 梁芳, 徐锡伟, 等. 升降轨道ASAR雷达干涉揭示的巴姆地震(M_W6.5)3D同震形变场[J]. 遥感学报, 2006, 10(4): 489-496.
[9] Wang R, Xia Y, Grosser H. The 2003 Bam (SE Iran) earthquakes: precise source parameters from satellite radar interferometry[J]. Geophys J Int, 2004, 112(6): 1-6.
[10] 孙建宝, 梁芳, 沈正康, 等. 汶川M_S8.0地震InSAR形变观测及初步分析[J]. 地震地质, 2008, 30(3): 789-795.
[11] 姚鑫, 张永双. 基于差分干涉雷达的汶川地震同震形变特点[J]. 地质力学学报, 2009, 15(2): 151-161.
[12] Yuan X, Liu Y, A Yuen D. An efficient system for creating synthetic InSAR images from simulations[J]. Pure Appl Geophys, 2008, 165: 671-691.
[13] Okada Y. Surface deformation due to shear and tensile faults in a half-space[J]. Bulletin of the Seismological Society of America, 1985, 75(4): 1135-1154.
[14] Wang R, Martin F L, Roth F. Computation of deformation induced by earthquakes in a multi-layered elastic crustFORTRAN programs EDGRN/EDCMP[J]. Computers & Geosciences, 2003, 29: 195-207.
[15] 穆东. 干涉合成孔径雷达成像技术研究[D]. 南京: 南京航空航天大学, 2001.
[16] 何平, 许才军. 卫星轨道误差对SAR干涉处理的影响研究[J]. 大地测量与地球动力学, 2009, 28(5): 54-57.
[17] Zebker H A. Atmospheric effects in interferometric systhetic aperture radar surface deformation and topographic maps[J]. Journal of Geophysical Research, 1997, 102(B4): 7547-7563.
[18] Xu Z W, Wu J, Wu Z S. A survey of ionospheric effects on space-based radar[J]. Waves Random Media, 2004, 14: S189-S273.
[19] Sandwell D T, Myer D, Mellors R. Accuracy and resolution of ALOS Interferometry: vector deformation maps of the Father′s Day intrusion at Kilauea[J]. IEEE Transactions on Geoscience and Remoting Sensing, 2008, 46(11): 3524-3534.
[20] Tarayre H, massonnet D. Atmospheric propagation heterogeneities revealed by ERS-1 interferometry[J]. Geophysical Research Letters, 1996, 23(9): 989-992.
[21] 刘艳华. 合成孔径雷达干涉(InSAR)测量技术测量地表变形的理论与算法研究[D]. 北京: 中国测绘科学研究院, 2007.
[22] 陈艳玲, 黄珹, 王小亚. 电离层对SAR干涉测量的影响综述[J]. 地球物理学进展, 2010, 25(3): 823-830.
[23] Rosen P A, Hensley S, Zebker H A. Surface deformation and coherence measurements of Kilauea Volcano, Hawaii, from SIR-C radar interferometry[J]. Journal of Geophysical Research, 1996, 101(E10): 23109-23125.