基于Sentinel-1A卫星SAR数据地震同震形变场反演——以意大利佩鲁贾地震为例
|
摘要
目的验证Sentinel-1A卫星SAR数据在地震导致的地表形变测量方面的可行性。方法利用获取到的意大利佩鲁贾市地震前后两景Sentinel-1A卫星SAR数据并使用ENVI软件反演其地震同震形变场,得到雷达视线(LOS)方向地表形变图。结果反演结果表明:(1)地震导致的雷达视线方向地表沉陷和抬升最大值分别约为0.31 m和0.32 m,地震中心地表形变量约为0.1 m;(2)地震中心附近形成"西北-东南"走向长约30km的断层;(3)由地表形变剖面图可知:位于地震中心以北、以西和以南3个方向的地表有所沉陷,地震中心以东地区地表被抬升。结论 Sentinel-1A数据在地震地表形变场反演具有可行性,可作为研究地震机理和防灾减灾参考数据,但实际精确形变值还需要GPS测量数据验证。
Objective To verify the feasibility of Sentinel-1 A satellite SAR data in the measurement of surface deformation caused by earthquakes.Menthods In this paper,we used the SAR data of Sentinel-1 A satellite before and after the Perugia earthquake and the ENVI software to retrieve the co-seismic deformation field,and obtained the surface deformation map in the direction of radar line of sight(LOS).Results The inversion results showed that:(1)the maximum values of surface subsidence and uplift caused by the earthquake in the radar line of sight direction respectively were about 0.31 m and 0.32 m,respectively and the surface shape variable of earthquake center was about 0.1 m;(2)near the center of earthquake,a fault with a strike length of about 30 km was formed in the northwest-southeast direction;(3)from the surface deformation profile,it can be seen that the surface subsides in three directions:north,west and south of the seismic center,and the surface is uplifted to the east of the seismic center.Conclusion The sentinel-1 A data is feasible in seismic surface deformation field inversion,which can be used as reference data for earthquake mechanism and disaster prevention and reduction,but the actual exact deformation value needs to be verified by GPS.
Objective To verify the feasibility of Sentinel-1 A satellite SAR data in the measurement of surface deformation caused by earthquakes.Menthods In this paper,we used the SAR data of Sentinel-1 A satellite before and after the Perugia earthquake and the ENVI software to retrieve the co-seismic deformation field,and obtained the surface deformation map in the direction of radar line of sight(LOS).Results The inversion results showed that:(1)the maximum values of surface subsidence and uplift caused by the earthquake in the radar line of sight direction respectively were about 0.31 m and 0.32 m,respectively and the surface shape variable of earthquake center was about 0.1 m;(2)near the center of earthquake,a fault with a strike length of about 30 km was formed in the northwest-southeast direction;(3)from the surface deformation profile,it can be seen that the surface subsides in three directions:north,west and south of the seismic center,and the surface is uplifted to the east of the seismic center.Conclusion The sentinel-1 A data is feasible in seismic surface deformation field inversion,which can be used as reference data for earthquake mechanism and disaster prevention and reduction,but the actual exact deformation value needs to be verified by GPS.
引文
[1] 李少青,杨武年,杨彦通,等.Sentinel-1A_IW雷达影像在地震形变信息提取中的应用[J].测绘,2016,39(2):56-59.
[2] Facts and figures of Sentinel-1A[EB/OL].(2016-04-25)[2017-03-21].http://www.esa.int/Our_Activities/Observing_the_Earth/Copernicus/Sentinel-1/Facts_and_figures.
[3] 杨久东,田辉,周春霞.基于D-InSAR技术的意大利拉奎拉地震同震形变分析[J].大地测量与地球动力学,2010,30(1):9-13.
[4] 中国地震台网中心.意大利6.6级地震[EB/OL].(2016-10-30)[2017-03-21].https://news.ceic.ac.cn/CC20161030144020.html.
[5] MASSONNET D,SOUYRIS J C.Imaging with synthetic aperture radar[M]//Imaging with synthetic aperture radar.EPFL Press,2008.
[6] MONTI GUARNIERI A,GUCCIONE P,PASQUALI P,et al.Multi-mode ENVISAT ASAR interferometry:techniques and preliminary results[J].Radar,Sonar and Navigation,IEEE Proceedings,2003,150(3):193-200.
[7] RICHARD M.GOLDSTEIN,CHARLES L.W.Radar interferogram filtering for geophysical applications[J].Geophysical Research Letters,1998,25(21):4035-4038.
[8] BARAN I,STEWART M P,KAMPES B M,et al.A modification to the Goldstein radar interferogram filter[J].IEEE Transactions on Geoscience & Remote Sensing,2003,41(9):2114-2118.
[9] GHULAM A,AMER R,RIPPERDAN R.A filtering approach to improve deformation accuracy using large baseline,low coherence DInSAR phase images[J].Geoscience & Remote Sensing Symposium,2010:3494-3497.
[10] Reigber A,Moreira J.Phase unwrapping by fusion of local and global methods[C]//Institute of Electrical and Electronics Engineers.IEEE International Geoscience and Remote Sensing Symposium Proceedings.Remote Sensing-A Scientific Vision for Sustainable Development:1997 vol.2.Singapore:IGARSS,1997:869-871.
[11] Holecz F,Moreira J,Pasquali P,et al.Height model generation,automatic geocoding and a mosaicing using airborne AeS-1 InSAR data[C]//Institute of Electrical and Electronics Engineers,IEEE International Geoscience and Remote Sensing Symposium Proceedings.Remote Sensing-A Scientific Vision for Sustainable Development:1997 vol.2.Singapore:IGARSS,1997:1929-1931.
[12] Goblirsch W,Pasquali P.Algorithms for calculation of digital surface models from the unwrapped interferometric phase[C]//Institute of Electrical and Electronics Engineers,International Geoscience and Remote Sensing Symposium:1996 vol.1.Lincoln,NE,USA:IGARSS,1996:656-658.
[2] Facts and figures of Sentinel-1A[EB/OL].(2016-04-25)[2017-03-21].http://www.esa.int/Our_Activities/Observing_the_Earth/Copernicus/Sentinel-1/Facts_and_figures.
[3] 杨久东,田辉,周春霞.基于D-InSAR技术的意大利拉奎拉地震同震形变分析[J].大地测量与地球动力学,2010,30(1):9-13.
[4] 中国地震台网中心.意大利6.6级地震[EB/OL].(2016-10-30)[2017-03-21].https://news.ceic.ac.cn/CC20161030144020.html.
[5] MASSONNET D,SOUYRIS J C.Imaging with synthetic aperture radar[M]//Imaging with synthetic aperture radar.EPFL Press,2008.
[6] MONTI GUARNIERI A,GUCCIONE P,PASQUALI P,et al.Multi-mode ENVISAT ASAR interferometry:techniques and preliminary results[J].Radar,Sonar and Navigation,IEEE Proceedings,2003,150(3):193-200.
[7] RICHARD M.GOLDSTEIN,CHARLES L.W.Radar interferogram filtering for geophysical applications[J].Geophysical Research Letters,1998,25(21):4035-4038.
[8] BARAN I,STEWART M P,KAMPES B M,et al.A modification to the Goldstein radar interferogram filter[J].IEEE Transactions on Geoscience & Remote Sensing,2003,41(9):2114-2118.
[9] GHULAM A,AMER R,RIPPERDAN R.A filtering approach to improve deformation accuracy using large baseline,low coherence DInSAR phase images[J].Geoscience & Remote Sensing Symposium,2010:3494-3497.
[10] Reigber A,Moreira J.Phase unwrapping by fusion of local and global methods[C]//Institute of Electrical and Electronics Engineers.IEEE International Geoscience and Remote Sensing Symposium Proceedings.Remote Sensing-A Scientific Vision for Sustainable Development:1997 vol.2.Singapore:IGARSS,1997:869-871.
[11] Holecz F,Moreira J,Pasquali P,et al.Height model generation,automatic geocoding and a mosaicing using airborne AeS-1 InSAR data[C]//Institute of Electrical and Electronics Engineers,IEEE International Geoscience and Remote Sensing Symposium Proceedings.Remote Sensing-A Scientific Vision for Sustainable Development:1997 vol.2.Singapore:IGARSS,1997:1929-1931.
[12] Goblirsch W,Pasquali P.Algorithms for calculation of digital surface models from the unwrapped interferometric phase[C]//Institute of Electrical and Electronics Engineers,International Geoscience and Remote Sensing Symposium:1996 vol.1.Lincoln,NE,USA:IGARSS,1996:656-658.