昆仑山MS81地震震源参数的多破裂段模拟研究

昆仑山MS81地震的已有研究结果在破裂带长度、破裂面方向、破裂面大小等震源破裂特征参数方面存在较大差异.本文采用D-InSAR技术首次获得昆仑山MS81地震干涉同震形变场，结合野外科学考察的实测值，进行了主破裂带InSAR视线向变化量的分解，通过对InSAR分解结果、野外科学考察、遥感解译等多源数据综合分析，重新划分了昆仑山地震的次级破裂段.进而通过对地震南北盘同震应变的分析，发现了昆仑山地震的南北两盘分别受挤压和拉张两种应力作用，研究表明多种岩石在拉张和压力作用下其最小主应力下的杨氏模量表现出非线性弹性特征，从而提出对昆仑山地震地表位移及震源特征参数分析时应考虑非线弹性介质导致的非线性弹性位移分布特征.基于上述原因，本文对Okada线弹性位错模型的算法进行了改进，提出了“多震源、非均一位错分量、多破裂段叠加”的线弹性模型，该模型模拟出的形变场干涉纹图较好地体现了地震形变场的分布特征，并由此获得了一套较为完整的地震发震断层的几何学特征参数，为破裂带长度、破裂面方向、破裂面大小等震源破裂特征参数研究提供了较好的解释.     

昆仑山断层围陷波的分析和研究

对2001年昆仑山口西Ms81级地震产生的断层带,布设了沿断层和横跨断层的两条人工地震测线.通过对观测资料的定量分析和处理,求得了昆仑山断层带内部的细结构.分析工作包括从S波震相开始的振幅谱计算、速度频散计算、群速度测量,并用面波频散方法反演S波速度结构,用振幅谱比的方法估计断层带的Q值.野外试验结果表明,S波震相与围陷波组的时间差随炮点与台站之间距离增大而增加,在断层带外的测点上观测到与断层带相关的场地效应.最后得出昆仑山断层带宽度为250m、速度结构为断层内低速的分层结构和Q值为15（断层内）和30（围岩）.虽然昆仑山口西地震的震级比美国加州Landers地震的震级（Ms76）大,且地震产生的破裂带长度长得多,但是这两个地震断层带的宽度却相差不大.     

2001年昆仑山口西8.1级地震同震形变场特征的初步分析a

利用差分干涉雷达测量技术获取的宏观震中区的同震形变场，结合对地震活动性、震源机制、野外考察等资料分析，对昆仑山口西8.1级地震同震形变场特征进行了研究. 结果表明:宏观震中位于库赛湖东北侧，宏观震中区发震断层可分为两个形变中心区域，其中西段长约42 km，东段长约48 km，整个发震断层主破裂段长90 km；由干涉形变条纹分布格局可清楚地判断出发震断层的左旋走滑特征；断层两盘变形特征不同，南盘变形程度明显大于北盘；宏观震中附近最大斜距向位移量为288.4 cm，最小斜距向位移量为224.0 cm，宏观震中发震断层最大左旋水平位错为738.1 cm，最小地面左旋水平位错为551.8 cm.      

D-InSAR技术应用于汶川地震地表位移场的空间分析

基于星载合成孔径雷达差分干涉测量技术(D-InSAR),利用7个条带共112景日本ALOS/PALSAR raw格式雷达数据,采用两通差分干涉处理模式,获取了2008年5月12日汶川MS8.0地震发震断层周围约450km×500km区域的同震形变干涉纹图。通过对干涉纹图的定性分析,确定了非相干带的分布范围,据此对相位连续条带和相位不连续条带采用不同的相位解缠方案,实现了7个条带的成功解缠,获得了数值化的干涉形变场图像,并通过形变等值线和跨断层形变剖面线等方法对干涉形变场的空间分布和演化特征进行了分析。结果表明:汶川地震造成的地表形变场沿映秀-北川断裂带分布,形变范围很大,但主要集中在发震断层南北两侧各约100km的近场区。其中断层附近由西向东宽约3015km,长约250km的区域为非相干带,是本次地震中变形最强烈并伴有地表破裂发生的区域,其形变梯度已超出InSAR测度能力。在非相干带两侧宽度各约70km,具有清晰可辨连续完整并向发震断层收敛的包络状干涉条纹区域是次一级形变区,距离发震断层越近,形变梯度和幅度越大,其视线向位移为北盘沉降,南盘抬升。相对于数据条带南北边缘,北盘最大累积沉     

2000年姚安MS65地震震源及近邻区的低QS

利用震源距23 km范围内观测的2000年姚安MS65地震余震记录，计算了震源及近邻区域的尾波规一直达S波在频率15～20 Hz之间的傅里叶谱振幅.结果显示谱振幅随震源距增大而增大， 在对谱振幅进行了震源辐射方向性校正之后, 才出现谱振幅随震源距衰减的现象.由此获得了震源及近邻区域S波的Q(f)值,可表示为QS(f)=89f098其值比由尾波得出的姚安地区的平均QC(f)值低得多，表明了震源破裂带的强烈非均匀性对QS(f)的重大影响.     

2000年姚安MS65地震震源及近邻区的低QS

利用震源距23 km范围内观测的2000年姚安MS65地震余震记录，计算了震源及近邻区域的尾波规一直达S波在频率15～20 Hz之间的傅里叶谱振幅.结果显示谱振幅随震源距增大而增大， 在对谱振幅进行了震源辐射方向性校正之后, 才出现谱振幅随震源距衰减的现象.由此获得了震源及近邻区域S波的Q(f)值,可表示为QS(f)=89f098其值比由尾波得出的姚安地区的平均QC(f)值低得多，表明了震源破裂带的强烈非均匀性对QS(f)的重大影响.     

根据2001年Mw 78可可西里强震InSAR同震测量结果反演东昆仑断裂两侧地壳弹性介质差异

利用2001年 Mw 78 可可西里强震InSAR同震测量结果，反演了青藏高原北部东昆仑断裂两侧地壳弹性介质差异.InSAR测量结果显示断层南侧的同震位移比北侧的大20％～30%.根据人工地震反射剖面建立岩石圈模型，以断层两侧杨氏模量差异和震源破裂深度为反演变量，通过有限元方法模拟实测得到的同震位移剖面.反演得到最佳断层破裂深度为20～22km，断层南侧杨氏模量相对北侧比值为81%～92%.结果表明，断层两侧弹性介质性质存在明显差异，由于构造运动作用，断层南部地壳不及北部地壳坚硬.前人利用地震层析成像和大地电磁测深等手段推断青藏高原内昆仑山断裂以南可可西里－羌塘地块地壳内广泛发育低速高导层，我们通过形变场力学分析得到与此相一致的结果.     

根据2001年Mw 78可可西里强震InSAR同震测量结果反演东昆仑断裂两侧地壳弹性介质差异

利用2001年 Mw 78 可可西里强震InSAR同震测量结果，反演了青藏高原北部东昆仑断裂两侧地壳弹性介质差异.InSAR测量结果显示断层南侧的同震位移比北侧的大20％～30%.根据人工地震反射剖面建立岩石圈模型，以断层两侧杨氏模量差异和震源破裂深度为反演变量，通过有限元方法模拟实测得到的同震位移剖面.反演得到最佳断层破裂深度为20～22km，断层南侧杨氏模量相对北侧比值为81%～92%.结果表明，断层两侧弹性介质性质存在明显差异，由于构造运动作用，断层南部地壳不及北部地壳坚硬.前人利用地震层析成像和大地电磁测深等手段推断青藏高原内昆仑山断裂以南可可西里－羌塘地块地壳内广泛发育低速高导层，我们通过形变场力学分析得到与此相一致的结果.     

基于线弹性位错模型及干涉雷达同震形变场反演1997年玛尼Mw75级地震参数－Ⅰ.均匀滑动反演

1997年11月8日西藏玛尼Mw75级地震是干涉雷达技术应用于地震观测以来的一次重要事件.Peltzer等认为该形变场体现了地壳的非线性弹性形变响应，挤压象限和拉张象限的杨氏模量比值为2.相同的情形并没有在其他地震同震形变场中发现，通常用线弹性理论就能够很好地解释观测数据.考虑到他们所用断层模型的简化程度和纯走滑约束等，本研究认为非线性弹性解释是牵强的.应用广泛使用的Okada线弹性位错模型，采用卫星观测得到的地表断层位置，去除倾滑为零的约束，基于非线性优化反演方法寻求拟合雷达观测的最佳断层几何参数和均匀滑动参数.结果表明，线弹性模型能够满意解释观测数据.断层在朝阳湖以东的段落最深达到165 km，随着断层向两侧延伸，深度逐渐减小；反演得到的断层倾     

基于线弹性位错模型及干涉雷达同震形变场反演1997年玛尼Mw75级地震参数－Ⅰ.均匀滑动反演

1997年11月8日西藏玛尼Mw75级地震是干涉雷达技术应用于地震观测以来的一次重要事件.Peltzer等认为该形变场体现了地壳的非线性弹性形变响应，挤压象限和拉张象限的杨氏模量比值为2.相同的情形并没有在其他地震同震形变场中发现，通常用线弹性理论就能够很好地解释观测数据.考虑到他们所用断层模型的简化程度和纯走滑约束等，本研究认为非线性弹性解释是牵强的.应用广泛使用的Okada线弹性位错模型，采用卫星观测得到的地表断层位置，去除倾滑为零的约束，基于非线性优化反演方法寻求拟合雷达观测的最佳断层几何参数和均匀滑动参数.结果表明，线弹性模型能够满意解释观测数据.断层在朝阳湖以东的段落最深达到165 km，随着断层向两侧延伸，深度逐渐减小；反演得到的断层倾     

孕震区震前D-InSAR干涉形变场动态演化图像分析

采用D-InSAR技术获取了1997年11月8日西藏玛尼地震前3期地表干涉形变场时空演化图像。通过对干涉形变场图像进行分析,得到了震前积累形变场时空演化特征与震源破裂面闭锁之间的对应关系。结果表明:玛尼地震前10个月孕震区地表形变场就开始出现了与发震断层性质一致的变化趋势———左旋扭动,且随着发震时间的临近,南盘形变加剧,形变中心向E进一步扩展。南盘变形量明显大于北盘,南盘为主动盘。到震前2个半月,断层两盘的闭锁进一步加剧,有随时破裂的可能。从震前1996年4月16日到震前2个半月之间,南盘局部累计扭曲变形量达344mm     

2004年7月11日西藏Mw62级地震震源破裂过程研究

基于有限断层模型，利用远场体波波形数据研究了2004年7月11日西藏MW62级地震的震源破裂过程.结果表明该地震是一个以倾滑为主的浅源正断层型地震，震源深度为125km，断层面走向152°，倾角44°，平均倾滑角-117°.破裂在震中处成核，然后以28km/s的平均速度向两侧传播，在震中以东偏北5km处达到最大滑动43cm.该地震主张力轴近E W方向，受浅部NNW SSE或N S向裂谷带控制，青藏高原南部的逆冲运动是引发这次地震的直接原因.     

基于D-InSAR技术的西藏改则地震同震形变场特征分析

2008年1月9日在西藏改则地区发生6.9级地震, 接着1月16日又发生6.0级余震, 表现为双震型。 本文利用差分干涉测量技术(D-INSAR), 通过对欧空局ASAR数据进行处理获取了其同震干涉形变场。 通过分析表明: 改则地震干涉形变场呈双贝壳状, 影响范围约33 km×30 km, 以北东向地震破裂带为分界线分为西北视线向沉降盘与东南视线向隆升盘, 最大视线向沉降形变量约53.2 cm, 最大视线向隆升形变量约11.3 cm。 西北沉降盘又存在东、 西两个形变中心, 推测西部形变中心受6.0级余震的控制, 东部形变中心受6.9级主震的控制。 宏观震中位置应位于左旋走滑改则—洞错断裂与依布茶卡—日干配错断裂的分阶部位(左阶), 构造应力场以张性拉伸为主, 导致地震破裂为典型的正断层破裂, 与此次地震的干涉形变场特征及哈佛大学震源机制解吻合。     

利用星载D-INSAR技术获取的地表形变场研究张北-尚义地震震源破裂特征

近 10年来 ,合成孔径雷达差分干涉测量技术 (D INSAR ,DifferentialInterferometricSyn theticApertureRadar,简称 :雷达差分干涉测量 )取得了令世人瞩目的成绩 ,已成为极具有潜力的空间对地观测新技术。但如何利用D INSAR测量结果来揭示震源破裂信息 ,成为目前被关注的问题。本文简单地论述了震源位错模型的基本原理 ,并以 1998年 1月 10日张北 -尚义地震为例 ,利用差分干涉形变场图像 ,采用弹性半空间介质中的震源位错模型 ,正演了发震断层某些几何学和运动学性质。结果表明 :张北 -尚义地震破裂面属左旋逆走滑断层 ,走向SEE NWW2 72° ,倾角 4 6° ;破裂方向为由SEE至NWW的单侧破裂 ;破裂面长 9km ,宽 8km ,深度 8km ;三个方向的位移矢量分别为 :2 90mm、5 6 0mm和 0mm。     

基于实测值分段线性内插模型的InSAR视线向同震位错分解-以青藏高原昆仑山MemＳsub8.1地震为例x,auto,auto,415px);}stylehf=httprztl.com business week launch

利用现场GPS定位的实测值，在尝试过最小二乘拟合并取得初步成果的基础上，进而采用更符合破裂带形态的线性立方插值函数，在破裂带主断面上，建立起一种理论与实测相结合的InSAR视线向(LOS)变化量的分解方程，得到了InSAR视线向位移分解的具有唯一性的解析解. 本文的解析法吸纳了GPS定点实测值的精度优势，利用InSAR全天候、准实时获取连续形变场的技术特点，通过数学近似，最终获得了主破裂带上连续变化的水平位错及垂直位错同震形变曲线.      

利用GPS和InSAR数据反演2011年日本东北MW9.0地震断层的同震滑动分布

本文首先对Envisat/ASAR数据进行干涉处理, 获取2011年日本东北M_W9.0地震的地表InSAR同震形变场; 然后通过对InSAR同震形变数据重采样方法的深入分析, 选择条纹率法结合干涉图的空间相干性对InSAR同震形变数据进行重采样; 最后基于弹性半空间位错模型, 联合InSAR与GPS形变数据, 采用最小二乘法反演发震断层的滑动分布. 研究结果表明: 考虑相干性的条纹率重采样方法, 更适用于形变场中存在除断层外的有限边界、 且形变场范围较大的InSAR数据重采样处理; 断层滑动主要发生在地表以下50 km范围内, 最大滑动量为49.9 m, 矩张量为4.89×1022 N·m, 所对应的矩震级为M_W9.1, 与地震学反演的结果比较吻合.      

2008年新疆乌恰Mw6.7地震震源机制与形变特征的InSAR研究

2008年10月5日新疆乌恰M_w6.7级地震发生在南天山、帕米尔高原及塔里木盆地交汇地带,基于地震波反演的震源机制解确定的震源深度存在较大差异.本文利用日本ALOS卫星的PALSAR图像,获得了本次地震的同震形变场,基于卫星视线向（LOS）和方位向（Azimuth）的形变,采用均匀弹性半无限位错模型和有界最小二乘（BVLS）算法,以网格矩形位错元法对发震断层的几何产状、滑移及分布进行了估算,结果表明本次地震以逆断破裂为主,断层面上最大位错量接近3.4 m,形变中心位于73.8040°E,39.5335°N,深度约5 km,震级估算为M_w6.6;地震发生在走向46°,倾角48°的断层上,发震断层长30 km,宽14 km,闭锁深度9 km,符合该地区浅源地震多发的构造特点,发震断层为乌合沙鲁断裂带.InSAR反演的滑移形变主要集中于地下2~7 km,表明乌恰地震为浅源地震,可能与该断层附近历史地震未完全释放的残余应力积累有关.同时,InSAR反演的断层位错分布呈现双破裂特征,震级分别为M_w6.5和M_w6.1,可能与本次地震的主震和余震相对应,也可能是由主震激发而产生的两组破裂.     

基于InSAR分析2015年尼泊尔MW7.8地震形变场特征及断层滑动分布反演

2015年4月25日尼泊尔发生了M_W7.8地震, 本文基于震前、 震后两景Sentinel-1A雷达影像, 采用D-InSAR两轨差分干涉法提取了此次地震的同震形变场。 结果显示, 同震形变场位于喜马拉雅造山带—主边界逆冲断裂(MBT)和主前锋逆冲断裂(MFT)附近, 形变场整体表现为自西北向往东南方向延伸近150 km的纺锤形包络状, 以大面积隆起抬升形变为主, 视线向最大隆升形变达1.18 m, 抬升区北侧存在一小沉陷区, 以InSAR观测值定位同震最大形变中心。 基于均匀介质弹性半空间模型(Okada模型)与InSAR观测数据反演了断层滑动分布。 反演结果表明该地震属于典型逆冲型地震, 发震断层为主喜马拉雅逆冲断裂(MHT), 同震破裂从主喜马拉雅逆冲断裂(MHT)向上沿着主前锋逆冲断裂(MFT)传递。 基于InSAR同震形变场局部形变细节, 结合震区地质背景、 断裂分布及断层运动特征, 获得了同震破裂拟出露地表迹线。     

基于InSAR同震形变场反演汶川Mw7.9地震断层滑动分布

通过综合分析2008年5月12日汶川地震野外地震地质考察的地表破裂带空间分布及分段资料,结合InSAR干涉形变场资料,构建了五段断层几何结构模型,该模型与野外地震地质考察结果在多数分段上基本一致;基于此五段断层模型,运用敏感性迭代拟合算法反演了汶川地震InSAR同震形变场,获得了断层滑动分布及部分震源参数.结果表明,基于余震精定位获得的地震断层倾角模型模拟的同震形变场与InSAR形变场吻合较好,且残差较小;反演的滑动分布主要集中于地下0～20 km,最大滑动量分别位于北川及青川等地区,最大可达到10 m;沿SW-NE走向,断层面的滑动方向主要以右旋兼逆冲形式为主,在汶川及都江堰地区以强烈的逆冲为主兼有一定右旋走滑分量,在北川及映秀地区以逆冲兼右旋运动为主,在平武及青川等地区则逐渐过渡为以右旋运动为主兼有一定的逆冲分量,其中汶川地区的平均滑动角为97°,北川地区的平均滑动角为119°,青川地区平均滑动角为138°.反演矩张量为7.7×10²⁰ N·m,矩震级达M_w7.9.     

基于InSAR同震形变场反演汶川M_w7.9地震断层滑动分布

通过综合分析2008年5月12日汶川地震野外地震地质考察的地表破裂带空间分布及分段资料,结合InSAR干涉形变场资料,构建了五段断层儿何结构模型,该模型与野外地震地质考察结果在多数分段上基本一致;基于此五段断层模型,运用敏感性迭代拟合算法反演了汶川地震InSAR同震形变场,获得了断层滑动分布及部分震源参数.结果表明,基于余震精定他获得的地震断层倾角模型模拟的同震形变场与InSAR形变场吻合较好,且残差较小;反演的滑动分布主要集中于地下0～20 km,最大滑动量分别位于北川及青川等地区,最大可达到10 m;沿SW-NE走向,断层面的滑动方向主要以右旋兼逆冲形式为主,在汶川及都江堰地区以强烈的逆冲为主兼有一定右旋走滑分量,在北川及映秀地区以逆冲兼右旋运动为主,在平武及青川等地区则逐渐过渡为以右旋运动为主兼有一定的逆冲分量,其中汶川地区的平均滑动角为97°,北川地区的平均滑动角为119°,青川地区平均滑动角为138°.反演矩张量为7.7×10~(20)N·m,矩震级达M_w7.9.