地震绕射波弱信号U-net网络提取方法

绕射波携带大量小尺度非均匀地质体信息,对于提高地震勘探分辨率具有重要意义.绕射波能量远小于反射波,在地震记录中常被强反射波掩盖,因此分离并单独成像绕射波,为探测小尺度非均匀地质体的关键问题.传统绕射波分离方法受限于理论模型假设,对陡倾角反射波去除效果不佳,且易对绕射波造成损伤.基于经典编码-解码框架下的U-net网络和注意力机制,本文提出了一种绕射波智能分离方法,通过编码器自动提取地震数据中的绕射波特征,再由解码器恢复绕射波,从而隐性去除反射波.该方法作为端到端的机器学习,训练后的U-net网络可自适应地分离绕射波.本文通过数值模拟数据与实际数据构建训练数据集,利用训练后的U-net网络分离绕射波,并将结果偏移成像.数值模型测试和实际资料应用表明,融合了注意力机制的U-net网络能够有效压制反射波能量,保留绕射波动力学特征,克服了传统绕射波分离方法难以去除陡倾角反射的局限性,其提取的绕射波弱信号特征较为完整,能够进一步提高地震成像分辨率,在小尺度断裂刻画上具有优势.     

地震绕射波分离方法研究进展

随着油气田勘探和开发的需求,利用地震绕射波成像探测地下不规则地质体边界和小尺度构造已经成为地震偏移成像的一个热点研究方向.从勘探地球物理学家发现地震绕射波开始,学者们致力于研究如何在常规反射地震记录中分离并提取出绕射波.过去半个多世纪,关于绕射波的研究取得了显著进展并成功应用到实际生产.本文总结了以往关于绕射波分离技术的研究成果,介绍绕射波的概念及其运动学和动力学特征,按照技术方法进行总结和分类并通过选择具有代表性的实例来向读者展示不同方法的分离效果.     

局部倾角滤波和预测反演联合分离绕射波

地震绕射波源于介质非连续性,从地震记录中将绕射波分离出来并进行成像,其结果对研究诸如碳酸盐岩缝洞储层这类复杂非均质储层具有重要意义.对炮集记录进行平面波分解,在地层倾角不大的假设下,反射波和绕射波同相轴在平面波分解剖面上存在较大的倾角差异.基于此,我们提出分步进行绕射波分离的方法:(1)利用局部倾角滤波方法将绕射波的较大倾角信息成分分离出来,此时,余下的部分包含有反射波和残留的低倾角绕射波信息;(2)利用频率-空间域预测反演方法从上述含有反射波和残留的低倾角绕射波信息中分离出残留绕射波成分;(3)将两次分离的绕射波信息相加得到最终的绕射波估计.用该方法能够得到相对完整的绕射波信息,有效地克服了靠单一的倾角差异进行绕射波分离时明显损失低倾角信息,从而影响绕射波成像结果横向分辨率这一问题.理论与实际资料试算验证了该方法的有效性.     

基于机器学习的时间反转散射体检测方法

地下小尺度散射体的检测和识别对于提高地震勘探的分辨率具有重要意义,目前业界普遍采用绕射波分离及成像方法检测地下散射体,而绕射波成像的质量主要取决于绕射波和反射波波场分离的程度.本文将被动源震源定位问题中常用的时间反转原理引入到地下散射体检测中,首先通过分析被动源和主动源模型反传波场的聚焦状态,验证了时间反转原理应用于地下散射体检测中的可行性;并引入机器学习中的朴素贝叶斯分类算法,给出适用于时间反转散射体检测的分类算法框架,计算模型中每个点成为散射体的概率,最终检测出地下散射体最有可能存在的位置.散射体模型和Sigsbee2a模型的试算结果证实了本文方法在不需对反射波和绕射波分离的情况下,即可实现对地下散射体的检测和定位,同时由于考虑了多次散射的影响,检测结果能准确反映地下散射体的位置.     

基于倾角-偏移距域道集的绕射波成像

地震绕射波是地下非连续性地质体的地震响应,绕射波成像对地下断层、尖灭和小尺度绕射体的识别具有重要的意义.在倾角域共成像点道集中,反射波同相轴表现为一条下凸曲线,能量主要集中在菲涅耳带内,绕射波能量则比较发散.由于倾角域菲涅耳带随偏移距变化而存在差异,因此本文提出一种在倾角-偏移距域道集中精确估计菲涅耳带的方法,在各偏移距的倾角域共成像点道集中实现菲涅耳带的精确切除,从而压制反射波.在倾角-偏移距域道集中还可以分别实现绕射波增强,绕射波同相轴相位校正,因此能量弱的绕射波可以清晰地成像.在倾角域共成像点道集中,反射波同相轴的最低点对应于菲涅耳带估计所用的倾角,因此本文提出一种在倾角域共成像点道集中直接自动拾取倾角场的方法.理论与实际资料试算验证了本文绕射波成像方法的有效性.     

绕射波的一组新数据

1985年3月3日智力中部城市互尔帕来索西部近海发生了一次7.6级地震,牡丹江地震台763长周期地震仪记录了一组绕射波(表)。经计算,震中距为158.9°,绕射波走时为18分19.5秒。为便于比较,根据《震相走时便查表》求出了绕射波走时方程: Tpdif=383.3+4.4Δ (33公里),r=1     

基于反稳相滤波的边缘绕射成像方法研究

绕射波是地下小尺度不连续体的地震响应,是对地下断层、尖灭点、河道边界等高精度成像的关键.为了充分利用绕射波信息,本文通过在传统的Kirchhoff偏移算子中引入反稳相滤波器,使满足Snell定律的镜面反射得到有效压制,从而凸显绕射能量;进一步,针对边缘绕射存在的极性反转现象,在成像过程中引入校正极性算子实现边缘绕射正确成像.在实现方法的基础上,通过理论模型和实际资料边缘绕射成像试处理发现:本文方法能有效改善断点及断层边界成像质量,提高地震解释精度;通过极性校正,能一定程度提高边缘绕射点成像分辨率.     

平面波预测滤波分离绕射波方法研究（英文）

常规地震资料处理通常将地震波场当作反射波处理,绕射波的存在会降低处理成果的分辨率,因此将其当做干扰噪音进行压制。然而,绕射波往往包含了非常重要的地质信息(裂缝、溶洞、断层等特殊构造),即使偏移将绕射波归位,因为绕射波能量较弱也会被反射同相轴遮盖。将反射波和绕射波进行分离并单独成像处理,将提高绕射目标构造的成像精度。根据在平面波域反射波拟线性、绕射波拟双曲的形态差异,本文利用平面波预测滤波技术对平面波炮记录进行波场分离。首先,采用平面波解构滤波器估算地震同相轴的局部倾角,利用局部倾角信息预测并提取反射波,然后从全波场中减去反射波,间接分离出绕射波,对绕射波成像可获得高分辨率的绕射目标体成像结果。通过对2D SEG盐丘模型进行波场分离与成像,分析表明:本文采用的平面波预测滤波技术解决了平面波解构滤波的极性反转问题,保持了绕射波原有的相位信息,保证了波场分离后的绕射波场保真度,成像结果表明用平面波预测滤波方法提取的绕射波进行偏移可改善非均质目标体的成像精度。     

频率域中的合成地震记录

波动理论的合成地震记录是地震资料处理和解释中的重要研究课题。 目前,制作波动理论的合成记录主要是用Trorey的克希霍夫方法,正在研究的有有限差分法和有限元素法。本文提出一种新的方法,在频率域中制作波动理论的合成记录,给出了与克希霍夫绕射波方程等价的频率域中的绕射波方程。在国产DJS-11型计算机上,应用两种绕射波方程实现了合成地震记录,并且进行了比较,认为频率域方法是有效的,可以用于波动方程偏移和地震资料解释,进一步工作有可能获得多层变速模型的波动理论记录。     

论折射法勘探中的绕射波射线轨迹

根据超声波三维地震模拟实验的结果,研究了绕射波射线轨迹的问题。对绕射波的运动学特征进行了讨论,并引用极值光程律进行了解释。实验和理论的分折指出,绕射点的位置经常按照一定的规律偏离震源--接收点联线的一侧。 本文还提出了“双曲线叠加法”,从而在实际较复杂的情况下,可以用它迅速地确定绕射点的空间位置及波的旅行时间。 文中提出绕射波的运动学特征,从中可以提供一些有关断层或绕射棱结构的情况。至于绕射波的动力学特征及机制将另文阐述。所获的结果及结论有助于物探资料的解释。     

用相关分析法提高波速测量精度

本文尝试利用相关分析法，计算回波之间的相关系数，以确定反射波之间的走时，从而提高波速测量精度，实测表明，这种方法是可行的。     

火山岩地震屏蔽层的转换波叠前时间偏移成像

在反射地震转换波资料处理中,准确求取共转换点一直是一个难题,采用叠前时间偏移技术能避免共转换点道集的抽取,而且能够使转换波归位到真正的反射点上,实现准确成像.本文针对火山岩地震屏蔽层的转换波成像问题,通过对转换波共近似转换点道集进行速度分析,建立了转换波叠前时间偏移的初始速度场,通过速度扫描和纵、横波速度比值扫描确定最佳的偏移速度场和纵、横波速度比值,实现了在火山岩高速层覆盖区域的转换波偏移成像.实际资料的成像结果表明,本文采用的近似转换点计算以及转换波叠前时间偏移方法是有效的.     

南海南部深部结构的复杂构造地震成像

南海南部海底崎岖、地下构造复杂、二维地震资料中多次波非常发育,在地震数据处理中很难获得令人满意的成像效果,为此本文采用了预测反褶积后F-K域滤波及叠前时间偏移后Radon变换去除多次波的组合方案,有效地衰减了多次波,突出了有效信号能量;采用速度分析点加密法获取了更为准确的速度,实现了陡倾角构造与海山的准确成像;采用百分...     

非对称型强飓风中的准平衡流特征分析Ⅱ:波动结构

利用非对称波分量分解和小波分析的方法,对准平衡动力模型下非对称强飓风内中尺度波动的空间结构和时间序列特征进行分析.结果表明,平衡流场中1波扰动占主要地位且具有涡旋Rossby波的典型结构特征,准平衡流各波数下扰动的空间分布反映了中尺度波动的混合性质;模式大气和准平衡垂直运动的全局功率谱中,超过信度检验的强波动信号中不仅包含分别表征重力波和涡旋Rossby波的高频和低频波动信号,还存在表征具有物理性质不可分特性的混合涡旋Rossby-重力波的中频波动.混合波的出现建立了不同频段波动之间的能量交换通道,其信号的变化对飓风系统的强弱变化具有一定的指示作用.非平衡垂直运动的波动强信号则主要集中在高频和低频区域,反映了在飓风强度变化情况下,与高频重力波有关的快波调整过程所引起的垂直扰动的振荡和传播.强垂直风切变对飓风内中尺度波动的切向和径向传播具有重要影响,当环境垂直风切变很强时,准平衡1波扰动在径向和切向方向上均呈"驻波"形态,随着环境垂直风切变的减弱,1波扰动以混合波波速逆基本气流传播.     

Effects of seismic wave propagation upon buried pipelines

The behaviour of long straight buried pipelines subjected to seismic wave propagation is investigated. Well-known relationships for determining upper bounds for the axial strain and curvature in the pipeline as well as relationships for relative displacement and rotation at the pipeline joints are discussed. The assumption that the seismic excitation can be modelled as a travelling wave having a shape which remains unchanged as it traverses the pipeline is examined in detail. It is shown that this assumption is unconservative when the effective propagation velocity of the seismic waves with respect to the pipeline is such that the actual time lag (separation distance between points divided by effective propagation velocity) is less than a ‘cross-over’ time lag. Cross-over time lags for 22 pairs of ground displacements recorded during the 1971 San Fernando Earthquake are presented in this paper. Finally, methods for estimating the propagation speed of the seismic waves along or with respect to the pipeline are discussed.     

RESEARCH ON TIME-DISTANCE INTERPRETATION OF MULTIPLES AND GHOSTS BASED ON MARINE SEISMIC DATA

Currently, the study on the active fault in the land areas is relatively mature, while there is still lack of detection and research on active faults in the sea areas. Marine exploration, which is different from land areas, has a prominent problem due to the existence of strong reflecting interfaces such as water surface and seafloor in the sea, thus the recording is often accompanied by interference of multiples on seafloor reflections. In addition, because of the characteristics of marine seismic exploration, the source exciting in the water and the geophone receiving in the water, ghost wave usually can be recorded simultaneously during the reflected wave propagation. This phenomenon makes it difficult to distinguish the effective waves and the noise, and has always plagued the data and seriously affects the quality of records. In the offshore and other regions of complex structures, such as inclined interfaces, it is difficult to eliminate the interference of multiples accurately by traditional multiples suppression methods, which are based on the horizontal interface assumption. This paper combines the sea area seismic data and its acquisition method, uses simplified model to simulate the multiples based on the time-distance analysis of multiples and their ghost wave in inclined interface. The time-distance characteristics of the multiples and their ghost waves from different interfaces(including the inclined interface)are obtained, and they are consistent with the actual records. The multiples time-distance simulation can help to distinguish the causes of reflected waves, summarize the multiple-wave time-distance characteristics from different interfaces(including inclined interfaces), and analyze the relationship between the characteristics of multiple waves and primary waves. In particular, this simulation has a significant effect on characterizing the internal multiples that are difficult to identify due to inconspicuous periodicity and the multiples of the inclined interface which present the phenomenon that the vertex of the time-distance curve is shifted. On this basis, relying on the time-distance analysis of ghost wave, we analyze the travel time difference characteristics between reflected waves and their accompanying ghost waves. The differences of the travel time characteristics of different orders ghost wave and reflected wave are summarized and the symmetry of the travel time difference between inclined interface and horizontal interface of ghost waves and reflected waves is analyzed. We simulate the distraction of the ghost wave event with the event of the reflected wave and analyze the influence of the ghost wave on the sea area seismic records. These results can improve the practical interpretation of seismic data. At last, the time-distance information is used to synthesize sea area seismic records, which can help us carry out the effective data processing and understand the characteristics of the time-distance and velocity of multiples in different interfaces and the layer artifact caused by multiples. This study combines the time-distance simulation of multiples and their ghost wave with conventional seismic data processing to analyze the pre-stack and post-stack features of multiple waves and their ghost waves in the seismic records of the sea area. The results of this study are conducive to the effective identification of multiples in seismic records in the sea, provide a theoretical basis for multi-wave suppression and prediction, and may facilitate the future study of sea-area seismic activity detection.     

Analysis of effects of active sources on observed phase velocity based on the thin layer method

In the wave field induced by active sources, the observed phase velocity of surface waves is influenced by both mode incompatibility (i.e. non-planar spread of surface waves is idealized as plane waves) and body waves. Effects of sources are usually investigated based on numerical simulations and physical models. Several methods have been proposed to mitigate the effects. In application, however, these methods may also have difficulties since the energy of the body waves depends on soil stratification and parameters. There are multiple modes of surface waves in layered media, among which the higher modes dominate the wave field for soils with the irregular shear velocity profiles. Considering the mode incompatibility and the higher modes, we derive analytical expressions for the effective phase velocity of the surface waves based on the thin layer stiffness method, and investigate the effects of the body waves on the observed phase velocity through the phase analysis of the vibrations of both the surface waves and the body waves. The results indicate that the effective phase velocity of the surface waves in layered media varies with the frequency and the spread distance, and is underestimated compared to that of the plane surface waves in the spread range less than about one wavelength. The oscillations that appeared in the observed phase velocity are due to the involvement of the body waves. The mode incompatibility can be ignored in the range beyond one wavelength, while the influence range of the body waves is far beyond one wavelength. The body waves have a significant influence on the observed phase velocity of the surface waves in soils with a soft layer trapped between the first and the second layers because of strong reflections.     

天然地震频率范围内首都圈地区近地表S波速度结构c

近地表沉积层的S速度结构是强地面震动模拟和地震灾害估计的重要参数,尤其是浅部的S波速度结构在工程上具有重要的应用意义．目前大部分资料来源于工程钻孔或工程地震探测,很少有地震波频率范围内的S波速度结构,或者深度达数百米的S波速度结构．通过对天然地震的井下摆波形记录的分析,提供了一种测量地震波频率范围深达数百米的S波速度的有效方法．收集了首都圈地区44个井下摆的近震记录,利用广义射线方法确认了直达S波及其在地表的反射波震相,并通过测量不同台站上两个震相的到时差,获得了首都圈地区浅层100——500m 深度范围的S波速度结构．研究发现,浅部100m 的平均S波速度低于300m/s．当深度增加到500m 时S波速增加到800m/s,平均速度梯度为0.8 （m/s）/m．研究结果表明,井下摆地震记录波形是研究沉积盆地浅层S波结构的重要资料,将为沉积盆地的强地面震动模拟提供重要基础参数．      

Feasibility of waveform inversion of Rayleigh waves for shallow shear-wave velocity using a genetic algorithm

Conventional surface wave inversion for shallow shear (S)-wave velocity relies on the generation of dispersion curves of Rayleigh waves. This constrains the method to only laterally homogeneous (or very smooth laterally heterogeneous) earth models. Waveform inversion directly fits waveforms on seismograms, hence, does not have such a limitation. Waveforms of Rayleigh waves are highly related to S-wave velocities. By inverting the waveforms of Rayleigh waves on a near-surface seismogram, shallow S-wave velocities can be estimated for earth models with strong lateral heterogeneity. We employ genetic algorithm (GA) to perform waveform inversion of Rayleigh waves for S-wave velocities. The forward problem is solved by finite-difference modeling in the time domain. The model space is updated by generating offspring models using GA. Final solutions can be found through an iterative waveform-fitting scheme. Inversions based on synthetic records show that the S-wave velocities can be recovered successfully with errors no more than 10% for several typical near-surface earth models. For layered earth models, the proposed method can generate one-dimensional S-wave velocity profiles without the knowledge of initial models. For earth models containing lateral heterogeneity in which case conventional dispersion-curve-based inversion methods are challenging, it is feasible to produce high-resolution S-wave velocity sections by GA waveform inversion with appropriate priori information. The synthetic tests indicate that the GA waveform inversion of Rayleigh waves has the great potential for shallow S-wave velocity imaging with the existence of strong lateral heterogeneity.     

高斯波包反射走时速度反演方法

扰动高斯波包理论指出,在Gabor域描述模型的扰动成分,且入射波场为短时宽带信号时,扰动波场可在时间域通过高斯波包算子描述.在此基础上通过拟合反射波的走时,提出一种速度反演方法.反射波走时残差利用地震道局部波形的互相关函数表示,以走时残差的二范数作为目标函数,优化目标函数实现对速度场的反演.基于一阶Born近似,利用扰动高斯波包理论推导出目标函数对速度场的梯度是本文理论部分的核心内容.梯度包括两部分:正传的背景波场与反传的扰动高斯波包之间的互相关,反传的背景波场和正传的扰动高斯波包之间的互相关.梯度表达式中背景波场和扰动波场均利用高斯波包算子模拟.计算梯度的具体算法中,如何模拟扰动波场,以及如何计算反射波的走时残差是两个要点,文中对此做了详细的讨论.数值实验进一步阐述了反演的实现策略,实验结果表明高斯波包反射走时速度反演方法和实现策略有效可行,并得到了理想的反演结果.