[发明专利]一种应用于大数据计算的航空重力数据物性层析计算方法

权利要求书

1.一种应用于大数据计算的航空重力数据物性层析计算方法，其特征在于：包括以下步骤：

A1：通过按照地表网格采集到的航空重力数据d₁，计算概率函数η_q：

其中，N为航空重力数据d₁中数据采集点的数量，Δg(x,y,z)为航空重力数据d₁中位于地表(x,y,z)采集点的重力异常测量值，Δg_u(x,y,z)为地表(x,y,z)采集点重力异常的格林函数；

A2：计算地下网格的密度m_q：

m_q＝Δρ·η_q (2)

其中，Δρ为迭代计算的步进量；

A3：按照地表网格计算重力正演数据d₂；

其中，Δg(x,y,z)为重力正演数据d₂中位于地表(x,y,z)采集点的重力异常正演值，G为万有引力常数，Q为地下网格剖分形成的小立方体的数量，(x_q,y_q,z_q)为地下网格剖分形成的某个小立方体的中心点坐标，Δρ_q是所述小立方体的密度，Δv_q是所述小立方体的体积；

A4：对重力正演数据与航空重力数据之间的差值d₂-d₁，应用公式(1)计算得到概率函数偏差η_q′；

A5：更新地下网格密度m_q：

m_q＝m_q+Δρ·η_q′ (4)

A6：重复步骤A3至A5，直到重力正演数据d₂与航空重力数据d₁均方误差恒定或最小，输出地下网格密度m_q。