[发明专利]基于改进迭代法的地空时间域电磁系统高精度下延拓方法有效
申请号: | 201810878106.X | 申请日: | 2018-08-03 |
公开(公告)号: | CN109085652B | 公开(公告)日: | 2019-12-06 |
发明(设计)人: | 关珊珊;姜曜;黎东升;嵇艳鞠 | 申请(专利权)人: | 吉林大学 |
主分类号: | G01V3/38 | 分类号: | G01V3/38 |
代理公司: | 21229 沈阳维特专利商标事务所(普通合伙) | 代理人: | 杨群<国际申请>=<国际公布>=<进入国 |
地址: | 130012 吉*** | 国省代码: | 吉林;22 |
权利要求书: | 查看更多 | 说明书: | 查看更多 |
摘要: | |||
搜索关键词: | 感应电动势 地表 迭代法 延拓 电磁数据 电磁系统 迭代公式 时间域 改进 修正 传感器测量 测量数据 基线校正 接收线圈 精度确定 深度计算 视电阻率 算法确定 鱼群 采样 预设 校正 叠加 验证 | ||
1.基于改进迭代法的地空时间域电磁系统高精度下延拓方法,其特征在于,包括如下步骤:
1)计算长接地导线源在地表z=0平面的感应电动势,接收点之间需等间距分布,利用FFT向上延拓理论计算z=h h>0平面的空中感应电动势;
2)采用鱼群智能方法确定改进迭代法中的三个参数:α、β、γ,具体确定过程为:
通过鱼群方法选择最优的α、β、γ三个参数值,首先完成相对误差目标函数的构造:
η=MAX{[un(x,y,z=h)-f(x,y,z=h)]/un(x,y,z=h)×100} (1)
其中f(x,y,z=h)为步骤1)通过FFT向上延拓出的z=h平面感应电动势,un(x,y,z=h)为向下延拓的地表感应电动势再进行FFT向上延拓得到的空中感应电动势,进行鱼群算法的初始化,其中包括α、β、γ三个参数值的初始化,通过预设最大尝试次数M、精度值Q、逐步迭代确定最优参数值,首先判断η<Q,若满足条件迭代结束,否则按照最大尝试次数M完成最优参数确定;
并且,改进迭代法计算公式为:
3)通过改进迭代法,下延拓出z=0平面感应电动势,与步骤1)计算出的地表理论感应电动势进行比较,验证方法的正确性;
4)采用接收线圈传感器测量z=h平面感应电动势,为保证下延拓精度,测线条数最少为五条,且保持平行;
5)将测得数据进行叠加、采样处理,判断处理后的数据是否为等间距,否则进行插值;
6)采用改进迭代法,下延拓出z=0平面感应电动势,计算视电阻率,视深度,并完成视电阻率-视深度成像。
2.按照权利要求1所述的基于改进迭代法的地空时间域电磁系统高精度下延拓方法,其特征在于:
所述步骤3)中,采用改进迭代法实现空中感应电动势的下延拓计算,进一步提高地-空电磁系统数据解释的精度,为了保证下延拓在高频时解的稳定性,不能采用FFT向下延拓方法,因此采用改进迭代法实现空中场的向下延拓,将步骤1)中延拓出的z=h平面感应电动势f(x,y,z=h)作为地表z=0平面的初始感应电动势,并将此次延拓出的z=h平面感应电动势un(x,y,z=h)与f(x,y,z=h)之间的差值作为误差项:
errorn=un(x,y,z=h)-f(x,y,z=h) (2)
对向下延拓出的地表感应电动势un(x,y,z=0)进行修正:
其中n为迭代次数,增加积分项和微分项可以在相同迭代次数时增加计算精度;
对un+1(x,y,z=0)进行向上延拓、计算误差项、修正,直到errorn小到可以忽略为止,将此时得到的un+1(x,y,z=0)作为最终的z=0平面感应电动势,完成空中感应电动势的向下延拓。
该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于吉林大学,未经吉林大学许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服】
本文链接:http://www.vipzhuanli.com/pat/books/201810878106.X/1.html,转载请声明来源钻瓜专利网。