Broadband synthetic aperture borehole radar interferometry

Trials in mines have established that wideband VHF borehole radars (BHR), working in the 10–100-MHz band, can be used to probe the rockmass between boreholes over ranges from <5 m to as much as 150 m with submeter resolution. There is evidence that ore bodies reflect these radar signals both specularly and diffusely, much as the ground/air interface does when overflown by synthetic aperture radar (SAR). In both SAR and BHR, multiple flight lines, together with diffuse reflections admit the possibility of developing interferometric 3D images of the object. This paper examines the possibility of imaging buried objects in three dimensions by interferometrically combining broadband VHF borehole radar profiles shot in adjacent pairs of boreholes. Broadbanding in BHR has the advantage of releasing the image from 2nπ phase ambiguities, but practically, interferometric borehole radar (InBHR) needs high signal-to-noise ratios (SNR) to avoid noise capture. This means that 3D InBHR is limited to ranges in wavelengths which are less than the rock's attenuation factor Q. Interferometric methods are developed which are capable of mapping ore bodies and other structures in three dimensions. Tangent plane migration methods are developed here in order to reconstruct surfaces that lie in the near-field of sparse interferometric arrays.     

瑞雷波多阶模式频散曲线稀疏正则化反演方法研究

目前瑞雷波多阶模式频散曲线反演中仅考虑数据的拟合,缺乏对模型的约束,不能很好地刻画地层间断面的问题,针对此问题,研究了瑞雷波多阶模式频散曲线稀疏正则化反演方法.正演模拟基于广义反射-透射系数法,数值计算上采用一种快速求根方法,与二等分方法相比,能够在很短的时间内达到最优的收敛效果;反演建模时采用L1范数正则化方法对模型...     

基于稀疏表示的地面磁共振信号提取方法

在多孔隙含水层中地面磁共振（surface nuclear magnetic resonance,SNMR）信号呈现多弛豫衰减特性,常规盲源分离方法和单指数拟合方法引起信号严重失真和信息缺失等问题.本文提出了基于稀疏表示的随机噪声背景下多弛豫SNMR信号的提取方法.根据SNMR信号的衰减特征,设计了精确刻画SNMR信号且与随机噪声不相关的离散衰减余弦冗余字典.其次,针对多弛豫SNMR信号稀疏度未知的问题,通过设置合理的残差比阈值控制迭代次数,改进了广义正交匹配追踪（generalized orthogonal matching pursuit,gOMP）算法,使得该方法应用于SNMR信号的提取时,具有更好的自适应性和普适性.再次,鉴于SNMR测量数据为多次独立重复采集的结果,提出了基于数据流的SNMR信号提取策略,在提高算法鲁棒性的同时,保证了信号提取结果的唯一性.最后,通过仿真和实测数据证明了基于gOMP算法的稀疏表示方法可以显著地提升多弛豫SNMR信号的提取质量,降低随机噪声对含水层反演结果的影响,提高SNMR探测能力.     

基于广义模型约束的时间域航空电磁反演研究

由于航空电磁具有海量数据,因此快速有效的成像和反演手段至关重要.本文针对层状介质模型推导与实现了广义模型约束条件下时间域航空电磁一维反演.从正则化反演的目标函数出发,通过改变模型约束项构造Lp范数反演和聚焦反演,进而通过改变模型求解域构造出基于小波变换的稀疏约束反演.针对不同反演方法目标函数的构建方式,本文进一步从数学原理上分析不同反演方法的预期效果,并通过理论模型和实测数据进行验证.结果表明L0.8范数反演、聚焦反演和基于小波变换的稀疏约束反演可以得到更符合地下层状介质陡变界面的反演结果.     

Inversion-based deblending using migration operators

In this paper, we compare the denoising- and inversion-based deblending methods using Stolt migration operators. We use Stolt operator as a kernel to efficiently compute apex-shifted hyperbolic Radon transform. Sparsity promoting transforms, such as Radon transform, can focus seismic data into a sparse model to separate signals, remove noise or interpolate missing traces. Therefore, Radon transforms are a suitable tool for either the denoising- or the inversion-based deblending methods. The denoising-based deblending treats blending interferences as random noise by sorting the data into new gathers, such as common receiver gather. In these gathers, blending interferences exhibit random structures due to the randomization of the source firing times. Alternatively, the inversion-based deblending treats blending interferences as a signal, and the transform models this signal by incorporating the blending operator to formulate an inversion problem. We compare both methods using a robust inversion algorithm with sparse regularization. Results of synthetic and field data examples show that the inversion-based deblending can produce more accurate signal separation for highly blended data.     

基于数据规则化和稀疏反演的三维表面多次波压制方法

表面多次波是海洋地震勘探中的主要问题.目前,二维数据驱动的表面多次波压制技术(SRME)已经比较成熟,并且已经成为工业界压制海洋表面多次波的主流方法.但是由于二维SRME算法没有考虑横测线方向上多次波的贡献,导致在处理实际三维海洋资料时存在比较大的误差.将二维SRME算法扩展到三维空间后可以得到三维SRME算法,但是由于目前实际采集的三维海洋资料的观测系统存在拖缆漂移,而且横测线方向采样过于稀疏,直接应用三维SRME算法无法准确预测表面多次波.本文提出的通过数据规则化配合稀疏反演的三维表面多次波压制方法能够解决这种实际资料和三维SRME算法之间的矛盾.本文通过研究数据规则化与反规则化技术,使得数据分布满足三维SRME的要求;通过研究稀疏反演技术,有效解决了横测线方向采样稀疏对于多次波预测的影响,三维实际海洋资料的应用结果验证了方法的有效性和可行性.     

The utilization of the double focal transformation for sparse data representation and data reconstruction

In many cases, seismic measurements are coarsely sampled in at least one dimension. This leads to aliasing artefacts and therefore to problems in the subsequent processing steps. To avoid this, seismic data reconstruction can be applied in advance. The success and reliability of reconstruction methods are dependent on the assumptions they make on the data. In many cases, wavefields are assumed to (locally) have a linear space–time behaviour. However, field data are usually complex, with strongly curved events. Therefore, in this paper, we propose the double focal transformation as an efficient way for complex data reconstruction. Hereby, wavefield propagation is formulated as a transformation, where one‐way propagation operators are used as its basis functions. These wavefield operators can be based on a macro velocity model, which allows our method to use prior information in order to make the data decomposition more effective. The basic principle of the double focal transformation is to focus seismic energy along source and receiver coordinates simultaneously. The seismic data are represented by a number of localized events in the focal domain, whereas aliasing noise spreads out. By imposing a sparse solution in the focal domain, aliasing noise is suppressed, and data reconstruction beyond aliasing is achieved. To facilitate the process, only a few effective depth levels need to be included, preferably along the major boundaries in the data, from which the propagation operators can be calculated. Results on 2D and 3D synthetic data illustrate the method's virtues. Furthermore, seismic data reconstruction on a 2D field dataset with gaps and aliased source spacing demonstrates the strength of the double focal transformation, particularly for near‐offset reflections with strong curvature and for diffractions.     

L1范数约束被动源数据稀疏反演一次波估计

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

基于压缩感知的Curvelet域联合迭代地震数据重建

由于野外采集环境的限制,常常无法采集得到完整规则的野外地震数据,为了后续地震处理、解释工作的顺利进行,地震数据重建工作被广泛的研究.自压缩感知理论的提出,相继出现了基于该理论的多种迭代阈值方法,如CRSI方法（Curvelet Recovery by Sparsity-promoting Inversion method）、Bregman迭代阈值算法（the linearized Bregman method）等.CSRI方法利用地震波形在Curvelet的稀疏特性,通过一种基于最速下降的迭代算法在Curvelet变换域恢复出高信噪比地震数据,该迭代算法稳定,收敛,但其收敛速度慢.Bregman迭代阈值法与CRSI最大区别在于每次迭代时把上一次恢复结果中的阈值前所有能量都保留到本次恢复结果中,从而加快了收敛速度,但随着迭代的进行重构数据中噪声干扰越来越严重,导致最终恢复出的数据信噪比低.综合两种经典方法的优缺点,本文构造了一种新的联合迭代算法框架,在每次迭代中将CRSI和Bregman的恢复量加权并同时加回本次迭代结果中,从而加快了迭代初期的收敛速度,又避免了迭代后期噪声干扰的影响.合成数据和实际数据试算结果表明,我们提出的新方法不仅迭代快速收敛稳定,且能得到高信噪比的重建结果.     

固体介质中孤立波的传播及演化特征

由于地球介质中广泛存在断裂、微裂缝等地质现象,实际地震资料中会出现形似孤立波的非线性地震现象。因此,对固体介质中孤立波的研究有利于解释这些非线性地震现象的形成机制。本文基于KdV方程,以雷克子波作为初始条件,采用伸缩子机理构建体力模型,利用有限差分的方法模拟了孤立波的演化过程。理论结果表明,非线性地震纵波可以从雷克子波逐渐演化成孤立波,而且地震波的初始振幅和频散系数对模拟结果也有重要影响。通过与实际资料的对比也能说明这种演化的可能性。同时根据方程系数矩阵中元素带状分布的特征,采用稀疏矩阵的存取方法,可以减小计算内存,提高计算效率。