《GPRMAX2D知识点汇总.doc》由会员分享，可在线阅读，更多相关《GPRMAX2D知识点汇总.doc（9页珍藏版）》请在金锄头文库上搜索。
1、GprMax是爱丁堡大学的Antonis Giannopoulos于1996年推出来的一种基于时域有限差分（FDTD）算法和理想匹配层（PML）边界吸收条件的探地雷达正演数值模拟软件，用于探地雷达成像研究。其中，GprMax2D是二维正演，GprMax3D为三维正演。该软件可以在Windows、Linux和MacOS三个平台上使用。本文主要针对Windows平台进行说明。一、软件获得该软件为免费软件，可以去GprMax官网下载。也可点此直接下载。二、准备工作软件无需安装，下载后用解压工具解压，找到Windows文件夹，直接双击GprMax2D.exe即可运行。但是为了方便批量模拟，建议把Win
2、dows文件夹下的文件（cygwin1.dll、GprMax2D.exe和GprMax3D.exe）全部复制到系统盘系统搜索路径下，例如：C:Windows。如果只进行二维正演，只复制cygwin1.dll和GprMax2D.exe即可。当然也可放在任意路径下，只要使用时包含所在路径就OK了。笔者习惯放在系统要目录下，即C:。另外，最好建立一个输入文件和模拟结果存放的专用文件夹，并且把tools文件夹下的文件（gprmax.m、gprmax2g.m、gprmax3g.m、gprmaxde.m和gprmaxso.m）全部复制到该专用文件夹下，注意文件夹名最好使用英文。如果只进行二维正演，只复制g
3、prmax.m和gprmax2g.m即可。笔者习惯使用D:GPR。正演结果需要用MATLAB进行绘图，因此需要安装有MATLAB软件。软件下载地址和安装方法此处不再说明，以后使用将假设读者已经成功安装MATLAB软件。三、软件使用一般进行数值模拟时通常都是若干个对比模型进行模拟，因此本文只介绍批量模拟的方法。如此一来，笔者将认为读者的软件和笔者一样放在C:，输入文件和gprmax.m、gprmax2g.m文件已经放在D:GPR。欲进行正演模拟需要先建立输入文件，输入文件的建立将在后面介绍，此处假设在D:GPR已经有若干输入文件，文件名分别为fname1.in，fname2.in，fname3.
4、in。打开记事本，写入以下两行直线间的内容，另存为*.bat文件，例如：GPR.bat。其中括号里的内容为本条语句的说明。-path C:;%path% (把软件所在路径包含进系统搜索路径)cd/d D:GPR (设置专用文件夹路径为活动路径)gprmax2d D:GPRfname1.in （按顺序进行批量正演模拟）gprmax2d D:GPRfname2.ingprmax2d D:GPRfname3.in-完成后保存并关闭，然后双击运行*.bat文件即可进行批量正演模拟四、建立输入文件输入文件是纯文本文件，可以用任何文字处理程序编辑，而且内容全部为英文。打开记事本，写入输入文件内容，另存为*
5、.in即可。输入文件必需包括模拟模型的所有必要信息，每一个有效命令行都必须以符号（#）开始，否则该命令行无效。所有命令的一般语法为：#命令名: 参数1 参数2 参数3 为了方便介绍GprMax2D的命令及参数，做了以下约定：f 代表实数，如1.5i 代表整数，如15c 代表单个字母，如ystr 代表字条串，如airfile 代表文件名，如test.in所有的长度单位均为m所有的时间单位均为s所有的频率单位均为Hz模型空间的坐标原点（0,0）在左下角输入文件中除了#analysis:和#end_analysis:命令外，先后顺序一般是没有要求的。GprMax2D 2.0版一共有32个命令：-#t
6、itle: str输出文件的标题，包含在输出文件中。#domain: f1 f2指定模型空间大小，单位是米。f1 x方向的大小f2 y方向的大小#dx_dy: f1 f2指定网格步长，单位是米。f1 x方向的步长，即xf2 y方向的步长，即y模型划分的网格数目是由#domain:和#dx_dy:两个命令共同决定的。网格划分越小模型越精确，但模拟的计算量也越在，需要权衡确定。#time_step_stability_factor: f1时间步长的稳定系数，取值范围（0，1，实际使用时间为f1t。#time_window: f1 or i1时窗大小，单位是秒。f1 为时窗，如f1=16e-9，时窗
7、为16nsi1 为时间步长数，时窗为i1t。#messages: c1是否开启屏幕信息。y(es)开启n(o)关闭，默认n。#number_of_media: i1介质数目，i1应大于10。系统分配了10种介质（自由空间即空气free_space，金属pec和8种用户自定义介质）的内存空间，当模型中使用介质数量超过10种时需要使用该命令。#nips_number: i1指定分配的空间。存储重要信息模型的数组大小可以在内部计算，但通常比实际需要的大，可以使用较小的数字节省空间。如果分配过小，系统会发出错误并提醒分配更多空间#media_file: file1要包含的介质文件的文件名（如果有必要还
8、包括文件路径）。介质文件格式为：# f1 f2 f3 f4 f5 f6 str1参数的具体含义参考#medium:命令。#geometry_file: file1存储模型文件（*.geo）。#medium: f1 f2 f3 f4 f5 f6 str1在输入文件中定义介质属性。f1 静态时介质的相对介电常数f2 在理论上的无限频率的相对介电常数f3 介质的弛豫时间（秒）f4 静态时介质的电导率（西门子/米）f5 相对磁导率的介质f6 介质的磁导率str1 介质的名称如果不想使用debye介质，可设置f3为0.0，此时系统将只使用f1和f4来描述介质介电性能。如果介质为非磁性，则设置f5为1.0
9、，f6为0.0。如果使用f3的值，该值应始终高于模型中使用的时间步长t。-该组命令不建议使用，除非是熟悉ABC的人用来定制和优化吸收边界条件。#abc_type: pml恢复默认的PML边界。#abc_order: i1ABC次序。取值范围1，2，3，默认为3。次序越低，性能越差。#abc_stability_factors: f1 f2 f3指定ABC的稳定系数。#abc_optimization_angles: f1 f2 f3优化ABC为特定的入射角。#abc_mixing_parameters: f1 f2 f3 f4 f5 f6改变使用ABC的时间和偏导数的离散化。默认使用box。#
10、pml_layers: i1指定PML使用的网格数量。默认为8。数量越多性能越佳。-#box: f1 f2 f3 f4 str1使用定义好的介质定义一个矩形模型。f1 矩形左下角x坐标，单位米f2 矩形左下角y坐标，单位米f3 矩形右上角x坐标，单位米f4 矩形右上角y坐标，单位米str1 使用的介质的名称#cylinder: f1 f2 f3 str1使用定义好的介质定义一个圆形模型。f1 圆心x坐标，单位米f2 圆心y坐标，单位米f3 圆半径R，单位米str1 使用的介质的名称#x_segment: f1 f2 f3 f4 f5 str1使用定义好的介质定义一个圆沿x轴裁剪的模型。f1 圆
11、心x坐标，单位米f2 圆心y坐标，单位米f3 沿x轴开始的位置f4 沿x轴结束的位置f5 圆半径R，单位米str1 使用的介质的名称#y_segment: f1 f2 f3 f4 f5 str1使用定义好的介质定义一个圆沿y轴裁剪的模型。f1 圆心x坐标，单位米f2 圆心y坐标，单位米f3 沿y轴开始的位置f4 沿y轴结束的位置f5 圆半径R，单位米str1 使用的介质的名称#triangle: f1 f2 f3 f4 f5 f6 str1f1 第一顶点x坐标，单位米f2 第一顶点y坐标，单位米f3 第二顶点x坐标，单位米f4 第二顶点y坐标，单位米f5 第三顶点x坐标，单位米f6 第三顶点y
12、坐标，单位米str1 使用的介质的名称模型的建立是有先后顺序的，后建立的模型会把先建立的模型覆盖住。-该组命令是数值模拟的核心。一个输入文件中可以出现多组数值模拟命令，但是数值模拟的控制命令必须放在#analysis:和#end_analysis:之间。#analysis: i1 file1 c1数值模拟开始的命令。i1 总步数file1 模拟结果输出文件名（*.out）c1 类型，取值范围a（文本文件）或b（二进制文件）#end_analysis:数值模拟结束的命令，无参数。#tx: f1 f2 str1 f3 f4发射天线的属性。f1 发射天线位置的x坐标，单位米f2 发射天线位置的y坐标，单位米str1 定义过的发射天线的名称f3 发射天线的延迟时间f4 发射天线的清除时间#rx: f1 f2接收天线的属性。f1 接收天线位置的x坐标，单位米f2 接收天线位置的y坐标，单位米#rx_box: f1 f2 f3 f4 f5 f6接收区域的属性。f1 接收区域左下角的x坐标，单位米f2 接收区域左下角的y坐标，单位米f3 接收区域右上角的x坐标，单位米f4 接收区域右上角的y坐标，单位米
