学位论文-实验四位场边缘识别程序设计实验.doc

《重磁资料处理与解释》实验四 位场边缘识别程序设计实验 专业名称：地球物理学 学生姓名： 学生学号： 指导老师：王万银、纪新林、纪晓琳、邱世灿 提交日期：2016-1-3 1、基本原理 地质目标体边缘时指断裂构造线、不同地质体的边界线，实际上是具有一定密度或磁性差异的地质体的边界线，在地质体的边缘附近，重、磁异常变化率较大，故所有的边缘识别方法均利用这一特点进行设计。现在有重、磁位场边缘识别方法分为数理统计、数值计算以及其他三大类。 数值类边缘识别方法均利用极大值位置或零值位置确定地质体的边缘位置，其依据的理论基础是二度体铅垂台阶模型的重力异常特征。在该模型边缘处重力异常总水平导数和解析信号振幅达到极大值、垂向导数达到零值。故可以利用这些特征位置来确定二度体铅垂台阶的边缘位置，确定倾斜二度体、不规则二度体及三度体边缘位置的理论均为二度体铅垂台阶模型理论的推广，但确定的边缘位置和真实的位置有一定的偏差。该偏差随着地质体边界形状、埋深、水平尺寸及物性差异等的变化而变化。因此，边缘识别结果是一种定性或半定量解释结果，与定量解释结果有一定区别，识别结果可作为边缘位置定量反演的初值。 (1)垂向导数：垂向导数方法利用零值位置确定地质体的边缘位置，重力异常可以直接使用，对磁力异常必须转化为磁源重力异常或化极磁力异常才可以使用 ( 1.1) （2）解析信号振幅：解析信号振幅也是利用极大值位置来确定地质体的边缘位置适用于重、磁力异常 (1.2) （3）总水平导数（THDR） (1.3) 2、输入/输出数据格式设计 依据上述原理，现在对上述各种边缘识别方法进行程序设计。 2.1 输入数据格式设计 本次实验给了正演的重力异常数据，为.GRD格式，均为实型变量。 例如： DSAA 201 201 -1000.000000 1000.000000 -1000.000000 1000.000000 5.549671E-01 23.539846 5.549671E-01 5.634658E-01 5.721339E-01 5.808522E-01 5.897312E-01 5.987253E-01 6.078691E-01 6.171604E-01 … 2.2 输出数据格式设计 计算结果输出数据格式与输入格式对应，格式为.GRD格式，均为实型变量。 例如： DSAA 201 201 -1000.000 1000.000 -1000.000 1000.000 -0.1465084 0.3190881 -3.6523044E-02 -3.3485338E-02 -3.3061244E-02 -3.2748722E-02 -3.2688729E-02 -3.2654848E-02 -3.2723978E-02 -3.2787599E-02 -3.2927759E-02 -3.3138681E-02 … 2.3 参数文件数据格式设计 将以上部分量保存在一个文件中，该文件名变量为cmdfile，字符串变量，长度不超过80，全路径名。在该文件中保存的参数如下： 输入数据文件名input_file，字符串变量，长度不超过80； 输出vdr数据文件名output_file_vdr，字符串变量，长度不超过80； 输出thdr数据文件名output_file_thdr，字符串变量，长度不超过80; 输出asm数据文件名output_file_asm，字符串变量，长度不超过80 factor_m：扩边比例因子，实型变量（1）； 总体设计 此次程序采用IPO结构设计，首先通过读取cmd文件，得到相关输入参数：输入数据文件名gravity.grd、输出vdr文件名field_vrd.grd、输出thdr文件名field_thdr.grd、输出asm文件名field_asm.grd、扩边比例因子factor_m；然后确定确定扩边网格的大小，扩边数据点号位置；再从观测面位场数据文件中读取数据。下一步，进行二维余弦扩边，将扩完边的数据进行快速二维傅里叶变换，转换到频率域；接下来在频率域求出在x,y,z方向的导数并反变换；最后求出VDR、THDR、ASM数据。最后去除扩边部分后输出。总体设计见表1。 输入参数：输入数据文件名gravity.grd、输出vdr文件名gravity_v