高阶交错网格有限差分弹性波场模拟的精度分析

交错网格波场数值模拟是目前地震正演中广泛使用的方法,为对比分析不同阶数的差分格式下产生的计算效率和精度差异,重新推导了弹性波方程的4种时间4阶、空间2N阶的差分公式及系数,并计算了他们的稳定性条件。利用这4种差分格式进行弹性波场数值模拟,对比分析了波场快照、合成地震记录及CPU时间。结果表明：时间4阶、空间6+6阶精度的交错网格有限差分方法在进行地震波场数值模拟时具有较高的计算精度和计算效率。     

西藏新生代裂谷系成因的探讨

通过区域重力场的三维密度扰动成像和地震层析成像研究,我们确定西藏新生代裂谷系反映为上中地壳低密度带,延深可达到42 km;同时,青藏高原下地壳也有低密度和低地震波速的物质蠕动流。和东非大裂谷等其他裂谷带不同,西藏新生代裂谷通常是多条平行裂谷组成的裂谷系,它的形成机制也是特殊的,可称为陆—陆俯冲型裂谷系。印度板块的陆—陆俯冲造成的正交方向拉张与西藏新生代裂谷系形成密切相关。地震变形空白区和陆—陆俯冲型裂谷带空间分布之间有一定的对应关系,即裂谷带的源头指向地震变形空白区。陆—陆俯冲型裂谷系的发育过程可分为以下四个环节:①陆—陆俯冲造成前沿带的地壳破裂和地震;②地震变形空白区地应力集中;③挤压地应力向俯冲带前方发散并且转化为張应力,造成前沿带正交方向的地壳破裂和地震;④大地震后应力释放,产生的回跳继续使地壳变形,每一次地震都促使裂谷的进一步发育。青藏高原的下地壳物质蠕动流对中上地壳产生的底辟作用,也促进了西藏新生代裂谷系的形成。地壳拆离面前端上中地壳的成倍加厚使温度升高,造成下地壳流向上挤出,从而使上中地壳張裂。特提斯大洋板块俯冲下去的残块在软流圈下沉也使软流圈上涌,也导致下地壳物质蠕动和西藏新生代裂谷系的形成。     

珠江口盆地惠州地区珠江组砂体上倾尖灭的地质-地球物理“逐级预测”方法

砂体上倾尖灭油气藏是岩性油气藏中重要类型之一,砂体上倾尖灭的刻画与描述是上倾尖灭圈闭有效性评价的关键。在以三角洲前缘沉积为主构成的三级层序内,三角洲沉积相关砂体上倾尖灭的地质-地球物理“逐级预测”描述方法对于认识上倾尖灭油气藏分布规律及其油气勘探具有重要意义。该方法由5个部分组成:在三级层序识别的基础上,运用高精度层序地层划分技术建立了四级层序构成的高频层序地层格架;在高频层序约束下,利用地震沉积学技术描述沉积相的平面分布;在地质模型指导下,利用频谱成像技术定性确定尖灭带的分布;利用相控储层反演半定量、定量预测砂体尖灭线;根据不同成因砂体的沉积特点,利用正演模拟技术开展尖灭线的外推预测。     

地基承载力确定的一种方法

根据三轴和直剪试验的特点并以Mohr强度理论为基础建立了Mohr包线的曲线形式。基于弹性理论导出了条形地基在均布荷载作用下塑性区方程,根据方程可以确定临塑荷载及塑性区的范围。与基于直线型Mohr-Coulomb准则建立起来的塑性区方程相比,该方程能反映土体抗剪强度增长随正应力增加逐渐趋缓的特点。算例分析表明：按两种方法计算得出的临塑荷载相差不大,塑性区连成一片时的荷载也基本一致;两种方法得出的塑性区的形状极为相似;随荷载增加,塑性区向下发展的速度比按传统方法得出的结果稍慢。     

某工程地基软弱夹层的补强设计与检测

对于地基主要受力层范围内局部存在软弱夹层的复杂地质条件,通常的处理方法或挖除换填或采用复合地基,针对具体工程探讨了的不同补强处理设计和施工方法的适用性,利用以往的研究成果,提出了按沉降控制进行补强设计的新思路,并针对处理后的深层复合土层不便直观和开挖检验的难题,结合多道瞬态面波法检验补强加固效果,实现了预期的目的,为类似工程实践提供了有益的借鉴。     

不耦合装药爆破对硬岩应力场影响的数值分析

爆破地震勘探石油是一种重要的方法,但爆破地震效应与爆破参数、地质条件等密切相关。采用动力有限元软件ANSYS/LS-DYNA,对柱状炸药与药孔壁之间为空气或其他介质以及空隙间距变化时碳酸盐岩岩石中爆炸应力波的传播规律和爆炸地震波能量的衰减特性进行了数值模拟研究,得到了不耦合装药爆炸时岩石应力、振动速度的衰减规律以及与不耦合系数、间隙介质的关系,分析了不耦合效应对爆炸地震波能量的影响。研究表明,不耦合或耦合不好时会使岩体中爆炸应力波的强度大大降低;耦合状态对岩体应力及速度的衰减系数和衰减指数影响较大;在空隙中注水或灌满泥浆会改善它们的耦合关系,增大下传的爆破能量。所得成果可为我国西南地区优选适合碳酸盐地层地震勘探的激发因素提供技术途径和方法。     

含有部分脱胶的浅埋圆夹杂对SH波的散射

采用复变函数法研究了含有部分脱胶的浅埋圆柱形弹性夹杂对SH波的散射与地震动问题。在圆形弹性夹杂中构造一个满足脱胶部分应力自由的驻波函数,将其展开为含有一个待定系数的Fourier级数。在半空间中,介质应满足脱胶部分应力自由、公共边界处位移和应力连续的边界条件,从而建立起求解该问题的无穷代数方程组。最后,给出了地表位移的数值结果。给果表明,入射波参数、脱胶位置、弹性夹杂参数以及埋深对地表位移都有一定程度的影响。     

ADI＋内推法：3D傍轴波方程有限差分法的精确解

为3D地震偏移傍轴波方程开发了一种精确的有限差分法解。因为偏移问题中的变数是复值而不是实数，其高阶空间导数的虚数部分不可忽略，所以常规偏移使用的变方向隐合(ADI)方案会产生误差。通过下述方法保护3D傍轴波外推器的精度：①保留这些高阶项，使它在常规偏移中不产生视方位各向异常；②向下外推时滤除非物理损耗波。精确解的实施包括两个步骤。     

基于贝叶斯理论的逐次迭代非线性AVA反演方法

常规AVA三参数反演方法均基于横波速度与纵波速度之比γ为常数这一假设条件,且常被近似地取为0.5。然而在许多情况下γ并不为常数,而是在横向与纵向都渐变。若一概假定γ等于0.5,反演出的岩性参数势必要偏离真实值,因此有必要合理地选择γ。笔者基于贝叶斯理论,提出逐次迭代非线性AVA的反演方法。该方法把γ看成横向与纵向都渐变的反演初始背景,通过给定初始模型计算初始背景γ,并采用逐次迭代的策略求解该反演问题,解决了关于γ的选取问题以及由于引入变γ值而带来的非线性问题,提高了AVA三参数反演结果的精确度。     

双复杂介质条件下的反Q滤波偏移延拓算子研究

地下介质的非弹性吸收衰减效应造成地震资料中深层信号频带窄、能量弱、成像困难,因此需要对地层的非弹性吸收衰减进行补偿。在已有的补偿方法基础上,基于单程波方程,采用与品质因子和频率相关的复速度代替介质真实速度,将粘滞传播算子在频率-空间域采用连分式展开,得到相移项和补偿项,针对连分式的误差,在频率-波数域进行补偿校正,并用逐步累加法解决起伏地表问题。该方法提高了成像精度,有效补偿了由于地层非弹性吸收造成的振幅和频率损失。模型试算和实际资料处理表明,该方法具有一定的适应性和实用性。