一种基于海洋磁力异常小波变换的居里面深度计算方法技术

本发明专利技术公开了一种基于海洋磁力异常小波变换的居里面深度计算方法，首先提取并规则化实测海洋磁力异常数据；然后对磁异常数据进行小波变换；再通过功率谱的计算获得磁性体的的顶深和平均中心深度；最后根据顶深和平均中心深度估算居里面的深度。利用磁力异常小波4阶变换，实现压制浅部异常，突出深部异常功能，进一步利用功谱法计算磁性体的顶深和平均中心深度，最终获得居里面深度，在一定程度上提高了居里面深度计算的准确性。

【技术实现步骤摘要】
一种基于海洋磁力异常小波变换的居里面深度计算方法
本专利技术涉及海洋磁力勘探
，更具体的说是涉及一种基于海洋磁力异常小波变换的居里面深度计算方法。
技术介绍
居里面也称为居里等温面，是指地壳岩石中铁磁性矿物因温度达到居里面而变为顺磁性界面，一般可将顺磁性理解为无磁性，则居里面可视为岩石圈磁性层的界面。居里面深度的计算方法目前主要有两种方法：地热法和谱分析法。地热法通过大地热流值、地温梯度等进行估算，由于难以获取研究区域的足够的地热数据，因此地热法在实际应用中并不采用，更多的是通过有限的地热点对利用谱分析法求得的居里面深度进行校正比较。谱分析方法通过利用磁异常数据进行反演计算。主要有两种：谱峰法和矩谱法。谱峰法是最早被SpectorandGrant在1970年提出的基于统计理论模型的“等效论”磁异常数据处理方法，他们指出磁异常的对数径向平均傅里叶谱与磁性层顶深有关，并且谱在频率域波数轴会出现峰位，这与磁性层的厚度有关。BhattacharyyaandLeu(1975a；1975b；1977)在谱峰法的基础上，提出了矩谱法直接计算磁性本文档来自技高网...

【技术保护点】
1.一种基于海洋磁异常小波变换的居里面深度计算方法，其特征在于，包括以下步骤：/n步骤1：在实测海洋磁力数据中提取实测海洋磁力异常数据，并进行规则化，获得磁异常数据；/n步骤2：对所述磁异常数据进行小波变换，获得小波域磁异常数据；/n步骤3：根据所述小波域磁异常数据计算磁性体的顶深z

【技术特征摘要】
1.一种基于海洋磁异常小波变换的居里面深度计算方法，其特征在于，包括以下步骤：
步骤1：在实测海洋磁力数据中提取实测海洋磁力异常数据，并进行规则化，获得磁异常数据；
步骤2：对所述磁异常数据进行小波变换，获得小波域磁异常数据；
步骤3：根据所述小波域磁异常数据计算磁性体的顶深zt；
步骤4：计算所述磁性体的平均中心深度z0；
步骤5：根据所述顶深zt和所述平均中心深度z0计算居里面深度：
hc＝zb＝2z0-zt；
其中hc为居里面深度，zb为磁性体底部深度。


2.根据权利要求1所述的一种基于海洋磁异常小波变换的居里面深度计算方法，其特征在于，所述步骤1中提取所述实测海洋磁力异常数据和进行所述规则化的公式为：
[x,y]＝meshgrid(xx,yy)；
△T(x,y)＝griddata(xx,yy,△T,x,y)；
其中，xx和yy分别为原始测量点经度和纬度；x和y分别为经网格化处理的经度和纬度；△T为获得的所述磁异常数据；meshgrid为对xx和yy进行网格化处理；griddata是依据测量值对网格化后的经纬度进行插值处理。


3.根据权利要求2所述的一种基于海洋磁异常小波变换的居里面深度计算方法，其特征在于，所述步骤2中进行所述小波变换的公式为：



其中分别为空间水平、垂直和对角线方向的高频细节部分；A0为小波变换系数；An为小波分解的主系数；f(x，y)为磁异常时间域表达；n为小波分解的阶数；h、m、d分别为空间水平、垂直和对角线方向的小波分解阶数。


4.根据权利要求3所述的一种基于海洋磁异常小波变换的居里面深度计算方法，其特征在于，所述步骤3中计算所述顶深zt的具体过程为：
步骤31：计算所述小...

【专利技术属性】
技术研发人员：田先德，孔敏，孙杨，刘志杰，王风帆，舒雨婷，余佳，
申请(专利权)人：国家海洋信息中心，
类型：发明
国别省市：天津;12