一种基于平面地质图的背斜构造三维建模方法技术

本发明专利技术公开了一种基于平面地质图的背斜构造三维建模方法，属于地质领域，用于解决背斜构造需要人为花费大量时间采集数据和实地勘测才将背斜的空间形态呈现出来，此类方法费时费力，且无法将背斜上地貌特征一一呈现的问题，背斜构造三维建模方法具体包括：栅格数据处理；背斜褶皱的智能判定；实体空间剖分与建模；剥蚀量计算；本发明专利技术是将背斜构造的空间形态以及背斜上地貌特征快速且准确进行呈现。态以及背斜上地貌特征快速且准确进行呈现。态以及背斜上地貌特征快速且准确进行呈现。

【技术实现步骤摘要】
一种基于平面地质图的背斜构造三维建模方法


[0001]本专利技术属于地质领域，涉及背斜构造三维建模技术，具体是一种基于平面地质图的背斜构造三维建模方法。

技术介绍

[0002]背斜指岩层发生褶曲时，形状向上凸起者，在一般平地上，背斜的地层上半部受到侵蚀变平，会形成中间古老，两侧较新的地层排列方式，背斜与向斜具有重要的实践意义，由于背斜岩层向上拱起，且油、气的密度比水小，所以背斜常是良好的储油、气构造，背斜因岩体向上拱起，有支撑作用，且不易渗水，所以适合在下方挖隧道；在现有技术中，背斜构造需要人为花费大量时间采集数据和实地勘测，从而才能够将背斜的空间形态呈现出来，此类方法费时费力，且无法将背斜上剥蚀等特征一一呈现，为此，我们提出一种基于平面地质图的背斜构造三维建模方法。

技术实现思路

[0003]针对现有技术存在的不足，本专利技术目的是提供一种基于平面地质图的背斜构造三维建模方法。
[0004]本专利技术所要解决的技术问题为：如何将背斜构造的空间形态以及背斜上地貌特征快速且准确进行呈现。
[0005]本专利技术的目的可以通过以下技术方案实现：一种基于平面地质图的背斜构造三维建模方法，背斜构造三维建模方法具体包括：步骤S1：栅格数据处理；步骤S2：背斜褶皱的智能判定；步骤S3：实体空间剖分与建模；步骤S4：剥蚀量计算。
[0006]进一步地，所述栅格数据处理过程具体如下：步骤S11：提取地形等高线；根据地形地质图上的地形特征提取地形等高线，并保存为矢量格式，比如mapgis的wl格式或者shapefile的线文件存储，并赋予高程属性；步骤S12：构造数字高程模型；通过步骤S11获取的地形等高线矢量文件，结合用户给定的网格参数，进行网格化处理，生成数字高程模型，并保存为geotiff格式；步骤S13：提取岩层属性和产状；根据地形地质图给定的地质信息，将每一个独立的地质单元进行矢量化，并赋予地层属性，并保存为mapgis的wp或者shapefile的面文件；对地形地质图中给定的产状位置进行矢量化，并保存为mapgis的wt或者
shapefile的点文件，赋予对应的产状和岩层属性；步骤S14：提取地质界线；根据地形地质图给定的地质信息，将岩层接线进行矢量化，并赋予界线属性。
[0007]进一步地，网格参数为像元地面分辨率；界线属性为界线两侧的岩层属性。
[0008]进一步地，所述背斜褶皱的智能判定过程具体如下：步骤S21：空间要素的初始化；建模范围：Bbox=｛
‘
west
’
：west，
’
east
’
：east，
’
south
’
：south，
’
north
’
：north｝；数字高程模型：Dem.tif；地质单元：units=｛u_1：｛
‘
name
’
：name_1，
’
era
’
：era，
’
points
’
：points_1｝，
……
｝；地质界线：boundaries=｛b_1：｛
‘
name
’
：name_1，
’
points
’
：points_1，
’
left
’
：u_m，
’
right
’
：u_n｝，
……
｝；地质产状：occurence=｛u_1：[[x，y，dip，angle]，
……
]，
……
｝；步骤S22：背斜关键标志；步骤S23：核部与两翼确定；根据背斜构造的定义，根据获得的地质界线可作为背斜褶皱的基本界线要素，位于集合bdy3的首位的界线为背斜两翼的最外侧界线ly=[b01，b10]，位于bdy3中心的地质界线为背斜核部c=[b0n，bn0]。
[0009]进一步地，所述步骤S22的分析过程具体如下：步骤S221：地质界线走向一致性；遍历步骤S21获得的地质界线结果，计算所有的地质界线的走向，并筛选走向一致的地质界线；bdy1=[b1，b2，b3
……
，bn]；步骤S222：地质单元时代
‑
新
‑
老
‑
新排列唯一性；遍历步骤S21获得的走向一致性地质界线bdy1，检测界线两侧地质单元时代，确定地质单元时代排列遵循新
‑
老
‑
新的变化规则，去除不符合约束的地质界线，获得结果；bdy2=[b01，b02，b03，
……
b0n]；步骤S223：地质单元排列对称性；遍历步骤S222获取的结果bdy2，检测地质界线两侧的地质单元在空间上是否呈现排列对称性，去除不符合约束条件的地质界线后得到构造背斜的界线要素；bdy3=[b01，b02，b03，
……
b0n，bn0
……
，b30，b20，b10]。
[0010]进一步地，b01右侧和b10左侧的岩层倾向呈近似轴对称，且属于同一岩层单位，b01和b10位背斜褶皱建模的边界约束。
[0011]进一步地，所述实体空间剖分与建模的工作过程具体如下：步骤S31：背斜约束边界计算；根据两翼的地质界线ly，读取b01和b10边界约束地层的产状，分别读取b01和b10地质界线的K01、L01和K10、L10两点的坐标；并根据b01.left和b10.right的地层产状构造平面K01L01Q1P1和平面
K10L10Q1P1，K10L10Q1P1平面与地质界线b01左侧的岩层B倾向、倾角一致，K01L01Q1P1平面与地质界线b02右侧的岩层B倾向、倾角一致；平面K01L01Q1P1和平面K10L10Q1P1相交于直线P1Q1，计算Q1P1点的坐标Q1（Q1x，Q1y，Q1z）和P1（P1x，P1y，P1z），比较Q1点和P1点高程，取较高值得点作为模型顶点；面K01L01Q1P1和面K10L10Q1P1构成背斜的外围约束边界；步骤S32：岩层约束界面构造；按照步骤S31获取背斜外围约束边界的方法，读取地质界线b02和b20的K2、L2和M2、N2坐标，读取b02右侧和b20左侧的岩层产状；分别根据倾向和倾角构造面K02L02Q2P2、K20L20Q2P2，两个面相交于直线P2Q2，分别计算岩层C的约束界面控制点坐标Q2（Q2x，Q2y，Q2z）和P2（P2x，P2y，P2z）；运用同样的方法，根据bdy3中剩余地质界线的参数和空间属性，直至完成岩层N约束界面控制点坐标Qn（Qnx，Qny，Qnz）和Pn（Pnx，Pny，Pnz）；步骤S33：立方体剖分与剥蚀地层三维重建；根据设定的最小剖分单元立方体边长，以地质界线b01的L01作为起点，依据数字高程模型dem.tif及地层单元units和地质界线bdy3本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种基于平面地质图的背斜构造三维建模方法，其特征在于，背斜构造三维建模方法具体包括：步骤S1：栅格数据处理；步骤S2：背斜褶皱的智能判定；步骤S3：实体空间剖分与建模；步骤S4：剥蚀量计算。2.根据权利要求1所述的一种基于平面地质图的背斜构造三维建模方法，其特征在于，所述栅格数据处理过程具体如下：步骤S11：提取地形等高线；根据地形地质图上的地形特征提取地形等高线，并保存为矢量格式，同时赋予高程属性；步骤S12：构造数字高程模型；通过步骤S11获取的地形等高线矢量文件，结合用户给定的网格参数，进行网格化处理，生成数字高程模型，并保存为geotiff格式；步骤S13：提取岩层属性和产状；根据地形地质图给定的地质信息，将每一个独立的地质单元进行矢量化，并赋予地层属性，并保存为mapgis的wp或者shapefile的面文件；对地形地质图中给定的产状位置进行矢量化，并保存为mapgis的wt或者shapefile的点文件，赋予对应的产状和岩层属性；步骤S14：提取地质界线；根据地形地质图给定的地质信息，将岩层接线进行矢量化，并赋予界线属性。3.根据权利要求2所述的一种基于平面地质图的背斜构造三维建模方法，其特征在于，网格参数为像元地面分辨率；界线属性为界线两侧的岩层属性。4.根据权利要求1所述的一种基于平面地质图的背斜构造三维建模方法，其特征在于，所述背斜褶皱的智能判定过程具体如下：步骤S21：空间要素的初始化；建模范围：Bbox＝{
‘
west
’
：west，
’
east
’
：east，
’
south
’
：south，
’
north
’
：north}；数字高程模型：Dem.tif；地质单元：units＝{u_1：{
‘
name
’
：name_1，
’
era
’
：era，
’
points
’
：points_1}，
……
}；地质界线：boundaries＝{b_1：{
‘
name
’
：name_1，
’
points
’
：points_1，
’
left
’
：u_m，
’
right
’
：u_n}，
……
}；地质产状：occurence＝{u_1：[[x，y，dip，angle]，
……
]，
……
}；步骤S22：背斜关键标志；步骤S23：核部与两翼确定；根据背斜构造的定义，根据获得的地质界线可作为背斜褶皱的基本界线要素，位于集合bdy3的首位的界线为背斜两翼的最外侧界线ly＝[b01，b10]，位于bdy3中心的地质界线为背斜核部c＝[b0n，bn0]。5.根据权利要求4所述的一种基于平面地质图的背斜构造三维建模方法，其特征在于，
所述步骤S22的分析过程具体如下：步骤S221：地质界线走向一致性；遍历步骤S21获得的地质界线结果，计算所有的地质界线的走向，并筛选走向一致的地质界线；bdy1＝[b1，b2，b3
……
，b...

【专利技术属性】
技术研发人员：吴勇庆，曾晓剑，李洋，
申请(专利权)人：广东省粤东地质工程勘察有限公司，
类型：发明
国别省市：