共查询到19条相似文献，搜索用时 218 毫秒
1.
2.
利用地震波正向传播方程对属于波形线性反演问题近似求解方法的地震数据偏移成像进行重新推导,得到了适合散射地震数据的散射偏移成像方法和适合反射地震数据的反射偏移成像方法.以地震波传播的散射理论为出发点,首先根据描述一次散射波正向传播的线性方程研究建立散射地震数据的偏移成像方法理论;利用高频近似对产生散射波场的地下速度扰动函数的空间变化进行近似,推导出地下反射率函数,再由散射波传播方程推导出基于反射率函数的反射波传播方程,然后根据描述一次反射波正向传播的线性方程研究建立反射地震数据的偏移成像方法理论.本文指出和修正了Claerbout偏移成像方法中的不足,提出的地震数据偏移成像方法是对当前偏移成像方法理论的完善,使反射地震数据偏移成像具有了更坚实的数学物理理论基础,得到的偏移成像结果相位正确、位置准确、分辨率提高. 相似文献
3.
4.
页岩的强各向异性特征挑战地震波传播数值模拟方法的精度极限, 特别是易引起频散的高频波(>100 Hz)传播的数值模拟.鉴于目前我国页岩气地震勘探主要以常规地震声波资料为主,本文首先介绍了一种VTI介质声波方程的任意偶数阶有限差分数值模拟方法,并讨论其稳定性条件和吸收边界条件.任意偶数阶的差分解可有效提高计算精度,压制数值频散噪声.针对页岩较强的各向异性特征,本文比较了不同模型的声波方程和VTI介质声波方程计算得到的地震响应.数值结果表明,各向异性对地震波的运动学(相位)和动力学(振幅)特性影响作用明显.因此,在页岩气地震勘探资料处理的各个环节必须充分考虑各向异性的影响,采取有别于常规油气勘探的处理流程和技术. 相似文献
5.
煤层陷落柱是煤田勘探开发中常见的一种典型的非均匀地质体.由于来自陷落柱的反射信号少、反射能量弱,使得基于反射波原理的常规地震成像方法难以有效识别陷落柱.本文以散射波理论为基础,采用数值模拟方法,研究了陷落柱的散射波场特征,研究表明地面接收的波场中含有来自陷落柱陡倾角界面的散射波场.通过共散射点道集波场的模拟,可以清晰地识别散射波,获得地下散射点和非均匀地质体的信息,判断散射点的位置,从而勾画出不均匀地质体的形态.采用等效偏移距假设抽取共散射点道集,在此基础上进行叠前偏移,对陷落柱成像;模拟与实际数据成像结果对比表明此方法能够合理地提取散射点的散射波场信息,对陷落柱形态及内部结构准确成像,是一种有效的煤田陷落柱成像方法. 相似文献
6.
海底冷泉广泛分布于全球大陆边缘海域,是未来能源勘探的重要领域.本文根据海底冷泉羽流特征,设计多个典型2D冷泉气泡散射地球物理模型并观测其数值模拟的地震响应特征.冷泉正演模拟的关键是建立典型地震海底冷泉散射模型,本文首先构建了平稳随机散点模型,设计了符合实际的气泡半径和气泡体积分数的参数模型,通过参数填充和计算得到了等效速度模型.由此建立的模型在适应性和灵活度上凸显出很大优势,极大地提升了地震响应的可靠程度.单炮地震记录显示冷泉羽流区散射能量最强,地震波场散射清晰,散射波场沿冷泉延伸方向具有明显的方向性.复杂冷泉气泡散射模型的正演模拟结果表明地震响应能够准确描述海底冷泉羽流,对冷泉羽流的研究和寻找海洋天然气资源具有重要意义. 相似文献
7.
针对深部隐伏的金属矿床,以铜陵某矿区为例,应用散射波正演模拟技术,对金属矿地震数据野外采集中涉及的偏移距、覆盖次数、道间距等参数进行了不同炮记录的模拟对比,辅助指导非均匀地质体金属矿地震勘探野外数据采集的观测系统设计;针对该地区实际采集的地震数据,应用基于等效偏移距的共散射点成像方法,对金属矿地震数据进行了处理,通过与常规反射波法处理结果对比,散射成像效果较好,说明散射波技术在金属矿地震勘探中应用有一定的可行性和适用性. 相似文献
8.
通过对波场的时间二阶积分运算以增强地震数据中的低频成分,提出了一种可有效减小对初始速度模型依赖性的地震数据全波形反演方法—时间二阶积分波场的全波形反演方法.根据散射理论中的散射波场传播方程,推导出时间二阶积分散射波场的传播方程,再利用一阶Born近似对时间二阶积分散射波场传播方程进行线性化.在时间二阶积分散射波场传播方程的基础上,利用散射波场反演地下散射源分布,再利用波场模拟的方法构建地下入射波场,然后根据时间二阶积分散射波场线性传播方程中散射波场与入射波场、速度扰动间的线性关系,应用类似偏移成像的公式得到速度扰动的估计,以此建立时间二阶积分波场的全波形迭代反演方法.最后把时间二阶积分波场的全波形反演结果作为常规全波形反演的初始模型可有效地减小地震波场全波形反演对初始模型的依赖性.应用于Marmousi模型的全频带合成数据和缺失4Hz以下频谱成分的缺低频合成数据验证所提出的全波形反演方法的正确性和有效性,数值试验显示缺失4Hz以下频谱成分数据的反演结果与全频带数据的反演结果没有明显差异. 相似文献
9.
地震槽波动力学特征物理-数学模拟及应用进展 总被引：1，自引：0，他引：1
复杂的地质构造和生产环境使得中国的煤业生产时常蒙受巨大的经济损失和人员伤亡,而地震槽波井下探测技术做为当前可探明煤层中局部小构造和异常体的有效技术,具备分辨率高和预测性强等特点,是煤矿生产中井下探测的一种有效的地震勘探方法.本文回顾了地震槽波勘探研究的发展概况,并基于地震槽波在煤层中传播的物理特征和影响煤层中地震槽波传播等主要因素,对地震槽波动力学特征方面已经开展的物理模拟和数值模拟工作进行了综合阐述,进而指出了地震槽波勘探存在的理论问题、应用难题和可能的发展方向. 相似文献
10.
时间采样以及空间采样问题是地震勘探基本问题,空间采样方式包括规则采样、非均匀(非规则)采样、随机采样等.从地震波场的带限性、空间域非均匀性、地下构造多样性出发,提出多采样率地震勘探技术(MrSET),尝试组合不同空间采样方式.多采样率地震勘探观测系统设计一般应满足5项准则,探讨了多采样率数据处理基本方法,并以合成数据和实际数据的处理效果为基础,总结了多采样率地震勘探技术的7个方面优势.在相同道密度前提下,与规则采样相比,多采样率地震勘探技术更能满足复杂地震波场、弱能量绕射波与散射波充分采样的需要,是一项能提高反射地震数据成像精度的地震勘探技术,经过进一步研究将有巨大潜力成为(陆上、海洋)节点地震勘探时代常规技术. 相似文献
11.
一般而言,由于地下非均质体的存在所产生的二次波源,由它再生成新的波场,叫散射波场。目前,将散射波作为有效波来成像,已开始在溶洞和裂缝等特殊地质体的识别中得到应用。但对于野外资料采集来说,地表复杂地区,如戈壁、砾石区和山前带,大量存在的散射波却是干扰波,它们的存在会严重影响资料的品质,而其研究与实际应用国内外还很少。因此,通过正演模拟,分析散射波的基本特征,在此基础上研究散射波的去噪方法显得十分必要。本文从地震波运动学时距关系出发,研究了反射波和散射波的几何特征;然后用有限差分正演,模拟了散射波场,用理论模型研究并测试了局域双曲线Radon变换散射波去噪新方法。对于实际炮集资料,分析了F-K滤波方法压制散射噪声的局限,采用局域双曲Radon变换有效地去除了炮集中存在的散射噪声,取得了较好的应用效果。 相似文献
12.
基于精细积分法的三维弹性波数值模拟(英文) 总被引：1，自引：0，他引：1
波动方程有限差分法是地震数值模拟中的一种重要的方法,对理解和分析地震传播规律、分析地震属性和解释地震资料有着非常重要的意义。但是有限差分法由于其离散化的思想,产生了不稳定性。精细积分法在有限差分法的基础上,在时间域采用解析解的表达形式,在空间域保留任意差分格式,发展成为半解析的数值方法。本文结合并发展了以往学者的成果,推导了任意精细积分法的三维弹性波正演模拟计算公式,并对其稳定性进行了数值分析。在计算实例中,实现了精细积分法二维和三维弹性波模型的地震正演模拟,对计算结果的分析表明,精细积分法反射信号走时准确,稳定性好,弹性波场相较于声波波场,弹性波波场成分更为丰富,包含了更多波型成分(PP-和PS-反射波、透射波和绕射波),这对实际地震资料的解释和储层分析有重要的意义。实践证明,该方法可直接应用到弹性波的地质模型的数值模拟中。 相似文献
13.
Diffraction, once considered noise rather than signal in seismic tomography, has recently been successfully used for delineating smaller objects whose size is comparable to or less than the wavelength of the incident seismic wave. But Most of the seismic diffraction tomography works are based on Born's and Rytov's approximation of weak scattering to linearize the expression for the scattered field due to an object and this is the basis of Fourier Diffraction Theorem. However, In this paper, we used the amplitude of the total scattered wave field without considering any approximation of weak scattering to obtain the synthetic scattered field data for forward modeling. These synthetic data are then inversed by Very Fast Simulated Annealing method which is a nonlinear inversion technique to obtain the unknown model parameters. This algorithm is applicable to both weak as well as strong scatterers. The simulated results give fairly good estimate for unknown model parameters of the scattering object within acceptable error limits. 相似文献
14.
15.
Reverse-time migration (RTM) is based on seismic numerical modeling algorithms, and the accuracy and efficiency of RTM strongly depend on the algorithm used for numerical solution of wave equations. Finite-difference (FD) methods have been widely used to solve the wave equation in seismic numerical modeling and RTM. In this paper, we derive a series of time–space domain staggered-grid FD coefficients for acoustic vertical transversely isotropic (VTI) equations, and adopt these difference coefficients to solve the equations, then analyze the numerical dispersion and stability, and compare the time–space domain staggered-grid FD method with the conventional method. The numerical analysis results demonstrate that the time–space domain staggered-grid FD method has greater accuracy and better stability than the conventional method under the same discretizations. Moreover, we implement the pre-stack acoustic VTI RTM by the conventional and time–space domain high-order staggered-grid FD methods, respectively. The migration results reveal that the time–space domain staggered-grid FD method can provide clearer and more accurate image with little influence on computational efficiency, and the new FD method can adopt a larger time step to reduce the computation time and preserve the imaging accuracy as well in RTM. Meanwhile, when considering the anisotropy in RTM for the VTI model, the imaging quality of the acoustic VTI RTM is better than that of the acoustic isotropic RTM. 相似文献
16.
二维频率空间域的数值模拟方法具有以下的优势:多炮模拟时,计算成本比时间域方法低;无累计误差;在地震反演中处理多震源模拟时,只需要有限的几个频率就可以得到好的反演结果.差分离散化形成的稀疏系数矩阵,需要求解一个巨大规模的线性方程组,最大瓶颈是需要海量的计算机内存,导致计算量庞大.本文在前人研究的基础上,采用嵌套剖分网格排序法,极大限度减少对计算机内存的需求,从而减少了计算量.针对弹性波数值模拟的特征,提出二维频率空间域弹性波多炮模拟的快速计算流程.数值模拟试验证明使用嵌套剖分排序法的弹性波多炮数值模拟比压缩存储法具有节省存储量、计算效率高等优势,为后续的二维频率空间域弹性波全波形反演奠定了很好的基础. 相似文献
17.
We propose a new method for removing sea-surface multiples from marine seismic reflection data in which, in essence, the reflection response of the earth, referred to a plane just above the sea-floor, is computed as the ratio of the plane-wave components of the upgoing wave and the downgoing wave. Using source measurements of the wavefield made during data acquisition, three problems associated with earlier work are solved: (i) the method accommodates source arrays, rather than point sources; (ii) the incident field is removed without simultaneously removing part of the scattered field; and (iii) the minimum-energy criterion to find a wavelet is eliminated. Pressure measurements are made in a horizontal plane in the water. The source can be a conventional array of airguns, but must have both in-line and cross-line symmetry, and its wavefield must be measured and be repeatable from shot to shot. The problem is formulated for multiple shots in a two-dimensional configuration for each receiver, and for multiple receivers in a two-dimensional configuration for each shot. The scattered field is obtained from the measurements by subtracting the incident field, known from measurements at the source. The scattered field response to a single incident plane wave at a single receiver is obtained by transforming the common-receiver gather to the frequency–wavenumber domain, and a single component of this response is obtained by Fourier transforming over all receiver coordinates. Each scattered field component is separated into an upgoing wave and a downgoing wave using the zero-pressure condition at the water-surface. The upgoing wave may then be expressed as a reflection coefficient multiplied by the incident downgoing wave plus a sum of scattered downgoing plane waves, each multiplied by the corresponding reflection coefficient. Keeping the upgoing scattered wave fixed, and using all possible incident plane waves for a given frequency, yields a set of linear simultaneous equations for the reflection coefficients which are solved for each plane wave and for each frequency. To create the shot records that would have been measured if the sea-surface had been absent, each reflection coefficient is multiplied by complex amplitude and phase factors, for source and receiver terms, before the five-dimensional Fourier transformation back to the space–time domain. 相似文献
18.
波动方程数值模拟方法是研究地震波场传播的一种重要手段,本文采用交错网格高阶有限差分方法分别对双程声波方程和双程弹性波方程进行了波场数值模拟,并且根据定位原理采用傅立叶有限差分算子进行了单程波方程数值模拟,在分析定位原理的基础上,对其计算过程稍作修改,将延拓到地面的波场直接由每个检波点接收,无需横向叠加过程,得到了单程声波方程共炮记录.基于不同波动方程的数值模拟结果表明,双程波方程结果包含直达波、多次波等干扰波,信噪比低;单程波数值模拟结果只包含了介质分界面的一次反射波,信噪比高,但对于大角度入射波误差较大,并且对于同一个地质模型而言,双程弹性波方程计算速度最慢,双程声波方程次之,单程声波方程计算速度最快.因此对于复杂地质模型,三种模拟方法可以取长补短,综合应用. 相似文献