共查询到18条相似文献，搜索用时 453 毫秒
1.
弹性波逆时偏移不受倾角和偏移孔径的限制,能够实现任意复杂构造的高精度多波成像,是目前最精确的多分量资料偏移成像方法之一.逆时偏移算法的核心是波场延拓,传统波场延拓以水平基准面为边界条件,基于固定采样步长进行规则网格剖分,采用阶梯近似法处理起伏地表和复杂构造界面时会产生台阶散射,严重影响起伏地表复杂构造的成像精度.基于无网格节点模型,定量分析了弹性波模拟中径向基函数有限差分法的频散关系和稳定性条件.基于此,提出一种基于QR径向基函数的高精度有限差分方法,并提出一种优化的起伏地表自适应节点剖分方法,推导了精确的无网格自由边界条件和弹性波无网格混合吸收边界条件,形成了新的基于无网格的起伏地表弹性波数值模拟方法.此外,本文将此无网格径向基函数有限差分方法应用于精确的纵横波场矢量分解公式,实现了起伏地表弹性波逆时偏移成像.通过对高斯山丘模型,起伏凹陷模型和起伏地表Marmousi-2模型进行数值试算,验证了本文方法的有效性和可行性. 相似文献
2.
采用曲线网格有限差分法描述复杂起伏地形(或不规则波阻抗界面)时,波场正演中可以避免因阶梯近似导致的虚假散射,进而波场逆时偏移可对起伏地表模型进行准确成像.文中以弹性波逆时偏移理论为基础,求解一阶速度-应力方程,推导出了弹性波正向传播和逆时传播的曲线网格差分格式,使用完全匹配吸收边界压制边界反射,采用互相关成像条件,实现了起伏层状介质中的波场逆时偏移.三层起伏、尖灭模型,以及起伏地表条件下的部分盐丘模型结果表明:曲线网格有限差分法逆时偏移法是一种高效、准确的逆时偏移法. 相似文献
3.
4.
随着能源和资源勘查开采工作的深入，地形强烈起伏的盆山耦合地区的地震资料处理解释技术正日益成为山地地震勘探面临的重要挑战.逆时偏移方法作为精确的地震偏移成像方法之一，能对地下结构进行高精度成像.逆时偏移的核心是地震波场延拓，由于传统的地震波场延拓技术往往基于水平地表条件，相应的方法在直接处理强地形起伏条件下的地震资料时往往存在一定的精度损失.本文引入一种精度无损的处理起伏边界的模型参数化方法：基于贴体网格的地形"平化"策略发展了与地形有关的地震波波动方程数值模拟方法，采用零延迟归一化互相关成像条件实现了起伏地表条件下的弹性波场逆时偏移成像.对工业界的标准Marmousi模型和盐丘模型进行改造，获得了相应起伏地形条件下的复杂几何模型，开展了起伏地表下的地震偏移成像数值试验.结果表明基于贴体网格"平化"策略的逆时偏移成像方法具有较高的灵活性，可适应不同类型起伏地表采集的地震资料，显示出该方法在地震勘探领域的良好应用前景. 相似文献
5.
尽管叠前逆时偏移成像精度高,但仅针对单一纵波的成像也可能形成地下介质成像盲区,由于基于弹性波方程的逆时偏移成像可形成多波模式的成像数据,因此弹性波逆时偏移成像可提供更为丰富的地下构造信息.本文依据各向同性介质的一阶速度-应力方程组构建震源和检波点矢量波场,再利用Helmholtz分解提取纯纵波和纯横波波场,使用震源归一化的互相关成像条件获得纯波成像,避免了直接使用坐标分量成像而引起的纵横波串扰问题.针对转换波成像的极性反转问题,文中提出一种共炮域极性校正方法.为有效节约存储成本,也提出一种适用于弹性波逆时偏移的震源波场逆时重建方法,在震源波场正传过程中,仅保存PML边界内若干层的速度分量波场,进而逆时重建出所有分量的震源波场.本文分别对地堑模型和Marmousi2模型进行了弹性波逆时偏移成像测试,结果表明:所提出的共炮域极性校正方法正确有效,基于波场分离的弹性波逆时偏移成像的纯波数据能够对复杂地下构造准确成像. 相似文献
6.
7.
8.
9.
地震资料含有各种类型多次波，而传统成像方法仅利用地震一次反射波成像，在地震成像前需将多次波去除.然而，多次波携带了丰富的地下结构信息，多次波偏移能够提供除反射波外的额外地下照明.修改传统逆时偏移方法，用包含一次反射波和多次波的原始记录代替震源子波，将SRME方法预测的表面多次波代替一次反射波作为输入数据，可将表面多次波成像.多次波成像的挑战和困难在于大量串扰噪声的产生，针对表面多次波成像中的成像噪声问题，将最小二乘逆时偏移方法与多次波分阶思想结合起来，发展可控阶数的表面多次波反演成像方法，有望初步实现高精度的表面多次波成像.在消除原始记录中的表面多次波后，通过逆散射级数方法预测得到层间多次波，将层间多次波作为逆时偏移方法的输入数据可将其准确归位到地下反射位置.数值实验表明，多次波成像能够有效地为地下提供额外照明，而可控阶表面多次波最小二乘逆时偏移成像方法几乎完全避免成像噪声. 相似文献
10.
针对鄂尔多斯盆地西缘黄土塬区复杂地表和复杂地下构造导致难以准确成像问题,采用浅层潜水波层析反演（DWT）速度建模技术,同时辅以中深层反射波层析成像技术,形成一套实用叠前深度偏移速度建模方法。首先生成一个基于钻井和解释信息的起始近地表速度,其次利用潜水波层析反演建立初始近地表模型,将其与常规处理获得的中深层速度模型进行匹配拼接,建立起初始的起伏地表全速度模型,然后在此基础上利用基于反射波的网格层析进行中深层速度建模,经过多轮次迭代,最终获得可靠的高精度速度模型。鄂尔多斯盆地西缘MJT工区地震资料的成像处理验证了这一套速度建模技术的有效性,地下构造成像更合理也更精确。 相似文献
11.
逆时偏移成像(RTM)常用来处理复杂速度模型,包括陡倾角及横向速度变化剧烈的模型.与常规偏移成像方法(如Kirchhoff偏移)相比,逆时偏移成像能提供更好的偏移成像结果,近些年逆时偏移成像越来越广泛地应用到勘探地震中,它逐渐成为石油地震勘探中的一种行业标准.电磁波和弹性波在动力学和运动学上存在相似性,故本文开发了基于麦克斯韦方程组的电磁波逆时偏移成像算法,并将其应用到探地雷达数据处理中.时间域有限差分(FDTD)用于模拟电磁波正向和逆向传播过程,互相关成像条件用于获得最终偏移结果.逆时偏移成像算法中,偏移成像结果受初始模型影响较大,而其中决定电磁波传播速度的介电常数的影响尤为重要.本文基于时间域全波形反演(FWI)算法反演获得了更为精确的地下介电常数模型,并将其反演结果作为逆时偏移成像的初始介电常数模型.为了验证此算法的有效性,首先构建了一个复杂地质结构模型,合成了共偏移距及共炮点探地雷达数据,分别应用常规Kirchhoff偏移算法及逆时偏移成像算法进行偏移处理,成像结果显示由逆时偏移成像算法得到的偏移结果与实际模型具有较高的一致性;此外本文在室内沙槽中进行了相关的物理模拟实验,采集了共偏移距及共炮点探地雷达数据,分别应用Kirchhoff和叠前逆时偏移成像算法进行处理,结果表明叠前逆时偏移成像在实际应用中能获得更好的成像效果. 相似文献
12.
Elastic reverse-time migration (RTM) can reflect the underground elastic information more comprehensively than single-component Pwave migration. One of the most important requirements of elastic RTM is to solve wave equations. The imaging accuracy and efficiency of RTM depends heavily on the algorithms used for solving wave equations. In this paper, we propose an efficient staggered-grid finite-difference (SFD) scheme based on a sampling approximation method with adaptive variable difference operator lengths to implement elastic prestack RTM. Numerical dispersion analysis and wavefield extrapolation results show that the sampling approximation SFD scheme has greater accuracy than the conventional Taylor-series expansion SFD scheme. We also test the elastic RTM algorithm on theoretical models and a field data set, respectively. Experiments presented demonstrate that elastic RTM using the proposed SFD scheme can generate better images than that using the Taylor-series expansion SFD scheme, particularly for PS images. FurH. thermore, the application of adaptive variable difference operator lengths can effectively improve the computational efficiency of elastic RTM. 相似文献
13.
The elastic reverse time migration approach based on the vector-wavefield decomposition generally uses the scalar product imaging condition to image the multicomponent seismic data. However, the resulting images contain the crosstalk artefacts and the polarity reversal problems, which are caused by the nonphysical wave modes and the angle-dependent reduction of image amplitudes, respectively. To overcome these two problems, we develop an amplitude-preserving elastic reverse time migration approach based on the vector-decomposed P- and S-wave seismic records. This approach includes two key points. The first is that we employ the vector-decomposed P- and S-wave multicomponent records to independently reconstruct the PP and PS reflection images to mitigate the crosstalk artefacts. The second is that we propose two schemes in addressing the issue of polarity reversal problem in the conventional PP image. One solution is to adopt the angle-dependent equation. Another one is to reconstruct an amplitude-preserving PP image with a separated scalar P-wave particle velocity, which has a clear physical meaning. Numerical examples using two-dimensional and three-dimensional models demonstrate that the proposed elastic reverse time migration approach can provide the images with better amplitude-preserving performance and fewer crosstalk artefacts, compared with the conventional elastic reverse time migration approach based on the scalar product imaging condition. 相似文献
14.
针对基于互相关成像条件的探地雷达（GPR）逆时偏移计算效率低、存储量大及易产生低频假象的不足，本文将激发振幅成像条件应用于GPR逆时偏移成像中.通过在源点电磁波场正向传播过程计算每个网格点的能量密度，并保存最大能量密度的时刻和相应的电磁波场值；在接收点电磁波场逆向传播过程提取每个网格点最大能量密度时刻及对应的电磁波场值，并利用保存的最大能量源点电磁波场及走时做归一化，从而获得了依赖反射系数成像剖面，避免了源点正向传播电磁波场的存储和重建.此外，为了提高电磁波场的模拟精度，采用了基于三角形剖分的时间域有限单元法（FETD）计算电磁波正向和逆向传播过程.最后通过模型试算表明：激发振幅成像条件相比于归一化互相关成像条件，成像结果低频噪声更弱，空间分辨率更高，计算效率提高了近2倍. 相似文献
15.
16.
17.
各向异性介质纯P波方程完全不受横波的干扰,在一定程度上可以减缓由于介质各向异性引起的数值不稳定,本文推导了具有垂直对称轴的横向各向同性(VTI)介质纯P波一阶速度-应力方程.由于纯P波方程存在一个分数形式的伪微分算子,无法直接采用有限差分法求解.针对该问题,本文采用伪谱法和高阶有限差分法联合求解波动方程,重点分析了混合法求解纯P波一阶速度-应力方程的稳定性问题,并给出了混合法求解纯P波方程的稳定性条件.数值模拟结果表明纯P波方程伪谱法和高阶有限差分混合法能够进行复杂介质的正演模拟,在强变速度、变密度的地球介质中仍然具有较好的稳定性. 相似文献