[发明专利]一种磁化率反演方法及系统

说明书

本发明涉及一种磁化率反演方法及系统，所述方法包括：S1，基于观测磁场强度和理论磁场强度得到空间域残差；S2，基于波数域残差通过二维离散傅里叶反变换表示所述空间域残差；S3，基于所述空间域残差得到波数域模型扰动，进一步基于傅里叶反变换得到空间域模型扰动；S4，基于所述空间域模型扰动得到磁化率反演模型。本发明降低了所需内存空间，提高了计算效率，进一步提高了反演效率。

技术领域

本发明涉及地球物理技术领域，特别是涉及一种磁化率反演方法及系统。

背景技术

磁测数据在地球物理勘探领域中有着广泛的应用。例如，在矿物勘探中,磁方法可以用于决定体参数；在石油勘探中，磁方法可以用来映射地下的沉积特性和缺点，用来控制沉积盆地的沉积背景。

物理参数的反演，例如计算磁化率和磁化强度，是使用磁场数据的主要科学问题。实际应用中，能够反演目标地理区域的磁化率，得到该目标地理区域的磁化率模型，通过该磁化率模型能够反映该目标地理区域内各个地理位置的磁化率，进而将磁化率符合预设要求的某地理位置选择出来，作为进一步的研究对象。

现有磁化率反演技术中，多采用广义逆迭代反演优化方法和阻尼广义逆迭代反演优化方法，但是现有反演技术都存在一定的缺点，包括：1)由于计算区域剖分单元的数量巨大，耗费大量的内存，且计算效率较低；2)测点必须有规则分布，不利于野外施工；3)反演方程组求解复杂，且收敛较慢；4)不能控制参数提高反演的分辨率；5)由于算法本身的特点，不利于将复杂的大型反演问题分解为有限多个小问题进行并行求解。

发明内容

有鉴于此，本发明的目的是提供一种磁化率反演方法及系统，以解决现有技术存在的问题。

为实现上述目的，本发明提供了如下方案：

一种磁化率反演方法，包括：

S1，基于观测磁场强度和理论磁场强度得到空间域残差；

S2，基于波数域残差通过二维离散傅里叶反变换表示所述空间域残差；

S3，基于所述空间域残差得到波数域模型扰动，进一步基于傅里叶反变换得到空间域模型扰动；

S4，基于所述空间域模型扰动得到磁化率反演模型。

优选地，所述S3包括：

S31，基于形函数对所述空间域残差积分，得到线性方程组；

S32，对所述线性方程组进行正则化处理得到正则化方程组；

S33，基于所述正则化方程组进行求解得到所述波数域模型扰动；

S34，对所述波数域模型扰动做傅里叶反变换得到所述空间域模型扰动。

优选地，所述S4包括：

S41，基于所述空间域模型扰动对模型矢量进行更新得到更新后的模型矢量；

S42，基于更新后的模型矢量得到平均均方误差，并对所述平均均方误差进行判断；当所述平均均方误差大于误差设定值时，返回至“S1”；当所述平均均方误差小于或等于误差设定值时，执行“S43”；

S43，基于所述空间域模型扰动得到磁化率反演模型。

优选地，所述形函数为线性插值形函数和二次插值形函数中任意一者。

本发明还提供了一种磁化率反演系统，包括：

空间域残差模块，用于基于观测磁场强度和理论磁场强度得到空间域残差；

波数域残差模块，用于基于波数域残差通过二维离散傅里叶反变换表示所述空间域残差；