[发明专利]复杂磁性体井中磁场三分量计算方法、装置、设备及介质

权利要求书

1.复杂磁性体井中磁场三分量计算方法，其特征在于，包括：

确定目标区域，构建内部包含所述目标区域的棱柱体模型；对棱柱体模型进行网格剖分，剖分成若干个小棱柱体单元，给每个小棱柱体单元的磁化率赋值；

采用三维离散卷积法，计算各小棱柱体单元在井中观测点处产生的井中磁场三分量的各个分量；

基于各小棱柱体单元在井中观测点处产生的井中磁场三分量的各个分量，得到整个棱柱体模型在井中观测点处产生的井中磁场三分量。

2.根据权利要求1所述的复杂磁性体井中磁场三分量计算方法，其特征在于，目标区域中每个小棱柱体单元的磁化率为常值，不同小棱柱体单元的磁化率取值不同，以此刻画任意磁化率分布、复杂形状的磁性体，位于空气部分的小棱柱体单元的磁化率值设为零，以此刻画起伏地形。

3.根据权利要求1或2所述的复杂磁性体井中磁场三分量计算方法，其特征在于，所有小棱柱体单元在井中观测点(x₀,y₀,z_n)处产生的井中磁场三分量包括9个分量，分别是b_x1(x₀,y₀,z_n)、b_x2(x₀,y₀,z_n)、b_x3(x₀,y₀,z_n)、b_y1(x₀,y₀,z_n)、b_y2(x₀,y₀,z_n)、b_y3(x₀,y₀,z_n)、b_z1(x₀,y₀,z_n)、b_z2(x₀,y₀,z_n)和b_z3(x₀,y₀,z_n)，计算公式分别为：

其中各井中观测点的水平坐标为x₀，y₀，各井中观测点的垂直坐标为z_n，n＝1,2,...,N，N为井中观测点个数，p＝1,2,...,P，P为棱柱体模型x方向小棱柱体单元的剖分个数，q＝1,2,...,Q，Q为棱柱体模型y方向小棱柱体单元的剖分个数，r＝1,2,...,R，R为棱柱体模型z方向小棱柱体单元的剖分个数，m_x(ξ_p,η_q,ζ_r)、m_y(ξ_p,η_q,ζ_r)和m_z(ξ_p,η_q,ζ_r)分别表示中心点坐标位置为(ξ_p,η_q,ζ_r)的小棱柱体单元的磁化强度的x分量、y分量和z分量；h₁(x₀-ξ_p,y₀-η_q,z_n-ζ_r)、h₂(x₀-ξ_p,y₀-η_q,z_n-ζ_r)、h₃(x₀-ξ_p,y₀-η_q,z_n-ζ_r)、h₄(x₀-ξ_p,y₀-η_q,z_n-ζ_r)、h₅(x₀-ξ_p,y₀-η_q,z_n-ζ_r)、h₆(x₀-ξ_p,y₀-η_q,z_n-ζ_r)、h₇(x₀-ξ_p,y₀-η_q,z_n-ζ_r)、h₈(x₀-ξ_p,y₀-η_q,z_n-ζ_r)、h₉(x₀-ξ_p,y₀-η_q,z_n-ζ_r)为9个加权系数。