探地雷达小波域三维波动方程偏移

阐述了矩阵多分辨分析理论中的标准形式与非标准形式,并以Hilbert算子为例,说明了算子多分辨表示的压缩效果,为小波域偏移算法奠定了理论基础.从三维雷达波动方程出发,利用爆炸反射原理和浮动坐标变换,推导出三维探地雷达波动方程差分格式,并通过方程分裂算法及小波多分辨算法,在小波域求解波场外推矩阵,进而得到探地雷达小波域三维波动方程偏移算法,在此基础上,开发了探地雷达小波域偏移处理程序,并把该程序应用于三个球体空洞的3-D正演结果及实际的雷达数据中,通过对比偏移处理前后的雷达资料,得知该三维偏移算法能使3-D正演剖面中的反射波归位、绕射波收敛,极大地提高了雷达剖面的分辨率,有利于探地雷达资料的地质解释.     

基于小波熵的地下管线探地雷达频率波数偏移成像研究

地下管线实际形状与准确位置的确定是探地雷达图谱解释的主要内容之一.制约管线探地雷达反射归位精度最主要的因素是电磁波在地层中的传播速度,如何准确确定电磁波传播速度是探地雷达偏移处理能否成功的关键.为便于对探地雷达图谱特征进行定量解释,开展了不同大小、埋深与重叠等情况管线的正演模拟与模型实验.结合小波熵理论与频率波数偏移成像方法,提出了一种基于小波熵的探地雷达频率波数偏移成像方法,将其应用于地下管线正演模拟与模型实验探地雷达信号偏移处理.研究结果表明:探地雷达图像经基于小波熵的频率波数偏移处理后,电磁波速度的估计值与真实值误差可控制在2％以内;随着管线管径大小与埋置深度的增加,绕射双曲线的开口弧度不断增大,管与管之间绕射波交叉干扰越严重;探地雷达信号经基于小波熵的频率波数偏移处理后可大大提高图像的分辨率,压制了多次波的干扰,使绕射双曲线更收敛,能量更聚焦,从而有利于地下管线实际形状与具体位置的准确提取.研究结果对城市地下管线探地雷达图谱解译工作提供技术支持.     

相控阵探地雷达数据傅氏有限差分偏移研究

偏移归位是探地雷达资料处理中的一种重要方法。本文利用傅里叶有限差分法对相控阵探地雷达样机采集的实测数据进行偏移归位。异常点绕射波形得到有效的归位,异常体空间位置与实际测量的结果基本一致并且处理速度快,不仅表明了该方法的有效性,证明相控阵探地雷达技术具有一定研究价值,也为仪器的进一步完善提供了思路。     

Kirchhoff偏移法在探地雷达正演图像处理中的应用

本文首先通过对波动方程的分析,得出了声波波动方程和雷达波波动方程形式一致性,从而说明了把广泛应用在地震数据处理中的偏移技术引入到GPR资料处理中的可行性;其后说明了时域有限差分法（FDTD）的原理,并用它合成了几种常见的雷达正演剖面;最后利用Kirchhoff积分偏移法对正演所得的雷达剖面进行偏移处理,通过对比偏移处理前后的雷达正演剖面,可知Kirchhoff积分偏移法能使雷达正演剖面中的反射波的归位,绕射波收敛,从而大大地提高了雷达正演剖面的分辨率,更好地指导GPR剖面的地质解释和验证偏移方法的有效性.     

衰减雷达波有限元偏移

高频雷达波在地球介质中有较强的衰减,反演中不可忽略．为此文中首先给出了含衰减项的雷达波的有限元方程及其偏移理论．用有限差分法或有限元法可正演合成雷达波资料,加入一定的扰动后用含衰减项的雷达波有限元方程做偏移,实例结果表明,考虑衰减项的偏移结果能使界面更好地归位,这为提高探地雷达地质解释的分辨率提供了可能性,为逐渐地实现符合雷达波自身动力学特点的处理系统奠定了基础．     

基于波动方程有限差分算法的接收函数正演与偏移

针对接收函数正演与偏移, 本文采用波动方程有限差分算法. 借鉴成熟的勘探地震学方法, 引入等效速度概念, 建立接收函数转换波与地震勘探反射波的等效走时方程, 实现了基于波动方程有限差分算法的接收函数正演与偏移. 数值计算表明, 波动方程有限差分叠后偏移方法可以对点绕射和穹隆构造模型实现高精度成像. 本文利用数值计算讨论了波动方程有限差分叠后偏移与Kirchhoff叠后偏移对于接收函数偏移的适用性, 还对偏移过程中速度模型的误差进行了分析.     

混合域单程波传播算子及其在偏移成像中的应用

以地震波的单程波传播方程为基础,利用算子近似展开的方法推导出了当前波动方程叠前深度偏移方法研究中广泛使用的裂步Fourier、Fourier有限差分和广义屏传播算子的一般形式及其近似式．讨论了它们间差异、相互关系以及他们的特点,最后给出了基于裂步Fourier、Fourier有限差分和广义屏传播算子的偏移成像方法时Marmousi模型的偏移成像结果,以说明它们间的优劣与计算效率．     

机载探地雷达数值模拟及逆时偏移成像

机载探地雷达可以用于人类无法到达的危险地区、植被严重覆盖的地下目标体探测,然而由于机载探地雷达的特殊性,影响机载探地雷达探测效果的因素包括天线的极化方向、天线的飞行高度以及地表粗糙度等.为了研究这些影响因素与探测效果之间的关系,用三维时间域有限差分模拟电磁波的传播过程,以沙漠地区地下空洞掩体的机载探地雷达探测为实例,分别模拟了不同天线极化方向、天线高度及地表粗糙度情况下的机载探地雷达剖面,分析了各因素对机载探地雷达探测地下空洞目标体的影响.天线极化方向与目标体走向垂直更有利于地下目标体探测;天线距离地表越近,可以获得更高分辨率的雷达剖面;沙漠地表起伏越大,雷达剖面中的散射杂波能量越强,浅部地下目标体信号容易被掩盖.为了消除起伏地形造成的散射杂波,提出用逆时偏移成像技术对共炮集机载探地雷达数据进行偏移成像,成像结果优于基尔霍夫偏移成像结果.     

机载探地雷达的进展以及数值模拟

机载探地雷达可能解决危险环境或广域条件下的近地表探测问题,用于解决环境、生态或军事方面的问题.然而由于种种原因,该技术的发展却显得比较慢.为了推进该技术的发展,本文介绍了目前世界范围内机载探地雷达的进展,并利用时间域有限差分法对一些典型模型进行数值模拟,并用特定的偏移成像方法对模拟结果进行成像.目前存在的机载探地雷达主要有三种类型:第一种为将常规探地雷达天线悬挂在直升飞机上,第二种为针对机载探地雷达开发的雷达系统,第三种为具有探地能力的合成孔径雷达.数值模拟结果表明,不管是水平地面的情况下,还是起伏界面的情况下,机载探地雷达都能清楚探测一定深度范围内的地下目标.可见,机载探地雷达是存满希望的一种方法.     

SEG/EAGE盐丘和推覆体模型的波动方程三维叠前深度偏移成像

三维波动方程叠前深度偏移是复杂介质中进行构造成像、弹性参数反演的重要环节．由于其技术实现不仅涉及波场延拓理论的创新,而且需要大规模计算,因而研究难度较大．本文以实验效果的取得为目的,完整地实现了SEG／EAEG盐丘和推覆体模型的三维波动方程辛几何算法的叠前深度偏移成像计算．文中详细考察了所研制的波动方程三维叠前深度偏移软件系统及其对复杂地质构造的成像能力,具体包括：1)对于盐丘模型,文中讨论了成像参数的选择、地震子波对成像精度的影响、完成二维及三维叠前深度偏移的比较;2)对推覆体模型,文中进行了脉冲响应测试;3)由两个模型的成像结果可见本文的波动方程三维叠前深度偏移软件系统已具有适应强速度横向变化、复杂构造的成像能力。     

Reverse time migration of ground-penetrating radar with full wavefield decomposition based on the Hilbert transform

The conventional reverse time migration of ground-penetrating radar data is implemented with the two-way wave equation. The cross-correlation result contains low-frequency noise and false images caused by improper wave paths. To eliminate low-frequency noise and improve the quality of the migration image, we propose to separate the left-up-going, left-down-going, right-up-going and right-down-going wavefield components in the forward- and backward-propagated wavefields based on the Hilbert transform. By applying the reverse time migration of ground-penetrating radar data with full wavefield decomposition based on the Hilbert transform, we obtain the reverse time migration images of different wavefield components and combine correct imaging conditions to generate complete migration images. The proposed method is tested on the synthetic ground-penetrating radar data of a tilt-interface model and a complex model. The migration results show that the imaging condition of different wavefield components can highlight the desired structures. We further discuss the reasons for incomplete images by reverse time migration with partial wavefields. Compared with the conventional reverse time migration methods for ground-penetrating radar data, low-frequency noise can be eliminated in images generated by the reverse time migration method with full wavefield decomposition based on the Hilbert transform.     

叠前逆时偏移影响因素分析

反射地震勘探中的偏移成像技术是获取地下介质构造形态最有效的手段之一.在叠前深度域偏移方法中,目前工业界采用的方法包括基于射线理论的波动方程积分解法和基于波动理论的微分波动方程单程波解法,这两类方法难以处理地震波横向速度变化剧烈的高陡倾角构造成像问题.近年来勘探地震学研究领域发展起来的叠前逆时偏移采用了双程波求解微分波动方程的算法,这种方法具有相位准确、不受介质横向速度变化和高陡倾角构造的影响、成像精度高、可以利用回转波正确成像等优点,从理论上弥补了当前工业界常规地震偏移所面临的成像缺陷.然而,叠前逆时偏移成像方法从理论走向实用尚需解决如下问题:计算速度和数据存储空间的节省、初始速度模型的建立、震源子波的选择、数值模型边界条件的定义和假像的消除等等.对于计算速度和存储量大的问题,随着计算机硬件的快速发展,将会不断得到改善,同时可以采取一些计算技术和存储策略来加以缓解.本文主要针对初始速度模型的建立、震源子波的选择、数值模型边界条件的定义和假像的消除这些因素,利用简单模型进行了分析.对于反射波造成的传播路径上的假像,给出了一种振幅补偿滤波方法.对勘探地球物理学界给出的SEG/EAGE二维盐丘模型、Marmousi模型和本研究设计的崎岖海底模型进行了叠前逆时偏移成像,均取得了较好的成像效果.     

双平方根波动方程偏移速度分析

传统的剩余校正（RMO）偏移速度分析方法基于走时原理,在陡倾角和欠照明地区,因为不能得到充分的角度域信息而失效.本文将展示一种基于波场延拓理论的偏移速度分析方法,即波动方程偏移速度分析（WEMVA）.这种方法先利用成像优化方法获得剩余成像,再利用剩余成像反演剩余速度.此类方法继承了波动方程偏移方法的优点和缺点.波动方程偏移速度分析是一种线性反演方法,它要求对Born近似的展开序列作一阶截断.高阶部分的丢失必然带来巨大的截断误差,因此剩余成像必须也进行线性化,以适应大速度扰动和大延拓步长.因此,在此类算法中,剩余成像的获取和线性化是偏移速度分析的关键.在叠前偏移算子中,因为双平方根算子的数学表达式更为简洁,所以本文基于对波动方程偏移速度分析初步讨论,并通过模型验证其原理.     

利用矩阵分解理论的双程波方程叠前深度偏移方法

叠前深度偏移理论及方法一直是地震数据成像中研究的热点问题.业界对单程波叠前深度偏移方法和逆时深度偏移开展了深入的研究,但对双程波方程波场深度延拓理论及成像方法的研究还鲜有报道.本文以地表记录的波场值为基础,利用单程波传播算子估计波场对深度的偏导数,为在深度域求解双程波方程提供充分的边界条件,并提出利用矩阵分解理论实现双程波方程的波场深度外推.通过对强速度变化介质中传播波场的计算,与传统的单程波偏移方法相比,本文提出的偏移方法计算的波场与常规有限差分技术计算的波场相一致,证明了本方法计算的准确性.通过对SEAM模型的成像,在相同的成像参数下,与传统的单程波偏移算法和逆时深度偏移算法方法相比,本文提出的偏移方法能够提供更少的虚假成像和更清晰的成像结果.本文所提偏移算法具有深度偏移和双程波偏移的双重特色,推动和发展了双程波叠前深度偏移的理论和实践.     

利用偏移进行视反射率估计的初步研究

视反射率估计是地震数据处理解释中的一项重要内容,通常采用反演的方法得到.本文以地震偏移和地震线性反演理论相结合为基础,并利用保幅单程波传播算子和保幅波动方程叠前偏移算法以及成像空间中的角度域波动方程偏移成像和照明补偿等方法技术,提出了一种利用单程波波动方程偏移进行地下反射面视反射率估计方法,并进行了理论模型的数值试验.这种估计方法得到的视反射率估计是一种近法向入射的小角度反射率.     

基于控制照明的合成震源记录交互剩余偏移速度分析

提出了一种新的偏移速度分析方法——基于控制照明的合成震源记录交互剩余偏移速度分析方法．与其他类似偏移速度分析方法的不同点在于：（1） 叠前深度偏移采用基于波动理论的快速合成震源记录算法；（2）偏移方法采用平面波震源，与速度分析方法一致；（3）应用控制照明技术，避免了因横向变速而导致的平面波震源波场在传播过程中的畸变，从而减小了速度分析的误差；（4）实用的速度谱设计，使交互偏移速度分析可行且易于操作．模型和新疆实际资料的试算表明该方法是一种有效和实用的偏移速度分析方法．     

等效偏移距偏移方法在声波反射成像测井中的应用研究(英文)

声波反射成像测井能够获得井眼周围构造的重要信息,然而,由于接收到的反射波信号远小于井眼模式波、信噪比较低,且每次发射只有8道接收,因此应用常规的偏移成像处理方法成像效果不好。本文应用一种基于散射理论的等效偏移距偏移方法,对声波反射成像测井模拟数据及现场数据进行偏移处珲。结果表明,与常规叠后深度偏移方法相比,等效偏移距方法可以有效提高覆盖次数,对于低信噪比的声波反射测井资料可以获得较好的井旁构造成像效果。