[发明专利]一种重力及梯度异常正演方法

权利要求书

1.一种重力及梯度异常正演方法，其特征在于，包括以下步骤：

获取格架函数并储存该格架函数；

建立模型空间，对部分模型空间中的异常体剩余密度赋值；

确定观测点并确定该观测点周围选取的计算格架函数的范围，该范围形成一长方体，确定观测点与该长方体的相对位置关系；

根据所述相对位置关系并利用平移等效性同时调用格架函数，代数求和获得该长方体单元对观测点的剩余密度异常体；

计算整个模型空间对该观测点的重力异常体；

循环全部完成，得出工作区的磁张量异常体；

在建立模型空间，对部分模型空间中的异常体剩余密度赋值时，对模型空间中一部分进行选取以形成子空间，确定后计算子空间；获取格架函数为计算子空间重力异常体的格架函数并储存该格架函数；

其中，所述计算子空间包含将计算子空间剖分成网格，所述计算子空间只限于子空间内的网格点在观测点产生的异常体值；所述计算子空间中包含确定网格的间距、异常体在网格中的位置，并将计算子空间计算区平面的单位换算为实际里程，并对异常体的剩余密度赋值；

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

确定观测点与长方体的相对位置关系包含在计算子空间内，观测点每循环到一处，所述长方体在计算字空间内遍历循环一次，长方体在循环中每到一处，观测点与其有一个相对位置关系。