【摘要】 一种地球物理勘探中的重磁数据处理方法，其包括步骤：在波数域将重磁资料频谱分别与任意方向的n阶导数或积分频率响应函数、延拓回返响应函数及差分法下延算子频率响应函数相乘，得到下延后的频谱；再经过快速傅氏反变换，获得向下延拓后的三维重磁数据体。利用向下延拓之后的三维重磁数据体结合地震、地质、钻井资料进行对比、分析并提交下一步勘探的有利区带及目标。该方法不仅有更大的向下延拓深度并且有较好的抗噪能力，从而获得稳定的下半空间三维重磁资料分布图像。 【专利类型】发明授权 【申请人】杨辉 【申请人类型】个人 【申请人地址】100083 北京市海淀区学院路20号中国石油勘探开发研究院物探技术研究所 【申请人地区】中国 【申请人城市】北京市 【申请人区县】海淀区 【申请号】CN200810110500.5 【申请日】2008-06-13 【申请年份】2008 【公开公告号】CN101285896B 【公开公告日】2010-07-21 【公开公告年份】2010 【授权公告号】CN101285896B 【授权公告日】2010-07-21 【授权公告年份】2010.0 【IPC分类号】G01V7/06 【发明人】杨辉 【主权项内容】一种地球物理勘探中的重磁数据处理方法，该方法包括以下步骤：(1)首先利用重力仪或磁力仪在野外采集重磁资料，或者在室内对原有的重磁异常图进行数字化而得到重磁资料，然后对野外采集或室内数字化的重磁资料进行各种常规校正及改正；(2)利用Surfer软件中的Kriging插值方法将重磁资料网格化，形成规则网重磁资料fa(x，y)，其中x、y分别为重磁资料的平面网格点坐标；(3)进行镶边处理：利用低精度的重磁资料进行补充镶边；将规则网重磁资料fa(x，y)转变为其线数和点数扩大为2的整数幂的重磁资料fb(x，y)；(4)利用快速傅氏正变换(FFT)将镶边后的重磁资料fb(x，y)变换到波数域，获得fb(x，y)的频谱Fb(u，v)，Fb(u，v)＝FFT[fb(x，y)]，其中u、v分别为x和y方向的波数；(5)在波数域，将频谱Fb(u，v)与重磁场任意方向的n阶导数或积分频率响应函数相乘，增强局部或区域重磁异常，即：Fc(u，v)＝Fb(u，v)×qn(u，v)其中n为有理数，i为虚数；α、β、γ分别为所求导数或积分方向的三个方向余弦；n为导数或积分的阶数，即分辨率参数；(6)在波数域，根据向下延拓深度h，将增强后的局部或区域重磁异常Fc(u，v)与延拓回返频率响应函数相乘，抑制高频成分，使其满足向下延拓的条件，即：其中h为向下延拓深度，p为延拓回返次数；h、p为正有理数；(7)在波数域，在抑制高频成分的基础上，与差分法向下延拓频率响应函数相乘，得到向下延拓后的频谱Fe(u，v)Fe(u，v)＝Fd(u，v)×d(u，v)；(8)再经过快速傅氏反变换(IFFT)，获得向下延拓后的重磁资料fe(x，y)，fe(x，y)＝IFFT[Fe(u，v)]；(9)进行反镶边处理，即从fe(x，y)中提取原数据大小的向下延拓深度h后的重磁数据ff(x，y)；(10)重复步骤(6)，从而获得下半空间三维重磁数据体fg(x，y，z)，其中，x，y，z为重磁异常的空间网格点坐标；(11)利用向下延拓之后的三维重磁数据体结合地震、地质、钻井资料进行对比、分析并提交下一步勘探的有利区带及目标。FSB00000043037600011.tif, FSB00000043037600012.tif, FSB00000043037600021.tif, FSB00000043037600022.tif 【当前权利人】中国石油天然气股份有限公司 【当前专利权人地址】北京市东城区东直门北大街9号中国石油大厦 【引证次数】2.0 【自引次数】2.0 【家族引证次数】2.0 【家族被引证次数】21
