关于Gprmax正演模拟结果显示空白的原因分析
Gprmax正演模拟结果显示空白
gprMax正演模拟结果显示空白,与发射天线极化方向有关,电场方向应垂直剖面方向,与参数设置关系不是很大。二维 gprMax 显示空白与天线收发距有关。
文章目录
- Gprmax正演模拟结果显示空白
- 一、问题的提出
- 二、问题分析
- 三、问题解决
一、问题的提出
用Gprmax正演模拟时,偶然发现模拟出来的结果是空白的,结果如下:
但是波形图却有结果:
虽然波形不是很明显,但是可以看出部分结果。
这是我用来模拟的In文件
#title: 3Dkondon_100cm_20cm#material: 12.0 0.01 1.0 0.0 soil
#material: 1.0 0 1.0 0.0 air#domain: 4.000 2.100 2.100
#dx_dy_dz: 0.010 0.010 0.010
#time_window: 60e-9#box: 0.000 0.000 0.000 4.000 2.000 2.000 soil
#box: 1.750 1.000 0.750 2.250 1.500 1.250 air#waveform: ricker 1.0 500e6 my_ricker
#hertzian_dipole: z 0.100 0.950 2.000 my_ricker#rx: 0.140 0.950 2.000
#src_steps: 0.040 0.000 0.000
#rx_steps: 0.040 0.000 0.000#geometry_view: 0 0 0 4.000 2.100 2.100 0.010 0.010 0.010 PEC n
二、问题分析
首先怀疑是参数设置的问题,我修改了介质的介电常数,当异常体的介电常数与背景介电常数相差不大时可以得到正常的结果,当介电常数很大时,介质对电磁波能量吸收强,导致模拟出来的图片显示空白。
然而不同介质有固定的介电常数,不能随意更改。
拖动波形图,发现直达波振幅非常大,与反射波相差几个数量级。
更新:(2021年10月4日)
空白的原因与天线极化方向有很大关系,见下篇博客的分析。
更新:(2021年10月5日)
这篇文章用于模拟的模型如下,文章里面的图是测线在边缘模拟的结果,测线不在异质体上方,所以结果有点奇怪(更新一下,以免误导)
三、问题解决
去除直达波即可解决这个问题,发现去除直达波后双曲线反射弧即显示出来了。
去直达波的方法是直接令直达波振幅等于0。
上面分享的方法不一定是对的,但是有效果。
[^1]本文仅做学习交流。
关于Gprmax正演模拟结果显示空白的原因分析相关推荐
- ParaView绘制gprMax正演模拟的波场快照方法(1)
ParaView绘制gprMax正演模拟的波场快照方法(1) gprMax是一款优秀的基于时域有限差分方法(FDTD)的电磁波数值模拟脚本软件,其正演模拟的结果通过波场快照的形式可以直观的显示出来,通 ...
- gprMax 正演模拟中Ex、Ey、Ez三个分量之间的关系分析
gprMax 正演模拟中Ex.Ey.Ez三个分量之间的关系分析 在 GPR 应用中,电场分量通常是测得量.我们一般的正演模拟用哪个电场分量呢 文章目录 gprMax 正演模拟中Ex.Ey.Ez三个分量 ...
- gprMax电磁波正演模拟方法
文章首发于:https://blog.zhaoxuan.site/archives/37.html: 第一时间获取最新文章请关注博客个人站:https://blog.zhaoxuan.site. 目录 ...
- 声波正演c语言程序,二维频率域声波方程正演模拟
1. 概述 频率域波场正演相对于时间域数值模拟来说,有其自身的优势.首先,在多炮数值模拟情况下,频率域相对于时间域效率更高,每个频率成分的阻抗矩阵只需要计算一次,加入并行计算后可以极大地提高计算效率: ...
- 重力勘探正演模拟matlab,裴雪林, 郭万松 (1995) 高精度重力勘探技术在国内外的应用. 断块油田, 5, 8-11....
ABSTRACT: 本文典型地质模型三维重力正演模拟研究是在二维重力正演方法总结的基础上,对地质模型作三维重力正演探究.其中所建的地质模型以球体为主.分别导出了二维与三维重力布格异常,重力异常各阶导数 ...
- 基于Madagascar的二维地震声波波动方程正演模拟
最近在将SU写的地震勘探的程序迁移到Madagascar上,初步尝试,写了一个二维声波方程正演程序,很简单,也很基本,只能输出波场快照,没有吸收边界条件,贴出来,供大家参考.代码和脚本如下: #inc ...
- VOS客户端上服务器和软交换状态显示红灯的原因分析
VOS客户端上服务器和软交换状态灯显示红灯,实际后台服务都在运行的原因 登录VOS3000客户端 原因分析 VOS客户端连接到服务器是属于TCP连接,TCP连接一旦断开后,服务器和软交换的状态灯就会 ...
- 地震射线追踪与有限差分正演模拟小软件
翻看自己研究生时候用java写的小软件,功能还是很强大的 https://github.com/Rtoax/Seismic-Processing
- 基于CUDA的VTI介质有限差分正演模拟与逆时偏移及ADCIGs提取
简单明了"CUDA"."C语言"."nvcc"."VTI介质"."RTM"."全孔径接收& ...
最新文章
- JQuery中html()和val()的用法区别
- nginx安装-添加MP4播放模块
- 【杂题总汇】NOIP2013(洛谷P1967) 货车运输
- matlab中线性变压器,用MATLAB优化计算差动变压器的线性范围
- python画方波_python实现周期方波信号频谱图
- 細微之處看看mysql與sql server的一些差別
- 对SqlServer2008中的日志进行截断的方法
- python使用函数的头文件_Python3.6下在Python中调用C函数
- Python 进阶 —— 可变参数(*args, **kw)与参数收集的逆过程
- 从零开始撸一个Fresco之gif和Webp动画
- wpe手机中文汉化版下载_wpe中文版
- 单片机实验中用到的元器件(Multisim14)
- c语言给数组整体赋值,c语言给数组赋值有哪些形式
- 『.NET Core CLI工具文档』dotnet-publish
- HttpWatch使用教程
- whai is gradient vanishing and exploding ?
- 各种内部排序方法及其比较实验报告
- java 鼠标动眼睛动_java实现眼睛跟着鼠标转动的跟踪眼代码
- vba数据类型,运算符,内置函数,循环判断语句,窗体控件
- bzoj1615 【Usaco2008 Mar】The Loathsome Hay Baler 麻烦的干草打包机
热门文章
- bootstrap4高度占一半_减肥选对了碳水,意味着成功了一半
- 一些不错的sql语句
- [USACO09OCT]热浪Heat Wave
- java常见面试题——java常见笔试题
- Python3之requests模块
- 【异常】INFO: TopologyManager: EndpointListener changed ...
- 【转】SQLite3简介及在.Net程序中的使用(自增列的创建)
- java调用javafx_Java代码调用JavaFX的例子(大力推荐)
- mysql 表连接 on_MysQL中表连接中“using”和“on”之间的区别是什么?
- SQL必知必会-组合查询