1.一种二度体重力异常积分解数值模拟方法,其特征在于,所述方法包括:
根据观测点在笛卡尔坐标系下的位置信息,构建包含目标区域的二度体模型;所述观测点坐标包括x轴方向和z轴方向;所述二度体模型是将所述目标区域分成若干个矩形,所述二度体的截面形状用于确定每个矩形的空间域异常体剩余密度值;
根据所述二度体模型,确定对应的空间域重力异常表达式;所述空间域重力异常表达式包含空间域异常体剩余密度值;
对所述空间域重力异常表达式沿x轴方向进行一维傅里叶变换,得到波数域重力异常表达式;其中,所述波数域重力异常表达式包含波数域异常体剩余密度值;所述波数域异常体剩余密度值是对所述空间域异常体剩余密度值进行非均匀采样快速傅里叶变换得到的;
根据所述矩形x轴方向的尺寸信息,确定截止频率,根据所述截止频率确定波数域波数采样值范围,根据所述波数采样值范围和预先设置的采样点总数,在波数域进行非均匀采样,得到波数采样值;
将所述波数采样值代入所述波数域重力异常表达式,得到波数域重力异常场值;
通过对所述波数域重力异常场值进行一维非均匀快速傅里叶反变换,得到目标区域内任一点的空间域重力异常场值。
2.根据权利要求1所述的方法,其特征在于,所述空间域重力异常表达式为:
式中:
3.根据权利要求2所述的方法,其特征在于,所述波数域重力异常表达式为:
式中:
4.根据权利要求1所述的方法,其特征在于,所述非均匀采样快速傅里叶变换为:
式中:i为虚数单位,
其中,所述非均匀快速傅里叶变换的实现步骤为:
根据
式中:m表示过采样因子,
根据采样点的值和权重因子,计算新傅里叶变换基对应的傅里叶变换系数
采用均匀的快速傅里叶变换,得到:
式中,
5.根据权利要求1所述的方法,其特征在于,将所述目标区域分成若干个矩形包括:
将所述目标区域分成若干个矩形,每个矩形的x轴方向尺寸
6.根据权利要求1所述的方法,其特征在于,根据所述矩形沿x轴方向的尺寸信息,确定截止频率,包括:
根据所述矩形沿x轴方向的尺寸信息,确定截止频率为:
式中,
7.根据权利要求6所述的方法,其特征在于,根据所述截止频率确定波数域波数采样值范围,根据所述波数采样值范围和预先设置的采样点总数,在波数域进行非均匀采样,得到波数采样值,包括:
根据所述截止频率确定波数域波数采样值范围为
根据所述波数采样值范围得到波数采样值,包括:
在
在
式中,
8.一种二度体重力异常积分解数值模拟装置,其特征在于,所述装置包括:
二度体模型构建模块,根据观测点在笛卡尔坐标系下的位置信息,构建包含目标区域的二度体模型;所述观测点坐标包括x轴方向和z轴方向;所述二度体模型是将所述目标区域分成若干个矩形,所述二度体的截面形状用于确定每个矩形的空间域异常体剩余密度值;
空间域重力异常表达式确定模块,用于根据所述二度体模型,确定对应的空间域重力异常表达式;所述空间域重力异常表达式包含空间域异常体剩余密度值;
波数域重力异常表达式确定模块,用于对所述空间域重力异常表达式沿x轴方向进行一维傅里叶变换,得到波数域重力异常表达式;其中,所述波数域重力异常表达式包含波数域异常体剩余密度值;所述波数域异常体剩余密度值是对所述空间域异常体剩余密度值进行非均匀采样快速傅里叶变换得到的;
波数采样值确定模块,用于根据所述矩形沿x轴方向的尺寸信息,确定截止频率,根据所述截止频率确定波数域波数采样值范围,根据所述波数采样值范围和预先设置的采样点总数,在波数域进行非均匀采样,得到波数采样值;
波数域重力异常场值确定模块,用于将所述波数采样值代入所述波数域重力异常表达式,得到波数域重力异常场值;
空间域重力异常场值确定模块,用于通过对所述波数域重力异常场值进行一维非均匀快速傅里叶反变换,得到目标区域内任一点的空间域重力异常场值。
9.一种计算机设备,包括存储器和处理器,所述存储器存储有计算机程序,其特征在于,所述处理器执行所述计算机程序时实现权利要求1至7中任一项所述方法的步骤。
10.一种计算机可读存储介质,其上存储有计算机程序,其特征在于,所述计算机程序被处理器执行时实现权利要求1至7中任一项所述的方法的步骤。