SRME多次波衰减方法在海洋地震资料中的应用

多次波问题一直是海洋地震资料处理中的难题,基于波动方程的多次波压制方法SRME是近20年来发展起来的一种有效的多次波衰减方法.该方法的特点是基于波动理论,利用地震数据本身来预测多次波,而不需要地下介质的任何信息,即这种方法可以应对地下任何复杂的情况,因此应用范围较广.本文通过SRME方法在海洋地震资料LH油田的应用实例证明了该方法的有效性--SRME多次波衰减方法对海底有关的多次波的衰减非常有效,特别是对近道多次波模拟准确,未损伤有效波,一定程度上提高了地震数据的分辨率.     

地震多次波成像

地震资料含有各种类型多次波,而传统成像方法仅利用地震一次反射波成像,在地震成像前需将多次波去除.然而,多次波携带了丰富的地下结构信息,多次波偏移能够提供除反射波外的额外地下照明.修改传统逆时偏移方法,用包含一次反射波和多次波的原始记录代替震源子波,将SRME方法预测的表面多次波代替一次反射波作为输入数据,可将表面多次波成像.多次波成像的挑战和困难在于大量串扰噪声的产生,针对表面多次波成像中的成像噪声问题,将最小二乘逆时偏移方法与多次波分阶思想结合起来,发展可控阶数的表面多次波反演成像方法,有望初步实现高精度的表面多次波成像.在消除原始记录中的表面多次波后,通过逆散射级数方法预测得到层间多次波,将层间多次波作为逆时偏移方法的输入数据可将其准确归位到地下反射位置.数值实验表明,多次波成像能够有效地为地下提供额外照明,而可控阶表面多次波最小二乘逆时偏移成像方法几乎完全避免成像噪声.

     

福山凹陷浅表层火山岩对地震资料品质的影响浅析

根据福山凹陷实际钻井地层和地震剖面情况,建立地质模型,通过对二维模型的正演计算和照明分析计算,分析因浅表层、深层火山岩而产生的多次波对地震资料品质的影响程度,并提出针对性的措施。     

自适应虚同相轴方法压制地震层间多次波

地震数据中发育的层间多次波是影响速度分析和偏移成像的精度和可靠性的关键.通常情况下,层间多次波的动校正量、叠加速度和频率与一次波并无明显差异,从而对识别、预测和压制多次波带来了极大挑战.传统虚同相轴方法基于物理图像和定性公式,其预测的层间多次波振幅和相位精度难以满足实际需求,造成了其对匹配算法的过度依赖.本文针对传统虚同相轴方法的理论缺陷和计算精度问题,通过理论推导得到了新的自适应虚同相轴方法.相比于传统方法,自适应虚同相轴方法能够显著提高压制多次波能力,同时减少对匹配算法的依赖.本文给出了自适应虚同相轴方法的推导过程,并运用一维和二维模型算例验证了方法相较于传统虚同相轴方法的多次波预测精度优势.通过在PLUTO模型和实际陆地地震数据上的应用实例,证明了本文新研究的自适应虚同相轴方法对去除层间多次波,恢复并突出目标储层同相轴,提高地震成像分辨率的显著作用.

     

地震波衰减机理及能量补偿研究综述

随着中、浅层油气勘探程度的不断提高,深层油气藏的勘探已经成为当今勘探研究的主要对象,而深层油气勘探要获得突破关键是用于综合研究的地震剖面要有较高的信噪比和分辨率,深层与浅层能量均衡.因此,国内外对地震波衰减机理和补偿方法进行了大量研究.本文综述了基于岩石物理和散射理论的地震波的衰减机理,重点概括了地震波衰减的影响因素,指出衰减与孔隙度、频率、温度、压力、岩性等因素有关,最后分类介绍了各种地震波衰减的能量补偿方法.     

基于构建虚同相轴压制地震层间多次波

基于构建虚同相轴来估计层间多次波的方法是克希霍夫积分表示定理的一个延伸发展.本文通过构建虚地震同相轴巧妙地将散射点从地下移到了表面,利用表面的散射点来预测层间多次波.由于预测与实际地震记录中的层间多次波振幅存在偏差,采用了多道的L1范数匹配算法来实现真振幅的多次波压制.对简单模型和南黄海地质模型的数据处理实例验证了本文研究方法的正确性和有效性.     

基于地震剖面瞬时峰值能量频率的衰减分析与气藏预测方法

由于流体本身的粘滞性和摩擦性,使得地震波在含油气地层中传播时衰减形成低频阴影.通过对叠后地震资料做时频分解,分析地层衰减特征,是当前除了AVO技术之外,比较流行的直接识别流体和气藏的技术.本文基于波动方程正演,模拟了地震波对含油气地层的响应.然后采用连续小波变换分析不同尺度剖面的瞬时能量特征,提取分频剖面上每一个点的峰值能量对应的频率,得到瞬时峰值能量频率剖面.通过分析瞬时峰值能量频率剖面,分析了地震波在地下介质中传播时的衰减特征.最后应用此方法对实际地震剖面进行了分析,验证了其有效性.     

利用基于T-K算子的非线性能量衰减分析技术进行储层含油气检测

以地震资料为主体的储层含油气性检测技术在油气田勘探开发过程中发挥着重要的作用.基于含油气储层可以造成地震波高频能量衰减剧烈而低频能量衰减较弱的特性,结合小波变换与T-K算子对地震波不同频率能量进行非线性跟踪,通过定量表征地震记录高频非线性能量的衰减梯度,为储层含油气检测提供有效信息.基于T-K算子的非线性能量是不同频率...     

核电厂地震安全性评价衰减关系影响分析

地震危险性分析中的不确定性处理和表征,一直是核电厂厂址地震安全性评价中倍受关注的重要问题,尤其是日本福岛核事故后,无论是确定核电厂厂址的设计基准地震动,还是进行核电厂地震风险评价,都更加重视地震危险性分析中的不确定性.本文通过理论分析重点说明了衰减关系的不确定性,包括标准差和截断水平对核电厂地震安全性评价的影响,并在此基础上,通过算例和讨论说明了概率性方法截断水平的选取问题,探讨了现行确定性方法和概率性方法在截断水平选取上的差异.分析计算结果表明,在地震活动较弱的区域,概率性方法截断水平为3,确定性方法截断水平为0的现行做法是恰当的.但是,对于发震构造大震复发间隔较小的区域,为了使二者在超越概率方面协调,恰当提高确定性方法的截断水平更为合理.     

异常地层压力声波衰减分析及地层压力地震数据外推反演

澳大利亚西北大陆架卡拉汶盆地含有丰富的油气资源, 但其地质条件复杂, 普遍存在异常地层高压分布,特别是该区区域盖层-中深层巨厚的Muderong页岩层内存在局部异常高压带,给石油钻井带来巨大困难.本文根据卡拉汶盆地北部地区10余口井的实测地层压力资料、声波测井数据和联井地震剖面, 研究该页岩层内异常压力千赫兹尺度声波响应的衰减特征,及其与地层有效压力之间的关联计算模型.通过联井地震数据的井-震相关性分析进行井-震过渡,研究井旁地震数据在异常高压带的衰减特征,及其与地层有效压力之间的关联计算模型,从而实现井震资料联合进行地层压力地震外推反演.研究表明,Muderong页岩层地层有效压力与井中声波速度/阻抗品质因子和井旁地震阻抗品质因子具有良好的正向相关性,根据二者交汇数据拟合经验公式反算得到的各井有效压力结果基本反映出实测有效压力曲线的变化特征.卡拉汶盆地北部地区地震数据具有良好的井-震相关性,地震资料品质相对可靠,确保了声-压相关分析的井-震过渡.最后,通过地震阻抗品质因子与地层有效压力之间的关联计算模型,进行该区异常地层压力地震外推反演,预测结果与钻井证实的异常压力分布特征基本符合.本研究为卡拉汶盆地中深部区域盖层局部异常高压带的地层压力外推反演提供了一种有效途径.     

深地震反射剖面构造信息识别研究

随着我国深部探测计划的实施,近年来在中国大陆地区实施了多条深地震反射剖面,获得了大量揭示地壳上地幔精细结构的宝贵资料. 为了更充分地挖掘这些资料中所包含的深部构造信息,本文在研究和总结前人工作的基础上,提出一种快速识别深地震反射剖面构造格架的方法. 该方法通过数据预处理、强振幅提取、中值滤波、对象识别、连续性计算和连续性滤波实现对深地震反射剖面构造格架识别. 同时,该方法还可通过对象倾角计算对复杂区域进行属性分析.     

地震波衰减的物理机制研究

对近年来有关地震波衰减的物理机制进行评述和初步研究认为,对于地壳介质,造成地震波衰减的主要原因是地壳内存在大量裂隙,裂纹中饱含水或部分含水,地震波传播时引起裂隙中的流体运动,从而造成地震波的衰减,对于上地幔,扩散控制的位错阻尼机制可能是造成地幔地震波衰减的一个重要原因,而上地幔软流层部分熔融的存在也是造成地震波衰减不可忽视的原因．     

基于深度神经网络的地震强反射剥离方法

在储层预测工作中,储层弱反射信号淹没在强反射信号之中的情况非常常见,这不利于精确识别和描述储层结构.本文提出了一种基于深度神经网络的强反射剥离方法,用于辅助储层弱反射信号的检测工作.该方法在卷积模型的框架下将强反射预测问题分解为地震子波预测与强反射预测两个子优化问题,并采用AIDNN与U-Net两个深度神经网络分别求解...     

双相介质中地震波衰减的物理机制

High-frequency seismic attenuation is conventionally attributed to anelastic absorption. In this paper, I present three studies on high-frequency seismic attenuation and propose that the physical mechanism results from the interference of elastic microscopic multiple scattering waves. First, I propose a new theory on wave propagation in a two-phase medium which is based on the concept that the basic unit for wave propagation is a nano- mass point. As a result of the elasticity variations of pore fluid and rock framework, micro multiple scattering waves would emerge at the wavelength of the seismic waves passing through the two-phase medium and their interference and overlap would generate high- frequency seismic attenuation. Second, I present a study of the frequency response of seismic transmitted waves by modeling thin-layers with thicknesses no larger than pore diameters. Results indicate that high-frequency seismic waves attenuate slightly in a near-surface water zone but decay significantly in a near-surface gas zone. Third, I analyze the seismic attenuation characteristics in near-surface water and gas zones using dual-well shots in the Songliao Basin, and demonstrate that the high-frequency seismic waves attenuate slightly in water zones but in gas zones the 160-1600 Hz propagating waves decay significantly. The seismic attenuation characteristics from field observations coincide with the modeling results. Conclusions drawn from these studies theoretically support seismic attenuation recovery.     

Scattering and anelastic attenuation of seismic energy in northern Greece

The relative contribution of scattering (Q _s ^–1 ) and intrinsic (Q _i ^–1 ) attenuation to the totalS-wave attenuation for the frequencies of 1.5, 3.0, 6.0 and 12.0 Hz has been studied by applying the radiative energy transfer theory, Data of local earthquakes which occurred in northern Greece and were recorded by the permanent telementered network of the Geophysical Laboratory of the University of Thessaloniki have been used. The results show that in this area the scattering attenuation is dominant over all frequencies while intrinsic attenuation is significantly lower. The estimatedQ _s ^–1 andQ _i ^–1 values have frequency dependences off ^–0.72 andf ^–0.45, respectively. The frequency dependence ofQ _s ^–1 is the same as that of the codaQ _c ^–1 , obtained by applying the single scattering model, which probably implies that the frequency dependence of the coda wave attenuation is attributed to the frequency dependence of the scattering attenuation.Q _c ^–1 values are very close to scattering attenuation for short lapse times, (10–20 sec), and intermediate between scattering and intrinsic attenuation for the longer lapse times, (50–100 sec). This difference is explained as the result of the depth-dependent attenuation properties and the multiple scattering effects.     

呼和浩特—包头盆地岩石圈细结构的深地震反射探测

跨呼和浩特-包头盆地(以下简称"呼包盆地")完成的91.8 km长的深地震反射剖面,揭示了呼包盆地的岩石圈精细结构和断裂的深、浅构造特征.结果表明,本区地壳和岩石圈具有清晰的层状反射结构特征,其中,地壳厚度约45~48 km,岩石圈厚度约82~87 km.莫霍面在大青山之下出现约3.5 km的抬升,暗示大青山的隆升不是因为地壳物质增厚所致,即大青山可能不存在"山根".呼包盆地为南浅、北深的"箕状"断陷盆地,盆地沉积层最厚处位于大青山山前,其厚度约为7~8 km.鄂尔多斯北缘断裂和大青山山前断裂作为呼包盆地的南、北边界断裂,在剖面上均表现为由3~4条断裂组成的"Y"字形断裂构造,它们对呼包盆地的形成、地层沉积、基底变形和地震活动都有重要的控制作用.剖面揭示的岩石圈深断裂位于大青山山前断裂的下方,该断裂向上进入上地壳,向下切割中-下地壳、莫霍面,进入上地幔.深断裂的存在为深部热物质的上涌与能量强烈交换提供了通道,而上涌的软流层物质与岩石圈地幔发生交代和侵蚀作用导致岩石圈减薄.     

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

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

     

北京地区地壳精细结构的深地震反射剖面探测研究

长度100 km、NW向穿过三河—平谷8.0级地震区和北京地区主要断裂构造的深地震反射剖面,揭示了该区地壳精细结构图像和断裂的深浅构造特征.结果表明,该区地壳以TWT6～7 s左右的强反射带为界分为上地壳和下地壳,上地壳厚约18～21 km,下地壳厚约13～15 km.剖面TWT3～4 s以上,反射层位丰富,构造形态清晰,且在剖面上具有明显不同的构造特征;在三河—平谷地震区以西,剖面揭示了2～3组反射能量较强的反射震相和一系列错断基底面的断裂,在三河—平谷地震区以东,为一套自东向西倾伏的密集强反射层,这套反射具有典型的沉积盆地特征,盆地最深处约为8～9 km.剖面揭示的地壳深断裂倾角陡直,该断裂切割、扰动了下地壳物质和壳幔过渡带,向上延伸至上地壳,将地壳深部构造与浅部断裂联系在一起,构成了该区最主要的深浅构造特征.     

基于随机介质模型的统计分析在深反射地震中的应用

深反射地震勘探是揭示地球深部构造的一种地球物理探测方法.由于不同于沉积地层的反射特征,深反射地震资料的解释需要不同于浅层水平层状介质的解释方法.基于随机介质模型的统计分析为深反射地震解释提供了新的途径.本文首先介绍了地球深部结构中非均质性的统计特性研究成果,阐述了随机介质模型的基本概念及二维随机介质的建立方法,讨论了随...