首页 / 专利库 / 信号处理 / 稀疏变换 / 小波变换 / 离散小波反变换 / 一种1D位场异常曲线构建2D位场异常剖面的方法

一种1D位场异常曲线构建2D位场异常剖面的方法

阅读:1020发布:2020-07-07

专利汇可以提供一种1D位场异常曲线构建2D位场异常剖面的方法专利检索,专利查询,专利分析的服务。并且一种1D位场异常曲线构建2D位场异常剖面的方法,步骤为:第一步,将含有m个 采样 点,点距为△x的一维位场异常数据读取到一维数组f中;第二步,对f进行一维离散小波多尺度分解;第三步,分别提取每一尺度小波细节的高频分量,并进行一维 离散小波反变换 ;第四步,将尺度维等效为拟深度维,并按由小尺度到大尺度的顺序组建多尺度高频分量的位场异常数据,形成二维数据F;第五步,对二维数据F进行网格化及插值处理,然后将其绘制成二维位场异常剖面图,便实现了一维位场异常曲线构建二维位场异常剖面,获得的二维剖面具有异常信息展示直观、易于地质解释的优点。,下面是一种1D位场异常曲线构建2D位场异常剖面的方法专利的具体信息内容。

1.一种1D位场异常曲线构建2D位场异常剖面的方法,其特征在于,包括如下步骤:
第一步,将含有m个采样点,点距为△x的一维位场异常数据读取到一维数组f中;
第二步,对一维数组f进行一维离散小波多尺度分解正变换,一维离散小波多尺度分解正变换为:
式中,ψ称为基本小波;ψ*为ψ的共轭函数;a表示伸缩系数,反映特定基函数的宽度;b表示平移系数,指定沿x轴平移的位置,且 并令a0=2,b0=1,j∈Z,k∈Z;△x为点距;i=0,1,…,m-1,为采样序号;
第三步,分别提取每一尺度小波细节的高频分量,并进行一维离散小波多尺度分解反变换,一维离散小波多尺度分解反变换为:
式中,ψ称为基本小波;j∈Z,k∈Z;x=i·△x,为数据采样点坐标;△x为点距;i=0,
1,…,m-1,为采样序号;
第四步,将尺度维等效为拟深度维,并按由小尺度到大尺度的顺序组建多尺度高频分量的位场异常数据,形成二维数据F;
第五步,对二维数据F进行网格化及插值处理,然后将其绘制成二维位场异常剖面图,便实现了一维位场异常曲线构建二维位场异常剖面。
2.根据权利要求1所述的一种1D位场异常曲线构建2D位场异常剖面的方法,其特征在于,包括如下步骤:
第一步,将含有155个采样点,点距为2km的一维布格重异常数据读取到一维数组f中;
第二步,对一维数组f进行一维离散小波13尺度分解正变换,一维离散小波13尺度分解正变换为:
式中,ψ称为基本小波;ψ*为ψ的共轭函数;a表示伸缩系数,反映特定基函数的宽度;b表示平移系数,指定沿x轴平移的位置,且 并令a0=2,b0=1,j=0,1,…,12,k∈Z;△x=2km,为点距;i=0,1,…,154,为采样序号;
第三步,分别提取每一尺度小波细节的高频分量,并进行一维离散小波多尺度分解反变换,一维离散小波多尺度分解反变换为:
式中,ψ称为基本小波;j=0,1,…,12,k∈Z;x=i·△x,为数据采样点坐标;△x=2km,为点距;i=0,1,…,154,为采样序号;
第四步,将尺度维等效为拟深度维,并按由小尺度到大尺度的顺序组建13尺度高频分量的布格重力异常数据,形成二维数据F;
第五步,对二维数据F进行网格化及插值处理,然后将其绘制成二维布格重力异常剖面图,便实现了一维布格重力异常曲线构建二维布格重力异常剖面。

说明书全文

一种1D位场异常曲线构建2D位场异常剖面的方法

技术领域

[0001] 本发明属于位场数据处理技术领域,特别涉及一种1D位场异常曲线构建2D位场异常剖面的方法。

背景技术

[0002] 位场包括重场、磁场电场。位场的异常是来自不同深度、不同尺度场源体综合响应及共同叠加的结果。受位场勘探方法理论自身的限制,野外位场数据采集难以获得深度(维)信息,致使位场异常的地质解释不够直观,从而增大了位场资料地质解释工作的难度,并降低了地质解释的准确性。

发明内容

[0003] 为了克服上述现有技术的不足,本发明的目的在于提出了一种1D位场异常曲线构建2D位场异常剖面的方法,该方法具有异常信息展示直观,易于地质解释的优点。
[0004] 为了实现上述目的,本发明采用的技术方案是:一种1D位场异常曲线构建2D位场异常剖面的方法,包括如下步骤:
[0005] 第一步,将含有m个采样点,点距为△x的一维位场异常数据读取到一维数组f中;
[0006] 第二步,对一维数组f进行一维离散小波多尺度分解正变换,一维离散小波多尺度分解正变换为:
[0007]
[0008] 式中,ψ称为基本小波(也称为母小波);ψ*为ψ的共轭函数;a表示伸缩系数,反映特定基函数的宽度(也叫做尺度);b表示平移系数,指定沿x轴平移的位置,且并令a0=2,b0=1,j∈Z,k∈Z;△x为点距;i=0,2,…,m-1,为采样序号;
[0009] 第三步,分别提取每一尺度小波细节的高频分量,并进行一维离散小波多尺度分解反变换,一维离散小波多尺度分解反变换为:
[0010]
[0011] 式中,ψ称为基本小波(也称为母小波);j∈Z,k∈Z;x=i·△x,为数据采样点坐标,△x为点距;i=0,2,…,m-1,为采样序号;
[0012] 第四步,将尺度维等效为拟深度维,并按由小尺度到大尺度的顺序组建多尺度高频分量的位场异常数据,形成二维数据F;
[0013] 第五步,对二维数据F进行网格化及插值处理,然后将其绘制成二维位场异常剖面图,便实现了一维(1D)位场异常曲线构建二维(2D)位场异常剖面。
[0014] 本发明的有益效果是:
[0015] 本方法基于浅层场源体产生的位场异常频率相对较高,而深层场源体产生的位场异常频率相对较低,并且随深度增加,位场异常的频率在逐渐降低为假设前提,通过一维离散小波多尺度分解技术,将来自不同深度、不同场源的位场异常进行多尺度分解,然后把尺度维等效为拟深度维,并按由小尺度到大尺度的次序将位场异常进行组合,从而将一维(1D)位场异常曲线拓展构建出二维(2D)位场异常剖面,此二维(2D)剖面可直接用于地质解释;该方法具有异常信息展示直观,易于地质解释的优点。附图说明
[0016] 图1为实际一维(1D)布格重力异常曲线图。
[0017] 图2为本发明一维离散小波13尺度分解并重构获得的不同尺度高频分量布格重力异常曲线图。
[0018] 图3为本发明一维(1D)布格重力异常曲线13尺度分解构建获得的二维(2D)布格重力异常剖面图。

具体实施方式

[0019] 下面结合附图和实施例对本发明进一步详细说明。
[0020] 参见图1、2、3,一种1D位场异常曲线构建2D位场异常剖面的方法,包括如下步骤:
[0021] 第一步,将含有m个采样点,点距为△x的一维位场异常数据读取到一维数组f中,参见图1;
[0022] 第二步,对一维数组f进行一维离散小波多尺度分解正变换,一维离散小波多尺度分解正变换为:
[0023]
[0024] 式中,ψ称为基本小波(也称为母小波);ψ*为ψ的共轭函数;a表示伸缩系数,反映特定基函数的宽度(也叫做尺度);b表示平移系数,指定沿x轴平移的位置,且并令a0=2,b0=1,j∈Z,k∈Z;△x为点距;i=0,2,…,m-1,为采样序号;
[0025] 第三步,分别提取每一尺度小波细节的高频分量,并进行一维离散小波多尺度分解反变换,一维离散小波多尺度分解反变换为:
[0026]
[0027] 式中,ψ称为基本小波(也称为母小波);j∈Z,k∈Z;x=i·△x,为数据采样点坐标,△x为点距;i=0,2,…,m-1,为采样序号,参见图2;
[0028] 第四步,将尺度维等效为拟深度维,并按由小尺度到大尺度的顺序组建多尺度高频分量的位场异常数据,形成二维数据F;
[0029] 第五步,对二维数据F进行网格化及插值处理,然后将其绘制成二维位场异常剖面图,便实现了一维(1D)位场异常曲线构建二维(2D)位场异常剖面,参见图3。
[0030] 实施实例
[0031] 将含有155个采样点,点距为2km的实测一维布格重力异常数据为例说明实施步骤:
[0032] 第一步,将含有155个采样点,点距为2km的一维布格重力异常数据读取到一维数组f中,图1;
[0033] 第二步,对一维数组f进行一维离散小波13尺度分解正变换,一维离散小波13尺度分解正变换为:
[0034]
[0035] 式中,ψ称为基本小波(也称为母小波);ψ*为ψ的共轭函数;a表示伸缩系数,反映特定基函数的宽度(也叫做尺度);b表示平移系数,指定沿x轴平移的位置,且并令a0=2,b0=1,j=0,1,…,12,k∈Z;△x=2km,为点距;i=0,1,…,154,为采样序号;
[0036] 第三步,分别提取每一尺度小波细节的高频分量,并进行一维离散小波多尺度分解反变换,一维离散小波多尺度分解反变换为:
[0037]
[0038] 式中,ψ称为基本小波(也称为母小波);j=0,1,…,12,k∈Z;x=i·△x,为数据采样点坐标;△x=2km,为点距;i=0,1,…,154,为采样序号,图2;
[0039] 第四步,将尺度维等效为拟深度维,并按由小尺度到大尺度的顺序组建13尺度高频分量的布格重力异常数据,形成二维数据F;
[0040] 第五步,对二维数据F进行网格化及插值处理,然后将其绘制成二维布格重力异常剖面图,便实现了一维(1D)布格重力异常曲线构建二维(2D)布格重力异常剖面,图3。
[0041] 实例效果说明:
[0042] 图1是实际一维(1D)布格重力异常曲线,其中横坐标为采样点序号,纵坐标为布格重力异常值△g,该曲线沿测线方向(横向)上的异常特征表现较明显,但缺乏深度域的信息,无法明确判断出由什么样的地质体、在什么位置引起了这些布格重力异常。
[0043] 图2中每一条曲线代表不同尺度的布格重力异常曲线,其中横坐标为采样点序号,纵坐标为布格重力异常值△g。
[0044] 图3是由一维(1D)布格重力异常曲线经过一维离散小波多尺度分解及重构获得的二维(2D)布格重力异常剖面,其中横坐标为采样点序号,纵坐标为尺度(拟深度),色标代表二维剖面不同位置处的布格重力异常值△g,由于本发明将尺度维等效为拟深度维,从而拓展出深度域信息,从该剖面图的异常值大小及异常轮廓可有效识别出布格重力异常在二维剖面上的分布特征,并定性判断出引起布格重力异常的地质体是高密度体还是低密度体。
高效检索全球专利

专利汇是专利免费检索,专利查询,专利分析-国家发明专利查询检索分析平台,是提供专利分析,专利查询,专利检索等数据服务功能的知识产权数据服务商。

我们的产品包含105个国家的1.26亿组数据,免费查、免费专利分析。

申请试用

分析报告

专利汇分析报告产品可以对行业情报数据进行梳理分析,涉及维度包括行业专利基本状况分析、地域分析、技术分析、发明人分析、申请人分析、专利权人分析、失效分析、核心专利分析、法律分析、研发重点分析、企业专利处境分析、技术处境分析、专利寿命分析、企业定位分析、引证分析等超过60个分析角度,系统通过AI智能系统对图表进行解读,只需1分钟,一键生成行业专利分析报告。

申请试用

QQ群二维码
意见反馈