CN1877366A - 重磁延拓回返垂直导数目标优化处理技术 - Google Patents

重磁延拓回返垂直导数目标优化处理技术 Download PDF

Info

Publication number
CN1877366A
CN1877366A CN 200610098673 CN200610098673A CN1877366A CN 1877366 A CN1877366 A CN 1877366A CN 200610098673 CN200610098673 CN 200610098673 CN 200610098673 A CN200610098673 A CN 200610098673A CN 1877366 A CN1877366 A CN 1877366A
Authority
CN
China
Prior art keywords
technology
heavy
continuation
returning
vertical
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
Application number
CN 200610098673
Other languages
English (en)
Other versions
CN100434934C (zh
Inventor
杨辉
邹才能
张研
李明
李建忠
文百红
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China Petroleum and Natural Gas Co Ltd
Original Assignee
杨辉
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by 杨辉 filed Critical 杨辉
Priority to CNB2006100986730A priority Critical patent/CN100434934C/zh
Publication of CN1877366A publication Critical patent/CN1877366A/zh
Application granted granted Critical
Publication of CN100434934C publication Critical patent/CN100434934C/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Geophysics And Detection Of Objects (AREA)

Abstract

一种在压制重磁异常高频噪音的前提下,提高重磁异常分辨率的处理方法。特征是:首先将重磁不规则网格数据网格化,形成规则网数据;对规则网的重磁数据进行傅氏正变换,求出原始数据的频谱;进行频谱分析,分析信号、噪音的分布区间,确定滤波参数;利用解析法技术向上延拓;利用差分法技术向下回返;解析向上延拓和差分法向下回返过程重复p次;与广义垂直n次导数滤波算子相乘;经过傅氏反变换,求出处理后的数据;用绘图软件绘制处理后数据的平面图或者立体图。效果是:通过将延拓回返技术与广义垂直n次导数技术组合,获得最好的频率响应滤波器,在滤除高频噪音与提高分辨率之间获得平衡,突出不同目标层段异常,有效提高重磁异常分辨率。

Description

重磁延拓回返垂直导数目标优化处理技术
技术领域
本发明涉及地球物理勘探中的重力勘探和磁法勘探处理技术领域,特别是一种在压制重磁异常中的高频噪音的前提下,提高重磁异常分辨率的处理方法。
背景技术
地质勘探主要应用的地球物理勘探方法有重力勘探、磁法勘探、电法勘探、地震勘探等。
重力勘探是通过利用重力仪观测地下物质密度差异引起的重力异常以查明地下的地质构造和岩性异常体。磁法勘探则是利用磁力仪观测地下物体磁性差异引起的磁力异常以查明地下的地质构造和磁性异常体如火山岩等。
重磁勘探可以分为三个大的环节:野外重磁异常采集、重磁异常处理、重磁异常解释。
重磁异常处理的主要任务是根据位场理论和滤波原理,利用电子计算机等设备和相应的资料处理软件,对野外采集的重磁异常进行加工处理,以获得能反映地下地层结构和地层岩性变化的平面和剖面分布特征,供解释人员确定含油气盆地的坳隆结构、有利二级构造带以及有利目标。
重磁异常处理过程主要包括以下几个处理步骤:1、预处理;2、网格化;3、压制干扰;4、位场转换。
重磁异常处理的重要目的是提高资料的横向分辨率,突出地下场源体引起的异常,以更好的反映局部构造、断裂等地质现象。根据原始资料的精度以及要解决的不同勘探对象的地质问题,如区域构造、局部构造、深层构造、浅层构造等。可以采用垂直一次导数、垂直二次导数、垂直四次导数技术进行资料的处理。到目前为止,还未有介绍及使用垂直三次导数技术的报道,并且导数的次数仅为1、2或4的正整数。导数的次数为正整数有明确的物理意义。垂直导数是一个比较有效的异常增强方法之一,该种计算可以在空间域或频率域实现。空间域滤波的方法如Elikins、Rosenbach、王家林、徐世浙、王光杰、李明、杨辉、王宜昌等。频率域滤波的方法如王家林、姚长利、Gupta等。与空间域滤波方法相比,频率域滤波的主要优点在于用频率域乘积代替了空间域的褶积,具有实现方便、滤波运算比较简单,处理计算量小、计算速度快,无数据边部损失等,所以目前人们基本上都是在频率域中进行重磁异常的处理。
垂直n次导数理论上相当于一个高频放大器,它对高频放大作用特别明显,而实际重磁观测资料中存在误差及浅部影响往往是高频的,造成了重磁处理的不稳定性。为此,国内外学者进行了大量的研究,其中主要是采取串联一个低通滤波器压制高频干扰的方法,如运算时常常加一个圆滑因子等低通滤波器对高频放大作用进行抑制。低通滤波器因子选择非常重要,也往往较难确定,易造成对观测异常压制过大或对高频压制不够等缺点,选择一个好的低通滤波器尤为重要。
通过大量的理论模型试验,优选出延拓回返滤波作为低通滤波器,该方法通过多次的延拓回返来压制高频噪音,最终获取低频信息。
中国专利公开号:1667433,公开了一种位场高分辨率视深度滤波方法。该申请的操作过程是:
a、对视深度滤波在波数域中的算子FPDF乘上水平一次导数算子Fgx或垂直一次导数算子Fgz或垂直二次导数算子Fgzz,得到高分辨率视深度滤波在波数域中的算子FHPDF。
①采用水平一阶导数的高分辨率视深度滤波波数域的算子为:
F HPDF = ( i 2 πu ) [ exp ( - 2 πh u 2 + v 2 ) ] p ( 6 - exp ( - 2 πh u 2 + v 2 - 2 cos ( 2 πhu ) - 2 cos ( 2 πhv ) ] p + 1
②采用垂直一阶导数的高分辨率视深度滤波波数域的算子为:
F HPDF = ( 2 π u 2 + v 2 ) [ exp ( - 2 πh u 2 + v 2 ) ] p ( 6 - exp ( - 2 πh u 2 + v 2 - 2 cos ( 2 πhu ) - 2 cos ( 2 πhv ) ] p + 1
③采用垂直二阶导数的高分辨率视深度滤波波数域的算子为:
F HPDF = ( 4 π 2 ( u 2 - v 2 ) ) [ exp ( - 2 πh u 2 + v 2 ) ] p ( 6 - exp ( - 2 πh u 2 + v 2 - 2 cos ( 2 πhu ) - 2 cos ( 2 πhv ) ] p + 1
式中:h为延拓高度;u,v分别为x和y方向的波数;p为回返次数
b、对原始观测数据进行傅氏变换,然后与高分辨率视深度滤波在波数域中的算子FHPDF相乘,再经过傅氏反变换,得到高分辨率视深度滤波在2D/3D空间分布的数据体。
c、通过绘图软件可绘制2D断面图或者3D立体图、水平切片图、垂直切片图;用于研究局部或断裂构造在平面及空间的分布规律及确定场源的质心深度。
该专利申请的滤波方法原理是将位场数据下延到地面以下研究地下不同深度的位场数据的分布,采用了水平一次导数以及垂直一次导数、垂直二次导数方法处理,以提高分辨率,目的是确定局部构造及断裂的分布及质心深度。该滤波方法特别适用于大比例尺位场资料的处理以及确定位场场源的质心深度。该滤波方法在其适用范围和所要达到的目的与本发明不相同,其原理不同,实现的方法步骤也不相同。本发明的原理是在观测面上进行资料的目标处理,步骤方法是通过延拓回返压制噪音,针对不同的勘探深度的目标,通过不同次数的导数来提高分辨率,并且,从滤波器的角度将传统导数次数的正整数概念拓展到了正有理数,使导数滤波器成为一个广义的连续变化的滤波器。资料处理的目的是突出不同深度的勘探目标,实现资料的目标处理,适用的范围是大、中、小比例尺重磁异常均适用,效果也不相同。
发明内容
本发明的目的在于提供一种重磁延拓回返垂直导数目标优化处理技术,在压制高频噪音的前提下,提高重磁异常分辨率的处理方法。通过不同参数的选择,突出不同目标层段的异常。克服传统的垂直二次导数在对信号放大的同时,对高频噪音同时放大的不足,造成出现虚假异常,使地质解释人员在分析资料时,真假难辨,发生误判,出现与实际情况不符的解释、判断,造成勘探的失误。
本发明采用的技术方案是:
重磁延拓回返垂直导数目标优化处理技术,它是一种在压制重磁异常高频噪音的前提下,提高重磁异常分辨率的处理方法。根据重磁原始资料的精度以及要解决的不同勘探对象的地质问题,如区域构造、局部构造、深层构造、浅层构造等,采用广义垂直n次导数技术进行重磁异常处理,并通过多次的向上延拓和向下回返处理,压制高频噪音,突出不同目标层段的异常。重磁异常的处理需要在计算机上进行,重磁异常的处理过程依次含有以下9个步骤:
①首先将不规则网格观测的重磁异常网格化,形成规则网数据异常。
②对规则网的重磁异常进行傅氏正变换,得到原始异常的频谱。
③进行频谱分析:在频谱图上根据信号、噪音的分布区间,振幅大小及相位差异确定滤波参数,滤波参数包括导数的次数n、向上延拓和向下回返的高度h,以及向上延拓和向下回返的次数p。
④利用解析法技术进行向上延拓,向上延拓采用的公式是:
F up = exp ( - 2 πh u 2 + v 2 ) ;
式中:u,v分别为x和y方向的波数,高度h的取值范围在500——3000米之间。
⑤利用差分法技术向下回返,向下回返采用的公式是:
F down = ( 6 - exp ( - 2 πh u 2 + v 2 ) - 2 cos ( 2 πhu ) - 2 cos ( 2 πhv ) ) .
⑥重复解析法技术向上延拓和差分法技术向下回返的过程,重复p次;p为5-20的整数。
⑦在上述处理的基础上,与广义垂直n次导数滤波算子相乘,垂直n次导数滤波器的公式是:
F = ( 2 π u 2 + v 2 ) n ;
式中:n为0.1-5.0的正有理数。
n为0.1-5.0的正有理数本发明的创新点之一,拓展了传统导数次数为正整数的概念。效果是:通过广义导数次数的选择,突出更多不同深度地质体所产生的异常,实现重磁异常的目标处理。
⑧再经过傅氏反变换,获得重磁异常处理后的重磁信号。
⑨用绘图软件将处理后的重磁信号绘制成平面图或者立体图,然后对处理之后的重磁异常结合其它勘探方法获得的资料进行对比、分析及解释。
傅氏正变换是将时间域的数据通过傅氏变换变换到频率域。
傅氏反变换是将频率域的数据通过傅氏变换变换到时间域。
解析法向上延拓是用解析的方法,根据观测面上的实测重磁异常确定观测面以上空间的某一高度上重磁异常。解析法向上延拓是适定的。
差分法向下回返是用差分的方法,将解析法向上延拓的某一高度上重磁异常回返到原观测面上。差分法向下回返是不适定的,因此,在回返的过程中,高频成分得到了不同程度的衰减。
上述离散数据网格化和确定滤波参数的过程,本专业技术领域分析人员能熟练完成。
本发明方法适用于大、中、小比例尺重磁异常的处理。该方法原理是在地面上研究重磁异常的分布情况。与中国专利公开号:1667433的一种位场高分辨率视深度滤波方法相比较,在适用范围和所要达到的目的与本发明不相同,其原理不同,实现的方法步骤也不相同。
本发明的有益效果:本发明的优点在于在压制高频噪音的前提下,提高重磁异常分辨率,通过将向上延拓和向下回返技术与广义垂直n次导数技术的组合,获得了一个最好的频率响应滤波器,在滤除高频噪音与提高分辨率之间获到了一种最佳平衡。n为0.1-5.0的正有理数本发明的创新点之一,打破了传统导数次数只能为正整数的概念,通过广义导数次数n和延拓回返高度以及回返次数的组合选择,突出更多不同目标层段的异常,实现重磁异常的目标处理。
本发明的上述重磁延拓回返垂直导数目标优化处理技术有较大的改进,使其完全符合地质勘探科研及生产的需要并产生了好用及实用的效果。是一项新颖、进步、实用的新的重磁异常处理方法。
下面结合附图和具体实施例对本发明作进一步详细说明。
附图说明
图1重磁延拓回返垂直导数目标优化处理过程的计算机程序流程图
图2本发明与传统滤波器频谱对比图
图3本发明滤波器频谱特征图
图4济阳坳陷某地区重力资料传统二次导数计算结果图
图5济阳坳陷某地区重力资料本发明技术计算结果图
图6华北某地区重力资料本发明技术计算结果图
图7松辽盆地北部徐家围子断陷磁力资料本发明技术计算结果图
图8松辽盆地南部长岭断陷磁力资料本发明技术计算结果图
图9松辽盆地北部某断陷重力资料本发明技术计算结果图
具体实施方式
实施例1:参阅附图1。以济阳坳陷某地区重力资料为例,对传统的二次导数与本发明技术进行效果对比。本发明的方法为在室内处理时,按照发明内容所述的9个步骤对重力资料进行处理,即:
①首先将不规则网格观测的重力异常网格化,形成规则网数据;
②对规则网的重力数据进行傅氏正变换,得到原始异常的频谱;
③进行频谱分析:在频谱图上根据信号、噪音的分布区间,振幅大小及相位差异确定滤波参数,导数的次数n=2、向上延拓和向下回返的高度h=750m,以及向上延拓和向下回返的次数p=10;
④利用解析法技术进行向上延拓,向上延拓采用的公式是:
F up = exp ( - 2 πh u 2 + v 2 ) ;
⑤利用差分法技术向下回返,向下回返采用的公式是:
F down = ( 6 - exp ( - 2 πh u 2 + v 2 ) - 2 cos ( 2 πhu ) - 2 cos ( 2 πhv ) ) ;
⑥重复解析法技术向上延拓和差分法技术向下回返的过程要重复10次;
⑦在上述处理的基础上,与垂直2次导数滤波算子相乘,垂直2次导数滤波器的公式是:
F = ( 2 π u 2 + v 2 ) 2 ;
⑧再经过傅氏反变换,获得重力资料处理后的重磁信号;
⑨用绘图软件将处理后的重磁信号绘制成平面图或者立体图,然后对处理之后的重力资料结合其它资料进行对比、分析及解释。
参阅附图2。图2是本发明与传统滤波器频谱对比图。本发明的滤波器为一个主频较低(f0=0.178Hz)的窄带滤波器,对信号进行放大,对噪音进行压制,带通效果好,因此,处理结果图件异常规律性较强;传统的二次导数为一个主频较高(f0=0.318Hz)的宽带滤波器,对信号与噪音均放大,带通效果不好,因此,处理结果图件异常规律性较差,图面较乱。另外,本发明的滤波器滤波曲线较陡,在频率0.08-0.18Hz较传统二次导数滤波器数值要大,意味着在该频率段,本发明的滤波器具有较大的放大作用,更能突出传统二次导数突出不够较低频(如深层潜山构造等)的信号。
参阅附图3。图3本发明滤波器频谱特征图。本发明的滤波器随着导数阶次的增加,主频逐渐向高频移动,意味着通过导数阶次n的选择,可以突出不同目标层段的异常。
参阅附图4。图4为济阳坳陷某地区重力资料传统二次导数计算结果图。可以看到,图面较乱,存在南北向的半系统误差,两条北东向断裂不清,渤古1井、渤古2井位于重力低上,与地震及钻井资料不符。
参阅附图5。图5为济阳坳陷某地区重力资料本发明技术计算结果图。可以看到,图面较有规律,两条北东向断裂十分清楚,北面北东向断裂断面较缓;渤古1井位于重力高边部,渤古2处于高部位,与地震及钻井资料吻合
实施例2:参阅附图1。以下为华北某地区重力异常的延拓回返垂直导数目标优化处理。本发明的方法为在室内处理时,按照发明内容所述的9个步骤对重力资料进行处理,即:
①首先将不规则网格观测的重力异常网格化,形成规则网异常;
②对规则网的重力异常进行傅氏正变换,得到原始异常的频谱;
③进行频谱分析:在频谱图上根据信号、噪音的分布区间,振幅大小及相位差异确定滤波参数,导数的次数n=3、向上延拓和向下回返的高度h=1250m,以及向上延拓和向下回返的次数p=12;
④利用解析法技术进行向上延拓,向上延拓采用的公式是:
F up = exp ( - 2 πh u 2 + v 2 ) ;
⑤利用差分法技术向下回返,向下回返采用的公式是:
F down = ( 6 - exp ( - 2 πh u 2 + v 2 ) - 2 cos ( 2 πhu ) - 2 cos ( 2 πhv ) ) ;
⑥重复解析法技术向上延拓和差分法技术向下回返的过程要重复12次;
⑦在上述处理的基础上,与垂直3次导数滤波算子相乘,垂直3次导数滤波器的公式是:
F = ( 2 π u 2 + v 2 ) 3 ;
⑧再经过傅氏反变换,获得重力资料处理后的重磁信号;
⑨用绘图软件将处理后的重磁信号绘制成平面图或者立体图,然后对处理之后的重力资料结合其它资料进行对比、分析及解释。
参阅附图6。从图6上可以看到异常分布规律性较强,主要发育南北向、东西向异常,与该区地震资料及钻井资料吻合。另外,在工区的中部新发现了北西向两个有利的局部异常,推测是深部潜山的反映。在从前应用其它处理方法处理过该地区重力资料,处理后的解释认为本地不存在或存在较小的北西向局部异常。
实施例3:松辽盆地磁力异常的向上延拓和向下回返垂直导数目标优化处理过程。本发明的方法为在室内处理时,按照发明内容所述的9个步骤对磁力异常进行处理,即:
①首先将不规则网格观测的磁力异常网格化,形成规则网异常;
②对规则网的磁力异常进行傅氏正变换,得到原始异常的频谱;
③进行频谱分析:在频谱图上根据信号、噪音的分布区间,振幅大小及相位差异确定滤波参数,导数的次数n=2、向上延拓和向下回返的高度h=1500m,以及向上延拓和向下回返的次数p=8;
④利用解析法技术进行向上延拓,向上延拓采用的公式是:
F up = exp ( - 2 πh u 2 + v 2 ) ;
⑤利用差分法技术向下回返,向下回返采用的公式是:
F down = ( 6 - exp ( - 2 πh u 2 + v 2 ) - 2 cos ( 2 πhu ) - 2 cos ( 2 πhv ) ) ;
⑥重复解析法技术向上延拓和差分法技术向下回返的过程要重复8次;
⑦在上述处理的基础上,与垂直2次导数滤波算子相乘,垂直2次导数滤波器的公式是:
F = ( 2 π u 2 + v 2 ) 2 ;
⑧再经过傅氏反变换,得到磁力异常处理后的重磁信号;
⑨用绘图软件将处理后的重磁信号绘制成平面图或者立体图,然后对处理之后的磁力资料结合其它资料进行对比、分析及解释。
参阅附图7、附图8。图7为松辽盆地北部徐家围子断陷磁力资料的向上延拓和向下回返垂直二次导数目标优化处理结果图,经与该区50余口钻遇火山岩的深井对比后发现,该结果较好的反映了火山岩的分布情况。该区23口后验井资料对比表明,预测火山岩的成功率在80%以上。另外,推荐位于磁力资料梯度带上的探井优先上钻,徐深7、徐深9等井均获得了高产天然气流。图8为松辽盆地南部长岭断陷磁力资料的向上延拓和向下回返垂直二次导数目标优化处理结果图,利用该处理结果结合其它资料推荐上钻该断陷第一口风险探井长深1井,长深1井获天然气无阻流量百万方,松辽盆地南部首次取得了历史性突破。
我国的火山岩油气勘探刚刚开始,急需针对火山岩的勘探技术,该项技术对于我国松辽、准噶尔、四川等盆地的火山岩油气勘探有着非常重要的意义,是一项有重要实用价值的技术。
实施例4:参阅附图1。以下为松辽盆地北部某断陷重力异常的延拓回返垂直导数目标优化处理。本发明的方法为在室内处理时,按照发明内容所述的9个步骤对重力异常进行处理,即:
①首先将不规则网格观测的重力异常网格化,形成规则网异常;
②对规则网的重力异常进行傅氏正变换,得到原始异常的频谱;
③进行频谱分析:在频谱图上根据信号、噪音的分布区间,振幅大小及相位差异确定滤波参数,导数的次数n=2.5、向上延拓和向下回返的高度h=1000m,以及向上延拓和向下回返的次数p=8;
④利用解析法技术进行向上延拓,向上延拓采用的公式是:
F up = exp ( - 2 πh u 2 + v 2 ) ;
⑤利用差分法技术向下回返,向下回返采用的公式是:
F down = ( 6 - exp ( - 2 πh u 2 + v 2 ) - 2 cos ( 2 πhu ) - 2 cos ( 2 πhv ) ) ;
⑥重复解析法技术向上延拓和差分法技术向下回返的过程要重复8次;
⑦在上述处理的基础上,与垂直2.5次导数滤波算子相乘,垂直2.5次导数滤波器的公式是:
F = ( 2 π u 2 + v 2 ) 2.5 ;
⑧再经过傅氏反变换,获得重力资料处理后的重磁信号;
⑨用绘图软件将处理后的重磁信号绘制成平面图或者立体图,然后对处理之后的重力资料结合其它资料进行对比、分析及解释。
参阅附图9。该区重力异常处理的目的是确定火山岩的岩性。从图9可以看到钻遇中基性火山岩的的升深101、升深7、尚深1、尚深3等井位于重力高上,而钻遇酸性火山岩的升深2、徐深1井位于重力低上,火山岩岩性与重力异常间有较好的对应关系。结合其它资料预测了该区的火山岩岩性分布,位于重力低上的后验井升深203证实预测是正确的。

Claims (10)

1、一种重磁延拓回返垂直导数目标优化处理技术,其特征在于:在压制重磁异常高频噪音的前提下,提高重磁异常分辨率,根据重磁原始资料的精度以及要解决的不同勘探对象的地质问题,如区域构造、局部构造、深层构造或浅层构造,采用广义垂直n次导数技术进行重磁异常的处理,并通过多次的向上延拓和向下回返处理,压制高频噪音,突出不同目标层段的异常,重磁异常的处理需要在计算机上进行,重磁异常的处理过程依次含有以下步骤:
①首先将不规则网格观测的重磁异常网格化,形成规则网数据异常;
②对规则网的重磁异常进行傅氏正变换,获得原始异常的频谱;
③进行频谱分析:在频谱图上根据信号、噪音的分布区间,振幅大小及相位差异确定滤波参数,滤波参数包括导数的次数n、向上延拓和向下回返的高度h,以及向上延拓和向下回返的次数p;
④利用解析法技术进行向上延拓,向上延拓采用的公式是:
F up = exp ( - 2 πh u 2 + v 2 ) ;
式中:u,v分别为x和y方向的波数,高度h的取值范围在500--3000米之间;
⑤利用差分法技术向下回返,向下回返采用的公式是:
F down = ( 6 - exp ( - 2 πh u 2 + v 2 ) - 2 cos ( 2 πhu ) - 2 cos ( 2 πhv ) ) ;
⑥重复解析法技术向上延拓和差分法技术向下回返的过程,重复p次;p为5-20的整数;
⑦在上述处理的基础上,与广义垂直n次导数滤波算子相乘,广义垂直n次导数滤波器的公式是:
F = ( 2 π u 2 + v 2 ) n ,
式中:n为0.1-5.0的正有理数;
⑧再经过傅氏反变换,得到重磁异常处理后的重磁信号;
⑨用绘图软件将处理后的重磁信号绘制成平面图或者立体图,然后对处理之后的重磁异常结合其它勘探方法获得的资料进行对比、分析及解释。
2、根据权利要求1所述的重磁延拓回返垂直导数目标优化处理技术,其特征是:所述的p次,p为5-20。
3、根据权利要求2所述的重磁延拓回返垂直导数目标优化处理技术,其特征是:所述的p次,p为8-12。
4、根据权利要求1或2或3所述的重磁延拓回返垂直导数目标优化处理技术,其特征是:所述的n次,n为0.1-5.0的正有理数。
5、根据权利要求4所述的重磁延拓回返垂直导数目标优化处理技术,其特征是:所述的n次,n为1-4。
6、根据权利要求4所述的重磁延拓回返垂直导数目标优化处理技术,其特征是:所述的n次,n为0.5。
7、根据权利要求4所述的重磁延拓回返垂直导数目标优化处理技术,其特征是:所述的n次,n为1.5。
8、根据权利要求4所述的重磁延拓回返垂直导数目标优化处理技术,其特征是:所述的n次,n为2.5。
9、根据权利要求4所述的重磁延拓回返垂直导数目标优化处理技术,其特征是:所述的n次,n为3.5。
10、根据权利要求4所述的重磁延拓回返垂直导数目标优化处理技术,其特征是:所述的n次,n为3。
CNB2006100986730A 2006-07-12 2006-07-12 重磁延拓回返垂直导数目标优化处理方法 Expired - Fee Related CN100434934C (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CNB2006100986730A CN100434934C (zh) 2006-07-12 2006-07-12 重磁延拓回返垂直导数目标优化处理方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CNB2006100986730A CN100434934C (zh) 2006-07-12 2006-07-12 重磁延拓回返垂直导数目标优化处理方法

Publications (2)

Publication Number Publication Date
CN1877366A true CN1877366A (zh) 2006-12-13
CN100434934C CN100434934C (zh) 2008-11-19

Family

ID=37509850

Family Applications (1)

Application Number Title Priority Date Filing Date
CNB2006100986730A Expired - Fee Related CN100434934C (zh) 2006-07-12 2006-07-12 重磁延拓回返垂直导数目标优化处理方法

Country Status (1)

Country Link
CN (1) CN100434934C (zh)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102368094A (zh) * 2010-10-22 2012-03-07 中国石油天然气股份有限公司 一种地球物理勘探中的矩形窗带通保幅滤波数据处理方法
CN107291659A (zh) * 2017-05-16 2017-10-24 哈尔滨工程大学 平面地磁异常场一步向上延拓平面模量梯度场的递归余弦变换法
CN109902315A (zh) * 2017-12-07 2019-06-18 核工业北京地质研究院 一种圈定隐伏花岗岩岩体深部边界的方法
CN112464521A (zh) * 2020-10-28 2021-03-09 中国石油天然气集团有限公司 走滑断裂确定方法及装置
CN113686329A (zh) * 2021-08-27 2021-11-23 中国人民解放军国防科技大学 一种基于地磁数据的垂直高度位场测量方法

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
AU594781B2 (en) * 1986-06-16 1990-03-15 Stolar, Inc. Method for constructing vertical images of anomolies in geological formations
GB2417086B (en) * 2001-12-18 2006-08-16 Bhp Billiton Innovation Pty Method of processing marine magnetic gradient data and exploration methods using that data
AUPS306802A0 (en) * 2002-06-20 2002-07-11 Wmc Resources Ltd A data acquisition unit, system and method for geophysical data
GB2415511B (en) * 2004-06-26 2008-09-24 Statoil Asa Processing electromagnetic data

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102368094A (zh) * 2010-10-22 2012-03-07 中国石油天然气股份有限公司 一种地球物理勘探中的矩形窗带通保幅滤波数据处理方法
CN102368094B (zh) * 2010-10-22 2013-10-16 中国石油天然气股份有限公司 一种地球物理勘探中的矩形窗带通保幅滤波数据处理方法
CN107291659A (zh) * 2017-05-16 2017-10-24 哈尔滨工程大学 平面地磁异常场一步向上延拓平面模量梯度场的递归余弦变换法
CN107291659B (zh) * 2017-05-16 2020-09-25 哈尔滨工程大学 平面地磁异常场一步向上延拓平面模量梯度场的递归余弦变换法
CN109902315A (zh) * 2017-12-07 2019-06-18 核工业北京地质研究院 一种圈定隐伏花岗岩岩体深部边界的方法
CN109902315B (zh) * 2017-12-07 2023-06-16 核工业北京地质研究院 一种圈定隐伏花岗岩岩体深部边界的方法
CN112464521A (zh) * 2020-10-28 2021-03-09 中国石油天然气集团有限公司 走滑断裂确定方法及装置
CN112464521B (zh) * 2020-10-28 2024-05-28 中国石油天然气集团有限公司 走滑断裂确定方法及装置
CN113686329A (zh) * 2021-08-27 2021-11-23 中国人民解放军国防科技大学 一种基于地磁数据的垂直高度位场测量方法
CN113686329B (zh) * 2021-08-27 2023-07-25 中国人民解放军国防科技大学 一种基于地磁数据的垂直高度位场测量方法

Also Published As

Publication number Publication date
CN100434934C (zh) 2008-11-19

Similar Documents

Publication Publication Date Title
CN1252492C (zh) 大地电磁波电阻率测量方法及其仪器
CN108931814B (zh) 一种基于多属性融合的基岩裂缝预测的方法
CN104047597B (zh) 含油气泥页岩地层测井曲线标准化方法
CN1291241C (zh) 等速砂泥岩地层的储层预测油气的方法
CN100351650C (zh) 一种利用叠前地震波形反演构建虚拟井数据的方法
CN1797037A (zh) 一种地震波波阻抗反演的方法
CN107526101A (zh) 一种获取地震反射波的采集和处理方法
CN109799540B (zh) 基于地质信息约束下的火山岩型铀矿床磁化率反演方法
CN1877366A (zh) 重磁延拓回返垂直导数目标优化处理技术
CN102749643A (zh) 一种波动方程正演的瑞利面波频散响应计算方法及其装置
CN107515290A (zh) 岩石矿物组分含量定量计算方法
WO2023000257A1 (zh) 一种砂岩型铀矿成矿有利部位地质-地震三维预测方法
CN105738952A (zh) 一种水平井区储层岩石相建模方法
CN105301647A (zh) 区分灰质泥岩和砂岩的方法
Elkarmoty et al. Deterministic three-dimensional rock mass fracture modeling from geo-radar survey: A case study in a sandstone quarry in Italy
Shan et al. Spectral decomposition and a waveform cluster to characterize strongly heterogeneous paleokarst reservoirs in the Tarim basin, China
CN1188712C (zh) 利用振电信号结合地震波进行勘探和测试
CN1869734A (zh) 重力勘探数据处理变密度地形校正方法
Erokhin et al. CSP-method prospecting of fracture-cavernous reservoirs in the Bazhen formation of the Salym oilfield
CN104375187B (zh) 一种基于微测井和非地震资料设计大炮井深的方法
CN103278852B (zh) 利用地震数据体波形结构特征模型预测油气的方法
CN109669211A (zh) 海底岩性薄砂体储层预测方法和装置、存储介质及终端
Schwenk Constrained parameterization of the multichannel analysis of surface waves approach with application at Yuma Proving Ground, Arizona
CN102147480A (zh) 一种三维观测系统测试方法
Stucchi et al. High‐resolution reflection seismic survey at the Patigno landslide, Northern Apennines, Italy

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
ASS Succession or assignment of patent right

Owner name: CHINA PETROLEUM +. NATURAL GAS CO., LTD.

Free format text: FORMER OWNER: YANG HUI

Effective date: 20110620

C41 Transfer of patent application or patent right or utility model
COR Change of bibliographic data

Free format text: CORRECT: ADDRESS; FROM: 100083 NO. 20, XUEYUAN ROAD, HAIDIAN DISTRICT, BEIJING TO: 100007 CHINA PETROLEUM BUILDING, NO. 9, DONGZHIMEN NORTH STREET, DONGCHENG DISTRICT, BEIJING

TR01 Transfer of patent right

Effective date of registration: 20110620

Address after: 100007 Dongcheng District, Dongzhimen, China, North Street, No. 9 Oil Mansion, No.

Patentee after: China Petroleum & Natural Gas Co., Ltd.

Address before: 100083 Haidian District, Xueyuan Road, No. 20,

Patentee before: Yang Hui

CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20081119

Termination date: 20160712

CF01 Termination of patent right due to non-payment of annual fee