基于InSAR和GPS观测数据的尼泊尔地震发震断层特征参数联合反演研究

利用日本ALOS-2和欧空局Sentinel-1A卫星获得的尼泊尔地震同震形变场,结合GPS同震位移数据,联合反演了断层滑动分布特征和空间展布.结果表明:尼泊尔地震的同震形变场主要集中在150km×100km的范围内,且分为南北两个相邻的形变中心,南形变中心的视线向抬升量约为1.2m,北形变中心的视线向沉降量约为0.8m,均位于发震断层上盘.位于形变抬升区的KKN4和NAST两个GPS站,抬升量和南向运动量均达到了m级,而远离震区的其他GPS台水平和垂直观测量均在1cm以内.联合反演得到的断层位错分布主要集中在沿走向150km,沿倾向70km的范围内,最大滑动量为5.59m,平均滑动量为0.94m.断层面倾角在浅部约为7°,随着深度增加,倾角逐渐变大,到垂直深度20km时倾角接近12°;5月12日MW7.2级余震位于主震破裂区的"凹"型滑动缺损区域;主震破裂区的上边界与MBT空间位置十分吻合,主震破裂区主要集中的MBT以北50~60km处,垂直深度为8~9km,倾角为9°,继续向北时主震破裂面以10°~12°的倾角向深延伸,在18~20km可能与MHT交汇.因此,初步判定MBT为此次地震的发震断层.     

青藏块体东北缘弹性块体边界负位错反演与强震地点预测

利用青藏块体东北缘1993——1999年、1999——2001年和2001——2003年3个时段的GPS水平运动速度场成果,通过研究建立的弹性块体边界负位错模型反演,给出该区构造块体及其边界断裂运动与变形差异、构造应变场时空分布、高应变能积累部位以及强度的定性分析和定量描述. 在充分考虑区域构造及块体应变的前提下,研究与该区6级左右及6级以上地震地点预测相关的背景性前兆的共性特征.      

川滇菱形块体东边界地壳形变研究

根据地质资料和GPS资料。在断层分布密集的地区用模拟断层的位错替代多断层位错,建立了川滇菱形块体东边界的非震位错模型。利用GPS区域网资料．采用改进的蒙特卡罗法和遗传算法,联合反演了川滇菱形块体东边界的活动特性．得到如下认识：1）小江断裂是高应变积累闭锁区,鲜水河断裂次之,安宁河和则木河断裂应变积累程度最低。2）2001～2004年与1999～2001年相比．鲜水河北西段应变积累即锁定程度没有明显变化。南东段略有加强;安宁河略微减弱;则木河没有明显的变化;小江断裂各段的锁定程度均有显著的增强。3）各断裂的运动特点在空间上呈现出分段性。4）边界带的变形主要集中在川滇菱形块体一侧。     

地震波衰减的物理机制研究

对近年来有关地震波衰减的物理机制进行评述和初步研究认为,对于地壳介质,造成地震波衰减的主要原因是地壳内存在大量裂隙,裂纹中饱含水或部分含水,地震波传播时引起裂隙中的流体运动,从而造成地震波的衰减,对于上地幔,扩散控制的位错阻尼机制可能是造成地幔地震波衰减的一个重要原因,而上地幔软流层部分熔融的存在也是造成地震波衰减不可忽视的原因．     

走滑断层地震地表断裂位错估计方法研究

为对一次地震中可能造成的地表断裂位错作出较准确的估计,采用拟静力弹塑性有限元方法,分别对覆盖土层为粉质黏土和黏土情况下,走滑断层引发的地震地表断裂进行了数值模拟分析。根据历史震害数据回归拟合的震级M与基岩位错 的关系式以及数值计算结果,建立了震级M与地表位错 的关系式。公式中考虑了土层厚度H对地表位错 的影响,而不仅仅局限于根据震级M的大小通过统计公式来估算地表断裂位错 。结果表明,走滑断层引发的地震地表断裂位错不仅与震级的大小有关,还与土层厚度和土层性质有关;在相同震级作用下,随着土层厚度的增加,地表位错逐渐减小;在相同震级和相同土层厚度下,上覆土层为粉质黏土时产生的地表位错要大于上覆土层为黏土时产生的地表位错;根据拟合的公式估计出不同震级情况下可不考虑走滑断层影响的临界覆盖土层厚度值,有助于提高活动断裂地震危险性评估工作的可靠性。     

球体位错理论计算程序的总体设计与具体实现

本文详细介绍了球体位错理论计算程序的总体设计思想、 各类配套文件的具体内涵以及各类输出文件的物理含义, 同时介绍了程序的使用方法和注意事项, 为读者独立使用该程序提供参考。 球体位错理论计算程序主要由三部分组成: ① 位错格林函数计算程序, 基于具体的球对称地球模型提供离散的二维格林函数数值框架; ② 积分计算程序, 对离散的格林函数数值框架进行双二次样条插值运算, 并对四类独立位错源对应的格林函数进行适当组合, 从而计算出任意位置任意类型震源在地表产生的同震变形(含位移、 应变、 重力变化和大地水准面变化); ③ 辅助文件, 用于提供发震断层模型和计算点位信息。 一般情况下, 读者不需要理解位错格林函数计算程序和积分计算程序, 只需要对辅助文件提供的信息进行针对性改动, 就可以计算目标地震在目标观测站引起的同震变形。     

汶川地震形成位移场的空间分布

采用Yoshimitsu Okada及Steketee的断裂位错模型和汶川地震现有成果,通过坐标转换和合成矢量的方法将断裂位错模型用于映秀—北川断裂、灌县—安县断裂和北川—青川断裂组成的断裂系统的错动研究中,理论上计算龙门山近断裂地区的水平位移场(包括沿走向和垂直走向方向)和垂直位移场从震源到地表的分布.模型参数源于现有的研究成果和野外实地考察,计算得到的地表位移场与先前研究成果,包括GPS实测数据,具有的近似性表明了模型的正确性.但GPS只能测定地表水平和垂直变形,对于地下变形情况及其分布却无法描述;目前对于地下变形的研究主要基于对大量仪器记录的地震资料进行反演或通过野外观测进行推测;但是反演多集中于空间较大范围,这样虽可阐述断裂运动引起的大范围位移趋势,而对于震源附近空间介质位移的描述却略显粗糙;同时,野外观测误差较大.本文通过理论模型计算汶川地震中近断裂区域(距断裂50km)内的位移及分布.由计算发现在近断裂区域内垂直位移场和垂直走向方向位移场变化趋势一样,幅值都是从震源到地表逐渐减小;沿断裂走向方向的水平位移场从震源到地表逐渐变大,同时计算还表明位移场的变化在断裂上盘比下盘剧烈,余震分布主要集中于断裂上盘,这说明余震分布和位移剧烈程度存在某种相关性.     

球体位错理论在2011年日本强震中的应用研究

利用球体位错理论计算了2011年日本强震产生的远场同震位移与应变,并利用GPS远场数据修正了该强震的总地震矩。结果表明：1）利用球体位错理论计算得到的理论水平位移场显示,垂直于发震断层的广大区域同震位移较大,位移矢量总体都指向震中地区,震中距约5 000千米的地方亦产生了3 mm以上的同震水平位移。理论位移与远场GPS观测结果具有良好的一致性;2）比较两个独立断层模型对应的理论同震位移场发现,震源西部地区远场位移总体上只有1%～4%的微小差异,而东部广大海域的差异则达到同震信号的6%～15%,震中周围差异更大。该差异表明,相对于震源仅局域覆盖的日本本土GPS观测数据对2011年日本强震的断层滑动分布模型的约束能力有限;3）依据中国及邻区的远场GPS同震观测数据修正2011年日本强震的总地震矩,把该地震释放的总能量约束在(3.24～4.96)×10 ²² Nm,相应的矩震级为Mw8.97～9.10;4）2011年日本强震在华北地区产生的同震应变与该区的长期应力变化背景场大体相反,表明该强震使华北地区的地壳产生了松弛效应。     

远场GPS同震位移的变形特征分析——以2011年日本东北大地震为例

利用2011日本东北大地震的亚洲东部GPS远场水平位移数据,通过计算及比较分析4个独立震源解和倾/滑角效应,研究了远场地震变形的空间分布特征和规律。结果表明,对于低角度逆冲型地震来说,水平位移格林函数在远场位移场中占主导作用。基于此确定了远场GPS水平位移场能在反演中有效利用的空间范围,为联合远场GPS数据和近场GPS数据反演断层滑动分布作理论准备,并为远场数据反演位错Love数（h、l）提供理论依据。     

基于粘弹性球体地球模型的震后位移与重力变化计算软件

在考虑地球曲率、成层结构、可压缩性与自重的前提下,Tanaka等［1-2］提出一套较为完备的粘弹性球体位错理论,可计算全球任意位置由地震产生的同震与震后形变（含位移、重力变化、大地水准面变化）。Gao等［3］给出了与上述理论相匹配、界面友好的计算软件,能计算30个震后时间点对应的震后形变。本研究针对Gao等的软件进行改进,可计算震后任意时间点对应的震后形变。新软件由3个部分组成：1）与32个震后时间点相对应的32套离散格林函数数值框架;2）格林函数插值计算程序,可针对上述32套格林函数数值框架进行插值运算,输出任意震后时间点对应的格林函数数值结果;3）积分计算程序,调用上述格林函数数值结果,计算任意类型地震在地表任意位置产生的同震与震后形变。一般情况下,使用者只需按要求准备辅助文件,提供发震断层模型和观测站位置信息,以及震中周围地区地幔粘滞性因子,先后运行格林函数插值计算程序和积分计算程序,即可计算出目标地震在地表任意位置产生的同震与震后形变。本文基于粘弹球体位错理论与弹性球体位错理论,分别计算2011年日本MW9.0地震引起的远场同震位移,2套结果的高度一致性证明了新程序的正确性。最后,介绍需要注意的若干事项,便于使用者掌握该软件。