[发明专利]地球物理数据处理的方法和装置

说明书

本发明公开了一种地球物理数据处理的方法和装置，通过获取多个测点的位置坐标和原始地球物理数据；对所述多个测点的位置坐标和所述原始地球物理数据进行非规则傅里叶变换，获得傅里叶变换系数；选取并存储所述获得的傅里叶变换系数中的特定傅里叶变换系数；对所述特定傅里叶变换系数进行傅里叶反变换，获得重构的地球物理数据函数。通过上述公开的方法和装置，可以对非规则采样的地球物理数据进行全局性重构，在重构规则网地球物理数据的同时避免了原始地球物理数据有效信息的丢失，以便用于后续地球物理数据分析。

技术领域

本发明涉及地理勘探技术领域，尤其涉及一种地球物理数据处理方法和装置。

背景技术

重磁勘探是指在地球上进行的，并且以地球的重力场作为被探测物体的引力场，以地球的磁场作为被探测物体的磁化场，从而用重、磁力仪器观测被探测物体受地球的吸引所产生的重力异常及其被地磁场磁化后所产生的磁异常，来达到探测目的的勘探方法。由于重磁勘探应用广泛，成本不高，因此在勘探领域一般是其他勘探方法之前的首选方法。

在重磁勘探的实践中，由于受地形等多种客观因素影响，在一些测点上不能实际测量，从而造成实测点分布不均匀，观测值并不能均匀覆盖观测区域。但是异常的重磁数据处理，要求数据均匀分布。通常通过插值的方法进行数据重构，插值的方法是基于数据的局部连续性假设及观测点相邻空间的观测值来获取均匀的网格化数据，该网格化数据用于后续的数据处理和分析。然而，用插值的方法重构规则网重磁数据的过程中，由于未考虑到数据的全局性特征，可能会造成原始观测数据的丢失。

发明内容

基于上述现有技术的不足，本发明提出了一种地球物理数据处理的方法和装置，以实现将非规则点获得的地球物理数据重构为规则网地球物理数据的同时，保持原始数据的有效信息。

为了实现上述目的，现提出的方案如下：

本发明第一方面公开了一种地球物理数据处理的方法，包括以下步骤：

获取多个测点的位置坐标和原始地球物理数据；

对所述多个测点的位置坐标和所述原始地球物理数据进行非规则傅里叶变换，获得傅里叶变换系数；

选取并存储所述获得的傅里叶变换系数中的特定傅里叶变换系数；

对所述特定傅里叶变换系数进行傅里叶反变换，获得重构的地球物理数据函数。

可选地，在上述地球物理数据处理的方法中，所述对所述多个测点的位置坐标和所述原始地球物理数据进行非规则傅里叶变换，获得傅里叶变换系数，包括：

将所述多个测点的位置坐标和所述原始地球物理数据代入非规则傅里叶变换公式，获得傅里叶变换系数；

其中，所述非规则傅里叶变换公式为：

式中：x_j为测点的位置坐标，k为傅里叶变换系数F(k)所对应的频率，f(x_j)表示在测点x_j处的地球物理数据，N表示测点的个数，Δx_j为权重。

可选地，在上述地球物理数据处理的方法中，所述选取并存储所述获得的傅里叶变换系数中的特定傅里叶变换系数，包括：

选取并存储所述获得的傅里叶变换系数中能量最大的傅里叶变换系数[F(k)]_n，作为所述特定傅里叶变换系数；

其中，所述[F(k)]_n表示第n个存储的所述能量最大的傅里叶变换系数。

可选地，在上述地球物理数据处理的方法中，所述对所述特定傅里叶变换系数进行傅里叶反变换，获得重构的地球物理数据函数，包括：