CN106780741A - 基于多点产状数据的地质曲面快速拟合插值方法 - Google Patents

基于多点产状数据的地质曲面快速拟合插值方法 Download PDF

Info

Publication number
CN106780741A
CN106780741A CN201710073770.2A CN201710073770A CN106780741A CN 106780741 A CN106780741 A CN 106780741A CN 201710073770 A CN201710073770 A CN 201710073770A CN 106780741 A CN106780741 A CN 106780741A
Authority
CN
China
Prior art keywords
geological
occurrence
data
interpolation
plane
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
CN201710073770.2A
Other languages
English (en)
Other versions
CN106780741B (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.)
Tianjin University
PowerChina Chengdu Engineering Co Ltd
Original Assignee
Tianjin University
PowerChina Chengdu Engineering Co Ltd
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 Tianjin University, PowerChina Chengdu Engineering Co Ltd filed Critical Tianjin University
Priority to CN201710073770.2A priority Critical patent/CN106780741B/zh
Publication of CN106780741A publication Critical patent/CN106780741A/zh
Application granted granted Critical
Publication of CN106780741B publication Critical patent/CN106780741B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • 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

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Geometry (AREA)
  • Software Systems (AREA)
  • Remote Sensing (AREA)
  • Computer Graphics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)
  • Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)

Abstract

本发明涉及一种根据地质勘测数据进行计算机地质建模、模拟和分析的方法，其公开了一种基于多点产状数据的地质曲面快速拟合插值方法，充分利用离散地质点的勘探数据，实现对地质点在空间位置和产状走势上的拟合重构工作。该方法包括：A、根据所输入的地质数据，分析提取其中属性相似的地质点，建立地质数据原始数据点集；B、构建地质插值基准平面；C、对每一个原始数据集中的地质点，根据其空间位置和产状信息构建产状四边形；D、提取每一个产状四边形的四个顶点，建立产状数据点集；E、将原始数据点集与产状数据点集合并，形成地质插值点集；F、以地质插值基准平面和地质插值点集，运用基于平面弹性理论的曲面样条插值法，构建地质曲面。

Description

基于多点产状数据的地质曲面快速拟合插值方法
技术领域
本发明涉及一种根据地质勘测数据进行计算机地质建模、模拟和分析的方法，特别涉及一种基于多点产状数据的地质曲面快速拟合插值方法。
背景技术
工程地质勘探需要处理诸多裂隙、断层等结构面的编录与分析工作，这些结构面多以地质曲面进行表征。在工程地质现场通过勘探数据快速准确地构建这些地质曲面，是进一步分析地质结构面特征、开展地质评估和实施地质救险的重要前提。
高博青[1]等提出了一种新的网格划分方法，而B样条曲面对于地形建模则更适用[2,3,4]，传统的地质曲面逼近拟合方法主要面向大范围工程地质分析，其所使用的地质分析点数量大，且分析过程仅考虑地质点的空间位置，而不考虑产状走势，因此所构建的地质曲面仅能满足空间位置的逼近，而不能反映重要区域的产状走势变化。
因此，基于上述问题，本发明提供了一种基于多点产状数据的地质曲面快速拟合插值方法。
参考文献：
[1]高博青，李铁瑞，潘炜.自由曲面网格结构的网格划分方法[P].发明专利ZL201510777920.9，2015
[2]M·施特拉森布克-克莱克,H·迈尔.基于曲面区域的地形建模[P].发明专利ZL200710196359.0，2007
[3]钟登华，李明超，刘杰，等.水利水电工程地质信息的三维统一模型构建方法[P]，发明专利ZL200610013425.1，2006
[4]张明，厉海涛，崔树鑫.一种高自由度的参数化曲面建模方法[P].发明专利ZL201610362889.7，2016
发明内容
本发明所要解决的技术问题是：提出一种基于多点产状数据的地质曲面快速拟合插值方法，充分利用离散地质点的勘探数据，实现对地质点在空间位置和产状走势上的拟合重构工作。
本发明解决上述技术问题所采用的方案是：
基于多点产状数据的地质曲面快速拟合插值方法，包括以下步骤：
步骤A、根据所输入的地质数据，分析提取其中属性相似的地质点，建立地质数据原始数据点集；
步骤B、以原始数据点中地质点的空间中点为平面根节点，以它们的平均产状为平面方向，构建地质插值基准平面；
步骤C、对每一个原始数据集中的地质点，根据其空间位置和产状信息构建产状四边形；
步骤D、提取每一个产状四边形的四个顶点，建立产状数据点集；
步骤E、将原始数据点集与产状数据点集合并，形成地质插值点集；
步骤F、以步骤B所得的地质插值基准平面和步骤E所得的地质插值点集，运用基于平面弹性理论的曲面样条插值法，构建地质曲面。
作为进一步优化，步骤A中的输入地质数据为指定地质区域内的所有地质编录属性数据，其中包括地质点空间位置、产状走势、属性类别、建模标识。
作为进一步优化，步骤B中的地质插值基准平面，是以地质点集的平均产状方向作为平面方向，并由平均产状计算得出法向量；该基准平面通过所选定的地质点集的空间中心。
作为进一步优化，步骤C中的产状四边形，是以地质点的空间位置为形心，以地质点的产状方向为倾斜方向的四边形；该四边形可通过地质点空间位置和产状方向，并给定四边形边长来求得。
作为进一步优化，步骤F中的基于平面弹性理论的曲面样条插值法，是一类用于构造地质曲面的插值数学模型；该方法将计算域视为无限延伸平面，根据实际情况确定插值基平面，在基平面上一点给予一定的垂直位移时，基平面各处均会随着产生位移。
本发明的有益效果是：
将空间内的离散地质点进行分析拟合，并重构地质曲面；充分利用地质点的空间位置数据和产状走势数据，使重构生成的地质曲面既满足地质点空间位置的逼近，又满足地质点产状走势的匹配。
附图说明
图1为基于多点产状数据的地质曲面快速拟合插值方法流程图。
具体实施方式
本发明旨在提出一种基于多点产状数据的地质曲面快速拟合插值方法，充分利用离散地质点的勘探数据，实现对地质点在空间位置和产状走势上的拟合重构工作。
如图1所示，本发明中的基于多点产状数据的地质曲面快速拟合插值方法包括以下步骤：
步骤A，根据所输入的地质数据，分析提取其中属性相似的地质点，建立地质数据原始数据点集；
本步骤中，输入地质数据为指定地质区域内的所有地质编录属性数据，其中包括地质点空间位置、产状走势、属性类别、建模标识等具体信息。
步骤B，以原始数据点中地质点的空间中点为平面根节点，以它们的平均产状为平面方向，构建地质插值基准平面；
本步骤中，地质插值基准平面，是以地质点集的平均产状方向作为平面方向，并由平均产状计算得出法向量；该基准平面通过所选定的地质点集的空间中心。
步骤C，对每一个原始数据集中的地质点，根据其空间位置和产状信息构建产状四边形；
本步骤中，产状四边形是以地质点的空间位置为形心，以地质点的产状方向为倾斜方向的四边形；该四边形可通过地质点空间位置和产状方向，并给定四边形边长来求得。
步骤D，提取每一个产状四边形的四个顶点，建立产状数据点集；
步骤E，将原始数据点集与产状数据点集合并，形成地质插值点集；
步骤F，以步骤B所得的地质插值基准平面和步骤E所得的地质插值点集，运用基于平面弹性理论的曲面样条插值法，构建地质曲面。
本步骤中，基于平面弹性理论的曲面样条插值法，是一类用于构造地质曲面的插值数学模型；该方法将计算域视为无限延伸平面，根据实际情况确定插值基平面，在基平面上一点给予一定的垂直位移时，基平面各处均会随着产生位移。

Claims (5)

1.基于多点产状数据的地质曲面快速拟合插值方法，其特征在于，包括以下步骤：
步骤A、根据所输入的地质数据，分析提取其中属性相似的地质点，建立地质数据原始数据点集；
步骤B、以原始数据点中地质点的空间中点为平面根节点，以它们的平均产状为平面方向，构建地质插值基准平面；
步骤C、对每一个原始数据集中的地质点，根据其空间位置和产状信息构建产状四边形；
步骤D、提取每一个产状四边形的四个顶点，建立产状数据点集；
步骤E、将原始数据点集与产状数据点集合并，形成地质插值点集；
步骤F、以步骤B所得的地质插值基准平面和步骤E所得的地质插值点集，运用基于平面弹性理论的曲面样条插值法，构建地质曲面。
2.如权利要求1所述的基于多点产状数据的地质曲面快速拟合插值方法，其特征在于，步骤A中的输入地质数据为指定地质区域内的所有地质编录属性数据，其中包括地质点空间位置、产状走势、属性类别、建模标识。
3.如权利要求1所述的基于多点产状数据的地质曲面快速拟合插值方法，步骤B中的地质插值基准平面，是以地质点集的平均产状方向作为平面方向，并由平均产状计算得出法向量；该基准平面通过所选定的地质点集的空间中心。
4.如权利要求1所述的基于多点产状数据的地质曲面快速拟合插值方法，步骤C中的产状四边形，是以地质点的空间位置为形心，以地质点的产状方向为倾斜方向的四边形；该四边形可通过地质点空间位置和产状方向，并给定四边形边长来求得。
5.如权利要求1所述的基于多点产状数据的地质曲面快速拟合插值方法，步骤F中的基于平面弹性理论的曲面样条插值法，是一类用于构造地质曲面的插值数学模型；该方法将计算域视为无限延伸平面，根据实际情况确定插值基平面，在基平面上一点给予一定的垂直位移时，基平面各处均会随着产生位移。
CN201710073770.2A 2017-02-10 2017-02-10 基于多点产状数据的地质曲面快速拟合插值方法 Active CN106780741B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710073770.2A CN106780741B (zh) 2017-02-10 2017-02-10 基于多点产状数据的地质曲面快速拟合插值方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710073770.2A CN106780741B (zh) 2017-02-10 2017-02-10 基于多点产状数据的地质曲面快速拟合插值方法

Publications (2)

Publication Number Publication Date
CN106780741A true CN106780741A (zh) 2017-05-31
CN106780741B CN106780741B (zh) 2020-05-12

Family

ID=58955910

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710073770.2A Active CN106780741B (zh) 2017-02-10 2017-02-10 基于多点产状数据的地质曲面快速拟合插值方法

Country Status (1)

Country Link
CN (1) CN106780741B (zh)

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101950030A (zh) * 2010-08-30 2011-01-19 钟世航 一种用弹性波反射法作隧道施工地质预报确定地质体反射面产状和空间位置的方法
CN103646423A (zh) * 2013-12-24 2014-03-19 中国科学院地质与地球物理研究所 一种三维地质建模的方法及装置
CN103713332A (zh) * 2013-12-26 2014-04-09 广西华锡集团股份有限公司 一种裂隙岩体巷道轴线走向确定方法
CN104183017A (zh) * 2014-07-29 2014-12-03 浙江大学 基于地面三维激光点云实现地质体产状自动提取的方法
WO2015137916A1 (en) * 2014-03-10 2015-09-17 Landmark Graphics Corporation Modeling geologic surfaces using unilateral non-node constraints from neighboring surfaces in the stratigraphic sequence

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101950030A (zh) * 2010-08-30 2011-01-19 钟世航 一种用弹性波反射法作隧道施工地质预报确定地质体反射面产状和空间位置的方法
CN103646423A (zh) * 2013-12-24 2014-03-19 中国科学院地质与地球物理研究所 一种三维地质建模的方法及装置
CN103713332A (zh) * 2013-12-26 2014-04-09 广西华锡集团股份有限公司 一种裂隙岩体巷道轴线走向确定方法
WO2015137916A1 (en) * 2014-03-10 2015-09-17 Landmark Graphics Corporation Modeling geologic surfaces using unilateral non-node constraints from neighboring surfaces in the stratigraphic sequence
CN104183017A (zh) * 2014-07-29 2014-12-03 浙江大学 基于地面三维激光点云实现地质体产状自动提取的方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
周晓琴 等: "基于DSI 技术的高精度三维地质曲面构建方法", 《深圳大学学报理工版》 *
张敬博: "层状岩体三维建模与可视化方法研究", 《中国优秀硕士学位论文全文数据库 信息科技辑》 *
王孜越: "基于三维实景平台的水电工程野外地质编录填图系统研究", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》 *

Also Published As

Publication number Publication date
CN106780741B (zh) 2020-05-12

Similar Documents

Publication Publication Date Title
CN102609982B (zh) 空间地质数据非结构化模式的拓扑发现方法
CN103824329B (zh) 一种地质勘探三维可视化储量估算方法
CN102194253A (zh) 一种面向三维地质层面结构的四面体网格生成方法
CN103500468B (zh) 基于边界约束的空间曲面重建方法
CN101303414A (zh) 一种基于水平集的地层面及地质体生成方法
CN105913494A (zh) 多尺度裂缝精细地质建模及数值模拟方法和装置
CN105118091B (zh) 一种构建多精度非均匀地质网格曲面模型的方法和系统
CN104240300B (zh) 基于分布式并行的大规模点云复杂空间曲面重构方法
CN104167020B (zh) 基于约束Delaunay三角网的空间倒转曲面重建方法
CN102867332B (zh) 基于复杂边界约束的多级细分网格曲面拟合方法
CN103907118A (zh) 用于在储层模拟系统中进行粗化的系统和方法
CN110334365B (zh) 一种非均质压裂后储层流动数值模拟方法及系统
CN103198192A (zh) 一种用于虚拟三维零部件快速装配的吸附方法
CN109829217A (zh) 压裂性裂缝油藏产能模拟方法及装置
CN102298660A (zh) 一种离散元法边界建模的通用方法
CN102147934A (zh) 基于tin的地质界面三维形态分析方法
Valero-Lara et al. A non-uniform Staggered Cartesian grid approach for Lattice-Boltzmann method
CN112687001B (zh) 三维地质结构模型随机生成及不确定性分析方法
CN105653746B (zh) 基于嵌入式离散裂缝模型的压裂井建模及模拟方法
CN106780741A (zh) 基于多点产状数据的地质曲面快速拟合插值方法
CN108595782A (zh) 一种离散裂缝中基质与裂缝间的传质计算方法
CN105957141A (zh) 一种基于对称结构的三维花朵建模方法
CN104881548B (zh) 基于离散域力场的相邻工序模型几何演变的获取方法
Wu et al. Characterization and reconstruction of rough fractures based on vector statistics
CN104036024A (zh) 一种基于GACUC和Delaunay三角网的空间聚类方法

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
GR01 Patent grant
GR01 Patent grant