简介:传统上,有限差分的差分系数一般可以通过泰勒级数展开法或优化方法来极小化频散误差得到。基于泰勒级数展开的差分法在有限的波数范围内精度较高,但在这个范围之外会产生较强的数值频散;基于最小二乘的优化有限差分法能在更大的波数范围内达到较高的精度,并可以在较小的计算需求内获得全局最优解。本文将基于最小二乘的优化有限差分法从二维正演模拟推广到三维,形成了计算效率高、高精度范围宽、适合并行计算的三维声波优化有限差分方法。频散分析及正演模拟表明本文发展的有限差分方法可以很好地压制数值频散。最后,将本文发展的有限差分方法应用到三维逆时偏移的震源波场延拓和检波点波场延拓中,并结合有效边界存储策略与checkpointing技术在GPU集群上实现三维逆时偏移以提高计算效率、减少存储量。三维逆时偏移试算结果表明本文三维优化有限差分方法与传统的有限差分法相比可以获得更高精度的偏移成像结果。
简介:频率-波数域单程波算子能高效地模拟地震波在复杂介质中的传播,但是在描述波的大角度传播和速度横向扰动变化较大介质中传播的问题时仍然存在一定误差。这类误差是由于对单平方根算子使用Taylor展开式的近似程度不足所造成。为了进一步提高泰勒展开式的精确性,本文提出一种利用粒子群智能算法优化级数展开系数的高阶广义屏算子对单平方根算子的展开级数进行优化处理。新的偏移算法能在保持单程波偏移算法高效的前提下进一步提高偏移算子在大角度的成像精度和对强横向速度变化介质的适应性。通过脉冲响应实验,验证了基于粒子群算法优化级数的高阶广义屏算子能够提高常规的高阶广义屏算子的成像精度和成像角度。根据对二维SEG/EAGE盐丘模型的成像处理,基于粒子群算法优化级数的高阶广义屏算子对盐丘下面的断层取得了更高质量的成像,说明粒子群优化级数的高阶广义屏算子比常规的高阶广义屏算子具有更好的横向速度适应性。为了检验本文所提算法对实际资料的处理能力,我们利用常规的偏移处理技术和本文所提算法对一条海上二维数据进行了偏移成像处理,对比分析成像剖面发现本文所提算法描述了更加清晰的层位信息和更高质量的偏移剖面。本文所提算法能有效提高高阶广义屏偏移在广角度成像的能力,具有一定实际应用价值。
简介:常规长排列非双曲动校正公式是在VTI介质中得到的,它不能满足任意空间取向TI(ATI)条件下的扩展.本文以VTI介质中非双曲动校正公式为基础,基于我们推导得出的ATI介质中精确四次时差系数解析解和NMO速度解析解,给出ATI介质中长排列优化的非双曲动校正公式.通过与各向异性射线追踪方法计算所得出的"精确走时"结果对比,研究表明优化后的非双曲动校正公式能精确地描述任意强弱、ATI介质中随测线方位变化的走时曲线,可以用来替代耗时、多偏移距、多方位的射线追踪方法正演拟合ATI介质中长偏移距反射走时,为利用非双曲时距的各向异性参数反演提供理论基础性认识。
简介:频率空间域地震波数值模拟具有独特的优势:可以同时模拟多源的波传播、每个频率之间独立并行地计算、计算频带选择灵活、不存在累计误差、容易模拟粘弹性介质中地震波传播。但是该方法的最大瓶颈是对于计算机内存的巨大需求。我们使用压缩存储系数矩阵的方法,极大地减少了计算机内存的需求量。同时为了减少短差分算子的数值频散,引用了频率空间域25点弹性波波动方程的差分格式,并使用了最小二乘意义下求出的优化差分系数。为了克服边界反射,采用了最佳匹配层吸收边界条件。数值模拟试验证明:用压缩存储系数矩阵及优化差分系数的频率空间域25点差分格式进行弹性波正演模拟,可以减少数值频散,提高计算精度。使用较大的网格间距,降低计算机内存需求,并保持较高的计算效率。该正演方法为后续弹性波偏移和弹性参数反演提供较好的基础。
简介:为了探明介质的非完全弹性对瑞利型槽波传播的影响以及利用瑞利型槽波品质因子QR反演煤层碳化程度、裂隙、煤层厚度等煤层属性,在本文中,我们将复速度引入瑞利型槽波频散方程,计算了煤层瑞利型槽波的品质因子QR曲线;采用控制变量法,分析了瑞利型槽波品质因子QR随煤层厚度、煤层纵横波品质因子以及围岩纵横波品质因子的变化规律。研究表明:瑞利型槽波品质因子QR曲线与群速度曲线变化趋势一致;埃里相(Airy-phase)频率附近瑞利型槽波的品质因子QR最小,且埃里相频率随煤层厚度增加而减小;瑞利型槽波品质因子QR随煤层横波品质因子QS2的增大而增大。利用瑞利型槽波品质因子QR曲线可以实现对瑞利型槽波的吸收衰减进行补偿、对煤层厚度进行预测以及岩性成像。更多还原
简介:地震低频信息能够提高分辨率与成像精度,改善反演质量,甚至直接进行油气检测,需要对其进行有效保护与拓展。对于子波而言,缺失低频信息会导致主瓣幅度降低、第一旁瓣幅度增加,并出现次级旁瓣呈周期震荡衰减的现象;从合成地震记录和典型地质模型来看,低频缺失会产生假同相轴,造成分辨率提高的假象,且模型不同位置的特征存在一定差异;对缺失低频的模型数据进行波阻抗反演,会造成构造失真、岩性改变的假象,特别是高陡构造和薄互层。针对缺失低频的地震资料,本文还研究了基于压缩感知与稀疏约束的拓频方法,开发了相应的模块,并对实际CIP道集进行处理,取得了较好的应用效果。
简介:本文以中观孔隙结构的White模型为基础,构建了部分饱和孔隙介质模型,利用Biot方程的建立思路和Johnson推导的体变模量,推导了部分饱和孔隙介质中的纵波方程,并以平面波为例,求取了方程的衰减系数,分析了地震勘探频带范围内地震波的衰减特性。结果表明:在部分饱和孔隙介质中,地震波在低频段也会发生明显的衰减和频散现象,频率越大,衰减越大;且第二纵波的衰减比第一纵波更为明显;这一结论弥补了Biot理论在描述地震勘探频带范围内波的衰减现象的不足。文中还研究了孔隙度、饱和度和模型内径尺寸对第纵波衰减特性的影响机理,主要表现在在地震勘探频带范围内,波的衰减随孔隙度的增大而增大,随含油气饱和度的增大而减小,当孔隙内径尺寸小于二分之一外径尺寸时,波的衰减随内径尺寸的增大而增大,当内径尺寸大于二分之一外径尺寸时,波的衰减随内径尺寸增大而减小。
简介:大型、构造复杂的三维物理模型可用于模拟油气勘探。构造逼近实际地质状况的模拟具有制作技术难度大、质量控制严格等特点,可用于采集宽方位、多方位和全方位的地震数据,从而进行多种三维处理、解释方法验证。本文针对中国西部前陆盆地地表条件复杂地下构造复杂,导致成像不理想等问题,基于复杂的地下构造,设计制作了目前世界上模拟施工面积最大、构造最复杂的KS(塔里木盆地克深勘探工区)物理模型。本文的模型技术的进步主要涉及3个方面:模型的设计方法、模型的浇铸流程和数据采集,首次给出了物理模型的三维真实速度模型,定量分析了物理模型的制作精度,绝对误差小于3mm,可以满足方法试验的需要。该模型基于三维形态测量技术建立了三维真实速度模型,可作为方法试验的基础数据。因此,该模型可作为地震物理模拟技术的标准。
简介:储层的裂缝分析在油气勘探中已变得愈来愈重要。品质因子Q值是表征地下介质对地震波吸收衰减特性的一个重要参数,不仅能够反映介质内部的本质特征,同时也可用于裂缝识别。在储层地震属性和反演中,地震叠前资料比叠后资料包含更加丰富的储层信息,因此,利用叠前地震道集进行吸收衰减参数提取能够进一步提高衰减系数估算的准确性。本文提出了一种在改进S变换基础上利用叠前地震CMP道集资料提取品质因子Q的方法。改进的叠前Q提取方法,首先是利用时频分辨率可调的改进S变换对叠前CMP道集进行时频分析,推导了利用改进s变换的频谱比法计算公式,其次采用能谱密度比替代传统的振幅比来逐道求取估算Q值的频谱比的斜率;然后建立谱比斜率与炮检距之间的关系式,采用多道拟合的方法消除炮检距影响,从叠前地震道集资料中提取准确的品质因子Q。最后本文把该叠前Q值计算方法引入到前进潜山储层裂缝预测中,纵波吸收衰减的各向异性具有好的实际应用效果。
简介:常规声波测井解释方法建立在各向同性地层中的垂直井基础之上,对各向异性地层中的水平井和斜度井并不适用。在水平井和斜度井的钻井过程中,岩屑会溅到井壁并最终在井壁底部形成一层薄岩屑层。同时,高速层和各向异性也可能会影响声波时差的测量。本文采用结合了Hybrid-PML的交错网格时域有限差分方法(SGS-FDTD),研究阵列声波仪器在各向异性介质中的水平井和斜度井条件下,不同因素对测得时差所带来的影响。模拟结果揭示了声波时差在高速层,薄岩屑层,井斜角,层厚,各向异性等不同影响因素下的相应规律。我们发现当波长远小于井壁与高速层的距离时,能够准确测得目的层时差;然而在两者处于同一量级时的近场条件下,几何声学理论不再适用;当井眼底部存在岩屑层时,费马原理依旧适用,对真实时差测量没有影响;在各向异性地层环境下,井斜角的增大使得测得时差由垂直向时差逐渐趋近于水平向时差;对于特定源距的声系,仪器在目的层中所移动的距离超过一定距离时,方可获得薄层时差值。基于不同井斜角和不同各向异性大小均匀TI地层模型的模拟结果,绘制了时差估计图版,可以对任意井斜角和各向异性比的模型进行时差的定量估计。通过对不同声系和不同弹性参数的模拟实验,证明了这种时差估计方法对TI地层中的水平井和斜度井准确有效。
简介:南海的构造与演化与资源环境等关系密切是本文研究重点。本文针对南海的东北部构造及其块体构造方向,利用所采集的区域地震剖面,通过解析地震相与构造及其演化的关系,提出以下观点:(1)构造分区特点明晰,可划分为五个不同构造单元,构造单元之间既有联系又相互独立;(2)南海沉积盆地无论表现为拉张-弱挤压-强挤压的何种构造格局,其区域构造应力场是统一的;(3)首次发现反射地震剖面上显示出两个浅俯冲点。每个块体构造层呈手风琴风箱式折曲并向东聚敛,体现沉积盆地从发育、成长、结束、消亡不同阶段在南海的表现,其块体俯冲方向以及块体包络区域性倾伏方向均与区域应力场方向一致。更多还原
简介:希尔伯特—黄变换(Hilbert-Huangtransform,HHT)是一种新的适合非平稳和非线性信号的分析方法,由于地震信号一般呈现出非平稳与非线性特性,因此HHT非常适合地震信号的分析。本文首先介绍了HHT中关于经验模态分解(EmpiricalModeDecomposition,EMD)的实现过程,在此基础上分析了几种基于EMD获得本征模态函数(IntrinsicModeFunctions,IMF)来计算瞬时频率的算法,其中利用了两个采样间隔瞬时频率的平均来计算瞬时频率,较好地反映了地震信号频率成分随时间变化的特征。将该方法应用于四川东北部某地区海相碳酸盐岩地层三维地震叠后偏移数据处理,提取"三瞬"地震属性,与传统的希尔伯特变换提取的"三瞬"地震属性进行对比,结果表明基于HHT的"三瞬"地震属性结果具有更高的分辨率,IMF2的瞬时相位能够较好地刻画台地边缘生物礁相,IMF2的瞬时频率亦具有较好的分带性。将IMF2的"三瞬"地震属性与钻井等资料结合分析,能够更好地识别沉积相的分布。
简介:在多个区块的测井评价工作发现,孔隙结构直接影响储集层的品质和油气层的电阻率,是测井准确评价流体性质的关键。岩石物理资料表明不同区块内影响储层孔隙结构的因素不同,但效果是一致的,即孔隙结构的复杂程度控制着储层的储集能力和渗透能力。孔隙结构的复杂程度影响了储层中导电流体的分布和含量,从而控制了储层的电阻率。储层出现低阻油气层的内因均为复杂的孔隙结构(骨架导电及工程原因除外)。测井储层评价在分析控制储层孔隙结构复杂程度的地质因素及储层分类的基础上,针对不同类型储层采用不同的模型、参数和标准,可以有效的认识储层品质和识别不同类型储层的流体性质。
简介:基于岩石物理模型和广义各向异性Zoeppritz方法在频率域计算裂缝型反射体模型反射波方位地震AVO响应。反射体模型为粘弹各向异性、有限厚度的地质体,其地震反射波形序列包含如下动力学信息,即分界面处介质的波阻抗和非弹性差异、反射体内部波的各向异性传播、在传播路径上的频散与衰减,以及来自顶底界面的反射波的调谐与干涉等。计算表明,速度频散和衰减增顶界面反射波大入射角反射时的振幅,而减弱底界面反射振幅。对于固定人射角的方位地震方位地震响应,PP波反射特征表现为随方位角的增加反射波形序列延续时间变长,而PSV和PSH转换类型反射波的方位各向异性变化特征稳定且受储层厚度影响较小,表现为PSV波反射振幅随方位角增加而增加,PSH波在0°和90°方位无反射能量,在45°方位反射振幅最强。
简介:研究了基于Poynting矢量的角度域逆时偏移成像及成像幅度的校正方法。根据Poynting矢量进行波场角度分解,由此构建局部成像矩阵及局部照明矩阵。在局部成像矩阵中建立的角度域成像条件,有效地消除了低波数干扰,同时可在局部成像矩阵中进行角度域共成像点道集抽取、倾角估计等运算。利用局部照明矩阵进行了基于全波波动方程的时域照明分析,在局部照明矩阵中计算倾角域幅度校正因子。根据逆时偏移的像计算共倾角像,利用校正因子对各角度的像进行校正,进而实现对成像结果的校正,从而实现了一种高效的倾角域幅度校正方法。最后通过SEG/EAGE模型进行数值计算验证了文中所述的计算方法。
简介:岩石物理建模及脆性指数构建是影响脆性预测精度的两大重要环节。现有页岩模型对有机质的模拟争议较大,需要寻找合理的岩石物理理论来模拟富有机质页岩。同时,现有脆性公式种类繁多,各公式的适用性值得探究。本文利用Self-ConsistentApproximationandtheDifferentialEffectiveMedium(SCA+DEM)理论,通过模拟有机质与粘土的耦合性,构建各向异性富有机质岩石物理模型。与前人理论对比,初步验证了本模型的有效性;同时,基于模型构建脆性模板,分析物性参数对各脆性指数公式的影响。结果显示:各脆性公式对不同物性条件下地层的敏感性不同,基于杨氏模量构建的脆性指数对矿物含量的变化较敏感,而基于拉梅系数构建的脆性公式对孔隙度/孔隙流体敏感。应综合各脆性指数公式并结合地层物性信息,以达到最优的预测结果。