[发明专利]地球物理数据处理的方法和装置在审
| 申请号： | 201811088255.2 | 申请日： | 2018-09-18 |
| 公开（公告）号： | CN109188535A | 公开（公告）日： | 2019-01-11 |
| 发明（设计）人： | 徐亚;南方舟;黄松;胡祖志;杨辉;郝天珧 | 申请（专利权）人： | 中国科学院地质与地球物理研究所 |
| 主分类号： | G01V3/08 | 分类号： | G01V3/08;G01V3/38 |
| 代理公司： | 北京集佳知识产权代理有限公司 11227 | 代理人： | 李慧引;王宝筠 |
| 地址： | 100029 *** | 国省代码： | 北京;11 |
| 权利要求书： | 查看更多 | 说明书： | 查看更多 |
| 摘要： | |||
| 搜索关键词： | 傅里叶变换系数 地球物理数据 方法和装置 物理数据 原始地球 地球物理 位置坐标 数据处理 非规则 测点 重构 傅里叶反变换 傅里叶变换 有效信息 重构的 采样 存储 分析 | ||
1.一种地球物理数据处理的方法，其特征在于，包括以下步骤：
获取多个测点的位置坐标和原始地球物理数据；
对所述多个测点的位置坐标和所述原始地球物理数据进行非规则傅里叶变换，获得傅里叶变换系数；
选取并存储所述获得的傅里叶变换系数中的特定傅里叶变换系数；
对所述特定傅里叶变换系数进行傅里叶反变换，获得重构的地球物理数据函数。
2.根据权利要求1所述的方法，其特征在于，所述对所述多个测点的位置坐标和所述原始地球物理数据进行非规则傅里叶变换，获得傅里叶变换系数，包括：
将所述多个测点的位置坐标和所述原始地球物理数据代入非规则傅里叶变换公式，获得傅里叶变换系数；
其中，所述非规则傅里叶变换公式为：
式中：xj为测点的位置坐标，k为傅里叶变换系数F(k)所对应的频率，f(xj)表示在测点xj处的地球物理数据，N表示测点的个数，Δxj为权重。
3.根据权利要求1所述的方法，其特征在于，所述选取并存储所述获得的傅里叶变换系数中的特定傅里叶变换系数，包括：
选取并存储所述获得的傅里叶变换系数中能量最大的傅里叶变换系数[F(k)]n，作为所述特定傅里叶变换系数；
其中，所述[F(k)]n表示第n个存储的所述能量最大的傅里叶变换系数。
4.根据权利要求1所述的方法，其特征在于，所述对所述特定傅里叶变换系数进行傅里叶反变换，获得重构的地球物理数据函数，包括：
利用傅里叶反变换公式对所述特定傅里叶变换系数[F(k)]l进行傅里叶反变换，获得重构的地球物理数据函数g(xj)；
其中，所述傅里叶反变换公式为：其中l∈n，xj为测点的位置坐标，k为傅里叶变换系数F(k)所对应的频率，n为存储的所述特定傅里叶变换系数的个数。
5.根据权利要求4所述的方法，其特征在于，所述利用傅里叶反变换公式对所述特定傅里叶变换系数[F(k)]l进行傅里叶反变换，获得重构的地球物理数据函数g(xj)之后，还包括：
计算所述原始地球物理数据f(xj)与所述重构地球物理数据g(xj)的残差：δf(xj)＝f(xj)-g(xj)；
判断||δf(xj)||≤ε是否成立，所述ε为误差阈值；
如果判断出||δf(xj)||≤ε不成立，则将所述残差作为当前的地球物理数据f(xj)＝δf(xj)，并将所述当前的地球物理数据替换所述原始地球物理数据返回执行对所述多个测点的位置坐标和原始地球物理数据进行非规则傅里叶变换，获得傅里叶变换系数步骤。
6.根据权利要求2至5任一项所述的方法，其特征在于，所述原始地球物理数据，包括：一维地球物理数据，其中，所述非规则傅里叶变换公式的权重Δxj利用Δxj＝xj-xj-1计算得到；
或者，所述原始地球物理数据，包括：二维地球物理数据，其中，所述非规则傅里叶变换公式的权重Δxj利用计算得到；S(x)为采样矩阵，G(x)为高斯函数，m表示所述地球物理数据的维度。
该专利技术资料仅供研究查看技术是否侵权等信息，商用须获得专利权人授权。该专利全部权利属于中国科学院地质与地球物理研究所，未经中国科学院地质与地球物理研究所许可，擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作，请联系【客服】
本文链接：http://www.vipzhuanli.com/pat/books/201811088255.2/1.html，转载请声明来源钻瓜专利网。
