CN113111927A - 一种基于地球物理数据使用堆叠法的三维地质建模方法 - Google Patents
一种基于地球物理数据使用堆叠法的三维地质建模方法 Download PDFInfo
- Publication number
- CN113111927A CN113111927A CN202110353243.3A CN202110353243A CN113111927A CN 113111927 A CN113111927 A CN 113111927A CN 202110353243 A CN202110353243 A CN 202110353243A CN 113111927 A CN113111927 A CN 113111927A
- Authority
- CN
- China
- Prior art keywords
- dimensional
- data
- model
- stacking
- lithology
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 60
- 239000010438 granite Substances 0.000 claims abstract description 17
- 230000010354 integration Effects 0.000 claims abstract description 12
- 238000005553 drilling Methods 0.000 claims abstract description 11
- 239000011435 rock Substances 0.000 claims abstract description 11
- 238000004458 analytical method Methods 0.000 claims abstract description 10
- 238000000605 extraction Methods 0.000 claims abstract description 10
- 238000007781 pre-processing Methods 0.000 claims abstract description 9
- 238000003672 processing method Methods 0.000 claims abstract description 4
- 238000012360 testing method Methods 0.000 claims description 23
- 238000012549 training Methods 0.000 claims description 13
- 238000012545 processing Methods 0.000 claims description 7
- 238000010801 machine learning Methods 0.000 abstract description 11
- 238000011160 research Methods 0.000 abstract description 4
- 230000009286 beneficial effect Effects 0.000 abstract description 2
- 238000013528 artificial neural network Methods 0.000 description 18
- 238000012706 support-vector machine Methods 0.000 description 16
- 230000000694 effects Effects 0.000 description 7
- 230000033558 biomineral tissue development Effects 0.000 description 2
- 229910052500 inorganic mineral Inorganic materials 0.000 description 2
- 239000011707 mineral Substances 0.000 description 2
- 230000004075 alteration Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000001360 synchronised effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/21—Design or setup of recognition systems or techniques; Extraction of features in feature space; Blind source separation
- G06F18/214—Generating training patterns; Bootstrap methods, e.g. bagging or boosting
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/24—Classification techniques
- G06F18/241—Classification techniques relating to the classification model, e.g. parametric or non-parametric approaches
- G06F18/2411—Classification techniques relating to the classification model, e.g. parametric or non-parametric approaches based on the proximity to a decision surface, e.g. support vector machines
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
- G06T17/05—Geographic models
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Data Mining & Analysis (AREA)
- General Physics & Mathematics (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Evolutionary Biology (AREA)
- Evolutionary Computation (AREA)
- Bioinformatics & Computational Biology (AREA)
- General Engineering & Computer Science (AREA)
- Artificial Intelligence (AREA)
- Life Sciences & Earth Sciences (AREA)
- Geometry (AREA)
- Software Systems (AREA)
- Remote Sensing (AREA)
- Computer Graphics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明公开了一种基于地球物理数据使用堆叠法的三维地质建模方法，其特征在于，处理方法包括：本方法包括Stacking集成方法，所述Stacking集成方法包括三维特征提取方法和数据预处理方法，其中所述三维特征提取方法包括三维钻孔模型，通过所述三维钻孔模型分析得到的岩性样本，在通过岩性样本分析的结果得出不同的岩石种类，在通过这些不同的岩石种类得出玄武岩、花岗岩或者其他的岩性；本发明的有益效果：本方法基于机器学习中的Stacking模型，融合了RF,SVM,GBDT,NN,XGBoost等机器学习算法，基于原始的地质数据(三维坐标、岩性等)与地球物理数据(磁化率、剩余密度等)，训练出一个能够对研究区进行三维岩性分类的模型，以达到快速三维地质建模的目的。
Description
技术领域
本发明涉及矿产资源勘测领域，更具体地说，属于一种基于地球物理数据使用堆叠法的三维地质建模方法。
背景技术
近些年机器学习方法发展迅速，当模型的输入特征维度较高并且与岩性的映射关系为非线性的时候，机器学习算法比传统的地质统计学具有更好的预测效果。近些年来许多研究着重于将一些算法应用在三维地质建模的某个领域，如基于机器学习的钻孔数据隐式三维地质建模方法，基于机器学习与可靠度算法的围岩动态分级方法等。要使三维地质模型在输入变量较多的情况下相对准确，必须建立一种综合单个分类器的集成分类器来提高岩性分类的准确度，使用机器学习方法建立三维地质模型，以期减少工作量，提升模型准确率。
故本方法基于机器学习中的Stacking模型，融合了RF,SVM,GBDT,NN, XGBoost等机器学习算法，基于原始的地质数据(三维坐标、岩性等)与地球物理数据(磁化率、剩余密度等)，训练出一个能够对研究区进行三维岩性分类的模型，以达到快速三维地质建模的目的。
发明内容
针对现有应用及技术存在的不足，本发明提供了一种基于地球物理数据使用堆叠法的三维地质建模方法。该方法能够很好地解决传统方法中效率和模型准确性等问题。基于此方案，用户可以通过电脑编程输入程序，因此，用户可以放心的借助电脑终端进行同步操作，对于矿产资源勘查的准确性和效率提升问题具有十分重要的意义。
为实现上述目的，本发明提供了如下技术方案：
一种基于地球物理数据使用堆叠法的三维地质建模方法，处理方法包括：本方法包括Stacking集成方法，所述Stacking集成方法包括三维特征提取方法和数据预处理方法，其中所述三维特征提取方法包括三维钻孔模型，通过所述三维钻孔模型分析得到的岩性样本，在通过岩性样本分析的结果得出不同的岩石种类，在通过这些不同的岩石种类得出玄武岩、花岗岩或者其他的岩性，然后通过这些不同岩性分类组成三维密度模型和三维磁化率模型；所述数据预处理方法包括原始数据集，通过所述原始数据集得出标准化数据、处理数据类别不平衡和数据划分训练测试集，在通过所述所述原始数据集得出标准化数据、处理数据类别不平衡和数据划分训练测试集得出最终的数据集。
进一步的，通过所述三维密度模型和三维磁化率模型可以得出所述原始数据集。
进一步的，通过所述原始数据集分别得出SVM、GBDT、NN、RF、XGboost, 通过所述SVM、GBDT、NN、RF、XGboost组成测试模型。
进一步的，所述测试模型包括测试集和训练集，所述测试集和训练集得出RF、SVM、NN、RF、XGboost。
本发明的有益效果：
本发明通过一种基于地球物理数据使用堆叠法的三维地质建模方法解决了原有的地球物理数据勘探技术问题：
本方法基于机器学习中的Stacking模型，融合了RF,SVM,GBDT,NN, XGBoost等机器学习算法，基于原始的地质数据(三维坐标、岩性等)与地球物理数据(磁化率、剩余密度等)，训练出一个能够对研究区进行三维岩性分类的模型，以达到快速三维地质建模的目的。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1为本发明的技术流程图。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
如附图1所示的一种基于地球物理数据使用堆叠法的三维地质建模方法，处理方法包括：本方法包括Stacking集成方法，所述Stacking集成方法包括三维特征提取方法和数据预处理方法，其中所述三维特征提取方法包括三维钻孔模型，通过所述三维钻孔模型分析得到的岩性样本，在通过岩性样本分析的结果得出不同的岩石种类，在通过这些不同的岩石种类得出玄武岩、花岗岩或者其他的岩性，然后通过这些不同岩性分类组成三维密度模型和三维磁化率模型；所述数据预处理方法包括原始数据集，通过所述原始数据集得出标准化数据、处理数据类别不平衡和数据划分训练测试集，在通过所述所述原始数据集得出标准化数据、处理数据类别不平衡和数据划分训练测试集得出最终的数据集。
优选的，通过所述三维密度模型和三维磁化率模型可以得出所述原始数据集。
优选的，通过所述原始数据集分别得出SVM、GBDT、NN、RF、XGboost, 通过所述SVM、GBDT、NN、RF、XGboost组成测试模型。
优选的，所述测试模型包括测试集和训练集，所述测试集和训练集得出 RF、SVM、NN、RF、XGboost。
优选的，最终的Stacking集成学习器在测试集上的泛化效果总体上来说是优于其它五个子分类器的。其中五个子分类器中，支持向量机(SVM)与神经网络(NN)性能比较差，神经网络在花岗岩、玄武岩与其它岩性上的F1分数分别为0.56、0.77与0.61。梯度提升树(GBDT)与XGboost在单分类器中泛化性能比较好，XGboost在花岗岩、玄武岩与其它岩性上的F1分数分别为 0.80、0.89与0.83，而且在花岗岩与其它岩性类别上甚至优于stacking集成模型。对于花岗岩来说，stacking集成分类器的F1分数为0.80与梯度提升树(GBDT)和XGboost的分数相当，但是召回率达到了0.83；stacking集成学习器在花岗岩上的综合分数比较高，其中精确率、召回率与F1分数分别达到了0.85、0.96与0.90；在其它岩性类别上，Stacking集成方法的精确率也达到了0.92。以上的几个性能指标综合反映了最终的集成模型取得了比单个分类器更好的泛化效果，并且在与成矿密切相关的花岗岩与玄武岩类别上的分类与预测效果最为理想。
本发明通过一种基于地球物理数据使用堆叠法的三维地质建模方法解决了原有的地球物理重数据勘探技术问题：
本方法包括Stacking集成方法，所述Stacking集成方法包括三维特征提取方法和数据预处理方法，其中所述三维特征提取方法包括三维钻孔模型，通过所述三维钻孔模型分析得到的岩性样本，在通过岩性样本分析的结果得出不同的岩石种类，在通过这些不同的岩石种类得出玄武岩、花岗岩或者其他的岩性，然后通过这些不同岩性分类组成三维密度模型和三维磁化率模型；所述数据预处理方法包括原始数据集，通过所述原始数据集得出标准化数据、处理数据类别不平衡和数据划分训练测试集，在通过所述所述原始数据集得出标准化数据、处理数据类别不平衡和数据划分训练测试集得出最终的数据集，通过所述三维密度模型和三维磁化率模型可以得出所述原始数据集，通过所述原始数据集分别得出SVM、GBDT、NN、RF、XGboost,通过所述SVM、GBDT、 NN、RF、XGboost组成测试模型，进一步的，所述测试模型包括测试集和训练集，所述测试集和训练集得出RF、SVM、NN、RF、XGboost。
最终的Stacking集成学习器在测试集上的泛化效果总体上来说是优于其它五个子分类器的。其中五个子分类器中，支持向量机(SVM)与神经网络(NN) 性能比较差，神经网络在花岗岩、玄武岩与其它岩性上的F1分数分别为0.56、 0.77与0.61。梯度提升树(GBDT)与XGboost在单分类器中泛化性能比较好， XGboost在花岗岩、玄武岩与其它岩性上的F1分数分别为0.80、0.89与0.83，而且在花岗岩与其它岩性类别上甚至优于stacking集成模型。对于花岗岩来说，stacking集成分类器的F1分数为0.80与梯度提升树(GBDT)和XGboost 的分数相当，但是召回率达到了0.83；stacking集成学习器在花岗岩上的综合分数比较高，其中精确率、召回率与F1分数分别达到了0.85、0.96与0.90；在其它岩性类别上，Stacking集成方法的精确率也达到了0.92。以上的几个性能指标综合反映了最终的集成模型取得了比单个分类器更好的泛化效果，并且在与成矿密切相关的花岗岩与玄武岩类别上的分类与预测效果最为理想。
在本发明创造的描述中，需要说明的是，除非另有明确的规定和限定，术语“安装”、“相连”、“连接”应做广义理解，例如，可以是固定连接，也可以是可拆卸连接，或一体地连接；可以是机械连接，也可以是电连接；可以是直接相连，也可以通过中间媒介间接相连，可以是两个元件内部的连通。对于本领域的普通技术人员而言，可以通过具体情况理解上述术语在本发明创造中的具体含义。
需要说明的是，在本文中，术语“上”、“下”、“内”、“外”“前端”、“后端”、“两端”、“一端”、“另一端”等指示的方位或位置关系为基于附图所示的方位或位置关系，仅是为了便于描述本发明和简化描述，而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作，因此不能理解为对本发明的限制。
尽管已经示出和描述了本发明的实施例，对于本领域的普通技术人员而言，可以理解在不脱离本发明的原理和精神的情况下可以对这些实施例进行多种变化、修改、替换和变型，本发明的范围由所附权利要求及其等同物限定。
Claims (4)
1.一种基于地球物理数据使用堆叠法的三维地质建模方法，其特征在于，处理方法包括：本方法包括Stacking集成方法，所述Stacking集成方法包括三维特征提取方法和数据预处理方法，其中所述三维特征提取方法包括三维钻孔模型，通过所述三维钻孔模型分析得到的岩性样本，在通过岩性样本分析的结果得出不同的岩石种类，在通过这些不同的岩石种类得出玄武岩、花岗岩或者其他的岩性，然后通过这些不同岩性分类组成三维密度模型和三维磁化率模型；所述数据预处理方法包括原始数据集，通过所述原始数据集得出标准化数据、处理数据类别不平衡和数据划分训练测试集，在通过所述原始数据集得出标准化数据、处理数据类别不平衡和数据划分训练测试集得出最终的数据集。
2.根据权利要求1所述的一种基于地球物理数据使用堆叠法的三维地质建模方法，其特征在于：通过所述三维密度模型和三维磁化率模型可以得出所述原始数据集。
3.根据权利要求1所述的一种基于地球物理数据使用堆叠法的三维地质建模方法，其特征在于：通过所述原始数据集分别得出SVM、GBDT、NN、RF、XGboost,通过所述SVM、GBDT、NN、RF、XGboost组成测试模型。
4.根据权利要求1所述的一种基于地球物理数据使用堆叠法的三维地质建模方法，其特征在于：所述测试模型包括测试集和训练集，所述测试集和训练集得出RF、SVM、NN、RF、XGboost。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110353243.3A CN113111927A (zh) | 2021-04-01 | 2021-04-01 | 一种基于地球物理数据使用堆叠法的三维地质建模方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110353243.3A CN113111927A (zh) | 2021-04-01 | 2021-04-01 | 一种基于地球物理数据使用堆叠法的三维地质建模方法 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN113111927A true CN113111927A (zh) | 2021-07-13 |
Family
ID=76713423
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110353243.3A Pending CN113111927A (zh) | 2021-04-01 | 2021-04-01 | 一种基于地球物理数据使用堆叠法的三维地质建模方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113111927A (zh) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115937568A (zh) * | 2022-09-29 | 2023-04-07 | 中国地质大学（北京） | 一种玄武岩构造背景分类方法、系统、装置及存储介质 |
| CN116502145A (zh) * | 2023-03-15 | 2023-07-28 | 西南交通大学 | 基于改进ER规则的Stacking集成学习围岩分级方法、装置及介质 |
-
2021
- 2021-04-01 CN CN202110353243.3A patent/CN113111927A/zh active Pending
Non-Patent Citations (1)
| Title |
|---|
| JIA R等: "A stacking methodology of machine learning for 3D geological modeling with geological-geophysical datasets, Laochang Sn camp, Gejiu (China)", 《COMPUTERS & GEOSCIENCES》 * |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115937568A (zh) * | 2022-09-29 | 2023-04-07 | 中国地质大学（北京） | 一种玄武岩构造背景分类方法、系统、装置及存储介质 |
| CN115937568B (zh) * | 2022-09-29 | 2024-05-07 | 中国地质大学（北京） | 一种玄武岩构造背景分类方法、系统、装置及存储介质 |
| CN116502145A (zh) * | 2023-03-15 | 2023-07-28 | 西南交通大学 | 基于改进ER规则的Stacking集成学习围岩分级方法、装置及介质 |
| CN116502145B (zh) * | 2023-03-15 | 2024-01-30 | 西南交通大学 | 基于改进ER规则的Stacking集成学习围岩分级方法、装置及介质 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US11016214B2 (en) | Dolomite reservoir prediction method and system based on well and seismic combination, and storage medium | |
| CN103077720B (zh) | 一种说话人识别方法及系统 | |
| CN113111927A (zh) | 一种基于地球物理数据使用堆叠法的三维地质建模方法 | |
| CN111596978A (zh) | 用人工智能进行岩相分类的网页显示方法、模块和系统 | |
| CN104318109A (zh) | 基于支持向量机的页岩气储层识别方法 | |
| CN104047598A (zh) | 非均质古岩溶碳酸盐岩储层产能预测方法 | |
| CN110344824B (zh) | 一种基于随机森林回归的声波曲线生成方法 | |
| Zhu et al. | Rapid identification of high-quality marine shale gas reservoirs based on the oversampling method and random forest algorithm | |
| CN108952699A (zh) | 一种复杂地质钻进过程地层岩性智能识别方法 | |
| CN109388816A (zh) | 一种复杂岩性的分级识别方法 | |
| CN110245686A (zh) | 一种计算石英相对百分含量的岩性快速识别方法 | |
| CN103376468A (zh) | 基于神经网络函数逼近算法的储层参数定量表征方法 | |
| CN105277979B (zh) | 一种地震属性的优化方法和装置 | |
| CN109298464A (zh) | 致密砂岩储层成岩相测井识别方法及装置 | |
| CN107885719A (zh) | 基于人工智能的词汇类别挖掘方法、装置及存储介质 | |
| Guoqiang et al. | Construction of well logging knowledge graph and intelligent identification method of hydrocarbon-bearing formation | |
| Yin et al. | FMI image based rock structure classification using classifier combination | |
| CN115330007A (zh) | 基于深度融合模型的储层多分类预测方法 | |
| CN110907996B (zh) | 一种致密气藏自动识别方法 | |
| CN110956388A (zh) | 一种海上油气藏增产方案生成方法 | |
| CN112711068A (zh) | 一种砂岩中油气有效储层预测方法及装置 | |
| CN115808710A (zh) | 一种致密储层叠前地震数据特征提取方法 | |
| Ma | The Research of Stock Predictive Model based on the Combination of CART and DBSCAN | |
| CN104314558B (zh) | 用斯通利波能量损失度判别储层有效性的方法 | |
| CN116108392A (zh) | 一种基于改进随机森林算法的地质构造识别技术 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| RJ01 | Rejection of invention patent application after publication | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20210713 |