黏弹性介质中瑞雷波有限差分数值模拟与波场分析

瑞雷面波经常被用来反演地表浅层横波速度,受到越来越广泛的关注。对瑞雷波的研究一般都基于完全弹性介质,而实际地层更接近黏弹性介质,对黏弹性介质中的瑞雷面波进行模拟更具实际意义。本文采用广义标准线性体模型来描述黏弹性介质,并采用交错网格有限差分法对考虑水平自由表面的黏弹介质进行正演模拟,再与弹性介质中的结果进行对比分析。首先采用非线性最优化算法根据期望常数品质因子直接求取松弛时间来拟合常Q模型,并给出广义标准线性固体的具体算例,实施自由表面条件时采用声学-弹性边界近似法,通过剪切模量不变来考虑自由表面上、下横向应力保持连续的条件。对于非自由表面,采用非分裂的多轴卷积完全匹配层来吸收波场。然后对几种典型的数值模型进行正演模拟计算,数值解与解析解的对比验证了本文方法的准确性与有效性,正演结果的对比表明波场尤其是面波频散会受黏弹性影响,因此有必要在面波勘探中考虑黏弹性因素。      

利用远震研究四川盆地及其周缘瑞利面波相速度和方位各向异性

基于四川盆地及周边的245个宽频带台站2010年9月—2014年9月期间的远震记录,提取双台路径瑞利面波相速度频散资料,反演得到四川盆地20～120s的高分辨率瑞利面波相速度及各向异性空间分布.在丰富区域地球物理基础数据的同时,结合已有研究成果对地壳上地幔变形耦合进行探讨,结果表明短周期(20～30s)的相速度分布与四川盆地的地质构造特征相吻合,作为川滇地块、松潘—甘孜地块和四川盆地之间的边界——龙门山断裂带和鲜水河断裂带对上述三个地块上地壳的速度结构具有明显的控制作用;松潘—甘孜地块,特别是川滇地块中下地壳普遍表现为明显的低速异常,表明中下地壳相对软弱;而四川盆地的中下地壳整体呈现相对高速,表明四川盆地具有相对坚硬的中下地壳.研究区域东南角接近北扬子地块与南扬子地块的缝合部位,呈现高速异常.四川盆地南部和东南邻区不同周期均具有较强的各向异性,且快波方向较为一致,反映这些地区不同深度变形耦合较好.四川盆地西部、北部及东北部邻区,不同周期的各向异性快波方向变化较大,不同深度变形耦合较差.这些特征与绕喜马拉雅东构造结的物质流动被扬子地块的高速地壳阻挡的宏观认识基本一致.     

鄂尔多斯西缘北段的地壳结构和块体间变形关系

本文对喜马拉雅计划二期部分台站的远震波形数据进行接收函数提取,利用接收函数共转换点叠加方法获得阿拉善地块、鄂尔多斯地块以及银川—河套盆地下方0~80 km深度的速度间断面结构.结果表明:鄂尔多斯地块成层性好,地壳厚度为38~42 km,康拉德界面为18~22 km,阿拉善地区的Moho面深度为38~45 km.河套盆地地壳厚度约52 km,银川断陷盆地和贺兰山下方的Moho面最深为~55 km.鄂尔多斯西缘构造边界下方Moho面变化明显,且黄河断裂为深大断裂直接切割莫霍界面.根据本文的间断面成像结果我们进一步确定阿拉善地块与鄂尔多斯地块分属不同的大地构造单元.与此同时,我们推测贺兰山以西70~80 km范围内和鄂尔多斯地块西缘北段存在地壳增厚变形的可能.     

青藏高原东南缘热流估算及与地震活动相关性分析

青藏高原东南缘地区内部构造运动强烈,是地热资源发育与地震事件频发的活动地区.大地热流记录了发生在地球深部各种作用过程的热学信息,可以作为地质构造活动和地震活动研究的有效约束,但是大范围的热流数据测量很难实现,因此,本文根据居里面深度结合放射性元素分布等计算了青藏高原东南缘的大地热流分布.首先,通过地表放射性元素的分布计算出地表产热量的分布,然后,利用相关地热参数之间的关系迭代计算出该地区地壳上下层的热导率分布,最终估算出地表热流及地下不同深度处热流值的分布.本文结果表明:(1)青藏高原东南缘的大地热流位于44～108mW·m~(-2)之间,平均75mW·m~(-2),符合研究地区西南高、东北低的背景趋势,地壳内部热流值随深度的增加而降低.大部分地区地表热流异常与实际地热带分布相吻合,如川西、藏东南与滇西地区等地为地热高值区,川东和楚雄等地为热流低值区.(2)结合其他地球物理探测结果,总结了地壳内部热流与地震事件的联系:在地热梯度带地区,当两侧地层在一定深度范围内存在明显物性差异时,地震事件高发.     

背景噪声和地震面波联合反演渭河盆地 及邻区壳幔S波速度结构

利用陕西及邻区测震台网和中国地震科学台阵探测项目共257个宽频带台站记录的连续波形与远震数据,采用基于射线追踪的面波频散直接反演方法获得了渭河盆地及邻区地壳上地幔顶部S波速度结构,成像结果显示:1)渭河盆地顶部形成于新生代的沉积层造成其浅部显著的低速异常,盆地中、上地壳为低速结构,低速带延深至约25km深处,莫霍面相对两侧突变上隆,上地幔高速体侵入下地壳,可能与中—新生代上地幔基性-超基性铁镁质物质底侵有关。2)南鄂尔多斯块体地壳浅层东薄西厚的低速结构可能与块体遭受的整体掀斜、差异性抬升和强烈而不均匀的剥蚀有关。壳内不存在明显的低速异常,说明壳内低速体并没有贯穿整个鄂尔多斯地块,鄂尔多斯南段仍保留着稳定克拉通属性,至今还未遭受明显改造。3)秦岭造山带东、西深部结构存在差异,具有分段特征。造山带下地壳底部的低速异常,可能与造山带受青藏高原东北缘隆升和向外扩展等构造活动的影响有关,分析认为秦岭造山带存在青藏高原物质E流的下地壳流通道的可能性不大。     

武当地区构造类型及其控矿作用

文中主要讨论了武当地区构造类型及控矿作用。构造类型基本包括伸展滑脱构造、推覆构造和脆性构造,分别发生在三个不同的时期。伸展滑脱构造主要存在于前印支期,挤压推覆构造存在于印支期,脆性断裂出现在印支晚期。多重韧性滑脱面出现在武当群地层形成时的不同构造岩性界面上;不同层次顺层伸展滑脱面(DF1—DF5)常出现在堆垛的构造岩性面上,并优先富集Au、Ag矿;逆冲推覆构造使先存滑脱面变形,伴随褶皱作用使成矿元素流动、搬运、富集。     

2002年6月21～23日梅雨锋上的MCS分析

对2002年6月21～23日发生在黄海-长江流域梅雨锋上的一系列由中尺度对流系统(MCS)引起的强对流天气过程进行了云顶辐射亮温(TBB)和天气形势、气象要素等的诊断分析.结果表明梅雨锋上的MCS是造成这次长江流域产生暴雨甚至大暴雨级强降水的主要系统;MCS的强降水中心(区)与TBB≤-32℃的频率FM、TBB≤-52℃的频率Fs高值中心(区)关系密切,由FM和Fs可以相当明显地反映降水的总体特征,而且比使用TBB低值中心(区)来反映降水的总体特征能力要强;地面的梅雨锋、低层的切变线、对流层上层的高压这种上下层系统的配置十分有利于强对流天气的发生.     

半日潮特征海域平均大潮高潮面的计算方法

计算半日潮特征海域平均大潮高潮面的关键在于,一是应有足够长时段且满足精度的潮汐实测或预报数据;二是应尽可能准确选定大潮时刻高潮位。针对大潮时刻高潮位在时间选取上的困难,从月相物理意义出发,证明了当日月黄经差值约为0°时,大潮时刻对应月相朔;当日月黄经差值约为180°时,大潮时刻对应月相望。提出并实现了顾及半日潮龄值及逐时日、月黄经差的平均大潮高潮面计算方法。     

刚性桩复合地基中性面深度与桩土应力比计算

刚性桩复合地基中性面深度及桩土应力比的简化计算主要基于桩侧摩阻力线性分布假设,当桩身较长时,桩端侧摩阻力的计算值会远大于实际值,致使中性面深度及桩土应力比的计算结果与实际差别较大,故有必要对线性分布模式予以修正。据此将桩侧摩阻力分布简化为分段线性模式,考虑负摩阻力作用及桩上、下刺入变形,根据褥垫层-桩-土变形协调关系推导了刚性桩复合地基中性面深度、桩顶面桩土应力比、中性面桩土应力比计算公式。最后通过模型试验与工程实例验证,计算值与实测值吻合较好。