[发明专利]一种重磁数据处理方法及系统

摘要

本发明公开了一种重磁数据处理方法及系统。所述方法及系统将重磁数据的去噪、填充和扩边这三个问题统一考虑，基于凸集投影方法同时实现了重磁数据的去噪、填充和扩边处理。所述方法及系统首先对原始重磁异常数据补零，并通过计算归一化径向平均功率谱来确定迭代过程最终的截止波数，然后对重磁数据进行去噪、填充和扩边迭代直至达到预定的迭代次数，获得处理后的重磁数据。本发明提供的方法及系统原理简单、操作方便、数据处理效率高，同时具有较高的去噪、填充精度，且填充和扩边的衔接处光滑无畸变，非常适合重磁数据的同时去噪和扩充。

主权项

1.一种重磁数据处理方法，其特征在于，所述方法包括：获取重磁异常数据；将所述重磁异常数据中有数据缺失的位置补零，获得补零后重磁位场数据；计算所述补零后重磁位场数据的归一化径向平均功率谱；所述计算所述补零后重磁位场数据的归一化径向平均功率谱，具体包括：采用快速傅里叶变换方法求解所述补零后重磁位场数据的重磁位场数据功率谱；以所述重磁位场数据功率谱的中心为圆心，以基频的整数倍为半径作圆，生成多个环带；每个所述环带对应一个不同的频率；求解各个所述环带内所有点的功率谱值的平均值，获得多个平均功率谱值；以径向波数为横坐标，以各个所述环带内的所述平均功率谱值的对数为纵坐标，生成径向平均功率谱；对所述径向平均功率谱作归一化处理，生成所述重磁位场数据的归一化径向平均功率谱；根据所述归一化径向平均功率谱确定截止波数；所述根据所述归一化径向平均功率谱确定截止波数，具体包括：根据所述归一化径向平均功率谱的曲线特征，选取所述归一化径向平均功率谱下降预设比例a所对应的波数为截止波数c_K；根据所述截止波数对所述补零后重磁位场数据进行去噪、插值和扩边处理，获得处理后的重磁数据；所述根据所述截止波数对所述补零后重磁位场数据进行去噪、插值和扩边处理，获得处理后的重磁数据，具体包括：对所述补零后重磁位场数据作二维傅里叶变换，获得所述补零后重磁位场数据的频谱G(u,v)；根据所述截止波数确定低通滤波函数其中c_k表示第k次迭代的截止波数，1≤k≤K且满足c₁≤c₂≤…≤c_K和1＜c_k≤min(M,N)，M和N是插值扩边后重磁位场数据的尺寸；D(u,v)＝[(u‑M/2)²+(v‑N/2)²]^1/2；采用所述低通滤波函数T_k(u,v)对所述频谱G(u,v)进行低通滤波处理，获得低通滤波后留下的谱分量；对所述谱分量作二维傅里叶反变换，获得反变换后的重磁位场数据；利用所述反变换后的重磁位场数据重置所述补零后重磁位场数据，获得重置后的重磁位场数据；判断迭代次数k是否不大于预设迭代次数K，获得第一判断结果；若所述第一判断结果为所述迭代次数k不大于所述预设迭代次数K，令k＝k+1，返回所述根据所述截止波数确定低通滤波函数的步骤；若所述第一判断结果为所述迭代次数k大于所述预设迭代次数K，输出所述重置后的重磁位场数据作为处理后的重磁数据。