基于快速非支配排序遗传算法的VTI介质多分量叠前联合反演

在VTI介质中,由于引入了各向异性参数使得多分量多参数地震反演问题的非线性程度显著增加,因此采用传统的权重加权法构建单目标函数进行反演得到的反演结果往往并不理想.本文以反射率法为基础,结合快速非支配排序遗传算法研究了一种VTI介质的多分量叠前联合反演方法.该方法以反射率法为正演方程,应用互相关原理构建PP波和PSV波的多目标函数,进而采用快速非支配排序遗传算法全局寻优获得VTI介质的厚度、纵横波速度、密度和各向异性参数等多个参数.在正演的过程中,反射率法可以考虑几何扩散、吸收衰减、透射损失、多次波以及纵横波旅行时不匹配等地震波传播效应,更能精确地描述地震波在地下地层中的真实传播情况;在反演的过程中,快速非支配排序遗传算法可以在不引入权重系数的条件下同时优化多个目标函数,获得联合反演问题的Pareto最优解,既不添加权重系数影响又充分利用多分量地震数据.模型测试结果验证了该反演方法的有效性和可行性.     

北京及其邻区Q值分布特征的研究

地震波通过地球介质时,由于实际介质並非完全弹性,因此,地震波的振幅除了因波阵面的几何扩散而减少外,还有因介质内摩擦等因素造成的吸收而使地震波的振幅随距离而减小。研究地震波在地球介质中的衰减特性,对研究震源本身及其周围环境都很重要,归纳起来有以下几方面的意义: 1.用于解决一系列与地震波在实际地球介质中传播有关的地震问题;     

基于拟微分算子的横向各向同性介质黏声波方程

由于地层的各向异性与黏滞性,地震波在传播过程中表现出方向上的各向异性和振幅值的吸收衰减特征.若简单地将地下介质视为各向同性,实际地质结构的地震波响应可能会被曲解,甚至忽略了有效信息.为了得到深部储层的高精度成像,本文将标准线性固体模型(SLS)扩展到各向异性介质中,并结合各向异性拟声波方程,推导出时间空间域的各向异性介质黏声波方程组,该方程中用伪微分算子来表述地震波的衰减特征.在逆时偏移的过程中,为了解决高频不稳定问题,文中通过引入规则化算子来构建稳定的逆时传播算子.模型的正演模拟表明该方程既能准确的描述各向异性介质下的标量波的传播规律,又能体现地层的吸收衰减效应.     

傅里叶有限差分法保幅叠前深度偏移方法

地震数据中饱含有丰富的走时信息和振幅信息. 为解决传统偏移方法中几何扩散和入射角变化引起的振幅误差问题，本文提出了一种实用的波动方程保幅地震偏移方法. 该方法从全声波方程出发进行单程波保幅分解，得到直观、高效率的直接面对地震波传播波场的压力分量进行延拓的保幅偏移单程波方程，进而推导出一个含有6项的傅里叶有限差分法保幅偏移的算子方程；修改边界条件和成像条件，使修改后的边界条件和成像方程中考虑振幅补偿，从而从三方面补偿几何扩散损失和入射角变化对振幅的影响. 脉冲响应测试、单炮记录的数值试验以及SEG/EAGE盐丘模型的叠前偏移结果表明，该方法不但可以使散射能量聚焦、归位，提高成像精度；而且可以输出正确反映地下反射系数的振幅信息，为后续的地震属性分析(如AVO/AVA)提供更真实的地震信息.     

Q值简介

弹性波通过介质传播时,它的振幅随着传播距离或时间而衰减。有多种原因可以造成这种衰减,如波动振幅随着距离的增加而产生的几何扩散,在介质边界处的反射、折射以及介质的非弹性吸收等。所谓介质吸收是指由于介质的内阻使弹性波部分弹性能量耗散而转化为热能的现象。为了描述介质吸收特性的强弱,引入了一个无量纲的因子“Q”,称为介质的品质因数。“Q”值是用来度量介质中振动或波动能量的非弹性衰减率的物理量,是介质所固有的特性。对于     

一种基于正则化策略的稳定衰减补偿逆时偏移

地下介质通常具有黏滞性,地震波在地下介质的传播过程中将不可避免地伴随着振幅衰减和速度频散,进而影响地震成像的准确性和分辨率.基于黏滞介质的衰减补偿逆时偏移能沿地震波的传播路径恢复其所经历的振幅衰减和相位畸变,有效提升成像效果.然而,由于地层对地震波的吸收衰减效应呈指数变化,衰减补偿过程中高低频分量的非同步增长易导致补偿算法数值不稳定.为此,本文提出了一种基于正则化策略的稳定衰减补偿逆时偏移方法.该方法基于解耦的常Q分数阶拉普拉斯算子黏滞声波方程描述地震波在地下介质中的传播效应,将振幅正则化因子引入该方程的时间-波数域的解析解中,以确保在补偿过程中地震波场能稳定延拓.二维、三维合成数据以及实际资料的偏移算例均证实了该方法的可行性和有效性,所提出的方法能有效地处理衰减补偿中的不稳定问题,明显提升地震资料的偏移成像质量.     

各向异性衰减薄层地震响应特征研究

薄储层的叠前地震响应特征研究,特别是针对具有速度各向异性的含流体薄层,对储层描述具有十分重要的意义.文中基于波动方程数值模拟方法,正演得到各向同性弹性、各向同性衰减、速度各向异性、各向异性衰减模型的地震波场,并对比分析了四种模型的纵波(PP)和转换横波(PS)地震反射特征.研究结果表明:在衰减介质背景下,引入各向异性,PP和PS波的反射波振幅较弹性介质均减弱,且衰减因素对薄层振幅的影响强于各向异性.同时,VTI各向异性衰减在单频PS振幅曲线上表现出强差异性,而HTI各向异性衰减则会影响PP和PS波单频振幅曲线的极值点幅值和位置,通过分析单频振幅曲线的极值点振幅和极值点位置对各向异性衰减薄层预测有指导作用,尤其对平时较难分辨的VTI各向异性衰减薄层,单频分析方法的优势更明显.     

角度域叠前时间偏移振幅补偿方法

本文从单程波方程出发,推导出角度域叠前时间偏移成像幅值计算方法.文章首先在时间域实现了速度梯度引起的幅值透射损失补偿算法;借鉴波动方程偏移的反褶积成像条件,对地震波能量随传播距离的几何扩散进行部分补偿;针对观测系统中炮点覆盖引起的角度域非均匀照明提出了相应的振幅补偿方法.三维水平层状理论模型成像结果表明获得的角道集具有较好的AVA特性.     

频率域有限差分数值模拟地震波衰减和震电效应特征

地震波在地下含流体孔隙介质中传播时,会引起中观尺度的"局域流",进而产生地震波震电效应.基于Biot(1941)固结理论的准静态方程,在频率域中采用空间有限差分方法,正演模拟虚岩石物理岩样的地震波衰减和震电效应.与时间域虚岩石物理方法相比,该方法既可以直接求取任一频率下的地震波衰减和电势,便于应用于实际岩样的预测分析,也避免了讨论岩样外表面施加的力源函数表达式及时间剖分稳定性条件等问题.首先利用周期性层状介质模型验证了本文所描述方法的有效性,并进一步求取分析了周期性层状介质两种不同特征单元的渗流电流密度及电势,数值模拟结果表明由中观尺度"局域流"引起的震电效应电势振幅数量级在实验室测量范围之内,随后,分析研究了四种不同高渗介质占比值的地震衰减及震电效应特征.最后,将本文提出的震电效应数值计算方法推广至二维,并求取了二维斑块饱和模型的地震波衰减、速度频散、电势的振幅和相位角数值结果.     

基于Ganley理论含有衰减特性零偏VSP正演研究及修正

含有衰减的VSP资料能够提供丰富的地下岩层和岩性信息,本文根据Ganley理论对其进行了正演计算.引入品质因子Q值,考虑界面发生反射和透射,对震源位于地表及地下某一深度进行讨论,分别计算水平层状介质中的下行波和上行波.与不含有衰减特性的正演记录及其时频分析对比可以看出,高频成分的吸收、主频的降低、较多的能量衰减揭示了地...     

有效压力与饱和度变化的气藏多波时移地震AVO响应

天然气在开发过程中,储层有效压力和含气饱和度均会发生变化,研究有效压力和含气饱和度的变化对地震响应特征的影响,在基于时移地震的剩余气分布预测研究中具有重要意义。天然气和石油的声学性质有着明显的差异,油藏时移地震的研究成果不能直接应用于气藏,因此需要开展气藏的时移地震研究。利用Shapiro模型表征干岩石弹性模量随有效压力的变化,借助Batzle-Wang方程描述流体速度随压力的变化关系,联合Gassmann理论进行流体替代,表征饱和流体岩石速度随含气饱和度的变化,建立了饱和流体岩石速度随有效压力和饱和度变化的岩石物理模型。基于该模型,对不同含气饱和度和不同有效压力下的气藏储层模型进行了多波时移地震叠前振幅变化（AVO）模拟。结果表明多波时移地震AVO技术可以有效地区分有效压力变化和含气饱和度变化,为进一步开展气藏多波时移地震流体监测提供了理论参考依据。      

基于频变AVO反演的页岩储层含气性地震识别技术

结合地震岩石物理技术,研究了叠前频变AVO反演在四川盆地龙马溪组页岩储层含气性识别中的应用.首先,应用Backus平均理论将测井数据粗化为地震尺度储层模型,应用传播矩阵理论进行高精度地震正演及井震标定,分析页岩气储层地震响应特征.其次,基于岩心观测结果,应用Chapman多尺度裂缝理论设计页岩气储层理论模型,研究储层衰减、频散以及对应的地震反射特征.应用该理论模型测试频变AVO反演方法,计算结果表明：对于研究区地层结构和地震数据,区分流体类型的优势频率不是地震子波的主频,还受层间调谐干涉等储层结构因素控制,也进一步说明理论模型测试和标定的重要性.最后,将频变AVO反演技术应用到四川盆地龙马溪组页岩地层,计算得到的频散属性为页岩气储层含气性识别提供依据.     

基于南黄海海相油气勘探的地震采集技术研究

南黄海盆地中部隆起区前第三系油气储层主要为中一古生代海相碳酸盐岩地层,厚度大、分布广,具有良好的油气前景.由于厚度小于1000 m的新近系低速地层直接覆盖在高速碳酸盐岩地层之上,对地震波的向下传播起到了强烈的屏蔽作用,使之深部目的层地震资料品质差,影响了油气勘探的进程.根据区域地质、地球物理特征,进行了基于AV0(Am...     

基于流体替换技术的地震AVO属性气藏识别(英文)

传统上,油藏地球物理工程师是基于测井数据进行流体替换,计算油藏饱和不同流体时的弹性参数,并通过地震正演模拟分析油藏饱和不同流体时的地震响应,从而进行油气藏识别研究。该研究方案为油藏研究提供了重要的弹性参数和地震响应信息,但这些信息仅限于井眼位置。对于实际油藏条件,地下储层参数都是随位置变化而变化的,如孔隙度、泥质含量和油藏厚度等,因此基于传统流体替换方案得到的流体变化地震响应信息对于油气藏识别具有很大的局限性。研究通过设定联系油藏弹性参数与孔隙度、矿物组分等参数的岩石物理模型,并基于三层地质模型,进行地震正演模拟与AVO属性计算。得到油藏孔隙度、泥质含量和储层厚度变化时地震AVO属性,并建立了饱和水储层和含气储层对应AVO属性(包括梯度与截距)之间的定量关系。建立的AVO属性之间的线性关系可以实现基于地震AVO属性直接进行流体替换。最后,应用建立的流体替换前后AVO属性之间线性方程,对模拟地震数据直接进行流体替换,并通过流体替换前后AVO属性交汇图分析实现了气藏识别。     

基于反射系数谱理论的薄层多波AVO

AVO分析是目前地震勘探潜在油气储层的一个重要方法。体散射信息包含了层结构、岩性和孔隙流体信息,对地震勘探非常有用。但是基于 Zoeppritz方程的传统AVO分析之只包含了单层信息。薄层厚度定量解释对构造解释、储层描述和储层横向预测都非常重要。本文阐述的基于频率域弹性传播矩阵反射系数谱方法既考虑了层界面引起的振幅变化(Zoeppritz方程),也考虑了层内传播引起的振幅变化。因此该反射系数谱既包括单一层界面信息,也包括层内体散射信息。该反射系数谱是层厚和频率的连续函数,便于分析频率和层厚对反射系数谱的影响。可分析的薄层厚度可以无限小,直至消失。可分析的频率是任意的和连续的。这是对时间域反射系数做傅里叶变换无法实现的。地震波的传播是复杂的,各种波型是同时存在而且相互转换的,该反射系数谱考虑了各种波型在传播过程中的相互转换以及多次波。与比射线方法比更便于正演薄层多波多分量AVO响应。     

Study on attribute characterization for reservoir dynamic monitoring by seismic

Study on characterizing reservoir parameters dynamic variations by time-lapse seismic attributes is the theoretical basis for effectively distinguishing reservoir parameters variations and conducting time-lapse seismic interpretation,and it is also a key step for time-lapse seismic application in real oil fields. Based on the rock physical model of unconsolidated sandstone,the different effects of oil saturation and effective pressure variations on seismic P-wave and S-wave velocities are calculated and analyzed. Using numerical simulation on decoupled wave equations,the responses of seismic amplitude with different offsets to reservoir oil saturation variations are analyzed,pre-stack time-lapse seismic attributes differences for oil saturation and effective pressure variations of P-P wave and P-S converted wave are calculated,and time-lapse seismic AVO (Amplitude Versus Offset) response rules of P-P wave and P-S converted wave to effective pressure and oil saturation variations are compared. The theoretical modeling study shows that it is feasible to distinguish different reservoir parameters dynamic variations by pre-stack time-lapse seismic information,including pre-stack time-lapse seismic attributes and AVO information,which has great potential in improving time-lapse seismic interpreta-tion precision. It also shows that the time-lapse seismic response mechanism study on objective oil fields is especially important in establishing effective time-lapse seismic data process and interpreta-tion scheme.     

Seismic wavefield modeling based on time-domain symplectic and Fourier finite-difference method

Seismic wavefield modeling is important for improving seismic data processing and interpretation. Calculations of wavefield propagation are sometimes not stable when forward modeling of seismic wave uses large time steps for long times. Based on the Hamiltonian expression of the acoustic wave equation, we propose a structure-preserving method for seismic wavefield modeling by applying the symplectic finite-difference method on time grids and the Fourier finite-difference method on space grids to solve the acoustic wave equation. The proposed method is called the symplectic Fourier finite-difference (symplectic FFD) method, and offers high computational accuracy and improves the computational stability. Using acoustic approximation, we extend the method to anisotropic media. We discuss the calculations in the symplectic FFD method for seismic wavefield modeling of isotropic and anisotropic media, and use the BP salt model and BP TTI model to test the proposed method. The numerical examples suggest that the proposed method can be used in seismic modeling of strongly variable velocities, offering high computational accuracy and low numerical dispersion. The symplectic FFD method overcomes the residual qSV wave of seismic modeling in anisotropic media and maintains the stability of the wavefield propagation for large time steps.     

三维AVO解释技术在洪泽地区岩性油藏勘探中的应用

洪泽地区由于沉积的特点，储层横向变化快，油藏受构造、岩性、油源多因素控制。在对该区三维AVO属性体解释中，利用多元回归方法求取了横波曲线，分岩性和含油气性统计了纵、横波、泊松比参数分布规律，建立了本区的含油砂岩的正演模型，从而降低了AVO解释的多解性。通过井-震结合对四种AVO属性数据体进行了标定，并确定了各属性体应用范围，进而进行了储层和含油气检测。实践表明，该方法能有效地利用AVO属性数据体进行储层预测及油气检测，具有一定的推广价值。