[发明专利]物探中的椭圆窗方向带通保幅滤波数据处理方法有效
申请号: | 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 ′ = u cos α - v sin α v ′ = u sin α + v cos α ]]> 其中:u′、v′分别为x′、y′方向的波数;α为u轴与u′轴的夹角或坐标旋转角度;7)用下列公式计算椭圆窗函数R,使得滤波器在波数域u′Ov′坐标系为椭圆窗:R = [ ( | u ′ | - u ′ 0 ) 2 Δ u ′ 0 2 + ( | v ′ | - v ′ 0 ) 2 Δ v ′ 0 2 ] 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/,转载请声明来源钻瓜专利网。