Monitoring vertical ground deformation in the North China Plain using the multitrack PSInSAR technique
-
摘要: 针对PSInSAR技术监测大空间尺度地表垂直形变场的不足,提出通过基准转换与数据拼接、融合不同轨道多个独立片区(Frame)获取大空间尺度地表垂直形变场方法.由于不同轨道主影像选取的参考基准不一致、不同轨道垂直基线不同、不同片区相位解缠的参考PS不同等因素影响,使得PSInSAR获得的不同轨道地表垂直形变场的空间基准存在差异,提取的地表形变信息不具有可比性.利用两个相邻片区重叠区域的PS点集,根据迭代最近点算法提取两片区公共区域内同名点,并兼顾同名点距离的权,精确计算两片区的配准参数,实现两片区的"无缝"连接,最后使各独立片区拼接成一个连续的大空间尺度地表垂直形变场.本文推导了多轨PSInSAR技术获得的不同片区形变场基准转换与数据拼接的数学模型,并以华北平原(115.32°E—118.79°E, 36.81°N—40.58°N)为实验区,解析了3个不同轨道共12个片区PS目标,获取了研究区2006—2010年地表垂直形变速度场.分析表明:(1)研究区大范围处于下沉状态,形成以北京、廊坊、天津、沧州、泊头—德州等城市为中心向外扩展的地表沉降发展态势,几个沉降中心的平均沉降速率分别达到-34.7 mm/a、-26.3 mm/a、-64.2 mm/a、-34.6 mm/a和-37.7 mm/a;(2)地表沉降的空间分布特征表明,城市工业生产和生活抽取地下水是地表沉降的主要诱因,农业灌溉和油气开采是导致华北平原大范围地表沉降的另一个重要因素;(3)研究区最大沉降带沿北北东向展布,与区内断裂分布一致,表明沉降的空间分布受到断裂带控制;(4)利用地表同期一等水准测量成果检验了本研究结果,精度达到4.72 mm,表明本文提出的数据处理策略能够满足大时空尺度地表形变监测的需求.Abstract: In this paper, we present a method to reveal vertical deformation characteristics through reference transformation, data splicing and data fusion in multiple regions, which can overcome the shortage of traditional single-track PSInSAR techniques. In order to solve the problem that the ground deformation is not comparable, which is derived from various methods with reference difference in multi-tracks, vertical baseline difference in multi-tracks and reference PS difference in phase unwrapping, we calculate accurate transformation parameters and realize seamless connection of adjacent two frames through corresponding PS point analysis of the overlapped region and same point identification using the shortest distance algorithm based on proper weight. On the basis of algorithm research, we deduce the model to extract deformation information through data connection and reference transformation in multi-track PSInSAR technology, and then process the data collected in the North China Plain(115.32°E—118.79°E, 36.81°N—40.58°N) which include 3 different orbits and 12 frames, and obtain the vertical deformation of the study area for the period 2006—2010. The results indicate that the research area is in a sinking status, forming several subsidence centers including Beijing, Langfang, Tianjin, Cangzhou and Botou-Dezhou with average subsidence rates up to -34.7 mm/a, -26.3 mm/a,-64.2 mm/a, -34.6 mm/a and -37.7 mm/a, respectively. The spatial distribution of subsidence shows that groundwater extraction of industry and daily life are the main reasons, meanwhile agricultural irrigation and gas exploitation are the other primary causes. The major subsidence zones, trending NNE, coincide with the orientations of the fault zones, implying the subsidence distribution is controlled by the faults. These research results are verified by leveling measurements with errors about 4.72 mm, which shows that the data processing method in this paper can meet the requirement of large spatial scale deformation monitoring.
-
[1] Adam N, Gonzalez F R, Parizzi A, et al. 2011. Wide area persistent scatterer interferometry: algorithm and examples. Fringe 2011. Frascati, Italy: IEEE, 1857-1860.
[2] Beijing Institute of Surveying and Mapping. 1999. Code for Urban Survey (CJJ 8-99). Beijing:China Building Industry Press(in Chinese).
[3] Casu F, Manzo M, Pepe A, et al. 2008. SBAS-DInSAR analysis of very extended areas: first results on a 60000 km2 test site. IEEE Geoscience and Remote Sensing Letters, 5(3): 438-442.
[4] Chang C P, Yen J Y, Hooper A, et al. 2010. Monitoring of surface deformation in northern Taiwan using DInSAR and PSInSAR techniques. Terrestrial, Atmospheric & Oceanic Sciences, 21(3): 447-461.
[5] Cuenca M C, Hanssen R, Hooper A, et al. 2011. Surface deformation of the whole Netherlands after PSI analysis. // Fringe 2011. Frascati, Italy.
[6] Decriem J, árnadóttir T, Hooper A, et al. 2010. The 2008 May 29 earthquake doublet in SW Iceland. Geophysical Journal International, 181(2): 1128-1146.
[7] Earth. 2012. Land subsidence in the Fen Wei Basin. http://www.diqiuzazhi.com.
[8] Ferretti A, Prati C, Rocca F. 2001. Permanent scatterers in SAR interferometry. IEEE Transactions on Geoscience and Remote Sensing, 39(1): 8-20.
[9] Ferretti A, Prati C, Rocca F. 2002. Nonlinear subsidence rate estimation using permanent scatterers in differential SAR interferometry. IEEE Transactions on Geoscience and Remote Sensing, 38(5): 2202-2212.
[10] Funning G J, Bürgmann R, Ferretti A, et al. 2007. Creep on the Rodgers Creek fault, northern San Francisco Bay area from a 10 year PS-InSAR dataset. Geophysical Research Letters, 34(19): L19306.
[11] Ge D Q, Zhang L, Wang Y, et al. 2010. Merging multi-track PSI result for land subsidence mapping over extended area. // Proceeding of IEEE Geoscience and Remote Sensing Symposium (IGARSS 2010). Honolulu, Hawaii, USA, 3522-3525.
[12] Guo B Y,He M,Liu J D. 2011. Monitoring the ground subsidence in nantong city with the ps-insar technology.Journal of Geological Hazards and Environment Preservation.(in Chinese),22(4):103-107.
[13] Guo K Y,Yu Z J,Fang Z,et al.2004. Survey and Evaluation of groundwater resources and geological hazard in the Yangtze River Delta area. Seminar of Geological environment of coastal zones and city development.(in Chinese).2004,Tianjin.
[14] He Q C. 2005. Land Subsidence in the North China Plain(NCP). // Proc. the 7th International Symposium on Land Subsidence (in Chinese). 18-29.
[15] He Q C, Liu W B, Li Z M. 2006. Land Subsidence Survey and Monitoring in the North China Plain. Geological Journal of China Universities.(in Chinese), 12(2):195-209.
[16] Hooper A, Zebker H, Segall P, et al. 2004. A new method for measuring deformation on volcanoes and other natural terrains using InSAR persistent scatterers. Geophysical Research Letters, 31(23), doi: 10.1029/2004GL021737.
[17] Hooper A, Segall P, Zebker H. 2007. Persistent scatterer interferometric synthetic aperture radar for crustal deformation analysis, with application to Volcán Alcedo, Galápagos. Journal of Geophysical Research: Solid Earth (1978-2012), 112(B7), doi: 10.1029/2006JB004763.
[18] Jiang Z S. 2003.Horizontal strain field and tectonic deformation of china mainland revealed by preliminary gps result. Chinese J.Geophys.(in Chinese), 46(3): 352-358.
[19] Luo S M,Dong Y H,Han Y P,et al. 2012a. Comparison of Adjustment Results Both of Dynamics and Statics for Repeated Levelling Networks. Journal of Geodesy and Geodynamics.(in Chinese), 32(1):97-100+122.
[20] Luo S M, Du K F, Chang L, et al. 2014.Ground subsidence rates of Beijing area constrained by PS-InSAR analysis. Journal of Geodesy and Geodynamics. (in Chinese), 34(1): 43-46.
[21] Luo S M,Yang G H,Dong Y H,et al. 2012b. Displacement field characteristics on L'Aquila Earthquake from PSInSAR time series. Progress in Geophysics.(in Chinese), 27(4):1307-1315.
[22] Qu C Y, Shan X J, Song X G, et al.2011.The PSInSAR technique and its application to the study on crustal deformation of the Haiyuan fault zone. Chinese J.Geophys. (in Chinese), 54(4): 984-993.
[23] Yang G H, Jiang Z S, Liu G Y, Han Y P. 2007. Possible relation of horizontal movement field in north china to kunlun mountain MS8.1 earthquake. Journal of Geodesy and Geodynamics.(in Chinese), 27(2):10-15.
[24] Zhang Y G,Shuai P, Hu X K,et al. 2006.Study on evolution of deformation field in north china area with gps data. Journal of Geodesy and Geodynamics.(in Chinese), 26(1):36-41.
计量
- 文章访问数: 1731
- PDF下载数: 2371
- 施引文献: 0