[发明专利]基于网格点格架函数的重力异常及梯度异常的正演方法

权利要求书

1.基于网格点格架函数的重力异常及梯度异常的正演方法，其特征在于，包括以下步骤：

步骤1：开始；

步骤2：剖分模型空间，对异常体剩余密度赋值；

步骤3：分别计算重力异常和梯度异常的格架函数；

步骤4：判断观测点与长方体单元相对位置；

步骤5：利用对称互换性和平移等效性；

步骤6：调用重力异常和梯度异常的格架函数；

步骤7：代数求和获得该长方体对该观测点的重力异常以及梯度异常；

步骤8：得出整个模型体对该观测点的重力异常和梯度异常；

步骤9：计算平面内的循环全部完成；

步骤10：得出结果;

所述步骤2具体为：将计算空间剖分成网格，确定网格的间距、异常体在网格中的位置，并将计算空间计算区平面的单位换算为实际里程，并对异常体的剩余密度赋值;

步骤3具体为：确定计算区平面内观测点的坐标，确定计算空间内剖分所得长方体模型单元的坐标，观测点在计算区平面内循环，长方体单元在整个计算空间内循环，分别计算模型空间内所有网格点对第一个观测点的重力异常和梯度异常的格架函数，然后存储待后期计算调用;

步骤4具体为：在计算空间内，观测点每循环到一处，长方体模型单元在计算空间内遍历循环一次，长方体模型单元在循环中每到一处，观测点与其有一个相对位置关系;

步骤5-9具体为：将观测点和长方体单元的位置参数代入重力异常格架函数中，利用对称互换性与平移等效性，即可调出8个已存储的网格点对观测点格架函数，然后代数求和，获得该长方体单元对观测点的重力异常和梯度异常，将所有长方体单元对该观测点的重力异常和梯度异常求和，即为整个模型体对该观测点的重力异常和梯度异常，当计算空间内的循环全部完成，即可得出计算区的重力异常和所有梯度异常。