探地雷达正演模拟,基于时域有限差分方法,一

2024-06-10 04:28

本文主要是介绍探地雷达正演模拟,基于时域有限差分方法,一,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

     声明:本博客中的公式均是在Word中使用AxMath写好后截图使用的,欢迎引用,但请标注来源。

本系列会有四篇博客:

第一篇内容:

1、基础知识掌握

2、Maxwell方法差分求解原理

第二篇内容:

1、基于C++的TE波波动方程实现

2、边界问题的产生及处理

第三篇内容:

1、基于实际操作流程的GPR正演模拟

2、简单并行化加速GPR正演模拟

第四篇内容:

1、GPR剖面获取及成像

2、简单数据处理及分析

   探地雷达(Ground-Penetrating-Radar,GPR)是一种使用高频电磁波探测近地表地下结构的地球物理勘探方法,该方法对地表无破坏,操作简便,已被广泛应用于考古、城市道路建设、岩土工程调查等领域,有着非常良好的使用效果。但是,当前公开的GPR实际资料较少且难以形成庞大的知识库作为先验性息以备研究人员对实测剖面进行识别,所以,使用GPR正演模拟获得一套资料就显得尤为重要了。

        当前,基于Maxwell方程的GPR正演模拟主要使用的方法有:

        1、有限元方法(论文和开源软件都有)

        2、时域有限差分方法(GPRMax软件已经集成,但该软件没有界面,使用需要通过命令行程序执行)

        其中,时域有限差分方法(Finite-Difference-Time-Domain,FDTD)由K.S.Yee于1966年提出,是一种非常成熟且概念简单的正演模拟方法,因此,本系列博客决定基于FDTD实现GPR正演模拟。

一、时域有限差分方法

        首先从Maxwell旋度方程及其本构方程讲起:

      Maxwell旋度方程:

Maxwell旋度方程

本构关系:

电磁场向量本构关系

在上述公式中,E是电场强度(V/m),H是磁场强度(H/m),D是电通量密度(C/m^2),B是磁通量密度(Wb/m^2),J是电流密度(A/m^2),介电常数\varepsilon _0=8.85*10^{-12}F/m,磁导常数\mu_0=4\pi *10e^{-7}H/m 。

        将旋度方程转换为差分式是比较容易理解的,这里为了方便表示,使用了矩阵形式表示二维微分方程形式的TE波波动方程:

TE波波动方程矩阵形式

其中,

参照参考文献1中的Yee元胞结构:

Yee元胞结构

        通过Yee元胞离散后,就可以对上述波动方程公式进行空间取样。

        这时要注意:虽然文献1中的差分取样是中心差分,该方式的好处是可以得到二阶计算精度(具体可以通过Taylor公式计算),但是计算机中无法通过上述方式实现,只能是在规定点取样,所以,理论上是二阶精度,但是,计算中仍然是一阶精度。

        首先是系数矩阵问题,公式详见参考文献1P18-P19:

系数矩阵

由n状态更新到n+1状态就可以:

具体下标差分方程

上述内容为使用FDTD实现GPR正演模拟所需要的基础知识,重新梳理一遍:通过Yee元胞对Maxwell方程进行离散,得到电磁场的更新方程。

参考文献:葛德彪,闫玉波. 电磁波时域有限差分法第2版.西安:西安电子科技大学出版社,2005:133-137,37,118.

Shen, H, Y., Li, X, S., Duan, R, F., et al., (2023), Quality evaluation of ground improvement by deep cement mixing piles via ground-penetrating radar. Nature Communications,14,34-48. https://doi.org/10.1038/s41467-023-39236-4

冯德山,戴前伟,何继善等. 探地雷达GPR正演的时域有限差分实现(英文)[J].地球物理学进展,2006,(02):630-636.

李静,刘津杰,曾昭发等. 基于变换光学有限差分探地雷达数值模拟研究[J].地球物理学报,2016,59(06):2280-2289.

这篇关于探地雷达正演模拟,基于时域有限差分方法,一的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/1047208

相关文章

SpringBoot日志配置SLF4J和Logback的方法实现

《SpringBoot日志配置SLF4J和Logback的方法实现》日志记录是不可或缺的一部分,本文主要介绍了SpringBoot日志配置SLF4J和Logback的方法实现,文中通过示例代码介绍的非... 目录一、前言二、案例一:初识日志三、案例二:使用Lombok输出日志四、案例三:配置Logback一

Python实现无痛修改第三方库源码的方法详解

《Python实现无痛修改第三方库源码的方法详解》很多时候,我们下载的第三方库是不会有需求不满足的情况,但也有极少的情况,第三方库没有兼顾到需求,本文将介绍几个修改源码的操作,大家可以根据需求进行选择... 目录需求不符合模拟示例 1. 修改源文件2. 继承修改3. 猴子补丁4. 追踪局部变量需求不符合很

mysql出现ERROR 2003 (HY000): Can‘t connect to MySQL server on ‘localhost‘ (10061)的解决方法

《mysql出现ERROR2003(HY000):Can‘tconnecttoMySQLserveron‘localhost‘(10061)的解决方法》本文主要介绍了mysql出现... 目录前言:第一步:第二步:第三步:总结:前言:当你想通过命令窗口想打开mysql时候发现提http://www.cpp

Mysql删除几亿条数据表中的部分数据的方法实现

《Mysql删除几亿条数据表中的部分数据的方法实现》在MySQL中删除一个大表中的数据时,需要特别注意操作的性能和对系统的影响,本文主要介绍了Mysql删除几亿条数据表中的部分数据的方法实现,具有一定... 目录1、需求2、方案1. 使用 DELETE 语句分批删除2. 使用 INPLACE ALTER T

MySQL INSERT语句实现当记录不存在时插入的几种方法

《MySQLINSERT语句实现当记录不存在时插入的几种方法》MySQL的INSERT语句是用于向数据库表中插入新记录的关键命令,下面:本文主要介绍MySQLINSERT语句实现当记录不存在时... 目录使用 INSERT IGNORE使用 ON DUPLICATE KEY UPDATE使用 REPLACE

CentOS 7部署主域名服务器 DNS的方法

《CentOS7部署主域名服务器DNS的方法》文章详细介绍了在CentOS7上部署主域名服务器DNS的步骤,包括安装BIND服务、配置DNS服务、添加域名区域、创建区域文件、配置反向解析、检查配置... 目录1. 安装 BIND 服务和工具2.  配置 BIND 服务3 . 添加你的域名区域配置4.创建区域

mss32.dll文件丢失怎么办? 电脑提示mss32.dll丢失的多种修复方法

《mss32.dll文件丢失怎么办?电脑提示mss32.dll丢失的多种修复方法》最近,很多电脑用户可能遇到了mss32.dll文件丢失的问题,导致一些应用程序无法正常启动,那么,如何修复这个问题呢... 在电脑常年累月的使用过程中,偶尔会遇到一些问题令人头疼。像是某个程序尝试运行时,系统突然弹出一个错误提

电脑提示找不到openal32.dll文件怎么办? openal32.dll丢失完美修复方法

《电脑提示找不到openal32.dll文件怎么办?openal32.dll丢失完美修复方法》openal32.dll是一种重要的系统文件,当它丢失时,会给我们的电脑带来很大的困扰,很多人都曾经遇到... 在使用电脑过程中,我们常常会遇到一些.dll文件丢失的问题,而openal32.dll的丢失是其中比较

python中字符串拼接的几种方法及优缺点对比详解

《python中字符串拼接的几种方法及优缺点对比详解》在Python中,字符串拼接是常见的操作,Python提供了多种方法来拼接字符串,每种方法有其优缺点和适用场景,以下是几种常见的字符串拼接方法,需... 目录1. 使用 + 运算符示例:优缺点:2. 使用&nbsjsp;join() 方法示例:优缺点:3

Mysql中深分页的五种常用方法整理

《Mysql中深分页的五种常用方法整理》在数据量非常大的情况下,深分页查询则变得很常见,这篇文章为大家整理了5个常用的方法,文中的示例代码讲解详细,大家可以根据自己的需求进行选择... 目录方案一:延迟关联 (Deferred Join)方案二:有序唯一键分页 (Cursor-based Paginatio