CN103513288A - 一种二维网格数据的补偿方向滤波方法 - Google Patents
一种二维网格数据的补偿方向滤波方法 Download PDFInfo
- Publication number
- CN103513288A CN103513288A CN201210211223.3A CN201210211223A CN103513288A CN 103513288 A CN103513288 A CN 103513288A CN 201210211223 A CN201210211223 A CN 201210211223A CN 103513288 A CN103513288 A CN 103513288A
- Authority
- CN
- China
- Prior art keywords
- node
- east
- filtering
- orientation
- south
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Landscapes
- Image Processing (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明涉及地球物理勘探的重力、磁力勘探数据处理技术,是一种二维网格数据的补偿方向滤波方法,将原始数据进行网格化处理,在窗口内划分方位角和每个方位对应的窗口内的节点分布,设置每个方位上背景场滤波时节点系数和剩余场补偿滤波求取时的节点系数,计算窗口内滤波方位的背景场值,循环迭代完成补偿方向滤波。本发明滤波的方向性强,满足了重磁处理提取方向异常信息的需要,同时也能满足压制某方向异常信息的需要。本发明推动了重磁勘探数据处理解释技术的进步,同时也为其它学科的二维网格数据的方向滤波处理提供了新方法。
Description
技术领域
本发明涉及地球物理勘探的重力、磁力勘探数据处理领域,属于数字滤波技术,是一种二维网格数据的补偿方向滤波方法。
背景技术
在磁力(包括航磁)资料采集中,由于测线间的一致性差异等原因引起的沿线的线性异常是具有方向的,这些方向性异常是地质研究的干扰(即数据噪音),因而,在重磁勘探资料处理解释中需要提取具有方向性的异常信息,以满足解释断层、构造的需要,需要在资料处理解释前进行处理,即利用方向滤波方法压制重磁数据中方向性的重磁噪音数据,使得滤波处理后的重磁数据更好地反映地质信息是需要解决的问题。
目前地球物理勘探的重力、磁力勘探数据处理中常用的频率域椭圆方向滤波是在频率域内实现,它的作用是提取出沿选择方向的方向性异常。这种方法常见的步骤是:1.将重力(或磁力)数据网格化,2.对重力网格数据进行傅立叶变换,3.选择滤波方向,选取椭圆方向滤波器,4.对傅立叶变换后的数据进行椭圆方向滤波,获得滤波后的频谱数据,5.对滤波后的频谱数据进行反傅立叶变换,即获得椭圆方向滤波后的重力(或磁力)异常网格数据。这种方法不具有迭代性和局部场补偿作用,而且方向性差。
地质出版社,2005.6,P194-204公开了一种《重力场与重力勘探》方法,利用空间域滑动窗口加权平均法建立一种圆滑滤波器,对数据具有圆滑滤波作用,可获得区域场,但是这种方法的方向滤波方法性能差,方向性异常往往被圆滑,无法提取重磁勘探资料处理解释中具有方向性的异常信息,不能利用方向滤波方法压制重磁数据中方向性的重磁噪音数据。
发明内容
本发明目的在于为重磁数据处理提供一种具有较好方向滤波功能,可提取方向性重磁异常数据或压制具有方向性的噪音数据的二维网格数据的补偿方向滤波方法。
本发明通过以下步骤实现:
1)进行野外重力或磁力勘探,采集得到原始数据并进行网格化处理;
2)运算窗口采用7行7列窗口,迭代次数取1~10;
步骤2)所述的运算窗口大小为7行7列。
3)在窗口内划分方位角和每个方位对应的窗口内的节点分布;
步骤3)所述的窗口内划分方位角是划分为南和北方位、东和西方位、北东和西南方位、北西和东南方位、北北东和南南西方位、北东东和南西西方位、北北西和南南东方位、北西西和南东东方位16个。
步骤3)所述的窗口内每对方位的节点分布相同,按以下分布:
南和北这对方位的节点分布是一组,窗口节点的行序为从下到上,窗口节点的列的序为从左到右,坐标组行序,列序为:南和北方位的窗口内节点分布为(1,4)、(2,4)、(3,4)、(4,4)、(5,4)、(6,4)、(7,4);
东和西方位的节点分布是一组,东和西方位的窗口内节点分布为(4,1)、(4,2)、(4,3)、(4,4)、(4,5)、(4,6)、(4,7);
北东和南西方位的节点分布是一组,北东和南西方位的窗口内节点分布为(1,1)、(2,2)、(3,3)、(4,4)、(5,5)、(6,6)、(7,7);
北西和南东方位的节点分布是一组,北西和南东方位的窗口内节点分布为(7,1)、(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,1)、(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,1)、(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)按照以下公式计算窗口内滤波方位的背景场值gbj:
gbj=(1/∑(k(i,j)))×(∑(k(i,j)×g(i,j)))
式中:
k(i,j)为节点(i,j)的系数值,g(i,j)为节点(i,j)的重力或磁力值;
步骤5)所述的滤波方位的背景场值gbj,当东西向滤波时,先计算计算点(4,4)的背景场值gbj(4,4):
gbj(4,4)=(1/42)×(3×g(4,1)+5×g(4,2)+7×g(4,3)+12×g(4,4)+7×g(4,5)+5×g(4,6)+3×g(4,7));
式中g(4,1)、g(4,2)、g(4,3)、g(4,4)、g(4,5)、g(4,6)、g(4,7)分别表示计算窗口内节点(4,1)、(4,2)、(4,3)、(4,4)、(4,5)、(4,6)、(4,7)的重力值或者是磁力值;
依次移动运算窗口,计算工区平面网格的每一个计算点的背景场值;
其它各方位滤波,替换所选坐标节点和对应节点的系数、按各方位的系数加权平均计算。
6)按照以下公式计算窗口内滤波方位的剩余场补偿滤波值gbch:
gbch=(1/∑(kbch(i,j)))×(∑(kbch(i,j)×gbch(i,j)))
式中:
kbch(i,j)为节点(i,j)的系数值kbch,g(i,j)为节点(i,j)的重力或磁力值;
步骤6)所述的剩余场补偿滤波值gbch,东西向滤波时:
gbch(4,4)=gbj(4,4)-((1/7)×(gbj(1,4)+gbj(2,4)+gbj(3,4)+gbj(4,4)+gbj(5,4)+gbj(6,4)+gbj(7,4)));
然后,计算计算点的补偿方向滤波输出值gf(4,4),
gf(4,4)=g(4,4)-gbj(4,4)+gbch(4,4);
再计算该点的剩余异常值gs(4,4),
gs(4,4)=g(4,4)-gf(4,4);
依次移动运算窗口,计算每一个计算点的补偿方向滤波值gf1(i,j)、剩余异常值gs1(i,j)。
7)循环迭代直到完成迭代运算次数,完成补偿方向滤波。
所述的循环迭代是将第一次计算得到的各个节点的补偿方向滤波值gf计为gf1(i,j),将第一次计算得到的各个节点的剩余异常值gs计为gs1(i,j),对gs1(i,j)进行下一次的补偿方向滤波运算,即把gs1(i,j)视为步骤5)中的原始数据g(i,j),进行迭代,计算步骤5)、步骤6),获得新的g’f(i,j)、g’s(i,j),则第二次迭代的结果为:
gf2(i,j)=gf1(i,j)+g’f(i,j).
gs2(i,j)=gs1(i,j)-g’f(i,j).
式中(i,j)表示计算点的节点在工区坐标,g’f、g’s分别表示本次迭代中计算得到的计算点的补偿方向滤波、剩余异常值。
本发明对实测数据处理,获得了多次迭代后的补偿方向滤波异常和补偿方向滤波剩余异常,滤波的方向性强,满足了重磁处理提取方向异常信息的需要,同时也能满足压制某方向异常信息的需要。本发明推动了重磁勘探数据处理解释技术的进步,同时也为其它学科的二维网格数据的方向滤波处理提供了一种新方法。
附图说明
图1是本发明二维网格的窗口选取和16方位节点分布图;
图2是本发明背景场滤波节点系数和补偿场滤波节点系数分布图;
图3是原始磁异常平面图;
图4是本发明南北向滤波磁异常平面图;
图5是本发明南北向滤波剩余磁异常平面图。
具体实施方式
以下结合实例详细说明本发明。
本实例通过以下步骤实现:
1)进行野外航空磁力勘探,采集得到航空磁力原始数据并进行网格化处理,见图3;
2)运算窗口采用7行7列窗口(见图1),迭代次数取2;
步骤2)所述的运算窗口大小为7行7列,见图1。
3)在窗口内划分方位角和每个方位对应的窗口内的节点分布,见图1;
步骤3)所述的窗口内划分方位角是划分为南和北方位、东和西方位、北东和西南方位、北西和东南方位、北北东和南南西方位、北东东和南西西方位、北北西和南南东方位、北西西和南东东方位16个,见图1。
步骤3)所述的窗口内每对方位的节点分布相同,按以下分布:
南和北这对方位的节点分布是一组,见图1的1,窗口节点的行序为从下到上,窗口节点的列的序为从左到右,坐标组行序,列序为:南和北方位的窗口内节点分布为(1,4)、(2,4)、(3,4)、(4,4)、(5,4)、(6,4)、(7,4);
东和西方位的节点分布是一组,见图1的2,东和西方位的窗口内节点分布为(4,1)、(4,2)、(4,3)、(4,4)、(4,5)、(4,6)、(4,7);
北东和南西方位的节点分布是一组,见图1的3,北东和南西方位的窗口内节点分布为(1,1)、(2,2)、(3,3)、(4,4)、(5,5)、(6,6)、(7,7);
北西和南东方位的节点分布是一组,见图1的4,北西和南东方位的窗口内节点分布为(7,1)、(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,1)、(3,2)、(3,3)、(4,3)、(4,4)、(4,5)、(5,5)、(5,6)、(5,7);
北北西和南南东方位的节点分布是一组,见图1的7,北北西和南南东方位的窗口内节点分布为(1,5)、(2,5)、(3,5)、(3,4)、(4,4)、(5,4)、(5,3)、(6,3)、(7,3);
北西西和南东东方位的节点分布是一组,见图1的8,北西西和南东东方位的窗口内节点分布为(5,1)、(5,2)、(5,3)、(4,3)、(4,4)、(4,5)、(3,5)、(3,6)、(3,7)。
4)设置南和北方位上背景场滤波时节点系数和剩余场补偿滤波求取时的节点系数;
图2是背景场滤波节点系数和补偿场滤波节点系数分布图。图中1为东-西向,2为南-北向,3为北西-东南向,4为东北-西南向,5为北西西-南南东向,6为北北西-南南东向,7为北东东-西西南向,8为北北东-南南西向。
步骤4)所述的节点系数为:
南北向滤波时,背景场滤波(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,节点系数分布见图2的1;
5)按照以下公式计算窗口内滤波方位的背景场值gbj:
gbj=(1/∑(k(i,j)))×(∑(k(i,j)×g(i,j)))
式中:k(i,j)为节点(i,j)的系数值,g(i,j)为节点(i,j)的重力或磁力值;
步骤5)所述的滤波方位的背景场值gbj,南北方向滤波时,先计算计算点(4,4)的背景场值gbj(4,4):
gbj(4,4)=(1/42)×(3×g(1,4)+5×g(2,4)+7×g(3,4)+12×g(4,4)+7×g(5,4)+5×g(6,4)+3×g(7,4));
式中g(1,4)、g(2,4)、g(3,4)、g(4,4)、g(5,4)、g(6,4)、g(7,4)分别表示计算窗口内节点1,4)、(2,4)、(3,4)、(4,4)、(5,4)、(6,4)、(7,4)的磁力值;
依次移动运算窗口,计算工区平面网格的每一个计算点的背景场值。
6)按照以下公式计算窗口内滤波方位的剩余场补偿滤波值gbch:
gbch=(1/∑(kbch(i,j)))×(∑(kbch(i,j)×gbch(i,j)))
式中:
kbch(i,j)为节点(i,j)的系数值kbch,g(i,j)为节点(i,j)的重力或磁力值;
步骤6)所述的剩余场补偿滤波值gbch,南北向滤波时:
gbch(4,4)=gbj(4,4)-((1/7)×(gbj(4,1)+gbj(4,2)+gbj(4,3)+gbj(4,4)+gbj(4,5)+gbj(4,6)+gbj(4,7)));
然后,计算计算点的补偿方向滤波输出值gf(4,4),
gf(4,4)=g(4,4)-gbj(4,4)+gbch(4,4);
再计算该点的剩余异常值gs(4,4),
gs(4,4)=g(4,4)-gf(4,4);
依次移动运算窗口,计算每一个计算点的补偿方向滤波值gf1(i,j)、剩余异常值gs1(i,j)。
7)循环迭代直到完成2次迭代运算,完成补偿方向滤波。
所述的循环迭代是将第一次计算得到的各个节点的补偿方向滤波值gf计为gf1(i,j),将第一次计算得到的各个节点的剩余异常值gs计为gs1(i,j),对gs1(i,j)进行下一次的补偿方向滤波运算,即把gs1(i,j)视为步骤5)中的原始数据g(i,j),进行迭代,计算步骤5)、步骤6),获得新的g’f(i,j)、g’s(i,j),则第二次迭代的结果为:
gf2(i,j)=gf1(i,j)+g’f(i,j).
gs2(i,j)=gs1(i,j)-g’f(i,j).
式中(i,j)表示计算点的节点在工区坐标,g’f、g’s分别表示本次迭代中计算得到的计算点的补偿方向滤波航磁异常(见图4)、补偿方向滤波剩余磁异常平面图,见图5。
对比图3、图4和图5,可见采用本方法实施南北向补偿方向滤波,获得了南北向方向异常图(见图4)和消除南北向条带磁异常现象的剩余磁异常(见图5),滤波分离的干净、效果明显。本发明对实测航磁数据进行了实验处理,获得了明显的滤波效果。
Claims (8)
1.一种二维网格数据的补偿方向滤波方法,特点是通过以下步骤实现:
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)循环迭代直到完成迭代运算次数,完成补偿方向滤波。
2.根据权利要求1所述的方法,特点是步骤2)所述的运算窗口大小为7行7列。
3.根据权利要求1所述的方法,特点是步骤3)所述的窗口内划分方位角是划分为南和北方位、东和西方位、北东和西南方位、北西和东南方位、北北东和南南西方位、北东东和南西西方位、北北西和南南东方位、北西西和南东东方位16个。
4.根据权利要求1所述的方法,特点是步骤3)所述的窗口内每对方位的节点分布相同,按以下分布:
南和北这对方位的节点分布是一组,窗口节点的行序为从下到上,窗口节点的列的序为从左到右,坐标组行序,列序为:南和北方位的窗口内节点分布为(1,4)、(2,4)、(3,4)、(4,4)、(5,4)、(6,4)、(7,4);东和西方位的节点分布是一组,东和西方位的窗口内节点分布为(4,1)、(4,2)、(4,3)、(4,4)、(4,5)、(4,6)、(4,7);
北东和南西方位的节点分布是一组,北东和南西方位的窗口内节点分布为(1,1)、(2,2)、(3,3)、(4,4)、(5,5)、(6,6)、(7,7);
北西和南东方位的节点分布是一组,北西和南东方位的窗口内节点分布为(7,1)、(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,1)、(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,1)、(5,2)、(5,3)、(4,3)、(4,4)、(4,5)、(3,5)、(3,6)、(3,7)。
5.根据权利要求1所述的方法,特点是步骤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;
其它各个方位的背景场滤波时节点系数和剩余场补偿滤波求取时的节点系数设置采用东西向滤波和南北向滤波设置方法类推。
6.根据权利要求1所述的方法,特点是步骤5)所述的滤波方位的背景场值gbj,当东西向滤波时,先计算计算点(4,4)的背景场值gbj(4,4):
gbj(4,4)=(1/42)×(3×g(4,1)+5×g(4,2)+7×g(4,3)+12×g(4,4)+7×g(4,5)+5×g(4,6)+3×g(4,7));
式中g(4,1)、g(4,2)、g(4,3)、g(4,4)、g(4,5)、g(4,6)、g(4,7)分别表示计算窗口内节点(4,1)、(4,2)、(4,3)、(4,4)、(4,5)、(4,6)、(4,7)的重力值或者是磁力值;
依次移动运算窗口,计算工区平面网格的每一个计算点的背景场值;
其它各方位滤波,替换所选坐标节点和对应节点的系数、按各方位的系数加权平均计算。
7.根据权利要求1所述的方法,特点是步骤6)所述的剩余场补偿滤波值gbch,东西向滤波时:
gbch(4,4)=gbj(4,4)-((1/7)×(gbj(1,4)+gbj(2,4)+gbj(3,4)+gbj(4,4)+gbj(5,4)+gbj(6,4)+gbj(7,4)));
然后,计算计算点的补偿方向滤波输出值gf(4,4),
gf(4,4)=g(4,4)-gbj(4,4)+gbch(4,4);
再计算该点的剩余异常值gs(4,4),
gs(4,4)=g(4,4)-gf(4,4);
依次移动运算窗口,计算每一个计算点的补偿方向滤波值gf1(i,j)、剩余异常值gs1(i,j)。
8.根据权利要求1所述的方法,特点是步骤7)所述的循环迭代是将第一次计算得到的各个节点的补偿方向滤波值gf计为gf1(i,j),将第一次计算得到的各个节点的剩余异常值gs计为gs1(i,j),对gs1(i,j)进行下一次的补偿方向滤波运算,即把gs1(i,j)视为步骤5)中的原始数据g(i,j),进行迭代,计算步骤5)、步骤6),获得新的g’f(i,j)、g’s(i,j),则第二次迭代的结果为:
gf2(i,j)=gf1(i,j)+g’f(i,j);
gs2(i,j)=gs1(i,j)-g’f(i,j);
式中(i,j)表示计算点的节点在工区坐标,g’f、g’s分别表示本次迭代中计算得到的计算点的补偿方向滤波、剩余异常值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210211223.3A CN103513288B (zh) | 2012-06-21 | 2012-06-21 | 一种二维网格数据的补偿方向滤波方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210211223.3A CN103513288B (zh) | 2012-06-21 | 2012-06-21 | 一种二维网格数据的补偿方向滤波方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103513288A true CN103513288A (zh) | 2014-01-15 |
CN103513288B CN103513288B (zh) | 2016-06-08 |
Family
ID=49896299
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210211223.3A Active CN103513288B (zh) | 2012-06-21 | 2012-06-21 | 一种二维网格数据的补偿方向滤波方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103513288B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103577607A (zh) * | 2013-11-20 | 2014-02-12 | 哈尔滨工程大学 | 一种基于地磁异常数据形态特征的边界补偿方法 |
CN105891889A (zh) * | 2016-03-31 | 2016-08-24 | 中国石油天然气集团公司 | 一种重力异常边界增强的方法及装置 |
CN113627051A (zh) * | 2021-07-23 | 2021-11-09 | 中国地质科学院地球物理地球化学勘查研究所 | 一种重力异常场分离方法、系统、存储介质和电子设备 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4570122A (en) * | 1983-03-30 | 1986-02-11 | Mobil Oil Corporation | Magnetic exploration with reduction of magnetic data to the equator |
US20040004480A1 (en) * | 2002-07-08 | 2004-01-08 | Von Wimmerspeg Udo | Precision grid survey apparatus and method for the mapping of hidden ferromagnetic structures |
CN101661115A (zh) * | 2008-08-29 | 2010-03-03 | 中国石油集团东方地球物理勘探有限责任公司 | 基于标准格架的快速三维重力、磁力物性反演的方法 |
CN102338890A (zh) * | 2010-10-22 | 2012-02-01 | 中国石油天然气股份有限公司 | 一种地球物理勘探中的圆形窗带通保幅滤波数据处理方法 |
CN102338884A (zh) * | 2010-10-22 | 2012-02-01 | 中国石油天然气股份有限公司 | 物探中的椭圆窗方向带通保幅滤波数据处理方法 |
-
2012
- 2012-06-21 CN CN201210211223.3A patent/CN103513288B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4570122A (en) * | 1983-03-30 | 1986-02-11 | Mobil Oil Corporation | Magnetic exploration with reduction of magnetic data to the equator |
US20040004480A1 (en) * | 2002-07-08 | 2004-01-08 | Von Wimmerspeg Udo | Precision grid survey apparatus and method for the mapping of hidden ferromagnetic structures |
CN101661115A (zh) * | 2008-08-29 | 2010-03-03 | 中国石油集团东方地球物理勘探有限责任公司 | 基于标准格架的快速三维重力、磁力物性反演的方法 |
CN102338890A (zh) * | 2010-10-22 | 2012-02-01 | 中国石油天然气股份有限公司 | 一种地球物理勘探中的圆形窗带通保幅滤波数据处理方法 |
CN102338884A (zh) * | 2010-10-22 | 2012-02-01 | 中国石油天然气股份有限公司 | 物探中的椭圆窗方向带通保幅滤波数据处理方法 |
Non-Patent Citations (1)
Title |
---|
郭志宏: "航磁及梯度数据正反演解释方法技术实用化改进及应用", 《中国优秀博硕士学位论文全文数据库(博士) 基础科学辑》, no. 11, 15 November 2006 (2006-11-15), pages 1 - 81 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103577607A (zh) * | 2013-11-20 | 2014-02-12 | 哈尔滨工程大学 | 一种基于地磁异常数据形态特征的边界补偿方法 |
CN103577607B (zh) * | 2013-11-20 | 2017-06-20 | 哈尔滨工程大学 | 一种基于地磁异常数据形态特征的边界补偿方法 |
CN105891889A (zh) * | 2016-03-31 | 2016-08-24 | 中国石油天然气集团公司 | 一种重力异常边界增强的方法及装置 |
CN105891889B (zh) * | 2016-03-31 | 2018-05-04 | 中国石油天然气集团公司 | 一种重力异常边界增强的方法及装置 |
CN113627051A (zh) * | 2021-07-23 | 2021-11-09 | 中国地质科学院地球物理地球化学勘查研究所 | 一种重力异常场分离方法、系统、存储介质和电子设备 |
CN113627051B (zh) * | 2021-07-23 | 2024-01-30 | 中国地质科学院地球物理地球化学勘查研究所 | 一种重力异常场分离方法、系统、存储介质和电子设备 |
Also Published As
Publication number | Publication date |
---|---|
CN103513288B (zh) | 2016-06-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US10739482B2 (en) | Surface wave prospecting method and acquisition equipment | |
CN111164462B (zh) | 一种人工源面波勘探方法、面波勘探装置及终端设备 | |
CN104932010B (zh) | 一种基于近道镶边稀疏Radon变换的绕射波分离方法 | |
CN109164483A (zh) | 多分量地震数据矢量去噪方法及多分量地震数据矢量去噪装置 | |
CN104793253A (zh) | 基于数学形态学的航空电磁数据去噪方法 | |
CN105652320B (zh) | 逆时偏移成像方法和装置 | |
CN107561589A (zh) | 一种近地表横波层速度模型建立方法 | |
CN105068131A (zh) | 一种航磁数据调平方法 | |
CN102944905B (zh) | 一种基于方向小波分析的重磁异常处理方法 | |
CN105204065A (zh) | 一种拾取初至波的方法和装置 | |
CN107861156B (zh) | 绕射波的提取方法及装置 | |
CN102736108B (zh) | 基于样条拟合的真三维地震数据噪声压制方法 | |
CN103513288A (zh) | 一种二维网格数据的补偿方向滤波方法 | |
CN104391325A (zh) | 不连续非均质地质体检测方法和装置 | |
CN110058319B (zh) | 一种大地电磁数据采集方法、装置及终端设备 | |
CN102338884B (zh) | 物探中的椭圆窗方向带通保幅滤波数据处理方法 | |
CN104765069A (zh) | 一种压制同步激发采集邻炮干扰的方法 | |
CN112946760B (zh) | 基于正则化方法的未爆弹三维立体成像方法、装置及系统 | |
CN109143344A (zh) | 绕射波分离方法、装置和电子设备 | |
CN102608658A (zh) | 强能量干扰抑制方法 | |
CN106842322B (zh) | 一种二氧化碳驱油监控地震时差校正方法 | |
CN105891889A (zh) | 一种重力异常边界增强的方法及装置 | |
CN104155692B (zh) | 获取微地震的震源点的方法及装置 | |
CN111352158B (zh) | 地震信号增强方法及装置 | |
CN111337973B (zh) | 地震数据重建方法、系统 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |