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

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

相关文章

Linux换行符的使用方法详解

《Linux换行符的使用方法详解》本文介绍了Linux中常用的换行符LF及其在文件中的表示,展示了如何使用sed命令替换换行符,并列举了与换行符处理相关的Linux命令,通过代码讲解的非常详细,需要的... 目录简介检测文件中的换行符使用 cat -A 查看换行符使用 od -c 检查字符换行符格式转换将

SpringBoot实现数据库读写分离的3种方法小结

《SpringBoot实现数据库读写分离的3种方法小结》为了提高系统的读写性能和可用性,读写分离是一种经典的数据库架构模式,在SpringBoot应用中,有多种方式可以实现数据库读写分离,本文将介绍三... 目录一、数据库读写分离概述二、方案一:基于AbstractRoutingDataSource实现动态

Java中的String.valueOf()和toString()方法区别小结

《Java中的String.valueOf()和toString()方法区别小结》字符串操作是开发者日常编程任务中不可或缺的一部分,转换为字符串是一种常见需求,其中最常见的就是String.value... 目录String.valueOf()方法方法定义方法实现使用示例使用场景toString()方法方法

Java中List的contains()方法的使用小结

《Java中List的contains()方法的使用小结》List的contains()方法用于检查列表中是否包含指定的元素,借助equals()方法进行判断,下面就来介绍Java中List的c... 目录详细展开1. 方法签名2. 工作原理3. 使用示例4. 注意事项总结结论:List 的 contain

macOS无效Launchpad图标轻松删除的4 种实用方法

《macOS无效Launchpad图标轻松删除的4种实用方法》mac中不在appstore上下载的应用经常在删除后它的图标还残留在launchpad中,并且长按图标也不会出现删除符号,下面解决这个问... 在 MACOS 上,Launchpad(也就是「启动台」)是一个便捷的 App 启动工具。但有时候,应

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