一种二度体重力异常快速数值模拟方法、装置和设备

1.本技术涉及计算机技术领域，特别是涉及一种二度体重力异常快速数值模拟方法、装置和计算机设备。


背景技术：

2.重力勘探是地球物理勘探中发展相对比较成熟的方法，由于其施工简单、成本低、效率高以及勘探深度大等优点已经被广泛地应用在区域调查、矿产资源勘探和石油天然气资源勘查中。正演是反演的基础，因此学者们对于重力场的正演计算非常重视。现实中所有地质体都是三维的，正演以三维模型为基础相对更加符合实际情况，但三维模型占用内存大，计算效率低；在实际模拟中，通常将其走向方向的尺度要远比垂直其走向方向的尺度大的多的地质体，采用二度体来代替，能够大幅缩减计算成本，相应的反演算法也更加容易实现。
3.对于大规模任意密度分布的地质体的二维重力异常计算，通常采用数值方法。文献(wu,l.,tian,g.high-precision fourier forward modeling of potential fields.geophysics,2014,79(5):g59-g68.)提出了一种频率域二度体重力异常正演的gauss-fft方法，该方法能够很好抑制传统快速傅里叶变换的边界效应问题，但随着高斯点个数的增加，其计算量成倍增加，不能实现精度和效率的统一。文献(chen,x.，chen,l.,luo,t.,et al.fast forward modelling of gravity anomalies of two-dimensional bodies of arbitrary shape and density distribution,applied geophysics,2020,17(5-6):688-695.)提出了一种基于任意截面形状和密度分布的二度体重力异常快速正演算法，该方法引入toeplitz矩阵与向量相乘的快速算法实现了组合矩形二度体重力异常的计算，但随着网格数的增加，该方法所需内存较大，同时降低了计算效率。
4.因此，现有算法存在计算效率低、占用内存大的问题。


技术实现要素：

5.基于此，有必要针对上述技术问题，提供一种能够提高计算效率、减小内存占用的二度体重力异常快速数值模拟方法、装置和计算机设备。
6.一种二度体重力异常快速数值模拟方法，所述方法包括：
7.根据待探测的目标体的分布，构建目标区域的二度体模型；所述二度体模型将所述目标区域沿x轴方向均匀剖分，沿z轴方向非均匀剖分，形成多个长方形网格，所述长方形网格的中心点为场点；
8.将所述目标体进行网格剖分，以网格的中心点为源点，并设定源点所在长方形的密度值，得到目标体密度分布模型；所述目标体密度分布模型的网格间隔与所述二度体模型对应部分的间隔一致；
9.根据所述二度体模型和所述目标体密度分布模型，得到任一场点上的重力异常公式；
10.对所述重力异常公式进行离散，得到任一观测高度上的重力异常积分公式；所述重力异常积分公式只在所述目标体密度分布模型区域进行积分；所述重力异常积分公式中包括源点的重力异常核函数单元积分系数以及源点密度值；
11.根据所述重力异常积分公式，得到任一观测高度上任一层所述目标体密度分布模型产生的重力异常一维卷积结果；
12.根据所述重力异常一维卷积结果，在所述目标体密度分布模型的垂向进行累加，得到整个目标体密度分布模型产生的重力异常。
13.在其中一个实施例中，还包括：根据所述二度体模型和所述目标体密度分布模型，得到任一场点上的重力异常公式为：
[0014][0015]
其中，gz表示重力异常，g表示万有引力常数，(x,z)表示场点坐标，(x
′
,z
′
)表示源点坐标，ρ(x
′
,z
′
)表示源点密度值。
[0016]
在其中一个实施例中，还包括：对所述重力异常公式进行离散，得到任一观测高度上的重力异常积分公式为：
[0017][0018]
其中，(xi,z0)表示二度体模型中编号为(i,0)的小长方形的中心坐标，(x
′m,z
′n)表示目标体密度分布模型中编号为(m,n)的小长方形的中心坐标；ρ(x
′m,z
′n)表示目标体密度分布模型中编号为(m,n)小长方形的密度值；i＝0,1,
…nx-1；m＝0,1,
…nxs-1；n＝0,1,
…nzs-1，n
x
表示所述二度体模型在x方向上剖分的网格数，n
xs
表示所述目标体密度分布模型在x方向上剖分的网格数，n
zs
表示所述目标体密度分布模型在z方向上剖分的网格数；f(x
i-x
′m,z
0-z
′n)表示重力异常核函数单元积分系数。
[0019]
在其中一个实施例中，还包括：根据所述重力异常积分公式，得到任一观测高度上任一层所述目标体密度分布模型产生的重力异常表示式：
[0020][0021]
其中，重力异常核函数单元积分系数表示为：
[0022][0023]
δx，δz分别表示x方向和z方向的网格间隔。
[0024]
在其中一个实施例中，还包括：通过一维快速离散卷积算法实现所述重力异常核函数单元积分系数和所述源点密度值的快速相乘。
[0025]
在其中一个实施例中，还包括：将所述重力异常核函数单元积分系数展开为向量形式，得到重力异常核函数单元积分系数向量
[0026][0027]
将所述源点密度值展开为向量形式，得到源点密度扩展向量
[0028]
通过一维快速离散卷积算法实现所述重力异常核函数单元积分系数和所述源点密度值的快速相乘：
[0029][0030]
其中，表示任一观测高度上第n层所述目标体密度分布模型产生的重力异常，f和f-1
分别表示一维离散傅里叶正反变换算子；表示提取矩阵的前n
x
个元素。
[0031]
在其中一个实施例中，还包括：根据所述重力异常一维卷积结果，在所述目标体密度分布模型的垂向进行累加，得到整个目标体密度分布模型产生的重力异常为：
[0032][0033]
其中，gz(xi,z0)表示观测高度z0上得到的整个目标体密度分布模型产生的重力异常。
[0034]
一种二度体重力异常快速数值模拟装置，所述装置包括：
[0035]
二度体模型构建模块，用于根据待探测的目标体的分布，构建目标区域的二度体模型；所述二度体模型将所述目标区域沿x轴方向均匀剖分，沿z轴方向非均匀剖分，形成多个长方形网格，所述长方形网格的中心点为场点；
[0036]
目标体密度分布模型构建模型，将所述目标体进行网格剖分，以网格的中心点为源点，并设定源点所在长方形的密度值，得到目标体密度分布模型；所述目标体密度分布模型的网格间隔与所述二度体模型对应部分的间隔一致；
[0037]
重力异常一维卷积结果确定模块，用于根据所述二度体模型和所述目标体密度分布模型，得到任一场点上的重力异常公式；对所述重力异常公式进行离散，得到任一观测高度上的重力异常积分公式；所述重力异常积分公式只在所述目标体密度分布模型区域进行积分；所述重力异常积分公式中包括源点的重力异常核函数单元积分系数以及源点密度值；根据所述重力异常积分公式，得到任一观测高度上任一层所述目标体密度分布模型产生的重力异常一维卷积结果；
[0038]
垂向叠加模块，用于根据所述重力异常一维卷积结果，在所述目标体密度分布模型的垂向进行累加，得到整个目标体密度分布模型产生的重力异常。
[0039]
一种计算机设备，包括存储器和处理器，所述存储器存储有计算机程序，所述处理器执行所述计算机程序时实现以下步骤：
[0040]
根据待探测的目标体的分布，构建目标区域的二度体模型；所述二度体模型将所述目标区域沿x轴方向均匀剖分，沿z轴方向非均匀剖分，形成多个长方形网格，所述长方形网格的中心点为场点；
[0041]
将所述目标体进行网格剖分，以网格的中心点为源点，并设定源点所在长方形的密度值，得到目标体密度分布模型；所述目标体密度分布模型的网格间隔与所述二度体模型对应部分的间隔一致；
[0042]
根据所述二度体模型和所述目标体密度分布模型，得到任一场点上的重力异常公式；
[0043]
对所述重力异常公式进行离散，得到任一观测高度上的重力异常积分公式；所述重力异常积分公式只在所述目标体密度分布模型区域进行积分；所述重力异常积分公式中包括源点的重力异常核函数单元积分系数以及源点密度值；
[0044]
根据所述重力异常积分公式，得到任一观测高度上任一层所述目标体密度分布模型产生的重力异常一维卷积结果；
[0045]
根据所述重力异常一维卷积结果，在所述目标体密度分布模型的垂向进行累加，得到整个目标体密度分布模型产生的重力异常。
[0046]
一种计算机可读存储介质，其上存储有计算机程序，所述计算机程序被处理器执行时实现以下步骤：
[0047]
根据待探测的目标体的分布，构建目标区域的二度体模型；所述二度体模型将所述目标区域沿x轴方向均匀剖分，沿z轴方向非均匀剖分，形成多个长方形网格，所述长方形网格的中心点为场点；
[0048]
将所述目标体进行网格剖分，以网格的中心点为源点，并设定源点所在长方形的密度值，得到目标体密度分布模型；所述目标体密度分布模型的网格间隔与所述二度体模
型对应部分的间隔一致；
[0049]
根据所述二度体模型和所述目标体密度分布模型，得到任一场点上的重力异常公式；
[0050]
对所述重力异常公式进行离散，得到任一观测高度上的重力异常积分公式；所述重力异常积分公式只在所述目标体密度分布模型区域进行积分；所述重力异常积分公式中包括源点的重力异常核函数单元积分系数以及源点密度值；
[0051]
根据所述重力异常积分公式，得到任一观测高度上任一层所述目标体密度分布模型产生的重力异常一维卷积结果；
[0052]
根据所述重力异常一维卷积结果，在所述目标体密度分布模型的垂向进行累加，得到整个目标体密度分布模型产生的重力异常。
[0053]
上述二度体重力异常快速数值模拟方法、装置、计算机设备和存储介质，通过根据待探测的目标体的分布，构建目标区域沿x轴方向均匀剖分，沿z轴方向非均匀剖分的二度体模型，以剖分形成的长方形网格的中心点为场点；将目标体进行网格剖分，以网格的中心点为源点，并设定源点所在长方形的密度值，得到目标体密度分布模型，其中，目标体密度分布模型的网格间隔与二度体模型一致；根据二度体模型和目标体密度分布模型，得到任一场点上的重力异常公式；对重力异常公式进行离散，得到任一观测高度上的重力异常积分公式，其中，重力异常积分公式只在目标体密度分布模型区域进行积分；根据重力异常积分公式，得到任一观测高度上任一层目标体密度分布模型产生的重力异常一维卷积结果；根据重力异常一维卷积结果，在目标体密度分布模型的垂向进行累加，得到整个目标体密度分布模型产生的重力异常。本发明对目标区域和待探测的目标体分别构建了网格模型，在计算重力异常时只需对目标体网格进行积分，极大地降低了内存需求，减少了积分系数计算，在保证原有精度不变的前提下提高了计算效率。
附图说明
[0054]
图1为一个实施例中二度体重力异常快速数值模拟方法的流程示意图；
[0055]
图2为另一个实施例中二度体重力异常快速数值模拟方法的流程示意图；
[0056]
图3为一个具体实施例中截面为长方形二度体模型示意图；
[0057]
图4为一个具体实施例中本发明方法和解析解计算的地面重力异常；
[0058]
图5为一个具体实施例中本发明方法与解析解计算的重力异常相对误差；
[0059]
图6为一个具体实施例中本发明方法与传统算法随着网格剖分大小的时间变化曲线；
[0060]
图7为一个实施例中二度体重力异常快速数值模拟装置的结构框图；
[0061]
图8为一个实施例中计算机设备的内部结构图。
具体实施方式
[0062]
为了使本技术的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本技术进行进一步详细说明。应当理解，此处描述的具体实施例仅仅用以解释本技术，并不用于限定本技术。
[0063]
在一个实施例中，如图1所示，提供了一种二度体重力异常快速数值模拟方法，包
括以下步骤：
[0064]
步骤102，根据待探测的目标体的分布，构建目标区域的二度体模型。
[0065]
待探测的目标体为引起二度异常的地质体，如层状矿体、脉状矿体、向斜、背斜等。
[0066]
二度体模型场点坐标包括x轴方向和z轴方向，二度体模型将目标区域沿x轴方向均匀剖分，沿z轴方向非均匀剖分，形成多个长方形网格，长方形网格的中心点为场点。沿z轴方向非均匀剖分具体指的是，z轴方向根据场变化的快慢进行剖分，在场变化快的地方网格间隔减小，在场变化慢的地方网格适当稀疏，在保证精度的前提下提高了计算效率。
[0067]
步骤104，将目标体进行网格剖分，以网格的中心点为源点，并设定源点所在长方形的密度值，得到目标体密度分布模型。
[0068]
目标体密度分布模型的网格间隔与二度体模型对应部分的间隔一致。
[0069]
步骤106，根据二度体模型和目标体密度分布模型，得到任一场点上的重力异常公式。
[0070]
由于源点是基于另外建立的目标体密度分布模型，目标体密度分布模型的网格数远远少于整个目标区域的网格。因此，在进行重力异常的计算时只需对目标体网格进行积分，极大地降低了内存需求，减少了积分系数计算，在保证原有精度不变的前提下提高了计算效率。
[0071]
步骤108，对重力异常公式进行离散，得到任一观测高度上的重力异常积分公式。
[0072]
重力异常积分公式只在目标体密度分布模型区域进行积分。重力异常积分公式中包括源点的重力异常核函数单元积分系数以及源点密度值。
[0073]
步骤110，根据重力异常积分公式，得到任一观测高度上任一层目标体密度分布模型产生的重力异常一维卷积结果。
[0074]
步骤112，根据重力异常一维卷积结果，在目标体密度分布模型的垂向进行累加，得到整个目标体密度分布模型产生的重力异常。
[0075]
重力异常一维卷积结果为某一层目标体模型产生的重力异常，在累加时垂向多层重力异常时，只需累加目标体垂向剖分的网格层数，即可得到整个组合模型产生的重力异常，进一步减小计算量，提高计算效率。
[0076]
上述二度体重力异常快速数值模拟方法中，通过根据待探测的目标体的分布，构建目标区域沿x轴方向均匀剖分，沿z轴方向非均匀剖分的二度体模型，以剖分形成的长方形网格的中心点为场点；将目标体进行网格剖分，以网格的中心点为源点，并设定源点所在长方形的密度值，得到目标体密度分布模型，其中，目标体密度分布模型的网格间隔与二度体模型一致；根据二度体模型和目标体密度分布模型，得到任一场点上的重力异常公式；对重力异常公式进行离散，得到任一观测高度上的重力异常积分公式，其中，重力异常积分公式只在目标体密度分布模型区域进行积分；根据重力异常积分公式，得到任一观测高度上任一层目标体密度分布模型产生的重力异常一维卷积结果；根据重力异常一维卷积结果，在目标体密度分布模型的垂向进行累加，得到整个目标体密度分布模型产生的重力异常。本发明对目标区域和待探测的目标体分别构建了网格模型，在计算重力异常时只需对目标体网格进行积分，极大地降低了内存需求，减少了积分系数计算，在保证原有精度不变的前提下提高了计算效率。
[0077]
在其中一个实施例中，还包括：根据二度体模型和目标体密度分布模型，得到任一
场点上的重力异常公式为：
[0078][0079]
其中，gz表示重力异常，g表示万有引力常数，(x,z)表示场点坐标，(x
′
,z
′
)表示源点坐标，ρ(x
′
,z
′
)表示源点密度值。
[0080]
在其中一个实施例中，还包括：对重力异常公式进行离散，得到任一观测高度上的重力异常积分公式为：
[0081][0082]
其中，(xi,z0)表示二度体模型中编号为(i,0)的小长方形的中心坐标，(x
′m,z
′n)表示目标体密度分布模型中编号为(m,n)的小长方形的中心坐标；ρ(x
′m,z
′n)表示目标体密度分布模型中编号为(m,n)小长方形的密度值；i＝0,1,
…nx-1；m＝0,1,
…nxs-1；n＝0,1,
…nzs-1，n
x
表示二度体模型在x方向上剖分的网格数，n
xs
表示目标体密度分布模型在x方向上剖分的网格数，n
zs
表示目标体密度分布模型在z方向上剖分的网格数；f(x
i-x
′m,z
0-z
′n)表示重力异常核函数单元积分系数。
[0083]
可以看出本方法只需对异常体网格进行积分，极大地降低了内存需求，减少了积分系数计算，在保证原有精度不变的前提下提高了计算效率。
[0084]
在其中一个实施例中，还包括：根据重力异常积分公式，得到任一观测高度上任一层目标体密度分布模型产生的重力异常表示式：
[0085][0086]
其中，重力异常核函数单元积分系数表示为：
[0087][0088]
δx，δz分别表示x方向和z方向的网格间隔。
[0089]
在其中一个实施例中，还包括：通过一维快速离散卷积算法实现重力异常核函数单元积分系数和源点密度值的快速相乘。
[0090]
在其中一个实施例中，还包括：将重力异常核函数单元积分系数展开为向量形式，得到重力异常核函数单元积分系数向量
[0091][0092]
将源点密度值展开为向量形式，得到源点密度扩展向量
[0093][0094]
通过一维快速离散卷积算法实现重力异常核函数单元积分系数和源点密度值的快速相乘：
[0095][0096]
其中，表示任一观测高度上第n层目标体密度分布模型产生的重力异常，f和f-1
分别表示一维离散傅里叶正反变换算子；表示提取矩阵的前n
x
个元素。
[0097]
在其中一个实施例中，还包括：根据重力异常一维卷积结果，在目标体密度分布模型的垂向进行累加，得到整个目标体密度分布模型产生的重力异常为：
[0098]
[0099]
其中，gz(xi,z0)表示观测高度z0上得到的整个目标体密度分布模型产生的重力异常。
[0100]
可以看出只需累加异常体垂向剖分的网格层数，即n
zs
层，则可得到整个组合模型产生的重力异常，进一步减小计算量，提高计算效率。
[0101]
在另一个实施例中，如图2所示，提供了一种二度体重力异常快速数值模拟方法，包括：
[0102]
模型构建：对研究区域进行网格剖分，并确定水平和垂向小长方形个数；对异常体进行剖分，并确定异常体在研究区域的相对位置，水平和垂向剖分长方形个数；
[0103]
给出二度体重力异常计算公式：根据模型参数对二度体重力异常计算公式进行离散；
[0104]
计算核函数单元积分系数：对离散公式进行单元分析，得到核函数单元积分系数；
[0105]
一维快速离散卷积：采用一维快速离散卷积实现核函数与密度的快速相乘；
[0106]
输出重力异常：对各层模型的重力异常累加得到整个组合模型产生的重力异常。
[0107]
应该理解的是，虽然图1、2的流程图中的各个步骤按照箭头的指示依次显示，但是这些步骤并不是必然按照箭头指示的顺序依次执行。除非本文中有明确的说明，这些步骤的执行并没有严格的顺序限制，这些步骤可以以其它的顺序执行。而且，图1、2中的至少一部分步骤可以包括多个子步骤或者多个阶段，这些子步骤或者阶段并不必然是在同一时刻执行完成，而是可以在不同的时刻执行，这些子步骤或者阶段的执行顺序也不必然是依次进行，而是可以与其它步骤或者其它步骤的子步骤或者阶段的至少一部分轮流或者交替地执行。
[0108]
在一个具体实施例中，目标区域有一个截面为规则矩形的二度体模型如图3所示，模拟区域范围为：x方向从-1000m到1000m，z方向从0m到1000m(z轴垂向向下为正)，水平网格间隔均为10m，垂向网格间隔为5m，水平和垂向网格个数为200
×
200，截面为矩形的异常体分布范围为：x方向从-100m到100m，z方向从400m到600m，剩余密度为1000kg/m3，计算地面测线的重力异常。
[0109]
本发明方法利用fortran语言编程实现，运行程序所用的个人电脑配置为：cpu-inter core i7-8700，主频为3.2ghz，运行内存为24.00gb。图4为本发明方法和解析解计算的地面的重力异常，可以看出本发明方法和解析解计算的重力异常吻合很好，图5为解析解和本发明方法计算的相对误差，可以看出相对误差都小于1.8
×
10-5，可以看出本方法计算精度较高。
[0110]
为了进一步对本发明的计算效率进行验证，采用本发明方法和传统的一维离散卷积方法进行对比。图6即给出的本发明方法和传统一维离散卷积方法随着网格剖分大小的时间变化曲线，当网格剖分大小为10000
×
10000时，传统方法需要8s而本发明方法仅需2s；当网格剖分大小为40000
×
40000时，传统方法需要194s，本发明方法需12.56s，本发明方法计算效率提高了15.4倍；可以看出随着剖分网格个数的增加，本算法的优势随之更大，计算效率更高。
[0111]
在一个实施例中，如图7所示，提供了一种二度体重力异常快速数值模拟装置，包括：二度体模型构建模块702、目标体密度分布模型构建模型704、重力异常一维卷积结果确定模块706和垂向叠加模块708，其中：
[0112]
二度体模型构建模块702，用于根据待探测的目标体的分布，构建目标区域的二度体模型；二度体模型将目标区域沿x轴方向均匀剖分，沿z轴方向非均匀剖分，形成多个长方形网格，长方形网格的中心点为场点；
[0113]
目标体密度分布模型构建模型704，将目标体进行网格剖分，以网格的中心点为源点，并设定源点所在长方形的密度值，得到目标体密度分布模型；目标体密度分布模型的网格间隔与二度体模型对应部分的间隔一致；
[0114]
重力异常一维卷积结果确定模块706，用于根据二度体模型和目标体密度分布模型，得到任一场点上的重力异常公式；对重力异常公式进行离散，得到任一观测高度上的重力异常积分公式；重力异常积分公式只在目标体密度分布模型区域进行积分；重力异常积分公式中包括源点的重力异常核函数单元积分系数以及源点密度值；根据重力异常积分公式，得到任一观测高度上任一层目标体密度分布模型产生的重力异常一维卷积结果；
[0115]
垂向叠加模块708，用于根据重力异常一维卷积结果，在目标体密度分布模型的垂向进行累加，得到整个目标体密度分布模型产生的重力异常。
[0116]
重力异常一维卷积结果确定模块706还用于根据二度体模型和目标体密度分布模型，得到任一场点上的重力异常公式为：
[0117][0118]
其中，gz表示重力异常，g表示万有引力常数，(x,z)表示场点坐标，
[0119]
(x
′
,z
′
)表示源点坐标，ρ(x
′
,z
′
)表示源点密度值。
[0120]
重力异常一维卷积结果确定模块706还用于对重力异常公式进行离散，得到任一观测高度上的重力异常积分公式为：
[0121][0122]
其中，(xi,z0)表示二度体模型中编号为(i,0)的小长方形的中心坐标，(x
′m,z
′n)表示目标体密度分布模型中编号为(m,n)的小长方形的中心坐标；ρ(x
′m,z
′n)表示目标体密度分布模型中编号为(m,n)小长方形的密度值；i＝0,1,
…nx-1；m＝0,1,
…nxs-1；n＝0,1,
…nzs-1，n
x
表示二度体模型在x方向上剖分的网格数，n
xs
表示目标体密度分布模型在x方向上剖分的网格数，n
zs
表示目标体密度分布模型在z方向上剖分的网格数；f(x
i-x
′m,z
0-z
′n)表示重力异常核函数单元积分系数。
[0123]
重力异常一维卷积结果确定模块706还用于根据重力异常积分公式，得到任一观测高度上任一层目标体密度分布模型产生的重力异常表示式：
[0124][0125]
其中，重力异常核函数单元积分系数表示为：
[0126][0127]
δx，δz分别表示x方向和z方向的网格间隔。
[0128]
重力异常一维卷积结果确定模块706还用于通过一维快速离散卷积算法实现重力异常核函数单元积分系数和源点密度值的快速相乘。
[0129]
重力异常一维卷积结果确定模块706还用于将重力异常核函数单元积分系数展开为向量形式，得到重力异常核函数单元积分系数向量
[0130]
将源点密度值展开为向量形式，得到源点密度扩展向量
[0131]
通过一维快速离散卷积算法实现重力异常核函数单元积分系数和源点密度值的快速相乘：
[0132][0133]
其中，表示任一观测高度上第n层目标体密度分布模型产生的重力异常，f和f-1
分别表示一维离散傅里叶正反变换算子；表示提取矩阵的前n
x
个元素。
[0134]
垂向叠加模块708根据重力异常一维卷积结果，在目标体密度分布模型的垂向进行累加，得到整个目标体密度分布模型产生的重力异常为：
[0135][0136]
其中，gz(xi,z0)表示观测高度z0上得到的整个目标体密度分布模型产生的重力异常。
[0137]
关于二度体重力异常快速数值模拟装置的具体限定可以参见上文中对于二度体
重力异常快速数值模拟方法的限定，在此不再赘述。上述二度体重力异常快速数值模拟装置中的各个模块可全部或部分通过软件、硬件及其组合来实现。上述各模块可以硬件形式内嵌于或独立于计算机设备中的处理器中，也可以以软件形式存储于计算机设备中的存储器中，以便于处理器调用执行以上各个模块对应的操作。
[0138]
在一个实施例中，提供了一种计算机设备，该计算机设备可以是终端，其内部结构图可以如图8所示。该计算机设备包括通过系统总线连接的处理器、存储器、网络接口、显示屏和输入装置。其中，该计算机设备的处理器用于提供计算和控制能力。该计算机设备的存储器包括非易失性存储介质、内存储器。该非易失性存储介质存储有操作系统和计算机程序。该内存储器为非易失性存储介质中的操作系统和计算机程序的运行提供环境。该计算机设备的网络接口用于与外部的终端通过网络连接通信。该计算机程序被处理器执行时以实现一种二度体重力异常快速数值模拟方法。该计算机设备的显示屏可以是液晶显示屏或者电子墨水显示屏，该计算机设备的输入装置可以是显示屏上覆盖的触摸层，也可以是计算机设备外壳上设置的按键、轨迹球或触控板，还可以是外接的键盘、触控板或鼠标等。
[0139]
本领域技术人员可以理解，图8中示出的结构，仅仅是与本技术方案相关的部分结构的框图，并不构成对本技术方案所应用于其上的计算机设备的限定，具体的计算机设备可以包括比图中所示更多或更少的部件，或者组合某些部件，或者具有不同的部件布置。
[0140]
在一个实施例中，提供了一种计算机设备，包括存储器和处理器，该存储器存储有计算机程序，该处理器执行计算机程序时实现上述方法实施例中的步骤。
[0141]
在一个实施例中，提供了一种计算机可读存储介质，其上存储有计算机程序，计算机程序被处理器执行时实现上述方法实施例中的步骤。
[0142]
本领域普通技术人员可以理解实现上述实施例方法中的全部或部分流程，是可以通过计算机程序来指令相关的硬件来完成，所述的计算机程序可存储于一非易失性计算机可读取存储介质中，该计算机程序在执行时，可包括如上述各方法的实施例的流程。其中，本技术所提供的各实施例中所使用的对存储器、存储、数据库或其它介质的任何引用，均可包括非易失性和/或易失性存储器。非易失性存储器可包括只读存储器(rom)、可编程rom(prom)、电可编程rom(eprom)、电可擦除可编程rom(eeprom)或闪存。易失性存储器可包括随机存取存储器(ram)或者外部高速缓冲存储器。作为说明而非局限，ram以多种形式可得，诸如静态ram(sram)、动态ram(dram)、同步dram(sdram)、双数据率sdram(ddrsdram)、增强型sdram(esdram)、同步链路(synchlink)dram(sldram)、存储器总线(rambus)直接ram(rdram)、直接存储器总线动态ram(drdram)、以及存储器总线动态ram(rdram)等。
[0143]
以上实施例的各技术特征可以进行任意的组合，为使描述简洁，未对上述实施例中的各个技术特征所有可能的组合都进行描述，然而，只要这些技术特征的组合不存在矛盾，都应当认为是本说明书记载的范围。
[0144]
以上所述实施例仅表达了本技术的几种实施方式，其描述较为具体和详细，但并不能因此而理解为对发明专利范围的限制。应当指出的是，对于本领域的普通技术人员来说，在不脱离本技术构思的前提下，还可以做出若干变形和改进，这些都属于本技术的保护范围。因此，本技术专利的保护范围应以所附权利要求为准。