CN106875484B - 一种基于三维地形的地质堆积体快速拟合建模方法 - Google Patents

一种基于三维地形的地质堆积体快速拟合建模方法 Download PDF

Info

Publication number
CN106875484B
CN106875484B CN201710073627.3A CN201710073627A CN106875484B CN 106875484 B CN106875484 B CN 106875484B CN 201710073627 A CN201710073627 A CN 201710073627A CN 106875484 B CN106875484 B CN 106875484B
Authority
CN
China
Prior art keywords
geological
depth
accumulation body
node
dimensional terrain
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
Application number
CN201710073627.3A
Other languages
English (en)
Other versions
CN106875484A (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 CN201710073627.3A priority Critical patent/CN106875484B/zh
Publication of CN106875484A publication Critical patent/CN106875484A/zh
Application granted granted Critical
Publication of CN106875484B publication Critical patent/CN106875484B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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
    • 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/20—Finite element generation, e.g. wire-frame surface description, tesselation

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Geometry (AREA)
  • Software Systems (AREA)
  • Computer Graphics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Remote Sensing (AREA)
  • Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

本发明涉及一种根据地质勘测数据进行计算机地质建模、模拟和分析的方法，其公开了一种基于三维地形的地质堆积体快速拟合建模方法，用以根据野外地质勘测获得的地质堆积体信息进行计算机处理、重构堆积体模型。该方法包括以下步骤：A、在三维地形实景中标记地质堆积体的地表出露范围；B、指定深度控制点，给予其距离地表的预测深度；C、在多段线内部构建插值点，形成Delaunay划分网格面；D、针对网格面中的每一个节点，计算其深度；E、将网格曲面展示到三维地形实景之中，形成地质堆积体的下包络曲面，并获得堆积体实体。该方法能够在野外现场快速实时推求地质堆积体的分布、形状、分界面等地质信息，为工程地质现场勘探提供重要参考信息与快速建模工具。

Description

一种基于三维地形的地质堆积体快速拟合建模方法
技术领域
本发明涉及一种根据地质勘测数据进行计算机地质建模、模拟和分析的方法，特别涉及一种基于三维地形的地质堆积体快速拟合建模方法。
背景技术
工程地质野外勘探过程中，地表堆积体的绘制和记录是地质填图的重要工作。准确记录并分析地质堆积体，能够对区域地质的分布情况做出完整的描述，并对地质灾害的预测。传统的地质勘测手段仅能通过地质堆积体的地表出露边界估算堆积体范围，通过少许采样点预测堆积体地下分布情况。预测结果的准确性很大程度上取决于地质工作者的经验，所获取的地质信息和构建的地质模型很难被用于工程地质三维综合分析之中。
Delaunay三角化平面[1,2]是均分曲面的主要方法；而非均匀有理B样条(NURBS)和不规则三角网(TIN)则是曲面建模的主要方法[3,4]，因此，基于以上情况和研究，本发明通过建立地面三角网，提出了一种基于三维地形的地质堆积体快速拟合建模方法。
参考文献：
[1]鲁才，李林超，胡光岷.基于矢量闭合的三维建模方法[P].发明专利ZL201410449265.X，2014
[2]王颀，赵洲，孙健等。一种基于大规模地震数据的三维地质曲面重建方法[P].发明专利ZL201610397592.4，2016
[3]钟登华，李明超，刘杰等。水利水电工程地质信息的三维统一模型构建方法[P].发明专利ZL200610013425.1，2006
[4]李明超，钟登华，秦朝霞，等.基于三维地质模型的工程岩体结构精细数值建模[J]，岩石力学与工程学报,2007,26(9):1893-1898
发明内容
本发明所要解决的技术问题是：提出一种基于三维地形的地质堆积体快速拟合建模方法，用以根据野外地质勘测获得的地质堆积体信息进行计算机处理、重构堆积体模型。
本发明解决上述技术问题所采用的方案是：
一种基于三维地形的地质堆积体快速拟合建模方法，包括以下步骤：
步骤A、根据野外地质现场测绘，在三维地形实景中标记地质堆积体的地表出露范围，以封闭多段线表示；
步骤B、指定封闭多段线内部的深度控制点，给予其距离地表的预测深度；
步骤C、在多段线内部构建插值点，形成Delaunay划分网格面；
步骤D、针对网格面中的每一个节点，根据以下规则计算其深度：堆积体边界处的深度为0，内部指定点处的深度为给定值；其余节点根据其与边界和每个指定点处的距离插值计算出对应深度；插值完成后的网格面从边界到内部指定点处渐进加深，形成平滑下凹曲面；
步骤E、将网格曲面展示到三维地形实景之中，形成地质堆积体的下包络曲面，并获得堆积体实体。
作为进一步优化，步骤A中的三维地形实景，是工程地质野外编录区域所对应的带实景贴图的三维地形模型，其是由陆摄或航摄手段获取的影像和高程数据叠加处理形成的能够准确表征野外地形环境的三维数字场景。
作为进一步优化，步骤A中的地质堆积体的地表出露范围，是地质工作者根据现场勘探信息，在三维地形实景上圈定的地质堆积体边界，该边界是首尾闭合的封闭多段线，多段线内属于堆积体。
作为进一步优化，步骤A中的野外地质现场测绘，是地质工作者在野外现场，根据现场实际地质情况，利用平板电脑在三维地形实景上快速勾划出地质堆积体边界的过程；三维地形实景和地质堆积体边界是快速拟合建模方法的重要数据来源。
作为进一步优化，步骤B中的预测深度，是地质工作者根据现场地勘现象和地质经验所预测的堆积体最大深度，在垂直方向上，该深度以上部分都属于地质堆积体。
作为进一步优化，步骤D具体包括：
步骤D1：获取堆积体边界的封闭多段线，并定义其上的点深度为0；
步骤D2：获取堆积体内部的若干深度控制点，并定义其深度为用户输入值；
步骤D3：遍历网格面中的每一个节点，计算每一个节点在每一个深度控制点下的深度值，包括：寻找距离该节点最近的边界；作深度控制点通过该节点至该边界的线段；以该节点距离控制点距离和上述线段距离的比值作为权重参数，根据控制点深度计算获得该节点处的深度；
步骤D4：将每一个节点在每一个深度控制点下的深度值按距离每个控制点的距离设定权重，综合所有控制点的深度值，计算获得该节点的综合深度；
步骤D5：针对所有节点重复D3-D4的步骤，计算每个节点的综合深度，并生成能反映各部位不同深度值的平滑下凹网格曲面。
作为进一步优化，步骤E中的堆积体实体是堆积体的下包络曲面与三维实景地形面交切形成的三维实体，用于表征地质堆积体在三维实景中的位置与形态。
本发明的有益效果是：
在数字三维地形实景的支持下，该方法能够在野外现场快速实时推求地质堆积体的分布、形状、分界面等地质信息，为工程地质现场勘探提供重要参考信息与快速建模工具。
附图说明
图1为基于三维地形的地质堆积体快速拟合建模方法流程图；
图2为插值拟合示意图。
具体实施方式
本发明旨在提出一种基于三维地形的地质堆积体快速拟合建模方法，用以根据野外地质勘测获得的地质堆积体信息进行计算机处理、重构堆积体模型。
参照图1，本发明中的基于三维地形的地质堆积体快速拟合建模方法包括以下步骤：
步骤A，根据野外地质现场测绘，在三维地形实景中用直线标记地质堆积体的地表出露范围，直线首位相连，形成多段线；标注完毕后，检验多段线是否封闭，若不封闭则重新选择；如果封闭则可进行下一步；
步骤B，指定封闭多段线内部的若干深度控制点，给予其距离地表的预测深度；若不选择自定义多点深度，则需赋值中心点深度；
步骤C，在多段线内部构建插值点，形成Delaunay划分网格面；
步骤D，针对网格面中的每一个节点，根据以下规则计算其深度：堆积体边界处的深度为0，内部指定点处的深度为给定值；其余节点根据其与边界和每个指定点处的距离插值计算出对应深度；插值完成后的网格面从边界到内部指定点处渐进加深，形成平滑下凹曲面；
步骤E，将网格曲面展示到三维地形实景之中，形成地质堆积体的下包络曲面，并获得堆积体实体。
本发明一个实施例如下：
(1)获取初始条件。在三维实景中选择要研究的区域，以封闭的多段线表示，多段线与地表吻合；系统自动判别多段线是否封闭，如不封闭无法进行下一步计算；系统可以自动捕捉研究区域的地质信息。
(2)建立直线方程，以图2为例，A、B、M、X横纵坐标已知，利用MX线上两点确定MX线方程：
Ax1+By1+C＝0
利用AB线上两点确定AB线方程：
Ax2+By2+C＝0
(3)求辅助交点。联立两直线方程：
Figure BDA0001223577330000041
得到直线相交交点O(x0,y0)。
(4)计算插值权重，获取插值深度。分别计算MO距离Dm，MX距离D1，设置中心点深度为最大深度，并计算其他点的权值，最终用地表高程和权值深度做差得到地表点的下凹深度：
Figure BDA0001223577330000042
其中，z0为区域内各点下凹面处高程；z为区域内各点地表高程；D1为各点到中心点实际距离；Dm为各点到中心点最大距离；Hmax为中心点最大深度。
(5)建立插值曲面。计算得到Delaunay划分网格面，可将计算得到的模型从工程地质野外数字测绘编录系统中导出。

Claims (5)

1.一种基于三维地形的地质堆积体快速拟合建模方法，其特征在于，包括以下步骤：
步骤A、根据野外地质现场测绘，在三维地形实景中标记地质堆积体的地表出露范围，以封闭多段线表示；
步骤B、指定封闭多段线内部的深度控制点，给予其距离地表的预测深度；所述预测深度，是地质工作者根据现场地勘现象和地质经验所预测的堆积体最大深度，在垂直方向上，该深度以上部分都属于地质堆积体；
步骤C、在多段线内部构建插值点，形成Delaunay划分网格面；
步骤D、针对网格面中的每一个节点，根据以下规则计算其深度：堆积体边界处的深度为0，内部指定点处的深度为给定值；其余节点根据其与边界和每个指定点处的距离插值计算出对应深度；插值完成后的网格面从边界到内部指定点处渐进加深，形成平滑下凹曲面，具体包括步骤D1-D5：
步骤D1：获取堆积体边界的封闭多段线，并定义其上的点深度为0；
步骤D2：获取堆积体内部的若干深度控制点，并定义其深度为用户输入值；
步骤D3：遍历网格面中的每一个节点，计算每一个节点在每一个深度控制点下的深度值，包括：寻找距离该节点最近的边界；作深度控制点通过该节点至该边界的线段；以该节点距离控制点距离和上述线段距离的比值作为权重参数，根据控制点深度计算获得该节点处的深度；
步骤D4：将每一个节点在每一个深度控制点下的深度值按距离每个控制点的距离设定权重，综合所有控制点的深度值，计算获得该节点的综合深度；
步骤D5：针对所有节点重复D3-D4的步骤，计算每个节点的综合深度，并生成能反映各部位不同深度值的平滑下凹网格曲面；
步骤E、将网格曲面展示到三维地形实景之中，形成地质堆积体的下包络曲面，并获得堆积体实体。
2.如权利要求1所述的一种基于三维地形的地质堆积体快速拟合建模方法，其特征在于，步骤A中的三维地形实景，是工程地质野外编录区域所对应的带实景贴图的三维地形模型，其是由陆摄或航摄手段获取的影像和高程数据叠加处理形成的能够准确表征野外地形环境的三维数字场景。
3.如权利要求1所述的一种基于三维地形的地质堆积体快速拟合建模方法，其特征在于，步骤A中的地质堆积体的地表出露范围，是地质工作者根据现场勘探信息，在三维地形实景上圈定的地质堆积体边界，该边界是首尾闭合的封闭多段线，多段线内属于堆积体。
4.如权利要求1所述的一种基于三维地形的地质堆积体快速拟合建模方法，其特征在于，步骤A中的野外地质现场测绘，是地质工作者在野外现场，根据现场实际地质情况，利用平板电脑在三维地形实景上快速勾划出地质堆积体边界的过程；三维地形实景和地质堆积体边界是快速拟合建模方法的重要数据来源。
5.如权利要求1所述的一种基于三维地形的地质堆积体快速拟合建模方法，其特征在于，步骤E中的堆积体实体是堆积体的下包络曲面与三维实景地形面交切形成的三维实体，用于表征地质堆积体在三维实景中的位置与形态。
CN201710073627.3A 2017-02-10 2017-02-10 一种基于三维地形的地质堆积体快速拟合建模方法 Active CN106875484B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710073627.3A CN106875484B (zh) 2017-02-10 2017-02-10 一种基于三维地形的地质堆积体快速拟合建模方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710073627.3A CN106875484B (zh) 2017-02-10 2017-02-10 一种基于三维地形的地质堆积体快速拟合建模方法

Publications (2)

Publication Number Publication Date
CN106875484A CN106875484A (zh) 2017-06-20
CN106875484B true CN106875484B (zh) 2020-03-20

Family

ID=59167007

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710073627.3A Active CN106875484B (zh) 2017-02-10 2017-02-10 一种基于三维地形的地质堆积体快速拟合建模方法

Country Status (1)

Country Link
CN (1) CN106875484B (zh)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107860375A (zh) * 2017-11-03 2018-03-30 广州地理研究所 一种基于三维激光扫描技术的滑坡灾害体积快速提取方法
CN111127645B (zh) * 2019-12-25 2023-02-10 中国电建集团中南勘测设计研究院有限公司 三维地质覆盖层界面建模方法、装置、设备及存储介质
CN115165014A (zh) * 2022-07-08 2022-10-11 中国地质科学院岩溶地质研究所 一种计算岩溶洼地尾渣库中尾渣堆体积的方法
CN117876621A (zh) * 2024-03-07 2024-04-12 贵州省第一测绘院(贵州省北斗导航位置服务中心) 一种基于高分辨率遥感图像和地形数据的国土测绘方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101582173A (zh) * 2009-06-24 2009-11-18 中国石油集团川庆钻探工程有限公司 复杂地质构造块状模型构建方法
CN105184864A (zh) * 2015-08-14 2015-12-23 中国能源建设集团安徽省电力设计院有限公司 一种用于天然地基换填工程量计算的场地地层三维地质结构模型生成方法
CN106097446A (zh) * 2016-06-07 2016-11-09 中国石油集团川庆钻探工程有限公司地球物理勘探公司 一种基于大规模地震数据的三维地质曲面重建方法

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB2512372B (en) * 2013-03-28 2020-07-29 Total Sa Method of modelling a subsurface volume
FR3016983A1 (fr) * 2014-01-28 2015-07-31 IFP Energies Nouvelles Procede de construction d'un maillage volumique pour la modelisation de structures geologiques

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101582173A (zh) * 2009-06-24 2009-11-18 中国石油集团川庆钻探工程有限公司 复杂地质构造块状模型构建方法
CN105184864A (zh) * 2015-08-14 2015-12-23 中国能源建设集团安徽省电力设计院有限公司 一种用于天然地基换填工程量计算的场地地层三维地质结构模型生成方法
CN106097446A (zh) * 2016-06-07 2016-11-09 中国石油集团川庆钻探工程有限公司地球物理勘探公司 一种基于大规模地震数据的三维地质曲面重建方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
An effective method for 3D geological modeling with multi-source data integration;Qiang Wu 等;《Computers & Geosciences》;20050228;第31卷(第1期);全文 *
基于网格的三维地质体建模方法研究;王威 等;《岩土力学》;20120430;第33卷(第4期);全文 *
自然邻点插值方法在三维地质建模中的应用;郭艳军等;《解放军理工大学学报(自然科学版)》;20091231;第10卷(第6期);全文 *

Also Published As

Publication number Publication date
CN106875484A (zh) 2017-06-20

Similar Documents

Publication Publication Date Title
CN112381937B (zh) 一种基于钻孔和复杂地质剖面的多源地质数据耦合建模方法
Enge et al. From outcrop to reservoir simulation model: Workflow and procedures
CN106875484B (zh) 一种基于三维地形的地质堆积体快速拟合建模方法
CN106846476B (zh) 基于三维实景和赤平投影的岩石块体稳定性快速评价方法
CN111079217B (zh) 一种基于bim的岩土工程综合勘察信息解译方法及系统
CN103824329B (zh) 一种地质勘探三维可视化储量估算方法
CN116152461B (zh) 地质建模方法、装置、计算机设备及计算机可读存储介质
Zhong et al. Enhanced NURBS modeling and visualization for large 3D geoengineering applications: an example from the Jinping first-level hydropower engineering project, China
CN112862967B (zh) 用于建立三维地质模型的方法和装置
CN104504754B (zh) 一种油气储层多点统计建模的方法及装置
CN111161415B (zh) 一种岩溶区地下水岩溶发育通道三维空间模拟方法
FR3043227A1 (zh)
CN111383336A (zh) 三维地质模型构建方法
Cosso et al. Surveying and mapping a cave using 3d laser scanner: the open challenge with free and open source software
CN116486025A (zh) 一种基于大数据云计算技术的城市地质数据处理平台
CN112231818B (zh) 一种基于bim的山地城市三维地质建模方法
CN115690354A (zh) 一种基于三维实景数值分析的浅埋隧道施工动态控制方法
CN110765665B (zh) 一种地学动态建模方法及系统
CN115100373A (zh) 地表与地下环境三维一体化表征方法、装置、介质及设备
CN109949415B (zh) 一种三维地表与地质体模型拓扑一致建模的系统及方法
Wu et al. The 3D modelling techniques of digital geological mapping
Ji et al. An automated method to build 3D multi-scale geological models for engineering sedimentary layers with stratum lenses
CN111951393B (zh) 水平构造台地的三维模型构造方法及装置
Shen et al. Three-dimensional modeling of loose layers based on stratum development law
CN113887046B (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
GR01 Patent grant
GR01 Patent grant