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

摘要

一种基于大数据计算的航空重力数据物性层析处理方法，包括以下步骤：对地下网格剖分形成的小立方体总数进行切分，构建RDD，并持久化到内存；计算概率函数；计算地下网格的密度；迭代计算地表网格重力正演数据和地下网格密度；输出地下网格密度。本发明将重力层析反演问题由传统的欠定方程求解方法，转变成为正演迭代计算方法，计算步骤简单、计算过程稳定、计算结果准确度高；采用大数据分布式并行计算框架Spark，使上述迭代计算和复杂多路处理过程可以充分借助内存数据存储实现，达到了接近实时的处理效果。

主权项

1.一种基于大数据计算的航空重力数据物性层析处理方法，其特征在于：包括以下步骤：B1：对地下网格剖分形成的小立方体总数进行切分，构建RDD，并持久化到内存；B2：通过按照地表网格采集到的航空重力数据d1，计算概率函数ηq：其中，N为航空重力数据d1中数据采集点的数量，Δg(x,y,z)为航空重力数据d1中位于地表(x,y,z)采集点的重力异常测量值，Δgu(x,y,z)为地表(x,y,z)采集点重力异常的格林函数；B3：计算地下网格的密度mq：mq＝Δρ·ηq(2)B4：迭代计算地表网格重力正演数据d2和地下网格密度mq；所述步骤B4中计算地表网格重力正演数据d2和地下网格密度mq的迭代步骤如下：B4.1：对RDD进行MAP算子，按照地表网格计算各自的重力正演数据d2；其中，Δg(x,y,z)为重力正演数据d2中位于地表(x,y,z)采集点的重力异常正演值，G为万有引力常数，Q为地下网格剖分形成的小立方体的数量，(xq,yq,zq)为地下网格剖分形成的某个小立方体的中心点坐标，Δρq是所述小立方体的密度，Δvq是所述小立方体的体积；B4.2：对RDD进行Reduce算子，合并重力正演数据d2；B4.3对重力正演数据与航空重力数据之间的差值d2‑d1，应用公式(1)计算得到概率函数偏差ηq′；B4.4更新地下网格密度mq：mq＝mq+Δρ·ηq′(4)B4.5：重复步骤B4.1至B4.4，直到重力正演数据d2与航空重力数据d1均方误差恒定或最小；B5：输出地下网格密度mq。