n维广义Radon变换

本文提出了一种在任意变换曲线族情况下的广义Ｒａｄｏｎ变换及反投影公式，研究了在何种变换曲线情况下，其反变换能唯一恢复原出数。事实上，并非任意变换曲线簇作广义Ｒａｄｏｎ变换都能由象函数唯一恢复原出数，同时探讨了广义Ｒａｄｏｎ变换比常规Ｒａｄｏｎ变换地震数据处理中具有更广泛的应用。     

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

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

分区连续函数的Radon变换的高精度反演法

Ｒａｄｏｎ变换及其反演是ＣＴ技术的数学基础，常用的反演算法虽然在函数比较光滑的情况下都能获得较满意的结果，但当函数有较强的奇性时精度就比较差．对函数在每一子区域中是光滑，但在不同于区域交界处有间断的情况，本文给出了有高精度的反演算法．     

平面Radon变换的反演公式

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

一类广义Radon变换及反投影

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

一维波动方程的奇性反演与小波

利用线性化的技巧及Ｇ．Ｂｅｙｌｋｉｎ引进的奇性反演的概念，使波动方程非均匀背景场的反演问题有了实质性的进展，但在进行线性化简时，往往会将数值小但奇性高的项略去，因而使反演结果失真，本文利用小波变换这一工具，在化简时保留了奇性的主要部分，从而使反演所得的结果从奇性分析的观点看来更为。     

人工源极低频电磁法三维LBFGS反演（英文）

人工源极低频电磁法因其具有场源长度大,观测范围广的特点,使得极低频电磁波的传播会受到电离层和位移电流的影响。现阶段针对该方法的三维正反演研究尚处于起步阶段。作为一项探索性的尝试,本文首先给出了电离层、空气以及地下介质耦合情况下一维极低频背景电磁场的计算方案,并对电离层影响下背景电磁波的传播特征进行了分析、总结。通过将之与求解二次电场的交错网格有限差分数值模拟算法整合,实现了人工源极低频电磁法的三维正演。针对人工源极低频探测中可能遇到的近区、过渡区数据反演问题,本文进而采用针对该方法的三维有限内存Broyden-Fletcher-GoldfarbShanno(Limited-memory BFGS,LBFGS)带源反演算法,实现了对全区张量阻抗数据的直接反演。文中详细介绍了目标函数梯度计算这一LBFGS反演中的核心问题。合成数据反演算例结果表明在LBFGS反演中,选择恰当的近似Hessian矩阵能够有效提高反演效率。高低阻异常同时存在下的反演模型响应告诉我们张量阻抗反对角元素对恢复地下电性结构的贡献远大于主对角元素。与常规标量数据反演相比,张量数据反演在异常体的恢复和背景电阻率的控制方面具有明显的优势。     

n—D Radon变换与波动方程偏移

本文把Ｒａｄｏｎ变换公式推广到任意ｎ维的情况。同时结合ｎ维Ｒａｄｏｎ变换和摄动理论提出了一种既能用于地面资料又能用于ＶＳＰ资料的偏移方法。     

基于非均匀Fourier变换的地震数据重建方法研究

不规则采样地震数据会对地震数据的多道处理造成严重影响,将非均匀Fourier变换和贝叶斯参数反演方法相结合,对不规则空间带限地震数据进行反演重建.对每一个频率依据最小视速度确定出重建数据的带宽,然后从不规则地震数据中估计出重建数据的空间Fourier系数.将不规则地震数据重建视为信息重建的地球物理反演问题,运用贝叶斯参数反演理论来估计Fourier系数.在反演求解时,使用共轭梯度算法,以保证求解的稳定性,加快解的收敛速度.理论模型和实际资料处理验证了本方法的有效性和实用性.     

CDSN地震事件数据格式的自动变换

在取得ＣＤＳＮ数据之后，所遇到的第一个问题就是数据的格式变换，本文介绍了ＣＤＳＮ数据格式变换的全自动方法，采用Ｃ＿ＦＯＲＭ命令无需提供输入文件名和输出文件名就可以进行格式变换。     

Radon变换的抽样取值与小波

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

高分辨率Radon变换方法及其在地震信号处理中的应用

Radon变换方法在地震资料处理中广泛采用,在地震同相轴识别和估计方面具有良好效果．无论是倾斜叠加,还是广义Radon变换方法,一般采用最小二乘反演方法实现．目前,在提高反演算法的效率和分辨率方面仍值得研究．本文从倾斜叠加的定义出发,阐明Radon变换分辨率问题的来源和解决办法．采用最小二乘反演方法研究高分辨率抛物线Radon变换和双曲Radon变换时,给出稀疏约束预条件共轭梯度法求解的高分辨率Radon变换的实现方法,同阻尼最小二乘方法相比,分辨率和精度明显提高,文中给出了模型算例．根据有效波和多次波NMO后剩余时差不同,采用高分辨率抛物线和双曲Radon变换可以压制多次波,分别给出了方法原理,最后给出应用实例．研究表明,稀疏约束预条件共轭梯度法可以有效实现高分辨率Radon变换;数值算例表明,算法计算效率和精度较高,可以更好地实现多次波压制．     

一类广义Radon变换的反演:迭代重建再投影的外插算法

本文对不完全数据下的一类广义Radon变换的反演问题,建议用迭代重建再投影的外插算法重建图像。数值模拟结果表明,迭代重建再投影的外插算法对重建图像是有效的。     

各向异性Radon变换及其在多次波压制中的应用

即使采用分辨率很高的双曲Radon变换,对速度各向异性发育介质及长偏移距情况下的地震数据,其Radon域内能量仍不收敛.为了克服此难题,我们在Radon变换的积分路径中考虑了非双曲走时的影响,通过引入非双曲时差公式中的各向异性非椭圆率η参数,可以准确描述出长偏移距条件下来自同一层位的时距曲线,并推导了由偏移距、慢度、非椭圆率三参数控制的积分曲线正反变换公式,我们称之为各向异性Radon变换.离散化求解时,各向异性Radon变换是时变的,频率域快速算法已不适用,本文采用了最优相似系数加权Gauss-Seidel迭代算法,保持其计算精度的同时也有较高的计算效率.将此方法应用在模型数据以及实际长偏移距海上地震数据的多次波压制处理中,收到了较好的处理效果.     

Application of the continuous wavelet transform of gravity and magnetic data to estimate sub‐basalt sediment thickness

The Continuous Wavelet Transform was recently proposed for the interpretation of gravity and magnetic potential data. We utilize the Continuous Wavelet Transform of gravity and magnetic data to address one of the most common issues in exploration geophysics: mapping of sub‐basaltic sedimentary strata. We observe that the magnetic response of the basaltic layer is dominant in a three‐layer case of a basalt‐sediment‐basement, whereas the gravity signal is dominated by the base of the sediment. Thus the Continuous Wavelet Transform of the magnetic data is related to the thickness of the basalt and the Continuous Wavelet Transform of the gravity data is related mostly to the bottom of the sediment. These observations are demonstrated with a synthetic model and a few field examples. Derived depths using Continuous Wavelet Transform are in good agreement with known vertical cross‐sections. Therefore, Continuous Wavelet Transform analysis of both gravity and magnetic data offers a possibility for primary information of sub‐basaltic sediment thickness, which can provide a basis for further detailed modelling.     

Time-domain hyperbolic Radon transform for separation of P-P and P-SV wavefields

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.     

加权抛物Radon变换叠前地震数据重建

基于部分动校正（NMO）后反射同相轴在CMP道集上的抛物线走时近似,给出了加权抛物Radon变换叠前地震数据重建方法（WPRT）. WPRT通过在迭代过程中引入变化着的权系数,拓展和改进了传统抛物Radon变换方法,使其可同时完成不规则采样的规则化和空道及近偏移距道重建,且有更高的计算效率. 文中给出了应用WPRT进行近偏移距和中偏移距的空地震道重建及数据规则化的算法实现. 理论模型和实际地震资料的地震数据重建结果显示了本文算法的优点.     

基于改进线性Radon变换的井孔地震上下行波场分离方法

反射波场分离是井孔地震资料处理中极其重要的一个环节,波场分离的质量直接影响成像结果的精度.不管是VSP还是井间地震资料,其反射波时距曲线都近似直线型,根据这一特征,本文提出一种改进的线性Radon变换方法来进行井孔资料的反射波上下行波场分离.该方法基于频率域线性Radon变换,通过引入一个新的变量λ来消除变换算子对频率的依赖性,避免了求取每一频率分量对应的不同变换算子,显著降低了计算成本;文中在求解该方法对应的最小二乘问题时,引入了发展较为成熟的高分辨率Radon变换技术来进一步提高波场分离的精度.采用本文方法进行井孔地震资料的上下行波场分离可以在保证分离精度的前提下有效地提高计算效率.根据上下行波在λ-f域内分布的特殊性,设计简单的滤波算子就可实现上下行波场的分离.最后通过合成数据试算以及实际资料处理（VSP数据和井间地震数据）验证了该方法的可行性和有效性.     

Fast hyperbolic deconvolutive Radon transform using generalized Fourier slice theorem

The hyperbolic Radon transform has a long history of applications in seismic data processing because of its ability to focus/sparsify the data in the transform domain. Recently, deconvolutive Radon transform has also been proposed with an improved time resolution which provides improved processing results. The basis functions of the (deconvolutive) Radon transform, however, are time-variant, making the classical Fourier based algorithms ineffective to carry out the required computations. A direct implementation of the associated summations in the time–space domain is also computationally expensive, thus limiting the application of the transform on large data sets. In this paper, we present a new method for fast computation of the hyperbolic (deconvolutive) Radon transform. The method is based on the recently proposed generalized Fourier slice theorem which establishes an analytic expression between the Fourier transforms associated with the data and Radon plane. This allows very fast computations of the forward and inverse transforms simply using fast Fourier transform and interpolation procedures. These canonical transforms are used within an efficient iterative method for sparse solution of (deconvolutive) Radon transform. Numerical examples from synthetic and field seismic data confirm high performance of the proposed fast algorithm for filling in the large gaps in seismic data, separating primaries from multiple reflections, and performing high-quality stretch-free stacking.