裂隙弱度参数的波场特征

本文将Schoenberg线性滑移裂隙介质理论应用到HTI(LSD)介质的弹性波波场数值模拟之中.主要针对描述裂隙状态的2个裂隙参数--法向弱度EN和切向弱度ET对波场的影响进行了有限元模拟.模拟结果揭示了EN主要影响快波而ET主要影响慢波的波场特征,通过EN和ET变化对应的裂隙含水程度说明了EN对裂隙中的流体敏感而E...     

城市密集建筑群对沉积谷地地震动放大效应的影响

采用有限元模拟方法建立了建筑群-沉积谷地二维模型,并在土体截断边界上施加粘弹性人工边界,在频域与时域中对比分析此体系和单独沉积谷地的地震反应,观察地震时沉积谷地与建筑群之间的动力相互作用规律。分析表明,沉积谷地中建筑群对谷地本身的地震反应具有显著影响。入射波频率较低时,由于共振效应的存在,在部分区域处建筑群-沉积谷地体系的地表位移响应幅值会大于单独沉积谷地,但随着入射波频率的增加,建筑群的存在又会对地震反应产生明显的减弱效果;建筑群对谷地的影响还与建筑高度和建筑间距有关,且不同位置处的响应也存在很大差异。计算结果可为沉积谷地中设防烈度的设置以及工程抗震设计提供部分理论依据。     

衬砌隧道中移动轴向激励作用下两相多孔介质动力响应

为了研究衬砌隧道中地铁振动在饱和地层中的传播情况,采用移动轴向激励模拟技术,建立了隧道-衬砌-两相介质的动力分析模型;利用波函数展开法、傅里叶变换法等,推导了频域内衬砌隧道移动轴向激励作用下两相多孔介质动力响应的解析解,并给出了两相多孔介质临界速度的经验公式;通过离散快速傅里叶逆变换得到了时-空域内两相介质的动力响应。结果表明：无衬砌隧道中移动的轴向常激励作用下,两相介质临界速度只与介质的剪切模量和密度有关,且数值接近其剪切波速的1.1倍;对于衬砌隧道,介质的临界速度随着衬砌剪切模量的增大而增大,随着衬砌密度的增大而减小,衬砌对振动向介质中的传播有一定的削弱作用,衬砌剪切模量与介质剪切模量相差越大,削弱越明显;动力响应频率越接近激振频率,其幅值越大,所对应的临界速度越小。     

饱和冻土中弹性体波传播特性影响参数研究

基于Leclaire对饱和双相孔隙弹性介质Biot模型的扩展,研究含有两种不同固相组分的三相多孔弹性介质中体波的传播特性。以饱和冻土为例,分析了各相体积分数、颗粒形状,接触参数等因素对波动方程中惯性参数、黏性参数、刚度参数的影响;对该三相介质模型进行了退化,分析了孔隙中只含液态水或固态冰时体波的特性;以饱和冻土为例,通过数值计算,探讨了饱和冻土中体波的相速度和衰减系数与胶结参数、接触参数、频率、饱和度、孔隙率等参数的关系。结果表明：与一般的饱和土不同,饱和冻土中存在5种体波,即3种纵波和2种横波;5种体波均具有弥散性和衰减性,且P1波、S1波弥散性和衰减性远小于P2、P3、S2波;胶结参数、饱和度、孔隙率对5种体波的传播特性影响显著,接触参数对传播特性影响较小。     

非均匀介质地震AVO响应模拟及分析

地震AVO （amplitude versus offset）技术是一项利用振幅信息研究岩性、检测油气的地震勘探技术。常规方法基于Zoeppritz方程计算模型界面处的反射系数,而实际地下非均匀介质中地震反射特征不仅与入射角度、物性差异有关,还与入射波频率、地层厚度、薄互层结构等因素有关。为此,应用传播矩阵理论充分考虑与这些因素有关的调谐干涉等传播效应,针对实际数据计算非均匀地下介质的高精度合成地震记录,对比Zoeppritz方程、Shuey二项近似方程、Shuey三项近似方程以及传播矩阵算法的模拟效果。研究发现：在小角度入射时Zoeppritz方程、Shuey二项近似方程和Shuey三项近似方程的反射振幅和波形基本一致,大角度入射时Zoeppritz方程与Shuey三项近似方程接近;Zoeppritz算法的模拟结果在小角度入射和浅层情况下与传播矩阵算法差别较小,而在大角度入射和深层情况下与传播矩阵差别较大,说明层间多次波的调谐干涉等传播效应不可忽略。     

基于解耦传播的波场分解方法在VTI介质弹性波逆时偏移中的应用

由于弹性波逆时偏移更符合实际情况,而且转换波的成像结果有更高的分辨率,因此弹性波逆时偏移越来越受到重视。然而,弹性波逆时偏移需要多波多分量数据,为了减少成像结果中的串扰假象,在逆时偏移过程中进行P波和S波分解就变得非常必要。结合基于向量的激发振幅成像条件,我们把基于解耦传播的波场分解方法应用到弹性波逆时偏移中,并对比了其在各向同性介质和横向各向同性（VTI）介质中的应用效果。结果说明,该方法可以在各向同性介质中完全分解P波和S波,并保留波场的向量信息。虽然在VTI介质中有较小的分解残余,但是该分解残余不会在逆时偏移结果中产生明显的串扰;因此,这种波场分解方法可以应用于各向同性介质和VTI介质弹性波逆时偏移。该方法是在时间空间域实现的,可以在波场传播过程中直接对P波和S波进行分离,应用方便,计算效率高。与利用Helmholtz分解的弹性波逆时偏移相比,该方法避免了在PS波成像结果中的极性反转问题。在复杂Hess VTI模型的逆时偏移结果中,高速岩体和断层的成像清晰,甚至是两个低速薄夹层也能较好成像;这说明该方法对复杂介质具有较好的适应性。PS波成像结果中的各向异性体成像清晰,说明各向异性介质弹性波逆时偏移可以对传统逆时偏移不能很好成像的构造进行成像。     

龙门山和相邻地域航磁场特征与汶川大地震

本文应用化极、水平及垂向导数、向上延拓、视磁化强度填图及磁性界面反演等方法处理了龙门山及相邻地域最新的航空磁测数据,分析了龙门山及相邻地域的航磁异常展布特征.研究结果表明:1)龙门山造山带与其东、西两侧可划分为三个磁异常区:松潘—甘孜磁异常区、龙门山负磁异常带、四川盆地磁异常区;三个区、带的地壳介质磁性结构存在明显差异.2)根据该区航空磁异常场的分布特征分别研究了,松潘—甘孜地块、龙门山造山带和四川地块的磁场特征.3)除识别前人识别的断层外,还推断鲜水河ES延伸甘洛—雷波北断裂作为四川盆地与滇西的界带.4)航空磁异常,磁性体上、下界面及磁源体深度的空间分布特征与汶川MS8.0大地震及芦山地震发生相关.     

裂缝诱导双相HTI介质地震波场错格伪谱法模拟与波场特征分析

裂缝诱导的双相具有水平对称轴的横向各向同性(HTI)介质模型是由一组平行排列的垂直裂缝嵌入到统计各向同性的流体饱和多孔隙岩石中而组成的,它综合考虑了裂缝型储层岩石的各向异性和孔隙性.高精度的地震波场数值模拟技术是研究该介质中地震波传播规律的主要方法.本文结合错格伪谱法和时间分裂法,求解描述该介质中地震波传播的一阶速度-应力方程.模拟了单层和双层模型中的地震波场,并对其进行了特征分析.研究结果表明:错格伪谱法能有效消除标准网格伪谱法波场模拟结果中出现的数值伪影现象,与时间分裂法结合能够获得稳定的、高精度的模拟结果;裂缝诱导双相HTI介质中的地震波场兼具裂缝各向异性介质和双相介质中传播的地震波的波场特征.     

复杂山地随机介质GMM-ULTI法射线追踪

对复杂山地介质的非均质性以及介质中地震波运动学特征进行深入研究,对于提高复杂山地区域地震勘探的效果有着重要的理论意义和实际价值.为了研究复杂山地非均质性和该介质中地震波的一些运动特性,提出了一种复杂山地随机介质的建模方法和一种新的射线追踪算法.与常规算法相比,复杂山地随机介质的生成方法采用更贴近实际介质特点的梯度介质作为背景介质,并在模型生成过程中加入地形修正步骤;新提出的GMM-ULTI射线追踪算法,充分融合群推进法、迎风思想、走时插值法的优势,采用先计算走时后追踪射线路径的两步策略完成射线追踪.算法分析与计算实例表明:复杂山地随机介质的生成方法能灵活、精细且更贴近实际地刻画复杂山地介质的非均质特点;新射线追踪算法兼顾精度和效率、能无条件稳定且灵活地适应复杂山地随机介质的特点;同时基于对几个模型试算结果的分析也得出了复杂山地随机介质中的地震波的一些传播规律.     

VTI介质纯P波混合法正演模拟及稳定性分析

各向异性介质纯P波方程完全不受横波的干扰,在一定程度上可以减缓由于介质各向异性引起的数值不稳定,本文推导了具有垂直对称轴的横向各向同性(VTI)介质纯P波一阶速度-应力方程.由于纯P波方程存在一个分数形式的伪微分算子,无法直接采用有限差分法求解.针对该问题,本文采用伪谱法和高阶有限差分法联合求解波动方程,重点分析了混合法求解纯P波一阶速度-应力方程的稳定性问题,并给出了混合法求解纯P波方程的稳定性条件.数值模拟结果表明纯P波方程伪谱法和高阶有限差分混合法能够进行复杂介质的正演模拟,在强变速度、变密度的地球介质中仍然具有较好的稳定性.