基于InSAR与GPS观测的汶川同震垂直形变场的获取

单新建 屈春燕 郭利民 张国宏 宋小刚张桂芳 温少妍 汪池升 徐小波 刘云华

1)中国地震局地质研究所，地震动力学国家重点实验室，北京 100029

2)中国地震局地震研究所，武汉 430071

3)新疆地震局，乌鲁木齐 830011

4)深圳大学，深圳 518060

0 引言

地表同震形变是三维形变，地表每一点都是三维形变变化。对于走滑断层，以水平形变变化为主，对于逆冲或正断层，以垂直形变变化为主。目前的合成孔径雷达干涉测量技术(In-SAR)、全球定位系统(GPS)、水准观测等手段都以某一类形变观测为主。因此，利用多手段观测数据联合解算地震同震形变场各个分量，对深刻理解断层运动方式、形变特征以及分析地震发生机理有着十分重要的作用。

InSAR是一种侧视形变测量技术，它获取的形变量是地表东西(EW)、南北(SN)、垂直(UP)3个形变分量在雷达视线向(LOS)上的投影，它观测到的地表形变量是一维的，并不能反映地表真实的形变量，这就是InSAR干涉测量中的视线向模糊问题。由InSAR观测值很难分离出地表水平和垂直方向的形变量。因此，如何从InSAR视线向形变场中分解出3个方向的形变值是一个重要研究内容。

目前在解算同震地壳三维形变场方面可以归结出以下几种方法:1)多视线向InSAR解算三维形变的方法。主要利用升降轨3个不同入射角的雷达卫星数据，获取同一地震同震形变场不同视线向形变分量，来实现三维形变场的解算(Wright et al.，2004;Hu et al.，2008)。该方法的缺点是同一地震不同视线方向的雷达数据未必在地震前后都完备;2)利用2个不同视线向InSAR形变量，结合offset偏移量方法获取的方位向形变场，解算同震三维形变场的方法(Yuri et al.，2001)。由于偏移量法是针对大形变尺度提出的方法，两者的结合在形变精度方面难以保证;3)利用2个不同视线向InSAR形变量，结合SN向模型形变场分量来解算三维形变的方法。由于InSAR观测结果对地表SN方向形变分量不敏感，可以利用初始断层模型模拟出的SN向形变分量，实现对三维形变场的解算(汪驰升等，2009)。但该方法除了受断层走向、断层形变方向限制外，由于采用了模型形变分量，使观测结果过早地掺杂了理论模拟量，可能造成观测量的不真实;4)GPS水平观测数据与InSAR视线向观测数据联合求解，获得地表同震三维形变场(Sverrir et al.，2002;Samsonov et al.，2006)。由于InSAR技术具有近场地表连续覆盖、视线向精度高的特点，已成为研究地震同震形变的重要方法之一，而GPS具有大区域、离散分布、水平精度高的观测优势，已广泛应用于各种规模尺度的地壳运动观测与构造形变研究。该方法将二者结合起来，充分利用了GPS的水平精度高以及InSAR视线向观测精度高且对垂直敏感的优势，较好地将离散的GPS观测点与地表连续观测的InSAR结合到一起。其缺点是要求GPS观测点密度较高，且需要选取适合的插值方法。

汶川地震所处的发震断层——龙门山断裂带为大规模的逆掩推覆构造变形带，而汶川地震发震断层以逆冲为主(邓起东等，1994;张培震等，2008)，断层两盘同震形变有较大的垂直运动分量(Xu et al.，2009)。以往对汶川地震同震形变场的认识大多基于GPS水平观测(国家重大科学工程“中国地壳运动观测网络”项目组，2008;Shen et al.，2009;Wang et al.，2011)、InSAR视线向观测(单新建等，2009;Qu et al.，2010)以及断层附近局部区域水准观测(Wang et al.，2008)，这些观测对认识汶川地震发震断层形变性质有积极作用。由于汶川地震以逆冲性质为主，同震垂直形变场的获取对理解龙门山断裂带的三维形变特征和空间分布有重要意义。本文采用InSAR与GPS相结合方法，以汶川地震为例，利用InSAR获取的斜距向同震形变场，结合GPS水平同震观测，通过校正以及插值等处理，获取了汶川地震垂直连续形变场。

1 形变观测数据获取与收集

1.1 汶川地震InSAR斜距向同震形变场

利用InSAR技术，采用地震前后日本ALOS/PALSAR数据，提取了汶川地震7个条带的地表同震形变场(图1)，从东向西依次是471、472、473、474、475、476和477条带。可以看出，整个同震形变场覆盖了汶川地震震区400km×500km的范围。考虑到卫星为升轨观测，蓝色负值区域为远离卫星，黄色正值区域为靠近卫星，即断层上盘以向东运动为主，位移为－110～－120cm。断层下盘以向西运动为主，位移为120～130cm。

图1 汶川地震InSAR和GPS水平位移观测结果Fig.1 Horizontal displacement of Wenchuan earthquake measured by InSAR and GPS.

由于InSAR形变观测结果是视线向的，是水平(EW和SN向)和垂直向形变在视线向的投影，因此，InSAR观测到的只是一维地表形变量，也就很难从汶川同震形变场中识别和分析出垂直形变信息和特征。而汶川地震发震断层以逆冲为主，断层两盘近场形变应有较大的垂直运动分量，断层近场垂直形变对认识汶川地震发震断层形变特征有重要作用。

1.2 汶川地震GPS同震形变观测数据

Wang等(2011)收集了汶川地震较为全面的GPS同震观测数据，给出了38个连续观测、435+12个流动观测，以及33个三角GPS测量结果。此外，还给出了42个水准点观测结果。从图1可以看出，约有57%，共计295个GPS观测点位于InSAR观测区域内，这就保证了有较多的GPS观测值对InSAR条带进行校正，以及对GPS水平观测结果进行插值处理。此外，42个水准观测点位也全部位于InSAR观测范围内。

2 汶川地震垂直形变场研究思路

图2 汶川地震同震垂直形变场获取流程Fig.2 Flowchart process deriving the vertical coseismic deformation of the Wenchuan earthquake.

整个汶川地震同震垂直形变场获取流程和思路见图2。考虑到GPS水平观测精度高，且点位密度足够多，以及InSAR连续形变场对垂直形变敏感等优势，采用了GPS水平观测量与InSAR视线向形变联合解算的方法，求取垂直形变分量。总体处理方法按以下步骤:1)采用Biharmonic样条插值对GPS水平观测量——水平位移量大小和方向进行插值，然后再分解为EW和SN 2个方向的形变场;2)对InSAR视线向形变场进行重采样，降低其空间分辨率;3)InSAR与GPS观测是在不同的坐标系下进行的，有必要利用GPS观测点对InSAR观测值进行参考点校正，便于以后的三维形变场的合成与分解。考虑到InSAR条带解缠是分断层上、下盘分别进行的，因此，校正时也分上、下盘分别进行，校正完后再拼接成一个完整的汶川地震InSAR视线向形变场;4)采用三维形变场分解模型，结合以上获取的EW和SN向形变场，以及校正后的InSAR视线向形变场，获取汶川地震同震垂直形变场。

3 数据处理流程

3.1 InSAR视线向形变量的校正

由于InSAR观测到的形变场是相对每个条带的边缘点独立解缠的，而GPS观测结果是相对于欧亚板块，因此，有必要利用GPS观测点对InSAR观测值进行参考点校正，这样就把InSAR与GPS观测值放在同一参考系下，便于以后的三维形变场的合成与分解。在校正前，首先对InSAR同震视线向形变场进行重采样，降低其空间分辨率。考虑到InSAR观测结果的点密度极高，而GPS观测点位相对稀少，因此，为降低计算耗时，对每个InSAR条带进行了重采样，重采样后每个条带的网格点为512×1 024;其次，由于龙门山断裂带附近的形变梯度极大，造成InSAR观测严重失相关。极震区InSAR条带解缠是分断层上、下盘分别进行的，因此，用GPS对InSAR进行参考点校正时，也必须分上、下盘分别进行。在GPS点位选取过程中，如果该GPS点位的观测值＞2倍的残差，则不采用该点作为校正点。这样每个InSAR条带在断层上盘或者下盘仍能保持足够多的GPS观测点来求取精度较高的InSAR与GPS的偏移量。图3给出了校正前、后InSAR视线向形变场的对比图。

经过对比分析图3a和b可以看出，经过校正后的形变场条带之间的一致性得到了改善。如476号条带的北段，校正后形变值有所降低，与相邻的475号条带的形变值更为接近。图像左上角471号条带，校正前形变值偏大，校正后形变值与472号条带接近。

图3 校正前(a)与校正后(b)InSAR视线向形变场对比图Fig.3 Comparison of uncorrected(a)and corrected(b)InSAR LOS deformation.

图4 475号条带形变剖面校正前(a)与校正后(b)InSAR与GPS观测值对比图Fig.4 Profiling comparison of uncorrected(a)and corrected(b)InSAR-and GPS-measured deformations along track No.475.

图4给出了以475号条带为例，由南向北的InSAR和GPS观测形变校正前、后的剖面对比。475号条带范围共包含GPS观测点61个，其中上盘19个，下盘42个。首先将GPS观测到的EW，SN，UP 3个方向的分量投影到视线向，即求取GPS观测在视线向的形变量，然后与InSAR观测值进行比较。断层上盘GPS与InSAR的标准差为44.4，下盘为12.2。按＜2倍的标准差来挑选GPS点位，筛选后，上盘有18个GPS点位，下盘有40个GPS点位。上盘GPS与InSAR观测值得平均差值为9.93cm，而下盘GPS与InSAR观测值得平均差值为－11.49cm，这2个值也就是InSAR观测值需要校正的偏移常量。校正后，上盘GPS与InSAR观测值的平均差值为－0.005cm，标准差降为19.09，下盘GPS与InSAR观测值的平均差值为－0.004cm，标准差降为5.1。可以看出，校正后InSAR与GPS一致性得到了改善。

3.2 GPS水平向观测值插值方法及结果

考虑到GPS同震水平形变变化的特性，采用了Biharmonic样条插值对GPS水平向观测形变矢量(水平位移量大小和方向)进行插值。采用Biharmonic样条插值考虑到3个原因:1)GPS水平观测的矢量方向和大小，在断层上盘或下盘地表变化应为连续有序变化;2)对于形变来讲，发震断层为间断面，即断层两盘的GPS运动方向相反，而该方法可以较好地解决断层形变间断面问题，符合发震断层两盘的形变物理特征;3)GPS观测点位不规则且离散分布，而该方法对控制点分布和数量没有要求。

Biharmonic样条插值具有以下3个特点:1)Biharmonic样条曲面插值是一个最小曲率插值，这是因为该方法既可利用各点的数值又可以利用斜率产生插值曲面，由于曲面满足Biharmonic方程，因此插值结果局部平滑性好。2)插值产生的结果是以各数据点为中心的Green函数的线性组合，Biharmonic方程在不同维空间中的解就是不同维的Green函数。3)Green函数的个数少于数据点的个数，从而使得插值结果不会与不精确数据点相匹配(Sandwell et al.，1987;王亚涛等，2007;李元金等，2011;陆鹏，2011)。

对于观测到的N个离散观测点形变数据Pi，i=1，…，N，Biharmonic样条插值问题可以转化为方程组(1)的求解:

式(1)中:▽4为Biharmonic算子;φ(P)为二维Green函数，

P表示平面位置，ω(P)表示平面位置P处的形变值。满足式(1)的通解为

其中，系数aj满足线性方程组(3):

在方程组中，当i=j时，令该项为0，如ajφPi-Pj( )=0。对于公式(3)，当输入N个形变观测点Pi(i=1，…，N)的坐标Pixi，yi( )及形变值ω(Pi)时，可以得到N个形变观测点Pi组成的线性方程组。然后求解线性方程组，得系数ai，i=1，…，N。有了ai，就可以利用公式(2)求取待插值点的形变值。

依据GPS观测水平精度较高优势，利用Biharmonic样条插值，对GPS水平量大小和方向进行插值，获得高精度的水平形变大小和方向矢量场，然后再分解为EW和SN 2个方向的形变场。图5给出了插值后得到的GPS的EW向和SN向形变场图像。

图5给出的形变场图像，相比离散的GPS观测点能较好地反映出地表连续形变场的变化特征。图5a中，红黄色阶为正值，表示地壳向东运动，蓝紫色阶为负值，表示地壳向西运动。可以看出，龙门山发震断层的逆冲分量由南向北逐步递减，其主要逆冲形变量集中在映秀—连山坪、茶坪—北川县城—南坝区段。同时，在映秀镇—连山坪一带出现了东向位移极大区域(白色区域)，断层上盘东向位移量达到585cm，而断层下盘向西位移量大约为100cm。而在茶坪—北川县城—南坝区段，断层下盘向西位移量达到－280cm，断层上盘向东位移量为290cm。发震断层南段上盘向东运动占绝对优势，而发震断层中北段，则转换为发震断层上、下盘对冲位移为主。

在图5b中，红黄色阶为正值，表示向北运动，蓝紫色阶为负值，表示向南运动。可以看出，发震断层南段向北运动明显，有一定的右旋走滑分量。在桥头一带，北向位移分量达到185cm，在青川南侧，北向位移分量为125cm。

3.3 三维形变观测量的分解

InSAR获取的形变量并不能代表真实的地表形变，而是地表东西(EW)、南北(SN)和垂直(UP)3个方向形变分量在雷达视线LOS向上的投影，这也就给出了4者之间的关系。如图6所示，θ为卫星中心入射角，α为卫星飞行方向(Heading方向)与北向夹角(顺时针)，α-3π/2为方位视线方向(Azimuth Look Direction，简记为ALD)，即距离向与北向的夹角。

根据InSAR测量与3D形变量的几何关系，同时约定目标远离雷达时LOS向形变量dLOS为负(视线向拉伸)，靠近雷达时dLOS为正(视线向缩短)，则视线向形变量(dLOS)与地表真实三维形变量(dU，dN，dE)之间的关系为

由公式(4)可以看出，只要知道dLOS、dU、dN和dE中的3个量，就可以求解出第4个量。

4 汶川同震垂直形变场获取与分析

图5 由GPS观测值插值得到的EW向(a)和SN向(b)形变场Fig.5 East-westward(a)and north-southward(b)deformation after interpolation based on GPS measurements.白色区域的形变范围为:240～585cm

图6 InSAR视线向与三维分量的几何关系Fig.6 The geometry between InSAR line of sight observation and the geographical coordination.

图5给出的由GPS插值得到的EW向和SN向形变场dN和dE，以及图3b得到的校正后的InSAR视线向形变场dLOS，利用公式(4)，解算出汶川地震同震垂直形变场dU。图7给出了解算出的汶川地震同震垂直形变场，覆盖范围为400km×500km。

图7 解算出的汶川地震同震垂直形变场Fig.7 The derived Wenchuan earthquake coseismic vertical deformation.白色区域的垂直形变范围为:320～550cm

由图7可以看出，汶川地震同震垂直形变场与图1的斜距向同震形变场，以及图5的EW向和SN向形变场有不同的形变特征。从分布上来看，有以下特征:1)较大的垂直形变并不是沿整个发震断层两侧都有，而是分段分布，主要集中分布在映秀－北川发震断裂的南段、中段和北端;2)断层两侧的垂直形变衰减较快，垂直形变变化量＞30cm的区域主要分布在断裂带附近，距发震断层的距离两盘之和不超过50km，这可能与断层倾角接近直立有关;3)在发震断层南段，断层两侧垂直形变极不对称，主要分布在汶川县城至都江堰之间30km×40km的断层上盘区域内，最大抬升在映秀镇至连山坪一带，抬升达到5.5m。而在断层下盘，沉降量较低，只有约－0.6cm;在中段，垂直形变表现为较强的反对称性。在茶坪—北川—南坝一带，断层上盘最大抬升区域在茶坪东侧，抬升量为255cm。北川县城北部的白泥抬升量为185cm，断层下盘永庆沉降量为－215cm。断层上盘桥头西侧的抬升量为215cm，断层下盘响岩的沉降量为－190cm。在发震断层北端，垂直形变呈对称分布，青川一带有较大面积的抬升，最大抬升量为120cm。Shen等(2009)采用GPS和InSAR数据对汶川同震形变场和断层滑动量分布进行了研究，结果表明，在映秀镇、北川和南坝3个地区断层滑动量出现极大值，对应地造成地表震害也最为严重。这一结果与本文垂直形变场分布基本一致。

474号条带范围内包含了88个GPS、水准测量点，去除＞3倍标准差的观测值，还剩下80个垂直观测值，可以用来与解算出的垂直形变值进行对比。图8所示为解算出的垂直形变与GPS、水准垂直观测形变的剖面对比。可以看出，InSAR解算出的垂直形变与GPS、水准观测在总体趋势上是一致的。较大的垂直形变主要在断层两侧附近区域。在断层下盘，远离断层区域，两者有4cm左右的差异。在断层上盘远离断层区域，两者有5cm左右的差异，但表现均较为稳定。而在断层附近，因形变起伏较大，两者差异较大。断层附近水准MB23标识点观测到的最大垂直位移为124.99cm，而从InSAR解算出的垂直位移为101.96cm，两者相差23.03cm。另外，在31.7°N附近，InSAR解算出的垂直形变为－28cm，而此处水准MB11标志点的垂直位移为3.16cm，相差31.9cm。断层附近两者相差较大的主要原因是由于在断层附近InSAR失相关严重，造成不能获取或获取精度不够。

图8 解算垂直形变与GPS、水准垂直观测形变的剖面对比图Fig.8 Comparison among the derived vertical deformation，GPS measurements，leveling observation in profiles.

5 结论

InSAR与GPS结合，能够在观测优势方向和空间域上互相补充，不但可提高观测精度，而且可获取丰富的地表形变场信息。InSAR获取的是一维视线向形变量，是水平形变和垂直形变在视线方向的投影。但当有较多的GPS水准观测数据时，两者可以分解出较为可靠的垂直形变信息。研究得到了以下认识:

(1)采用Biharmonic样条插值先对GPS水平形变的大小和方向进行插值，更符合地表形变变化有序等物理特征。考虑到断层上盘或者下盘的形变矢量的大小和方向变化呈有序稳定状态，先采用Biharmonic样条插值对地表水平形变矢量的大小和方向进行插值，然后再进行EW和SN向分解，较好地克服了其他插值方法在断层形变间断面附近出现紊乱现象，比直接对EW和SN向形变进行插值更接近真实的形变特征。

(2)InSAR与GPS观测的坐标框架不同，两者联合解算，必须进行坐标参考点校正。InSAR观测到的形变参考点是相对图像自身的，是近场的，而GPS观测结果是相对于欧亚板块，是远场的。因此，两者观测值之间存在一个偏差常量。在利用两者数据进行三维形变场合成或分解时，必须利用远场的GPS观测值对InSAR参考点进行校正。

(3)解算出的汶川地震同震垂直形变场既符合目前对发震断层形变特征的研究成果，又取得了一些新的认识。断层两侧垂直形变衰减较快，横跨断裂带形变量＞30cm的宽度不超过50km;沿断层垂直形变高值区分布不均匀，主要集中分布在发震断裂的南段、中段和北端。这3段各具特色，南段断层两侧垂直形变极不对称，主要以上盘剧烈抬升为主，最大抬升区域在映秀镇至连山坪一带，抬升量达到5.5m。中段表现为较强的反称性，断层一侧抬升另一侧沉降。该段上盘最大抬升区域在茶坪东侧，抬升量为255cm，下盘最大沉降量在永庆，沉降量为－215cm。北端垂直形变呈对称分布，表现为“狮子尾”图像特征，挂在发震断层最北端，最大抬升量为120cm，分布在青川北侧。

邓起东，陈社发，赵小麟.1994.龙门山及其邻区的构造和地震活动及动力学［J］.地震地质，16(4):389—403.

DENG Qi-dong，CHEN She-fa，ZHAO Xiao-lin.1994.Tectonics，seismicity and dynamics of Longmenshan Mountains and its adjacent regions［J］.Seismology and Geology，16(4):389—403(in Chinese).

国家重大科学工程“中国地壳运动观测网络”项目组.2008.GPS测定的2008年汶川MS8.0地震的同震位移场［J］.中国科学(D辑)，38(10):1195—1206.

Working Group of Crustal Movement Observation Network of China(CMONOC).2008.Coseismic displacement field of the MS8.0 Wenchuan earthquake in 2008 measured by GPS surveys［J］.Science in China(Ser D)，38(10):1195—1206(in Chinese).

李元金，罗立民，张鹏程，等.2011.基于校正靶特征与Biharmonic样条插值的XRII图像扭曲校正［J］.东南大学学报(自然科学版)，41(6):1213—1218.

LI Yuan-jin，LUO Li-min，ZHANG Peng-cheng，et al.2011.Distortion correction of XRII image based on calibration grid characteristic and Biharmonic interpolation［J］.Journal of Southeast University(Natural Science Edition)，41(6):1213—1218(in Chinese).

陆鹏.2011.基于Biharmonic样条的离散数据三维地形生成方法［J］.广西水利水电，2011(6):64—73.

LU Peng.2011.3D terrain generation out of discrete data by Biharmonic spline method［J］.Gunagxi Water Resource＆ Hydropower Engineering，2011(6):64—73(in Chinese).

单新建，屈春燕，宋小刚，等.2009.汶川MS8.0地震InSAR同震形变场观测与研究［J］.地球物理学报，52(2):496—504.

SHAN Xin-jian，QU Chun-yan，SONG Xiao-gang，et al.2009.Coseismic surface deformation caused by the Wenchuan MS8.0 earthquake from InSAR data analysis［J］.Chinese J Geophys，52(2):496—504(in Chinese).

汪驰升，单新建，张国宏，等.2009.基于ASAR升降轨数据解算于田MS7.3地震3D同震形变场［J］.地震，29(B10):105—112.

WANG Chi-sheng，SHAN Xin-jian，ZHANG Guo-hong，et al.2009.3D coseismic deformation field of the Yutian MS7.3 earthquake calculated from ASAR ascending and descending data［J］.Earthquake，29(B10):105—112(in Chinese).

王亚涛，董兰芳，倪奎.2007.基于Biharmonic样条插值的图像渐变算法及实现［J］.中国图像图形学报，12(12):2189—2194.

WANG Ya-tao，DONG Lan-fang，NI Kui.2007.Image morphing algorithm based on Biharmonic spline interpolation and its implementation［J］.Journal of Image and Graphics，12(12):2189—2194(in Chinese).

张培震，徐锡伟，闻学泽，等.2008.2008年汶川8.0级地震发震断裂的滑动速率、复发周期和构造成因［J］.地球物理学报，51(4):1066—1073.

ZHANG Pei-zhen，XU Xi-wei，WEN Xue-ze，et al.2008.Slip rates and recurrence intervals of the Longmenshan active fault zone and tectonic implications for the mechanism of the May 12 Wenchuan earthquake，2008，Sichuan，China［J］.Chinese J Geophys，51(4):1066—1073(in Chinese).

Hu Jun，Li Zhi-wei，Ding Xiao-li，et al.2008.Two-dimensional co-seismic surface displacements field of the Chi-Chi earthquake inferred from SAR image matching［J］.Sensors，8:6484—6495.

Qu Chun-yan，Shan Xin-jian，Zhang Guo-hong，et al.2010.Coseismic displacement field of the Wenchuan MS8.0 earthquake in 2008 derived using differential radar interferometry［J］.Journal of Applied Remote Sensing，4(1):043516(19 March 2010).doi:10.1117/1.3386043.

Samsonov S，Tiampo K.2006.Analytical optimization of a DInSAR and GPS dataset for derivation of three-dimensional surface motion［J］.IEEE Geosci Remote Sens Lett，3(1):107—111.

Sandwell D T.1987.Biharmonic spline interpolation of GEOS－3 and SEASAT altimeter data［J］.Geophysical Research Letters，14(2):139—142.

Shen Zheng-kang，Sun Jian-bao，Zhang Pei-zhen，et al.2009.Slip maxima at fault junctions and rupturing of barriers during the 2008 Wenchuan earthquake［J］.Nature Geoscience，2:718—724，Published online:27 September 2009.doi:10.1038/ngeo636.

Sverrir G，Freysteinn S.2002.Three-dimensional surface motion maps estimated from combined interferometric synthetic aperture radar and GPS data［J］.Journal of Geophysical Research，B10(107):1029.

Wang Qi，Qiao Xue-jun，Lan Qi-gui，et al.2011.Rupture of deep faults in the 2008 Wenchuan earthquake and uplift of the Longmen Shan［J］.Nature Geoscience，4:634—640.doi:10.1038/nego1210.

Wang Qing-liang，Cui Du-xin，Zhang Xi，et al.2009.Coseismic vertical deformation of the MS8.0 Wenchuan earthquake from repeated levelings and its constraint on listric fault geometry［J］.Earthquake Science，22:595—602.

Wright T J，Parsons B E，Lu Z.2004.Toward mapping surface deformation in three dimensions using InSAR ［J］.Geophys Res Lett，31:L01607.doi:10.1029/2003GL018827.

Xu Xi-wei，Wen Xue-ze，Yu Gui-hua，et al.2009.Coseismic reverse-and oblique-slip surface faulting generated by the 2008 MW7.9 Wenchuan earthquake，China［J］.Geology，37(6):515—518.

Yuri F，Mark S.2001.The complete(3－D)surface displacement field in the epicentral area of the 1999 MW7.1 Hector Mine earthquake，California，from space geodetic observation ［J］.Geophysical Research Letters，28(16):3063—3066.