断裂区带变形分析方法及应用

第３２卷第３期　西北地震学报　ＶｏＩ．３２　Ｎｏ．３　２０１０年９月　Ｎ０ＲＴＨＷＥＳＴＥＲＮ　ＳＥＩＳＭＯＬＯＧＩＣＡＬ　ＪＯＵＲＮＡＬ　Ｓｅｐｔ．，２０１０　断裂区带变形分析方法及应用　杨　博，周　伟，陈阜超，韩月萍　（中国地震局第一监测中心，天津３００１８０）　摘要：应用ＧＰＳ观测技术专门研究活动断裂区带的地壳形变及其服务于地震预测的数据处理方　法目前还不够完善与充分。本文针对这一问题并结合实际的需求，从场的角度提出了适合于断裂　区带活动特征分析的数据处理方法，并以川滇菱形块体东边界带为例，依据１９９９—２００４年的ＧＰＳ　资料对其进行了简要分析。结果表明：①经过这种处理与描述可使人们从多角度更清楚全面地了　解断裂及两侧的相对活动，及在空间上的演化特征；②走滑运动等值线显示川滇菱形块体东边界两　侧相对活动量最大达１５　ｍｍ／ａ，分布宽度约４００　ｋｍ，然而活动量的２／３只分布在断层两侧近１００　ｋｍ的范围内，清晰而定量地显现了高剪切应变的积累与存在的空间；③走滑运动梯度显示川滇形　块体东边界带南段变形大于北段；④张压性运动结果显示断裂两侧没有明显的差异变化。　关键词：ＧＰＳ：数据处理方法；梯度变化；活动断裂区带；ＪＩＩ滇块体东边界带　中图分类号：Ｐ２２７　文献标识码：Ａ　文章编号：１０００—０８４４（２０１０）０３～０２１５—０５　Ａｎａｌｙｓｉｓ　Ｍｅｔｈｏｄ　ｏｆ　ｔｈｅ　Ｃｒｕｓｔａｌ　Ｄｅｆｏｒｍａｔｉｏｎ　Ｏｉｌ　Ｆａｕｌｔ　Ｚｏｎｅ　ａｎｄ　Ｉｔｓ　Ａｐｐｌｉｃａｔｉｏｎ　ＹＡＮＧ　Ｂｏ，ＺＨＯＵ　Ｗｅｉ，ＣＨＥＮ　Ｆｕ—ｃｈａｏ，ＨＡＮ　Ｙｕｅ—ｐｉｎｇ　（Ｆｉｒｓｔ　Ｃｒｕｓｔａｌ　Ｍｏｎｉｔｏｒｉｎｇ　ａｎｄ　Ａｐｐｌｉｃａｔｉｏｎ　Ｃｅｎｔｅｒ，ＣＥＡ，Ｔｉａｎｊｉｎ　３００１８０，Ｃｈｉｎａ）　Ａｂｓｔｒａｃｔ：Ｔｈｅ　ｒｅｓｅａｒｃｈ　ｏｎ　ｃｒｕｓｔａｌ　ｄｅｆｏｒｍａｔｉｏｎ　ｅｓｐｅｃｉａｌｌｙ　ｉｎ　ａｃｔｉｖｅ　ｆａｕｌｔ　ｚｏｎｅ　ｕｓｉｎｇ　ＧＰＳ　ｔｅｃｈｎｉｑｕｅ　ａｎｄ　ｔｈｅ　ｒｅｌａｔｅｄ　ｄａｔａ　ｐｒｏｃｅｓｓｉｎｇ　ｍｅｔｈｏｄ　ｆｏｒ　ｅａ￣ｈｑｕａｋｅ　ｐｒｅｄｉｃｔｉｏｎ　ａｒｅ　ｎｏｗ　ｎｏｔ　ｐｅｒｆｅｃｔ　ａｎｄ　ｅｎｏｕｇｈ．Ａｉｍｅｄ　ａｔ　ｔｈｉｓ　ｐｒｏｂｌｅｍ　ａｎｄ　ｆｏｒ　ａｃｔｕａｌ　ｗｏｒｋ　ｒｅｑｕｉｒｅｍｅｎｔ，ｔｈｅ　ｄａｔａ　ｐｒｏｃｅｓｓｉｎｇ　ｍｅｔｈｏｄ　ｆｏｒ　ｔｈｅ　ｃｈａｒａ—ｅｔｅｒｉｓｔｉｃ　ａｎａｌｙｓｉｓ　ｏｆ　ｃｒｕｓｔａｌ　ｄｅｆｏｒｍａｔｉｏｎ　ｉｎ　ｆａｕｌｔ　ｚｏｎｅ　ｉｓ　ｓｕｇｇｅｓｔｅｄ　ｆｒｏｍ　ｆｉｅｌｄ　ｖｉｅｗ　ｉｎ　ｔｈｅ　ｐａｐｅｒ．Ｂａｓｅｄ　ｏｎ　ｔｈｅ　ＧＰＳ　ｄａｔａ　ｆｒｏｍ　Ｓｉｃｈｕａｎ—Ｙｕｎｎａｎ　ｒｅｇｉｏｎ　ｉｎ　１　９９９—２００４，ｔｈｅ　ｃｒｕｓｔａｌ　ｍｏｖｅｍｅｎｔ　ｃｈａｒａｃｔｅｒｓ　ｉｎ　ｅａｓｔｅｒｎ　ｂｏｕｎｄａｒｙ　ｆａｕｌｔ　ｚｏｎｅ　ｏｆ　Ｓｉｃｈｕａｎ—Ｙｕｎｎａｎ　ｂｌｏｃｋ　ａｒｅ　ａｎａｌｙｚｅｄ　ｂｙ　ｕｓｉｎｇ　ｔｈｉｓ　ｍｅｔｈｏｄ．Ｔｈｅ　ｒｅｓｕｌｔｓ　ｓｈｏｗ　ａｓ　ｆｏｌｌｏｗｓ：（１）Ｔｈｅ　ｒｅｌａ—　ｔｉｖｅ　ｍｏｖｅｍｅｎｔ　ａｎｄ　ｓｐａｔｉａｌ　ｖａｒｉａｔｉｏｎ　ｉｎ　ｂｏｔｈ　ｓｉｄｅｓ　ｏｆ　ｔｈｅ　ｆａｕｌｔ　ｚｏｎｅ　ｃａｎ　ｂｅ　ｓｈｏｗｎ　ｃｌｅａｒｌｙ　ｉｎ　ｄｉｆｆｅｒｅｎｔ　ｖｉｅｗ　ａｎｇｌｅｓ　ｂｙ　ｔｈｉｓ　ｄａｔａ　ｐｒｏｃｅｓｓｉｎｇ　ｍｅｔｈｏｄ．（２）Ｔｈｅ　ｄｅｆｏｒｍａｔｉｏｎ　ｗｉｄｔｈ　ｏｎ　ｂｏｔｈ　ｓｉｄｅｓ　ｏｆ　ｔｈｅ　ｆａｕｌｔ　ｚｏｎｅ　ｉｓ　ａｂｏｕｔ　４００　ｋｍ　ｗｉｔｈ　ｔｈｅ　ｍａｘｉｍｕｍ　ｍｏｖｅｍｅｎｔ　ｖａｌｕｅ　ｏｆ　１５　ｍｍ／ａ．Ｈｏｗｅｖｅｒ，２／３　ｏｆ　ｔｈｅ　ｍｏｖｅｍｅｎｔ　ｄｉｓｔｒｉｂｕｔｅ　ｉｎ　ａ　ｂｅｌｔ　ａｒｅａ　ａｂｏｕｔ　１００　ｋｍ　ｗｉｄｔｈ　ａｌｏｎｇ　ｔｈｅ　ｆａｕｌｔｓ．Ｔｈｉｓ　ｃｌｅａｒｌｙ　ｓｈｏｗｓ　ｔｈｅ　ｄｉｓｔｒｉｂｕｔｉｏｎ　ｏｆ　ｓｌｉｐ　ｍｏｖｅｍｅｎｔ　ａｎｄ　ｓｈｅａｒｉｎｇ　ｓｔｒａｉｎ　ａｃｃｕｍｕｌａｔｉｏｎ　ｆｒｏｍ　ｔｈｅ　ｍｏｔｉｏｎ　ｉｓｏｇｒａｍ．（３）Ｔｈｅ　ｄｅｆｏｒｍａｔｉｏｎ　ｏｎ　ｓｏｕｔｈｅｒｎ　ｐａｒｔ　ｉｓ　ｍｏｒｅ　ｔｈａｎ　ｔｈａｔ　ｏｎ　ｎｏｒｔｈｅｒｎ　ｐａｒｔ　ｏｆ　Ｅａｓｔｅｒｎ　ｂｏｕｎｄａｒｙ　ｚｏｎｅ　ｏｆ　Ｓｉｃｈｕａｎ——Ｙｕｎｎａｎ　Ｂｌｏｃｋ　ｓｈｏｗｅｄ　ｂｙ　ｔｈｅ　ｇｒａｄｉｅｎｔ　ｖａｒｉａ－－　ｔｉｏｎ　ｍａｐ　ｏｆ　ｓｌｉｐ　ｍｏｖｅｍｅｎｔ．（４）Ｉｎ　ｐｅｒｐｅｎｄｉｃｕｌａｒ　ｄｉｒｅｃｔｉｏｎ　ｔｏ　ｔｈｅ　ｆａｕｌｔ　ｔｈｅｒｅ　ｉｓ　ｎｏｔ　ｃｌｅａｒｌｙ　ｄｉｆｆｅｒｅｎｔ　ｄｅ—　ｏｆｒｍａｔｉｏｎ　ｂｅｔｗｅｅＴ１　ｔｗｏ　ｓｉｄｅｓ　０ｆ　ｔｈｅ　ｚｏｎｅ　Ｋｅｙ　ｗｏｒｄｓ：ＧＰＳ；Ｄａｔａ　ｐｒｏｃｅｓｓｉｎｇ　ｍｅｔｈｏｄ；Ｇｒａｄｉｅｎｔ　ｖａｒｉａｔｉｏｎ；Ａｃｔｉｖｅ　ｆａｕｌｔ　ｚｏｎｅ；Ｅａｓｔｅｒｎ　ｂｏｕｎｄａｒｙ　ｚｏｎｅ　ｏｆ　Ｓｉｃｈｕａｎ—－Ｙｕｎｎａｎ　ｂｌｏｃｋ　０　引言　ＧＰＳ技术的成熟与应用普及使地壳运动与地球　收稿日期：２００９－０４—２３　基金项目：“十一五”国家科技支撑课题（２００６ＢＡＣ０３～０１—０１）　作者简介：杨博（１９８５一），男（汉族），实习生，主要从事ＧＰＳ应用技术研究　２１６　西北地震学报　可　动力学的研究生机勃勃。在我国自１９９８年国家重　大科学工程中国地壳运动观测网络建成与运行以　来，产生了大量的观测数据，相应的应用研究、方法　研究结果不断涌现。主要有时间序列方面的研　究¨　，形变场方面的研究ｌ３　Ｊ、形变场与地震关系　方面的研究　Ｊ、负位错反演与地震预测方面的研　利用ＧＰＳ数据处理软件所给出的运动等结果　往往是ＩＴＲＦ参考框架下的结果，这样的结果在分　析全球各板块之间的差异运动时是很恰当的，但在　站　的　分析局部区域内部的差异运动时则不能一目了然，　所以在分析中国大陆地壳运动时又出现了以欧亚板　差　＿＝＝．　方　块为参考基准的数据处理方法、以中国大陆整体无　１●●●Ｊ　ｊＥ　建　立　第３２卷　误　究　、数值模拟方面的研究　等。然而，针对以断　裂区带为对象的研究从空间的角度来观察则不多　见。尽管一些学者开展了对断裂带的活动习性的分　段研究¨　及跨断层形变剖面研究Ｌ１¨等，但随着认　识的深入与实际的需求显示出不够完善，主要是因　旋转为参考的基准以及区域整体无旋转基准等数据　处理方法　”　，使人们分析起来显得清晰、直观、方　一　＝　：　便，突现了彼此之间的差异变化。本文以这些工作　．　．一　为基础，但没有按照文献［１２］、［１３］所给出的数据　Ａ　处理流程进行处理，而是直接利用观测运动求解研　为在空间上的变化信息或没有得到全面而又充分地　体现，或不够一目了然，或简化的成分较重。本文针　对断裂区带活动的特点和ＧＰＳ数据资源，试图在前　人工作的基础上争取进一步完善形变信息表述模式　及实现方法。　究区的整体运动并剔除就可满足分析的要求。这样　可使比较复杂的数据处理变得相对简单些。　首先要求解研究区的整体运动并剔除，由此来　获得研究区的相对差异运动。根据欧拉模型　１　数据处理与实现方法　一ｓｉｎ（　）ｓｉｎ（Ａ）　一（　。　ＣＯＳ（Ａ）　［：二］＝ｒ［一ｓｉ　；　（Ａ）　一　Ａ）ｃ。　三ｉ］一［　：］　ｃ２　为了简明起见，对所有测站的误差方程用矩阵　表示：　Ｖ＝ＡＸ＋Ｌ　（３）　析的一个基础。由于我们的研究对象是断裂带（断　层）的活动习性及其在空间上的演化，而断裂带的　走向是任意方向的，故需要进行投影或转换计算。　在此基础上依据最小二乘法即可组成法方程　ＮＸ＋Ｗ＝０　（４）　根据断裂带两侧块体主要表现为走滑活动、正断层　活动或逆断层活动，可通过旋转将直角坐标某一轴　平行于断裂带的走向，另一轴垂直于断裂带的走向，　对其求解可获得研究区的欧拉矢量　Ｘ＝一Ｎ　Ｗ　（５）　再利用式（１）即可解算任一测站的　、　，进而可获　得相应的差异运动：　＇ｌｆ／ｅ　ｅ　对Ｙｅ　、Ｙｎ　实施计算。首先，结合图１及坐标旋转公　式获得平行及垂直于断裂带的运动分量：　．ｖ　断裂带走向　ｒ　］＝［ｖｅ　—－—Ｌ　］　（６）　／　一式中　、　为研究区内任一点（Ａ，　）的欧拉东向运　动与北向运动；　、Ｖ　为相应的东向观测运动与北向　观测运动；　、　为相应的东向差异运动与北向差异　ｖ　运动；ｒ为地球平均半径。　通过式（６）我们可获得研究区相对差异运动，　以此作为先验速率重新计算以获取差异运动结果　（　）及相应的误差（ｍ　、ｍ　和Ｐ），并在此基础上　图１　坐标系及坐标变换示意图　Ｆｉｇ．１　Ｃｏｏｒｄｉｎａｔｅｓ　ｓｙｓｔｅｍ　ａｎｄ　ｉｔｓ　ｔｒａｎｓｆｏｌ￣ｎａｔｉｏｎ　绘制相应的图件。然而这一结果与图件仅是我们分　ｍ　∞　＝　第３期　十一　杨博等：断裂区带变形分析方法及应用　２１７　、＞　Ｃ　』‘　　＿　ｅ　，　２　，　ｍ　Ｓ　Ｏ　（、，７）　，　３。等值线的绘制可以用ｓｕｒｆｅｒ软件等，在此不再赘　述。除此之外，考虑到信息的全面性，进一步计算绘　Ｌ　ｎｎ　Ｖｅｅｓｉｎ　＋　ＭＣＯＳ　ＯＬ　２　Ｐ　ｍ　制相应的梯度图。根据梯度的定义　式中　和　：　分别为平行与垂直于断裂带的运动分　＋　量；　为旋转角，顺时针为正。　Ｓ　２　ｎ　ｇｒａｄ　＝　ｉ＋　（１５）　对　和　：　其进行精度评定，以考察其信息的显　著性。根据式（７）及误差传播定律，有　２　ｎ　　—２　Ｃ　Ｏ　Ｓ　ｍ　＝±￣／ｓｉｎ　・ｍ：＋ｃｏｓ　ｏ／・ｍ：＋２ｃｏｓＳ　ｎ　　Ｏ／・ｓｉｎ　・ｍ　（９）　式中ｍ　、ｍ　（及Ｐ）可从ＧＰＳ前期数据处理中所给出　）　的运动结果中获取，协方差ｍ　虽未直接给出，但可　由下式计算得到：　／＇／ｂ　＝ｍ　。ｍ　・Ｐ　（１０）　在实际分析工作中为了方便与直观，常常将形　变量、误差及其它与分析有关的要素（如断裂带、甚　至地形地貌等）绘制成图件，为统一坐标系，需要对　上述转换的两种运动量及其相应的误差换算到原坐　标系下。由图１可知，平行于断层的运动量其在Ｖｅ　ｏ　：　坐标系下的坐标为（　：　、０），垂直于断层运动量　的坐标则为（０、　：　）。顾及式（７），则平行于断层运　动量在原坐标系下的分量为　ＩＶｅｌ　ＶｅｅＣＯＳ（一　ｆ１１１　＝　ｓｉｎ（一　）　其相应误差变换的计算式为　』ｍｅ　ｍｅｅｅＯＳ（一　（１２）　【ｍ　ｌ＝ｍ　ｓｉｎ（一　）　同理，垂直于断裂带的运动量及其误差变换的　计算式为　ＩＶｅ２　一　ｍ　ｎ（一　ｆ　１３）　【　２＝ＶｎｎＣＯＳ（一Ｏ／）　Ｊｍ以　一ｒｎｎ．ｓｉｎ（一　’　（１４）　【Ｌｍ　＝ｍ：ｍ　＝ｍ　　ＣＯＳ（～Ｏ（～　）Ｌ）　考虑到图件的美观、简洁与信息的全面，在绘图　时误差量的正负应取决于位移（运动）的矢量方向，　即误差分解量的正负随位移（运动）分解量的正负　即可。这样不仅减少了图件中的不必要因素，还可　直接地观察位移（运动）的显著性（用不同的颜色或　深浅度表示，图２、３）。通常绘图时虽然给出矢量的　比例长度，但定量观察时仍觉得不十分方便。由于　这里已经是同方向的矢量（或平行、或垂直），所以　还可在同一类图里绘上描述大小的等值线，见图２、　ｌｇｒａｄｕ（　）Ｊ＝√（　）‘＋（　）‘　（１６）　可分别计算平行和垂直断裂带空间任一点的梯度　值。其中解析式／／，（　，Ｙ）的获取可参照文献［１２］所　给出的方法，也可根据其它数值逼近理论获取。得　到梯度值后可绘制等值线图。　２川滇东边缘断裂带及周围活动简析　川滇地区是我国地震活动最为频繁和强烈的地　区之一，中国地壳运动网络在此布设了较密集的　ＧＰＳ流动测站。现以１９９９—２００４年的观测资料为　例进行分析。由于计算是以ＩＴＲＦ为参照的运动结　果，不利于直接分析与观察，故根据式（１—６）采用　区域无旋转方法获得研究区的相对差异运动结果，　如图２所示。　图２川滇地区１９９９—２００４年相对于区域　无旋转基准的运动结果　Ｆｉｇ．２　Ｔｈｅ　ｃｒｕｓｔａｌ　ｍｏｖｅｍｅｎｔ　ｏｎ　Ｓｉｅｈｕａｎ—Ｙｕｎｎａｎ　ｒｅｇｉｏｎ　ｒｅｌａｔｉｖｅ　ｔｏ　ｎｏ・・ｎｅｔ－－ｒｏｔａｔｉｏｎ　ｒｅｆｅｒｅｎｃｅ　ｆｒａｍｅ　ｉｎ　ｌ９９９—２００４．　由图２可知，该地区水平运动在空间上的变化　是很有序的，最大形变区带分布在川滇菱形块体的　２１８　西北地震学报　第３２卷　东边界上，其它地域的差异活动较弱，说明东边界是　该区近代地壳构造活动最活跃部位。若观察东边界　东西两侧的相对运动，则该断裂及附近地域的差异　兼有较大的剪切变形的成分。进一步观察还可知该　地区的运动中既含有显著的变形分量，也含有刚性　分量。边界带西侧的点位运动方向恰好与青藏高原　物质进入该区的流动方向一致，向南运动在丽江处　性运动均在１０　ｍｍ／ａ以内，活动方式为左旋，向两　侧延伸这种差异运动在扩大，最大可达１５　ｍｍ／ａ。　它反映了该边界的活动并不纯粹具有走滑特性，还　９Ｒ。　１Ｏ０。　１ｎ，。　１　０４。Ｆ　部分向南东运动，部分继续向南，到达龙陵附近遇到　阻隔转向西。这是该区运动的基本特征。　９８。　１００　１０２　１０４。　（ａ）相对运动及其等值线　（ｂ）变化梯度　图３　平行于断裂带的相对运动及其等值线和变化梯度　Ｆｉｇ．３　Ｔｈｅ　ｒｅｌａｔｉｖｅ　ｍｏｖｅｍｅｎｔ　ｔｏ　ｂｅ　ｐａｒａｌｌｅｌ　ｔｏ　ｆａｕｌｔ，ａｎｄ　ｉｔｓ　ｉｓｏｇｒａｍ　ａｎｄ　ｃｈａｎｇｅ　ｇｒａｄｉｅｎｔ　为进一步了解形变变化在空间上与断层活动的　关系，利用上述有关公式给出了相关的结果，见图　３、４（梯度单位：ｒａｍ／１０　ｋｍ・ｙ　）。　图３（ａ）展现了平行于川滇块体东界断裂带的　２）。仔细观察可知，这种阻挡在南端也并非均等，　表现东西两端大中间小的格局，这可能与红河断裂　带的活动有关。　从地震的孕育角度看，这种在断裂带附近南北　运动（即走滑性运动）与等值线；图３（ｂ）为变化梯　度。可看出，沿鲜水河、则木河、小江断裂形成了一　向背离运动是断层具有走滑运动的体现，但由断裂　带附近区域运动等值线及梯度变化可知这并非只是　断层上的孤立运动，在断裂带及附近地域同时也有　剪切应变的存在与积累，具备地震孕发的条件。一　个正负等值线的分界带，等值线的展布方向基本上　与该带平行且形变梯度大，表明该带是东西两块体　相对运动的分界带。同时还显示该带的形变活动存　在于以断层为中心的具有一定宽度的空间上。东侧　向东形变的变化量逐渐变小（图３（ｂ）），至昭通时　基本上完成了较快的形变变化，再向东形变量的增　加已经很缓慢了。这说明昭通东侧介质刚性特征相　般地高形变、高梯度的地域在中短期内往往是中强　地震活动的主要场所，而在高形变、高梯度的背景下　存在低形变、低梯度的地域则有可能是强震活动的　场所，所以断裂带的中北段应更值得重视。　图４（ａ）为垂直于川滇块体东界断裂带的相对　运动（即张压性运动）及其等值线；图４（ｂ）是其运　对突出些，其西侧除断层走滑活动外其余属形变活　动。该带西侧的较快形变变化部分主要位于该带附　近，其它空间的形变变化相对宽缓，形变最大值位于　动变化的梯度。据此我们很难辨别垂直于断裂带的　运动与断层有什么关联性，因为不但运动的大小在　炉霍至丽江一线。尽管东西两块最大运动差达ｌ５　ｍｍ／ａ，但其中近１０　ｍｒｎ／ａ的变化分布在断裂带及其　附近大约１度的空间跨距上。西侧块体向南的运动　断裂带附近大体相当，而且其梯度变化也较缓慢。　确切地说，运动大小的空间分布与断裂带的活动没　有直接关系。在研究区北部表现为“西大东小”的　变化特征，变化最大的部位在离开断裂带的西侧块　体上（图４（ｂ）），表明该局部区域是北部挤压变形　向南变得越来越小，至其南端时已变成北向运动，这　表明物质的逃逸在南端受到了阻挡而改道向西（图　第３期　杨博等：断裂区带变形分析方法及应用　２１９　最大的区域（炉霍以北）。向南运动的大小与方向　也在发生变化，至龙陵及其西部地区时反向运动达　到最大。但变化最大的部位并不在此。在炉霍一丽　江一楚雄一玉溪一线存在一个较明显的北西向运动　变化梯度条带，同时也是运动性质差异变化的条带。　然而该形变条带并非是断裂构造活动带（至少目前　还未发现）。这可能是由该区不同断层共同控制造　成了这种复杂的形变形态。事实上，近几十年该条　带上已多次发生强烈地震，或许与此形变背景有某　种联系。此外，从表面上看研究区南部似乎存在着　相对拉张运动，参考图２和图３（ａ），看出是由于来　自北部物质的逃逸与推挤在南端受到阻挡而形成的　东西向挤出所致。　（ａ）相对运动及其等值线　（ｂ）变化梯度　图４垂直于断裂带的相对运动及其等值线和变化梯度　Ｆｉｇ．４　Ｔｈｅ　ｒｅｌａｔｉｖｅ　ｍｏｖｅｍｅｎｔ　ｔｏ　ｂｅ　ｐｅｒｐｅｎｄｉｃｕｌａｒ　ｔｏ　ｆａｕｌｔ　ａｎｄ　ｉｔｓ　ｉｓｏｇｒａｍ　ａｎｄ　ｃｈａｎｇｅ　ｇｒａｄｉｅｎｔ　３结语　综上所述，获得研究区相对差异运动（图２）可　带构造活动为主对川滇地区形变给出了分析，但还　较为粗浅，意在方法的应用。尽管如此仍可使我们　对断裂带构造活动及在空间上的关系有了更加清晰　的认知。可将本方法用于不同时段运动、形变等进　使我们对该区运动的概貌有较好的了解。然而在实　际工作中除了研究运动场的整体变化外，还常常要　结合断裂带（层）分析其走滑、张压活动的特征与习　性，及其相应的形变在空间上变化分布，并由此进一　步判定活动是否属于正常还是进入运动亏损或闭锁　状态，以及在此基础上确定相应的空间位置等。本　行对比演化分析上，因为当断裂带（层）活动发生变　化时，并不一定都体现在断裂带（层）的空问某部位　上，还可能体现在断裂带（层）的某时段活动上。　致谢：本项工作是在杨国华研究员指导下完成　的，并在本文的形成过程中又提出了详细的建议与　修改意见，在此谨表感谢。　文所给的方法有利于这样的分析。应该指出，在变　形机制的解释上不应离开研究区的整体运动（如图　２），否则有可能得出不客观的看法。文中虽以断裂　二＿　第３期　卜＿　＿［参考文献］　…　减分析［Ｊ］．公路，２００４，（４）：２４－２８．　…　２２５　蔡汉成等：交通荷载作用下边坡动力响应数值分析　Ｍ　Ｐｏｐｅｓｅｕ１．Ａ　ｓｕｇｇｅｓｔｅｄ　ｍｅｔｈｏｄ　ｆｏｒ　ｒｅｐｏｒｔｉｎｇ　ｌａｎｄｓｌｉｄｅ　ｒｅｍｅｄｉａｌ　［８］　樊秀峰，简文彬．交通荷载作用下边坡振动特性分析［Ｊ］．岩土　力学，２００６，２７（增）：１１９７－１２０１．　ｍｅａｓｕｒｅ［Ｊ］．Ｂｕｌｌｅｔｉｎ　ｏｆ　Ｅｎｇｉｎｅｅｒｉｎｇ　Ｇｅｏｌｏｇｙ　ａｎｄ　ｔｈｅ　ｅｎｖｉｒｏｎ－　ｍｅｎｔ，２００１，６０（１）：５８－６２，　［９］　郑黎明．铁路边坡岩体的振动特性和机理［Ｊ］，岩石力学与工　程学报，１９９４，１３（１）：６９—７８．　刘飞禹，蔡袁强，徐长节．交通荷载作用下软基加筋道路变形　机制研究［Ｊ］．岩土力学，２００９，３０（１）：１２９－１３４　姚海林，卢正，罗海宁．交通荷载作用下ｋｅｌｖｉｎ地基上不平整　［１Ｏｊ　张友葩，高永涛，王杰林．动荷载下边坡的失稳分析［Ｊ］．北京　科技大学学报，２００３，２５（２）：１１０—１１６．　路面的动力响应分析［Ｊ］．岩土力学，２００９，３０９（４）：８９１－８９６．　章东强，谢伟乎，于艳丽．交通荷载引起的环境振动实侧与分　金亮星，李建．高速铁路路堑岩质边坡的动力响应分析［Ｊ］．　振动工程学报，２００８．２１（４）：３５４—３５８．　析［Ｊ］．武汉理工大学学报，２００４，２６（９）：５７—５９．　查文华，洪宝宁．交通荷载下低路堤路基的动力响应［Ｊ］．江苏　大学学报，２００８，２９（３）：２６４－２６８．　［１２］　简文斌，胡志忠，樊秀丽，等．边坡对循环荷载的响应研究　［Ｊ］．岩石力学与工程学报，２００８，２７（】２）：２５６２－２５６７．　［１３］　叶四桥，唐红梅，李明．汽车荷载对滑坡稳定性的影响分析　［Ｊ］．重庆建筑大学学报，２００６，２８（５）：１０６一１０９．　［１４］　蒯行成，吾朝辉，李永红．车辆振动瞬态分析及动荷载［Ｊ］．中　南公路工程，２００４，２９（３）：４１－４４．　查文华，洪宝宁，徐毅．交通荷载下低路堤高速公路路面路基　振动测试与分析［Ｊ］．公路工程，２００７，３２（４）：１１３．１１７．　刘奉喜，刘建坤，房建宏，等车辆荷载作用下盐路公路振动衰　业　业　业业业韭尘业誊　业　妊　坐业　韭　业　业逝坐　业坐・　业业－：ｇｔ￣　业逝业－＆ｅ业　（上接２１９页）　［参考文献］　敬少群，吴云，乔学军，等．ＧＰＳ时间序列及其对昆仑山口西　８．１级地震的响应［Ｊ］．地震学报，２００５，２７（４）：３９４－４０１．　顾国华，张晶．中国地壳运动观测网络基准站ＧＰＳ观测的时间　序列结果［Ｊ］．大地测量与地球动力学，２００２，２２（２）：６１－６７．　张跃刚，帅平，胡新康，等．从ＧＰＳ观测看华北地区的形变场　演化［Ｊ］．大地测量与地球动力学，２００６，２６（１）：３６－４１．　江在森，张希，崔笃信，等．青藏块体东ｊＥ缘近期水平运动与变　形［Ｊ］．地球物理学报，２００１，４４（５）：６３６－６４４．　殷海涛，李杰，张玲，等．基于ＧＰＳ观测网的山东地区地壳运动　特征分析［Ｊ］．西北地震学报，２００８，３０（３）：２７６－２８１．　杨国华，江在森，刘峡，等．华北地区水平运动场动态特征及其　与昆仑山８．１级地震的可能关系［Ｊ］．大地测量与地球动力　学，２００７，２７（１）：６—１０．　［７］　杨国华，王敏，韩月萍．华北中北部地壳运动与张北地震［Ｊ］．　中国地震，２００１，１７（３）：３０１—３１１．　［８］　张希，江在森，王琪，等．青藏块体东北缘弹性块体边界负位错　反演与强震地点预测［Ｊ］．地震学报，２００５，２７（６）：６２０－６２９．　［９］　刘峡，傅容珊，杨国华，等．用ＧＰＳ资料研究华北地区形变场和　构造应力场［Ｊ］．大地测量与地球动力学，２００６，２６（３）：３３－３８．　［１Ｏ］　黄立人，郭良迁，马青．华北ＧＰＳ网复测结果的初步分析　［Ｍ］．北京：地震出版社，２００１，２１８－２２７．　杨国华，王琪，王敏，等．云南地区现今地壳水平运动与变形　特征［Ｊ］．大地测量与地壳动力学，２００３，２３（２）：７－１４．　［１２］　杨国华，江在森，武艳强，等．中国大陆整体无净旋转基准及　其应用［Ｊ］．大地测量与地球动力学，２００５，２９（４）：６－１Ｏ．　［１３］　胡新康，王倩，马青，等。区域无整体旋转基准的研究与应用　［Ｊ］．大地测量与地球动力学，２００７，２７（６）：５２－６０．