CN111951396B - 知识驱动的基岩面三维模型构建方法及装置 - Google Patents

知识驱动的基岩面三维模型构建方法及装置 Download PDF

Info

Publication number
CN111951396B
CN111951396B CN202010794872.5A CN202010794872A CN111951396B CN 111951396 B CN111951396 B CN 111951396B CN 202010794872 A CN202010794872 A CN 202010794872A CN 111951396 B CN111951396 B CN 111951396B
Authority
CN
China
Prior art keywords
bedrock
point
obtaining
drilling
value
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
CN202010794872.5A
Other languages
English (en)
Other versions
CN111951396A (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.)
Nanjing Normal University
Original Assignee
Nanjing Normal University
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 Nanjing Normal University filed Critical Nanjing Normal University
Priority to CN202010794872.5A priority Critical patent/CN111951396B/zh
Publication of CN111951396A publication Critical patent/CN111951396A/zh
Application granted granted Critical
Publication of CN111951396B publication Critical patent/CN111951396B/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
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T3/00—Geometric image transformations in the plane of the image
    • G06T3/40—Scaling of whole images or parts thereof, e.g. expanding or contracting
    • G06T3/4007—Scaling of whole images or parts thereof, e.g. expanding or contracting based on interpolation, e.g. bilinear interpolation

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Geometry (AREA)
  • Software Systems (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Remote Sensing (AREA)
  • Computer Graphics (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

本发明公开了一种知识驱动的基岩面三维模型构建方法及装置，方法具体包括步骤：首先，基于研究区地质图，获取基岩出露边界线；其次，基于基岩出露边界线，裁剪获得基岩埋藏区域Dem和出露区Dem；然后，基于钻孔点数据，根据地貌类型，构建虚拟钻孔和基岩埋深等值线；最后，基于基岩埋深等值线和钻孔点数据，插值得到各地貌单元的基岩面三维模型，与基岩出露区域的Dem合并后，得到研究区的完整基岩面三维模型。本发明实现了一种在钻孔数据相对稀疏的情况下，基于地学知识规律相对精确地构建基岩面三维模型的方法。本发明与现有技术相比，在一定程度上提高了钻孔稀疏区域的基岩面建模精度。

Description

知识驱动的基岩面三维模型构建方法及装置
技术领域
本发明涉及三维地质建模技术领域，尤其涉及一种知识驱动的基岩面三维模型构建方法及装置。
背景技术
基岩面是指地球陆地表面疏松物质(土壤和底土)底下的坚硬岩层面，它将第四系地质体和基岩地质体分隔开来，约束着第四系松散层的空间形态和分布范围。基岩面建模对于第四系松散层三维建模、地下资源评价、地下空间规划和利用等具有重要的研究意义和实用价值。
早期，主要根据地球物理方法观测和分析各种地球物理场的变化来探测地层岩性、地质构造等，并通过数据反演构建基岩面模型。该方法虽然适用范围广，但精度较低，且解析结果及建模质量因人而异。后来，随着工程地质钻孔数据量的不断积累，逐步形成了基于钻孔的基岩面建模方法。在钻孔数据丰富，且分布相对均匀的情况下，该方法建模精度较高；然而，该方法过分依赖钻孔数据，导致钻孔稀疏区域的基岩面模型精度较低。
发明内容
发明目的：本发明针对现有技术存在的问题，提供一种在提高钻孔稀疏区域的建模精度的知识驱动的基岩面三维模型构建方法及装置。
技术方案：本发明所述的知识驱动的基岩面三维模型构建方法包括：
(1)基于研究区的地质图，获取基岩出露区边界线，并对研究区的DEM数据进行裁剪，获得基岩埋藏区域的DEM数据dem1和基岩出露区的DEM数据dem2；
(2)根据研究区的地貌单元图层，将基岩埋藏区域划分为岗地、河谷、湖泊和平原这四类地貌单元类型，并存入地貌单元集合GU中；
(3)从地貌单元集合GU中获取任一地貌单元gui，若gui为平原，则执行步骤(6)；否则，执行步骤(4)；
(4)获取gui的钻孔集合RD，并根据钻孔集合RD和gui的地貌单元类型计算gui待构建埋深等值线的最大埋深值Hmax和最小埋深值Hmin，之后执行步骤(5)；
(5)基于Hmax和Hmin，构建地貌单元gui的埋深等值线，执行步骤(7)；
(6)基于Dem1，构建地貌单元gui的虚拟钻孔集合VD；
(7)基于埋深等值线、虚拟钻孔集合VD、钻孔集合RD和Dem1，构建地貌单元gui的基岩面三维模型bmi；
(8)重复执行步骤(3)至步骤(7)，直至所有地貌单元的基岩面三维模型构建完成，并存储至基岩面模型集合BM中；
(9)将Dem2和基岩面模型集合BM中所有栅格合并，生成研究区域完整的基岩面三维模型。
进一步的，步骤(4)具体包括：
(4-1)获取gui内的所有钻孔数据到集合RD＝{rdj|j＝1,2,…,nrd}，其中，rdj表示gui内第j个钻孔，nrd为钻孔个数；
(4-2)获取埋深值最大的钻孔rdm，将其埋深值记为Hm，根据下式计算最大埋深值Hmax；
式中，d为预设埋深线间距，符号[]为取整符号；
(4-3)若gui为岗地，则直接将Hn赋值为零；否则，获取埋深值最小的钻孔rdn，将其埋深值记为Hn；
(4-4)基于Hn采用下式，计算最小埋深值Hmin；
(4-5)判断是否满足Hmin＝d，若是，则执行步骤(5)；否则执行步骤(4-6)；
(4-6)根据钻孔rdm和钻孔rdn分别到地貌单元边界的最短距离调整最小值Hmin，并执行步骤(5)。
进一步的，步骤(4-6)具体包括：
(4-6-1)分别获取钻孔rdm、钻孔rdn到地貌单元边界的最短距离Dism、Disn；
(4-6-2)根据下式计算基岩埋深为Hmin-d的等值线到地貌单元边界的距离Dis；
式中，Hn、Hm分别为钻孔点rdn、rdm的基岩埋深值,符号||为取绝对值符号；
(4-6-3)若Dis≥0，则令Hmin＝Hmin-d，然后执行步骤(4-5)；否则，执行步骤(5)。
进一步的，步骤(5)具体包括：
(5-1)基于地貌单元gui的钻孔集合RD，获取用于构建基岩埋深为Hmin的等值线ch1的钻孔rdu和rdv；
(5-2)若gui为岗地，则基于地质图，将gui的基岩出露线作为埋深基准线ch0，否则，将gui地貌单元边界线作为埋深基准线ch0；
(5-3)分别计算钻孔点rdu、rdv到基岩埋深基准线ch0的最短距离Disu、Disv；
(5-4)根据下式，计算基岩埋深为Hmin的等值线ch1到基岩埋深基准线ch0的距离Dis1-0；
式中，Hu、Hv分别为钻孔点rdu、rdv的基岩埋深值,符号||为取绝对值符号；
(5-5)基于基岩埋深基准线ch0和距离Dis1-0，执行缓冲区分析，获取缓冲区边界线BLine0；
(5-6)根据地貌单元gui的范围裁剪BLine0，生成基岩埋深为Hmin的等值线ch1；
(5-7)循环执行(5-1)至(5-6)，依次构建基岩埋深为Hmin+d，Hmin+2d，…，Hmax的基岩埋深等值线，直至地貌单元gui的所有基岩埋深等值线构建完成。
进一步的，步骤(5-1)具体包括：
(5-1-1)获取gui的钻孔集合RD，并根据埋深值，将钻孔集合RD分为子集合RD1和RD2，其中子集合RD1中钻孔的埋深值均不大于Hmin，子集合RD2中钻孔的埋深值均大于Hmin；
(5-1-2)若RD1和RD2均不为空，则获取RD1中埋深值最大的钻孔，记为rdu，其基岩埋深记为Hu，获取RD2中埋深值最小的钻孔，记为rdv，其基岩埋深记为Hv，然后执行步骤(5-2)；否则，执行步骤(5-1-3)；
(5-1-3)若RD1为空集，获取RD2中埋深值最小的两个钻孔作为rdu、rdv；若RD2为空集，则获取RD1中埋深值最大的两个钻孔作为rdu、rdv。
进一步的，步骤(6)具体包括：
(6-1)获取地貌单元gui的钻孔集合RD，以钻孔集合RD为基础数据通过反距离加权插值方法得到栅格面DDem；
(6-2)从基岩埋藏区域的DEM数据dem1中获取地貌单元gui的DEM数据DEMi；
(6-3)构建地貌单元gui的外包矩形Rec，并在矩形Rec内部生成10×10网格，将位于gui的网格交叉点存储到点集合PR中；
(6-4)从点集合PR中获取任一点prp，获取点prp在DDem上的埋深值dr，以及点prp在DEMi上的高程值hr，计算dr和hr的差值difr，并作为点prp的Z坐标值；重复执行本步骤，直至完成集合PR中所有点的Z坐标值计算；
(6-5)读取地貌单元gui的边界线到集合BL＝{blk|k＝1,2,…,nbl}，其中，blk表示第k条边界线，nbl为边界线个数；
(6-6)从集合BL中获取任一边界线blk，将边界线blk离散为点集合PLk＝{plk,l|l＝1,2,…,npk}，其中，plk,l表示边界线blk上第l个点，npk为边界线blk上的点个数；若blk为基岩出露边界线，则执行步骤(6-7)，否则执行步骤(6-8)；
(6-7)从点集合PLk中获取任一点plk,l，获取点plk,l在DEMi上的高程值hl，作为点plk,l的Z坐标值；重复执行该步骤直至完成集合PLk中所有点的Z坐标值计算，然后执行步骤(6-9)；
(6-8)从点集合PLk中获取任一点plk,l，获取点plk,l在DDem上的埋深值dl，以及点plk,l在DEMi上的高程值hl，计算dl和hl的差值difl，作为点plk,l的Z坐标值；重复执行该步骤直至完成集合PLk中所有点的Z坐标值计算，然后执行步骤(6-9)；
(6-9)循环执行步骤(6-6)到(6-8)，直至集合BL中所有边界线被遍历，得到地貌单元gui所有边界线的点的Z坐标值；
(6-10)将包含有Z坐标值的集合PR中的点和所有集合PLk中的点存储到集合VD中，作为地貌单元gui的虚拟钻孔集合。
进一步的，步骤(7)具体包括：
(7-1)若地貌单元gui为平原，执行步骤(7-5)；否则将基岩埋深等值线上的点存储到集合PC中；
(7-2)从集合PC中获取任一点pcb，获取点pcb在地貌单元gui的DEM数据DEMi上的地表高程值hp，并基于基岩埋深等值线，获取点pcb的埋深值dp，计算高程值hp与埋深值dp的差值difp，作为点pcb的Z坐标值；重复执行该步骤直至完成集合PC中所有点的Z坐标值计算；
(7-3)从钻孔集合RD中获取任一钻孔rdj，读取rdj的高程属性得到高程值hrd，以及埋深属性得到埋深值drd；
(7-4)获取点rdj的X坐标值和Y坐标值，生成点P，计算高程值hrd与埋深值drd的差值difrd，作为点P的Z坐标值，将点P存储到集合PC中；
(7-5)循环执行步骤(7-3)至步骤(7-4)，直至集合RD中所有点数据都被处理；
(7-6)将地貌单元gui的虚拟钻孔集合VD中所有的点都存储到集合PC中；
(7-7)以集合PC为基础数据通过反距离加权插值方法得到栅格面bmi，栅格面bmi即为地貌单元gui的基岩面三维模型。
有益效果：本发明与现有技术相比，其显著优点是：本发明从不同地貌类型的基岩面埋深及基岩面起伏规律出发，形成一种知识驱动的基岩面三维建模方法，提高了钻孔稀疏区域的基岩面三维建模精度。
附图说明
图1是本实施例中采用的实验数据((a)钻孔分布图，(b)DEM，(c)基岩出露图层，(d)地貌单元图层)；
图2本发明的流程示意图；
图3是本发明岗地单元构建基岩埋深等值线的子流程图；
图4是本发明构建虚拟钻孔点的子流程图
图5是本实施例中构建的基岩面模型；
图6是本实施例中构建的基岩面模型三维示意图。
具体实施方式
下面对本发明技术方案进行详细说明，本实例的实验区采用的是南京市(图1)。区内地貌类型也较为齐全，钻孔数据虽相对丰富，但分布不均。下面结合附图，并通过描述一个具体的实施例，来进一步说明。
如图2所示，本实施例提供的知识驱动的基岩面三维模型构建方法包括：
(1)基于研究区的地质图，获取基岩出露区边界线，并对研究区的DEM数据进行裁剪，获得基岩埋藏区域的DEM数据dem1和基岩出露区的DEM数据dem2。
(2)根据研究区的地貌单元图层，将基岩埋藏区域划分为岗地、河谷、湖泊和平原这四类地貌单元类型，并存入地貌单元集合GU中。
(3)从地貌单元集合GU中获取任一地貌单元gui，若gui为平原，则执行步骤(6)；否则，执行步骤(4)。
(4)获取gui的钻孔集合RD，并根据钻孔集合RD和gui的地貌单元类型计算gui待构建埋深等值线的最大埋深值Hmax和最小埋深值Hmin，之后执行步骤(5)。
该步骤具体包括：
(4-1)获取gui内的所有钻孔数据到集合RD＝{rdj|j＝1,2,…,nrd}，其中，rdj表示gui内第j个钻孔，nrd为钻孔个数；
(4-2)获取埋深值最大的钻孔rdm，将其埋深值记为Hm，根据下式计算最大埋深值Hmax；
式中，d为预设埋深线间距，本实施例中，设置d＝5，符号[]为取整符号；
(4-3)若gui为岗地，则直接将Hn赋值为零；否则，获取埋深值最小的钻孔rdn，将其埋深值记为Hn；
(4-4)基于Hn采用下式，计算最小埋深值Hmin；
(4-5)判断是否满足Hmin＝d，若是，则执行步骤(5)；否则执行步骤(4-6)；
(4-6)根据钻孔rdm和钻孔rdn分别到地貌单元边界的最短距离调整最小值Hmin，并执行步骤(5)。
其中，步骤(4-6)具体包括：
(4-6-1)分别获取钻孔rdm、钻孔rdn到地貌单元边界的最短距离Dism、Disn；
(4-6-2)根据下式计算基岩埋深为Hmin-d的等值线到地貌单元边界的距离Dis；
式中，Hn、Hm分别为钻孔点rdn、rdm的基岩埋深值,符号||为取绝对值符号；
(4-6-3)若Dis≥0，则令Hmin＝Hmin-d，然后执行步骤(4-5)；否则，执行步骤(5)。
(5)基于Hmax和Hmin，构建地貌单元gui的埋深等值线，执行步骤(7)。
如图3所示，该步骤具体包括：
(5-1)基于地貌单元gui的钻孔集合RD，获取用于构建基岩埋深为Hmin的等值线ch1的钻孔rdu和rdv；
其中，步骤(5-1)具体包括：
(5-1-1)获取gui的钻孔集合RD，并根据埋深值，将钻孔集合RD分为子集合RD1和RD2，其中子集合RD1中钻孔的埋深值均不大于Hmin，子集合RD2中钻孔的埋深值均大于Hmin；
(5-1-2)若RD1和RD2均不为空，则获取RD1中埋深值最大的钻孔，记为rdu，其基岩埋深记为Hu，获取RD2中埋深值最小的钻孔，记为rdv，其基岩埋深记为Hv，然后执行步骤(5-2)；否则，执行步骤(5-1-3)；
(5-1-3)若RD1为空集，获取RD2中埋深值最小的两个钻孔作为rdu、rdv；若RD2为空集，则获取RD1中埋深值最大的两个钻孔作为rdu、rdv。
(5-2)若gui为岗地，则基于地质图，将gui的基岩出露线作为埋深基准线ch0，否则，将gui地貌单元边界线作为埋深基准线ch0；
(5-3)分别计算钻孔点rdu、rdv到基岩埋深基准线ch0的最短距离Disu、Disv；
(5-4)根据下式，计算基岩埋深为Hmin的等值线ch1到基岩埋深基准线ch0的距离Dis1-0；
式中，Hu、Hv分别为钻孔点rdu、rdv的基岩埋深值,符号||为取绝对值符号；
(5-5)基于基岩埋深基准线ch0和距离Dis1-0，执行缓冲区分析，获取缓冲区边界线BLine0；
(5-6)根据地貌单元gui的范围裁剪BLine0，生成基岩埋深为Hmin的等值线ch1；
(5-7)循环执行(5-1)至(5-6)，依次构建基岩埋深为Hmin+d，Hmin+2d，…，Hmax的基岩埋深等值线，直至地貌单元gui的所有基岩埋深等值线构建完成。
(6)基于Dem1，构建地貌单元gui的虚拟钻孔集合VD。
如图4所示，该步骤具体包括：
(6-1)获取地貌单元gui的钻孔集合RD，以钻孔集合RD为基础数据通过反距离加权插值方法得到栅格面DDem；
(6-2)从基岩埋藏区域的DEM数据dem1中获取地貌单元gui的DEM数据DEMi；
(6-3)构建地貌单元gui的外包矩形Rec，并在矩形Rec内部生成10×10网格，将位于gui的网格交叉点存储到点集合PR中；
(6-4)从点集合PR中获取任一点prp，获取点prp在DDem上的埋深值dr，以及点prp在DEMi上的高程值hr，计算dr和hr的差值difr，并作为点prp的Z坐标值；重复执行本步骤，直至完成集合PR中所有点的Z坐标值计算；
(6-5)读取地貌单元gui的边界线到集合BL＝{blk|k＝1,2,…,nbl}，其中，blk表示第k条边界线，nbl为边界线个数；
(6-6)从集合BL中获取任一边界线blk，将边界线blk离散为点集合PLk＝{plk,l|l＝1,2,…,npk}，其中，plk,l表示边界线blk上第l个点，npk为边界线blk上的点个数；若blk为基岩出露边界线，则执行步骤(6-7)，否则执行步骤(6-8)；
(6-7)从点集合PLk中获取任一点plk,l，获取点plk,l在DEMi上的高程值hl，作为点plk,l的Z坐标值；重复执行该步骤直至完成集合PLk中所有点的Z坐标值计算，然后执行步骤(6-9)；
(6-8)从点集合PLk中获取任一点plk,l，获取点plk,l在DDem上的埋深值dl，以及点plk,l在DEMi上的高程值hl，计算dl和hl的差值difl，作为点plk,l的Z坐标值；重复执行该步骤直至完成集合PLk中所有点的Z坐标值计算，然后执行步骤(6-9)；
(6-9)循环执行步骤(6-6)到(6-8)，直至集合BL中所有边界线被遍历，得到地貌单元gui所有边界线的点的Z坐标值；
(6-10)将包含有Z坐标值的集合PR中的点和所有集合PLk中的点存储到集合VD中，作为地貌单元gui的虚拟钻孔集合。
(7)基于埋深等值线、虚拟钻孔集合VD、钻孔集合RD和Dem1，构建地貌单元gui的基岩面三维模型bmi。
该步骤具体包括：
(7-1)若地貌单元gui为平原，执行步骤(7-5)；否则将基岩埋深等值线上的点存储到集合PC中；
(7-2)从集合PC中获取任一点pcb，获取点pcb在地貌单元gui的DEM数据DEMi上的地表高程值hp，并基于基岩埋深等值线，获取点pcb的埋深值dp，计算高程值hp与埋深值dp的差值difp，作为点pcb的Z坐标值；重复执行该步骤直至完成集合PC中所有点的Z坐标值计算；
(7-3)从钻孔集合RD中获取任一钻孔rdj，读取rdj的高程属性得到高程值hrd，以及埋深属性得到埋深值drd；
(7-4)获取点rdj的X坐标值和Y坐标值，生成点P，计算高程值hrd与埋深值drd的差值difrd，作为点P的Z坐标值，将点P存储到集合PC中；
(7-5)循环执行步骤(7-3)至步骤(7-4)，直至集合RD中所有点数据都被处理；
(7-6)将地貌单元gui的虚拟钻孔集合VD中所有的点都存储到集合PC中；
(7-7)以集合PC为基础数据通过反距离加权插值方法得到栅格面bmi，栅格面bmi即为地貌单元gui的基岩面三维模型。
(8)重复执行步骤(3)至步骤(7)，直至所有地貌单元的基岩面三维模型构建完成，并存储至基岩面模型集合BM中。
(9)将Dem2和基岩面模型集合BM中所有栅格合并，生成研究区域完整的基岩面三维模型，如图5和图6所示。
以上所揭露的仅为本发明一种较佳实施例而已，不能以此来限定本发明之权利范围，因此依本发明权利要求所作的等同变化，仍属本发明所涵盖的范围。

Claims (7)

1.一种知识驱动的基岩面三维模型构建方法，其特征在于该方法包括：
(1)基于研究区的地质图，获取基岩出露区边界线，并对研究区的DEM数据进行裁剪，获得基岩埋藏区域的DEM数据dem1和基岩出露区的DEM数据dem2；
(2)根据研究区的地貌单元图层，将基岩埋藏区域划分为岗地、河谷、湖泊和平原这四类地貌单元类型，并存入地貌单元集合GU中；
(3)从地貌单元集合GU中获取任一地貌单元gui，若gui为平原，则执行步骤(6)；否则，执行步骤(4)；
(4)获取gui的钻孔集合RD，并根据钻孔集合RD和gui的地貌单元类型计算gui待构建埋深等值线的最大埋深值Hmax和最小埋深值Hmin，之后执行步骤(5)；
(5)基于Hmax和Hmin，构建地貌单元gui的埋深等值线，执行步骤(6)；
(6)基于Dem1，构建地貌单元gui的虚拟钻孔集合VD；
(7)基于埋深等值线、虚拟钻孔集合VD、钻孔集合RD和Dem1，构建地貌单元gui的基岩面三维模型bmi；
(8)重复执行步骤(3)至步骤(7)，直至所有地貌单元的基岩面三维模型构建完成，并存储至基岩面模型集合BM中；
(9)将Dem2和基岩面模型集合BM中所有栅格合并，生成研究区域完整的基岩面三维模型；
步骤(7)具体包括：
(7-1)若地貌单元gui为平原，执行步骤(7-5)；否则将基岩埋深等值线上的点存储到集合PC中；
(7-2)从集合PC中获取任一点pcb，获取点pcb在地貌单元gui的DEM数据DEMi上的地表高程值hp，并基于基岩埋深等值线，获取点pcb的埋深值dp，计算高程值hp与埋深值dp的差值difp，作为点pcb的Z坐标值；重复执行该步骤直至完成集合PC中所有点的Z坐标值计算；
(7-3)从钻孔集合RD中获取任一钻孔rdj，读取rdj的高程属性得到高程值hrd，以及埋深属性得到埋深值drd；
(7-4)获取点rdj的X坐标值和Y坐标值，生成点P，计算高程值hrd与埋深值drd的差值difrd，作为点P的Z坐标值，将点P存储到集合PC中；
(7-5)循环执行步骤(7-3)至步骤(7-4)，直至集合RD中所有点数据都被处理；
(7-6)将地貌单元gui的虚拟钻孔集合VD中所有的点都存储到集合PC中；
(7-7)以集合PC为基础数据通过反距离加权插值方法得到栅格面bmi，栅格面bmi即为地貌单元gui的基岩面三维模型。
2.根据权利要求1所述的知识驱动的基岩面三维模型构建方法，其特征在于：步骤(4)具体包括：
(4-1)获取gui内的所有钻孔数据到集合RD＝{rdj|j＝1,2,…,nrd}，其中，rdj表示gui内第j个钻孔，nrd为钻孔个数；
(4-2)获取埋深值最大的钻孔rdm，将其埋深值记为Hm，根据下式计算最大埋深值Hmax；
式中，d为预设埋深线间距，符号[]为取整符号；
(4-3)若gui为岗地，则直接将Hn赋值为零；否则，获取埋深值最小的钻孔rdn，将其埋深值记为Hn；
(4-4)基于Hn采用下式，计算最小埋深值Hmin；
(4-5)判断是否满足Hmin＝d，若是，则执行步骤(5)；否则执行步骤(4-6)；
(4-6)根据钻孔rdm和钻孔rdn分别到地貌单元边界的最短距离调整最小值Hmin，并执行步骤(5)。
3.根据权利要求2所述的知识驱动的基岩面三维模型构建方法，其特征在于：步骤(4-6)具体包括：
(4-6-1)分别获取钻孔rdm、钻孔rdn到地貌单元边界的最短距离Dism、Disn；
(4-6-2)根据下式计算基岩埋深为Hmin-d的等值线到地貌单元边界的距离Dis；
式中，Hn、Hm分别为钻孔点rdn、rdm的基岩埋深值,符号||为取绝对值符号；
(4-6-3)若Dis≥0，则令Hmin＝Hmin-d，然后执行步骤(4-5)；否则，执行步骤(5)。
4.根据权利要求1所述的知识驱动的基岩面三维模型构建方法，其特征在于：步骤(5)具体包括：
(5-1)基于地貌单元gui的钻孔集合RD，获取用于构建基岩埋深为Hmin的等值线ch1的钻孔rdu和rdv；
(5-2)若gui为岗地，则基于地质图，将gui的基岩出露线作为埋深基准线ch0，否则，将gui地貌单元边界线作为埋深基准线ch0；
(5-3)分别计算钻孔点rdu、rdv到基岩埋深基准线ch0的最短距离Disu、Disv；
(5-4)根据下式，计算基岩埋深为Hmin的等值线ch1到基岩埋深基准线ch0的距离Dis1-0；
式中，Hu、Hv分别为钻孔点rdu、rdv的基岩埋深值,符号||为取绝对值符号；
(5-5)基于基岩埋深基准线ch0和距离Dis1-0，执行缓冲区分析，获取缓冲区边界线BLine0；
(5-6)根据地貌单元gui的范围裁剪BLine0，生成基岩埋深为Hmin的等值线ch1；
(5-7)循环执行(5-1)至(5-6)，依次构建基岩埋深为Hmin+d，Hmin+2d，…，Hmax的基岩埋深等值线，直至地貌单元gui的所有基岩埋深等值线构建完成。
5.根据权利要求4所述的知识驱动的基岩面三维模型构建方法，其特征在于：步骤(5-1)具体包括：
(5-1-1)获取gui的钻孔集合RD，并根据埋深值，将钻孔集合RD分为子集合RD1和RD2，其中子集合RD1中钻孔的埋深值均不大于Hmin，子集合RD2中钻孔的埋深值均大于Hmin；
(5-1-2)若RD1和RD2均不为空，则获取RD1中埋深值最大的钻孔，记为rdu，其基岩埋深记为Hu，获取RD2中埋深值最小的钻孔，记为rdv，其基岩埋深记为Hv，然后执行步骤(5-2)；否则，执行步骤(5-1-3)；
(5-1-3)若RD1为空集，获取RD2中埋深值最小的两个钻孔作为rdu、rdv；若RD2为空集，则获取RD1中埋深值最大的两个钻孔作为rdu、rdv。
6.根据权利要求1所述的知识驱动的基岩面三维模型构建方法，其特征在于：步骤(6)具体包括：
(6-1)获取地貌单元gui的钻孔集合RD，以钻孔集合RD为基础数据通过反距离加权插值方法得到栅格面DDem；
(6-2)从基岩埋藏区域的DEM数据dem1中获取地貌单元gui的DEM数据DEMi；
(6-3)构建地貌单元gui的外包矩形Rec，并在矩形Rec内部生成10×10网格，将位于gui的网格交叉点存储到点集合PR中；
(6-4)从点集合PR中获取任一点prp，获取点prp在DDem上的埋深值dr，以及点prp在DEMi上的高程值hr，计算dr和hr的差值difr，并作为点prp的Z坐标值；重复执行本步骤，直至完成集合PR中所有点的Z坐标值计算；
(6-5)读取地貌单元gui的边界线到集合BL＝{blk|k＝1,2,…,nbl}，其中，blk表示第k条边界线，nbl为边界线个数；
(6-6)从集合BL中获取任一边界线blk，将边界线blk离散为点集合PLk＝{plk,l|l＝1,2,…,npk}，其中，plk,l表示边界线blk上第l个点，npk为边界线blk上的点个数；若blk为基岩出露边界线，则执行步骤(6-7)，否则执行步骤(6-8)；
(6-7)从点集合PLk中获取任一点plk,l，获取点plk,l在DEMi上的高程值hl，作为点plk,l的Z坐标值；重复执行该步骤直至完成集合PLk中所有点的Z坐标值计算，然后执行步骤(6-9)；
(6-8)从点集合PLk中获取任一点plk,l，获取点plk,l在DDem上的埋深值dl，以及点plk,l在DEMi上的高程值hl，计算dl和hl的差值difl，作为点plk,l的Z坐标值；重复执行该步骤直至完成集合PLk中所有点的Z坐标值计算，然后执行步骤(6-9)；
(6-9)循环执行步骤(6-6)到(6-8)，直至集合BL中所有边界线被遍历，得到地貌单元gui所有边界线的点的Z坐标值；
(6-10)将包含有Z坐标值的集合PR中的点和所有集合PLk中的点存储到集合VD中，作为地貌单元gui的虚拟钻孔集合。
7.一种知识驱动的基岩面三维模型构建装置，包括处理器及存储在存储器上并可在处理器上运行的计算机程序，其特征在于：所述处理器执行所述程序时实现权利要求1-6中任意一项所述的方法。
CN202010794872.5A 2020-08-10 2020-08-10 知识驱动的基岩面三维模型构建方法及装置 Active CN111951396B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010794872.5A CN111951396B (zh) 2020-08-10 2020-08-10 知识驱动的基岩面三维模型构建方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010794872.5A CN111951396B (zh) 2020-08-10 2020-08-10 知识驱动的基岩面三维模型构建方法及装置

Publications (2)

Publication Number Publication Date
CN111951396A CN111951396A (zh) 2020-11-17
CN111951396B true CN111951396B (zh) 2024-02-27

Family

ID=73333074

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010794872.5A Active CN111951396B (zh) 2020-08-10 2020-08-10 知识驱动的基岩面三维模型构建方法及装置

Country Status (1)

Country Link
CN (1) CN111951396B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113393576B (zh) * 2021-05-27 2024-07-05 南京师范大学 基于地质图的松散层三维模型构建方法及装置

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109003330A (zh) * 2018-07-02 2018-12-14 南京师范大学 一种基于基岩边界约束的三维地层建模方法
CN110197013A (zh) * 2019-05-14 2019-09-03 南京泛在地理信息产业研究院有限公司 基于Morphing的河床基岩面建模方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP3788378B2 (ja) * 2002-03-27 2006-06-21 ソニー株式会社 ３次元モデル生成システム及び方法、並びにコンピュータ・プログラム

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109003330A (zh) * 2018-07-02 2018-12-14 南京师范大学 一种基于基岩边界约束的三维地层建模方法
CN110197013A (zh) * 2019-05-14 2019-09-03 南京泛在地理信息产业研究院有限公司 基于Morphing的河床基岩面建模方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
虚拟钻孔控制的三维地质体模型构建方法;林冰仙;周良辰;闾国年;;地球信息科学学报;20131015(05);全文 *

Also Published As

Publication number Publication date
CN111951396A (zh) 2020-11-17

Similar Documents

Publication Publication Date Title
AU2015312396B2 (en) Sediment transport simulation with parameterized templates for depth profiling
US20130338978A1 (en) Generating facies probablity cubes
CN103116192A (zh) 一种缝洞型碳酸盐岩油藏储集体建模方法
CN110056346B (zh) 一种基于趋势变化函数的油藏三维原始含水饱和度模拟方法
CN107329171A (zh) 深度域储层地震反演方法及装置
CN108986213B (zh) 一种基于叠置技术的三维地层建模方法
CN109138975B (zh) 一种基于时移测井数据的求解相渗特征曲线的新方法
CN110955982A (zh) 一种变质岩储层渗透率计算方法、装置及计算机存储介质
CN111475920A (zh) 一种深水盆地古水深的获取方法、系统、电子设备及存储介质
CN110501757A (zh) 点线面三维逐点井深设计方法
CN109884700A (zh) 多信息融合地震速度建模方法
CN106646603B (zh) 一种实际时移地震资料处理差异的可靠性判断方法
CN108150160A (zh) 一种地层中欠压实作用超压的求取方法
Tikhomirova et al. A review on methods of oil saturation modelling using IRAP RMS
CN111951396B (zh) 知识驱动的基岩面三维模型构建方法及装置
CN110344822B (zh) 含钙质致密储层含水饱和度的确定方法
CN110298473A (zh) 一种超采区地下水库选址及确定调蓄库容的方法
CN103543478A (zh) 地质形态插值的km方法
CN112100878A (zh) 一种利用沉积正演模拟建立时间域层序地层剖面的方法
CN113775326A (zh) 可动水饱和度的评价方法、装置、电子设备及介质
Ghassemi et al. Three-dimensional simulation of the Home Island freshwater lens: preliminary results
CN116629023A (zh) 一种融合地层正演模拟和地质统计学的三维油藏表征新方法
Szabó et al. Rock densities in the Pannonian basin-Hungary
CN115327627A (zh) 多信息融合的致密砂岩气分流河道表征方法及装置
CN112649880B (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