快速3D抛物Radon变换地震数据保幅重建

为解决3D AVO地震数据快速保幅重建问题,在传统3D抛物Radon变换的基础上提出一种3D快速高阶抛物Radon变换方法.该方法将传统抛物Radon变换与正交多项式相结合,通过正交多项式系数描述地震数据AVO信息,确保重建后的地震数据具有良好保幅效果.同时,该方法引入新变量λ_x=q_xf和λ_y=q_yf,通过对q_xf和q_yf的整体采样,消除了3D高阶抛物Radon变换算子对频率的依赖,使变换算子的求逆过程仅需计算一次,大大节省计算时间.理论模型和实际地震资料的处理结果表明,该方法重建效率高,保幅效果良好.

     

3D高阶抛物Radon变换地震数据保幅重建

本文结合传统3D抛物Radon变换（PRT）和AVO数据正交多项式拟合,给出了3D高阶抛物Radon变换方法（HOPRT）.该变换增加了描述AVO数据变化的梯度信息和曲率信息,拓展了传统3D抛物Radon变换方法,使其在具有AVO特征的数据重建中具有更高的准确度,从而提高AVO分析的可靠性.文中给出了3D高阶抛物Radon变换进行地震数据保幅重建的流程.理论模型和实际地震资料的重建结果显示了本文方法的优点.     

多项式Radon变换

Radon变换是数据处理中广泛应用的一种方法技术.本文介绍了次数为“2”的多项式Radon变换.讨论了多项式正反Radon变换的公式、实现方法和有关计算参数的选择.通过理论模型试算,对多项式、线性、抛物线三种Radon变换进行了比较.对实际资料进行了多项式Radon变换处理,给出了消除地震反射记录中线性干扰的算例.     

多项式Radon变换

Radon变换是数据处理中广泛应用的一种方法技术.本文介绍了次数为“2”的多项式Radon变换.讨论了多项式正反Radon变换的公式、实现方法和有关计算参数的选择.通过理论模型试算，对多项式、线性、抛物线三种Radon变换进行了比较.对实际资料进行了多项式Radon变换处理，给出了消除地震反射记录中线性干扰的算例.     

3D高阶抛物Radon变换在不规则地震数据保幅重建中的应用

3D地震数据不规则采样缺失重建是地震勘探数据处理流程中的重要问题.本文提出了一种基于具有保幅特性的非均匀高阶抛物Radon变换（NHOPRT）地震数据重建方法.在最小二乘反演方程中引入Delaunay三角网格剖分来计算空间不规则加权系数,从而获得最接近完整规则数据的高阶抛物Radon变换域系数.在用SVD求解反演方程过程中,利用高阶抛物Radon变换算子在频率域为指数函数,具有线性可分解特性,将二维空间的高阶抛物Radon变换算子分解为两个独立的一维空间变换算子,减小了变换算子的矩阵大小,从而很大程度地提高了计算效率.理论模型和实际地震数据重建测试证明了本文方法的有效性以及实用性.

     

3D高阶抛物Radon变换在不规则地震数据保幅重建中的应用

3D地震数据不规则采样缺失重建是地震勘探数据处理流程中的重要问题.本文提出了一种基于具有保幅特性的非均匀高阶抛物Radon变换（NHOPRT）地震数据重建方法.在最小二乘反演方程中引入Delaunay三角网格剖分来计算空间不规则加权系数,从而获得最接近完整规则数据的高阶抛物Radon变换域系数.在用SVD求解反演方程过程中,利用高阶抛物Radon变换算子在频率域为指数函数,具有线性可分解特性,将二维空间的高阶抛物Radon变换算子分解为两个独立的一维空间变换算子,减小了变换算子的矩阵大小,从而很大程度地提高了计算效率.理论模型和实际地震数据重建测试证明了本文方法的有效性以及实用性.     

数据严重缺少的Radon变换的反演

在许多实际问题中,由于客观条件的限制经常会遇到数据不全的Ｒａｄｏｎ变换的反演问题．在地球物理学的一些领域中更会遇到数据严重缺少的Ｒａｄｏｎ变换的反演问题──有限个点源及有限个方向投影问题．本文讨论此类反演问题,并给出具体的数值计算方法．     

数据严重缺少的Radon变换的反演

在许多实际问题中，由于客观条件的限制经常会遇到数据不全的Ｒａｄｏｎ变换的反演问题．在地球物理学的一些领域中更会遇到数据严重缺少的Ｒａｄｏｎ变换的反演问题──有限个点源及有限个方向投影问题．本文讨论此类反演问题，并给出具体的数值计算方法．     

λ-f域加权抛物Radon变换地震数据重建方法研究

大多数地震处理技术都要求地震数据具备完整性和规则性,然而,由于诸多因素的影响,地震勘探所采集到的资料普遍存在数据缺失,需要对地震数据进行重建.本文在传统Radon变换重建的基础上提出一种λ-f域加权抛物Radon变换的地震数据重建方法.通过引入新变量λ,消除了Radon变换算子对频率的依赖,使得Radon变换算子及算子的逆只需计算一次,显著提高了计算效率.同时,在λ-f域抛物Radon变换迭代计算过程中引入变化的权系数,更好地实现了λ-f域的能量聚焦.理论模型及实际数据试算表明,文中方法对地震数据重建精度较高,单道对比吻合较好.     

Radon变换的抽样取值与小波

Radon变换在实际应用中必须考虑抽样取值的问题,通信理论中常用的Shannon取样定理在地球物理学等领域中并不适用,因为它要求函数是频谱有限的.本义介绍利用小波理论建立新的取样定理的方法,这一方法的优点在于可以根据具体问题的要求选取适当的基本空间V_0,从而获得较佳的抽样插值效果.     

一类曲线图象重建卷积反投影算法

本文给出了一类平面曲线 s=4(x,y)·coxθ+B(x,y)·sinθ上的 Radon 变换卷积反投影算法,并进行了数值模拟,其重建结果与真实数据相吻合,其中 A(x,y)和 B(x,y)满足一定条件.     

Rayleigh-Wave Dispersive Energy Imaging Using a High-Resolution Linear Radon Transform

Multichannel Analysis of Surface Waves (MASW) analysis is an efficient tool to obtain the vertical shear-wave profile. One of the key steps in the MASW method is to generate an image of dispersive energy in the frequency-velocity domain, so dispersion curves can be determined by picking peaks of dispersion energy. In this paper, we propose to image Rayleigh-wave dispersive energy by high-resolution linear Radon transform (LRT). The shot gather is first transformed along the time direction to the frequency domain and then the Rayleigh-wave dispersive energy can be imaged by high-resolution LRT using a weighted preconditioned conjugate gradient algorithm. Synthetic data with a set of linear events are presented to show the process of generating dispersive energy. Results of synthetic and real-world examples demonstrate that, compared with the slant stacking algorithm, high-resolution LRT can improve the resolution of images of dispersion energy by more than 50%.     

一类广义Radon变换及反投影

本文提出了一类广义Ｒａｄｏｎ变换及反投影公式。这种变换形式比常规Ｒａｄｏｎ变换形式具有更广泛的意义及应用。数值计算结果表明，反投影恢复性较好。     

局域双曲线Radon变换压制散射波研究与应用

一般而言,由于地下非均质体的存在所产生的二次波源,由它再生成新的波场,叫散射波场。目前,将散射波作为有效波来成像,已开始在溶洞和裂缝等特殊地质体的识别中得到应用。但对于野外资料采集来说,地表复杂地区,如戈壁、砾石区和山前带,大量存在的散射波却是干扰波,它们的存在会严重影响资料的品质,而其研究与实际应用国内外还很少。因此,通过正演模拟,分析散射波的基本特征,在此基础上研究散射波的去噪方法显得十分必要。本文从地震波运动学时距关系出发,研究了反射波和散射波的几何特征;然后用有限差分正演,模拟了散射波场,用理论模型研究并测试了局域双曲线Radon变换散射波去噪新方法。对于实际炮集资料,分析了F-K滤波方法压制散射噪声的局限,采用局域双曲Radon变换有效地去除了炮集中存在的散射噪声,取得了较好的应用效果。     

平面Radon变换的反演公式

本文借助球面平均法和Ｈｉｌｂｅｒｔ空间算子理论,给出了平面上Ｒａｄｏｎ变换的反演公式,此公式是具体的、构造性的,便于数值计算,或者进一步探讨重建图象的性质。     

广义Radon变换的正交展开和反演

本文具体给出了广义Ｒａｄｏｎ变换在函数空间Ｌ＾２（Ｒ＾ｎ,Ｗｎ）上的正交函数展开形式,由此得到了广义Ｒａｄｏｎ变换的反演公式,从而推广了Ｌｏｕｉｓ和Ｄａｖｉｓｏｎ的工作。     

沿圆曲线的Radon变换数值解

研究具有紧支集且在支集内连续的二元函数沿上半圆曲线的Radon变换反演问题。基于对投影函数的Fourier变换,反演问题可以归结为具有弱奇性及震荡核的Abel积分方程的求解。我们证明了当圆曲线中心及半径在一定范围内变化时,在已知沿上半圆曲线的Radon变换情况下,这个积分方程的解具有唯一性,并给出了消除Abel积分方程弱奇性的数值方法。在考虑投影数据噪声的情况下,给出了多次加权改善系数矩阵条件数稳定的数值方法,并通过数值模拟验证所提出方法的有效性。     

求解Radon变换改进Fourer算法的误差分析

本文将对求解Ｒａｄｏｎ变换的改进Ｆｏｕｒｉｅｒ算法进行误差分析,证明了在Ｌ２范数下改进Ｆｏｕｒｉｅｒ算法是收敛的且具有Ｏ（１／ｑ）敛阶,其中２ｑ为像素点矩陈阶数。