Gprmax 三维地质雷达建模及在 paraview 中的可视化（1）

Gprmax 运行三维模型速度很慢，建议在有GPU加速的环境下运行，有关GPU加速的方法，在我的上一个博客link中有说明。

文章目录

Gprmax 三维地质雷达建模及在 paraview 中的可视化（1）
- 一、in文件
- 二、模型可视化
- 三、仿真结果

一、in文件

这里，我们想建立一个长 4 m ，宽 2 m，深 2 m 的研究区域，填充物质为土壤，查阅资料得知土壤的相对介电常数为 12，电导率为 0.01，在深度为 1 m 的地下埋有水平放置的金属管线，半径为 20 cm ，管壁厚度 1 cm ，采用 gprmax 内置的理想金属 pec 材料。假设地质雷达天线频率为 500 MHz，时窗设为 60 ns ，源和接收点每次移动 1 cm，测线在研究区域窄边 1 m 处，垂直管线走向，网格大小为 1 cm。按照上述要求，我们编写的 in 文件如下：

#title: 3Dpec_100cm_20cm#material: 12.0 0.01 1.0 0.0 turang
#material: 1.0 0 1.0 0.0 air#domain: 4.000 2.100 2.000
#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 turang
#cylinder: 2 1.00 0 2 1.00 2.00 0.200 pec
#cylinder: 2 1.00 0 2 1.00 2.00 0.190 air#waveform: ricker 1.0 500e6 my_ricker
#hertzian_dipole: z 0.100 2.000 1.000 my_ricker#rx: 0.200 2.000 1.000
#src_steps: 0.100 0.000 0.000
#rx_steps: 0.100 0.000 0.000#geometry_view: 0 0 0 4.000 2.100 2.000 0.010 0.010 0.010 PEC n#python:
for i in range(1,61):print('#snapshot: 0 0 0 4.000 2.100 2.000 0.010 0.010 0.010 {} snapshot{}'.format((i/10)*10e-9,i))
#end_python:

2、管线倾斜的情况

#title: 3Dpec_100cm_20cm#material: 12.0 0.01 1.0 0.0 turang
#material: 1.0 0 1.0 0.0 air#domain: 4.000 2.100 2.000
#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 turang
#cylinder: 1.5 1.00 0 2.5 1.00 2.00 0.200 pec
#cylinder: 1.5 1.00 0 2.5 1.00 2.00 0.190 air#waveform: ricker 1.0 500e6 my_ricker
#hertzian_dipole: z 0.100 2.000 1.000 my_ricker#rx: 0.200 2.000 1.000
#src_steps: 0.100 0.000 0.000
#rx_steps: 0.100 0.000 0.000snapshot: 0 0 0 4.000 2.100 2.000 0.010 0.010 0.010 58e-9 snapshot1#geometry_view: 0 0 0 4.000 2.100 2.000 0.010 0.010 0.010 PEC n#python:
for i in range(1,61,5):print('#snapshot: 0 0 0 4.000 2.100 2.000 0.010 0.010 0.010 {} snapshot{}'.format((i/10)*10e-9,i))
#end_python:

对于倾斜的管线，只用设置 cylinder 参数即可。

#title: 3Dpec_100cm_20cm#material: 12.0 0.01 1.0 0.0 turang
#material: 1.0 0 1.0 0.0 air#domain: 4.000 2.100 2.000
#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 turang
#cylinder: 2 0.70 0 2 1.30 2.00 0.200 pec
#cylinder: 2 0.70 0 2 1.30 2.00 0.190 air#waveform: ricker 1.0 500e6 my_ricker
#hertzian_dipole: z 0.100 2.000 1.000 my_ricker#rx: 0.200 2.000 1.000
#src_steps: 0.100 0.000 0.000
#rx_steps: 0.100 0.000 0.000#snapshot: 0 0 0 4.000 2.100 2.000 0.010 0.010 0.010 58e-9 snapshot1#geometry_view: 0 0 0 4.000 2.100 2.000 0.010 0.010 0.010 PEC npython:
for i in range(1,61):print('#snapshot: 0 0 0 4.000 2.100 2.000 0.010 0.010 0.010 {} snapshot{}'.format((i/10)*10e-9,i))
end_python:

二、模型可视化

将 vti 文件导入 paraview 中可视化三维模型。

第一个 in 文件对应的模型

第二个 in 文件对应的模型

第三个 in 文件对应的模型

三、仿真结果

1、一条测线的雷达扫描图

2、不同时刻的波场图

3、波场图

[^1]本文仅供学习交流。

Gprmax 三维地质雷达建模及在 paraview 中的可视化相关推荐

3DMAX在三维GIS建模中的应用与优化
3DMAX在三维GIS建模中的应用与优化黄名华 (广西基础地理信息中心广西南宁 530023) [摘要] 将一个城市或某一区域的地理空间数据制作成三维的可视化地理信息系统,在当今有着非常重要的现 ...
中国有完全自主的三维几何建模引擎和几何约束求解器吗？
工业软件,就是"工业大脑".三维CAD系统,是工业核心数据的来源,属于核心的研发设计类工业软件.三维 CAD 软件的两大底层核心技术:三维几何建模引擎和几何约束求解器,目前均面临严 ...
Python三维地址建模教程【Gempy】
Gempy 是一个开源 Python 库,用于生成完整的 3D 结构地质模型.该库是从界面.断层和层方向创建地质模型的完整开发,它还关联地质层序列以表示岩石侵入和断层顺序. 地质建模算法基于通用协同克 ...
航飞原始影像外方位元素_【技术】无人机倾斜摄影建模技术在虚拟现实中的应用...
(如有侵权,请联系删除) 摘要针对于虚拟现实平台中构建三维场景的费时费力问题,基于无人机倾斜摄影建模技术构建三维模型,利用 3DS Max 建模软件进行模型优化,并结合 Unity 3D 引擎构建 ...
利用计算机绘制地质图的思路和方法,基于平面地质图的三维地质建模方法研究...
摘要: 平面地质图综合了地质野外勘察工作成果与地质专家知识,揭示了地区的岩石.地层和地质构造等信息,是人们了解区域地质最易获取和最直接的数据源.研究基于平面地质图的三维地质体建模方法,能够有效地解决缺 ...
三维地质建模数据处理
三维地质建模计算在地质工程.地球物理.矿产勘查等领域获得了广泛的应用,常用软件包括GOCAD.Surpac.XModel.DMine等.通过三维地质建模,既可以表达空间几何对象,也可以表现空间属性分布 ...
《流浪地球》地下城怎么建？三维地学建模探秘地下空间
"上天.入地.下海"是人类畅想科技发展的主题,众多科幻小说.电影都描摹了未来奇妙的地下世界.在电影<流浪地球>中,随着地球自转停止.太阳远离,地表环境恶化,人类只得居住 ...
GOCAD 2009 完整版功能强大的三维地质建模软件
为什么80%的码农都做不了架构师?>>> GOCAD 2009 完整版功能强大的三维地质建模软件 GOCAD(Geological Object Computer Aided ...
python实现3d扫描_三维激光扫描建模仪（基于树莓派）
MAKER:mfx2/译:趣无尽 Cherry(转载请注明出处) 只需要将物体放入到一个盒子里,就可以得到它的 3D 模型数据,是不是很方便呢? 今天我们就用树莓派制作一台三维激光扫描建模仪.能够将物 ...

Gprmax 三维地质雷达建模及在 paraview 中的可视化