CN104899928B - 基于稀疏钻孔点的三维地质建模方法 - Google Patents
基于稀疏钻孔点的三维地质建模方法 Download PDFInfo
- Publication number
- CN104899928B CN104899928B CN201510254716.9A CN201510254716A CN104899928B CN 104899928 B CN104899928 B CN 104899928B CN 201510254716 A CN201510254716 A CN 201510254716A CN 104899928 B CN104899928 B CN 104899928B
- Authority
- CN
- China
- Prior art keywords
- model
- drilling
- data
- stratum
- geological
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 17
- 238000005553 drilling Methods 0.000 claims abstract description 41
- 230000035772 mutation Effects 0.000 claims abstract description 8
- 238000012360 testing method Methods 0.000 claims abstract description 7
- 238000005065 mining Methods 0.000 claims abstract description 3
- 238000005259 measurement Methods 0.000 claims description 5
- 238000011160 research Methods 0.000 claims description 5
- 238000002790 cross-validation Methods 0.000 claims description 3
- 238000011156 evaluation Methods 0.000 claims description 3
- 238000010561 standard procedure Methods 0.000 claims description 3
- 238000007689 inspection Methods 0.000 claims 1
- 230000015572 biosynthetic process Effects 0.000 abstract description 3
- 238000010276 construction Methods 0.000 abstract description 3
- 239000003245 coal Substances 0.000 abstract description 2
- 241001074085 Scophthalmus aquosus Species 0.000 abstract 1
- 238000010586 diagram Methods 0.000 abstract 1
- 239000000284 extract Substances 0.000 abstract 1
- 238000012986 modification Methods 0.000 abstract 1
- 230000004048 modification Effects 0.000 abstract 1
- 238000005755 formation reaction Methods 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 230000007812 deficiency Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 229910052500 inorganic mineral Inorganic materials 0.000 description 1
- 239000011707 mineral Substances 0.000 description 1
Landscapes
- Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
一种基于稀疏钻孔点的三维地质建模方法，属于煤矿地质的三维地质建模方法。建模方法：1、构建钻孔数据库，提取相应底层信息，产生完整钻层与不完整钻层的对比图，使用固定点的插值算法；2、基于勘探线剖面图，根据地层走势的勘探资料，设定横向虚拟钻孔的标记点，并以此坐标作为插值的输入参数形成整个钻孔；3.根据已构建的模型，运用交叉检验作为一种统计思路检验模型的精确度；4.以某一地层的等值线图为参照，对比模型的相应地层数据，统计相似性，直到模型修改到合理的可信度区间。模型构建过程在稀疏钻孔点的大背景下，做到了融入其他非钻孔资料，虚拟缺失钻孔数据，针对可能构建模型的误差，判别突变构造方法，使模型的构建质量提升。
Description
技术领域
本发明涉及一种煤矿地质的三维地质建模方法，特别是一种基于稀疏钻孔点的三维地质建模方法。
背景技术
地质体的真三维建模成为地质研究的热门领域，尤其是矿山地质体。基于钻孔数据的三维地质建模技术最为常见，三维地质建模较为直接、有效、完整的反映了地质地层信息。但是，由于历史原因，传统的钻孔数据相对稀少且勘探的目的迥异，三维地质建模的数据质量不足，而重新设计勘探钻孔又耗时耗资过多，所以现有稀疏钻孔点的建模不能满足矿山地质体研究的需要。
发明内容
本发明的目的是要提供一种基于稀疏钻孔点的三维地质建模方法，解决现有稀疏钻孔点的建模不能满足矿山地质体研究的需要的问题。
本发明的目的是这样实现的：该三维地质建模方法：
1、构建钻孔数据库，提取相应底层信息，产生完整钻层与不完整钻层的对比图，并规划在缺失的钻孔下纵向虚拟延伸，具体是使用设定点的插值算法；
2、基于勘探线剖面图，地质勘探线剖面图是地质勘探线剖面空间要素按一定比例尺的缩影，是现有最常见的矿区勘探资料之一，根据地质工程人员研究地层走势的勘探资料，设定横向虚拟钻孔的标记点，并以标记点的坐标作为插值的输入参数形成整个钻孔；
3、根据实测和虚拟的钻孔，利用地质工程人员的地质学知识划分地层，然后关联各层所属的钻孔钻层，然后选取插值算法形成一定粒度的地层模型，进而组成整个地质体模型；
4、针对以上构建的模型，运用交叉检验作为一种统计思路检验模型的精确度，
其中，为通过模型预测的结果，yi为用标准方法测定的结果，ym为yi的平均值，M为建模用的样品数；交叉验证均方根误差(RMSECV)和预测均方根误差(RMSEP)也是一个整体评价指标，是检验模型稳定性的一个统计指标；当RMSEP远大于RMSECV时，建模样品的代表性差、模型信息拟和不够或过拟和；当RMSEP远小于RMSECV时,验证样品代表性差；
5、以某一地层的等值线图为参照，对比模型的相应地层数据，统计相似性，并结合步骤4中的结果，推断相应位置是否具有突变构造，如果有则修正模型的源数据区，并进行相同步骤的再次验证，直到模型修改到合理的可信度区间，在采用该方法时，先要对钻孔数据库信息进行整理；基于整理的信息设定需要纵向虚拟延长的钻孔，并根据剖面数据和地层数据，设计横向的虚拟钻孔；最后基于交叉或者非钻孔资料，用基于统计的突变地质构造的判别不断修正完善模型的可信度低的区域。
上述方案针对处理稀疏钻孔点造成的数据质量不足，进而模型可信度低的问题，采用两种思维进行处理，一种是基于正向思维，即通过利用非钻孔数据或者利用合理的虚拟和插值方法弥补源数据不足；另外一种是基于反向思维，即在构建好的模型上通过交叉检验或者异构数据源的统计思路去修正模型，以提高建模质量。解决了现有稀疏钻孔点的建模不能满足矿山地质体研究的需要的问题，达到了本发明的目的。
主要优点如下：
1.模型构建过程在稀疏钻孔点的大背景下，做到了融入其他非钻孔资料，并通过一定方式合理地虚拟了缺失钻孔数据。
2.针对可能构建模型的误差，不依赖旧有的勘探资料，提出了基于统计的判别突变构造方法，使模型的构建质量提升。
附图说明：
图1是本发明的完整钻层与不完整钻层对比图。
图2是本发明的横向钻孔扩展示意图。
具体实施方式
实施例1：该三维地质建模方法：
1、构建钻孔数据库，提取相应底层信息，产生完整钻层与不完整钻层的对比图，并规划在缺失的钻孔下纵向虚拟延伸，具体是使用设定点的插值算法；
2、基于勘探线剖面图,地质勘探线剖面图是地质勘探线剖面空间要素按一定比例尺的缩影，是现有最常见的矿区勘探资料之一，根据地质工程人员研究地层走势的勘探资料，设定横向虚拟钻孔的标记点，并以此坐标作为插值的输入参数形成整个钻孔；
3、根据实测和虚拟的钻孔，利用地质工程人员的地质学知识划分地层，然后关联各层所属的钻孔钻层，然后选取插值算法形成一定粒度的地层模型，进而组成整个地质体模型；
4、针对以上构建的模型，运用交叉检验作为一种统计思路检验模型的精确度，
其中，为通过模型预测的结果，yi为用标准方法测定的结果，ym为yi的平均值，M为建模用的样品数；交叉验证均方根误差(RMSECV)和预测均方根误差(RMSEP)也是一个整体评价指标，是检验模型稳定性的一个统计指标；当RMSEP远大于RMSECV时，建模样品的代表性差、模型信息拟和不够或过拟和；当RMSEP远小于RMSECV时,验证样品代表性差；
5、以某一地层的等值线图为参照，对比模型的相应地层数据，统计相似性，并结合步骤4中的结果，推断相应位置是否具有突变构造，如果有则修正模型的源数据区，并进行相同步骤的再次验证，直到模型修改到合理的可信度区间。
在采用该方法时，先要对钻孔数据库信息进行整理；基于整理的信息设定需要纵向虚拟延长的钻孔，并根据剖面数据和地层数据，设计横向的虚拟钻孔；最后基于交叉和非钻孔资料，用基于统计的突变地质构造的判别不断修正完善模型的可信度低的区域。
图1为完整钻层与不完整钻层对比图；图中梯形线是最深地层数(左刻度)，呈递增(因为数据集是以层数排序的)。锯齿线为其最深下底到整个模型下界的距离差(右刻度)，整体呈下降趋势。最下一行数值(下刻度)是钻孔的标识号，即每个钻孔有最深层数和下界距离两个数据。
图2为横向钻孔扩展示意；图中方点为各地层与钻孔的交叉点，或者说标记点。虚垂线为虚拟钻孔，实垂线为实际钻孔，A点是预估的地层分界点，作为虚拟钻孔的插值信息输入点。
Claims (1)
1.一种基于稀疏钻孔点的三维地质建模方法，其特征是：
(1)、构建钻孔数据库，提取相应底层信息，产生完整钻层与不完整钻层的对比图，并规划在缺失的钻孔下纵向虚拟延伸，具体是使用设定点的插值算法；
(2)、基于勘探线剖面图,地质勘探线剖面图是地质勘探线剖面空间要素按一定比例尺的缩影，是现有最常见的矿区勘探资料之一，根据地质工程人员研究地层走势的勘探资料，设定横向虚拟钻孔的标记点，并以标记点的坐标作为插值的输入参数形成整个钻孔；
(3)、根据实测和虚拟的钻孔，利用地质工程人员的地质学知识划分地层，然后关联各层所属的钻孔钻层，然后选取插值算法形成一定粒度的地层模型，进而组成整个地质体模型；
(4)、针对以上构建的模型，运用交叉检验作为一种统计思路检验模型的精确度，
其中，为通过模型预测的结果，yi为用标准方法测定的结果，M为建模用的样品数；交叉验证均方根误差(RMSECV)和预测均方根误差(RMSEP)是一个整体评价指标，是检验模型稳定性的一个统计指标；当RMSEP远大于RMSECV时，建模样品的代表性差、模型信息拟和不够或过拟和；当RMSEP远小于RMSECV时,验证样品代表性差；
(5)、以某一地层的等值线图为参照，对比模型的相应地层数据，统计相似性，并结合步骤(4)中的结果，推断相应位置是否具有突变构造，如果有则修正模型的源数据区，并进行相同步骤的再次验证，直到模型修改到合理的可信度区间，在采用该方法时，先要对钻孔数据库信息进行整理；基于整理的信息设定需要纵向虚拟延长的钻孔，并根据剖面数据和地层数据，设计横向的虚拟钻孔；最后基于交叉检验或者非钻孔资料，用基于统计的突变地质构造的判别不断修正完善模型的可信度低的区域。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510254716.9A CN104899928B (zh) | 2015-05-18 | 2015-05-18 | 基于稀疏钻孔点的三维地质建模方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510254716.9A CN104899928B (zh) | 2015-05-18 | 2015-05-18 | 基于稀疏钻孔点的三维地质建模方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN104899928A CN104899928A (zh) | 2015-09-09 |
| CN104899928B true CN104899928B (zh) | 2018-12-04 |
Family
ID=54032571
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510254716.9A Active CN104899928B (zh) | 2015-05-18 | 2015-05-18 | 基于稀疏钻孔点的三维地质建模方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104899928B (zh) |
Families Citing this family (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106296817B (zh) * | 2016-08-08 | 2017-07-07 | 西安科技大学 | 一种基于地质数据的工作面煤层三维建模方法 |
| CN106558100B (zh) * | 2016-10-25 | 2019-08-30 | 中国电建集团华东勘测设计研究院有限公司 | 一种基于钻孔数据的地层模型自动建模方法 |
| CN106709990A (zh) * | 2016-11-25 | 2017-05-24 | 北京市地震局 | 适用于地质结构的三维模型建立的方法和装置 |
| CN112307541A (zh) * | 2020-10-27 | 2021-02-02 | 深圳市勘察研究院有限公司 | 一种城市地下空间岩土信息化综合集成数字交付方法 |
| CN112562080B (zh) * | 2021-01-07 | 2023-06-02 | 南京市测绘勘察研究院股份有限公司 | 一种基于钻孔数据的地质结构降维模型建模方法 |
| CN113808265B (zh) * | 2021-08-12 | 2024-06-18 | 贵州省地质调查院(贵州省地质矿产勘查开发局地质科学研究所) | 大尺度区域三维地质模型建模剖面的构建方法 |
| CN115619788B (zh) * | 2022-12-19 | 2023-03-21 | 武汉市测绘研究院 | 一种三维地质模型质量的自动化定量评价方法 |
Family Cites Families (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101727678A (zh) * | 2008-10-20 | 2010-06-09 | 车德福 | 复杂褶皱真三维建模方法 |
| CN101882171B (zh) * | 2010-05-21 | 2012-05-09 | 中交第二公路勘察设计研究院有限公司 | 一种交互式隧道及围岩体三维模型快速构建方法 |
| US20140152659A1 (en) * | 2012-12-03 | 2014-06-05 | Preston H. Davidson | Geoscience data visualization and immersion experience |
| CN103886641B (zh) * | 2014-04-18 | 2017-01-04 | 重庆市勘测院 | 山地城市区域地质三维模型构建集成方法 |
-
2015
- 2015-05-18 CN CN201510254716.9A patent/CN104899928B/zh active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN104899928A (zh) | 2015-09-09 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104899928B (zh) | 基于稀疏钻孔点的三维地质建模方法 | |
| AU2009260453B2 (en) | Heterogeneous earth models for a reservoir field | |
| CN109388817B (zh) | 一种储层裂缝三维建模方法 | |
| US9260948B2 (en) | Multi-level reservoir history matching | |
| US9416630B2 (en) | Grid calculation for fluid flow modeling | |
| US10641923B2 (en) | Method for characterizing and exploiting a subterranean formation comprising a network of fractures | |
| CN106894761B (zh) | 利用时间域地震体的地质导向模型的地质导向方法 | |
| CN104992468A (zh) | 缝洞型碳酸盐岩油气藏三维地质建模方法 | |
| CN110579802B (zh) | 一种天然气水合物储层物性参数的高精度反演方法 | |
| BRPI0820174B1 (pt) | Método para formar um modelo geológico de uma região da terra e método para perfurar um furo em uma região da terra | |
| US9176255B2 (en) | Permeability prediction systems and methods using quadratic discriminant analysis | |
| CN110073246B (zh) | 与质量控制有关的改进的方法 | |
| CA2911107C (en) | Local layer geometry engine with work zone generated from buffer defined relative to a wellbore trajectory | |
| CN103088803B (zh) | 鉴别岩溶塌陷致塌因素的方法 | |
| CN107633556B (zh) | 一种定量获得三维矿床地质模型不确定性的方法 | |
| Ray et al. | Characterizing and modeling natural fracture networks in a tight carbonate reservoir in the Middle East: A methodology | |
| CN117852416A (zh) | 基于数字化地质模型的多模态注浆预控分析方法及系统 | |
| CN111680380A (zh) | 基于地质力学特征空间展布的全三维压裂设计方法 | |
| CN107103165B (zh) | 一种计算油田储层陡倾角剪切裂缝线密度的方法 | |
| Ozkaya | Modeling finite-fracture networks in a partially fractured reservoir in the Middle East | |
| Sarkheil et al. | The fracture network modeling in naturally fractured reservoirs using artificial neural network based on image loges and core measurements | |
| CN113742962B (zh) | 一种基于水平井的页岩储层属性三维建模方法 | |
| CN115494552A (zh) | 一种斜井的测井曲线校正方法及其应用 | |
| Lyu et al. | Classification, modeling and characterization of marine carbonate paleokarst reservoirs in Tahe Oilfield, Tarim Basin, China | |
| CN112394416A (zh) | 非均质断裂控制储层预测方法及装置 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |