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

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

相关文章

Java中实现线程的创建和启动的方法

《Java中实现线程的创建和启动的方法》在Java中,实现线程的创建和启动是两个不同但紧密相关的概念,理解为什么要启动线程(调用start()方法)而非直接调用run()方法,是掌握多线程编程的关键,... 目录1. 线程的生命周期2. start() vs run() 的本质区别3. 为什么必须通过 st

C#之List集合去重复对象的实现方法

《C#之List集合去重复对象的实现方法》:本文主要介绍C#之List集合去重复对象的实现方法,具有很好的参考价值,希望对大家有所帮助,如有错误或未考虑完全的地方,望不吝赐教... 目录C# List集合去重复对象方法1、测试数据2、测试数据3、知识点补充总结C# List集合去重复对象方法1、测试数据

SpringBoot读取ZooKeeper(ZK)属性的方法实现

《SpringBoot读取ZooKeeper(ZK)属性的方法实现》本文主要介绍了SpringBoot读取ZooKeeper(ZK)属性的方法实现,强调使用@ConfigurationProperti... 目录1. 在配置文件中定义 ZK 属性application.propertiesapplicati

MyBatis设计SQL返回布尔值(Boolean)的常见方法

《MyBatis设计SQL返回布尔值(Boolean)的常见方法》这篇文章主要为大家详细介绍了MyBatis设计SQL返回布尔值(Boolean)的几种常见方法,文中的示例代码讲解详细,感兴趣的小伙伴... 目录方案一:使用COUNT查询存在性(推荐)方案二:条件表达式直接返回布尔方案三:存在性检查(EXI

Java调用C#动态库的三种方法详解

《Java调用C#动态库的三种方法详解》在这个多语言编程的时代,Java和C#就像两位才华横溢的舞者,各自在不同的舞台上展现着独特的魅力,然而,当它们携手合作时,又会碰撞出怎样绚丽的火花呢?今天,我们... 目录方法1:C++/CLI搭建桥梁——Java ↔ C# 的“翻译官”步骤1:创建C#类库(.NET

Java 枚举的基本使用方法及实际使用场景

《Java枚举的基本使用方法及实际使用场景》枚举是Java中一种特殊的类,用于定义一组固定的常量,枚举类型提供了更好的类型安全性和可读性,适用于需要定义一组有限且固定的值的场景,本文给大家介绍Jav... 目录一、什么是枚举?二、枚举的基本使用方法定义枚举三、实际使用场景代替常量状态机四、更多用法1.实现接

java String.join()方法实例详解

《javaString.join()方法实例详解》String.join()是Java提供的一个实用方法,用于将多个字符串按照指定的分隔符连接成一个字符串,这一方法是Java8中引入的,极大地简化了... 目录bVARxMJava String.join() 方法详解1. 方法定义2. 基本用法2.1 拼接

java连接opcua的常见问题及解决方法

《java连接opcua的常见问题及解决方法》本文将使用EclipseMilo作为示例库,演示如何在Java中使用匿名、用户名密码以及证书加密三种方式连接到OPCUA服务器,若需要使用其他SDK,原理... 目录一、前言二、准备工作三、匿名方式连接3.1 匿名方式简介3.2 示例代码四、用户名密码方式连接4

springboot项目中使用JOSN解析库的方法

《springboot项目中使用JOSN解析库的方法》JSON,全程是JavaScriptObjectNotation,是一种轻量级的数据交换格式,本文给大家介绍springboot项目中使用JOSN... 目录一、jsON解析简介二、Spring Boot项目中使用JSON解析1、pom.XML文件引入依

IDEA中Maven Dependencies出现红色波浪线的原因及解决方法

《IDEA中MavenDependencies出现红色波浪线的原因及解决方法》在使用IntelliJIDEA开发Java项目时,尤其是基于Maven的项目,您可能会遇到MavenDependenci... 目录一、问题概述二、解决步骤2.1 检查 Maven 配置2.2 更新 Maven 项目2.3 清理本