有序量最优分割法在计算中的改进

本文论述了“最优分割法”中,段的直径计算公式的推导分解过程,采用最简便的数学形式在计算机上进行运算,并对整个计算流程做了合理安排,充分利用了中间性数据,使计算机内存占用量及计算时间得到了较大节省。     

白兰高速公路沿线黄土场地震陷区判定与评价

采用DSD 160型电磁式振动三轴试验系统,通过对甘肃白银—兰州高速公路(国道109线)沿线黄土的震陷试验,探讨了该场地黄土在动荷载作用下的震陷特性。利用黄土震陷灾害预测方法和该地区黄土场地地震危险性分析结果,对该公路沿线黄土场地的震陷区进行了判定和评价。这些研究结果可作为该高速公路工程进行抗震设计和地震变形分析的基本依据,对黄土地区其它工程也有一定的参考价值。     

黄土场地地震液化研究

通过往返加荷动三轴试验，对兰州某民用机场扩建工程场地的饱和黄土和砂土进行了液化试验。在试验结果的基础上，运用抗液化剪应力判别方法、地震危险性分析计算结果以及根据有限元一维模型计算求得的该场地的地震剪应力，对该场地饱和黄土和砂土的地震液化进行了综合判断。结果表明，在未来遭受到50年10％超越概率的地震作用时，该场地的饱和黄土比饱和砂土更容易发生液化。     

中国境内天山的岩石圈速度结构——P波与S波接收函数的联合反演

通过合成地震图的计算,本文系统地分析研究了S波接收函数的运动学及动力学特征,并在接收函数非线性复谱比反演方法的基础上,进一步发展了基于贝叶斯理论的P波和S波接收函数的非线性联合反演方法. 利用本文发展的P波和S波接收函数的非线性联合反演方法和沿独库公路布设的由51个宽频带流动地震台站组成的地震台阵记录的远震体波波形数据得到了中国境内天山300 km深度范围地壳上地幔的P波和S波速度结构,得到了中国境内天山的岩石圈结构和造山动力学模型的一些新认识.     

剪切波速判别方法对黄土液化判别适用性研究

最近的研究结果表明饱和黄土甚至高含水率的黄土也具有很高的液化势和流态破坏势,而工程设计时关于黄土液化的判别,尚无统一标准,亦无经验可循.用Seed-Idriss简化判别法对兰州某民用机场的饱和黄土和砂土进行了液化判别,并以此结果为依据检验了现有剪切波速判别式对黄土液化判别的适用性.结论显示该场地的黄土剪切波速值过大,《岩土工程勘察规范》与Andrus和Stokoe剪应力比法都不适用.     

基于贝叶斯理论的接收函数与环境噪声联合反演

基于Bayes反演理论(Tarantola,1987,2005),在接收函数非线性复谱比反演方法基础上(刘启元等,1996),本文讨论了接收函数与地震环境噪声Rayleigh波相速度频散的联合反演.本文采用修正后的快速广义反射/透射系数方法(Pei et al., 2008,2009) 计算Rayleigh波相速度频散, 并引入地壳泊松比的全局性搜索.数值检验表明:(1)接收函数与环境噪声的联合反演能够有效地解决反演结果对初始模型依赖的问题,即使对地壳速度结构仅有非常粗略的初始估计(例如,垂向均匀模型),本文方法仍能给出模型参数的可靠估计;(2)由于环境噪声与接收函数在频带上的适配性明显优于地震面波,接收函数与环境噪声的非线性联合反演能更好地约束台站下方近地表的速度结构;对于周期范围为2~40s的环境噪声相速度频散,利用本文方法能够可靠推测台站下方0~80 km深度范围的S波速度结构, 其浅表速度结构的分辨率可达到1 km; (3)本文方法能够可靠地估计地壳泊松比,泊松比的全局性搜索有助于合理解释接收函数和环境噪声的面波频散数据.利用本文方法对川西台阵KWC05台站观测的接收函数与环境噪声的联合反演表明,该台站下方地壳厚度为44 km,上地壳具有明显的高速结构,24~42 km范围的中下地壳具有低速结构.该台站下方地壳的平均泊松比为0.262,壳内低速带的泊松比为0.27.     

汶川地震区地壳上地幔三维S波速度结构初步研究

2008年5月12日汶川MW7.9地震发生在龙门山断裂带。龙门山断裂带及其邻域的地壳上地幔三维速度结构的研究对于理解汶川大地震的动力学背景具有重要的意义。2006年10月至2009年10月,在国家重大基础研究项目(973)的支持下,中国地震局地质研究所地震动力学国家重点实验室在川西地区(26°~32°N,100°~105°E)布设了由297台宽频带数字地震仪组成的流动观测台阵(简称川西台阵)。根据川西台阵记录的环境噪声和远震波形数据,利用噪声成像技术和接收函数方法,我们研究了川西地区(29°~32°N,100°~105°E)地壳上地幔100km深度范围内的三维S波速度结构。本文得到的结果为研究川西高原和四川盆地的地壳结构提供了新的高分辨率观测证据。我们的结果表明:1)观测台阵覆盖的川滇地块、松潘-甘孜地块和四川盆地的地壳上地幔S波速度结构具有显着差异,龙门山断裂和鲜水河断裂带,作为地块间的边界断裂带,对两侧地壳结构具有明显的控制作用。2)观测台阵覆盖区域的地壳厚度存在明显差异,川滇地块的地壳厚度为60~64km,松潘-甘孜地块的地壳厚度为52~56km,四川盆地前陆的地壳厚度为46~52km,沿龙门山断裂带松潘-甘孜地块和四川盆地形成镶嵌结构,汶川地震震中处南北两侧的壳幔边界存在约6km的断错。3)四川盆地前陆低速特征表明相应区域存在厚度8~10km的沉积盖层,松潘-甘孜地块和川滇地块的中下地壳具有大面积分布的S波低速区,松潘-甘孜地块地壳平均泊松比高达0.29~0.31,汶川地震余震绝大多数分布在低速区上方的高速介质区域内,而四川盆地的中下地壳呈现整体性的高速特征,以汶川地震的震中为界,龙门山断裂带北段和南段的S波速度结构显示了明显的速度分段特征,其北段的S波速度总体上高于南段。4)本文给出的研究区地壳三维S波速度结构表明,川西高原中下地壳较为软弱,而四川盆地中下地壳的强度应明显高于松潘-甘孜地块,意味着四川盆地坚硬中下地壳可以阻挡松潘-甘孜地块向东的逃逸;另一方面,川西高原和川滇地块的中下地壳虽然均存在大面积的S波低速区,但松潘-甘孜地块内的地壳速度结构相对来说较为复杂,并形成了高、低速相间的结构特征,表明在四川盆地的阻挡作用下,该地块形成了折皱变形的结构。5)与S波低速区相应,松潘-甘孜地块和川滇地块中下地壳应处于部分熔融的状态,这对该区域存在中下地壳通道流(Channelflow)的推断是一个支持;但是,松潘-甘孜地块内是否存在中下地壳通道流仍有待进一步的深入研究。6)接收函数方位各向异性的偏振分析表明,以汶川地震震中为界,龙门山断裂西南侧处于挤压状态,而其东北侧的主压应力方向与断层走向大体平行,推断先存应力场可能驱动了汶川地震逆冲破裂之后沿龙门山断裂向北东方向的走滑破裂。     

川西龙门山及邻区地壳上地幔远震P波层析成像

本文利用川西地震台阵记录到的远震P波走时数据和非线性层析成像算法,获得龙门山地区400 km深度范围内的三维P波速度结构.为了适应川西地区复杂的地质结构,本文的层析成像方法采用了快速行进三维走时计算算法和Tarantola非线性反演算法.我们的结果揭示了川滇地块、松潘-甘孜地块和四川盆地三个不同地块构造差异及该区深部动力学特征.本文的研究表明：1）研究区地壳上地幔P波速度结构具有较为明显的分区特征,松潘-甘孜地块和川滇地块岩石圈速度较低,四川盆地岩石圈速度较高,四川盆地的岩石圈厚度从南250 km向北逐渐减薄至100 km.松潘-甘孜地块上地幔存在地幔上涌的特征.2）川滇地块和四川盆地仅是垂直接触关系,而在龙门山地区四川盆地前缘存在减薄的现象,并伴随松潘-甘孜地块上地幔低速物质有侵入四川盆地岩石圈下方的特征,这显示了四川盆地与松潘-甘孜地块和川滇地块的动力学关系的差异.3）以映秀为界,龙门山断裂带被从松潘-甘孜侵入的低速异常分为南北两段：龙门山南段和龙门山北段,汶川大地震及其余震序列均分布在龙门山断裂带的北段.在青藏高原向东挤压和地幔上涌的双重作用下造成松潘-甘孜地块隆升,由于汶川处于龙门山北段的最南端,应力容易在此集中.这些因素可能是汶川M_S8.0地震的基本动力学背景.本文的结果不支持四川盆地的俯冲及层间流动的动力学模型.     

根据接收函数反演得到的首都圈地壳上地幔三维S波速度结构

利用2002~2003年中国地震局地质研究所台阵实验室以唐山大震区为中心布设的40个流动宽频带地震台站和首都圈数字台网的33个宽频带台站的远震数据,采用接收函数非线性反演方法得到其中72个宽频带台站下方60 km深度范围内的S波速度结构.根据得到的各台站下方地壳上地幔的S波速度结构,并综合刘启元等（1997）用接收函数非线性反演方法得到的延怀盆地15个宽频带流动台站下方的地壳上地幔S波速度结构模型,给出了39°N～41°N,114°E~119.5°E区域内沿不同走向、不同深度S波速度分布.由于综合了利用首都圈数字地震台网的宽频带台站以及流动地震台阵的观测数据,本文给出了较前人同类研究空间分辨率更好的结果.结果表明: (1)研究区的速度结构,特别是怀来以东的速度结构十分复杂.在10~20 km深度范围内,研究区地壳具有高速和低速异常块体的交错结构.研究区中上地壳速度结构主要被与张渤地震带大体重合的NW向高速条带和穿越唐山大震区的NE向高速条带所控制,而其中下地壳的速度结构主要为延怀—三河—唐山地区上地幔隆起所控制.(2)研究区内存在若干壳内S波低速体,它们主要分布在唐山,三河及延怀盆地等地区.在这些地区,壳内低速体伴随着壳幔界面的隆起和上地幔顶部速度结构的横向变化.(3)地表断层分布与地壳速度结构分区有较好的相关性,表明断层对不同块体有明显的控制作用.其中,宝坻断裂,香河断裂和唐山断裂均为超壳断裂.(4)首都圈内大地震的分布与壳内低速体及上地幔顶部的速度结构有密切关系.对于唐山大地震的成因,仅考虑板块作用引起的水平应力场是不够的,有必要充分重视由于上地幔变形引起的地壳垂直变形和上地幔物质侵入造成的热效应.     

汶川M_S 8.0地震:龙门山断裂两侧地壳各向异性的初步研究

利用中国地震局地质研究所地震动力学国家重点实验室在川西地区布设的大规模密集流动宽频带地震台阵记录的远震P波波形数据,分别采用波形互相关和加权叠加方法研究了地壳各向异性。作为初步结果,得到了龙门山断裂两侧4个宽频带流动台站的接收函数莫霍面Ps转换震相快波方向和分裂时间延迟。结果表明:1)波形互相关方法总体上优于方位加权叠加方法,它不但可以给出快波方位,而且可以给出快慢波的时间延迟;方位加权叠加方法测量结果存在不确定性,其原因在于难以确定介质各向异性对称轴的方位;同时采用上述2种不同方法研究地壳各向异性参数有助于判断测量结果的可靠性;2)四川盆地内快波偏振方位基本一致,表明四川盆地地壳整体性较好,横向非均匀变形较弱;3)以汶川地震主震的震中区为界,松潘-甘孜地块北侧快波偏振方位与龙门山断裂近于平行,表明在四川盆地坚硬地壳的阻挡作用下,龙门山断裂附近松潘-甘孜地块北侧可能存在中下地壳软弱物质沿断层向NE方向的扩张变形,而其南侧处于正向挤压的状态。该结果有助于解释汶川MS8.0地震的单侧破裂过程及其余震发育的特征