[发明专利]一种基于数据空间快速计算的联合反演方法

说明书

本发明公开了一种基于数据空间快速计算的联合反演方法，包括如下步骤：S1：建立重磁电震联合反演目标函数；S2：建立电阻率、速度、密度和磁化率的初始模型；S3：正演计算和雅克比矩阵的求取；S4：计算交叉梯度偏导数和拉格朗日算子；S5：计算获得模型空间迭代模型；S6：计算获得数据空间迭代模型；S7：对获得的数据空间迭代模型进行正演计算，并与观测数据进行拟合差求取，如果拟合差增大、达到最大迭代次数或者拟合差达到阈值，则停止循环反演，输出该数据空间联合反演模型，否则返回S3。该基于数据空间快速计算的联合反演方法计算速度快、占用内存小。

技术领域

本发明涉及多物性联合反演的处理技术领域，特别提供了一种基于数据空间快速计算的联合反演方法。

背景技术

综合地球物理方法已经广泛应用于深部地质探测和矿产、水、石油、地热等资源勘查。而联合反演是其综合地球物理方法数据处理的主流方法，国外科学家已经提出了二种或多种物性参数之间结构耦合的联合反演方法，国内也实现了二种或多种物性参数之间的结构耦合联合反演方法。

目前，对开展大数据量的联合反演数据处理已经成为重要的发展方向。现有的联合反演多采用在模型空间下进行反演计算，但是当反演物性参数种类从二个增加到多个或者反演区域网格剖分数据量过大时，对计算机的计算速度和内存占用都将成为考验，将会严重制约了综合地球物理方法的发展。

发明内容

鉴于此，本发明的目的在于提供一种基于数据空间快速计算的联合反演方法，以解决现有的联合反演方法计算速度慢、占用内存大等问题。

本发明提供的技术方案是：一种基于数据空间快速计算的联合反演方法，包括如下步骤：

S1：建立重磁电震联合反演目标函数，所述联合反演目标函数如式(1)：

约束条件：τ(m)＝0

其中：

C_d＝diag[C_d1,C_d2,C_d3,C_d4]，C_m＝diag[C_m1,C_m2,C_m3,C_m4]，

α＝[α₁,α₂,α₃,α₄]，