简介：坐标变换法通过将物理空间的曲网格映射为计算空间的矩形网格,将起伏地表转化为水平地表,同时将物理空间的波动方程转化为计算空间的波动方程,在计算空间完成数值模拟,坐标变换的方法对处理起伏自由边界具有较好的适应性和应用效果。本文在传统坐标变换方法的基础上,根据计算区域速度差异采用不同的网格大小和采样时间步长,提出了一种基于时空双变网格的起伏地表坐标变换正演模拟方法。在编程实现算法的基础上,通过典型模型波场模拟试算结果分析可知：（1）变网格方法与常规方法波场模拟误差在0.5%左右;（2）变网格方法计算效率视不同的变网格区域面积及变网格大小可提高几倍量级,在本文模型和计算参数下提高约5倍。（3）在满足模拟精度及频散条件要求下,变网格方法较全局细网格算法能显著节约计算内存。为此,针对起伏地表数值模拟,本文方法具有较高的模拟计算精度和一定的适应性。

简介：逆时偏移被认为是处理陡倾角断层、回转波和逆掩推覆构造最有效的成像方法。考虑到当今勘探对地下复杂构造体、薄互层和隐蔽储层的高分辨率成像需求，井间、井中地震资料成为地面地震资料的一个很重要的补充。但由于地面地震、井中地震和井间地震存在很大的频率差异，本文将相应的变网格算法应用到偏移中。另外，也将Lanczos滤波引入到波场外推计算中，从而可以很好地消除变网格引起的人为反射。试算表明，对比地面、井中和井间的偏移，可以看到微构造、陡反射层都得到了更好地刻画和描述。基于变网格的逆时偏移在储层预测和岩性解释方面有广阔的应用前景。

简介：地震波场数值模拟方法对理解和分析地震波的传播规律具有着重要的意义。弹性波动方程能够模拟地下介质的实际情况，为偏移和成像提供有效的依据。在弹性波波场数值模拟中，旋转交错网格数值模拟（RSM）修改了标准交错网格数值模拟（SSM）方法，将同类的参数定义在同样的节点上，拓宽了稳定性条件的约束，但在低速区会出现较严重的频散。变阶数差分方法是自适应空间算子长度方法的一种变化和推广。它以理论频散误差研究为基础，结合实际波场传播的情况进行误差计算，对不同速度匹配不同的差分阶数。本文研究了变阶数旋转交错网格数值模拟（VRSM），即是籽变阶数方法应用到RSM中，它可以很好地解决RSM在低速区域的数值频散问题，以及减少不必要的时间损耗；同时讨论了旋转交错网格的理论频散特性，并基于波场分离的方法分析了实际波场传播的频散误差，将原方法的应用范围由声波推广到剪切波，由理论值推广到时变值。在数值模拟试验中，VRSM将被应用于水平层模型和Overthrust模型。通过阶数分配以及相应波场传播效果和计算时间的分析，验证了该方法应用于复杂介质波场模拟中的实用性和有效性。实验的结果表明VRSM能够合理分配不同速度所对应的差分阶数，能保证计算的精确性，并合理控制计算的时间。

简介：Gabor变换和S变换是常用的时频分析工具。根据测不准原理,它们的时频分解结果无法在时间域和频率域同时具有很高的分辨率。为了提高非平稳信号时频分解结果的分辨率,本文提出瞬时频率分布函数（IFDF）并利用它表达非平稳信号。当非平稳信号时频成分的分布满足测不准原理对信号可分辨的要求时,瞬时频率分布函数的支集和短时Fourier变换的小波脊支集是同一个集合。利用IFDF的该特征,本文提出一种迭代算法（Sparse-STFT）实现了信号的稀疏时频分解。该算法在每次迭代过程中利用残留信号的短时Fourier变换结果的脊支集更新信号的时频成分,每次迭代得到的时频成分的叠加结果即为最终的稀疏时频分解结果。文中的数值实验证明了Sparse-STFT可以有效地提高非平稳信号时频分解结果的分辨率。最后,本文将该方法应用于地震数据面波的压制中,取得了理想的处理结果。

简介：尾波是持续时间比主要波型长得多的多重散射波，它包含了直达波之外的部分有用信号。在油气田的开采过程中，由储层物性参数的微小变化而导致的储层速度的微小变化对初至的影响很小，无法通过初至变化直接观测，但是因为尾波为多重散射波，故储层速度的微小变化会在尾波的传播过程中被放大，因此利用尾波可以观测到这种变化。本文通过实验与数值模拟研究了尾波干涉方法在储层微小时移差异监测方面的有效性。合成地震数据是基于部分MarmousiII模型采用有限差分波动方程正演方法计算获得，为了模拟时移地震中的储层波速变化，在模型中选取波速发生微小变化的目标区，计算速度扰动前后的合成地震数据，进而观测尾波变化。研究表明利用尾波干涉可检测储层物性的微小变化，为时移地震信号监测提供新的方法和手段，从而提高时移监测的准确度，也为开发生产工作提供可靠的指导。

简介：通过地震数据获取裂缝储藏中流体的性质并对流体类型进行识别,是地震勘探岩性反演的重要问题之一。由于地震波的速度、储层的密度等弹性参数对某些流体不具有很强的敏感性,使只依赖振幅信息进行流体识别的传统AVO方法面临困境。作为传统叠前振幅反演的一个拓展,频变AVO（FDAVO）技术进一步考虑了振幅对频率的依赖关系,将这种依赖关系与地下裂缝结构、流体填充对应起来,能带来更丰富的流体信息。利用该技术,本文提出了一种基于地震数据参数化Chapman模型的贝叶斯反演新方法（BIDCMP）,它包含两步算法,即,FDAVO反演储层的非弹性属性和贝叶斯框架下的流体识别。首先,通过匹配观测数据和模型数据,构造差函数反演裂缝储层非弹性参数。随后,在贝叶斯框架下,使用马尔科夫随机场（MRF）作为先验模型,联合多参数场识别流体。本方法在计算过程中,除综合考虑了弹性参数场、测井资料等常规信息外,还特别地加人了第一步中反演得的非弹性参数的约束,从而充分利用了流体粘性差异,最后在最大后验概率（MAP）准则下输出最佳岩性一流体识别结果。分别对合成地震记录和模拟岩性—流体剖面验证本文方法的有效性,结果证明本文方法获得的流体识别结果准确可信。

简介：提高地震资料的分辨率是地震数据处理流程中的重要环节,对后续的精细构造解释起到重要作用。传统的提高分辨率方法大都假设地震资料是稳态的并且噪声水平不随空间发生变化,而实际情况不满足这一假设,导致提高分辨率处理后的效果达不到预期要求。针对这一问题,本文提出了一种基于时频二次谱的提高地震资料分辨率方法。首先,文中提出了基于S变换的时频二次谱,并结合模型论述了时变子波和反射系数在时频二次谱中的特征及其可分离性;其次,依据时变子波和反射系数在时频二次谱中的特征差异,构建了二维滤波器在地震记录的时频二次谱中提取时变子波的振幅谱;再次,文中研究了噪声环境中时变提高分辨率算子设计方法,并提出了依据时频谱能量强弱相对关系自适应确定频带拓宽范围的时变提高分辨率算子设计,进行提高分辨率处理;最后,文中对该方法进行了模型和实际数据的试处理,并与传统谱模拟方法和Q补偿方法的处理结果进行了对比分析,对比结果表明：本方法不需要估计Q值,提高分辨率能力不受震源子波频带的限制,在兼顾信噪比的前提下能够充分提高不同时间局部的地震数据的分辨率。