海底天然气水合物OBS多分量地震正演模拟（英文）

根据Ecker的水合物沉积物的三种微观模式,计算含水合物沉积层和含游离气沉积物的弹性模量,分析对比了水合物的不同微观模式、不同水合物饱和度以及不同游离气饱和度对沉积物弹性模量的影响;从纵横波分离的弹性波动方程出发,采用交错网格空间有限差分方法模拟地震波在海底天然气水合物沉积地层的传播,得到纵、横波的海底地震（OBS）共接收点道集。数值算例表明,当水合物作为流体的一部分或胶结颗粒骨架时,仅纵波记录上存在BSR;当水合物胶结颗粒接触,纵、横波记录上均存在BSR。并且,OBS会接收到上行纵波和上行横波在海底界面形成的转换波,干扰横波记录上BSR的识别。     

低信噪比浅水区OBS测线初至拾取研究

浅水区OBS(Ocean Bottom Seismograph)测线往往具有信噪比较低,空间假频发育严重的问题,导致OBS测线初至拾取困难,进而影响反演速度模型的可靠性.本文通过直达波校正,正演模拟校正和射线互逆实现对浅水区OBS测线初至的准确拾取和质控.直达波校正消除主要的时钟漂移量.通过多道地震与OBS的联合处理,充分利用多道地震浅层速度模型较准确的特点,通过正演模拟构建OBS道集,弥补了实际采集的OBS道集浅层空间假频严重的问题,正演模拟法基本消除了剩余时钟漂移量,校正后的OBS道集符合层析成像的要求,同时也解决了近偏移距初至拾取的问题.利用射线互逆的方法,通过不同OBS道集组合不同偏移距组的对比验证,解决了中远偏移距初至拾取对比验证问题.三步法基本解决低信噪比浅水区OBS测线初至拾取的问题,为速度反演获得准确的模型奠定了基础.     

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

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

三维波动方程有限差分正演方法

三维地震资料的处理和解释都需要有效的三维正演模型予以验证．本文提出一种在ｘ－ｔ域实现快速、高精度的有限差分正演方法,采用了独特的“平行四边形网格”,并用Ｐ－Ｒ交替差分格式使三维波动方程可以局部地分裂成二维求解方程,从而有效地减少运算量．     

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

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

三维波动方程有限差分正演方法

三维地震资料的处理和解释都需要有效的三维正演模型予以验证．本文提出一种在ｘ－ｔ域实现快速、高精度的有限差分正演方法，采用了独特的“平行四边形网格”，并用Ｐ－Ｒ交替差分格式使三维波动方程可以局部地分裂成二维求解方程，从而有效地减少运算量．     

高精度有限差分地震波正演方法

三维地震模拟不仅可以更精确地研究地震波的传播规律,而且是三维地震资料处理和解释的工具。本文采用精细积分法用于求解波动方程.对波动方程在空间采用差分格式,时间域采用积分法求解析的方法。文中详细论述了精细积分法的数值方法,并给出了计算公式。理论分析和数值算例的结果表明了用这种混合方法得到的解与精确解十分吻合,比有限差分法具有更高的精度。文中给出的地震波正演模拟算例说明了该方法适用于复杂地表和复杂构造地质体。     

傅里叶有限差分法三维波动方程正演模拟

傅里叶有限差分(FFD)法兼有相位屏法和隐式有限差分法二者的优势,能够处理复杂地质构造中的波传播问题,但在三维情形下,算子的双向分裂会引起明显的方位各向异性误差.本文用Fourier变换计算双向分裂过程中的高阶交叉项,消除了方位各向异性误差.该方法充分利用了FFD法在双域实现的算法结构,明显减少了由于引入误差校正所带来的计算量.将该方法应用于修改后的三维French模型的地震正演问题,并将得到的叠后记录、单炮记录同全波有限差分法的模拟结果进行对比,结果证实了该方法对一次反射波具有较高的模拟精度,在内存需求和计算效率方面则具有更大的优势.     

用于声波正演模拟的时空域高精度交错网格有限差分方法

地震正演模拟是逆时偏移和全波形反演中的核心问题之一,因为它们都需要高效、高精度地模拟波场正向和反向传播。为了提高数值模拟的精度,人们广泛采用高阶有限差分方法,但是大多数方法仅在空间上具有更高的精度,在时间上只有二阶精度。首先系统介绍时空域高精度交错网格有限差分方法的基本原理,然后利用模型验证方法的有效性,结果表明:时空域高精度交错网格有限差分方法拥有比常规交错网格有限差分方法更低的数值频散。      

基于模型正演的折射静校正方法研究与应用

随着地震勘探区域地质条件复杂化,静校正在资料处理中的作用日显重要.本文模拟研究了折射渡在不同地质条件中的特性,完成了正演模型下的折射波静校正量计算,并将静校正后记录与理论模型记录对比,分析该方法的准确性,总结折射波静校正方法的技术难点,提出一种更准确优越的解决方法.最后,运用该方法选取最合理的参数,并配合先进的处理软件...     

塔里木地区勘探地震正演模拟研究

以塔里木前陆盆地地质背景为模型,采用三维任意差分精细积分方法并行算法实现了库车地区三维正演模拟.三维任意差分精细积分方法通过时间域采用局部积分半解析方法求得波动方程的递推算子,与常规的差分法相比计算精度有较大提高;文中计算稳定性根据实际算例进行分析,采用稳定因子约束,得到较好的稳定性;边界条件采用改进的自适应吸收边界,并通过串行程序并行化,大大减少了三维正演模拟的耗时,完成了大数据量三维正演模拟.     

基于约束Markov决策过程的初至自动识别技术

随着地震数据采集技术的进步,地震数据量日益增加,全自动、高精度的地震初至走时拾取技术受到了更加广泛的关注.本文将初至拾取看作特征空间内带约束的Markov决策过程,在奖励函数空间,按一定准则全局寻优获得积累奖励值最大的路径,从而达到在高维空间自动拾取初至信息的目的.同时,状态值函数中包含与距离相关的折扣因子γ,使Markov决策过程拾取初至能够考虑地震数据的横向连续性,并且回避地震数据中的坏道信息.在此基础上,本文方法进一步引入受空间几何信息约束的动作（Actions）和转移概率（Transitions Probability）,从而降低了对起始状态和折扣因子选取的难度,让地震数据初至走时拾取更加准确和自动化.实际数据测试结果表明,在初至能量较弱（信噪比较低）情况或浅层存在相邻较近复杂波形时,本文提出的约束Markov算法仍能准确地进行初至走时的自动拾取,并且具有一定的质量监控能力,让拾取结果更有物理意义.

     

基于Lorenz规范的空间-波数域大地电磁三维正演

大地电磁勘探方法由于其成本低、施工简单、探测深度广等优点,广泛应用于矿产资源普查、油气勘探和深部构造研究等领域.如何提高大规模三维大地电磁数值模拟的精度和效率一直是研究热点.本文基于空间-波数域方法,实现了基于Lorenz规范的空间-波数域三维大地电磁数值模拟.基于二次场计算原理,引入Lorenz规范,将Maxwell方程组转化为关于二次场矢量位的亥姆霍兹方程;利用水平方向二维傅里叶变换,将空间域三维偏微分方程转换为多个波数下相互独立的常微分方程,方程采用二次插值有限单元法计算,得到定带宽线性方程组,方程计算量小、并行性好,采用追赶法求解,提高了算法效率;引入压缩算子,用迭代法逐次逼近真实解.充分利用了空间-波数域方法数值精度高、内存需求少、效率高的特点.设计棱柱体模型验证了算法的正确性、分析了算法的收敛性,说明算法对不同频率、不同电导率对比度模型均具有很好的适应性.利用Dublin（DTM1）模型进行三维大地电磁数值模拟,结果表明：在满足精度要求的前提下,空间-波数域算法比空间域算法占用内存少、耗时短;相比基于Coulomb规范的空间-波数域算法,基于Lorenz规范的空间-波数域方法耗时更短、占用内存更少,效率提高至少3倍以上,体现了新方法的优势.

     

基于稀疏表达的OBS去噪方法

海底地震仪（OBS）采集数据的去噪处理是开展OBS震相分析及后续处理反演的基础.本文结合曲波（Curvelet）变换及压缩感知提出一种稀疏化表达的OBS去噪方法,并通过与小波变化对比等探讨去噪效果.曲波变换具有抛物尺度及识别线性异常的优点,可以稀疏地表示OBS数据,再结合压缩感知思想对稀疏表达OBS数据进行去噪处理和重构.通过对变换后的系数进行基于L₁范数的冷却阈值迭代滤波,获得最优的变换系数,本文指出基于曲波变换的冷却阈值迭代法能够很好地对OBS数据去噪.对比小波和曲波两种变换在相同迭代次数下对理论模型数据进行去噪处理,表明曲波变换得到的结果信噪比更高.利用本文方法对渤海地区采集的OBS数据进行去噪处理获得了更加清晰连续的震相,噪声压制效果更明显,为震相拾取及后续速度模型反演奠定了良好的基础.

     

基于平均导数优化方法的VTI介质频率空间域正演

本文提出了一种新的基于平均导数优化方法(average-derivative optimal method,简称ADM)的二维VTI介质qP波波动方程频率空间域二阶9点格式,这种新算法将二维VTI介质qP波波动方程中中心空间导数项的差分近似表示为正交方向上3个网格点的加权平均形式.通过最小二乘优化方法求取空间导数项和加速度项的加权优化系数从而使数值频散达到极小化,每个波长所需要的网格点数在1%的误差范围内仅为3.57个网格点数,而VTI介质常规9点差分格式在相同的误差范围内则需要约12个网格点数,新方法的计算精度明显提高.复杂BP2007 2D VTI海洋标准模型数值模拟结果也验证了本文VTI介质9点ADM算法的有效性和准确性.     

薄层时频特征的正演模拟

基于正演模拟的薄层时频特征响应分析是薄层定性识别与厚度定量预测的基础与关键.本文基于波场延拓理论,利用深度域相移法对不同厚度薄层进行了正演模拟.采用广义S变换对零偏移距地震道的反射复合波进行瞬时频谱分析,发现当薄层顶底反射系数极性相反时,复合波的干涉作用具有升频降幅作用,而极性相同时具有升幅降频作用;理论推导和实验分析均证明主峰值频率与薄层双程旅行时间厚度存在定量解析关系,这为薄层厚度定量预测提供了重要的技术手段;由于峰值频率对反射系数大小与极性变化不敏感,因此,峰值频率方法与时域振幅方法相比更具稳定性.     

基于二次场的可控源电磁法三维有限元-无限元数值模拟

本文提出了一种新的混合有限元-无限元三维可控源电磁法（CSEM）问题快速高精度正演模拟算法.首先从电场双旋度方程出发,推导了水平电偶极子源的二次场边值问题,采用无限元代替截断边界条件和有限元离散内部计算区域的新策略,达到减小计算区域的目的,基于并行直接求解技术,实现多源CSEM问题的快速精确求解.其次,通过层状解析模型测试,一方面验证了新算法的正确性,另一方面通过与其他三种已知CSEM问题求解策略进行对比,表明了本文提出的基于二次场有限元-无限元算法具有离散区域小、求解速度快和计算精度高等优点.最后,通过3D模型计算,清晰直观地模拟了场源阴影效应,为野外数据的处理与解释提供指导.

     

等厚薄互层时频特征的正演模拟

本文基于波动理论,利用深度域相移法对不同互层数、不同单层厚度的等厚薄互层进行正演模拟.采用广义S变换分析零偏移距地震道反射复合波的瞬时频谱.研究发现当单层厚度大于3/16波长时,等厚薄互层时域特征为中、高频等幅振荡,瞬时频谱为与单层厚度成反比的单一峰值梳状函数;当单层厚度小于3/16波长时,震源子波的有限带宽对瞬时频谱的高频频率峰有压制作用,导致等厚薄互层中部的时域振幅和频域最大幅度均为低值,薄互层的结构特征不易分辨,说明宽高频的震源子波是提高薄互层分辨能力的关键.此外,理论推导和实验分析均证明:当地层厚度大于1/8波长时,其峰值频率与薄互层单层厚度存在定量解析关系,这为等厚薄互层单层厚度定量预测提供了重要的技术手段.     

Forward modeling of tight sandstone permeability based on mud intrusion depth and its application in the south of the Ordos Basin

Permeability is an important index in reservoir evaluation,oil and gas accumulation control,and production efficiency.At present,permeability can be obtained th...