一种二维网格数据的补偿方向滤波方法技术

本发明专利技术涉及地球物理勘探的重力、磁力勘探数据处理技术，是一种二维网格数据的补偿方向滤波方法，将原始数据进行网格化处理，在窗口内划分方位角和每个方位对应的窗口内的节点分布，设置每个方位上背景场滤波时节点系数和剩余场补偿滤波求取时的节点系数，计算窗口内滤波方位的背景场值，循环迭代完成补偿方向滤波。本发明专利技术滤波的方向性强，满足了重磁处理提取方向异常信息的需要，同时也能满足压制某方向异常信息的需要。本发明专利技术推动了重磁勘探数据处理解释技术的进步，同时也为其它学科的二维网格数据的方向滤波处理提供了新方法。

【技术实现步骤摘要】
【专利摘要】本专利技术涉及地球物理勘探的重力、磁力勘探数据处理技术，是，将原始数据进行网格化处理，在窗口内划分方位角和每个方位对应的窗口内的节点分布，设置每个方位上背景场滤波时节点系数和剩余场补偿滤波求取时的节点系数，计算窗口内滤波方位的背景场值，循环迭代完成补偿方向滤波。本专利技术滤波的方向性强，满足了重磁处理提取方向异常信息的需要，同时也能满足压制某方向异常信息的需要。本专利技术推动了重磁勘探数据处理解释技术的进步，同时也为其它学科的二维网格数据的方向滤波处理提供了新方法。【专利说明】
本专利技术涉及地球物理勘探的重力、磁力勘探数据处理领域，属于数字滤波技术，是。
技术介绍
在磁力(包括航磁)资料采集中，由于测线间的一致性差异等原因引起的沿线的线性异常是具有方向的，这些方向性异常是地质研究的干扰(即数据噪音)，因而，在重磁勘探资料处理解释中需要提取具有方向性的异常信息，以满足解释断层、构造的需要，需要在资料处理解释前进行处理，即利用方向滤波方法压制重磁数据中方向性的重磁噪音数据，使得滤波处理后的重磁数据更好地反映地质信息是需要解决的问题。目前地球物理勘探的重力、磁力勘探数据处理中常用的频率域椭圆方向滤波是在频率域内实现，它的作用是提取出沿选择方向的方向性异常。这种方法常见的步骤是:1.将重力(或磁力)数据网格化，2.对重力网格数据进行傅立叶变换，3.选择滤波方向，选取椭圆方向滤波器，4.对傅立叶变换后的数据进行椭圆方向滤波，获得滤波后的频谱数据，5.对滤波后的频谱数据进行反傅立叶变换，即获得椭圆方向滤波后的重力(或磁力)异常网格数据。这种方法不具有迭代性和局部场补偿作用，而且方向性差。地质出版社，2005.6，P194_204公开了一种《重力场与重力勘探》方法，利用空间域滑动窗口加权平均法建立一种圆滑滤波器，对数据具有圆滑滤波作用，可获得区域场，但是这种方法的方向滤波方法性能差，方向性异常往往被圆滑，无法提取重磁勘探资料处理解释中具有方向性的异常信息，不能利用方向滤波方法压制重磁数据中方向性的重磁噪音数据。
技术实现思路
本专利技术目的在于为重磁数据处理提供一种具有较好方向滤波功能，可提取方向性重磁异常数据或压制具有方向性的噪音数据的二维网格数据的补偿方向滤波方法。本专利技术通过以下步骤实现:I)进行野外重力或磁力勘探，采集得到原始数据并进行网格化处理；2)运算窗口采用7行7列窗口，迭代次数取I?10 ;步骤2)所述的运算窗口大小为7行7列。3)在窗口内划分方位角和每个方位对应的窗口内的节点分布；步骤3)所述的窗口内划分方位角是划分为南和北方位、东和西方位、北东和西南方位、北西和东南方位、北北东和南南西方位、北东东和南西西方位、北北西和南南东方位、北西西和南东东方位16个。步骤3)所述的窗口内每对方位的节点分布相同，按以下分布:南和北这对方位的节点分布是一组，窗口节点的行序为从下到上，窗口节点的列的序为从左到右，坐标组行序，列序为:南和北方位的窗口内节点分布为(I，4)、(2,4)、(3,4)、(4，4)、(5，4)、(6，4)、(7，4)；东和西方位的节点分布是一组，东和西方位的窗口内节点分布为(4，I)、(4,2),(4，3)、(4，4)、(4，5)、(4，6)、(4，7)；北东和南西方位的节点分布是一组，北东和南西方位的窗口内节点分布为(1，I)、(2，2)、(3，3)、(4，4)、(5，5)、(6，6)、(7，7)；北西和南东方位的节点分布是一组，北西和南东方位的窗口内节点分布为(7，I )、(6，2)、(5，3)、(4，4)、(3，5)、(2，6)、(1，7)；北北东和南南西方位的节点分布是一组(图1的5)，北北东和南南西方位的窗口内节点分布为(1，3)、(2，3)、(3，3)、(3，4)、(4，4)、(5，4)、(5，5)、(6，5)、(7，5)；北东东和南西西方位的节点分布是一组(图1的6)，北东东和南西西方位的窗口内节点分布为(3，I)、(3，2)、(3，3)、(4，3)、(4，4)、(4，5)、(5，5)、(5，6)、(5，7)；北北西和南南东方位的节点分布是一组，北北西和南南东方位的窗口内节点分布为(1，5)、(2，5)、(3，5)、(3，4)、(4，4)、(5，4)、(5，3)、(6，3)、(7，3)；北西西和南东东方位的节点分布是一组，北西西和南东东方位的窗口内节点分布为(5，I)、(5，2)、(5，3)、(4，3)、(4，4)、(4，5)、(3，5)、(3，6)、(3，7)。4)设置每个方位上背景场滤波时节点系数和剩余场补偿滤波求取时的节点系数；步骤4)所述的节点系数为:东西向滤波时，背景场滤波时节点系数和剩余场补偿滤波求取时的节点系数为:(4，1)、(4，2)、(4，3)、(4，4)、(4，5)、(4，6)、(4，7)节点的背景场滤波时节点系数自左向右分别为 3、5、7、12、7、5、3 ;剩余场偿滤波求取时的(1，4)、(2，4)、(3，4)、(4，4)、(5，4)、(6，4)、(7，4)节点系数分别为1、1、1、1、1、1、1 ；南北向滤波时，背景场滤波(1，4)、(2，4)、(3，4)、(4，4)、(5，4)、(6，4)、(7，4)节点系数分别为 3、5、7、12、7、5、3，剩余场偿滤波(4，1)、(4，2)、(4，3)、(4，4)、(4，5)、(4，6)、(4，7)节点的系数分别为1、1、1、1、1、1、1 ；其它各个方位的背景场滤波时节点系数和剩余场补偿滤波求取时的节点系数设置采用东西向滤波和南北向滤波设置方法类推。5)按照以下公式计算窗口内滤波方位的背景场值gw:gbJ = (I/ Σ (k (i, j))) X (Σ (k (i, j) Xg(i, j)))式中:k (i, j)为节点(i，j)的系数值，g(l, j)为节点(i，j)的重力或磁力值；步骤5)所述的滤波方位的背景场值gw，当东西向滤波时，先计算计算点(4，4)的背景场值gw (4，4):gbJ (4，4)= (1/42)X (3Xg(4, I)+5Xg(4，2)+7Xg(4，3)+12Xg(4，4)+7Xg(4，5) +5Xg(4, 6)+3Xg(4, 7))；式中g(4，I)、g(4，2)、g(4，3)、g(4，4)、g(4，5)、g(4，6)、g(4，7)分别表示计算窗 P内节点(4，I)、(4，2)、(4，3)、(4，4)、(4，5)、(4，6)、(4，7)的重力值或者是磁力值；依次移动运算窗口，计算工区平面网格的每一个计算点的背景场值；其它各方位滤波，替换所选坐标节点和对应节点的系数、按各方位的系数加权平均计算。6)按照以下公式计算窗口内滤波方位的剩余场补偿滤波值gbc;h:gbch= (1/ Σ (kbch (i, j))) X (Σ (kbch (i, j) Xgbch(i, j)))式中:kbch (i, j)为节点(i，j)的系数值kbc;h，g(i, j)为节点(i，j)的重力或磁力值；步骤6)所述的剩余场补偿滤波值gbc;h,东西向滤波时:gbch (4,4) =gbJ (4,4)-( (1/7) X (gbJ(1, 4) +gbJ (2，4) +gbJ (3，4) +g本文档来自技高网...

【技术保护点】
一种二维网格数据的补偿方向滤波方法，特点是通过以下步骤实现：1）进行野外重力或磁力勘探，采集得到原始数据并进行网格化处理；2）运算窗口采用7行7列窗口，迭代次数取1～10；3）在窗口内划分方位角和每个方位对应的窗口内的节点分布；4）设置每个方位上背景场滤波时节点系数和剩余场补偿滤波求取时的节点系数；5）按照以下公式计算窗口内滤波方位的背景场值gbj：gbj＝（1/∑（k（i,j）））×（∑（k（i,j）×g(i,j)））式中：k（i,j）为节点（i,j）的系数值，g(i,j)为节点（i，j）的重力或磁力值；6）按照以下公式计算窗口内滤波方位的剩余场补偿滤波值gbch：gbch=（1/∑（kbch（i,j）））×（∑（kbch（i,j）×gbch(i,j)））式中：kbch（i,j）为节点（i,j）的系数值kbch，g(i,j)为节点（i，j）的重力或磁力值；7）循环迭代直到完成迭代运算次数，完成补偿方向滤波。

【技术特征摘要】

【专利技术属性】
技术研发人员：刘云祥，
申请(专利权)人：中国石油天然气集团公司，中国石油集团东方地球物理勘探有限责任公司，
类型：发明
国别省市：