基于gprMax和matlab的复杂地下模型正演----以地下空洞为例

概要

gprMax 非均匀介质建模，如地下空洞、富水等等

整体架构流程

1.利用python或matlab生成h5文件
2.导入gprMax的in文件中
3.为h5文件的格个网格填充介电常数
4.gprMax正演开始
5.matlab或python处理out文件（一般只做增益即可）

技术细节

以探测公路空洞为例
matlab生成h5文件，代码如下

close all
clear
clc

% 此程序是写gprmax随机介质建模的几何文件
% 创建任意随机介质
% 创建地下层状多相随机介质

% 给定网格单元边长
dx = 0.002;
dy = 0.002;
dz = 0.002;

% 网格单元个数，注意x,y,x的方向
nex = 1;
ney = 1000;
nez = 2000;

% 设置保存的文件名
file_name = 'kongdong.h5';
if exist(file_name, 'file') == 2
    % 如果文件存在，删除文件
    delete(file_name);
end
% 创建一个HDF5文件
h5create(file_name,'/data',[nex ney nez]);
dx_dy_dz = [dx dy dz];

% 创建一个充满随机数的大矩阵
data = zeros(nex,ney,nez);

% 第一层填充0
% 第一层默认为0
data(:,1