首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到19条相似文献,搜索用时 109 毫秒
1.
井间高分辨率纵横波层析成像研究井间油藏   总被引:3,自引:1,他引:2       下载免费PDF全文
瞿辰  杨文采  于常青 《地球物理学报》2010,53(12):2944-2954
采用视速度-偏振井间波场分离法分离得到纵波和横波波场,使用跨孔地震走时层析成像级联方法计算井间纵波和横波速度结构以及泊松比空间分布.综合纵横波速度、泊松比和测井参数,分析两井间岩性和储层特性.数值模型试验表明该方法垂直分辨率达到主频波长的1/4,并能准确分辨倾斜断层和垂直断层.处理相同纵波资料结果表明本文方法分辨率远远高于以往的井间走时层析成像方法,比井间声波层析成像分辨率更高.使用该方法处理垦71区井间资料,得到同台同源纵波和横波速度结构和泊松比分布.泊松比的低/高变化与测井参数指示的砂/泥岩层基本吻合.根据泊松比分布,参考岩石地震物理学测试数据,可以区分井间介质中的泥岩层、砂岩层,以及砂岩层中的饱水和饱油区,确定储层的连通性,圈定井间的剩余油藏.  相似文献   

2.
井间地震资料精细解释方法研究与应用   总被引:1,自引:0,他引:1       下载免费PDF全文
井间地震资料分辨率很高,但是缺乏成熟的解释方法,本文针对井间地震资料的特点,借鉴地面地震解释技术对井间地震资料开展了精细解释方法研究.首先在深度域资料特点分析的基础上,利用区域速度和层析成像速度将深度域井间地震资料转换到时间域进行研究,并通过粗化降频建立起井间地震资料与地面地震资料之间的相似性联系,得到对井间地震低频反射特征的认识,结合井旁道合成地震记录标定建立起井间地震资料的反射特征,准确的解释了井间地震的反射层位.采用地震瞬时属性分析和时频属性分析得到了对沉积环境的认识,进一步通过波阻抗反演更好的认识井间储层特征.在实际应用中综合多种方法进行解释,并结合测井资料对比落实了井间砂体展布特征,得到对井间储层的精细描述结果,解决了两口井之间储层分布范围、横向连通性及厚度变化等问题,取得了较好的应用效果.  相似文献   

3.
在井间地震初至走时层析成像中,随着相邻地质体速度差的增大,使得射线分布不均匀,以及网格剖分不合适,导致层析成像结果不理想.物理和数值模型的井间走时层析成像表明:当速度差超过300%时,层析结果畸变较大;在300%~150%之间时,层析结果较好;低于150%时,层析结果好.在此基础上,提出了井间多尺度初至走时层析成像方法,即对同一模型采用多种网格剖分来同时进行层析成像,以获得研究区域的速度图像.数学和物理模型的井间多尺度走时层析结果表明:该方法很好地兼顾了层析成像的分辨率和精度,极大地改善了井间地震层析成像的质量.即使速度差超过30%,其多尺度的层析结果仍然较好.因此,这种方法具有实际应用的潜力.  相似文献   

4.
张超  姚华建  童平  刘沁雅  雷霆 《地球物理学报》1954,63(11):4065-4079
伴随层析成像(Adjoint Tomography)通过求解全波方程来准确模拟地震波在复杂介质中的传播,并利用波形信息来反演地下结构,是新一代的高分辨率成像方法.其中3-D伴随层析成像需要庞大的计算资源,而2-D反演相对更具计算效率.面波和远震体波是研究地壳上地幔速度结构的重要方法,它们对S波速度及Moho面的敏感度不同,通过联合反演,可以得到更为准确的S波速度结构及Moho面.通过两种数据的高度互补性,本文提出基于伴随方法的线性台阵背景噪声面波和远震体波联合成像方法,同时约束台阵下方S波速度结构及Moho面形态.我们将该方法应用到符合华北克拉通岩石圈典型结构特征的理论模型上,测试结果表明联合反演方法优势明显,相比于面波伴随层析成像,能获得更高分辨率的S波速度结构,同时能精准约束Moho面形态.相比于体波伴随层析成像,联合反演能有效压制高频假象,降低波形反演过程中的非线性化程度.本研究有望提供一种更为高效精准的线性台阵成像方法,搭建联合伴随层析成像理论框架,提升岩石圈成像分辨率,并为后续其他类型波形数据的引入提供思路和方法.  相似文献   

5.
井间地震层析成像的现状与进展   总被引:16,自引:5,他引:16  
综述了井间地震层析成像研究的现状,给出了小波变换域井间地震层析成像方法的最新进展,并对井间地震层析成像研究给予展望。  相似文献   

6.
地面地震、VSP和井间地震等多尺度地震资料在油气勘探开发中发挥着重要的作用,但不同的观测尺度和资料采集处理方式限制了多种资料之间的有效匹配和综合应用.通过对同一地质体模型采用波动方程开展正演模拟分析,进一步开展逆时偏移有利于更好地建立起多尺度地震资料之间的联系与差异.分析表明,地面地震、VSP和井间地震等不同尺度下的观测资料都能在对应的范围内准确反映地质体的地震反射特征,而观测尺度和观测方式的不同则导致多尺度地震资料在分辨率、波场信息丰富程度和观测范围等方面存在差异;正演模拟地震道的时频分析进一步表明多尺度地震资料在时间域和频率域都存在着较好的联系,特别是井间地震和VSP等井中地震资料能够提供更为合理、可靠的高频信息,在提高地面地震资料分辨率方面具有明显的优势,为开展多尺度地震资料的综合研究与应用提供了基础.  相似文献   

7.
使用广义最小残量方法迭代求解三维频率-空间域声波方程,反演时使用多尺度、多重网格的策略,探讨了如何快速实现高分辨率的三维频率-空间域迭代法声波全波形速度反演.通过对理论模型进行三维频率-空间域迭代法声波全波形反演数值试验,证实该方法的计算速度快、计算效率高,反演所得速度的分辨率高.从而为基于多尺度迭代求解三维频率-空间域声波方程的全波形速度反演成像打下方法基础.  相似文献   

8.
地面电磁法在实际工作中在很大程度上受到了探测深度和分辨率的限制.为了克服这些缺点,出现了井间电磁法.井间电磁法指的是在两个(或多个)钻孔中分别发射或接收电磁波信号,利用电磁波信号进行成像并探测井间物理性质的地球物理方法.由于发射机和接收机可以分别放置在很深钻孔中,其具有大透距、大探测深度的特点,因而广泛应用于工程环境物探、矿产勘查、石油勘探等中.针对不同的应用,产生了各具特点的一些特殊方法,包括井间无线电波成像、跨孔雷达、井间电磁成像.井间无线电波成像仪目前只测量电场强度数据,工作频率低,一般是单频的电磁波,频率范围通常在1 kHz至10 MHz.由于缺少走时数据修正射线路径,井间无线电波成像主要是进行基于直射线追踪的衰减层析成像.井间无线电波成像既可用于工程与环境地球物理也可用于找矿.跨孔雷达是钻孔雷达的一种探测方式,用高频电磁脉冲探测两个井孔间介电常数和电导率的变化.跨孔雷达层析成像也叫地质雷达CT,既可进行走时成像,还可进行衰减成像.一般来说,地质雷达CT的电磁波工作频率较高,中心频率通常在10 MHz和1 GHz之间,因此在分辨率指标上占有优势,跨孔雷达主要用于工程与环境地球物理.井间电磁成像采用更低的频率,测量复电磁信号,适合油气储集层监测,是一种地球物理前沿技术.经过在多个试验区初步试验表明,井间电磁成像是油藏研究的有效手段,可用于分析剩余油分布,寻找油气富集区,进而达到提高钻探高效井成功率和提高采收率的目的.本文详细介绍对比了这三种方法在理论和实践中的一些特点,并对未来的发展进行了展望.  相似文献   

9.
基于非规则网格声波正演的时间域全波形反演   总被引:2,自引:2,他引:0       下载免费PDF全文
全波形反演是地震资料处理中速度建模的有力工具,相比层析成像等速度建模方法它能够得到速度场的更高频成分.本文给出了基于声波方程格子法正演的时间域全波形反演方法,该方法用非规则、非结构化的三角网格来离散计算区域及模型参数,能实现网格粒度与反演分辨率在空间上的自动匹配,内存需求少,计算效率高;采用L-BFGS优化方法,以分频段变网格的方式实施多尺度反演.以二维Overthrust模型进行了速度反演数值测试,显示了该方法的高效性和潜力.  相似文献   

10.
利用井间地震层析成像方法探测隐伏断层   总被引:1,自引:0,他引:1  
针对目前隐伏断层探测中存在的技术难题,提出了在地面反射地震勘探初步确定的断点附近,用井间地震层析成像方法探测隐伏断层及其上断点的埋深。首先介绍井间地震技术的应用现状和井间地震层析成像的基本原理,并结合实例讨论了井间初至走时层析成像的野外数据采集和室内资料处理方法。芦花台隐伏断层的井间地震层析成像结果表明,采用弯曲射线追踪方法反演得到的跨芦花台隐伏断层的井间纵波速度结构图像与钻孔联合剖面揭示的断层特征相吻合。   相似文献   

11.
基于对数目标函数的跨孔雷达频域波形反演   总被引:2,自引:1,他引:1       下载免费PDF全文
波形反演在探地雷达领域的应用已有十余年历史,但绝大部分算例属于时间域波形反演.频率域波形反演由于能够灵活地选择迭代频率并可以使用不同类型的目标函数,因而更加多样化.本文的频率域波形反演基于时间域有限差分(FDTD)法,采用对数目标函数,可在每一次迭代过程中同时或者单独反演介电常数和电导率.文中详细推导了频率域波形反演的理论公式,给出对数目标函数下的梯度表达式,并使用离散傅氏变换(DFT)实现数据的时频变换,能够有效地减少大模型反演的内存需求.在后向残场源的时频域转换过程中,提出仅使用以当前频点为中心的一个窄带数据,可以消除高频无用信号的干扰,获得可靠的反演结果.为加速收敛,采用每迭代十次则反演频率跳跃一定频带宽度的反演策略.实验证明适当的频率跳跃能够在不降低分辨率的基础上有效地提高反演效率.通过两组不同情形下合成数据反演的分析对比,证明基于对数目标函数的波形反演结果准确可靠.最后,将该方法应用到一组实际数据,得到较好的反演结果.  相似文献   

12.
全球地磁感应测深数据三维反演   总被引:2,自引:0,他引:2       下载免费PDF全文
全球地磁感应测深能获得地幔转换带及下地幔上部的导电结构.但目前稀疏的地磁台站分布及部分台站的观测数据稳定性较差,影响了三维反演对地下电性结构的分辨力和反演可靠性.为此,区别于传统的L2-范数反演方法,本文提出并实现了基于L1-范数的地磁测深响应三维反演技术.在反演中,利用L1-范数度量数据预测误差,降低"飞点"数据的影响,将相关系数较小的C-响应估计也纳入反演数据中.三维正演模拟采用球坐标系下的交错网格有限差分法,反演采用有限内存拟牛顿法.文中利用指数概率密度分布函数构造非高斯噪声的合成数据,并采用棋盘模型对反演方法的可靠性进行了验证.之后,我们将本文提出的三维反演方法用于全球129个地磁观测台站的C-响应数据反演,结果表明在地幔转换带深部,中国东北地区为高导电异常,南欧和北非则均为高阻;夏威夷在900km以下为高导;菲律宾海及以东地区在转换带表现为明显的高阻,这些结果与前人研究结果一致.由于采用了更多的台站数据,我们的反演结果还发现一些新的异常:南美洲南端,转换带表现为明显的高导;澳大利亚东南部,地幔转换带深部,也存在一个明显的高导异常,这些异常分布和地震层析成像的低速区一致.因此,L1-范数三维反演能够充分利用全球C-响应数据信息,提高地磁测深对地球深部电性结构的分辨能力,更好的研究全球地幔电性结构.  相似文献   

13.
In this paper we propose a 3D acoustic full waveform inversion algorithm in the Laplace domain. The partial differential equation for the 3D acoustic wave equation in the Laplace domain is reformulated as a linear system of algebraic equations using the finite element method and the resulting linear system is solved by a preconditioned conjugate gradient method. The numerical solutions obtained by our modelling algorithm are verified through a comparison with the corresponding analytical solutions and the appropriate dispersion analysis. In the Laplace‐domain waveform inversion, the logarithm of the Laplace transformed wavefields mainly contains long‐wavelength information about the underlying velocity model. As a result, the algorithm smoothes a small‐scale structure but roughly identifies large‐scale features within a certain depth determined by the range of offsets and Laplace damping constants employed. Our algorithm thus provides a useful complementary process to time‐ or frequency‐domain waveform inversion, which cannot recover a large‐scale structure when low‐frequency signals are weak or absent. The algorithm is demonstrated on a synthetic example: the SEG/EAGE 3D salt‐dome model. The numerical test is limited to a Laplace‐domain synthetic data set for the inversion. In order to verify the usefulness of the inverted velocity model, we perform the 3D reverse time migration. The migration results show that our inversion results can be used as an initial model for the subsequent high‐resolution waveform inversion. Further studies are needed to perform the inversion using time‐domain synthetic data with noise or real data, thereby investigating robustness to noise.  相似文献   

14.
Elastic waves, such as Rayleigh and mode‐converted waves, together with amplitude versus offset variations, serve as noise in full waveform inversion using the acoustic approximation. Heavy preprocessing must be applied to remove elastic effects to invert land or marine data using the acoustic inversion method in the time or frequency domains. Full waveform inversion using the elastic wave equation should be one alternative; however, multi‐parameter inversion is expensive and sensitive to the starting velocity model. We implement full acoustic waveform inversion of synthetic land and marine data in the Laplace domain with minimum preprocessing (i.e., muting) to remove elastic effects. The damping in the Laplace transform can be thought of as an automatic time windowing. Numerical examples show that Laplace‐domain acoustic inversion can yield correct smooth velocity models even with the noise originating from elastic waves. This offers the opportunity to develop an accurate smooth starting model for subsequent inversion in the frequency domain.  相似文献   

15.
多震源编码技术可以提高全波形反演的计算效率,但同时会引入串扰噪声使反演结果质量降低. 全变分约束可以有效地压制层内噪声并突出模型界面,其与多震源技术的结合,能在大大提高弹性波全波形反演效率的同时提高反演质量. 本文提出了一种高效的动态多震源全波形反演策略,可以在离散串扰噪声的同时保证照明的均匀性. 根据残留串扰噪声的分布特征,构建了与之匹配的基于各向异性全变分约束的弹性波全波形反演方法. 为了减少周期跳跃效应,将基于稀疏约束的低频重构算法应用于弹性波地震记录,给出了利用快速梯度投影算法求解各向异性全变分约束的全波形反演流程. 模型数据测试结果表明本文的方法不仅能有效地抑制多震源方法导致的串扰噪声,同时也能降低观测数据中的噪声对反演结果的影响.  相似文献   

16.
Waveform inversion met severe challenge in retrieving long‐wavelength background structure. We have proposed to use envelope inversion to recover the large‐scale component of the model. Using the large‐scale background recovered by envelope inversion as new starting model, we can get much better result than the conventional full waveform inversion. By comparing the configurations of the misfit functional between the envelope inversion and the conventional waveform inversion, we show that envelope inversion can greatly reduce the local minimum problem. The combination of envelope inversion and waveform inversion can deliver more faithful and accurate final result with almost no extra computation cost compared to the conventional full waveform inversion. We also tested the noise resistance ability of envelope inversion to Gaussian noise and seismic interference noise. The results showed that envelope inversion is insensitive to Gaussian noise and, to a certain extent, insensitive to seismic interference noise. This indicates the robustness of this method and its potential use for noisy data.  相似文献   

17.
In this paper, we focus on the influences of various parameters in the niching genetic algorithm inversion procedure on the results, such as various objective functions, the number of the models in each subpopulation, and the critical separation radius. The frequency–waveform integration (F–K) method is applied to synthesize three-component waveform data with noise in various epicentral distances and azimuths. Our results show that if we use a zero-th-lag cross-correlation function, then we will obtain the model with a faster convergence and a higher precision than other objective functions. The number of models in each subpopulation has a great influence on the rate of convergence and computation time, suggesting that it should be obtained through tests in practical problems. The critical separation radius should be determined carefully because it directly affects the multi-extreme values in the inversion. We also compare the inverted results from full-band waveform data and surface-wave frequency-band (0.02–0.1 Hz) data, and find that the latter is relatively poorer but still has a higher precision, suggesting that surface-wave frequency-band data can also be used to invert for the crustal structure.  相似文献   

18.
傅磊  刘四新 《地球物理学报》2016,59(12):4464-4472
本文提出了一种初至纵波(P波)与瑞雷面波的交叉梯度联合反演策略.通过对初至P波进行全波形反演可以获得近地表P波速度结构;通过对仅含瑞雷面波信息的地震数据转换到频率-波数域进行加窗振幅波形反演(Windowed-Amplitude Waveform Inversion,w-AWI)可获得近地表横波(S波)速度结构.在二者反演的目标函数中均加入P波速度和S波速度的交叉梯度作为正则化约束项,使得在反演过程中P波速度和S波速度相互制约,相互约束,从而实现对地震初至P波与瑞雷面波的联合反演.数值模拟结果表明交叉梯度联合反演可以提高S波速度反演分辨率,而P波速度反演结果并没有得到提高.实际资料的反演结果表明,交叉梯度联合反演能够获得更加可信的近地表速度结构.  相似文献   

19.
欧洋  高文利  李洋  王宇航 《地球物理学报》2019,62(10):3843-3853
为了避免使用不合理初始辐射场强和方向性因子带来的误差,研究了估计辐射参数的井间电磁波层析成像技术.通过时域有限差分法模拟表明,天线长度与波长的比值、钻孔充填情况、钻孔周围介质的物性均会影响偶极天线的初始辐射场强或方向性因子;为此结合已知的分层资料,将它们设为未知参数,并设定初始辐射场强与发射点位置相关,方向性因子随射线角度而变化;采用正则化反演方法,由钻孔资料建立了模型方差目标函数,使得反演结果与钻遇的地质特征保持一致.通过理论模型试验和实例应用分析表明,相对于传统射线层析成像方法,估计辐射参数的正则化层析成像技术有助于提高反演的准确性.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号