山前褶皱带模型正演模拟及波场分析

巴布亚褶皱带是巴布亚盆地的主要油气聚集带,已累计发现40亿桶油当量,但其勘探程度依旧很低.该构造带地表大部分区域为灰岩所覆盖,地形起伏剧烈,地下构造复杂,逆掩推覆构造广泛发育,地震资料采集处理效果很差.为此,在构建逼近巴布亚褶皱带实际复杂地表复杂构造模型的基础上,进行基于波动方程的波场地震正演模拟和地震照明模拟,分析褶皱带复杂地震资料采集和处理的影响因素.结果表明,波场正演模拟方法能够很好地适应复杂地表和复杂构造模型,模拟结果与实际地震资料具有很好的一致性,地面不同位置激发得到的单炮记录和地震照明分布具有很大的差异,探讨了影响复杂构造成像的主要因素.针对工区实际地下目标,提出了地震采集观测系统优化设计和后续处理建议,为复杂地表复杂构造区地震资料采集处理提供了借鉴和参考依据.     

塌陷采空区的正演模拟及波场分析

老窑塌陷采空区造成的透水事故严重威胁着大型矿井的生产与安全,探测出塌陷采空区已成为当前三维地震的一项重要研究课题.针对老窑采空区不同的塌陷类型,建立多种不同内部结构的塌陷采空区地质模型。并采用有限差分法数值模拟技术对其进行正演模拟,分析出不同塌陷采空区在生成时间剖面上的波场特征.结合一些典型的工程实例,总结出塌陷采空区在地震资料上的判别标志以及识别特征,为矿井提供出准确的塌陷采空区位置.     

非均匀各向异性弹性波场正演

推导了用于非均匀各向异性弹性波场正演的伪谱法基本公式，对特征值法边界修正方程进行了全面的理论推导，给出了二维和三维问题固体，流体的特征变量及各类边界修正方程，选用2个模型，模拟了2．5维弹性波场，对S波分裂的偏振图像和时间延迟作了较详细的分析。     

基于射线追踪的微地震多波场正演模拟

对于微地震正演模拟,本文以射线追踪的原理为基础,对两点间的射线追踪问题进行了研究,应用二分算法、改进二分算法和微变网格算法对水平层状匀速模型、弯曲层状匀速模型和复杂地质模型进行射线追踪,使得计算效率和适用范围都得到了很大的改善.文中对每种算法误差范围和计算效率进行了对比验证,对于不同的地质模型,选用合适的算法才能在计算速度和精度上得到双重保证,最后正演模拟了多波三分量记录.在模型建立上引入了超薄层概念,并在前人模拟的直达波、透射波、反射波基础上拟了折射波,使正演模拟的多波场信息更丰富.文中的应用实例及模型结果表明:与二分法相比,改进二分法能够对弯曲界面进行射线追踪,并能保证结果的精度.弯曲层状模型中,改进二分法与微变网格法相比计算速度有显著提高,能够应用到资料的反演中.     

基于一步法波场延拓的正演模拟和逆时偏移成像

通常工业界实现逆时偏移算法时采用有限差分数值方法模拟地震波场,波场模拟常常受稳定性条件限制,且易产生数值频散,成像精度降低.本文引入了一步法波场延拓方法,首先构建声波传播算子,借助Chebyshev多项式和Jacobi-Anger展开式近似传播算子中的e指数项,进而实现波场递推,该方法时间步长的选取不受稳定性条件限制而且不存在空间频散现象.本文将一步法波场延拓方法用于逆时偏移成像的波场模拟,并提出双缓冲区存储策略,在不增加计算量的前提下,大幅降低了逆时偏移方法的波场存储量.波场模拟和逆时偏移成像测试表明,本文提出的一步法波场延拓方法模拟地震波场精度高,消除了频散影响,可在较大时间步长的情况下实现高精度波场模拟;提出的基于一步法波场延拓的逆时偏移方法成像质量好;基于双缓冲区存储策略的逆时偏移成像方法存储成本低.

     

离散介质波场正演模拟与弱化性影响分析

从流变力学角度建立的离散介质本构关系同时考虑了介质的弹性、粘性与弱化性,与损伤力学建立的本构关系比较表明弱化因子具有损伤因子的力学意义。本文采用交错网格伪谱法正演模拟了离散介质中波的传播,并以平面波解分析了体波的频散和介质的吸收系数。高频和低频情况下波场数值模拟的结果是:饱和多孔介质中体波速度随弱化因子的增加而降低,介质吸收系数随弱化因子的增加而增加,弱化性对不同体波频散的影响程度不同。     

深水陡坡带绕射多次波压制方法研究

在深水陡坡带和复杂海底构造区域,即使采用长电缆,由绕射产生的自由表面相关多次波因为传播方向的异变,仍然难以避免在自由表面发生下行反射的反射点位置位于电缆长度范围之外,造成预测多次波的近偏移距信息不足或者负偏移距信息缺失等问题,以致不能准确地预测和压制陡坡带的绕射多次波.本文在常规反馈迭代自由表面多次波压制方法的基础上做出相应改进,将邻炮的炮记录依照炮点、检波点互换原理引入到反馈迭代循环当中,补充绕射多次波预测所需要的近偏移距和负偏移距信息.经过Sigsbee2b模型和实际资料测试,达到较好压制深水陡坡带的自由表面相关绕射多次波的目地.     

地震多级散射波正演模拟方法

地震波在穿越地下散射体群时会产生多级散射波,分析其地震响应特征,可推断散射体的分布情况和性质。本文从二维标量波动方程出发,结合地震散射理论和波恩近似理论,推导了多级散射波方程。在此基础上,采用高阶有限差分法对双点散射体模型和复杂散射体模型进行数值模拟,分析了多级散射波的传播规律和波场特征,并通过抽取多级散射记录和各级散射记录的单道记录与参考单道记录的对比,验证了本文推导散射波方程的准确性。      

非规则极坐标系弹性波场正演模拟（英文）

在研究隧道、井孔等特殊地质环境下的波场传播时,极坐标正演模拟比传统的笛卡尔坐标正演模拟精度更高。而且在实际情况下,极坐标环境大多是不规则的,为了解决这一问题,提出了一种非规则极坐标系正演模拟方法提高模拟精度。该方法首先将具有不规则表面的空间在非规则极坐标系下进行网格剖分,并通过变换将不规则表面映射为规则表面,然后在非规则极坐标系下计算波场。在求解非规则极坐标系下的波动方程时,利用全交错网格对其进行差分离散求解。此外,利用人工吸收边界、极坐标系自由边界和圆周边界条件对虚假边界反射进行吸收。选取了三个模型来验证该方法,结果显示,相比于常规极坐标系正演模拟,本论文非规则极坐标系弹性波场正演模拟方法在对具有不规则表面的极坐标空间进行正演模拟时,可以更获得准确、稳定的波场模拟结果。非规则极坐标系弹性波全波形反演结果进一步证明了本文提出的弹性波正演模拟方法可以准确地模拟不规则极空间中的波场。     

远震波场正演模拟方法及应用

地震波场数值模拟是壳幔结构成像和深部探测的重要基础.?经典的远震波场数值模拟主要基于一维地球模型,包括解析法、半解析法和数值法等.?这些算法能够高效地计算理论地震图,但由于将地球假设为一维层状介质,难以考虑介质的横向非均匀性对地震波场的影响.?近年来随着计算机性能的不断提高,三维地震波场数值模拟方法得到快速发展,并被广...     

地震波场射线类正演模拟方法对比

射线类正演方法以其高效,灵活的特点,在复杂地区的勘探中发挥着非常重要的作用.常规射线正演方法只考虑地震波的运动学特征,精度较低,且对复杂构造适应性较差;高斯射线束正演方法作为常规射线正演方法的一种改进,同时考虑地震波的运动学和动力学特征,且无需两点射线追踪,兼具计算效率和精度;菲涅尔束正演模拟进一步从波动理论角度对高斯射线束的有效半宽度进行了限制,改进了波场模拟的精度和稳定性.对三种射线类正演模拟方法进行理论解剖和分析,在高斯射线束正演方法的基础上,用第一菲涅尔带半径约束高斯束束宽实现了菲涅尔束正演算法.基于理论模型和实际模型,对三种射线类正演方法和波动方程有限差分法进行对比分析,结果表明高斯射线束正演方法具有较高的计算精度;菲涅尔束正演方法更加稳定,精度更高,更接近于波动方程有限差分方法.高斯射线束和菲涅尔束均具有较高的计算效率,适合于对波场模拟精度不高但对运算效率有较高要求的应用环境.     

标量波Laplace-Fourier域有限差分正演模拟方法

<正>演是反演的基础,有效正确的正演差分格式可以保证反演结果的精度和效率.本文通过声波波动方程域间转换提出标量地震波Laplace-Fourier域数值模拟方法,并推导了同时引入衰减因子和频率的Laplace-Fourier域标量波方程的9点法有限差分格式和Laplace-Fourier域对应的加入PML(perfectly matched layer)吸收边界条件的差分格式,并通过模型试算验证了本文提出的Laplace-Fourier域正演方法的有效性和准确性,通过与时间域正演方法得到的地震记录比较,可以看出该方法能满足正演数值模拟的要求,为下一步进行Laplace-Fourier域标量波全波形反演奠定了基础.     

基于波场分离的接收函数正演及偏移

针对理论接收函数的正演和偏移,本文采用基于波场分离的有限差分数值模拟的算法,借鉴成熟的地震勘探处理技术,建立了基于波场分离的一阶速度-应力方程,实现接收函数阵列的正演和偏移.利用叠前逆时偏移的理论对接收函数进行成像.数值计算的结果表明,基于波场分离理论的正演能够得到接收函数的纯P波阵列和纯S波阵列,利用叠前逆时偏移的方法能够进行高精度成像.通过对比纯P波阵列、纯S波阵列和全波场阵列的成像结果发现,利用接收函数的纯S波阵列进行叠前逆时偏移能够获得较高精度的成像结果.     

三维孔洞储层建模及其地震波场正演模拟

本文研究的内容有二,一是针对塔河油田六区东中、下奥陶统碳酸盐岩孔洞储层发育的特征,通过精细的层位标定、构造解释与时深转换,搭建深度域等时地层框架模型,然后采用协克里金算法,将波阻抗反演体作为协变量,对孔洞储集体的三维孔隙度分布进行确定性建模;二是应用二维横向各向同性介质中的弹性波波动方程,基于交错网格差分格式,考虑初始...     

广角反射波场特征研究及正演模拟分析

本文主要讨论了广角反射波场特征及影响广角反射波波形识别的几个主要因素，通过高斯射线正演模拟分析，以及对实际广角反射资料(WAR)进一步处理与研究，得到了一些十分有意又的结论，这些结论可以作为一些特殊地质构造地区进行广角反射地震勘探的理论依据。     

声波介质一次散射波场高斯束Born正演

Born正演是一种常用的地震波场正演模拟方法,也是线性化地震反演的理论基础.在实际应用时,Born正演通常结合常规的地震射线方法进行实现.为了克服常规地震射线方法的弊端,并且保证地震波场的模拟精度和计算效率,本文提出了一种基于高斯束的一阶散射波场Born正演方法.该方法分为两个环节：首先,我们利用高斯束的走时和振幅信息将地下散射点处的反射率映射为地表束中心位置处的局部平面波;然后,我们利用逆倾斜叠加将局部平面波转化为接收点处的时空域散射波场.在具体的实施过程中,我们提出一种以wavelet-bank方式实现的局部平面波合成方法,同现有的算法相比,可以在保持计算精度的同时,大大减少计算时间;此外,我们还利用最速下降法优化了高斯束的迭代循环过程,进一步提高了Born正演的计算效率.两个模型的应用效果证明,本文所提出的高斯束Born正演方法可以精确、高效的实现声波介质一次散射波场的正演模拟,为三维大规模地震波场的正演问题提供了一种切实可行的实现方案.

     

基于AVO正演模拟的自由表面多次波压制方法

针对SRME和Radon变换压制自由表面多次波的近道多次波残留和伤害有效波等问题,主要采用基于AVO属性分析的多次波剔除法对近道多次波进行压制,该方法在迭代更新过程中,受近道多次波影响容易出现更新误差,同时迭代更新降低计算效率.本文对基于AVO属性分析的多次波剔除法进行改进,利用中远道地震数据进行约束,通过AVO正演模...     

基于单程波正演模拟的多尺度地震资料特征分析

多尺度地震资料在油藏地球物理中发挥着重要的作用,通过开展单程波正演模拟并对模拟结果进行偏移,有利于深入剖析多尺度地震资料之间的联系与差异.正演模拟表明,地面地震资料、VSP资料和井间地震资料三种尺度的地震资料都是相同地下地质体对应于该观测尺度下的地震波响应特征,差异在于由于观测方式的不同,资料所包含的波场信息丰富程度不一样,反映地质体信息的范围不一样,在分辨能力上具有一定的差异性.地面地震资料具有较强的横向连续性但纵向分辨率低,VSP资料具备比较精确的时深关系,井间地震资料分辨率高但是覆盖范围小,对三种多尺度地震资料开展综合研究有利于更好的揭示地下地质信息.     

黏弹性VTI介质频率空间域准P波正演模拟

有限差分方法是波场数值模拟的一个重要方法,时间域有限差分计算方法因按时间片递推计算,每个时间片的舍入误差会累积到下一片中,当时间片较多,最终会导致累积误差太大.而频率域计算是按频率片对空间网格进行整体求解方程组,其计算误差分配到了每个网格点上,并且各个频率片之间是独立计算的,因此不存在累计误差,而且在频率-空间域更易于模拟黏弹各向异性介质的衰减,同时采用25点优化差分算子克服了常规差分算子的数值频散.本文首先推导了黏弹性VTI介质的频率-空间域准P波波动方程,在此基础上,利用25点优化差分算子构造了差分格式,实现了频率空间域的正演模拟.并以断层模型为例,模拟并对比分析了利用弹性VTI介质波动方程差分格式和利用黏弹性VTI介质波动方程差分格式模拟的井间地震准P波波场记录,得出黏弹性介质地震记录由于衰减振幅相对较小,主频相对较低的结论.     

弱各向异性介质弹性波的准各向同性近似正演模拟

准各向同性(QI)近似可用于弱各向异性介质的正演模拟.本文通过运用QI方法的零阶和一阶近似,计算了VTI介质模型的地震记录.得出的地震记录与标准各向同性射线理论(IRT)和基于伪谱法的三维地震正演模拟得出的地震记录作了比较,可以认为是精确的合成地震记录.