基于PS—InSAR技术丽江断裂形变监测研究.docVIP

基于PS—InSAR技术丽江断裂形变监测研究 　　摘要：基于PSInSAR技术，利用欧空局28景ENVISAT ASAR数据，运用小基线集干涉处理方法，提取并分析了丽江断裂2003～2010年的平均形变速率。研究结果显示：（1）丽江断裂具左旋走滑形变特征；（2）年平均LOS向形变速率在3～5 mm/a之间，与地质资料较吻合。该研究为丽江地区地震预测提供一种可靠的断裂带活动数据，但丽江断裂进入城区后受城市化活动影响，断裂形变与地面沉降信息共存，较难分离。 　　关键词：PSInSAR；丽江断裂；形变速率；LOS向 　　中图分类号：P315727文献标识码：A文章编号：1000-0666（2013）03-0319-05 　　0前言 　　以丽江市区为界，丽江—金河断裂分为南西和北东两段。该断裂带晚第四纪以来表现出较强的活动性。在丽江—永胜公路12 km处可见该断裂东支断层的地质剖面，断层明显切错古近系上部丽江组砾岩。断面产状NE28°/NW47°，其中发育侧伏向SW，倾伏角35°～40°的擦痕，判断该断层的最新错动为左旋走滑—正断性质；对取自该断层的样品的热释光绝对年龄测定，其最新活动发生在距今3～21 ka。沿断裂带地震活动频繁，在区内长120 km区段内，自1481年以来共记录M≥40地震20余次，较强的两次地震为1751年5月25日剑川6级地震以及1996年2月3日丽江7级地震；这两次地震的等震线长轴方向均为NE向，受丽江—小金河断裂控制，形成一条横贯研究区东南部的地震活动带（苏生瑞等，2004；韩竹军等，2004；周光全等，1997）。 　　对研究区断裂带方向及活动性（速率）的测量，为定量化研究断裂构造对地应力场的影响以及揭示断裂带动力学机制提供必要的数据。虽然王永安等（2004）、牛安福等（2004）、李忠华等（1998）在丽江地区及其邻区开展过地应力及水准测量工作，但由于测点分散、测量深度相差较大，测量值和方向都较为分散。永久散射体雷达干涉测量（Permanent ScatterInterferometric Synthetic Aperture Radar，简称PSInSAR）技术的优点在于解决时间失相关的同时，计算并消除大气、地形影响与轨道误差，保证干涉处理的正常进行，并可以获取毫米级精度的地表形变信息。PSInSAR作为空间探测的新技术，为高精度监测断层微形变速率提供了时间域的连续动态??测和空间域的大量测点，为定量化计算研究区的应力场提供了必要的、精确的数据，在地震、地质灾害监测方面得到了广泛应用（Ferretti et al，2000，2001；张景发等，2006；陈国浒等，2010；郭炳跃等，2011；屈春燕等，2011；许才军等，2011）。本文应用PSInSAR技术解析了丽江断裂微形变速率及特征。 　　1PSInSAR数据处理方法111数据选取利用PSInSAR技术监测地壳形变需要同一地区大量的影像对，以保证较高的计算精度。目前常用的卫星数据有：欧洲ENVISAT、加拿大的RADARSAT、日本ALOS及德国的TerrSARX卫星。ENVISAT卫星由欧空局2002 年发射，搭载的雷达传感器ASAR工作波段波长为56 cm（C波段）。考虑天气状况、轨道参数及数据可获得性等因素，我们获取丽江地区28景ENVISAT卫星降轨IS2模式，VV极化SLC数据（表1），外部DEM数据采用美国宇航局喷气推进实验室SRTM3′数据（约90 m分辨率）。ENVISAT卫星的轨道数据采用DORIS精密轨道，以提高初始干涉测量的轨道精度。1地震研究136卷第3期刘智荣等：基于PS-InSAR技术的丽江断裂形变监测研究表1ASAR数据（垂直基线距由DORIS精密轨道计算） 　　Tab1ASAR Data（vertical baseline distance 　　2数据处理方法 　　因为活动断裂所处的地区地势往往较复杂（图1），东部与西南地区植被???响较严重，而且活动断裂的地表形变并不是严格的呈现线性形变特征，可能在某些时间内会有突然的变形作用。综合以上因素，本文应用斯坦福大学开发的Standford Method for Persistent Scatter（StaMPS）（Hooper，2006，2008）对丽江断裂形变特征进行提取。StaMPS方法的最大优势是可以在植被区或山区获取比较可靠的PS点，特别适用于活动断裂所处的地质与地理条件，而且它不需要先验的形变模型，就可以捕捉非线性形变信息，计算精度也可以满足活动断裂形变监测的要求。 　　22丽江断裂形变特征分析 　　由于雷达卫星干涉测量存在视线向模糊问题（洪顺英，2010），基于PSInSAR获取的年平均形变速率图反映的是卫星LOS向的形变测量结