RGIS2010数据处理反演-RGIS2006重力数据整理

2010年05月05日数据预处理数据处理反演rgis2006区域重力信息系统RGIS2006培训提纲RGIS软件应用流程1RGIS重磁数据预处理2RGIS重磁处理与反演334RGIS2006实测重力资料整理RGIS2010-应用流程RGIS软件应用流程数据预处理频率域平面数据处理空间域平面数据处理剖面数据处理数据处理反演培训提纲RGIS软件应用流程1RGIS重磁数据预处理2RGIS重磁处理与反演334RGIS2006实测重力资料整理RGIS2010-平面数据处理RGIS实测重力资料整理数据预处理RGIS2010-平面数据处理RGIS实测重力资料整理数据预处理——网格化文件数据要求被网格化的数据文件是ASCII明码格式后缀：“*.txt、*.dat、*.xls”文件，至少包含三列数据。为平面数据处理做准备RGIS2010-平面数据处理RGIS实测重力资料整理数据预处理——样条光滑输入数据格式为网格化后的数据为绘制等值线图而准备RGIS2010-平面数据处理RGIS实测重力资料整理数据预处理——网格数学计算输入数据格式为网格化后的数据两个网格化数据必须行列数、行列距一致RGIS2010-平面数据处理RGIS实测重力资料整理输入数据格式为网格化后的数据后缀：“*.txt、*.dat、*.xls”文件数据预处理——数据列操作RGIS2010-平面数据处理RGIS实测重力资料整理数据预处理——空白区填补输入数据格式为网格化后的数据（*.grd）为空间域平面数据处理准备RGIS2010-平面数据处理RGIS实测重力资料整理数据预处理——还原空白区输入数据格式为网格化数据（*.grd）第一个为填补了空区的数据第二个为原始有空区的数据培训提纲RGIS软件应用流程1RGIS重磁数据预处理2RGIS重磁处理与反演334RGIS2006实测重力资料整理RGIS2010-平面数据处理RGIS2010-平面数据处理为了获得所探测目标地质体（也称异常体或场源）的信息，获得可靠的地质解释，通常需要进行相关的数据处理，以突出目标体的异常特征。运行平面数据处理模块RGIS2010-平面数据处理平面数据处理模块RGIS2010-平面数据处理RGIS2010-平面数据处理RGIS2010-平面数据处理用于平面数据处理的数据格式为网格化数据（*.grd)RGIS2010-平面数据处理-延拓RGIS2010-平面数据处理——延拓向上延拓可压制浅而小的地质体的场，突出大地质体的场（包括浅部和深部的大地质体）。向下延拓有分解的作用（向下延拓要慎用），向上延拓则有合并场的作用。要注意延拓高度与反映的深度没有简单的对应关系。合适的延拓高度要通过试验确定。有时，想突出特点不同的几类信息，需要几个高度的延拓资料，分别突出不同的信息，否则没有必要延拓数个高度。参数：数据扩边的设置和延拓参数设置RGIS2010-平面数据处理—延拓RGIS2010-平面数据处理——延拓（频率域）实际扩边最好大于最小扩边数的2倍按照2的整数幂进行扩边延拓高度的设置延拓高度的单位和数据点距单位相同RGIS2010-平面数据处理——延拓（空间域）以网格数据的网格距为单位的RGIS2010-平面数据处理——导数水平导数（方向导数）主要用于突出走向垂直于求导方向的断裂、脉岩的位置和宽大地质体的边界线求导方向应依据实际和推断的构造方向来确定。当存在几组构造线时，应垂直每组构造线分别求导，不应不加分析地一律求取0°、45°、90°、135°4个方向的导数水平一导极大值近似的确定垂向断裂位置垂向一导突出局部异常。利用垂向一导、二导零值线划分磁性体边界需注意。RGIS2010-平面数据处理——导数（频率域）如果若只求垂向导数，则要置求导方向偏角为0只求水平方向上的导数时，需要置求导方向倾角为0参数：参数设置RGIS2010-平面数据处理——导数（空间域）计算跨度是指参与导数计算的两个网格点之间的距离，以点距数为单位，为偶数。计算跨度越大，压制高频成分越强，根据具体情况确定RGIS2010-平面数据处理——导数（空间域）RGIS2010-平面数据处理——水平总梯度模、解析信号水平总梯度模（水平方向导数模）、解析信号（总梯度模）是可以用来确定线性构造如断裂或岩性接触面的位置。将起伏地形上的观测资料换算为水平地形条件下的观测资料，然后反演。RGIS2010-平面数据处理——曲化平这里四个方向的扩边数是指计算数据进行曲化平前的扩边数。扩边大小不小于30个点距。四个方向的扩边数必须相等。曲化平平面高度值不能小于地形高程数据的最大值。即曲化平平面高度必须位于观测曲面所有测点之上RGIS2010-平面数据处理——正则化滤波水平几何尺度是指被滤掉异常的最大宽度。因子越大，被滤掉的微小异常越多，为了确定最好的滤波程度，可以试用不同的参数进行对比，从而确定最合适的几何尺度因子。补偿因子越小，滤波作用越强，压制高频成分效果越好。指数因子越大，补偿圆滑滤波器越接近理想滤波器，滤波作用越强，压制高频成分效果越好。RGIS2010-平面数据处理——补偿圆滑滤波把圆滑点设为中心点，求窗口内其它点异常的平均值，并付给圆滑点设置滤波窗口大小，窗口越大压制高频成分的滤波效果越好。此处所给定的滤波窗口的行列数必需是奇数。RGIS2010-平面数据处理——滑动平均把圆滑点设为中心点，求窗口内其它点异常的平均值，并付给圆滑点该对话框显示了趋势分析中趋势面的阶数、系数个数和系数值、趋势值统计参数和偏差统计参数。对于异常趋势分析效果的好坏，首先取决于全区区域异常能否用一个多项式来描述；其次取决于多项式函数的阶次。当所取阶次过高时，则区域异常中必然含有较多的局部异常成分。RGIS2010-平面数据处理——趋势分析用来凸现异常中可能存在的线性构造特征，压制干扰因素产生的异常，使得测区的线性异常信息更加突出。RGIS2010-平面数据处理——线性增强RGIS2010-平面数据处理——线性回归在一定面积内，用高程与布格重力异常或自由空间异常进行回归分析，可以大致了解地壳均衡补偿的程度将布格重力异常和地形两个的物理量进行相关分析，用来反映测区内地形起伏对重力异常的影响情况。布格重力异常数据地形数据RGIS2010-平面数据处理——相关分析在斜磁化影响严重时，才需进行化极处理大面积化极需用变磁倾角化极方法。化极只能使磁异常范围与磁性体的关系更为密切。但不能依据化极资料精确确定磁性体边界。RGIS2010-平面数据处理——化极包括：磁异常化极、低纬度化极、变纬度化极按高斯坐标网格化输入默认值90度N行方位角为减90度实际工区的倾角、偏角RGIS2010-平面数据处理——化极自定义坐标时，输入网格数据行与北的夹角行方位角减90度实际工区的倾角、偏角N行方位角N行方位角RGIS2010-平面数据处理——磁异常化极在赤道附近的低磁纬度地区，观测所得的磁场形态要远比在高磁纬度地区相同地质条件下复杂得多。RGIS2010-平面数据处理——低纬度化极当测区纬度跨度比较大，一般来说达到2度时常规化极已不能满足较高的精度要求，此时需要作变纬度化极处理。RGIS2010-平面数据处理——变纬度化极RGIS2010-平面数据处理——磁源重力异常磁异常可转换为重力异常，当测区同时具有重、磁资料时，为研究重、磁异常的相关性，识别同源异常和基本同源异常，可进行该项换算。有时为避开斜磁化条件下的反演，也进行该项换算。RGIS2010-平面数据处理——磁源重力异常假设剩余密度RGIS2010-剖面数据处理运行剖面数据处理模块RGIS2010-剖面数据处理测点点距相同剖面测点数必须为奇数数据为三列，列之间为空格00.000-208.7821000.000-208.816……20000.000-210.524第一列为剖面距离，第二列为测点高程，第三列为测点重/磁异常保存格式为ASCII格式,后缀为“*.dat”或“*.txt”数据格式RGIS2010-重磁异常反演——频率域解析延拓延拓高度：向上延拓输入负值，向下延拓输入正值。磁重转换磁化倾角：磁性体被磁化的方向。泊松比值：当泊松比值为1时得同源伪重力异常，当泊松比值取实际值时，得同源剩余密度的重力异常。RGIS2010-重磁异常反演——频率域水平分离场异常距离：欲分离的两个异常之间的距离，单位为点距。异常幅值比：是欲分离的两个异常的幅度比。先对原始数据求水平导数，根据极值确定异常之间的距离。异常幅值比大于1时分离效果好。水平分离功能只适用于两个形态相似的异常的分离。RGIS2010-重磁异常反演——频率域水平分离场延拓高度：向上延拓输入负值，向下延拓输入正值。RGIS2010-重磁异常反演——频率域提取区域/局部异常操作步骤：先进行对数功率谱计算，根据对数功率谱曲线，找出背景场源与局部场源的两个频谱段，利用鼠标左键在对数功率谱曲线图分别确定能代表这两个频谱段的四个点，两个频谱段的分析及四个点的确定是关键，应注意其代表性。选取点时一定要注意确保P1P2P3P4。第1、2点对应于低频段（深源或异常体）第3、4点对应于高频段（干扰或浅部）RGIS2010-重磁异常反演——频率域原始数据计算功率谱选点提取区域场提取局部场RGIS2010-重磁异常反演——频率域去高频滤波一般可取1.5，实际应用时最好要经过多次试验，以使其有针对性。RGIS2010-重磁异常反演——频率域RGIS2010-剖面数据处理——空间域培训提纲RGIS软件应用流程1RGIS实测重力资料整理2RGIS重磁数据预处理334地球物理（重磁电）三维数据可视化RGIS重磁处理与反演RGIS2010-重磁异常反演运行重磁反演解释模块RGIS2010-重磁异常反演剖面磁源深度计算用于通过剖面异常粗略计算，对于孤立异常，结果基本可靠应用方法：异常场响应的低频段快速衰减，干扰场的频率响应下降较缓；整条剖面分析或其中一段进行频谱分析；整条剖面分析，直接点击绘制功率谱按钮，计算整条剖面；对其中一段进行功率谱计算，需要选择剖面段：•第一个点：shift+鼠标左键选取•第二个点：shift+鼠标右键选取•第一点和第二点之间的点数应大于8个点RGIS2010-重磁异常反演剖面磁源深度计算RGIS2010-重磁异常反演剖面磁源深度计算低频数：从功率谱曲线上读出的源场响应的低频段的起始频数，不能为零，为大于等于1的整数。低频数：从功率谱曲线上读出的源场响应的低频段的终止频数，大于低频数的整数。深度值：计算所得场源体的顶面埋深，单位和剖面距离单位相同。RGIS2010-重磁异常反演平面磁源深度计算读入数据位网格数据异常形态简单最好是孤立异常体应用鼠标在功率谱曲线上选择低频段斜率相近的频段。RGIS2010-重磁异常反演3D磁源深度计算参数设置：窗口宽度：2的整数幂，16，32滑动距离：连续计算点之间的距离，单位是点距低频数：异常体功率谱起始频数高频数：异常体功率谱终止频数低频数和高频数需要多次试验选择合理的结果。计算结果：所得结果反应磁性层或磁性体顶面埋深，深度单位同读入数据坐标的单位相同，km或m。RGIS2010-重磁异常反演3D磁源深度计算均匀磁化磁异常曲线能够全面反映异常计算结果反应磁性体顶面平均深度磁性体形态较为简单计算结果的单位和计算数据横纵坐标的单位相同RGIS2010-重磁异常反演密度界面反演根据重力异常来计算地下密度界面的起伏情况。单位为：km数据网格化时，单位选择kmRGIS2010-重磁异常反演磁性界面反演根据磁异常来计算地下磁性界面的起伏情况。单位为：km数据网格化时，单位选择kmRGIS2010-重磁异常反演磁性界面反演根据磁异常来计算地下磁性界面的起伏情况。单位为：km数据网格化时，单位选择kmRGIS2010-重磁异常反演2.5D重磁剖面反演提供了成熟的二度半棱柱体模型重磁联合反演数据格式为明码数据（三列）：X，h，g/dtRGIS2010-重磁异常反演2.5D重磁剖面反演正演：不读入实测剖面，通过直接建立模型正演地质体重磁异常场反