探地雷达多尺度波形反演方法研究

来源 :吉林大学 | 被引量 : 0次 | 上传用户:panzx777
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
探地雷达(GPR)因其操作方便、工作效率高、具有描绘复杂近地表结构能力等特点,而被诸多行业青睐并应用于环境,水文和工程地球物理中。常规的探地雷达成像主要依赖于对雷达信号进行常规的处理和偏移等,通过利用雷达反射波信息,获得能够反映出地下异常体和地层界面位置的雷达剖面图。随着经济社会的发展,探地雷达的应用范围越来越广,对探地雷达的技术要求也越来越高,探地雷达成像不满足于仅获得界面位置,研究逐渐向能够获得地下具体电性参数的反演方法发展。传统的探地雷达层析成像方法基于射线追踪理论,如走时反演和第一周期最大振幅反演。基于上述射线技术提供的常规层析成像结果往往只能具有有限的分辨率,主要是因为在反演中仅利用了雷达数据所包含的小部分信息(即初至时间和第一周期最大振幅),且该方法仅能应用于跨孔探地雷达数据。为了实现提高分辨率、反演地面采集方案的探地雷达数据、实现能够为地下介质的物性参数提供定量解释的目的,波形反演技术被引入到了探地雷达领域并逐渐成为了研究的热点。在地球物理领域,波形反演技术最早被应用于地震技术探测与开发中,形成了一套基于波动方程的波形反演方案。目前,探地雷达波形方面的研究多集中在基于传统目标函数的跨孔雷达时间域全波形反演,利用雷达波的全部信息进行反演。虽然利用全部信息反演,可以获得高分辨率的反演结果,但同时在反演过程中也面临着较高的非线性,使得反演陷入局部最小。为了解决传统全波形反演非线性问题,最常用的手段是为波形反演提供一个相对准确的初始模型,但由于地下是未知的,如何提供初始模型,避免陷入局部最小,是波形反演研究方面的一个热点。本文采用基于麦克斯韦方程组的探地雷达矢量矩阵波形反演方案,首先实现了传统时间域全波形反演技术,利用对模型参数求偏导并引入虚拟向量源的概念,整理了针对于探地雷达波形反演的基础理论,给出了梯度和步长公式。通过采用时间域高阶有限差分的CPML吸收边界来实现对模型的正演模拟,在迭代过程中采用共轭梯度法来实现目标函数的收敛。在利用合成数据对该套反演理论进行验证后,本文反演了两套在不同地区采集的实际数据。在反演过程开始之前,要对雷达信号进行必要的预处理,并对实际数据进行三维数据转换二维数据和源子波估计等。在实现利用时间域全波形反演技术对跨孔雷达数据反演之后,本文将时间域波形反演技术拓展到了地面采集雷达数据反演中。通过滤波器将观测数据与正演数据的频率成分进行分解,在迭代过程中按照从低频到高频的规律,依次对不同的频率成分进行反演,构建了一套在时间域迭代、在频率域分解的多频段多尺度反演策略。多尺度策略能够有效避免在利用不准确的初始模型进行反演时,所面临的高非线性问题,降低对初始模型的依赖性。通过多频段多尺度策略,本文在时间域对地面采集的探地雷达合成数据进行反演,分别单独反演了介电常数和电导率。为了探讨更多利用波形反演技术实现探地雷达地面数据反演的可能性,本文提出了相位波形反演方法,并采用相位积分多尺度反演策略反演地面采集的探地雷达合成数据。在构建探地雷达相位波形反演目标函数时,通过傅立叶变(FFT)换分别获得实际观测数据和正演数据的振幅谱和相位谱,使用实际观测数据的振幅谱代替合成数据振幅谱,即对观测数据和正演数据的振幅归一化以消除振幅的影响,实现反演相位的目的。相位波形反演可以恢复与全波反演相同的雷达信息。相位波形反演不需要知道源子波的振幅,且在处理实际数据时振幅是难以估计的;相位积分多尺度反演策略允许较差的初始模式,这对于缺少先验信息的波形反演是非常重要的。通过相位多尺度反演方法,本文成功反演了地面采集方案的探地雷达合成数据,分别实现对介电常数和电导率的反演。在合成数据模拟中,相位多尺度反演方法可以对异常体和地层的位置信息准确还原,且反演得到的电性参数较为准确。本文将包络波形反演概念引入到探地雷达领域,通过希尔伯特变换(Hilbert transform)计算获得实际观测数据和正演数据的包络算子。在地震波形反演领域,地震记录的包络起伏和衰减被认为携带着超低频(即源谱中最低频率以下的频率,ULF)信号,可以用来估计长波速度结构。然后,包络反演可以用于恢复介质的低波数分量(平滑背景),以便可以减少波形反演的初始模型依赖性。本文使用希尔伯特变换的属性,我们在时域中分离相位和幅度信息。通过求导的方法并引入虚拟向量源的概念,详细推导了包络波形反演的梯度公式,并给出了最优步长。通过对信号的频谱进行分析,确定包络信号具有还原信号低频信息的能力,包络波形反演有更低的非线性,对初始模型参数的选择并不敏感,且具有抑制局部最小的能力。通过合成数据验证了包络目标函数有能力反演低频缺失跨孔雷达数据。在确定了包络目标函数抑制局部最小的能力之后,本文将包络目标函数与不依赖源子波的跨孔雷达波形反演相结合,构建了基于包络目标函数的不依赖源子波波形反演理论。利用对模型参数求偏导并引入虚拟向量源的概念,获得了基于包络目标函数不依赖源子波波形反演的梯度和步长公式,实现了介电常数和电导率的同步反演。根据褶积波场的滤波特性,通过改变反演过程中正演源子波的中心频率,制定了包络不依赖源子波波形反演的多尺度反演策略。避免了在反演过程中需要源子波估计这一复杂而又耗时的过程,在处理实际数据时,通过源子波估计获得的每一道波形往往区别较大,通常采用抽取其中相近的几道源子波平均的方法获得最终估计源子波,这为反演带来了不确定性,对收敛带来不利的影响。而包络不依赖源子波波形反演既避免了源子波估计这一过程,又抑制了常规不依赖源子波波形反演的高非线性。本文通过合成数据验证了包络不依赖源子波的波形反演方法,并成功利用该方法反演了实际数据。在文章的最后,本文针对不同的目标函数分别计算了它们的目标函数二维分布图。通过观察二维分布图,可以了解每一种方法对于不同物性参数初始模型的反演能力,初始模型参数选择的规律,每一种方法的特点。通过比较不同方法的目标函数二维分布图的形态,光滑程度与复杂程度,可以了解不同方法的非线性程度,对于构建多尺度反演方法有一定指导意义。
其他文献
农田灌溉水有效利用系数综合反映了农田灌溉水平,是国民经济和社会发展以及实现水利现代化的一项重要标准,同时也是国家实行最严格水资源管理制度“三条红线”的主要考核指标。
近年来国内外已研发出多种营养功能性面条。本文对近年来研制的营养功能性面条进行了归纳总结,分别对果蔬类、杂粮类、水产类、其他类等4大类营养功能性面条的开发应用和研究
本文通过对陕西省科技创新与金融创新耦合程度的现状及对接瓶颈进行分析的基础之上,提出促进陕西省科技创新与金融创新耦合模式发展的建议。旨在通过对陕西省科技创新与金融
现代新儒家第三代代表人物对儒学现代处境的评估,在上个世纪八十年代以后,已摆脱了其前辈之悲 观情绪的影响,变得乐观起来,以至于认为儒学似已走出现代困境,由"一阳来复"开始
研究目的:高盐(high salt, HS)饮食不仅是高血压的致病因素,还能独立于血压水平造成靶器官损伤。越来越多的证据显示免疫系统在高血压及其靶器官损伤发生中发挥重要作用,其中
本文以八十年代提出的“社会主义新人”形象塑造为研究对象,通过考察历史转折时期文坛多种力量的矛盾冲突与复杂纠缠,对这一富有历史特色而极具张力的文艺主张进行“重返现场
目的探究气压治疗与康复训练治疗脑卒中肩手综合征Ⅰ期患者的临床效果。方法选择2015-10—2016—10间南雄市人民医院进行治疗的60例脑卒中肩手综合征Ⅰ期患者,随机分为2组,其
<正>在后补贴时代,打车软件仅仅一个打车功能显然无以为继。综合来看,加强司机端的技术和产品优化,乘客端生活服务内容的扩展、更高级的商务体验才是比拼的要点。免费当道的
肌纤维是肌肉的基本组成物质。研究表明,肌纤维类型是决定肉品质的一个很重要的因素,因此,了解肌纤维类型的变化规律及调控机制至关重要。肌纤维的类型在生长过程中不断发生
王钦(以下简称“王”):甘阳老师,您好,很高兴您能接受我们的采访。我们就从刚刚结束的第三期“通识教育讲习班”谈起吧。我发现这次选择的阅读篇目中有章太炎的《馗书》和梁漱溟的