收藏本站
收藏 | 手机打开
二维码
手机客户端打开本文

跨孔雷达全波形层析成像反演方法的研究

吴俊军  
【摘要】:跨孔雷达全波形反演是一种使用全波形信息反演两钻孔之间地下信息的层析成像技术。常规的层析成像反演大部分采用射线追踪方法,其中基于初至时的射线追踪方法可以反演出速度剖面,基于最大振幅的层析成像可以反演出衰减剖面。在一定假设情况下,速度和衰减反演图像可以转换成介电常数和电导率图像。地下介质的物性参数是影响介电常数和电导率的主要因素,包括孔隙度,含水率,含盐率,金属矿含量等。但常规射线追踪方法有许多不足,究其原因主要是该方法仅采用了一小部分的信号信息。例如:射线追踪方法仅能解决大于信号主波长的异常体,分辨率近似为第一菲涅尔带的直径。该方法不能提供稳定的低速异常体反演结果。当目标体仅能在有限的几个方向被提供“照明”时,这些弱点就尤为明显。 波形反演方法的应用非常广泛,如地震技术探测与开发、医学成像、无损测试、隧道和地雷探测等。近几十年来,波形反演发展了一些能够准确成像的方法,以此来处理包括声波、弹性波、雷达波、微波、光波和x光波在内的不同波场数据,其中菲涅尔量、程函方程、散射场及全波形技术等波形估计成像反演方法在地震中得到了应用。另外在正演方面还包括在时间域和频率域的有限元方法和有限差分方法。 跨孔雷达全波形反演作为全波形反演的组成部分,目前仍处在高速发展的初期。全波形反演能够提供小于波长的分辨率,在理想状态下,能够达到二分之一到三分之一波长。应用全波形反演方法在分辨率上比应用射线方法能提高一个量级。由于全波反演方法考虑了记录道信号內的全部信息,故其能够反演出多种形态的异常,其中包括在传播中因能量衰减而消失的信息。全波反演方法能够精确地重建单个异常体的位置、尺寸以及电磁特性;能够确定背景场中亚波长尺寸的相邻的异常体的位置,尺寸;能够表明非均匀层状模型中两个对比度较大的界面的电磁参数的平均变化情况或随机变化情况,包括非均匀层状模型中的充水的隧道位置、尺寸和电导率以及亚波长尺寸的管道的位置、尺寸和电导率。 本文全面推导了全波形跨孔雷达层析成像反演方法,该方法利用雷达全幅度相位信息能够反演出地下高分辨率的介电常数和电导率图像,背景场介质可以是物性参数大、对比度非均匀的介质。本文通过基于局域网的分布式并行算法,有效的解决了巨量数据正演计算,从而实现了全波形跨孔层析成像方法在普通计算机上全面有效的运行,而之前这种成像方法仅能在超级计算机上进行计算。 在本文中首先建立了基于UPML吸收边界的TE模式下的FDTD二维正演算法。由于跨孔雷达测量方式中发射天线与接收天线极化方式为垂直于地表的z方向,因此选用FDTD正演算法时需要使用有Ex、Ez、及Hy的TE模式。 本文中,反演算法采用最速梯度法求解,通过应用包括时间维度在内的全波场信息与残场逆向传播的全波场信息乘积来计算梯度方向,且通过求取以步长为自变量的目标函数的极值确定步长公式。由于介电常数与电导率在量级上的区别,因此在确定迭代步长时需要对两种电性参数进行分别计算以提高收敛速度,这里需要分别使用不同的稳定因子。文中对介电常数及电导率分别采用对数化处理,并推导出新的梯度公式。采用参数对数化处理能够很好的提高收敛速度及稳定性,并且拓宽了异常与周边介质物性参数对比度。 在计算梯度时,需要对包含时间维度在内的电场和残场全波场信息进行储存运算,这需要超过几十G的内存空间以及多达数百次的FDTD正演。国外同行在计算该类别全波反演时均使用超级计算机,本文建立了一个基于局域网的小型PC机群,组织了超过20核CPU参与计算,实现了跨孔成像全波反演非超级计算机算法,使得该方法能够为普通PC用户提供服务。 通过均匀背景下直径为1m的圆柱形目标体的全波反演,了解了实测数据(合成数据)波形、逆时残场波形的形态及特点,分析了梯度求解中每个阶段各组成部分的形态特征。在步长计算中,给定了合适的稳定因子κ ε及κ σ。 本文中将源所在的点及水平方向左右各五点设置为背景场值,并在源所在点左右10点内做滑动平均,保证介质相对介电常数的连续性。解决了反演中由于源附近介电常数“跳动”引起的波形畸变问题,提高了算法稳定性。 在不降低分辨率的情况下,通过对反演模型采用抽稀的办法,将数据减小至原内存需求的六分之一,以解决大型模型对计算机内存的巨大需求。 通过构建的矩形目标体模型可以发现全波反演能够精确的反演出矩形目标体的棱边位置、电性参数等。而对于倾斜或垂直异常体,能够准确的反演出目标体上下表面,但对垂直目标体的中间区域不能得到很好的成像。通过分析对复杂模型加入实际地表的模型,可以了解当需要考虑地表影响时,在反演初始模型中应当将地表所在位置加入计算。通过在模型中加入背景噪声,可以发现垂直方向的小尺度噪声并不能影响目标体成像。通过建立水平相关长度及垂直相关长度各异的随机介质模型可以发现,在100MHz天线频率下,介质介电常数3-8之间时,能够清晰的明显的分辨出0.5m以上的水平薄层异常体,而对垂直薄层异常体,全波反演不能很好的反演出原始模型。 本文结合跨孔及VRP测量方式进行全波反演,全面分析了不同测量方式下,不同电场分量的波场及残差波场的情况,针对建立的复杂模型,通过结合跨孔及VRP观测系统的全波反演能够有效地提高圆柱形目标体的水平方向分辨率。 在使用介电常数及电导率同步反演时,本文提出将第一步单介电常数反演作为初始模型,有效的解决了同步反演收敛速度较慢的问题。并通过模型试验可以看出在同步反演时,电导率的反演效果要比介电常数反演效果好一些。另外,对存在薄层目标体或异常体时,不应当结合VRP测量方式,采用VRP测量方式会降低跨孔测量模式下水平薄层的识别能力。 在数值反演中,源的波形是给定的,而在实际测量情况下,天线发射波形往往并不知道,因此子波的波形估计非常有必要。本文分别介绍了将源作为未知量进行反演的直接估计方法,以及采用频率域反褶积进行波形估计的方法,并应用反褶积方法迅速的完成了子波波形估计,尽管造成了子波相位上与理论值的差异,但使用该子波估计方法能够有效的全波反演成像。


知网文化
【相似文献】
中国期刊全文数据库 前20条
1 王静,何挺,李玉环;基于高光谱遥感技术的土地质量信息挖掘研究[J];遥感学报;2005年04期
2 屈永华;王锦地;;遥感地表参数反演中的一种不确定性知识处理方法[J];北京师范大学学报(自然科学版);2007年03期
3 王艳姣;张培群;董文杰;张鹰;;基于BP人工神经网络的水体遥感测深方法研究(英文)[J];Marine Science Bulletin;2007年01期
4 许年春;赵明阶;吴德伦;;节理岩体应力波反演模型研究[J];岩土力学;2007年12期
5 郭洁滨;黄冲;王海光;孙振宇;马占鸿;;基于高光谱遥感技术的不同小麦品种条锈病病情指数的反演[J];光谱学与光谱分析;2009年12期
6 王晓玥;;富营养化水体中氮磷物质与光谱反射率的相关性分析[J];杭州师范大学学报(自然科学版);2009年06期
7 谢国清;鲁韦坤;李蒙;周伟民;余凌翔;李发荣;;基于光谱特征和HSV变换融合MODIS影像的滇池叶绿素a浓度监测[J];生态学杂志;2010年06期
8 沈晔娜;吕军;陈军华;陈丁江;史一鸣;;水源区河流非点源污染物入河量计算的水质方程反演方法[J];环境科学;2010年08期
9 李法礼;王鸿燕;曹学伟;任恩明;;有限元大地测量反演影响因素分析[J];科技信息;2010年23期
10 陈勇;韩波;;时间推移地震反演的连续模型与算法[J];地球物理学报;2006年04期
11 陆明华;郑忠;杨建礼;;叠前敏感弹性参数反演在MXZ地区的应用[J];石油天然气学报;2008年02期
12 邬国锋;崔丽娟;纪伟涛;;基于时间序列MODIS影像的鄱阳湖丰水期悬浮泥沙浓度反演及变化[J];湖泊科学;2009年02期
13 张荣忠;袁诚;;地震波形反演技术的应用[J];天然气工业;2009年02期
14 荀尚培;翟武全;范伟;;MODIS巢湖水体叶绿素a浓度反演模型[J];应用气象学报;2009年01期
15 杨煜;李云梅;王桥;乐成峰;孙德勇;黄昌春;;富营养化的太湖水体叶绿素a浓度模型反演[J];地球信息科学学报;2009年05期
16 盖利亚;刘正军;张继贤;;基于CHRIS数据的三峡坝区水体叶绿素浓度遥感监测[J];测绘科学;2010年04期
17 李娜;吕建升;Altermann W;;光谱分析在植被重金属污染监测中的应用[J];光谱学与光谱分析;2010年09期
18 梁亮;杨敏华;张连蓬;林卉;;小麦叶面积指数的高光谱反演[J];光谱学与光谱分析;2011年06期
19 张妍,薄立群,路兴昌;长春净月潭地物波谱数据库的研究与开发[J];遥感信息;2002年02期
20 马明东;江洪;刘世荣;朱春全;刘跃建;王金锡;;森林生态系统立地指数的遥感分析[J];生态学报;2006年09期
中国重要会议论文全文数据库 前10条
1 李细兵;范小平;;微动技术中运用邻里算法反演地表速度结构[A];中国地震学会第14次学术大会专题[C];2012年
2 李志伟;胥颐;郝天珧;刘劲松;;利用非线性方法反演琉球-台湾-吕宋地区的岩石层P波速度结构[A];中国科学院地质与地球物理研究所2007学术论文汇编(第六卷)[C];2008年
3 李帝铨;王光杰;底青云;王妙月;王若;;基于遗传算法的CSAMT最小构造反演[A];中国科学院地质与地球物理研究所2008学术论文汇编[C];2009年
4 徐升;张鹰;;长江口水深遥感反演模型研究[A];第十五届全国遥感技术学术交流会论文摘要集[C];2005年
5 史玉峰;陈健;;雷达遥感技术在森林资源管理与监测中的应用[A];第二届中国林业学术大会——S3 森林经理与林业信息化的新使命论文集[C];2009年
6 彭涛;李耀刚;武威;郭印;;深埋型岩溶勘察中的地球物理勘探技术[A];中国地球物理·2009[C];2009年
7 毛志华;黄海清;;CMODIS资料提取叶绿素a浓度的方法研究[A];第六届成像光谱技术与应用研讨会文集[C];2006年
8 张世荣;刘润泽;;瀑布沟水电站防渗墙质量无损检测技术[A];2004水利水电地基与基础工程技术——中国水利学会地基与基础工程专业委员会2004年学术会议论文集[C];2004年
9 曹俊兴;朱介寿;;跨孔电磁波层析成象的应用问题研究[A];1995年中国地球物理学会第十一届学术年会论文集[C];1995年
10 李小超;张玉君;崔志文;王克协;;Biot介质中井外爆炸点源与跨孔爆炸点源激发的井孔声场[A];中国地球物理学会第二十四届年会论文集[C];2008年
中国博士学位论文全文数据库 前10条
1 吴俊军;跨孔雷达全波形层析成像反演方法的研究[D];吉林大学;2012年
2 邵广周;多阶模式瑞利波频散特征与反演研究[D];长安大学;2009年
3 李翠琳;基于非线性贝叶斯理论的多模态界面波频散曲线反演研究[D];中国海洋大学;2011年
4 王凌;苹果树花期叶/冠N、P营养状况的卫星遥感反演研究[D];山东农业大学;2012年
5 郭荣文;贝叶斯MT反演的非线性和不确定度分析[D];中南大学;2011年
6 刘伟;植被覆盖地表极化雷达土壤水分反演与应用研究[D];中国科学院研究生院(遥感应用研究所);2005年
7 吕杰;基于机器学习和辐射传输模型的农作物叶绿素含量高光谱反演模型[D];中国地质大学(北京);2012年
8 丛丕福;海洋叶绿素遥感反演及海洋初级生产力估算研究[D];中国科学院研究生院(遥感应用研究所);2006年
9 张彦林;基于3S技术的山东省森林蓄积量定量估测研究[D];北京林业大学;2008年
10 李瑞;北方农牧交错带草地植被动态研究[D];北京林业大学;2008年
中国硕士学位论文全文数据库 前10条
1 王世龙;HZ地区地震资料反演研究[D];成都理工大学;2012年
2 韩波;大地电磁二维拟牛顿反演研究[D];中国地质大学;2012年
3 李伟;多波叠前AVA反演中的改进模拟退火方法研究[D];成都理工大学;2012年
4 赵文轲;坑(井)—地电磁测深视参数反演研究[D];中国地质大学(北京);2010年
5 李子;基于GST-MQ配点法的突发水污染事故反演模型研究[D];清华大学;2010年
6 金希;高分辨率SAR影像裸土信息提取及土壤含水量反演初探[D];浙江大学;2011年
7 范渭亮;毛竹林参数遥感定量反演与软件实现[D];浙江农林大学;2010年
8 程丽娟;矿区植被生化组分高光谱遥感定量反演[D];山东科技大学;2011年
9 刘艳;基于MODIS数据的积雪深度反演研究[D];武汉大学;2005年
10 王景然;基于地形校正月表矿物含量反演研究[D];吉林大学;2012年
中国重要报纸全文数据库 前4条
1 通讯员 夏瑛 宋娟 记者 曹颖;江苏开展反演路面状况研究[N];中国气象报;2011年
2 记者 江然 通讯员 林长城 张春桂;福建三项气象课题获省科技奖[N];中国气象报;2010年
3 记者 周敏飞 黄跃 特约通讯员 李家山;青年运河200员工日夜抢修[N];湛江日报;2008年
4 通讯员 陈爱萍;义长灌域管理局科研项目通过专家鉴定[N];巴彦淖尔日报(汉);2010年
 快捷付款方式  订购知网充值卡  订购热线  帮助中心
  • 400-819-9993
  • 010-62982499
  • 010-62783978