基于门控循环单元主导的深度网络半监督地震衰减补偿方法

地震波在地下传播时存在能量衰减,呈现非平稳性,导致地震资料的分辨率降低,不利于储层刻画.为了解决这个问题,需要对地震记录进行衰减补偿来消除非平稳性.考虑到GRU(Gated Recurrent Unit)即门控循环单元在处理时间序列时具有长时和短时记忆的优势,本文构建了以GRU为主导的网络结构(简称GRU网络), 提出了基于GRU网络的地震衰减补偿方法.实际地震数据处理中,测井资料往往非常有限,导致标签数据较少.为解决小样本问题,本文借鉴自编码器的思想,除井旁道和测井合成地震记录的标签外,将其余地震记录作为无标签数据引入训练,半监督地学习非平稳数据到平稳数据的非线性映射,实现了地震记录的衰减补偿.最后利用含噪合成地震记录和实际地震资料测试了本文提出的方法,证明了其有效性.

     

基于广义Curvelet变换和相空间特征聚类的VSP波场分离方法

垂直地震剖面(Vertical Seismic Profiling, VSP)蕴含大量的地层地质信息, 是连接地震反射数据与测井资料之间的桥梁.VSP数据分辨率高, 资料丰富, 包括上行波、下行波、转换波等多种类型的波场, 常被用于地层反射界面标定、介质衰减参数反演等.其中, 如何有效地分离波场是将VSP资料应用于油气勘探的关键之一.广义Curvelet变换是一种具有多尺度、多角度的相空间变换, 能够进行相空间域波场特性的提取.本文聚焦于VSP波场的上下行波分离问题, 提出一种基于广义Curvelet变换和相空间特征聚类方法.首先, 采用广义Curvelet变换对VSP波场进行特征提取, 构造融入相空间角度信息的特征数据集; 然后, 利用K-means聚类算法在相空间对波场特征进行分类; 最后, 对分类结果进行反变换, 完成VSP上下行波波场的自适应分离.为了验证方法的有效性, 本文将所提出的方法用于合成数据和实际VSP数据的波场分离, 并与常用的基于F-K变换的波场分离方法进行对比.处理结果表明, 本文方法角度分辨率高、抗噪能力强, 波场分离结果的保真保幅性好, 这为后续的成像与地层的特征分析提供了重要基础资料.

     

基于多特征融合自编码器的无监督地震相分类研究

地震相分类是地震数据解释中的一个重要步骤, 是地震数据与沉积相的连接工具.为了提高地震相分类精度和减少对有限人工标签的依赖, 本文提出了一种基于多特征融合自编码器的无监督地震相分类方法.首先, 提出了一种混合卷积和变分编码的多特征融合自编码器, 实现了地震数据中表征地震相的大量隐含特征提取.其次基于非负矩阵分解和K均值聚类实现了主特征分量分解和地震相聚类.实际地震数据应用结果和指标分析表明, 本文方法提取的隐含特征趋于正态分布, 且主特征分量中蕴含了不同地震相类别的响应, 从而可以获得更准确的地震相分类结果.在渤海湾盆地东营凹陷古近系沙河街组湖相沉积中, 清晰划分出了六类沉积微相的边界, 有利于揭示三角洲沉积环境演变.

     

基于目标区的局部角度域方向照明分析及其在偏移成像中的应用

叠前深度偏移已成为近年来地震偏移成像领域广泛应用的技术之一,其中复杂盐丘下方目标区域结构的成像始终是偏移成像中较难解决的问题.局部角度域波场分解可以提供在空间和方向上的双重局部化信息,因而这一技术被广泛地应用在方向照明分析、成像振幅校正等方面.本文在采用局部指数标架小波束进行角度域方向照明分析的基础上,研究采集系统对复杂盐丘下部层状结构成像质量的影响;同时通过分析目标区对应于相对地表采集系统分布的照明能量,确定叠前炮集数据中对盐下目标区成像起重要贡献的部分,并将这些数据用来对目标区成像,从而达到提高盐下结构成像质量的目的.本文将二维SEG-EAGE盐丘模型中盐下反射体作为目标区,分别计算其照明分析和由照明能量分布确定的炮集数据对目标区的偏移成像结果,同时通过与成像振幅校正前后的成像结果对比,说明了该方法对提高目标区结构成像质量的有效性.     

基于时空局域化dreamlet单程波算子的观测系统沉降法偏移

Dreamlet偏移成像目的是探索一类能够对地震波场和单程波传播算子同时分解和压缩的理论和方法,也即实现在压缩域的传播与成像.地震数据在时间和空间的非平稳性质,决定了要实现地震数据的有效稀疏表示,分解方法必须在时间和空间上同时具有局域化性质.Dreamlet由时间和空间局部分解原子的张量积构成,可以看作一种脉冲-小波束形式的波场分解原子.时空局域化的dreamlet单程波传播算子在对波场沿深度方向延拓时,地震数据在时间轴上总是向同一方向流动.随着深度的增加,部分用于成像浅层结构的数据归位至其空间位置后被dreamlet算子丢弃,波场的有效记录时间变短,每一步用于波场延拓的计算量也相应下降.为了充分发挥这一优势,本文介绍dreamlet传播算子与观测系统沉降法偏移相结合的理论与方法.观测系统沉降法偏移每一步都将记录到的所有数据向下延拓,沉降后的波场等效于把源和检波器都放置在目标深度所能接收的反射数据.Dreamlet观测系统沉降过程只保留用于成像观测系统下部地质结构的有效数据,自动丢弃已经用于成像观测系统上部而对下部成像没有贡献的信号.本文通过二维SEG/EAGE叠后和Marmousi叠前数据算例展示了dreamlet传播算子应用于观测系统沉降法偏移的这一特点.数值算例结果显示,在不影响成像质量的前提下,该偏移方法能够有效减少传播数据量,为发展一种快速高效的偏移方法提供了新的思路.

     

广义S变换域有色噪声与信号识别方法

为研究广义S变换域有色噪声与信号识别问题,首先推导白噪声在S变换和广义S变换下平均功率谱的统计规律,白噪声的平均功率谱与频率呈线性关系. 然后从理论上定性分析一般的有色噪声功率谱的分布特性,得出在广义S变换域有色噪声的功率谱服从自由度为2的χ2分布的规律, 并用Monte Carle方法进行了验证. 在上述工作基础上,本文提出了利用噪声和有效信号统计特性的差别,在广义S变换域中区分噪声和信号的方法,并用模型数据证明了方法的有效性.     

利用VSP资料直达波的包络峰值处瞬时频率提取介质品质因子

本文采用一个具有4个待定参数的函数去逼近震源子波,利用黏弹介质中单程波传播理论推导出了地震子波包络峰值处瞬时频率(EPIF)和品质因子之间的解析关系;同时,为提高瞬时频率的估计精度和抗噪性能,在小波域中发展了一种计算瞬时频率的方法,并在此基础上提出了估算VSP资料Q值的方法,简称小波域包络峰值处瞬时频率法(WEPIF)．合成的VSP数据衰减估计结果表明,与对数谱比法和中心频率偏移法相比较,WEPIF法受到界面反射波影响相对较小、计算结果稳定、精度相对较高．将WEPIF法用于某气田的单炮零偏6级VSP资料Q值估计,结果表明,吸收强弱与储层含气性高低有良好的对应关系．     

基于熵谱特征定量计算薄层厚度的方法研究

讨论了基于地震资料的熵谱特征定量计算薄互层厚度的原理。首先对谱均衡后的地震数据作归一化处理,然后对处理结果作最大熵谱估计,若估计结果的谱分辨率精度不能满足要求,则先对归一化的地震记录叠加适当的白噪声,然后再作最大熵谱,这样可以得到较高的估计精度。另外,文中对流程中关键步骤之一的最大熵谱估计进行了改进,并给出了基于相关技术用于从谱上提取周期的方法,提高了谱周期的估计精度。文中用合成地震数据证明了其有效性。     

叠前地震资料噪声衰减的小波域方法研究

讨论了叠前地震记录中随机噪声的模拟方法,采用定性分析与Monte Carlo仿真相结合的方法,研究了在物理小波为基本小波的变换域该类随机噪声局部小波功率谱的统计特性,得出：其局部小波功率谱服从χ22分布.在此基础上,利用假设检验给出了在给定置信度下确定有效信号在小波变换域内能量分布空间的方法,通过在有效信号能量分布空间重建信号来压制噪声.给出了模型及实际地震资料算例,证明了本文提出的方法的有效性.     

地震资料处理中小波函数的选取研究

本文给出了常见地震子波的一个模拟公式,可以很好地模拟零相位及混合相位子波,在一定意义上也可以近似模拟最大相位及最小相位子波．模拟出的子波加上适当的修正项后满足允许条件,可用作小波函数．与Ｍｏｒｌｅｔ小波类似,在实际应用中这些修正项在一定条件下可以略去,文中对Ｍｏｒｌｅｔ小波作了改造,使其能更好地适应于地震资料处理．研究了反射波能量及噪声等干扰波在时间－尺度域的分布特征与所选基本小波的关系．提出用地震子波（或与地震子波相近的函数）作为基本小波,对地震资料进行去噪及分频解释的方法．最后用实例证明方法的有效性．