GprMax2D ——英国建筑研究机构（BRE）示例 #2

BRE example #2

此示例表示的问题是来自第一个简单钢筋示例的更复杂的建模方案。该问题的示意图如图所示。遵循与前一示例中相同的逻辑步骤，可以确定模型的基本参数。手头的问题包括两种不同类型的材料 - 混凝土和湿砂。它们的介电特性的特定值没有特别限定。因此，假设湿砂的相对介电常数是Er = 20，并且混凝土的参数将是前面例子中使用的参数。湿砂的Er = 20的值可能相当高了，但是，这里我们假设它被水饱和。使用Ricker子源在F = 900兆赫兹，通过简单的计算表明的空间步长Δ l = 2 .5毫米，接近于经验法则设定的极限， 但它应该足够小，以获得相当准确的模型。

模型的大小设置为2600x800毫米，转换为1040x320个单元。时间窗口设置为12纳秒，并且沿着2400毫米的扫描线得115个GPR迹线。因此可以创建此模型的输入文件并读取

#medium: 6.0 0.0 0.0 0.01 1.0 0.0 concrete

#medium: 20.0 0.0 0.0 0.1 1.0 0.0 wet_sand

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

#domain: 2.5 0.5

#dx_dy: 0.0025 0.0025

#time_window: 12e-9

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

#box: 0.0 0.0 2.5 0.45 wet_sand

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

Use a cylinder of free space and then put a slab of concrete to cut it in half.

#cylinder: 1.05 0.3 0.2 free_space

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

#box: 0.0 0.3 2.5 0.45 concrete

#cylinder: 0.25 0.375 0.0125 pec

#cylinder: 0.45 0.375 0.0125 pec

#cylinder: 0.65 0.375 0.0125 pec

#cylinder: 0.85 0.375 0.0125 pec

#cylinder: 1.05 0.375 0.0125 pec

#cylinder: 1.25 0.375 0.0125 pec

#cylinder: 1.45 0.375 0.0125 pec

#cylinder: 1.65 0.375 0.0125 pec

#cylinder: 1.85 0.375 0.0125 pec

#cylinder: 2.05 0.375 0.0125 pec

#cylinder: 2.25 0.375 0.0125 pec

#triangle: 0.25 0.3 0.65 0.3 0.45 0.1 free_space

#triangle: 1.45 0.3 1.85 0.3 1.65 0.1 free_space

#box: 2.05 0.1 2.25 0.3 free_space

#cylinder: 2.15125 0.10125 0.1 wet_sand

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

#line_source: 1.0 900e6 ricker MyLineSource #analysis: 115 bre2.out b

#tx: 0.0875 0.4525 MyLineSource 0.0 12e-9

#rx: 0.1125 0.4525

#tx_steps: 0.02 0.0

#rx_steps: 0.02 0.0 #end_analysis:

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

#geometry_file: bre2.geo #title: BRE Model 2 #messages: y

检查对象构造命令在输入文件中出现的顺序很有趣。 例如，为了构造半圆形空隙，在构造湿砂的背景矩形之后使用#cylinder：命令。 湿砂矩形的垂直高度可以包括混凝土板的高度。 在第一种情况下，当引入板坯时，湿砂的性质将被混凝土的性质所取代。 在第二种情况下，可以将板坯放置在湿砂矩形上方，覆盖用于初始化模型的介质的值，该模型始终是自由空间。 第一种方法是优选的，因为它确保如果在板的坐标上产生错误，则不存在具有自由空间属性的任何节点。

在下图中，显示了存储在几何文件中信息的图像表示。观察该图像，模型的不同尺寸相对于其示意图是显而易见的。首先，模型的宽度比原理图中的宽度多100毫米。这对于确保ABC和第一个源点之间有足够的空间是必要的。其次，模型的一部分代表自由空间，第三，模型的高度小于示意图。在模型中没有必要包括所有600毫米的湿砂。因此，允许目标和ABC之间的一些空间可以忽略湿砂占据的其余空间。这样做是因为假设反射仅在模型内产生。如果用其他介质终止湿砂是重要的，则应将其包括在模型中。在所有四个示例中都使用了相同的方法，并且将节省大量的计算机内存和运行时间。

该示例的GprMax2D仿真结果如下图所示。模拟的GPR轨迹使用灰度图像格式显示。要创建此图像，GprMax2D数据需导入MATLAB。