[发明专利]一种获取地下三维密度结构的方法

权利要求书

1.一种获取地下三维密度结构的方法，其特征在于，包括如下步骤：

确定建模区域的纵坐标、横坐标、深度及二维2D剖面间距；

根据先验信息及建模区域，利用朴素贝叶斯分类方法进行重力反演建模，得到初始地质模型，具体为：

根据所述先验信息确定出模型参数矢量；其中，所述模型参数矢量包括：地质体的埋深、厚度、密度及波速度；

根据所述建模区域、先验信息确定出核函数；其中，所述核函数表示观测数据与模型参数之间的函数关系；

根据所述模型参数矢量及核函数确定出观测向量数据，所述观测向量数据构成训练样本；

将所述观测向量数据划分为不同的类别，并计算出每个类别的观测向量数据在所述训练样本中的频率；

根据所述核函数、每个类别的观测向量数据在所述训练样本中的频率，计算出模型区域中地质体的最大似然解；

使用分类器鉴别满足所述地质体最大似然解的模型参数；

根据所述分类器的鉴别结果生成地质模型，并计算出所述地质模型的重力异常值和实际重力异常观测值的拟合方差；

其中，当所述拟合方差＜设定的正数ε时，由分类器的鉴别结果反演出的地质模型参数作为反演结果，得到初始地质模型；

当所述拟合方差≥设定的正数ε时，修正模型参数矢量，重复上述步骤；

采用离散体模拟方法将所述初始地质模型构建成2D地质模型；

将所述2D地质模型拼接拟合成3D地质模型；

对所述3D地质模型进行可视化和结构解译处理，得到地下三维密度结构。

2.根据权利要求1所述的方法，其特征在于，所述分类器为隐马尔科夫链分类器。

3.根据权利要求1所述的方法，其特征在于，所述设定的正数ε为0.001-0.03。

4.根据权利要求1所述的方法，其特征在于，所述先验信息包括地球物理信息和地质信息；

所述地球物理信息包括重力场信息、电磁信息、地震信息及岩石物理信息；

所述地质信息包括地层年代信息、钻孔资料信息及区域地质背景信息。

5.根据权利要求1所述的方法，其特征在于，所述采用离散体模拟方法将所述初始地质模型构建成2D地质模型的步骤，具体为：

采用离散体模拟方法，结合先验信息对所述初始地质模型进行重力反演，得到2D地质模型。

6.根据权利要求1所述的方法，其特征在于，将所述2D地质模型拟合成3D地质模型的步骤，具体为：

采用人机交互法对所述2D地质模型的每个2D剖面进行修改，叠加、拼接拟合出3D地质模型。

7.根据权利要求1所述的方法，其特征在于，

采用GOCAD三维可视化模拟软件对3D地质模型进行可视化和结构解译处理，得到地下三维密度结构；或

采用Encom PA地质资料分析、模拟软件对3D地质模型进行可视化和结构解译处理，得到地下三维密度结构。