CN106780741A - 基于多点产状数据的地质曲面快速拟合插值方法 - Google Patents
基于多点产状数据的地质曲面快速拟合插值方法 Download PDFInfo
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 32
- 238000004458 analytical method Methods 0.000 claims abstract description 8
- 238000006073 displacement reaction Methods 0.000 claims description 4
- 239000002245 particle Substances 0.000 claims 8
- 241001269238 Data Species 0.000 claims 1
- 238000004088 simulation Methods 0.000 abstract description 3
- 238000005457 optimization Methods 0.000 description 4
- 238000004364 calculation method Methods 0.000 description 2
- 238000013178 mathematical model Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
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中的基于平面弹性理论的曲面样条插值法,是一类用于构造地质曲面的插值数学模型;该方法将计算域视为无限延伸平面,根据实际情况确定插值基平面,在基平面上一点给予一定的垂直位移时,基平面各处均会随着产生位移。
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)
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 |
-
2017
- 2017-02-10 CN CN201710073770.2A patent/CN106780741B/zh active Active
Patent Citations (5)
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)
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 |