基于PS—InSAR技术的丽江断裂形变监测研究
开篇：润墨网以专业的文秘视角，为您筛选了一篇基于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向的形变测量结果。从图3分析，拉市海附近有比较明显的LOS向沉降形变，年平均形变速率最大可达-96 mm/a，根据其所处的地理环境，推测如此明显的LOS向形变最有可能与拉市海周围的地面沉降有关，而并非受断层活动影响。另一个明显的LOS向沉降形变位于蛇山公园西南侧，形变影响范围只限于局部，没有形成比较明显的线性分析特征，推测其可能也与地面沉降有关。因此，PSInSAR测量的地表形变场不仅仅与活动断裂有关，还与地面沉降影响有关，是一个多因素共同作用的结果。
为了分析丽江断裂的形变特征，沿丽江断裂可见空间上不连续的形变梯度带，尤其在狮子山与打那子附近的形变界线最为明显。但是，狮子图3叠加在DEM上的平均速率
Fig3The mean velocity superimposed
on the DEM in LOS direction图4平均速率标准偏差
Fig4The standard deviation of mean
velocity in LOS direction山南侧沿丽江断裂NE向存在形变界线外，其西侧沿NNE方向似乎也存在一条长2～3 km的形变界线一直延伸至象山（图4）。通过仔细验证，发现狮子山以东、以北地区主要为丽江古城遗址，而且古城遗址与狮子山附近的形变界线基本吻合；丽江古城自2003年以来就一直处于受保护状态，地表变形相对稳定，而其它城区可能因地下水开采引起地面沉降。所以，狮子山南侧沿丽江断裂NE向的形变界线并不能完全代表丽江断裂的活动，它应该是活动断裂与地面沉降的综合结果；而狮子山西侧沿NNE向的形变界线可能主要与地面沉降有关。所以，丽江断裂进入城区后由于城市化活动影响，PSInSAR观测到的形变界线模糊或存在多解性。
打那子附近的形变界线可能与丽江断裂的活动是紧密相关的，因为此形变界线附近的人为活动影响弱，基本上可以排除城市化造成的地面沉降影响。垂直于丽江断裂走向进行形变速率剖面线分析（图5）可知，丽江断裂NW盘与SE盘的LOS向年变速率相差约3～5 mm/a，NW盘表现为相对远离卫星方向运动，南东盘表现为相对靠近卫星方向运动，反映了丽江断裂左旋走滑的特征。图5丽江断裂形变剖面线
Fig5Deformation section line of the Lijiang FaultPSInSAR揭示的2003～2010年丽江断裂的年变速率约为3～5 mm/a，与向宏发等（2002）给出的历史滑动速率基本相当，表明丽江断裂目前地壳活动相对稳定。
3结论
基于PSInSAR形变探测原理，运用小基线集干涉处理方法提取、分析了丽江断裂2003～2010年的平均形变速率。结果表明：丽江断裂具左旋走滑形变特征，年平均LOS向滑动速率可达3～5 mm左右，目前处于相对稳定状态。丽江城区局部以及拉市海附近存在较明显的地面沉降，最大沉降速率可达-96 mm/a。
本文的研究为丽江地区地震预测增加了可靠的断裂活动数据，有助于推进PSInSAR技术在活动断裂形变监测中的进一步应用。但是丽江断裂进入城区后，受城市化干扰因素影响，PS-InSAR测量获得的地表形变信息中存在地面沉降的干扰，目前尚没有好的方法进行分离。
感谢欧空局（ESA）提供的ASAR科学数据服务（编号C1P5623）。
参考文献：
陈国浒，刘云华，单新建2010PSInSAR技术在北京采空塌陷区地表形变测量中的应用探析[J].中国地质灾害与防治学报，21（2）：59-63
郭炳跃，何敏，刘建东2011利用PSInSAR技术监测南通市区地面沉降[J].地质灾害与环境保护，22（4）：103-107
韩竹军，向宏发，虢顺民，等2004滇西北丽江盆地北部区第四纪时期的左旋剪切拉张[J].科学通报，50（4）：356-362
洪顺英2010基于多视线向DInSAR技术的三维同震形变场解算方法研究及应用[D].北京：中国地震局地质研究所
李忠华，苏有锦，刘祖荫1998丽江70级地震滇西实验场区跨断层短基线短水准测量的变化特征[J].中国地震，14（2）：140-149
牛安福，张晶，薛艳2004丽江地震后形变场变化与中强地震之间的关系的研究[J].地震，24（增刊）：61-68
屈春燕，单新建，宋小刚，等2011基于PSInSAR技术的海原断裂带地壳形变初步研究[J].地球物理学报，54（4）：984-993
苏生瑞，王运生，王士天2004丽江地震区应力场研究[J].地质论评，50（1）：57-63
王永安，刘强，王世芹2004丽江70级地震前跨断层形变累积率的变化特征[J].地震研究，27（1）：61-65
向宏发，徐锡伟，虢顺民，等2002丽江—小金河断裂第四纪以来的左旋逆推运动及其构造地质意义[J].地震地质，24（2）：188-198
许才军，何平，温扬茂2011利用PSInSAR研究意大利Etna火山的地表形变[J].武汉大学学报：信息科学版，36（9）：1 012-1 016
张景发，龚利霞，姜文亮2006PSInSAR技术在地壳长期缓慢形变监测中的应用[J].国际地震动态，（6）：1-6
周光全，张建国，周瑞琦，等1997丽江70级地震的地震地质构造背景分析[J].地震研究，20（1），92-100