混合Radon变换地震噪声压制的应用

如何有效地从地震资料中剔除相干噪声一直是地震勘探十分关心的问题,尤其是多次波的压制和消除。笔者采用相似函数压制端点效应和截断效应的混合Radon变换从地震剖面中一次性分离面波等线性噪声,再应用自适应滤波技术在Radon域识别并剔除多次波,并使用双曲Radon反变换计算获得有效反射地震信号。理论模型及实测资料试算结果表明,该方法对多次波压制效果良好,同时,应用混合Radon变换可简化数据中噪声压制处理流程。     

基于稀疏反演算法的高分辨率Radon变换及其在多次波压制中的应用

针对目前双曲Radon变换Radon域能量团能量泄漏的问题,提出一种基于稀疏约束反演算法的高分辨率Radon变换进行多次波的压制. 该方法在常规双曲Radon变换的基础上加入稀疏约束条件,应用共轭导向梯度( CGG)反演算法进行求解,能有效地提高Radon域速度谱的聚焦效果,有利于分离一次波和多次波. 模型和实际资料试算证明了该方法的有效性和实用性.     

基于低频约束的高分辨率Radon变换多次波压制方法研究

抛物线Radon变换是地震资料处理中常用的一种多次波压制方法,常规Radon变换由于采用最小二乘算法的原因,对不同的曲率值采用相同的加权值进行求解,分辨率较低。为提高Radon域数据的分辨率,应采用对不同曲率值应用不同的加权值的方法。一般的思路是利用前一次迭代的结果得到加权矩阵,该算法需要迭代进行,运算量较大。这里拟采用低频约束的方式提高Radon变换的分辨率,即利用前一个频率的运算结果对下一个频率的计算进行约束,该算法分辨率高,计算速度快。模拟数据和实际数据的测试表明,本文方法在多次波压制处理中,可以快速有效地完成多次波去除。     

基于F-K和Radon变换的多次波衰减方法

多次波干扰是地震勘探中最重要的问题之一,其衰减技术也是地震数据处理的难题之一。本文主要介绍了多次波的形成、分类及其特征,并在 ProMax 系统上对合成地震记录的多次波进行了压制:首先采用合适的速度对时-空域的地震记录进行 NMO 校正;然后,分别用 F -K 变换和 Radon 变换方法将其变换到 F -K 域和 Radon域,对多次波所占据的主要区域进行充零衰减;最后,通过对应的反变换将其变换到时-空域,从而达到多次波衰减的目的。通过在 ProMax 系统中的应用表明:不同方法对多次波的衰减效果各有优劣,组合滤波方法能达到更好的效果。     

952南区多道地震多次波压制方法

多次波的压制是震资料处理中的难点和重点。特别是海上地震资料记录时间长,多次波发育。本文针对南沙海域地震地质条件。提出二维皮方法及其辅助办法,详细了二维Ｆ－Ｋ滤波的原理及演变过程,并通过对南沙海的实际资料自理得出了满意的结果。     

各向异性 Ｒadon 变换在深水地震多次波压制中的应用

在深海多缆、长拖缆地震勘探的各种噪音影响中,多次波对深层的有效反射影响难以消除。笔者介绍了一种高分辨率各向异性 Ｒadon 变换方法,在 Ｒadon 变换积分路径中增加非椭圆率 eta 参数,将传统 Ｒadon 变换积分曲线由慢度、偏移距两参数改为慢度、偏移距及非椭圆率三参数,更好地刻画了长偏移距下反射非双曲时差。通过对某海域实际采集的多缆长偏移距地震数据进行压制多次波处理,与传统方法相比取得了较为理想的效果,为后续深层构造成像提供了基础数据。     

线性Radon域地震干涉层间多次波预测方法

为避免时空域地震干涉技术预测层间多次波过程中出现假象干扰,影响地震资料处理。笔者结合线性Radon变换与地震干涉技术,在线性Radon域中利用地震干涉技术预测层间多次波。通过数值模拟验证了线性Radon域地震干涉层间多次波预测方法的有效性,并将该方法应用到复杂速度模型中。结果表明,线性Radon域地震干涉技术可以有效减少时空域地震干涉技术预测层间多次波存在的假象,压缩数据体积大小,提高计算效率,并能够适用于复杂地质环境。     

高分辨率抛物线拉冬变换多次波压制技术

对于孔径过小,采样不充分的地震数据,抛物线拉冬变换的多次波压制效果不理想,且存在空间假频问题。综合Sacchi、Mauricio、Todd Mojesky等对抛物线拉冬变换压制多次波的分析,阐明了高分辨率抛物线拉冬变换原理。通过理论模型试验和实际地震数据的分析证明,高分辨率抛物线拉冬变换可以突破离散采样和有限孔径的限制,实现对地震数据中的多次波及空间假频的快速、有效压制。     

基于地震信号相似性的S变换域多次波压制技术

多次波压制是海洋地震数据处理的关键,如何保护一次反射波的同时高效准确地识别并压制多次波是地震资料处理的重点。笔者在共中心点(CMP)叠加压制多次波方法的基础上,基于CMP道集地震信号与叠加道地震信号的相似性特征在S变换域进行多次波压制。首先,利用多次波速度动校正并对CMP道集数据减去叠加道数据完成第一次多次波减除,此时,波场中残余的多次波以随机噪声的形式存在;然后,对减除多次波的地震波场应用一次波速度动校正,并以一次波速度叠加道为初始模型对CMP道集S变换谱进行相似性滤波完成第二次残余多次波去除。针对理论数据和实际数据进行了有效性测试,结果表明本文所提方法能够较好地压制多次波,提高反射波的成像精度。     

基于Shearlet变换的非局部均值地震噪声压制

在地震勘探中,由于野外地震数据采集环境及仪器性能本身的限制,采集到地震信号中不可避免地会混入较强的噪声,极大影响后续处理、解释工作。而近几年,多尺度几何分析方法以其独特优势成为压制噪声的研究热点,本文提出在Shearlet域中引入非局部均值算法对地震噪声进行压制,该算法首先对地震信号进行非下采样Shearlet变换,然后采用非局部均值法对分解后系数子集进一步处理,并采用8个Sobel算子近似表示全方向结构,对权重函数进行改进,最后对系数进行Shearlet反变换,得到去噪后的地震信号。实验结果表明相比于传统非局部均值法,该联合算法能有效地压制随机噪声,同时对弱同相轴具有更好的保护作用,在地震资料处理中具有良好的实用性。     

Attenuation of multiple waves using the continuous wavelet transform (CWT)

This paper describes how the continuous wavelet transform is used to filter multiple waveforms in both time and frequency domains. It is well suited to process the stationary signals, and it shows the signal in both time and frequency scales. This new approach was tested first on synthetic data and then on real data. The results obtained on both cases were good. The method consists of identifying the multiples on which we apply a normal move out using the multiple velocity law. The multiples will be aligned and the primary reflections will not be aligned. This operation allows locating the multiples in the time-scale domain. We compute the continuous wavelet transform (CWT for short) in order to focus on the patterns relative to seismic events. To filter the multiples, we define a zone with frequency and time bounds. These bounds are deduced from the projection of the seismic trace. Then an automatic mask is applied to the pattern to be isolated. Filtering in time–frequency domain is done by keeping only the wavelet coefficients that are outside the mask and assigning zero to the coefficients larger than a threshold amplitude inside the defined zone. The mask shape does not matter, which is not the case in classical filtering, where both the window size and shape play a key role. The mask is defined from three parameters: time, frequency, and the wavelet coefficients. To go back to the time domain, one has to compute the wavelet transform inverse of the trace. This procedure is repeated for all traces. To reset the traces to their initial positions, we apply the dynamic correction inverse with the same velocity law as the multiples. It turns out that the attenuation of multiples by the CWT works fine, in particular, the two identified multiples were quasi eliminated (Fig. 10).     

对比抛物Radon正变换几种矩阵求解方法

抛物Radon变换法(Parabolic Radon Transform)在地震资料处理中有广泛的应用。PRT可对不同频率的地震数据解耦处理,这一特点使得抛物Radon变换的计算效率比双曲Radon变换有数量级上的提高。在频率域求解时,需要对每一个频率成份求解同样大小的线性方程组。求解抛物Radon正变换的计算方法主要有Levinson递推法、共轭梯度法、Cholesky分解法和直接矩阵求逆法。最小平方抛物Radon正变换所形成的矩阵具有Toeplitz结构,可采用Levinson递推法进行计算。高分辨率抛物Radon正变换所形成矩阵的Toeplitz结构被破坏,一般采用共轭梯度法或Cholesky分解法进行求解。这里详细推导了复Toeplitz矩阵的Levinson递推算法,并分别对求解方程的四种方法进行了讨论,最后给出抛物Radon正变换求解的数值算例,并对所给出的四种方程求解方法的计算效率及计算精度进行了对比。     

抗噪性加权抛物线Radon变换的地震道重建

叠前CMP道集中,走时为双曲线的地震反射同相轴经过部分动校正( NMO)之后变为近似抛物线,当道集中缺失某些地震道时,可利用抛物线Radon变换通过迭代计算进行地震道的重建和恢复. 笔者基于抛物线Radon变换的基本原理,利用加权抛物线Radon变换的计算方法,进行叠前地震道的恢复和重建,并且在数据域对其加权系数进行改进和优化,模型数据和实际资料的试算结果表明,改进后的计算方法即使在较强的随机噪声下依然可以稳建同相轴,且避免了常规方法在重建过程中产生虚假同相轴的缺点,与原始的加权抛物线Radon变换的计算方法相比,新方法具有压制随机噪声的优点,即具有较好的抗噪性,计算精度也有所提高,且算法稳健、高效,具有良好的应用前景.     

Evaluation of geodynamic activity of the Dead Sea transform fault by radon gas concentrations

Twelve radon lines of dosimeters (detectors) were placed across the Jordan Valley active fault, which is a segment of the active Dead Sea transform fault system. Each line of the dosimeters shows one or more peaks of radon anomaly concentrations. Some of these peaks prove the intersection of the fault trace with these lines in areas where the fault plane is inferred. In other lines, the peaks correspond to the arrangement of faults in areas of pull-apart basins (sag ponds) or pressure ridges, formed due to the left or right step of the fault. Sag ponds usually show low radon emanations, because they are the place for the accumulation of very fine sediments, which decreases their porosity and hence the upward migration of the radon gas. The northern part of the Jordan Valley relatively shows high radon emanation, which could be attributed to the presence of a seismic gap in the upper Jordan valley.     

基于L_1范数的多次波自适应减方法及应用分析

多次波问题在地震勘探中普遍存在。自由表面相关多次波压制(SRME)方法是目前多次波压制方法的主流,使用该方法的重要步骤之一是将由反馈迭代法预测得到的地震多次波经匹配后从原始数据中减去。基于L2范数的多次波自适应减方法有其适用范围,仅在某些情况下才有好的处理结果。这里基于迭代重加权最小二乘法(IRLS算法)的混合L1/L2范数来近似L1范数解,同时结合模型数据和实际数据进行多次波压制处理,并与基于L2范数的自适应减方法进行对比分析。结果显示,本方法不仅有效压制了多次波,而且还相对更好地保持了有效波的能量,这表明本方法可以在不同情况下实现更为普遍的多次波压制。     

一种土壤氡浓度测量的校正方法及其应用

通过对土壤氡浓度测量校正模型的应用和分析,初步建立了根据残留气体的土壤氡浓度校正方法,并在实际应用中论证了方法的可参考性及实用性。     

利用小波变换确定电磁测深曲线拐点的应用研究

电磁测深数据是地下地质体对电磁波传输特性的综合反映。由于地下情况的复杂性和地质体对电磁波特性反映的迟钝性,往往由电磁波携带的有用信息非常弱小和复杂,因此,要想通过正演对比和建模反演的方式来准确解释电磁测深资料和正确建立反演模型是相当困难的。电磁波在均匀介质中传播时,其幅度按指数规律单调递减,而随着地下地质体的变化将会改变这一规律。通过寻找测深曲线的拐点,就可以确定地下地质体的基本形态和分界情况,从而为资料定性解释和反演建模提供参考。然而在实际应用中,要直接求出观测曲线的拐点是不可能的。小波变换方法被称为“数学显微镜”方法,用小波变换来处理观测信息后,在不同尺度的小波变换结果中,分析结果是一种非常实用可行的方法。