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

权利要求书

1.复杂地质体井中重力三分量计算方法，其特征在于，包括：

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

计算所有小棱柱体单元在井中观测点处产生的井中重力。

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

3.根据权利要求1或2所述的复杂地质体井中重力三分量计算方法，其特征在于，采用三维离散卷积法计算所有小棱柱体单元在井中观测点(x₀,y₀,z_n)处产生的井中重力，包括3个分量，分别是g_x(x₀,y₀,z_n)、g_y(x₀,y₀,z_n)、g_z(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方向小棱柱体单元的剖分个数；ρ(ξ_p,η_q,ζ_r)表示中心点坐标位置为(ξ_p,η_q,ζ_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)分别表示为重力分量g_x(x₀,y₀,z_n)、g_y(x₀,y₀,z_n)、g_z(x₀,y₀,z_n)的加权系数。