探地雷达反问题的小波-同伦混合反演算法 Wavelet-Homotopy Hybrid Algorithm for Inverse Problem of Ground Penetrating Radar.pdf

计 算 物 理 第26卷第1期 V01．26．No．I CHINESE OFCOMPUTArrIONALPHYSICS 2009年1月 JOURNAL Jan．，2009 [文章编号】 1001．246XI2009)01．0087—07 探地雷达反问题的小波一同伦混合反演算法 李 壮h2， 韩 波1， 陈 勇1 (I．哈尔滨工业大学数学系，黑龙江哈尔滨150001；2．琼州学院计算机系，海南五指山572200) [摘 要】 构造一种求解探地雷达反问题的小波一同伦混合反演算法．利用小波将反问题分解到不同尺度上， 在最大尺度上采用同伦反演方法．求得次级尺度的初始近似解，在其余尺度上进行多次迭代修正，以获得全局 最小点．算法结合小波反演和同伦反演的优点，数值结果表明其具有较强的全局搜索能力． [关键词】 探地雷达；Maxwell方程；小波；同伦；混合反演 【中图分类号】P315 【文献标识码] A 0 引言 探地雷达(Ground Radar，简称GPR)是一种针对地下或物体内不可见目标的电磁探测技术，其 Penetrating 电磁模型的理论基础是Maxwell方程边值问题．在工程物探实际中，地下介质参数往往是未知的，探地雷达 反问题除完成常规电磁计算(正演)外，还要对未知介质参数进行反演计算．由于反问题具有不适定、非线性 和计算量大的特点，寻求稳定抗噪、收敛域宽而且计算量较小的实用算法，就成为反演计算中最为重要的研 究内容之一． 小波多尺度反演利用小波分解算法，将原始反问题分解在不同的尺度上，先在大尺度上迭代反演，得到 一个较好的参数估计，再将这个估计作为次一级尺度上的初始猜测，继续迭代修正，直到反演出原问题的全 局极小点．在小波多尺度反演的研究中，Dahlk…等针对带有线性算子的热传导方程逆问题，研究了如何将自 适应小波方法应用于反问题的一些思想；Cohen【2o等提出了基于小波一伽辽金离散格式求解反问题的两种数 值方法，“u旧1等用小波多尺度反演方法数值求解了一类椭圆型方程的分布参数识别问题；Michel¨1提出了一 种使用尺度函数和小波函数求解积分微分方程的小波多尺度正则化反演方法，并将其用于反演地球密度分 布参数．尽管取得一些成果，但尚存某些重要问题需待解决，其中之一即是，如何处理在最大尺度上对全局 极小的求解? 同伦方法是代数拓扑学中的一个基本概念，是求解非线性算子方程的一种快速稳定的收敛方法．其基 Jegenl8’等人将其应用于求解地震勘探及电磁勘探反问题，在国内由韩波【91首次应用到石油测井反演．同伦 反演做为一种高效实用的大范围收敛方法，在收敛范围和抗噪性方面表现出相对优势，但也存在一些问题， 如在多解或目标函数存在大量局部极小的情况下，如何保证其全局收敛性? 演与同伦反演优点、并能克服它们各自弱点的新型混合反演格式，结果表明，新算法能够在较短时间内得到 问题的最优近似解，是有效可行的． 1 基本理论 1．1 GPR反问题模型 三维情况下关于电场层和磁场日的GPR模型为 [收稿日期]2007—09—17；[修回日期】2008一Ol一21 【基金项目]国家自然科学基金及海南省自然科学基金(10701)资助项目 [作者简介】李壮(1967一)，男，黑龙江鸡西，副教授，博士，从事非线性反演算法研究 万方数据万方数据 88 计 算 物 理 第26卷 rotE：一P箬，rot日=s薯+扭“ (1) 其中平面(菇，)，)平行于地面，z轴垂直于地面并指向地下．t为时间，e，盯，／1分别为介电常数、电导率和磁导 F‘m～， C·m～．