GprMax2D ——建模示例

混凝土中的单钢筋

这个简单问题的原理图如图所示。对于半径为r = 25毫米的完美导弹钢筋，需要一个探地雷达模型，该钢筋位于厚度为600毫米，深度为75毫米的混凝土板中。

第1步

首先，应确定所涉媒体的构成参数。由于可以在不需要指定其参数的情况下对完美的导体和自由空间进行建模，因此仅需要定义混凝土的本构参数。

让我们假设具体可以使用 Er = 6.0 和 σ = 0.01 （S / m）建模。而且，这可以被认为是非磁性的。因此， μr = 1.0 和 σ * = 0.0。0因此，输入文件中应包含：

#medium: 6.0 0.0 0.0 0.01 1.0 0.0 concrete

具体的介质标识符concrete 可以在模型中使用。

第2步

确定源类型和激发频率。这些应该是众所周知的，包含在问题的规范中。我们假设可以使用一个ricker源来模拟GPR天线，并且需要f = 900 MHz 的中心频率。为了控制所述nu-merical dispersion（数值色散）的空间的步长Δ x和Δ y应适当地选择。根据rule of thumb（经验法则），模型中感兴趣的最高频率的波长应该通过至少10个单元来确定，这个最感兴趣的频率不是用于激发脉冲的中心频率。通过检查源波形的幅度谱可以很容易地确定这一点，如图所示，很明显，脉冲在比中心频率更高的频率处包含大量能量。通常，应该在计算中使用的最高频率的良好估计是脉冲中心频率的三到四倍。因此，我们假设在脉冲频谱中具有重要意义的最高频率是fm = 3f = 2700 MHz，混凝土中2700 MHz的波长为：

因此λ ≈4.5 cm。考虑到经验法则，空间步长可以设置为 4 .5 毫米。

第3步

为了首先设置用于空间步长的值，判定该模型具有由 Δ x =Δ y=Δ l。这是一个显而易见的选择，因为不需要在特定方向上有较小的步长。如果我们选择为 Δ l = 4 .5毫米，然后，钢筋的半径 r ≈ 6个单元。这不会被认为是非常粗糙的，但如果该值增加到 10 ，它将改善模型。考虑到这一点，空间步骤设置为

 Δ x =Δ y =Δ l = 2 .5 毫米

Δl= 2.5毫米的选择完全在经验法则所描述的限度内。

第4步

由于  Δ x =Δ y = 2 .5 毫米的时间步长 Δ t 的计算是

等到  Δ t = 5.8967 皮秒（ps）。对于 8 纳秒的时间窗口，模型将需要大约 1357 次迭代。

第5步

让我们进一步假设需要42.5厘米长的扫描，包括41个GPR迹线。为了确定模型的大小，我们首先确定来自板坯终端的所有反射是否重要。考虑到它的深度和时间窗口并非如此，所以混凝土板可以建模为半空间。因此，600×300毫米的模型应该足够了。 根据Δx和Δy的值，可以转换为240x120个单元。 因此在输入文件中我们添加：

#domain: 0.6 0.3

#dx_dy: 0.0025 0.0025

#time_window: 8.0e-9

从y方向的300毫米开始，让我们使用50作为自由空间（在平板的顶部），因此可以通过添加到输入文件来定义平板：

#box: 0.0 0.0 0.6 0.25 concrete

然后引入圆筒以便用其特性覆盖平板的特性。 因此我们在输入文件添加：

#cylinder: 0.3 0.175 0.025 pec

第6步

需要指定由41个GPR迹线组成的长度为42.5厘米的扫描。如果扫描的中心在钢筋上，并且我们假设发射器和接收器之间的间隔为50毫米，则第一个源应为（0.075,0.2525），第一个接收器应为（0.125,0.2525）。源和接收器的高度设置为距接口2.5毫米。源和接收器在x 方向（扫描方向）上移动的步长是10 毫米。

从版本2开始，在引入将计算扫描的分析步骤之前，必须包含源定义。因此，源定义的命令将在文件中插入：

#line_source：1.0 900e6 ricker MyLineSource

要计算扫描，请在输入文件中插入以下分析步骤

#analysis：41 rebar.sca b

#tx：0.075 0.2525 MyLineSource 0.0 8.0e-9

#rx：0.125 0.2525

#tx_steps：0.01 0.0

#rx_steps：0.01 0.0

#end_analysis：

输出文件名为rebar.sca，将以二进制格式存储。

该模型现已完成，但我们可以在输入文件中添加更多命令：

#title: Model of rebar in concrete

#messages: y

#geometry_file: rebar.geo

以获取几何文件并包含模型的标题。最终输入文件如下：

#medium: 6.0 0.0 0.0 0.01 1.0 0.0 concrete

---------------------------------------------------------------------------

#domain: 0.6 0.3

#dx_dy: 0.0025 0.0025

#time_window: 8.0e-9

---------------------------------------------------------------------------

#box: 0.0 0.0 0.6 0.25 concrete

#cylinder: 0.3 0.175 0.025 pec

---------------------------------------------------------------------------

#line_source: 1.0 900e6 ricker MyLineSource

#analysis: 41 rebar.sca b

#tx: 0.075 0.2525 MyLineSource 0.0 8.0e-9

#rx: 0.125 0.2525

#tx_steps: 0.01 0.0

#rx_steps: 0.01 0.0

#end_analysis:

---------------------------------------------------------------------------

#title: Model of rebar in concrete

#messages: y

#geometry_file: rebar.geo

绘制线条是为了分隔输入文件中的不同步骤和部分。 文件可以包含任何其他注释，如果更改命令的顺序，则不会产生任何影响。

使用GprMax2D和上面的输入文件，已经获得了模拟的GPR扫描， 结果如图所示，同时使用存储在几何文件中的信息构建模型空间的图像表示。