没有合适的资源?快使用搜索试试~ 我知道了~
基于多偏移距探地雷达数据的包络-波形反演.docx
1.该资源内容由用户上传,如若侵权请联系客服进行举报
2.虚拟产品一经售出概不退款(资源遇到问题,请及时私信上传者)
2.虚拟产品一经售出概不退款(资源遇到问题,请及时私信上传者)
版权申诉
0 下载量 160 浏览量
2023-02-23
20:07:46
上传
评论 2
收藏 880KB DOCX 举报
温馨提示
试读
14页
基于多偏移距探地雷达数据的包络-波形反演.docx
资源推荐
资源详情
资源评论
1. 引言
探地雷达(Ground Penetrating Radar, GPR)是一种用于研究浅地表电磁特性的无损检测
技术,主要基于电磁波传播原理,通过激发并记录高频电磁波(20~1000 MHz)在地下介质
中(地下障碍物或不连续地质体)的反射、绕射等传播特性来探测电性结构分布。目前被广
泛应用于土木工程场地勘查
[1]
、地雷探测
[2]
、环境与水文地质研究
[3,4]
、农业评估
[5,6]
、考古
调查
[7,8]
、极地考察
[9]
以及航空航天
[10,11]
等前沿工程与科学研究领域。
探地雷达的分辨率和探测深度极易受到脉冲源频率和地下电导率分布的限制,其中,
高频电磁脉冲具有较高的分辨率,但其在地下介质中的衰减更快,因此穿透深度较小。根
据测量方式的差异,探地雷达可进一步被划分为共偏移距和多偏移距两种类型。相较于共
偏移距的测量方式,多偏移距测量可以灵活调整发射和接收天线之间的距离,捕获到的绕
射波信息较为可靠,有助于实现精确的速度分析以及对地下电性参数的定量估计。因此开
展基于多偏移距探地雷达数据的研究是十分必要的。
目前,探地雷达数据的解释工作已经从单纯推断目标体空间位置、形态、尺寸以及层
界面等信息,上升到对介质电性参数(介电常数、电导率)的反演估计和对近地表结构的精
细刻画。现行的反演方法主要可分为基于射线理论的层析成像方法和基于全波场理论的全
波形反演方法。其中,全波形反演(Full-Waveform Inversion, FWI)最初起源于时域地震成
像,后被推广到频率域,其同时兼顾了波场的运动学和动力学特征(振幅和相位),是一种
高分辨率、高精度的成像方法。全波形反演已在地震勘探领域获得广泛应用,但由于地面
雷达测量方式对地下介质的照明能力有限,基于地面雷达数据的全波形反演具有挑战性。
Lavoué 等人
[12]
采用拟牛顿算法针对多偏移距地面雷达数据开展频率域 FWI,同时重构了介
电常数和电导率的 2 维分布。Feng 等人
[13,14]
将井间地震 FWI 和 GPR FWI 集成到一个联合
程序中,用于对不同地球物理参数(P 波速度、介电常数和电导率)模型的定量成像;Feng
等人
[15]
将多尺度反演策略和吉洪诺夫正则化(Tikhonov regularization, TV)约束方法结合应用
于时域 GPR FWI 过程,实现了介电常数和电导率的同步反演;王珣等人
[16]
在时间域采用
改进的全变差正则化策略,开展了 GPR 双参数多尺度同步反演研究;Huai 等人
[17]
在时间
域提出了一种基于模型的层剥离 FWI,有效增强了模型深部的重构质量和反演分辨率。
然而,当利用梯度类最优化算法求解基于 Born 近似框架下的全波形反演问题时,一
个严重偏离正解的初始模型会在波形匹配过程中产生“周波跳跃”问题,使得目标函数倾向
陷入到局部极值
[18]
。因此,全波形反演高度依赖精确的初始模型来确保反演过程的正确收
敛。低频成分对于恢复长波背景速度结构进而构建大尺度初始模型至关重要,但是实际采
集的数据中所携带的低频信息十分有限,这就使得研究低频缺失情况下的全波形反演方法
显得尤为重要。最初在地震勘探领域提出并发展起来的包络反演是用于重构低频成分比较
知名的一类方法。Bozdağ 等人
[19]
最早构建了基于瞬时相位和包络的全波形反演目标函数,
该种方法旨在利用 Hilbert 变换来实现波形相位和振幅信息的有效分离,从而降低了反问题
的非线性程度。Chi 等人
[20,21]
和 Wu 等人
[22]
构建了基于包络残差的全波形反演目标函数,
并利用伴随状态法求取了梯度算子,利用包络反演成功重建了地下介质的大尺度背景信
息。Wu 等人
[23]
提出了调制-褶积信号模型并在研究中赋予了包络数据以明确的物理意义。
刘新彤等人
[24]
研究了低频缺失下基于包络目标函数的全波形反演,在有效还原低频信息的
同时提供了对地下介质的定量解释,然而这一研究仅针对跨孔雷达的观测模式。与跨孔雷
达相比,地面雷达对地下介质的照明角度有限,这本身就增加了反问题的非线性和不确定
性;当地面雷达缺失低频成分的情况下,也会极大地制约全波形反演的稳定性和精确度。
因此,研究基于地面雷达数据的包络全波形反演具有重要意义。
此外,引用分频处理的多尺度反演策略也有利于缓解全波形反演的“周波跳跃”问题,
提高反演的稳定性。频率多尺度方法可按低频到高频依次反演多个离散频率或频率组,其
中,利用低频记录对局部极值的不敏感特性来重建平滑的大尺度背景结构,随后再利用高
频记录刻画局部精细结构
[25]
。Boonyasiriwat 等人
[26]
最初在地震勘探领域通过选择最优频带
和维纳滤波来实现更高效的时间域多尺度全波形反演;随后,这一策略被成功借鉴到探地
雷达的全波形反演中。Meles 等人
[27]
通过在模型更新过程中由低到高逐渐增加频率成分来
扩展数据的频率带宽,从而实现了对跨孔-井地数据的时-频域联合反演。Lavoué 等人
[12]
开
展了与频率采样有关的多参数全波形反演研究,通过数值实验证实了宽频带数据的同时反
演可以有效增强双参数同步重建质量。Li 等人
[28]
提出了一种基于跨孔雷达数据的 Laplace
域波形反演,并提出了频率跳跃的多尺度波形反演策略,实现了反演过程的快速收敛。
本文在现有的探地雷达时间域全波形反演框架下,回顾了地震勘探中信号包络提取和
常规包络反演方法的基本原理与实施流程,并将其成功推广到探地雷达领域,提出了基于
地面多偏移距雷达数据的包络-波形反演方法,实现了对地下介电常数模型的有效重构。采
用 3 层结构介电常数模型仿真模拟了数据低频信息缺失的情况,验证了包络-波形反演在增
强数据所携带的大尺度构造信息方面的能力与优势,有效克服了常规全波形反演在缺失低
频信息时无法正常工作的缺陷。此外,在所构建的包络-波形反演框架下,又有针对性地加
入了频率多尺度策略,在重建大尺度宏观背景构造的基础上,也有利于实现对小尺度结构
的精细刻画。本文还进行了抗噪性测试并验证了包络-波形反演在噪声条件下仍具有较强的
反演能力。
2. 包络-波形反演
2.1 探地雷达信号包络提取
一个解析信号[Math Processing Error]y~(t)是由一个实信号[Math Processing Error]y(t)
和该实信号所对应的希尔伯特变换[Math Processing Error]yH(t)组成的
[Math Processing Error]y~(t)=y(t)+jyH(t)
(1)
其中,[Math Processing Error]yH(t)可以表示为
[Math Processing Error]yH(t)=H{y(t)}=−1πP∫−∞+∞y(t)t−τdτ
(2)
其中,P 代表柯西主值。
解析信号[Math Processing Error]y~(t)也可以通过瞬时振幅[Math Processing Error]E(t)
以及瞬时相位[Math Processing Error]ϕ(t)来表示
[Math Processing Error]y~(t)=E(t)ejϕ(t)
(3)
其中,瞬时振幅[Math Processing Error]E(t)可表示为
[Math Processing Error]E(t)=ℜ{y~(t)}2+ℑ{y~(t)}2=y(t)2+yH(t)2
(4)
其中,[Math Processing Error]E(t)也称为实信号[Math Processing Error]y~(t)的包络。
值得注意的是,仅允许在时域而非频率域中将信号分离为上述相位和振幅的形式。
因此,探地雷达信号的包络也可以通过求取解析信号的振幅而获得,即
[Math Processing Error]Esyn=d(t)2+dH(t)2
(5)
其中,[Math Processing Error]Esyn 表示探地雷达合成数据的包络,d(t)表示探地雷达
合成数据, d
H
(t)则表示合成数据的希尔伯特变换。由希尔伯特变换的定义可知:信号包络
的提取实际上是一种非线性运算过程,这是因为 d
H
(t)相当于是将合成数据 d(t)的相位旋转
90°而获得的结果。
图 1(a)和图 1(b)分别给出了主频为 80 MHz 的源子波及其包络的波形与频谱数据;接
下来,利用高通滤波器剔除了低于 40 MHz 的全部以及 40~80 MHz 的部分低频成分,缺
失低频的子波及其包络的波形与频谱数据分别如图 1(c)和图 1(d)所示。从图 1(d)不难看
出:尽管源子波缺失了低频成分,但是其包络数据仍包含丰富的低频信息。
剩余13页未读,继续阅读
资源评论
罗伯特之技术屋
- 粉丝: 3651
- 资源: 1万+
下载权益
C知道特权
VIP文章
课程特权
开通VIP
上传资源 快速赚钱
- 我的内容管理 展开
- 我的资源 快来上传第一个资源
- 我的收益 登录查看自己的收益
- 我的积分 登录查看自己的积分
- 我的C币 登录后查看C币余额
- 我的收藏
- 我的下载
- 下载帮助
最新资源
- 基于matlab实现 powell算法 用matlab实现,使用方法内附详细说明.rar
- 基于matlab的手写字体识别程序,并对结果进行保存.rar
- 基于c语言指纹识别demo代码 包括了指纹图像方向图计算、频率计算、gabor滤波器增强,细化,特征点提取,特征点匹配.rar
- 基于c++NSGA-2思想的多目标优化程序,采用进化算法处理多目标实值优化问题.rar
- Linux系统中常用权限管理命令
- Coello Coello等人提出了MOPSO 该程序基于matlab实现针对测试函数matlab程序
- Linux系统中常用权限管理命令
- AIR-AP1815-K9-ME-8-5-182-0.tar For Cisco AP1815
- 实验七.zip
- ESP8266刷固件软件flash-download-tools-v3.6.5,AT固件,机智云固件
资源上传下载、课程学习等过程中有任何疑问或建议,欢迎提出宝贵意见哦~我们会及时处理!
点击此处反馈
安全验证
文档复制为VIP权益,开通VIP直接复制
信息提交成功