[发明专利]一种用于提取地球物理三维断裂信息的方法和计算设备

权利要求书

1.一种用于提取地球物理断裂信息的方法，其特征在于，方法包括：

根据二维谱矩和三维密度场数据构建出三维二阶对称谱矩；

根据所述三维二阶对称谱矩获得对应的协方差矩阵；

对所述协方差矩阵进行特征值分解，获得对角矩阵，根据所述对角矩阵的元素构建统计不变量；

由所述统计不变量构建地球物理三维断裂信息。

2.如权利要求1所述的用于提取地球物理三维断裂信息的方法，其特征在于，所述根据二维谱矩和三维密度场数据构建出三维二阶对称谱矩的方法包括：

三维密度场数据由重力场数据获得三维密度场函数，将所述三维密度场函数进行快速傅里叶变换得到三维密度场的傅里叶变换函数F(f_x，f_y，f_z)；

根据所述三维密度场的傅里叶变换函数计算得到三维密度场的三维功率谱密度G(f_x，f_y，f_z)；

根据所述三维密度场的三维功率谱密度计算得到三维密度场对应的三维二阶对称谱矩m_pqr。

3.如权利要求2所述的用于提取地球物理三维断裂信息的方法，其特征在于，所述三维密度场数据由重力场数据获得三维密度场函数，将所述三维密度场函数进行快速傅里叶变换得到三维密度场的傅里叶变换函数的具体方法包括：

三维密度场数据由重力场数据获得三维密度场函数；

将所述三维密度场函数进行快速傅里叶变换得到三维密度场的傅里叶变换函数F(f_x，f_y，f_z)；

式中，f_x、f_y和f_z分别表示x、y和z三个相互垂直方向上的波数；在由x、y和z轴组成的笛卡尔直角坐标系中，ρ表示密度场函数；(x，y，z)是空间坐标位置；e表示自然对数的底数；i表示虚数单位；dx、dy和dz表示x、y和z的微元。

4.如权利要求2所述的用于提取地球物理三维断裂信息的方法，其特征在于，所述根据所述三维密度场的傅里叶变换函数计算得到三维密度场的三维功率谱密度的具体方法包括：

由所述三维密度场的傅里叶变换函数F(f_x，f_y，f_z)得到三维密度场的傅里叶变换的共轭函数F^*(f_x，f_y，f_z)；

三维密度场的三维功率谱密度G(f_x，f_y，f_z)由所述三维密度场的傅里叶变换函数F(f_x，f_y，f_z)与三维密度场的傅里叶变换的共轭函数F^*(f_x，f_y，f_z)的乘积得到，即：

G(f_x，f_y，f_z)＝F(f_x，f_y，f_z)F^*(f_x，f_y，f_z)；

式中，f_x、f_y和f_z分别表示x、y和z三个相互垂直方向上的波数。