[发明专利]一种航空-地面-井中磁异常数据协同三维反演方法

权利要求书

1.一种空中-地面-井中磁异常数据协同三维反演方法，其特征在于：包括以下步骤：

a、利用航空化极磁数据d_a建立水平加权函数矩阵W_h，其对角线元素表达式为：

航空磁数据采样点个数为N₁，剖分的模型块体个数为M，W_h为M×M维矩阵，d_a(x,y,z₀)代表坐标为(x,y,z₀)处的航空化极磁数据，航测面高度h＝z₀；τ₁为水平加权强度因子，τ₁越大，则水平加权强度越大；通过上式可以看出，同一水平位置的块体的w_h数值相同，不随深度变化，因此，称其为水平加权函数；

b、利用井中化极磁数据d_w建立垂直加权函数矩阵W_z，其对角线元素表达式为：

井中磁数据采样点个数为N₂，剖分的模型块体个数为M，W_z为M×M维矩阵，Mo_w(x₀,y₀,z)代表坐标为(x,y,z₀)处的井中化极磁数据，井位坐标为(x₀,y₀)。τ₂为垂直加权强度因子，τ₂越大，则垂直加权强度越大；通过上式可以看出，同一深度位置的块体的w_z数值相同，与水平位置(x,y)无关，因此，称其为垂直加权函数；

c、综合以上水平加权函数和垂直加权函数，构建组合模型加权函数矩阵：

W_m＝W_hW_z

d、将地面磁数据d作为观测数据，引入组合模型加权函数W_m，建立目标函数：

其中，A为地面磁数据的正演算子矩阵，m为反演的模型磁化强度参数，d为地面磁数据，α为正则化因子，W_m为由航空、井中磁数据构建的组合模型加权函数矩阵；利用最优化算法可求解目标函数最优解，即反演得到模型磁化强度结果。