[发明专利]物探中的椭圆窗方向带通保幅滤波数据处理方法有效

专利信息
申请号: 201010516723.9 申请日: 2010-10-22
公开(公告)号: CN102338884A 公开(公告)日: 2012-02-01
发明(设计)人: 杨辉 申请(专利权)人: 中国石油天然气股份有限公司
主分类号: G01V1/28 分类号: G01V1/28
代理公司: 北京市中实友知识产权代理有限责任公司 11013 代理人: 刘天语;武尚修
地址: 100007 北京市*** 国省代码: 北京;11
权利要求书: 查看更多 说明书: 查看更多
摘要: 发明是物探技术中椭圆窗方向带通保幅滤波数据处理方法,采用插值方法将重磁资料网格化,利用低精度的重磁资料补充镶边,快速傅氏正变换,旋转坐标系使滤波器具有方向性,计算椭圆窗函数使得滤波器在波数域坐标系为椭圆窗,计算初始的汉宁窗滤波器并对初始的汉宁窗滤波器进行补偿:在波数域将重磁资料的频谱与带通滤波频率响应函数相乘,得到滤波后的波数域结果,傅氏反变换,反镶边处理,得到带通滤波后的重磁数据。本发明滤波器在通带内信号保幅,在阻带内信号完全抑制,能够压制重磁异常低频背景和高频噪声,突出中频有用信号,特别适合突出在波数域为椭圆窗的中频有用信号。
搜索关键词: 物探 中的 椭圆 方向 带通保幅 滤波 数据处理 方法
【主权项】:
1.一种物探中的椭圆窗方向带通保幅滤波数据处理方法,通过如下步骤实现:1)在工区利用重力仪或磁力仪在野外采集重力或磁力资料;2)对野外采集的重磁资料进行校正及改正;3)采用插值方法将重磁资料网格化,形成规则网重磁资料fa(x,y),其中x、y分别为重磁资料的平面网格点坐标;4)利用低精度的重磁资料补充镶边;将规则网重磁资料fa(x,y)转变为其线数和点数扩大为2的整数幂重磁资料fb(x,y);5)利用快速傅氏正变换FFT将镶边后的重磁资料fb(x,y)变换到波数域,获得fb(x,y)的频谱Fb(u,v),Fb(u,v)=FFT[fb(x,y)],其中u、v分别为x和y方向的波数;6)用下列公式进行波数域坐标旋转,从波数域uov坐标系旋转到新的波数域u′Ov′坐标系,从而使滤波器具有方向性:u=ucosα-vsinαv=usinα+vcosα]]>其中:u′、v′分别为x′、y′方向的波数;α为u轴与u′轴的夹角或坐标旋转角度;7)用下列公式计算椭圆窗函数R,使得滤波器在波数域u′Ov′坐标系为椭圆窗:R=[(|u|-u0)2Δu02+(|v|-v0)2Δv02]1/2·π]]>其中:u′0、v′0分别为中心波数,是波数域椭圆中心的位置坐标,Δu′0、Δv′0分别为通频带半宽度,是波数域椭圆的半轴长度,且u′0≥O,v′0≥0,Δu′0>0,Δv′0>0;8)用下列公式计算初始的汉宁窗滤波器:9)用下列公式对初始的汉宁窗滤波器进行补偿:其中:p为补偿次数,正整数;10)在波数域,将频谱Fb(u,v)与带通滤波频率响应函数相乘,得到滤波后的波数域的滤波结果Fc(u,v);11)经过快速傅氏反变换IFFT,获得椭圆窗带通滤波后的重磁资料fc(x,y),fc(x,y)=IFFT[Fc(u,v)];12)反镶边处理,从fc(x,y)中得到原数据大小的带通滤波后的重磁数据fd(x,y)。
下载完整专利技术内容需要扣除积分,VIP会员可以免费下载。

该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于中国石油天然气股份有限公司,未经中国石油天然气股份有限公司许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服

本文链接:http://www.vipzhuanli.com/patent/201010516723.9/,转载请声明来源钻瓜专利网。

×

专利文献下载

说明:

1、专利原文基于中国国家知识产权局专利说明书;

2、支持发明专利 、实用新型专利、外观设计专利(升级中);

3、专利数据每周两次同步更新,支持Adobe PDF格式;

4、内容包括专利技术的结构示意图流程工艺图技术构造图

5、已全新升级为极速版,下载速度显著提升!欢迎使用!

请您登陆后,进行下载,点击【登陆】 【注册】

关于我们 寻求报道 投稿须知 广告合作 版权声明 网站地图 友情链接 企业标识 联系我们

钻瓜专利网在线咨询

周一至周五 9:00-18:00

咨询在线客服咨询在线客服
tel code back_top