测井与全方位道集联合各向异性参数建模及成像

对于TTI介质各向异性成像,精确的各向异性参数估算与建模至关重要。在当前计算机硬件迅速发展及宽方位地震数据采集日益普及的情况下,成像必须考虑介质的各向异性。提出了基于测井数据与全方位角道集联合的各向异性参数建模方法,该方法先采用Walk-away VSP资料及测井数据可获得井位置点的各向异性参数δ,借助全方位共反射角道集,可获得与方位有关的各向异性时差信息。再采用全方位层析对各向异性速度、各向异性参数ε和δ进行迭代优化,建立精确的各向异性参数模型,可有效减少各向异性参数的多解性。通过实际地震数据对此方案进行了验证,该方法在裂缝探测及地震成像方面有较好的适用性。     

TTI介质各向异性参数优化提取方法

对传统基于井震联合网格层析的TTI介质各向异性参数建模结果的分析发现,在地层厚度变化剧烈的情况下,反演得到的参数模型会出现畸变。为了避免这种不合理的现象,提出了一种优化的TTI各向异性参数反演方法,即通过基于各向异性参数δ的井震误差插值方法将地层厚度信息引入各向异性参数反演流程,使反演结果更加精确。实际资料处理结果表明,优化方法可以有效避免地层厚度突变带来的各向异性参数畸变问题,模型参数不受地层厚度变化的影响,因而比传统方法更为精确、合理。     

基于角道集的井约束层析速度反演

本文研究的走时层析方法,是在角度域共成像点道集(ADCIGs,简称角道集)的基础上获取旅行时差,并通过井约束控制层析反演的精度。此法实现的关键步骤为:①利用叠前深度偏移剖面上拾取的层位界面生成初始层析速度模型;②按照选定的角度范围抽取共成像点道集,并求取各层位的旅行时差;③求取共成像点对应的灵敏度矩阵;④建立反演方程组,并加入正则化和井约束反演慢度更新量;⑤依据共成像点道集中的同相轴拉平程度判断所求速度的精度是否达到要求。理论模型试算及实际地震资料处理结果均表明:该方法易于实现,不必考虑复杂的反射问题,每一步迭代都是在前一步深度偏移的基础上进行速度场的更新,能够得到精度较高的速度场。     

高精度深度域层析速度反演方法

深度域层速度场是叠前深度偏移最重要的参数,直接决定了偏移成像的精度。传统层析速度建模方法由三部分组成,首先利用Dix公式或基于层位约束的Dix公式生成初始层速度模型,然后利用炮检距域道集对初始速度模型进行层析迭代,最后利用网格层析方法局部修饰速度模型。本文改进了传统层析速度建模流程：首先利用炮域波动方程相干反演方法建立初始速度模型,再利用角度域共成像点道集进行层析迭代,最后利用基于层位约束的网格层析方法局部修饰速度模型。渤海M区块三维实际数据试算证明本文对层析速度建模流程的改进是有效的,改进后的新流程可以大幅提高层速度场的精度。     

自动拾取的成像空间域走时层析速度反演

本文研究了一种成像空间域走时层析速度反演方法,利用波动方程双平方根算子叠前深度偏移提取的角度域共成像点道集(ADCIGs)作为速度分析道集,并基于剩余曲率自动拟合策略获取高精度的走时差。文中不仅给出ADCIGs道集剩余曲率自动拟合拾取的主要思路,而且给出成像空间域层析反演的实现。模型和实际资料试算验证了该方法具有较高的反演精度和计算效率。     

正交各向异性介质多参数建模方法研究

<正>交各向异性处理方法在各向异性特征明显的复杂裂缝地区的应用效果良好,但正交各向异性介质多参数建模方法尚不成熟。首先针对正交各向异性介质特性,通过分方位速度分析确定裂缝发育方向、提取方位各向异性参数,并采用融合公式建立正交各向异性参数初始模型;然后建立方位正交各向异性层析矩阵,制定先分方位、再联合的多参数反演策略;最后建立描述地震波速度方位特性的正交各向异性参数精细模型。三维模型数据试算和实际资料处理结果表明,所建立的模型参数分布合理,与真实地质情况吻合度高,有效改善了基于此模型的OVT域成像道集拉平效果,进而提升了断裂和缝洞储层成像效果,为准确预测断裂以及AVO精细分析提供了可靠的成果数据。     

初至波与反射波旅行时多尺度渐进联合层析成像

针对应用单-波反演难以更准确、更全面地重构地下速度信息及层析正反演中网格剖分大小相互矛盾的问题,本文采用Humphreys等^[1]的数据加权算法,对模型空间和数据空间进行加权,并在反演中采用多尺度渐进层析反演;通过比较模型上各种波的理论旅行时与实际拾取的旅行时残差,更新模型速度及界面深度参数。模型资料试算结果表明,应用初至波与反射波旅行时联合层析整体上要优于单-波旅行时层析的结果;多尺度渐进层析反演能有效地解决层析正反演网格剖分大小的矛盾;在-定范围内,随着反演网格尺度的逐渐减小,反演精度也会逐渐提高。实际资料试算结果表明该方法是可行的,可以为高精度地震勘探提供精细的速度模型。     

基于角道集剩余曲率分析的层析速度建模

中深层地下构造由于地震波场复杂、覆盖次数不足、信噪比低及有效反射较弱等原因造成速度建模精度低,进而影响了深度偏移质量。为此,本文提出一种中深层层析方法:利用角道集作为速度分析道集,推导了基于剩余曲率的层析建模方程,并将角道集的剩余曲率沿射线路径反投影进行速度更新。该方法既避免了常规反射层析在低信噪比叠前数据上同相轴拾取困难的问题,又无须在射线追踪中进行试射以拟合炮检点坐标,具有较高的反演精度。模型试算和实际数据处理结果均验证了方法的有效性。     

高斯波束角道集的旅行时层析速度分析

本文综合高斯波束叠前深度偏移准确、灵活、高效及其角道集能够直接反映成像点偏移深度与局部入射角之间的关系,避免多波至带来的假象干扰等优势,提出了基于高斯波束角道集的旅行时层析速度分析方法,推导了灵敏度矩阵和角道集剩余曲率之间的定量关系式,并利用自动拟合策略求取剩余曲率。模型和实际资料试算验证了该方法具有较高的反演精度和计算效率。     

基于岩石物理的页岩储层各向异性表征

针对罗家地区A井测井数据,应用Backus平均方法将模型从测井尺度粗化至地震尺度,应用各向异性传播矩阵模拟页岩含油储层的地震AVO响应,并结合实际地震资料进行井震标定;基于VTI介质反射系数理论公式进行VTI介质参数AVO反演,并应用井中岩石物理反演结果进一步计算页岩含油储层的裂缝空间分布。结果表明：①随着模型中各向异性参数ε数值的增大,远炮检距反射振幅逐渐减小;②由地震反演得到的各向异性参数和水平缝密度与岩石物理的计算结果一致,由反演得到的水平缝的空间分布与油气产量具有相关性,证明了反演方法的有效性。     

相似系数阈值滤波的数据驱动控制束偏移

高斯束偏移不仅具有接近波动方程偏移的精度,而且具有Kirchhoff偏移灵活、高效的特点。然而当实际地震采集数据中含有较强噪声时,易产生偏移假象而影响成像质量。为此,在传统高斯束偏移的基础上,根据有效信号和干扰信号在τ-p域中的相干性差异,发展了一种基于相似系数阈值滤波的数据驱动控制束偏移方法。采用数据驱动策略,在高斯束偏移成像过程中,先计算τ-p道集的相似系数,再通过设定相似系数阈值控制干扰信号,从而降低偏移剖面中的随机噪声;控制束偏移可以直接提取角度域共成像点道集,无需复杂的角度映射变换且具有更高信噪比。模型测试及实际资料处理结果表明：对于低信噪比数据,控制束偏移剖面的信噪比明显高于常规高斯束偏移,但会损失相对振幅的可靠性;尽管控制束偏移在τ-p道集的滤波过程增加了一定的计算量,但总体与常规高斯束偏移方法的计算效率相当;相似系数阈值参数选取十分关键,阈值较小时偏移噪声较强,但过大阈值也可能压制部分有效信息或产生偏移假象,选取合适的阈值参数才能得到较理想的偏移剖面。     

波动方程叠前深度偏移角度道集求取方法

针对常规直接转换法和直射线叠前时间偏移(PSTM)法提取的角度道集无法满足地质构造复杂地区的AVA分析与反演精度的问题,研究基于共炮检距道集波动方程保幅叠前深度偏移方法,提取角度域共成像点道集。首先,对地震数据进行叠前深度偏移得到炮检距域共成像点道集;其次,在频率-波数域对炮检距域共成像点道集利用快速插值映射法抽取角度域共成像点道集。通过与直接转换法和直射线PSTM法提取的角度道集比较,分析不同方法求得的角度道集信息的准确性和AVA特征。数值试算和实际资料处理结果表明:共炮检距道集波动方程保幅叠前深度偏移方法求取的角度道集的振幅相对保真,角度范围更广,更适用于AVA分析。     

TI介质角度域高斯束逆时偏移方法

高斯束逆时偏移兼具高斯束偏移灵活、高效和逆时偏移高精度的优势,具有面向目标成像的能力。本文将基于相速度的各向异性射线追踪算法引入高斯束逆时偏移中,并结合高斯束计算时的传播角度信息,实现了一种更为高效的TI介质角度域高斯束逆时偏移方法。模型试算表明：同传统基于弹性参数的各向异性算法相比,在保证成像精度的前提下,本文方法具有更高的计算效率,同时提取的角度域共成像点道集（Angle Domain Common Imaging Gather,ADCIG）不仅能够为后续偏移速度分析提供支撑,而且可以用于分角度叠加成像,压制成像噪声、提高成像质量。     

VTI介质角度域叠前深度偏移

研究了VTI介质角度域偏移方法,以各向同性双平方根方程的角度域偏移方法为基础,从VTI介质qP波频散关系出发,推导出VTI介质角度域偏移的波场延拓算子;在频率-波数域处理以横向均匀速度传播的波场,在空间域处理具有速度扰动特征的波场以提高波场延拓精度。模型试算和实际资料处理结果表明：各向同性偏移方法由于未考虑各向异性参数的影响,绕射波不能完全收敛,波场聚焦效果差,降低了成像剖面的分辨率和信噪比,不能对地质构造精确成像;VTI介质角度域偏移可对断层、盐丘、小尺度地质体精确成像。对于角度域偏移产生的角度域共成像点道集（ADCIG）而言,各向同性偏移的ADCIG同相轴无法校平,残留断点绕射波,波场无法正确聚焦,不能正确反映局部地质特征;VTI介质偏移的ADCIG同相轴较平直,角度范围更宽,波场归位准确,精度较高。因此,利用VTI介质角度域偏移方法可对复杂构造精确成像。     

VTI介质角度叠加逆时偏移

常规各向同性逆时偏移方法在处理各向异性问题时会导致层位成像不准确、同相轴能量不聚焦等问题,并且利用拉普拉斯滤波算子等简单滤波方法不能彻底去除成像结果中的强振幅低频噪声。针对这些问题,首先应用声波近似方程模拟qP波在各向异性VTI介质中的传播,在震源激发位置设定平滑的各向同性或椭圆各向异性震源环以压制qSV波的影响;在此基础上提出了VTI介质逆时偏移的角度叠加实现策略,利用角度域共成像点道集中不同角度成像结果的叠加压制低频噪声,提高成像精度。模型试算和应用实例偏移结果表明,该方法能够对高陡复杂构造准确成像,且能量均衡、低频噪声压制效果较好。     

基于成像域的多次波压制

常规多次波压制方法大都是在数据域基于滤波理论或者波动理论实现的,这些方法都不能有效地兼顾多次波压制精度和计算效率。为此,本文将波动理论与滤波理论相结合,基于成像角度域共成像点道集(ADCIGs)实现多次波压制。首先对成像域多次波映射进行了讨论,通过推导可知,成像域ADCIGs中多次波和一次波剩余时差不同,这一结论即为在成像域应用滤波方法压制多次波的理论基础;考虑到基于DSR方程的叠前深度偏移可以获得更多的角度信息,然后基于DSR方程提取ADCIGs,再引入适用于ADCIGs的正切平方估计高分辨率拉东变换分离一次波和多次波,实现多次波压制。由于本文方法将波动理论与滤波理论相结合,因此能同时兼顾多次波压制精度和计算效率,模型试算和实际资料处理验证了方法的有效性。     

弹性波逆时偏移的角度域共成像点道集提取方法

对于多分量地震勘探,由于介质中纵、横波的耦合作用,多波角道集存在能量串扰和低波数噪声干扰等问题。针对横波Poynting矢量构建过程中存在的应力串扰问题,基于解耦延拓方程及伪横波应力构建方法,对质点振动速度和应力进行解耦,实现了无能量串扰的纵波和横波Poynting矢量构建;在基于Poynting矢量提取PP波反射角和PS波反射角的基础上,提出了基于弹性逆时偏移的PP波和PS波角道集提取方法;针对PP波角度域共成像点道集存在的低波数噪声干扰问题,提出了角度域共成像点道集的角度衰减叠加策略,提高了偏移成像精度。应用模型数据验证了方法的正确性。     

基于角度域共成像点道集的层析速度反演方法研究

碳酸盐岩探区的复杂构造成像与速度反演是目前地震数据处理难点问题。碳酸盐岩构造复杂、目标层极深、断点发育,常规处理方法很难得到比较好的成像与速度反演结果。为此提出基于角度域共成像点道集(ADCIGs)、面向碳酸盐岩储层的层析速度分析方法,利用波动方程双平方根算子叠前深度偏移提取的角度域共成像点道集作为速度分析道集,能够更为准确地反映速度和深度的耦合关系,并且减少了假象等的干扰,使得速度分析结果精度更高。模型和实际结果表明,该方法能够较好地解决碳酸盐岩探区的速度反演问题。