重磁最新边缘识别实验报告.docVIP

位场边缘识别与程序设计 实验 学生姓名： 专业名称： 学生学号： 指导老师： 一 、基本原理 2 1.1 位场边缘识别方法研究进展概述 2 1.2 数值计算类边缘识别方法的研究及计算 2 1.2.1 垂向导数（VDR） 2 1.2.2 总水平导数（THDR） 2 1.2.3 解析信号振幅（ASM） 2 1.2.4 倾斜角（TA） 3 1.2.5 二阶导数类边缘识别 3 二、 输入/输出数据格式设计 3 2.1 主要变量设计 3 2.2位场边缘识别程序数据格式设计 4 三、 总体设计 5 四、测试结果 5 4.1 测试用例 5 4.2 测试结果 6 4.3 结果分析 8 五、结论与建议 9 5.1 结论 9 5.2 建议 9 附录：边缘识别程序源代码 10 一 、基本原理 1.1 位场边缘识别方法研究进展概述 地质目标体边缘是指断裂构造线、不同地质体的边界线，实际上是具有一定密度或磁性差异的地质体的边界线．在地质体的边缘附近，重、磁异常变化率较大，故所有的边缘识别方法均利用这一特点进行设计。现有重、磁位场边缘识别方法分为数理统计、数值计算和其他三大类。 1.2 数值计算类边缘识别方法的研究及计算 数值类边缘识别方法均利用极大值位置或零值位置确定地质体的边缘位置，其依据的理论基础是二度体铅垂台阶模型的重力异常特征．在该模型边缘处重力异常总水平导数和解析信号振幅达到极大值、垂向导数达到零值．故可以利用这些特征位置来确定二度体铅垂台阶的边缘位置．确定倾斜二度体、不规则二度体及三度体边缘位置的理论均为二度体铅垂台阶模型理论的推广，但确定的边缘位置与真实位置有一定偏差．该偏差随着地质体边界形状、埋深、水平尺寸及物性差异等的变化而变化．因此，边缘识别结果是一种定性或半定量解释结果，与定量解释结果有一定区别，识别结果可作为边缘位置定量反演的初值。 1.2.1 垂向导数（VDR） 垂向导数方法利用零值位置确定地质体的边缘位置，重力异常可直接使用，对磁力异常必须转换成磁源重力异常（假重力异常）或化极磁力异常才可以使用。 垂向导数方法研究历史较早，方法较成熟，应用频率较高．通过我们的试验和研究认为：随着地质体埋深的增大，垂向导数所确定的边缘位置偏离实际位置的距离越大。 1.2.2 总水平导数（THDR） 总水平导数是利用其极大值位置来确定地质体的边缘位置，适用于重力异常，对磁力异常必须换算成磁源重力异常或化极磁力异常才可以使用。 1.2.3 解析信号振幅（ASM） 解析信号振幅也是利用极大值位置来确定地质体的边缘位置，适用于重力异常和磁力异常。 1.2.4 倾斜角（TA） 倾斜角实质上是垂向导数和总水平导数的比值。由于倾斜角为一阶导数的比值，所以能很好地平衡高幅值异常和低幅值异常，起到边缘增强的效果。 1.2.5 二阶导数类边缘识别 为了增强边缘分辨能力，和克服当总水平导数等于０时倾斜角存在“解析奇点”，会使得计算结果不稳定。 输入/输出数据格式设计 依据上述原理，现对上述各种边缘识别方法实现进行程序设计。 2.1 主要变量设计 cmd_file:参数文件名 input_field_filename： 输入位场文件名 output_field_filename：输出转换后位场文件名 expan_2D_method： 2D扩边方法选择 Field_Deriv_method：边缘识别方法选择 num_area：计算区域大小选择 mpoint,nline：点数，线数 m0,m1,m2,m3：扩边后X方向上各端点 n0,n1,n2,n3：扩边后Y方向上各端点 Xmin,Xmax：X坐标最小，最大值 ymin,Ymax：Y坐标最小，最大值 Ori_Field：原始场值，二维数组，m3*n3 Deriv_Field：转换后后场值，二维数组，m3*n3 Deriv_operator_x：x方向转换因子，二维数组，m3*n3 Deriv_operator_y：y方向转换因子，二维数组，m3*n3 Deriv_operator_z：z方向转换因子，二维数组，m3*n3 2.2位场边缘识别程序数据格式设计 cmd文件，输入参数文件，格式如下：（cmd.txt） 原始位场文件，输入文件，grd格式(ASCII) 延拓后位场文件，输出文件，grd格式(ASCII) Grd格式如下： 总体设计 四、测试结果 4.1 测试用例 （1）观测面上的重力异常存放在“anomaly.grd”中，坐标单位为m。 （2）形体边界存放在“rectangle.bln”中， 坐标单位 4.2 测试结果 4.3 结果分析 1) 从图4.2.1-图4.2.7可以看出，VDR