关于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]本文仅做学习交流。