共查询到18条相似文献，搜索用时 80 毫秒
1.
复杂近地表散射衰减所有的地面观测波场,形成半随机半相干的近地表强散射噪音背景,弥漫整个炮集,淹没深层反射信号,是导致地震资料极低信噪比的主要原因.如何研究和评价近地表散射强度一直是石油勘探未解决的问题,这与起伏地表的粗糙度、近地表速度横向变化和结构倾角分布密切相关.基于前期复杂近地表边界元法波动方程数值模拟研究,本文提出一种复杂近地表散射振幅矩阵方法来分析近地表散射强度.首先对复杂近地表结构进行边界元配置方法离散,根据边界积分方程生成矩阵方程.我们不求解该矩阵方程(涉及海量计算),只是利用矩阵分析技术来解析矩阵方程中的散射振幅系数矩阵,研究复杂近地表结构对不同频率波场的散射强度.该方法利用边界元对近地表结构几何特征的精确表征,研究起伏地表和非规则地质分界面对地震波传播的影响,由基本解及其在边界上的法向导数经过高斯数值积分计算得到的散射振幅系数矩阵,不仅描述了任意两点之间的相互影响,同时还刻画了边界形状特征的影响,为评价不同地质结构的散射强度提供了可能性.作为初步评价手段,我们采用矩阵元素总和与矩阵维数之比作为表征散射振幅系数矩阵散射特征的标量复杂系数,通过理论和实际模型测试,形成了一套行之有效、计算快速的近地表复杂性分析方法. 相似文献
2.
3.
本文提出了一种计算不规则起伏地形中SH波散射的有效方法——局域边界元法.本方法基于传统边界元法,为计算复杂地表散射问题提供了一种更加高效的解决方案.根据地震波满足的边界积分方程中牵引力格林函数的特性,我们将自由边界分解成水平部分和起伏部分.通过公式推导,可将水平部分的位移由起伏部分的位移通过格林函数线性叠加表示,因此只需对起伏部分的位移进行直接求解,从而极大地减少了待求解的未知数个数,显著提高了计算效率.通过与半圆形山谷SH波平面波入射的解析解比较,验证了方法的正确性.数值模型比较显示,局域边界元模拟结果与传统边界元数值解完全吻合,但是大幅提高了计算效率.因此,局域边界元法可以作为模拟不规则地形中地震波散射的有效工具. 相似文献
4.
提出一种可对起伏地表采集的三维地震资料直接进行偏移成像的叠前时间偏移方法和流程.它用两个等效速度描述近地表和上覆层对地震波传播的影响,可对炮、检点不在同一水平面的三维地震资料直接进行叠前时间偏移处理.该方法不对近地表地震波传播做垂直出、入射假定,因此可适应高速层出露等不存在明显低、降速带情况.描述近地表和上覆层的两个等效速度参数可依据偏移道集的同相轴是否平直来确定,避免了确定近地表速度的困难;而对已知近地表速度的情况,则可进一步修正近地表速度,获得更好的成像效果.用三维起伏地表的理论数据和中国东部某工区实际数据验证了所发展方法和处理流程的有效性和实用性. 相似文献
5.
波动方程数值模拟是深入研究地震波传播规律的有效方法.有限差分法因其方法简单、精度高而得到了广泛的应用.但其缺点是不能准确模拟具有复杂几何形态的物性界面.因而当遇到起伏地表或复杂构造时，求解精度低.为了准确模拟起伏地形、复杂构造和复杂介质条件下的地震波场，本文采用有限元法模拟二维声波方程.用三角形单元模拟地形和速度界面；把单元内的场和波速均看作单元上的线性函数，以适应复杂介质压制边角散射；采用吸收边界条件去除来自截断边界上的反射；采用集中质量矩阵和集中阻尼矩阵使得显式时间递推无需对矩阵求逆，提高了计算效率.对模型的计算表明该方法正确有效. 相似文献
6.
模拟地震波传播的三维逐元并行谱元法 总被引：2，自引：0，他引：2
高效地震波场正演模拟对于复杂模型中地震波传播与成像研究至关重要.本文在谱元法原理框架内,对已有逐元谱元法改进,提出一种新的逐元并行谱元法求解三维地震波运动方程,并得到地震波场.逐元并行谱元法的核心思想在于在单元上进行质量矩阵与解向量的乘积运算,并将此运算平均分配至每一个CPU计算核心,此处理有利提升谱元法的并行计算效率.同时,根据Gauss-Lobatto-Legendre(GLL)数值积分点与插值点重合的特点,将稠密单元刚度矩阵的存储转化成单元雅克比矩阵行列式的值及其逆的存储,大幅减少谱元法计算内存开销.此外,在模型边界上利用逐元并行谱元法求解二阶位移形式完美匹配层(PML)吸收边界条件,消除边界截断而引入的虚假反射.通过逐元并行谱元法得到的数值解与解析解对比,以及实际地震波场模拟,数值结果证实了逐元并行谱元法用于地震波场模拟的高效性. 相似文献
7.
地球深部圈层及沉积盆地是一种分区非均匀介质系统,其中不规则地层边界(含起伏地表)对地震波的主要特征有显著影响,而地层的随机非均匀性则主要影响地震波的散射和衰减特征.为了精确刻画不规则地层边界对地震波的反射、透射效应以及非均质体散射引起的地震波衰减效应,全局广义R/T递推传播矩阵法(GGRTM)被提出并逐步发展成为继有限元和有限差分方法之后的另一种复杂介质高精度地震波传播半解析求解方法.在已有的此类方法中,不规则边界均匀地层GGRTM法的优势在于对不规则地层边界的反射和透射效应的准确模拟,而非均质地层薄板化GGRTM法则能准确描述非均质体散射对地震波衰减的影响.本文吸收这两种已有方法的优势,提出了一种考虑非均匀介质、不规则边界的全局广义R/T递推传播矩阵混合方法,并将其用于对边界不规则、层内非均质的复杂模型的二维SH波场模拟.随后在本文方法与边界元法对比研究的基础上讨论了方法的模拟精度.研究结果表明本文提出的混合法是一种解决复杂模型高精度地震模拟的有效方法. 相似文献
8.
产生地震动空间效应的主要原因是地质构造的变化会极大地改变复杂场地条件下到达不同观测点的地震波,因此地震波振幅和相位发生明显变化。本文利用有限元和边界元耦合的方法,通过建立雷克子波入射时的非线性响应模型求解复杂场地条件下的波动积分方程,进行复杂近地表构造中地震波场数值模拟,旨在分析复杂场地条件下入射方向不同的地震波的动力响应规律。研究结果显示近地表会极大地改变地震波传播特征,导致地震波传播趋向不一致性的因素是由于复杂场地条件重建了地震波的传播路径,尤其须指出的是非一致性地表位移响应源于软土层中低频子波的传播。在地震响应分析中导致系统不同动力响应的其他原因是主频和地震波速度结构的改变,同时地震波入射角度的改变会引起到达观测点的时间滞后,地表位移响应随着入射角度的增大而增加,地震动的空间效应随着入射角度的减小而愈趋明显。 相似文献
9.
10.
在地震波正反演研究中,考虑起伏地表和地震各向异性具有非常重要的理论意义和实际应用价值.本文在前人研究的基础上,将最短路径追踪算法引入到起伏地表各向异性介质模型的地震波走时计算中.模型剖分时,整体模型划分成正方形单元,起伏边界附近以不规则网格逼近,进而采用非规则节点布置实现非规则网格处的最短路径计算.追踪计算中采用Sena群速度近似公式,得到各向异性地震波的走时,实现了复杂地表情况下各向异性介质模型中地震波的射线追踪.理论模型计算结果显示,本文方法能够可靠地应用于复杂各向异性介质模型,具有较高的计算精度. 相似文献
11.
Michel Bouchon 《Pure and Applied Geophysics》1996,148(1-2):3-20
We review the application of the discrete wave number method to problems of scattering of seismic waves formulated in terms of boundary integral equation and boundary element methods. The approach is based on the representation of the diffracting surfaces and interfaces of the medium by surface distributions of sources or by boundary source elements, the radiation from which is equivalent to the scattered wave field produced by the diffracting boundaries. The Green's functions are evaluated by the discrete wave number method, and the boundary conditions yield a linear system of equations. The inversion of this system allows the calculation of the full wave field in the medium. We investigate the accuracy of the method and we present applications to the simulation of surface seismic surveys, to the diffraction of elastic waves by fractures, to regional crustal wave propagation and to topographic scattering. 相似文献
12.
矩形网格有限差分法在地震波传播数值模拟方面具有计算速度快的显著优势,但该方法在处理复杂边界问题上存在着效率低的严重缺陷.本文针对分段光滑曲线边界定义了尖点处的一种正则导数,给出了矩形网格情形分段光滑曲线网格边界点法向导数的一种插值计算方法.采用矩形网格有限差分法对复杂边界地球介质模型进行地震波场数值模拟,并采用波场系列快照技术揭示地震波在起伏地表和复杂介质中的传播规律.模拟结果表明：法向导数插值计算方法为矩形网格有限差分法处理复杂边界提供了有效途径,采用波场系列快照技术可以清晰地展现地震波在反射界面的反射和透射规律、在尖点的绕射规律以及在自由表面的直达波和多次反射规律. 相似文献
13.
The complexity of near surface intensifies the diversity of seismic wave fields, which makes study on near surface wavefields important in many aspects. The strong absorption of low velocity layer can affect the resolution of seismic data, and free boundary can cause surface wave. Considering the above problems, we focus on the Rayleigh wavefields simulation using finite-difference wave equation of higher-order staggered grids and PML boundary conditions. Free boundary, buried source and overlying low velocity layer are taken into consideration and point explosion source is adopted. Through some numerical simulation with different parameters, we quantitatively analyze relationship between wave intensity and source depth, as well as the energy variation with propagation and obtain some practical knowledge and conclusions. 相似文献
14.
15.
The boundary-volume integral equation numerical technique can be a powerful tool for piecewise heterogeneous media, but it is limited to small problems or low frequencies because of great computational cost. Therefore, a restarted GMRES method is applied to solve large-scale boundary-volume scattering problems in this paper to overcome the computational barrier. The iterative method is firstly applied to responses of dimensionless frequency to a semicircular alluvial valley filled with sediments, compared with the standard Gaussian elimination method. Then the method is tested by a heterogeneous multilayered model to show its applicability. Numerical experiments indicate that the preconditioned GMRES method can significantly improve computational efficiency especially for large Earth models and high frequencies, but with a faster convergence for the left diagonal preconditioning. 相似文献
16.
地震自发破裂模拟是震源动力学研究的重要内容，了解复杂的断层动力学破裂过程对深入认识震源特征和解释运动学反演结果具有重要意义.基于边界积分方程方法的破裂模拟已经被广泛使用，大多采用的是平面断层模型的结构化网格划分.由于实际的断层往往具有较为复杂的几何特征，为了更为灵活地刻画断层几何复杂性，我们建立断层模型的三角形网格离散方案，通过精确的解析解形式来计算断层各个单元之间的应力格林函数，联立滑动弱化摩擦准则和非奇异边界积分方程，对断层的自发破裂过程进行了模拟.在简单的平面断层模型下，将计算结果与前人的结果进行了对比，验证了方法的正确性与有效性.对于几种常见的复杂断层模型，例如弯折、阶跃、含障碍体断层等，我们模拟了其破裂过程并对计算结果进行了比较与分析.模拟结果表明，非结构化网格划分的边界积分方程方法能够很好地模拟平面矩形断层或由其组成的规则断层，同时也能成功地模拟具有复杂几何形状的不规则断层上的动力学破裂过程.本研究的结果显示了边界积分方程方法在模拟复杂断层系统的动力学破裂问题上具有较广阔的应用前景. 相似文献
17.
间断有限元方法(Discontinuous Galerkin method,简称DGM)在求解地震波动方程时具有低数值频散、网格剖分灵活等优点,因此,为适应数值模拟对模拟精度和复杂地质结构的要求,本文提出一种新的加权Runge-Kutta间断有限元(weighted Runge-Kutta discontinuous Galerkin,简称WRKDG)方法,用于求解三维D′Alembert介质中声波方程.本文不仅详细推导了其数值格式,特别地,根据常微分方程理论给出了满足数值稳定性条件的一般经验公式,并首次对该方法的数值频散和耗散进行了深入分析,且考虑了耗散参数对结果的影响.同时,我们也对该方法进行了精度测试,并分析了3D情形下WRKDG方法的并行加速比,结果表明3D WRKDG方法具有良好的并行性.最后,我们给出了包含均匀模型、非规则几何模型以及非均匀Marmousi模型在内的数值模拟算例.结果表明,该方法不仅计算准确,能与解析解很好地吻合,且能有效模拟包含球体在内的非规则模型及非均匀Marmousi模型中的衰减声波波场.数值模拟实验进一步验证了WRKDG方法在求解三维D′Alembert介质中声波方程时的正确性和有效性,并获得了对这种强衰减介质中波传播特征的规律性新认识. 相似文献