弹性波逆散射微扰论

本文将普遍声逆散射微扰论应用于弹性波层析成像问题,在Born变换下推出了以旋转角为补偿参数的各阶微扰重建公式,实现了对非均匀各向同性散射体内3个参数（质量密度ρ和两个Lamé系数λ,μ）的同时重建. 对于层析成像问题,在弹性波的传播过程中P波与SV波有耦合,但它们不会和SH波发生耦合,于是可以得到3个形式相对简单的标量方程. 在Born变换下,在散射波中引入微扰参数,将散射体的3个参数分别按该微扰参数展开,然后利用二维自由空间的Green函数分别得到散射的P波、SV波和SH波的积分表示. 最后,经一维傅氏变换后,得到Born变换下散射体3个参数的各阶微扰重建公式.     

角度域广义Radon变换多参数非迭代线性反演

基于Born近似的线性反演是当前流行的定量化地震成像与反演方法, 然而, 现有线性反演方法存在多种不可避免的问题, 比如迭代计算成本过高, 背景模型依赖, 以及多参数提取不确定/不稳定.为避免这些问题, 本文提出了角度域广义Radon变换(AD-GRT)多参数非迭代反演方法(2D声波双参数介质).AD-GRT是一种地震数据与角度域模型之间的积分变换, 该角度域模型能够保留地震散射数据的完整信息.基于该角度域模型, 地震数据可被精确重建, 即使背景模型是错误的.介质双扰动参数可在角度域模型中准确且稳定提取.为解决传统角度域离散方法存在的幅值震荡和不连续问题, 设计了两种有效的角度域离散单元分裂方法, 实现了平滑连续的角度域幅值提取.为提高AD-GRT反演的分辨率, 考虑了震源子波反滤波.利用三个数值模型验证了AD-GRT在角度域模型反演、数据重建、多参数提取方面的有效性.

     

A new scheme for numerical modelling of flow and transport processes in 3D fractured porous media

The objective of this work is to develop a new numerical approach for the three-dimensional modelling of flow and transient solute transport in fractured porous media which would provide an accurate and efficient treatment of 3D complex geometries and inhomogeneities. For this reason, and in order to eliminate as much as possible the number of degrees of freedom, the fracture network, fractures and their intersections, are solved with a coupled 2D–1D model while the porous matrix is solved independently with a 3D model. The interaction between both models is accounted for by a coupling iterative technique. In this way it is possible to improve efficiency and reduce CPU usage by avoiding 3D mesh refinements of the fractures. The approach is based on the discrete-fracture model in which the exact geometry and location of each fracture in the network must be provided as an input. The formulation is based on a multidimensional coupling of the boundary element method-multidomain (BEM-MD) scheme for the flow and boundary element dual reciprocity method-multidomain (BE-DRM-MD) scheme for the transport. Accurate results and high efficiency have been obtained and are reported in this paper.     

基于小波包变换的SAR干涉图去噪研究

基于SAR复干涉图中噪声满足加性噪声模型,提出对复干涉图实部和虚部分别应用小波包软阈值方法的去噪方案。以Daubeachies小波为小波包基函数,对添加了噪声的模拟SAR干涉图进行去噪实验,发现三级小波包分解的去噪效果明显优于一级和二级分解。并进一步对实际的含噪干涉图进行了去噪实验,与小波包软阈值方法直接对干涉图去噪的结果进行了比较。实验结果证实了本文方案的有效性。     

Application of wavelet transforms to the decoupling of free-roll decay data

This study presents the results of the decoupling of a free-roll decay test by wavelet transforms. A free-roll decay test was performed to determine the coefficients of damping terms in equations of motion. No motion of the model was restricted during the experiment. A slight yaw motion was found while the model was in the free-roll decay motion. It was necessary to extract a pure roll motion from the experimental data. Wavelet transforms were applied to the signals to extract the pure roll motion. The results were compared to those found through the Fourier transform. Discrete wavelet transforms were able to efficiently decouple the test signals, while the continuous wavelet transform and the Fourier transform could not.     

基于小波变换法定义的波群参数

通过小波变换分析的波能过程定义了两个描述波浪群性的参数,由数值模拟波浪和实测波浪资料对其与常用的波群参数进行对比分析,结果表明基于小波波能过程定义的群性参数是有效的,从而展示了小波变换用于在时频域上分析波群的能力。详细探讨了波浪记录长度对群性参数稳定性的影响,分析结果表明,波浪观测长度对于波群参数的影响较大,在考虑波浪群性的波浪模拟及分析时,建议模拟时间长度应在400~500个波以上。     

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

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

基于f-x域时频非凸正则化低秩矩阵近似的共偏移距道集去噪方法

随机噪声压制是地震数据处理的关键环节,而时频稀疏低秩近似算法逐道处理地震数据过程中无法利用信号的道间相干性。为此,将时频稀疏低秩近似与f-x域去噪结合,提出一种f-x域时频非凸正则化低秩矩阵近似算法。该算法对f-x域中每一单频分量作时频分解后,再对时频系数矩阵作低秩矩阵近似计算,能够利用信号和噪声的时频谱差异实现非平稳信号去噪处理。与共炮点道集和共中心点道集相比,共偏移距道集具有平缓甚至接近水平的同相轴结构,基本满足f-x域去噪的线性同相轴假设前提,建议将所提算法应用于共偏移距道集去噪处理。通过数值模拟和实际地震数据试算,证明本文方法能够有效压制随机噪声,同时保持有效信号不被损害。     

Applications of Lodrigues Matrix in 3D Coordinate Transformation

Three transformation models (Bursa-Wolf, Molodensky, and WTUSM) are generally used between two data systems transformation. The linear models are used when the rotation angles are small; however, when the rotation angles get bigger, model errors will be produced. In this paper, we present a method with three main terms: ① the traditional rotation angles θ , φ ,ψ are substituted with a , b, c which are three re-spective values in the anti-symmetrical or Lodrigues matrix; ② directly and accurately calculating the formula of seven parameters in any value of rotation angles; and ③ a corresponding adjustment model is established. This method does not use the triangle function. Instead it uses addition, subtraction, multiplication and division, and the complexity of the equation is reduced, making the calculation easy and quick.     

基于探地雷达信号处理的小波基选取研究

针对探地雷达信号处理和分析时小波基选取存在的问题,本文在分析探地雷达信号特点的基础上,首先从理论上讨论小波基的选取准则,然后再从实验角度进行对比、判别,认为在进行小波分解和重构时应该分别选择不同的小波基函数进行处理,这样可以保证重构信号的精确度,增强对信号的处理能力,从而也突破了以往分解与重构时都采用同一个小波基进行处理的做法.最后通过实际资料的处理,指出bior2.6小波基在进行雷达信号处理时效果最佳,不仅去噪彻底,而且能够保留有效信号的高频部分,提高信号的分辨率和信噪比,为后续解释工作打好了基础.