首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到18条相似文献,搜索用时 125 毫秒
1.
Radon变换是一种稀疏变换,被广泛应用于地震数据处理,其中线性Radon和抛物Radon最为常用.在实际地震数据中,直达波和面波的同相轴形态为线性,反射波为双曲型,单独使用线性Radon或抛物Radon变换时,不能确保所有同相轴在变换域的系数都是稀疏的,影响地震数据处理效果.本文提出的多路径Radon变换联合了线性Radon变换和抛物Radon变换,每个同相轴都有两种不同的变换参数来与积分路径相适应,能够兼顾不同形态的同相轴;然后利用最小二乘稀疏反演方法对不同形态同相轴匹配最佳积分路径,使其自动分离到两个不同的Radon域剖面且保持系数同时稀疏.从多路径Radon域剖面上,能够很容易地识别不同系数所对应的同相轴形态、时间截距以及速度,这些特征有利于提高利用Radon变换方法进行随机噪声压制、面波压制以及波型分离等技术的处理效果,该变换在模型数据和实际数据中的应用结果证明了本文方法有效性.  相似文献   

2.
熊登  赵伟  张剑锋 《地球物理学报》2009,52(4):1068-1077
高分辨率Radon变换存在计算效率和分辨率不能兼得的困境.时间域算法可以获得很高的分辨率,但计算效率非常低;频率域算法具有良好计算效率,但分辨率不理想.为此发展了混合域高分辨率抛物Radon变换,即对频率域抛物Radon变换引入时变的稀疏权.本文给出了一种新的混合域高分辨率抛物Radon变换实现方法,并将该算法应用于叠前数据衰减多次波.文中给出了Radon变换和衰减多次波的流程.理论和实际数据算例表明本文方法既有较高的分辨率又有很高的计算效率.  相似文献   

3.
多次波的存在会降低地震资料的信噪比,影响地震资料处理效果和后续的地质解释精度,压制多次波干扰是地震资料处理中的一个重要环节。利用有限差分方法正演模拟几个典型地质模型的地震波场,之后进行动校正,再通过抛物线Radon变换方法把t-x域一次反射波和多次波变换到τ-q域,由于存在速度的差异,Radon域记录中的一次反射波和多次波的能量互相分开,在Radon域对多次波的能量切除,保留反射波的有效信息,达到压制多次波的效果,提高资料的信噪比。通过对实际地震数据进行抛物线Radon变换结果也进一步验证了该方法在压制多次波中的应用效果。   相似文献   

4.
多方向正交多项式变换压制多次波   总被引:4,自引:3,他引:1       下载免费PDF全文
提出一种基于Radon 变换和正交多项式变换的多方向正交多项式变换压制多次波方法.抛物Radon变换对不同曲率方向的同相轴叠加,根据速度差异区分一次波和多次波,但Radon反变换会损伤振幅特性,不利于AVO分析.多方向正交多项式变换在Radon变换(某一曲率方向的零阶特性)的基础上,利用正交多项式变换进一步分析同相轴的高阶多项式分布特性,用正交多项式谱表征同相轴AVO特性;根据一次波和多次波速度差异和同相轴能量分布特征实现多次波压制.该方法的优点是仅用一个曲率参数就可描述同相轴剩余时差参数,提高了一次波和多次波的剩余时差分辨率.实验结果表明,该方法可以有效压制多次波并保留一次波AVO特性.  相似文献   

5.
多次波压制是地震数据处理的重要环节之一,而Radon变换方法以其优势被广泛应用于多次波处理环节中.本文提出利用时间延迟量Γ代替传统抛物Radon变换的曲率参数q,对地震记录同相轴进行积分运算.当使用时延Radon变换时,输入数据的道数作为重要参数直接影响Radon变换域的假频,而非传统抛物Radon变换空间方向的采样间隔.使用预条件算子通过重建间接增加输入数据的道数,消除Radon域的假频问题,提高成像精度,进而更好的实现多次波压制.利用本文方法对水平层状介质和复杂模型数据进行试算,结果表明,Radon域假频得到有效抑制,多次波压制效果明显.  相似文献   

6.
本文针对井间和3D VSP波场的线性特征,研究井孔地震波场线性高分辨率Radon变换算子,用于井孔地震波场分析与纵横波分离.在Radon变换原理分析基础上,采用基于柯西分布的高分辨率线性Radon变换对井孔数据进行Radon变换,其间通过对离散倾角叠加算子求取的研究,及对影响Radon能量收敛的重要参数阻尼因子算法的改进,使数据在Radon域以能量团的形式呈现,得到很好的收敛效果,基本解决了Radon域数据的一定程度的拖尾现象,消除了各能量团之间的平滑效应,采用柯西分布来规则化数据,提高了Radon域的分辨率,Radon域能量也收敛到一个点上,有利于上下行波或纵横波波场分离.最后通过反演结果和模型试算验证了该方法的可行性和稳定性.  相似文献   

7.
石颖  王维红 《地球物理学报》2012,55(9):3115-3125
基于波动方程预测的表面多次波压制方法可处理复杂地下介质的地震资料,但计算成本较高.基于滤波的多次波压制方法计算效率较高,但其成功应用仅局限于一次波和多次波有明显时差差别的地震数据,对来自速度逆转等复杂介质数据则较难获得满意的压制效果.本文将波动方程预测的反馈迭代法和滤波法有效结合,采用GPU(图形处理器)和CPU协同并行加速计算粗略预测表面多次波,随后在双曲Radon域比较分析原始数据和预测的多次波,设计合理有效的Butterworth型自适应滤波器,滤出原始数据Radon域中的多次波能量,进行Radon反变换后,在时空域将多次波从原始数据中减去,得多次波压制结果.文中对理论模拟的单炮数据、复杂的SMAART模型以及实际地震数据进行了计算,结果表明,结合基于波动方程预测和双曲Radon变换的方法有效突破了两种方法各自的局限性,可高效高精度地压制复杂地下介质的表面多次波.  相似文献   

8.
基于稀疏反演三维表面多次波压制方法   总被引:2,自引:1,他引:1       下载免费PDF全文
三维表面多次波压制是海洋地震资料预处理中的重要研究课题,基于波动理论的三维表面多次波压制方法(3DSRME)是数据驱动的方法,理论上来说,可有效压制复杂构造地震数据表面多次波.但该方法因对原始地震数据采集要求高而很难在实际资料处理中广泛应用.本文基于贡献道集的概念,将稀疏反演方法引入到表面多次波压制中,应用稀疏反演代替横测线积分求和,无需对横测线进行大规模重建,进而完成三维表面多次波预测,这样可有效解决实际三维地震数据横测线方向稀疏的问题.基于纵测线多次波积分道集为抛物线的假设,为保证预测后三维表面多次波和全三维数据预测的多次波在运动学和动力学特征上基本一致,文中对预测数据实施基于稳相原理的相位校正.理论模型和实际数据的测试结果表明,本文基于稀疏反演三维表面多次波压制方法可在横测线稀疏的情况下,有效压制三维复杂介质地震资料中的表面多次波,从而更好地提高海洋地震资料的信噪比,为高分辨率地震成像提供可靠的预处理数据保障.  相似文献   

9.
声波反射成像测井是利用阵列声波测井的反射波进行井旁构造成像的测井新方法,由于是从滑行波能量较强的阵列声波信号中来获得反射波,因此反射波的提取问题尤为重要,文中采用能使波场更加聚焦的高分辨率Radon变换来提取反射波.分析得出滑行波和反射波到时在阵列数据道情况下表现为线性或近似线性特征,因而采用线性Radon变换方法.首先分析了声反射成像测井不同数据道集的波场几何特征,得出滑行波和反射波在不同道集都表现有不同程度的视速度差异,为反射波提取提供了依据.传统的Radon变换采用的最小二乘法不能使同相轴在Radon域足够聚焦,结合最大熵原理和贝叶斯原理实现了线性高分辨率Radon变换,模型数据处理结果表明后者可以在较复杂条件下提取出期望的波场.采用高阶交错网格有限差分模拟得到软、硬地层声反射成像测井波场,处理结果表明,相对于最小二乘Radon变换,采用高分辨率Radon变换后可以使波场更加聚焦,能提取出高质量的反射波,并能适应不同性质地层,相比于目前应用较多的反射波提取方法,其结果更接近于理论波形.实际数据处理也表明了该方法的有效性和准确性.该方法在声反射成像测井的速度分析、数据道插值以及偶极反射波提取方面有一定潜力,需要进一步研究.  相似文献   

10.
带有多道相关的抛物线Radon变换法分离P-P、P-SV波   总被引:3,自引:0,他引:3       下载免费PDF全文
因为多波多分量地震勘探中P-P波和P-SV波通常混杂在一起,所以较好地分离P-P波和P-SV波能够提高数据处理和解释的质量.抛物线Radon变换法在分离P-P波和P-SV波时取得了一定效果,但是在离散叠加的计算过程中会带来假频,这些假频会干扰波场分离.本文针对这一问题,将多道相关算法引入抛物线Radon变换,发展了带有多道相关的抛物线Radon变换法.该方法利用叠加信号具有相似性的特点,依据多道相关中衡量多道信号相似性的能量比标准,对叠加过程加以控制,压制变换中出现的假频.本文用该方法对合成地震记录进行了波场分离,取得了较好的效果.  相似文献   

11.
变阻尼共轭梯度算法及其性能分析   总被引:2,自引:1,他引:1       下载免费PDF全文
为了提高反演的分辨率和计算效率,本文在传统阻尼共轭梯度法的基础上,提出了变阻尼共轭梯度算法.由于在最小二乘线性反演过程中,多数情况下都要计算偏导数矩阵,而偏导数矩阵列向量的长度大小决定了解向量在对应分量方向上前进的快慢,变阻尼共轭梯度算法的提出正是利用了偏导数矩阵的这一特点.从理论上讲,它要优于传统的固定阻尼共轭梯度法.最后通过计算实例证明了该算法计算精度高,稳定性好,收敛速度快.因此采用变阻尼共轭梯度算法进行地球物理反演是切实可行的.  相似文献   

12.
A time-domain hyperbolic Radon transform based method for separating multicomponent seismic data into P-P and P-SV wavefields is presented. This wavefield separation method isolates P-P and P-SV wavefields in the Radon panel due to their differences in slowness, and an inverse transform of only part of the data leads to separated wavefields. A problem of hyperbolic Radon transform is that it works in the time domain entailing the inversion of large operators which is prohibitively time-consuming. By applying the conjugate gradient algorithm during the inversion of hyperbolic Radon transform, the computational cost can be kept reasonably low for practical application. Synthetic data examples prove that P-P and P-SV wavefield separation by hyperbolic Radon transform produces more accurate separated wavefields compared with separation by high-resolution parabolic Radon transform, and the feasibility of the proposed separation scheme is also verified by a real field data example.  相似文献   

13.
In many situations, the quality of seismic imaging is largely determined by a proper multiple attenuation as preprocessing step. Despite the widespread application of surface-related multiple elimination (SRME) and estimation of primaries by sparse inversion (EPSI) for the removal of multiples, there still exist some limitations in the process of prediction and subtraction (SRME) or inversion (EPSI), which make the efficiency of multiple attenuation less satisfactory. To solve these problems, a new fully data-driven method called closed-loop SRME was proposed, which combines the robustness of SRME and the multi-dimensional inversion strategy of EPSI. Due to the selection of inversion approach and constraint, primary estimation by closed-loop SRME may fall into a local optimum during the solving process, which lowers the accuracy of deep information and weakens the continuity of seismic events. To avoid these shortcomings, we first modified the solving method for closed-loop SRME to an L1 norm-based bi-convex optimization method, which stabilizes the solution. Meanwhile, in the L1 norm constraint-based optimization process, the 3D sparsifying transform, being a 2D Curvelet-1D wavelet transform, is brought in as a 3D sparse constraint. In the 3D sparsifying domain, the data become sparser, thus making the result of optimization more accurate, the information of seismic events more continuous and the resolution higher. Examples on both synthetic and field data demonstrate that the method proposed in this paper, compared with the traditional SRME and closed-loop SRME, have an excellent effect on primary estimation and suppress multiples effectively.  相似文献   

14.
本文提出了一种基于模型空间压缩技术的大地电磁三维反演方法.该方法在传统大地电磁三维反演理论的基础上,通过小波变换将待反演的空间域模型参数映射到小波域进行反演,获得小波域更新模型后再通过小波逆变换得到空间域反演模型.由于小波变换具有压缩特性和多尺度分辨能力,本文反演方法可在一定程度上提高反演分辨率.为了提高反演效率,我们针对基于L_1范数的模型约束求解不易收敛的反演问题,提出了一种基于模型粗糙度的简单有效的预条件处理技术.为验证本文算法的有效性,本文首先对经典的"棋盘"模型进行三维反演测试.反演结果表明本文算法的反演效率与传统方法相当,但对于深部异常体具有更好的分辨能力.最后,我们通过对实测数据反演进一步验证本文算法的有效性.  相似文献   

15.
线性同相轴波场分离的高分辨率τ-p变换法   总被引:7,自引:5,他引:7       下载免费PDF全文
基于最小二乘τ-p变换和τ-p域模型稀疏分布的假设,本文给出高分辨率τ-p变换的推导及其模型空间域的离散采样公式,同时给出了保振幅线性同相轴波场分离的算法流程.在求解本文给出的高分辨率τ-p正变换时,由于待求解的矩阵不具备最小平方法所具有的Toeplitz结构,故采用Cholesky分解法进行计算.本文模拟了井间地震和阵列声波测井中的Stoneley上下行波的分离算法过程,高分辨率正反τ-p变换且滤波所得结果显示本文算法误差小和保振幅的特点.对于在τ-p域距离很近或时间域同相轴近于水平的线性波场,高分辨率算法的聚焦作用使得所分离波场畸变小,体现本文算法精度高的优点.理论模型试算表明本文给出的高分辨率τ-p变换线性波场分离算法具有稳定性、精度高和保振幅的特点.  相似文献   

16.
Apex shift hyperbolic Radon transform (ASHRT) is an extension of hyperbolic Radon transform (HRT). We have developed a novel sparsity-promoting framework for ASHRT by employing curvelet transform (CT) in the sparse inversion. RT-based seismic data processing can be considered as an optimization problem and a mixed norms inversion, therefore, objective function with CT can promote the sparsity of the transformed domain, which makes the sparse inversion more efficient. Compared with the conventional sparse inversion of ASHRT, the proposed method weights the sparse penalization, which indicates a sparser solution of ASHRT. We use synthetic and field data examples to demonstrate the performance of ASHRT. Compared to the conventional solution, the ours may lead to more accurately reconstructed results and have a better noise immunity.  相似文献   

17.
We present preconditioned non‐linear conjugate gradient algorithms as alternatives to the Gauss‐Newton method for frequency domain full‐waveform seismic inversion. We designed two preconditioning operators. For the first preconditioner, we introduce the inverse of an approximate sparse Hessian matrix. The approximate Hessian matrix, which is highly sparse, is constructed by judiciously truncating the Gauss‐Newton Hessian matrix based on examining the auto‐correlation and cross‐correlation of the Jacobian matrix. As the second preconditioner, we employ the approximation of the inverse of the Gauss‐Newton Hessian matrix. This preconditioner is constructed by terminating the iteration process of the conjugate gradient least‐squares method, which is used for inverting the Hessian matrix before it converges. In our preconditioned non‐linear conjugate gradient algorithms, the step‐length along the search direction, which is a crucial factor for the convergence, is carefully chosen to maximize the reduction of the cost function after each iteration. The numerical simulation results show that by including a very limited number of non‐zero elements in the approximate Hessian, the first preconditioned non‐linear conjugate gradient algorithm is able to yield comparable inversion results to the Gauss‐Newton method while maintaining the efficiency of the un‐preconditioned non‐linear conjugate gradient method. The only extra cost is the computation of the inverse of the approximate sparse Hessian matrix, which is less expensive than the computation of a forward simulation of one source at one frequency of operation. The second preconditioned non‐linear conjugate gradient algorithm also significantly saves the computational expense in comparison with the Gauss‐Newton method while maintaining the Gauss‐Newton reconstruction quality. However, this second preconditioned non‐linear conjugate gradient algorithm is more expensive than the first one.  相似文献   

18.
对于被动源地震数据,运用常规的互相关算法得到的虚拟炮记录中,不仅含有一次波反射信息,还包括了表面相关多次波.然而,通过传统的被动源数据稀疏反演一次波估计(EPSI)方法,可以求得只含有一次波,不含表面相关多次波的虚拟炮记录.本文改进了传统的被动源数据稀疏反演一次波估计问题的求解方法,将被动源稀疏反演一次波估计求解问题转化为双凸L1范数约束的最优化求解问题,避免了在传统的稀疏反演一次波估计过程中用时窗防止反演陷入局部最优化的情况.在L1范数约束最优化的求解过程中,又结合了2DCurvelet变换和小波变换,在2DCurvelet-wavelet域中,数据变得更加稀疏,从而使求得的结果更加准确,成像质量得到了改善.通过简单模型和复杂模型,验证了本文提出方法的有效性.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号