CN104899928A - 基于稀疏钻孔点的三维地质建模方法 - Google Patents

基于稀疏钻孔点的三维地质建模方法 Download PDF

Info

Publication number
CN104899928A
CN104899928A CN201510254716.9A CN201510254716A CN104899928A CN 104899928 A CN104899928 A CN 104899928A CN 201510254716 A CN201510254716 A CN 201510254716A CN 104899928 A CN104899928 A CN 104899928A
Authority
CN
China
Prior art keywords
model
borehole
data
models
stratum
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.)
Granted
Application number
CN201510254716.9A
Other languages
English (en)
Other versions
CN104899928B (zh
Inventor
陈伟
王惠
蔡斌
王珂
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China University of Mining and Technology CUMT
Original Assignee
China University of Mining and Technology CUMT
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by China University of Mining and Technology CUMT filed Critical China University of Mining and Technology CUMT
Priority to CN201510254716.9A priority Critical patent/CN104899928B/zh
Publication of CN104899928A publication Critical patent/CN104899928A/zh
Application granted granted Critical
Publication of CN104899928B publication Critical patent/CN104899928B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Geophysics And Detection Of Objects (AREA)
  • Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)

Abstract

一种基于稀疏钻孔点的三维地质建模方法,属于煤矿地质的三维地质建模方法。建模方法:1、构建钻孔数据库,提取相应底层信息,产生完整钻层与不完整钻层的对比图,使用固定点的插值算法;2、基于勘探线剖面图,根据地层走势的勘探资料,设定横向虚拟钻孔的标记点,并以此坐标作为插值的输入参数形成整个钻孔;3.根据已构建的模型,运用交叉检验作为一种统计思路检验模型的精确度;4.以某一地层的等值线图为参照,对比模型的相应地层数据,统计相似性,直到模型修改到合理的可信度区间。模型构建过程在稀疏钻孔点的大背景下,做到了融入其他非钻孔资料,虚拟缺失钻孔数据,针对可能构建模型的误差,判别突变构造方法,使模型的构建质量提升。

Description

基于稀疏钻孔点的三维地质建模方法
技术领域
本发明涉及一种煤矿地质的三维地质建模方法,特别是一种基于稀疏钻孔点的三维地质建模方法。
背景技术
地质体的真三维建模成为地质研究的热门领域,尤其是矿山地质体。基于钻孔数据的三维地质建模技术最为常见,三维地质建模较为直接、有效、完整的反映了地质地层信息。但是,由于历史原因,传统的钻孔数据相对稀少且勘探的目的迥异,三维地质建模的数据质量不足,而重新设计勘探钻孔又耗时耗资过多,所以现有稀疏钻孔点的建模不能满足矿山地质体研究的需要。
发明内容
本发明的目的是要提供一种基于稀疏钻孔点的三维地质建模方法,解决现有稀疏钻孔点的建模不能满足矿山地质体研究的需要的问题。
本发明的目的是这样实现的:该三维地质建模方法:
1、构建钻孔数据库,提取相应底层信息,产生完整钻层与不完整钻层的对比图,并规划在缺失的钻孔下纵向虚拟延伸,具体是使用设定点的插值算法;
2、基于勘探线剖面图,地质勘探线剖面图是地质勘探线剖面空间要素按一定比例尺的缩影,是现有最常见的矿区勘探资料之一,根据地质工程人员研究地层走势的勘探资料,设定横向虚拟钻孔的标记点,并以此坐标作为插值的输入参数形成整个钻孔;
3、根据实测和虚拟的钻孔,利用地质工程人员的地质学知识划分地层,然后关联各层所属的钻孔钻层,然后选取插值算法形成一定粒度的地层模型,进而组成整个地质体模型;
4、针对以上构建的模型,运用交叉检验作为一种统计思路检验模型的精确度,
RMSECV = 1 M - 1 Σ i = 1 N ( y i - y ^ i ) 2 , RMSEP = 1 M - 1 Σ i = 1 M ( y i - y ^ i ) 2
其中,为通过模型预测的结果,yi为用标准方法测定的结果,ym为yi的平均值,M为建模用的样品数;交叉验证均方根误差(RMSECV)和预测均方根误差(RMSEP)也是一个整体评价指标,是检验模型稳定性的一个统计指标;当RMSEP远大于RMSECV时,建模样品的代表性差、模型信息拟和不够或过拟和;当RMSEP远小于RMSECV时,验证样品代表性差;
5、以某一地层的等值线图为参照,对比模型的相应地层数据,统计相似性,并结合步骤4中的结果,推断相应位置是否具有突变构造,如果有则修正模型的源数据区,并进行相同步骤的再次验证,直到模型修改到合理的可信度区间。
上述方案针对处理稀疏钻孔点造成的数据质量不足,进而模型可信度低的问题,采用两种思维进行处理,一种是基于正向思维,即通过利用非钻孔数据或者利用合理的虚拟和插值方法弥补源数据不足;另外一种是基于反向思维,即在构建好的模型上通过交叉检验或者异构数据源的统计思路去修正模型,以提高建模质量。解决了现有稀疏钻孔点的建模不能满足矿山地质体研究的需要的问题,达到了本发明的目的。
主要优点如下:
1.模型构建过程在稀疏钻孔点的大背景下,做到了融入其他非钻孔资料,并通过一定方式合理地虚拟了缺失钻孔数据。
2.针对可能构建模型的误差,不依赖旧有的勘探资料,提出了基于统计的判别突变构造方法,使模型的构建质量提升。
附图说明:
图1是本发明的完整钻层与不完整钻层对比图。
图2是本发明的横向钻孔扩展示意图。
具体实施方式
实施例1:该三维地质建模方法:
1、构建钻孔数据库,提取相应底层信息,产生完整钻层与不完整钻层的对比图,并规划在缺失的钻孔下纵向虚拟延伸,具体是使用设定点的插值算法;
2、基于勘探线剖面图,地质勘探线剖面图是地质勘探线剖面空间要素按一定比例尺的缩影,是现有最常见的矿区勘探资料之一,根据地质工程人员研究地层走势的勘探资料,设定横向虚拟钻孔的标记点,并以此坐标作为插值的输入参数形成整个钻孔;
3、根据实测和虚拟的钻孔,利用地质工程人员的地质学知识划分地层,然后关联各层所属的钻孔钻层,然后选取插值算法形成一定粒度的地层模型,进而组成成整个地质体模型;
4、针对以上构建的模型,运用交叉检验作为一种统计思路检验模型的精确度,
RMSECV = 1 M - 1 Σ i = 1 N ( y i - y ^ i ) 2 , RMSEP = 1 M - 1 Σ i = 1 M ( y i - y ^ i ) 2
其中,为通过模型预测的结果,yi为用标准方法测定的结果,ym为yi的平均值,M为建模用的样品数;交叉验证均方根误差(RMSECV)和预测均方根误差(RMSEP)也是一个整体评价指标,是检验模型稳定性的一个统计指标;当RMSEP远大于RMSECV时,建模样品的代表性差、模型信息拟和不够或过拟和;当RMSEP远小于RMSECV时,验证样品代表性差;
5、以某一地层的等值线图为参照,对比模型的相应地层数据,统计相似性,并结合步骤4中的结果,推断相应位置是否具有突变构造,如果有则修正模型的源数据区,并进行相同步骤的再次验证,直到模型修改到合理的可信度区间。
在采用该方法时,先要对钻孔数据库信息进行整理。基于整理的信息设定需要纵向虚拟延长的钻孔,并根据剖面数据和地层数据,设计横向的虚拟钻孔。最后基于交叉或者非钻孔资料,用基于统计的突变地质构造的判别不断修正完善模型的可信度低的区域。
图1为完整钻层与不完整钻层对比图;图中梯形线是最深地层数(左刻度),呈递增(因为数据集是以层数排序的)。锯齿线为其最深下底到整个模型下界的距离差(右刻度),整体呈下降趋势。最下一行数值(下刻度)是钻孔的标识号,即每个钻孔有最深层数和下界距离两个数据。
图2为横向钻孔扩展示意;图中方点为各地层与钻孔的交叉点,或者说标记点。虚垂线为虚拟钻孔,实垂线为实际钻孔,A点是预估的地层分界点,作为虚拟钻孔的插值信息输入点。

Claims (1)

1.一种基于稀疏钻孔点的三维地质建模方法,其特征是:
(1)、构建钻孔数据库,提取相应底层信息,产生完整钻层与不完整钻层的对比图,并规划在缺失的钻孔下纵向虚拟延伸,具体是使用设定点的插值算法;
(2)、基于勘探线剖面图,地质勘探线剖面图是地质勘探线剖面空间要素按一定比例尺的缩影,是现有最常见的矿区勘探资料之一,根据地质工程人员研究地层走势的勘探资料,设定横向虚拟钻孔的标记点,并以此坐标作为插值的输入参数形成整个钻孔;
(3)、根据实测和虚拟的钻孔,利用地质工程人员的地质学知识划分地层,然后关联各层所属的钻孔钻层,然后选取插值算法形成一定粒度的地层模型,进而组成成整个地质体模型。
(4)、针对以上构建的模型,运用交叉检验作为一种统计思路检验模型的精确度,
RMSECV = 1 M - 1 Σ i = 1 N ( y i - y ^ i ) 2 , RMSEP = 1 M - 1 Σ i = 1 M ( y i - y ^ i ) 2
其中,为通过模型预测的结果,yi为用标准方法测定的结果,ym为yi的平均值,M为建模用的样品数;交叉验证均方根误差(RMSECV)和预测均方根误差(RMSEP)也是一个整体评价指标,是检验模型稳定性的一个统计指标;当RMSEP远大于RMSECV时,建模样品的代表性差、模型信息拟和不够或过拟和;当RMSEP远小于RMSECV时,验证样品代表性差;
(5)、以某一地层的等值线图为参照,对比模型的相应地层数据,统计相似性,并结合步骤4中的结果,推断相应位置是否具有突变构造,如果有则修正模型的源数据区,并进行相同步骤的再次验证,直到模型修改到合理的可信度区间。
CN201510254716.9A 2015-05-18 2015-05-18 基于稀疏钻孔点的三维地质建模方法 Active CN104899928B (zh)

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 true CN104899928A (zh) 2015-09-09
CN104899928B 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)

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106296817A (zh) * 2016-08-08 2017-01-04 西安科技大学 一种基于地质数据的工作面煤层三维建模方法
CN106558100A (zh) * 2016-10-25 2017-04-05 中国电建集团华东勘测设计研究院有限公司 一种基于钻孔数据的地层模型自动建模方法
CN106709990A (zh) * 2016-11-25 2017-05-24 北京市地震局 适用于地质结构的三维模型建立的方法和装置
CN112307541A (zh) * 2020-10-27 2021-02-02 深圳市勘察研究院有限公司 一种城市地下空间岩土信息化综合集成数字交付方法
CN112562080A (zh) * 2021-01-07 2021-03-26 南京市测绘勘察研究院股份有限公司 一种基于钻孔数据的地质结构降维模型建模方法
CN113808265A (zh) * 2021-08-12 2021-12-17 贵州省地质调查院(贵州省地质矿产勘查开发局地质科学研究所) 大尺度区域三维地质模型建模剖面的构建方法
CN115619788A (zh) * 2022-12-19 2023-01-17 武汉市测绘研究院 一种三维地质模型质量的自动化定量评价方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101727678A (zh) * 2008-10-20 2010-06-09 车德福 复杂褶皱真三维建模方法
CN101882171A (zh) * 2010-05-21 2010-11-10 中交第二公路勘察设计研究院有限公司 一种交互式隧道及围岩体三维模型快速构建方法
US20140152659A1 (en) * 2012-12-03 2014-06-05 Preston H. Davidson Geoscience data visualization and immersion experience
CN103886641A (zh) * 2014-04-18 2014-06-25 重庆市勘测院 山地城市区域地质三维模型构建集成方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101727678A (zh) * 2008-10-20 2010-06-09 车德福 复杂褶皱真三维建模方法
CN101882171A (zh) * 2010-05-21 2010-11-10 中交第二公路勘察设计研究院有限公司 一种交互式隧道及围岩体三维模型快速构建方法
US20140152659A1 (en) * 2012-12-03 2014-06-05 Preston H. Davidson Geoscience data visualization and immersion experience
CN103886641A (zh) * 2014-04-18 2014-06-25 重庆市勘测院 山地城市区域地质三维模型构建集成方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
宋永康等: "基于交叉检验法的地质预测模型优选", 《石油化工应用》 *
王媛妮等: "基于稀疏钻孔点的地层建模及可视化", 《测绘信息与工程》 *
陈兰珍、叶志华、赵静: "《蜂蜜近红外光谱检测技术》", 31 December 2012, 北京:中国轻工业出版社 *

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106296817A (zh) * 2016-08-08 2017-01-04 西安科技大学 一种基于地质数据的工作面煤层三维建模方法
CN106296817B (zh) * 2016-08-08 2017-07-07 西安科技大学 一种基于地质数据的工作面煤层三维建模方法
CN106558100A (zh) * 2016-10-25 2017-04-05 中国电建集团华东勘测设计研究院有限公司 一种基于钻孔数据的地层模型自动建模方法
CN106558100B (zh) * 2016-10-25 2019-08-30 中国电建集团华东勘测设计研究院有限公司 一种基于钻孔数据的地层模型自动建模方法
CN106709990A (zh) * 2016-11-25 2017-05-24 北京市地震局 适用于地质结构的三维模型建立的方法和装置
CN112307541A (zh) * 2020-10-27 2021-02-02 深圳市勘察研究院有限公司 一种城市地下空间岩土信息化综合集成数字交付方法
CN112562080A (zh) * 2021-01-07 2021-03-26 南京市测绘勘察研究院股份有限公司 一种基于钻孔数据的地质结构降维模型建模方法
CN112562080B (zh) * 2021-01-07 2023-06-02 南京市测绘勘察研究院股份有限公司 一种基于钻孔数据的地质结构降维模型建模方法
CN113808265A (zh) * 2021-08-12 2021-12-17 贵州省地质调查院(贵州省地质矿产勘查开发局地质科学研究所) 大尺度区域三维地质模型建模剖面的构建方法
CN115619788A (zh) * 2022-12-19 2023-01-17 武汉市测绘研究院 一种三维地质模型质量的自动化定量评价方法

Also Published As

Publication number Publication date
CN104899928B (zh) 2018-12-04

Similar Documents

Publication Publication Date Title
CN104899928A (zh) 基于稀疏钻孔点的三维地质建模方法
US10288544B2 (en) Method for characterizing the fracture network of a fractured reservoir and method for exploiting it
CN105785471A (zh) 一种矿井预开采煤层的冲击危险性评价方法
Close et al. Integrated workflows for shale gas and case study results for the Horn River Basin, British Columbia, Canada
US10641923B2 (en) Method for characterizing and exploiting a subterranean formation comprising a network of fractures
CN104992468A (zh) 缝洞型碳酸盐岩油气藏三维地质建模方法
CN104636980B (zh) 针对河道砂油藏类型油气汇集条件的地球物理表征方法
CN105005097A (zh) 利用重力、磁力、电磁、地震资料综合识别火成岩方法
CN103114840A (zh) 一种高-过高成熟页岩有机碳含量计算方法及装置
CN104749617A (zh) 一种多尺度裂缝储层正演模型建立方法
CN105223612A (zh) 一种基于地震信息的煤矿水害预测评价方法
CN103088803B (zh) 鉴别岩溶塌陷致塌因素的方法
Harris The role of geology in reservoir simulation studies
CN104181595A (zh) 用于复杂拉张构造体系断层伴生裂缝定量识别的新方法
CN105467438A (zh) 一种基于三模量的泥页岩地应力三维地震表征方法
CN107798189A (zh) 一种基于透明空间地球物理的精准动态防突方法
CN106033127B (zh) 基于横波速度变化率的地应力方位地震预测方法
CN102288996A (zh) 叠前afi反演定量预测含气概率的方法
CN114357750A (zh) 一种采空区充水状态评估方法
CN106154342B (zh) 一种确定洞穴充填物电阻率的方法
EP2661646B1 (en) Method and system of calculating a fault throw
CN111695228A (zh) 一种缝洞型碳酸盐岩油藏的多尺度裂缝建模方法
CN109358364B (zh) 一种建立地下暗河储集体地质模型的方法、装置及系统
Wilson et al. Fracture model of the Upper Freeport coal: Marshall County West Virginia pilot ECBMR and CO2 sequestration site
CN104167022A (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