VTI介质井间地震高斯束正演数值模拟

为了研究垂直横向各向同性(VTI)介质井间地震中,地震波的传播方向及波场特征,采用高斯束方法对复杂构造进行波场正演模拟。分析VTI介质中地震波的群相关系,并提出群相数据库思想来解决地震波在VTI介质中传播的方向和速度问题,改善了波场正演模拟的算法。同时高斯束正演方法不仅能够解决复杂构造的盲区问题,且其动力学特征更能反映各向异性介质对地震波传播的影响。较VTI介质频率-空间域弹性波正演模拟方法,正演结果一致且计算效率高,能有效分辨地震波场特征,为地震数据的处理与解释提供了借鉴。     

准P波方程紧致交错网格井间地震波场模拟及边界条件

研究井间地震波场的形成过程以及波场的传播机理、规律,对于指导实际井间地震勘探有着重要的意义.基于具有垂直对称轴的横向各向同性(VTI)介质中的一阶准P波方程,应用具有无条件稳定性质的紧致交错网格隐式差分格式求解该方程.重点研究了紧致交错网格求解该方程的完全匹配层(PML)吸收边界条件,在此基础上实现了VTI介质中一阶准P波方程的井间地震波场的正演模拟.数值算例表明:紧致交错网格能精准模拟VTI介质中准P波的传播过程,得到高精度的正演结果.一阶准P波方程能以足够的精度描述VTI介质中准P波特征.完全匹配层吸收边界能有效地解决人工边界问题,是一种高效的边界吸收算法.     

运用Hartley变换模拟地震波在正交各向异性介质中的传播

在地震传播理论中,地震波场的正演数值模拟一直是研究的热点。而在正演数值模拟方法的研究中,计算精度和计算效率是评价此方法的有效性及优越性的二个关键问题。伪谱法在计算精度与计算效率方面优越性十分明显,常用的有Fourier变换法和Hartley变换法;虽然Hartley变换法在求取导数时较Fourier变换法复杂,但由于Fourier变换法的计算同时涉及复数的实部与虚部,在计算速度和占用内存方面不如Hartley变换法。这里详细地阐述了利用Hartley变换求解正交各向异性介质波动方程的数值方法,并模拟了多种正交各向异性介质模型,对比分析了地震波在各向同性介质与正交各向异性介质中的传播差异。结论表明,该方法能正确、高效、直观地反应地震波在正交各向异性介质中的传播规律。     

VTI介质中弹性波方程正演的一阶混合吸收边界

波动方程正演是研究地震波在地下介质中的传播机理以及提高地震资料采集、处理、解释精度和效果的重要工具,截断边界的处理是波动方程正演的重要研究内容。针对在地震数值模拟中常用的三种(吸收)边界处理方法存在的不足与问题,提出了适用于VTI介质中弹性波动方程正演的一阶混合吸收边界条件,并给出了该方法在过渡区域的线性叠加系数公式。数值模拟实验表明:相同条件下,VTI介质一阶混合吸收边界方法不仅在吸收效果方面优于常规的PML边界条件,而且在存储量和计算效率方面也具有较大优势。     

各向异性介质储层数值模拟及波场特征研究

各向异性介质储层中波场特征的研究,有利于对多波资料的准确解释。在弹性波方程基础上,结合数值模拟技术与AVO技术,对多种介质模型和第二类含气砂岩储层模型进行数值模拟与波场特征描述。以表征介质各向异性程度的Thomsen参数为出发点,研究各参数对地震波场动力学和运动学特征的影响,并对第二类含气砂岩介质中模拟得到的正演记录,以及各向异性参数对波场的影响程度进行分析,进一步加深了对该类储层各向异性特征的认识,这为各向异性介质中的参数反演与地震资料解释提供了有力的参考依据。     

裂隙各向异性介质2.5维弹性波场数值模拟

研究地震波在各向异性介质中的传播规律是勘探地震学领域的一个重要内容。本文首先讨论了２．５维波动方程，其次，在２．５维空间，用有限元数值模拟方法模拟了裂隙各向异性介质波场，结果用波场切片做了显示。最后对模拟结果做了基本分析，这里讨论的模拟方法也适用于其它类型的得各向异性介。     

裂隙含流体、气体各向异性介质波场数值模拟

由目前普遍被勘探地震工作者承认的Ｈｕｄｓｏｎ本构方程出发，首先讨论了裂隙含流体、气体各向异性介质一维二、三分量波动方程，然后用有限元方法进行了波场数值模拟。该研究有益于加深对地震波在裂隙含流体、气体各向异性介质中传播规律的认识。     

裂隙含流体、气体各向异性介质波场数值模拟

由目前普遍被勘探地震工作者承认的Ｈｕｄｓｏｎ本构方程出发，首先讨论了裂隙含流体、气体各向异性介质一维二、三分量波动方程，然后用有限元方法进行了波场数值模拟。该研究有益于加深对地震波在裂隙含流体、气体各向异性介质中传播规律的认识。     

各向异性上覆层对裂缝储层地震AVAz影响的模拟及分析

为认识和研究上覆层为VTI时,垂直裂缝发育储层的地震反射AVAz响应特性,提出了相适应的正演模拟方法。该方法类比Shoenberg传播矩阵法中水平慢度分量的表达形式,引入了用入射角、方位角和VTI介质中相速度来表达的水平慢度分量作为输入,简化了各向异性反射率法的框架。利用该方法模拟和分析了各向异性层状模型的地震反射响应随方位的变化,讨论了入射角、上覆层各向异性程度对目标层顶底界面反射波AVAz特征的影响。结果表明,采用适当大的观测入射角进行AVAz分析更为合理,上覆层为VTI时会减弱目标层的AVAz特征。测井尺度理论模型的数值模拟说明了该方法的实用性。     

复杂介质的三维块状模型快速射线追踪

三维复杂介质的射线追踪,是目前在生产实际中应用最广泛的正演模拟技术。这里采用了块状结构的建模方法构造复杂三维地质体,以试射法作为射线追踪的基本算法。在试射过程中,引入了聚类的基本思想进行聚类加密,提高了射线追踪的效率。为了得到光滑连续的地层曲面,采用基于方位角的平滑处理对曲面进行平滑,可以进一步提高三维射线追踪的计算速度。同时,通过对迭代算法的改进,能够保证复杂模型的快速射线追踪。实际地质模型的射线追踪结果证明,改进后的射线追踪方法能够有效地提高复杂介质三维模型的射线追踪速度。     

济阳坳陷罗家地区各向异性页岩储层全波场地震响应模拟及分析

济阳坳陷罗家地区页岩储层的固有各向异性及水平层理发育使得储层呈现VTI（vertical transversely isotropy）各向异性,而VTI各向异性背景下垂直裂缝的发育使得储层进一步呈现等效正交各向异性特征。本文以正交各向异性介质作为罗家页岩油储层模型,在岩石物理建模中应用Backus平均理论将测井尺度的VTI各向异性粗化至地震尺度,并利用Schoenberg理论在VTI各向异性背景中引入垂直裂缝,进而得到正交各向异性等效介质模型,同时考虑裂缝尺度和流体填充等因素。之后,应用各向异性反射率法进行全波场地震正演模拟,计算正交各向异性页岩油储层AVAZ（amplitude versus azimuth）响应,通过振幅的方位特征进行储层裂缝识别。计算结果表明,PP波、PSV波和PSH波方位振幅响应各不相同,且方位振幅分布图的拟合形状可反映裂缝的发育方向,为页岩油储层中裂缝的地震识别提供依据。     

正交各向异性介质中地震记录的有限差分数值模拟

本文首先介绍了正交各向异性介质弹性波动方程的建立过程,然后讨论了高阶有限差分算法的实现原理、算法的稳定性和收敛性以及吸收边界条件的选取,并对正交各向异性介质作了正演模拟。通过波场快照、共炮点地震记录和ＶＳＰ模拟记录等正演分析,结果表明该算法在模拟各向异性介质中弹性波的传播是有效的     

基于Matcom 平台的VC + + 与Matlab 混合编程地震正演模拟

正演模拟是地震勘探的重要方法之一。为解决基于Matcom 4. 5 的VC + + 与Matlab 接口程序的地震正演模拟问题,通过用Matcom 4. 5 与VC + + 结合编程的新方法实现EDA 介质中的地震波正演模拟的实例,结果显示Matcom 实现VC + + 与Matlab 混合编程并提高了编程效率。     

TTI介质拟声波方程数值模拟

TTI介质qP波数值模拟方法因为考虑了倾角因素,可以比VTI介质qP波数值模拟方法更加准确地描述各向异性介质中地震波场的传播规律。文中用拟声波方程对TTI介质中的地震波场进行了高阶有限差分数值模拟,在改进衰减函数分布方式后,通过坐标变换,利用改进的完全匹配层(perfectly matched layer,PML)边界控制方程对波场边界进行吸收处理,取得了良好的效果;然后分析了拟声波方程数值模拟中的稳定性问题,并对波场中的伪横波进行压制。通过对不同模型的数值模拟,验证了文中使用的TTI介质拟声波波动方程的稳定性以及所采用的PML边界控制方程的可靠性和适用性。     

改进射线追踪算法的微震源反演

基于射线追踪理论,提出了一种改进的试射法路径追踪算法,可以较准确地模拟射线在层状介质中的传播,不但可计算出任意检波点能接收到的所有的反射波的旅行时及射线路径,而且具有快速、准确、适应复杂构造模型的优点。利用共扼梯度法,通过正演模型的修改迭代,研究微地震层状模型的反演算法。最后,对实际的检测资料进行反演,反演结果基本上和实际资料吻合。     

多尺度随机孔洞介质地震波正演及其逆时成像

常规地震勘探理论无法完全解释复杂多变的地质特征,较难达到高精度地震资料处理的目的,因此必须发展适应复杂非均匀介质的地震波场传播理论,以更好地刻画真实的复杂构造、岩性、裂隙油气藏的特点。根据实际火山岩样的随机特征,构建了相应的多尺度组合型随机孔洞介质模型,并进行了高精度地震波正演数值模拟和逆时成像研究,制作了在不同频率条件下的自激自收剖面和逆时偏移剖面,计算结果表明,构建的多尺度随机孔洞介质模型可以有效模拟真实的地质变化特征,同时指出准确区分由地层随机非均匀性引起的散射波场响应和干扰信号是实现高分辨率地震勘探的前提。     

TTI介质弹性波频率-空间域有限差分数值模拟

由周期性薄互层引起的VTI介质是研究比较广泛的一类各向异性介质。当VTI介质对称轴偏离垂向,本构坐标系与观测坐标系不重合时,会形成观测坐标系下的TTI介质。引入25点优化差分算子,推导出二维TTI介质频率域弹性波动方程;为压制边界反射,采用完全匹配层法吸收边界条件,并计算出优化差分系数;最后采用集中力源,模拟了弹性波在TTI介质中的传播过程。从波场快照和地面共炮记录可以看出,笔者采用的数值模拟算法能有效压制数值频散。TTI介质中的波场传播比较复杂,纵波传播相对稳定,横波波前的三分叉现象比较明显,并存在振幅奇异性。当VTI介质的对称轴偏转后,还会增加地面地震记录的复杂性。     

含直立裂隙介质的弹性波动方程正演模拟

从具有水平对称轴的横向各向同性（HTI）介质中的弹性波动方程出发,在交错网格空间中采用高阶差分算子对弹性波动方程进行差分离散,得到了HTI介质中地震波正演的高阶有限差分格式,研究并实现了PML吸收边界条件。在此基础上实现了HTI介质中弹性波方程的多波正演。数值算例表明,该方法能够精确模拟弹性波在复杂各向异性介质中的传播过程,得到高精度的正演记录。     

三维方位各向异性介质数值模拟及波场分析

采用交错网格技术将速度-应力方程中的速度对时间的导数转化为应力对空间的导数, 将弹性波动方程表示为与二阶双曲方程等价的一阶应力-速度公式, 以实现三维三分量地震波场模拟。对获得的波场快照、VSP记录图像和地面记录图像进行波场分析, 发现在波场图像中存在明显的拟P波、拟快横波和慢横波, 还出现了横波分裂、横波分裂盲区及波面三分叉等地震波在方位各向异性介质传播时产生的特殊现象。研究结果表明, 用交错网格三维模拟方法研究方位各向异性介质对方位各向异性介质中地震波传播反演能起到很强的辅助作用。     

利用高阶交错网格有限差分法数值模拟VTI介质井孔声场

横向各向同性（TI）介质是岩石地球物理中常见的一种现象,研究其井孔声场传播特征对声波测井理论以及为声波测井解释提供依据具有重要意义。针对具有垂直对称轴的横向各向同性（VTI）介质,根据柱坐标系条件下的弹性波波动方程,推导了速度-应力交错有限差分公式,采用时间二阶、空间十阶的交错有限差分算法对VTI介质中的井孔声场进行数值模拟。给出了在均匀介质中井孔声场不同时刻的波场快照,以及不同各向异性系数的VTI介质中的波场快照,计算了井轴上声源激发出的声波全波列波形。结果表明,在其他条件不变的条件下,VTI地层的各向异性系数的增大对横波的传播影响不大,但会使得纵波在纵向上的传播速度相对变小,径向上变化不大。各向异性系数的增大会使声波测井全波列首波信号时差变大,声波幅度略变小。