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;
B5:输出地下网格密度mq。
2.根据权利要求1所述的基于大数据计算的航空重力数据物性层析处理方法,其特征在于:所述步骤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均方误差恒定或最小。