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

权利要求书

1.一种磁化率反演方法，其特征在于，包括：

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

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

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

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

2.根据权利要求1所述的磁化率反演方法，其特征在于，所述S3包括：

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

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

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

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

3.根据权利要求1所述的磁化率反演方法，其特征在于，所述S4包括：

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

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

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

4.根据权利要求2所述的磁化率反演方法，其特征在于，所述形函数为线性插值形函数和二次插值形函数中任意一者。

5.一种磁化率反演系统，其特征在于，包括：

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

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

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

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

6.根据权利要求5所述的磁化率反演系统，其特征在于，所述模型扰动模块包括：

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

正则化单元，用于对所述线性方程组进行正则化处理得到正则化方程组；

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

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

7.根据权利要求5所述的磁化率反演系统，其特征在于，所述磁化率反演模块包括：

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

判断单元，用于基于更新后的模型矢量得到平均均方误差，并对所述平均均方误差进行判断；当所述平均均方误差大于误差设定值时，返回至所述空间域残差模块；

磁化率反演单元，用于当所述平均均方误差小于或等于误差设定值时，基于所述空间域模型扰动得到磁化率反演模型。

8.根据权利要求6所述的磁化率反演系统，其特征在于，所述形函数为线性插值形函数和二次插值形函数中任意一者。