Matlab实现的FEMIC的说明书

FEMIC程序是用来反演小回路频域电磁感应数据的。要启动代码，在Matlab命令窗口中输入start，然后点击“Enter”或“返回”按钮。然后会出现FEMIC的主界面，见图1。

它由几个输入区域组成，这几个区分别实现了：加载数据，反演过程控制和最终显示。

图1 主界面

下面对这些输入格式进行简介：

1 并行计算

用户可以指定反演使用多少个处理器。
如果用户提供的数字大于所使用计算机的容量，则将显示错误消息，指示该机器中允许的最大处理器数量如图2所示。

图2 并行计算输入错误的提示框

2 实测数据输入面板

输入数据面板由以下四块组成，如图3所示：

a 采集数据类型

指线圈配置，根据提供的数据输入文件自动分配。

b 数据输入文件

用户要求准备/上传扩展名为*.dat的数据输入文件，使用如图3所示的特定格式。数据文件应该只包含数值，由8列组成，排列如下：

• x-location：测深沿x方向的位置
• y-location：沿y的测深位置
• z-location：每个测深的仰角。当前程序版本假设是所有测深点的海拔均为零。
• Sounding_id：每个测深点的ID号。同一测深点的不同数据是相同的数字。
• Frequency：数据采集过程中使用的频率列表。
• STD误差：每次测量的标准偏差。
• Coil configuration：线圈配置，VDM输入1，HDM输入2，两者配置都有则输入3。
• 所有测深点测得的视电导率值列表

c 收发距偏移量

用户需上传*.dat文件，该文件由两列组成，如图3

• 数据采集期间的线圈间距。
• 给定线圈间距的相应频率值。

d 过滤/校准

用户需要上传电导率模型文件来校准实测电磁原始数据。该文件由以下四列组成：

• x-location：电导率值沿x方向的采样位置。
• y-location：电导率值沿y方向的采样位置。
• Depth：采样电导率值的深度。
• Conductivity：所有位置的采样电导率值列表。

要查看模型，点击“查看模型”按钮。为了运行校准和过滤处理，点击“运行过滤/校准”按钮。自动校准和过滤后的数据会替换原来的电磁原始数据。按下运行按钮后，将出现一条指示处理成功的消息。

图3 实测数据输入面板所需文件的格式

3.初始模型参数

这个模块允许用户利用任何可用的先验信息来输入参考模型。

用户可以简单地点击“创建分层模型”，为所有测深点创建一个均匀的起始模型。

初始模型对于所有测深点会自动为每个层/电导率设置约束权重为0。
约束权值是用于调整衡量反演模型到参考模型的偏离度的，取值范围为0 ~ 1。

约束权值为“0”表示没有约束，而如果值为“1”，则表示在反演执行过程中参考值将是硬约束。

用户若想在任意测深点处修改参考模型和约束权值，只需单击“添加图层”或插入图层，如图4所示。

图4 初始模型参数面板

 4. 模型/反演参数

在本面板，可直接输入反演参数，也可以加载*.dat文件来输入。该文件应该有一个包含以下参数的单列：

• Min. parameter：反演求解中使用的最小电导率值，标量
• Max. parameter：最大电导率值，标量。
• Bar:正障碍函数的参数，标量。
• Muh:水平正则化参数，标量或向量。
• Muv:垂直正则化参数，标量或向量。
• max iterations：最大反演迭代次数。
• err tol:变化误差的收敛准则。
• eca tol:对电导率变化的容限。
• Sx: x方向平滑度参数。
• Sy: y方向平滑度参数。
• Sz: z方向平滑度参数。

5 反演求解

本面板中，用户需选择：

• 反演问题的维数：一维、二维或三维反演。
• 可视化：用户可能会想添加灵敏度图，模型分辨率或反演得到的电导率模型内的最大DOI，通过选择相应的单选按钮。在显示常规反演结果后还会显示灵敏度和分辨率图。DOI估计是基于奥登伯格和李耀国教授（1999）的方法。

6 2D理论模型案例

我们模拟了一个电导率为0.02S/m的异常体，背景围岩均匀，电导率为0.01S/m（图5）。

图5 理论模型

 反演结果如图6所示。异常体恢复的很好，最大DOI估计匹配异常体的底部。

模型分辨率显示了模型域内由数据控制的区域(亮色区域)和反演约束(深蓝色区域)。这个案例的输入数据位于代码中的“synthetic example”文件夹中。

图6 a-测量数据；b-估计数据；c-反演模型；d-模型分辨率

 最后附上原始文档供参考。

原始文档资料下载