共查询到20条相似文献，搜索用时 171 毫秒
1.
基于反射波动方程,本文提出了一种估计地下反射率分布的地震数据最小二乘偏移方法.高频近似下,非齐次的一次反射波动方程的源项是由反射率与入射波场的时间一阶导数相互作用产生的.根据反射波动方程,利用线性最小二乘反演方法由地震反射数据重建出地下产生反射波的反射源,再结合波场正演计算出的地下入射波场,得到地下反射率分布的估计.在地下反射源的线性最小二乘反演重建中,我们采用迭代求解方法,并以地震波的检波器单向地下照明强度作为最小二乘优化问题中Hessian矩阵的近似. 相似文献
2.
非零偏VSP地震资料是一种多分量资料,处理非零偏VSP资料,弹性波叠前逆时深度偏移技术无疑是最适合的处理技术.本文从二维各向同性介质的弹性波波动方程出发,研究了对非零偏VSP资料进行叠前逆时深度偏移的偏移算法,讨论了逆时传播过程中的边值问题和数值频散问题及其相应的解决方案;采用求解程函方程计算得到地下各点的地震波初至时间作为成像时间,实现了非零偏VSP资料的叠前逆时深度偏移.最后进行了模型试算和非零偏VSP地震资料的试处理,结果表明该方法不受地层倾角限制,较适用于高陡构造地区或介质横向速度变化较大地区的非零偏VSP地震资料处理. 相似文献
3.
4.
地震干涉是近几年勘探地球物理领域研究的一个热点方向,它是波场重构的有力工具.然而,地震干涉往往引入虚假同相轴,影响波场重构的质量.为进一步分析虚假同相轴产生的原因并改善波场重构的质量,本文基于稳相分析,详细探讨子波主频、激发炮数、检波器埋深、检波器间距和地层倾角五个因素对波场重构的影响.模型结果表明,震检排列方式和地层倾角等因素通过改变有限震检范围内稳相点的位置来影响波场重构的质量.将垂直地震剖面VSP中的下行直达波与下行反射波进行地震干涉处理,可有效重构来自高陡反射面的下行反射波,将传统的VSP转化为单井地震剖面SWP.直接利用重构SWP波场进行成像,不仅扩大传统VSP的成像范围,而且避免常规勘探面临的静校正和近地表速度建模问题,为高陡构造成像提供一种新方法. 相似文献
5.
试图运用波场动力学途径,旨在增强观测系统分析评估,本文提出基于单程波方程的角度域照明分析方法.本文方法不仅可用于观测系统设计,对叠前成像资料的考量、叠前振幅补偿和AVA（Amplitude Versus Angle,振幅随角度变化）反演也必具重要的意义.文中通过目的层邻域单个绕射点的正、反传波场分析,给出地下目标的水平和倾斜界面随角度变化的成像照明强度.单程波方法可用于模拟地震波在复杂速度构造下的传播,并得到较准确的幅值和多次到达的波场,用以对绕射点正、反传分析时可在照明分析中简明地利用幅值信息和多次到达的波场.本文方法可显现复杂盖层下目标的照明情况.本文中也对照明能量随传播距离的几何扩散进行补偿,突显了角度域的照明均匀度和范围,从而使得所分析的结果可直接应用于AVA的分析及补偿.针对照明分析的特点,文中建议可用一种快速的单程波波场延拓策略——频率速度相关变步长波场深度延拓.数值实验结果表明了本文方法是可行且有效的. 相似文献
6.
预测多次波的逆散射级数方法与SRME方法及比较 总被引：1，自引：0，他引：1
逆散射级数方法和SRME方法是典型的两类基于波动方程的无需地下结构信息的自由表面多次波预测方法,本文详细讨论了这两类方法的基本思想、方法原理和实现思路,并对2维逆散射级数方法进行了降维简化,推导给出了1.5维逆散射级数衰减自由表面多次波的方法,减少了计算量,并降低了对3维规则观测系统的要求.利用经典SMAART模型数据测试比较了2维逆散射级数方法、1.5维逆散射级数方法及2维SRME方法预测自由表面多次波的效果及优缺点,比较了三者在方法实现、并行计算效率及对数据要求上的异同,并结合分析比较结果,给出了实际应用中方法选择的建议. 相似文献
7.
COMPUTATION OF ZERO-OFFSET VERTICAL SEISMIC PROFILES INCLUDING GEOMETRICAL SPREADING AND ABSORPTION*
Synthetic vertical seismic profiles (VSP) provide a useful tool in the interpretation of VSP data, allowing the interpreter to analyze the propagation of seismic waves in the different layers. A zero-offset VSP modeling program can also be used as part of an inversion program for estimating the parameters in a layered model of the subsurface. Proposed methods for computing synthetic VSP are mostly based on plane waves in a horizontally layered elastic or anelastic medium. In order to compare these synthetic VSP with real data a common method is to scale the data with the spherical spreading factor of the primary reflections. This will in most cases lead to artificial enhancement of multiple reflections. We apply the ray series method to the equations of motion for a linear viscoelastic medium after having done a Fourier transformation with respect to the time variable. This results in a complex eikonal equation which, in general, appears to be difficult to solve. For vertically traveling waves in a horizontally layered viscoelastic medium the solution is easily found to be the integral along the ray of the inverse of the complex propagation velocity. The spherical spreading due to a point source is also complex, and it is equal to the integral along the ray of the complex propagation velocity. Synthetic data examples illustrate the differences between spherical, cylindrical, and plane waves in elastic and viscoelastic layered media. 相似文献
8.
正演模拟技术在地震采集设计中的应用 总被引：1，自引：0，他引：1
随着地震勘探开发的不断深入和发展,地震勘探的主战场逐渐向复杂地区转移.复杂地区既指地表条件复杂的地区,也指地下地质构造和地层条件复杂的地区,这些都对地震勘探提出了新的挑战和更高的要求.地震正演模拟正是开展此类问题研究的一个重要手段和方法.目前市场上具有正演模拟功能的软件大多是根据射线理论采用射线追踪的方法来完成正演模拟的,这种方法不能很好地反映地震波的动力学特征,特别是在复杂地区难以得到正确的结果.本文利用高阶有限差分有效克服了常规有限差分算法求解波动方程的频散问题,并以高效的OpenMP并行计算模式进行了并行优化,较大程度上提高了正演计算的速度和精度;同时实现了二阶Higdon边界条件,改善了边界吸收效果;也在一定程度上提高了计算速度.塔中地区主要目的层埋藏深,逆断层发育,地震反射特征复杂,增加了地震勘探的难度.本文依据该地区的地质模型,利用波动方程正演技术论证了该地区的地震采集观测系统,为该区地震采集观测系统的设计提供了科学依据. 相似文献
9.
在地震勘探中,描述复杂介质的正演和反演问题通常包含许多反映介质不同特性的参数.同时获得这些参数对进行更准确的岩性描述和油藏预测具有重要的理论和现实意义.为了提高频率域黏弹性波动方程的零偏VSP多参数反演的精度,本文对多参数反演的可行性进行分析,明确了目标函数的敏感程度及参数之间的耦合情况,提出了一种基于走时约束的分频分步多参数反演策略.首先利用零偏VSP资料构建先验信息,然后分别利用高、低频数据进行两步反演,也就是"三个参数反演+五个参数反演"的过程,以提高反演的稳健性和精度.利用此方法可同时得到零偏VSP数据可靠的弹性波速度、密度和品质因子,为精确的时-深关系及含油气的解释和预测奠定基础,同时也可以为地面地震叠前反演提供可靠有效的约束,增强地面地震反演精度. 相似文献
10.
11.
偏移成像是VSP数据处理中的一个重要环节,常规的VSP成像方法通常利用VSP-CDP转换或Kirchhoff偏移,均存在保幅性差及成像精度低等问题,而波动方程叠前深度偏移被认为是对地下复杂构造进行成像的精确偏移方法.任意广角波动方程作为一种高精度的空间域单程波波动方程,同时由于只含有二阶偏导数项,易于数值实现,与其他单程波波动方程相比,具有更大的成像倾角,因此是偏移成像的有力工具之一.本文将AWWE推广应用到VSP数据成像中,实现了VSP时空域高角度单程波方程偏移.首先从三维标量任意广角波动方程出发,推导了完全匹配层吸收边界条件,在基本不增加计算量的前提下有效地压制了边界反射成像噪音,同时利用非线性反演算法优选参考速度来提高平方根算子的近似程度,从而提高高角度地层的成像精度.模型数值模拟实验验证了该方法的有效性,同时表明该方法在陡倾角构造情况下能取得很好的成像效果.最后对某地区实际观测的VSP资料进行了偏移成像,并与地面地震偏移结果进行了对比,显示出VSP波动方程偏移在成像分辨率上的优势. 相似文献
12.
受井中检波器串级数局限,垂直地震剖面(VSP)反射波成像范围窄,且不能对井中最浅接收点上方区域有效成像.虽然多次波成像能扩大成像范围,但在实际应用中尚有诸多困难和挑战.本文根据Wapenaar的地震干涉理论,基于上下行波场分离技术,研发了VSP地震干涉成像方法.该方法将VSP自由表面多次波重建为在地表震源位置激发(虚震源)接收的拟地面地震反射波,然后偏移成像,以达到对多次波间接成像的目的.通过数值模型实验,测试了VSP干涉成像的极限分辨率,并讨论了主要采集参数的影响,结果表明:该方法的垂向和水平极限分辨率分别达约10m和20m,且能分辨深度达6500m处的50m×100m溶洞;采用12至24道井中检波器串采集的VSP资料,其干涉成像结果显著优于VSP反射波成像,与相应的地面地震成像效果相当.将本文方法应用于新疆地区采集的VSP资料,结果表明:与VSP反射波成像相比,成像同相轴更加连续,成像范围显著扩大;与地面地震成像相比,成像结果相当,尤其在浅中部甚至更好.新方法不仅无需进行井中接收点静校正,且能显著增大成像范围,有利于成像同相轴的追踪对比、地震属性提取和地质解释,尤其对中国新疆地区深部缝洞型储层的成像,具有广泛的实际应用前景. 相似文献
13.
A deterministic near-surface correction procedure is developed for multicomponent VSP data, shot using directional sources and recorded using three-component receivers. The method is capable of removing unwanted effects of acquisition such as unequal source strengths or misorientations, but may also remove near-surface multiples and anisotropy. This is of considerable benefit for obtaining accurate and consistent estimates of subsurface anisotropy from different source combinations. Application of the technique is illustrated using a dataset from the Romashkino reservoir in Russia, where three or four different directional sources are used at the same source locations. The technique corrects for the large discrepancies which exist between the estimates obtained using different source combinations. Application of the technique to three wells in the survey region reveals a nearly isotropic subsurface, except for a few isolated zones of moderate to high (2 to 8%) anisotropy which lie close to the expected depth for the reservoir. Although there is no significant correlation with the production figures for each well, the qS1 polarization azimuth within the reservoir does vary at each well location, suggesting that this may be a more sensitive indicator of reservoir Drocess. 相似文献
14.
尽管叠前逆时偏移成像精度高,但仅针对单一纵波的成像也可能形成地下介质成像盲区,由于基于弹性波方程的逆时偏移成像可形成多波模式的成像数据,因此弹性波逆时偏移成像可提供更为丰富的地下构造信息.本文依据各向同性介质的一阶速度-应力方程组构建震源和检波点矢量波场,再利用Helmholtz分解提取纯纵波和纯横波波场,使用震源归一化的互相关成像条件获得纯波成像,避免了直接使用坐标分量成像而引起的纵横波串扰问题.针对转换波成像的极性反转问题,文中提出一种共炮域极性校正方法.为有效节约存储成本,也提出一种适用于弹性波逆时偏移的震源波场逆时重建方法,在震源波场正传过程中,仅保存PML边界内若干层的速度分量波场,进而逆时重建出所有分量的震源波场.本文分别对地堑模型和Marmousi2模型进行了弹性波逆时偏移成像测试,结果表明:所提出的共炮域极性校正方法正确有效,基于波场分离的弹性波逆时偏移成像的纯波数据能够对复杂地下构造准确成像. 相似文献
15.
Heather Schijns Suvi Heinonen Douglas R. Schmitt Pekka Heikkinen Ilmo T. Kukkonen 《Geophysical Prospecting》2009,57(6):997-1008
A traveltime inversion technique is applied to model the upper ∼40 m of the subsurface of a glaciated shield rock area in order to calculate static corrections for a multi-azimuth multi-depth walk-away vertical seismic profile and a surface seismic reflection profile. First break information from a seismic refraction survey is used in conjunction with a ray-tracing program and an iterative damped least-squares inversion algorithm to create a two-dimensional model of the subsurface. The layout of the seismic survey required crooked seismic lines and substantial gaps in the source and receiver coverage to be accounted for. Additionally, there is substantial topographical variation and a complex geology consisting of glaciofluvial sediment and glacial till overlying a crystalline bedrock. The resolution and reliability of the models is measured through a parameter perturbation technique, normalized χ2 values, root means square traveltime residuals and comparison to known geology. 相似文献
16.
Due to the complicated geophysical character of tight gas sands in the Sulige gasfield of China, conventional surface seismic has faced great challenges in reservoir delineation. In order to improve this situation, a large‐scale 3D‐3C vertical seismic profiling (VSP) survey (more than 15 000 shots) was conducted simultaneously with 3D‐3C surface seismic data acquisition in this area in 2005. This paper presents a case study on the delineation of tight gas sands by use of multi‐component 3D VSP technology. Two imaging volumes (PP compressional wave; PSv converted wave) were generated with 3D‐3C VSP data processing. By comparison, the dominant frequencies of the 3D VSP images were 10–15 Hz higher than that of surface seismic images. Delineation of the tight gas sands is achieved by using the multi‐component information in the VSP data leading to reduce uncertainties in data interpretation. We performed a routine data interpretation on these images and developed a new attribute titled ‘Centroid Frequency Ratio of PSv and PP Waves’ for indication of the tight gas sands. The results demonstrated that the new attribute was sensitive to this type of reservoir. By combining geologic, drilling and log data, a comprehensive evaluation based on the 3D VSP data was conducted and a new well location for drilling was proposed. The major results in this paper tell us that successful application of 3D‐3C VSP technologies are only accomplished through a synthesis of many disciplines. We need detailed analysis to evaluate each step in planning, acquisition, processing and interpretation to achieve our objectives. High resolution, successful processing of multi‐component information, combination of PP and PSv volumes to extract useful attributes, receiver depth information and offset/ azimuth‐dependent anisotropy in the 3D VSP data are the major accomplishments derived from our attention to detail in the above steps. 相似文献
17.
压制多次波是地震资料处理过程中的重要环节之一,多次波去除的效果直接影响到最后成像的质量．本文研究的逆散射级数法以波动方程和Born级数为基础,构造一个与自由表面多次波相关的子序列来预测并消除多次波,该方法不需要知道地下的速度结构信息,适用于各种复杂的地下情况．在逆散射级数法的实现过程中,地震子波和近道波场是必不可少的信息,本文采用能量最低法则和带限抛物Radon变换法来弥补原始地震资料中所缺失的这些信息．合成资料和实际地震资料的处理结果表明,逆散射级数法能够在去除自由表面多次波的同时保留有效波的信息,对于地下复杂结构的情况同样有效． 相似文献
18.
复杂地表的单程波动方程地震叠前正演 总被引：4，自引：0，他引：4
作者基于数学检波器和等时叠加原理,实现了复杂地表的单程波动方程地震叠前正演模拟。该方法采用虚拟的数学检波器接收地下的反射地震信号,灵活地将接收点布置在地表的任何地方,从而满足地表起伏的要求。此外,根据等时叠加原理, 该方法采用单程波动方程进行波场延拓和成像,计算简单快速。通过复杂正断层的数值模拟,得到了高信噪比的共炮集地震记录,并采用适用于起伏地形的深度偏移方法对该共炮集地震记录进行了叠前深度偏移,较好地实现了地震波的偏移归位,从而证明了这里提出的起伏地表的单程波动方程地震叠前正演方法是正确和有效的。 相似文献
19.
20.
J.W. SATTLEGGER J. ROHDE H. EGBERS G.P. DOHR P.K. STILLER J.A. ECHTERHOFF 《Geophysical Prospecting》1981,29(2):229-240
Two dimensional inverse modeling, a process to be applied after standard processing and interpretation, uses interfaces picked by the user. These interfaces are transformed into an approximate subsurface model. The subsurface model is represented by curved interfaces and interval velocities. The interfaces have to be unique functions of the line coordinate. Otherwise they may be arbitrarily curved and may begin or terminate anywhere along the section, e.g., at faults, pinchouts, salt domes and the like. Interval velocities may vary laterally along the section. The inverse modeling algorithm then modifies the model until traveltimes calculated from this model match the traveltimes observed as closely as possible in a least squares sense. The traveltimes corresponding to the model are obtained through ray tracing taking exact account of refraction. The traveltimes observed are the arrival times of single impulses before stacking contributing to the interfaces. These traveltimes are provided by ANAKON, a continuous interface analysis system. The comparison of INMOD results with those of well measurements and those of classical interval velocity computation from seismic data shows the accuracy of the method. Deviations of INMOD derived interface depths are within 2% of well data. 相似文献