关于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"."全孔径接收& ...
最新文章
- Set和存储顺序深入探讨、SortedSet排序的示例
- MATLAB表白利器
- 如何将Sql Server中的数据表导入到PowerDesigner中
- Thinking in Java,Fourth Edition(Java 编程思想,第四版)学习笔记(二)之Introduction to Objects...
- 【Android 界面效果9】9patch图片
- css3-响应式布局
- 对hash签名失败_详解Vue开发微信H5微信分享签名失败问题解决方案
- [css] 除了可以用js跟踪用户信息外,如果不用js,使用纯css怎么做呢?
- 会议 | 2018年全国知识图谱与语义计算大会(CCKS 2018)
- 【插件发布】JAVA微服务框架,Jeecg-P3-Base-System 1.0.0 插件开源发布
- MySQL if函数使用详解
- 吴恩达机器学习作业5.偏差和方差
- 更改后缀为.dex文件为.odex文件 让你的程序瘦身 运行更稳定
- 面试准备——springboot相关
- spring mvc 实现websocket(服务器主动消息推送)
- 非常经典的java编程题全集-共50题(1-10)
- Redis底层数据结构-SDS
- 奇兔recovery卡刷教程_如何使用recovery进行双清教程
- 万亿估值来了!3次大难不死,蚂蚁金服终于登顶世界第一
- 程序员的价值观与网络的复杂性
热门文章
- php mysql 框架 php5.2_简单快速安装Apache+PHP+MySql服务环境(四)—— 将php版本升级到7.2...
- Ubuntu 16.04 更新源
- 方阵和的行列式、方阵行列式的和
- http参数修改以及拦截
- LeetCode 295. Find Median from Data Stream
- 怎样练习一万小时(转)
- python格式化输出宽度对齐_python 中 字符串格式化示例,输出格式没有达到指定宽度,求帮忙!...
- 计算型存储:异构计算的下一个关键应用
- 抢鲜体验:openGauss从源码到主备
- 今晚直播:WLS/WAS故障基本分析介绍