【电磁】求解瞬变电磁TEM层状介质正演【含Matlab源码 2164期】

最新推荐文章于 2024-08-06 15:39:39 发布

Matlab领域

最新推荐文章于 2024-08-06 15:39:39 发布

阅读量881

点赞数 22

分类专栏： Matlab物理应用（初级版）文章标签： matlab

本文链接：https://blog.csdn.net/weixin_63266434/article/details/135729730

版权

Matlab物理应用（初级版）专栏收录该内容

277 篇文章 10 订阅

订阅专栏

本文介绍了使用Matlab求解瞬变电磁层状介质正演问题的方法，包括基于大定源回线的电磁场理论，提供部分源代码示例，并展示了不同模型的运行结果。文中还讨论了Matlab版本和相关研究参考文献。

摘要由CSDN通过智能技术生成

在这里插入图片描述

⛄一、获取代码方式

获取代码方式1：
完整代码已上传我的资源：【电磁】基于matlab求解瞬变电磁TEM层状介质正演【含Matlab源码 2164期】
点击上面蓝色字体，直接付费下载，即可。

获取代码方式2：
付费专栏Matlab物理应用（初级版）

备注：
点击上面蓝色字体付费专栏Matlab物理应用（初级版），扫描上面二维码，付费29.9元订阅海神之光博客付费专栏Matlab物理应用（初级版），凭支付凭证，私信博主，可免费获得1份本博客上传CSDN资源代码（有效期为订阅日起，三天内有效）；
点击CSDN资源下载链接：1份本博客上传CSDN资源代码

⛄二、瞬变电磁

1 引言
瞬变电磁法的主要工作装置有中心回线、重叠回线、大定源回线等, 其中又以大定源回线装置应用最为广泛。大回线源装置的工作方法是采用较长的发送回线, 利用阶跃波电流场源激励, 在大地产生过渡过程场, 断电瞬间在大地中形成涡旋交变电磁场。在回线内外一定的范围内, 可观测到这种由地下介质产生的二次感应电磁场, 随时间变化的衰减特性, 进而从测量得到的异常信号中分析出地下不均匀体的导电性能和位置, 以达到解决地质问题的目的。大回线装置与目标体具有最佳耦合、异常幅值大、形态简单、受旁侧地质体影响小的特点, 并对高阻层有很强的穿透能力, 对低阻层有较高的分辨能力。这种装置对铺设回线的要求不是很严格, 一旦铺好回线后, 不仅可采用多台接收机同时工作, 而且可以把接收线圈排成阵列, 发展成阵列式接收的观测系统。这种场源具有发射磁矩大、场均匀及随距离衰减慢等特点, 适合于密集点距采样, 精细探测。

大回线由于面积较大, 不能看作磁偶极子。因此, 研究大回线源激发的电磁场有着重要的意义。国内、外对大回线源产生电磁场的研究并不多, 翁爱华等首先求出任意一层电磁波的反射系数和透射系数, 然后据此计算了圆形发射回生的频率域电磁场;M.Poddar根据电磁学互易原理和垂直磁偶源产生的电磁场;求解出距形发射回线产生的频率域电磁场;Nagendra Partap Singh等利用超几何分布函数, 化简了圆形发射回线电磁场表达式中贝塞尔函数的乘积, 从而求解出频率域电磁场。

作者从基本的电磁场理论出发, 给出大回线发射源框内外不同位置的电磁响应表达式, 利用一种新的算法计算电磁响应, 并与传统的线性数值滤波法计算结果进行了对比分析 (对比所用的滤波权系数, 采用Anderson计算的801个系数) 。

2 理论基础
在这里插入图片描述
每个区间上的积分, 采用以下形式计算:

式中m是积分阶数;w是与横坐标x有关的权重向量。

通过对上述方程做简单的变换, 可以变换成与FHT算法一样的形式:
在这里插入图片描述

⛄三、部分源代码

close all
clc
clear
%%% 受Matlab计算精度所限，观测时间范围设置在10^(-7)至10（-2）s以内。
%%% 受Matlab计算精度所限，各层的测试电阻设置在500欧姆以下。
%%% FORTRAN语言才是电磁波场正反演的经典编程语言，MATLAB可用于简单测试
t=logspace(-7,-2,100);
%两层模型测试
[V1]=TEMforward4([500],[]);
[V2]=TEMforward4([100],[]);
[V3]=TEMforward4([500,100],[50]);
%三层模型测试
[V4]=TEMforward4([200],[]);
[V5]=TEMforward4([50],[]);
[V6]=TEMforward4([500],[]);
[V7]=TEMforward4([200,50,500],[50,100]);
%四层模型测试
[V20]=TEMforward4([500,100,1000,200],[50,200,500]);
%一层模型绘图
V11=[V1;V2];
figure(1);
loglog(t,V11,‘.-’);
legend(‘500欧姆’,‘100欧姆’);
xlabel(‘时间 / s’)
ylabel(‘电压 / V’)
%两层模型绘图
V=[V1;V2;V3];
figure(2);
loglog(t,V,‘.-’);
legend(‘第一层’,‘第二层’,‘两层’);
xlabel(‘时间 / s’)
ylabel(‘电压 / V’)
%三层模型绘图
V12=[V4;V5;V6;V7];
figure(3);
loglog(t,V12,‘.-’);
legend(‘第一层’,‘第二层’,‘第三层’,‘三层’);
xlabel(‘时间 / s’)
ylabel(‘电压 / V’)
%四层模型绘图
figure(4);
loglog(t,V20,‘.-’);
legend(‘{\rho}=[500,100,1000,200] H=[50,200,500]’);
xlabel(‘时间 / s’)
ylabel(‘电压 / V’)