EN中文

Journal of Geo-information Science >

2016 , Vol. 18 >Issue 10: 1322 - 1331

DOI: https://doi.org/10.3724/SP.J.1047.2016.01322

Orginal Article

The Research of Expression Method on Geological Fault Modeling

  • SUN Jingguang , 1 ,
  • GAO Tianpeng , 2, *
Expand
  • 1. School of Electronic and Information Engineering, LiaoNing Technical University, Huludao 125000, China
  • 2. School of Graduate, LiaoNing Technical University, Huludao 125000, China
*Corresponding author: GAO Tianpeng, E-mail:

Received date: 2015-09-21

  Request revised date: 2015-10-20

  Online published: 2016-10-25

Copyright

《地球信息科学学报》编辑部 所有
Fold

Abstract

The limitations of expressing faults, the interpolation dilemma and other problems exist in the 3D modeling and the visualization technique of geological faults. Through our research, the reason why these problems exist is that these methods mainly focus on the geometrical morphology of the geological faults and do not concern the formation mechanism of the geological faults. As we all know, no matter how complex the geologic bodies that containing the geological faults are, the formation mechanisms which they depend on are similar. Besides, the modeling of the geologic faults is composed of a series of processes. If we can model the faults through certain ways to express their formation mechanisms, then we may solve the above mentioned problems and limitations. Thus, we proposed an expression method to express the formation mechanism for modeling and visualizing the geological faults. Firstly, we analyzed the operands and operators needed to construct the geological fault expression, constructed the rules for the expression based on the theory of formation mechanism of the geological faults, and used the context-free grammar rules to describe the rules for fault expression based on the formation mechanism. Then, we formated the original fault exploration data to the standardize 2D data table, extracted the operators and operands from the 2D data table to construct the expression of geological faults. Finally, we calculated the expression fomular to generate the abstract faults model and achieved the three-dimensional geological modeling with the relevant exploration data based on the abstract faults model. Experiment is carried out based on DaLian River geological data which comprising a plurality of geological faults. Geological modeling results show that the expression method solves the limitations of fault expression, and avoids the interpolation dilemma accured in the geological layer fracture zone.

Key words： geological fault; expression of fault; rules of expression; abstract fault model; visualization of fault

Cite this article

SUN Jingguang , GAO Tianpeng . The Research of Expression Method on Geological Fault Modeling[J]. Journal of Geo-information Science, 2016 , 18(10) : 1322 -1331 . DOI: 10.3724/SP.J.1047.2016.01322

1 引言

断层是岩层或岩体顺破裂面发生明显的位移现象[1]。断层是在地壳中分布极为广泛的地质构造类型,大型断层控制和影响成矿作用,中小型断层控制矿床、矿体的产状,活动断层则直接影响工程施工和建筑物的稳定性。因此,研究地质断层的三维建模技术对分析地质断层结构、描述运动规律等具有重要的意义和实用价值。
地质三维建模及可视化的难点[2-3]主要体现在以下方面：① 地质断层数据类型多样、数据稀疏、获取艰难,导致断层建模数据提取复杂;② 断层通常以组合形式存在,导致形态多变,用现有的数据模型很难完整地描述与表达;③ 断层构造引起的地层不连续导致地层面重构[4-5]困境。针对上述问题,许多学者深入地研究了如下3类建模方法[3]：① 基于地层恢复的断层建模技术,又称为整体法建模技术,其首先假设断层不存在,根据所有的采样点数据进行三维地质建模获得没有断层的三维地质模型,然后提取断层数据,加入断层,恢复到地质体发生断裂后的真实模型。该方法有效地解决了断层数据较少、断距不大的地质体的断层建模。比较适合于断距较小、断层面较简单的断层建模。② 基于分区插值[6-8]的断层建模技术,又称为局部法建模技术,该技术首先根据断层数据把整个地质体分为不同的区块,然后分别对各个区块进行三维建模,最后按照一定的关系进行组合,形成一个完整的三维地质模型。该方法对采样点数据较多、断距较大的地质体能够准确地建模。③ 断层与地层统一构模技术[3]。根据数据特点,首先对每一个地层进行整体拟合,同时考虑断层线的约束作用,从而生成合理的地层层面;然后将断层面与地层面进行切割处理,调整交线上的点,并将其约束到断层面上;最后根据上、下层面及侧面形成完整的实体模型。
综上,局部法与整体法能够很好地应用于特定类型的断层,但具有明显的局限性问题（即不能准确模拟所有类型的断层）,而且局部法存在建模过程复杂,同时要求建模人员具有较好的地质学背景。断层与地层统一构模技术理论上能够适用于任意类型的断层,但含有断层的地质体往往形态复杂,而统一建模技术需要运用插值技术对包含断层的地层进行整体拟合,依然存在插值因境问题,插值困境即运用插值技术对地层断开裂处插值,而插值技术较适用于无断裂区域的地层,因而造成地层断裂区域插值的不准确性;另外,统一建模技术还需进行断层面与地层面的交割处理,具有一定的复杂性。
针对上述方法存在的问题,本文基于局部法建模技术的建模思路,提出从构建断层抽象模型入手,基于断层抽象模型进行地质断层三维建模及可视化表达。首先,分析建立表达式所需要的算子类型、对象类型以及表达式构建的方法与规则;然后,对断层数据进行规范化,形成构建表达式所需要的抽象断层表格数据,进而从中提取出构造断层表达式的算子序列和对象集合,形成断层表达式;最后,对表达式进行计算,形成断层抽象模型,结合勘探数据进行地质断层三维建模及可视化。采用本文提出的方法不仅能够解决整体法与局部法的局限性,而且避免了在地层断裂处进行插值以及地层面与断层面的交割处理,从而解决了插值困境并有效地降低了复杂度。
Tab.1 Operators of fault expression method

表1 断层表达式方法算子集合

算子 目数 结合性 优先级 类型
u(Up) 单 左 1 几何算子
d(Down) 单 左 1 几何算子
t(Translate) 单 左 1 几何算子
r(Rotate) 单 左 1 几何算子
p(Topology) 双 右 2 拓扑算子
=(Assignment) 双 左 2 混合算子

2 表达式方法建模的基本理论

（1）表达式方法
表达式[9]是说明计算的基本方式,其大部分特征都来自数学中的传统形式,主要包括运算对象（操作数）和算子（操作符）。具体描述如下：运算对象为算子作用的对象;算子指定对运算对象进行何种性质的操作;目数指算子作用于对象的个数;优先级指算子进行计算的顺序。
（2）表达式规则
表达式规则即规定表达式中的运算对象类型、算子集合及运算对象与算子的运算方式。运算对象的类型决定了运算对象的运算方式及哪些算子能够作用于运算对象。运算对象与算子的运算方式决定于算子的目数、优先级及结合性。
（3）上下文无关文法
上下文无关文法[10]（context-free grammar）可以用来对表达式方法规则的语法结构进行精确无误的形式描述,进而依据这种形式描述构造表达式文法规则分析器。若表达式经分析器分析为合法,则进行后续表达式的计算。一个上下文无关文法由终结符号集合、非终结符号集合、产生式集合和开始符号4个元素组成。

3 断层表达式规则建模思路

基于断层的形成机制[11]：原本连续的地质体在地壳运动的作用下,地质体的某一部分或某几部分沿破裂面发生明显位移或扭曲而形成断层。位移包括地质体沿破裂上升、下降、水平移动;扭曲是绕断裂面上某一轴旋转一定的角度。无论多么复杂的包含断层的地质体,具有相同的形成机制。从断层的形成机制出发来表达和模拟断层,可以适用于不同类型的断层,避免局部法建模技术和整体法建模技术的局限性。基于这种考虑,提出表达式方法的断层建模技术。
根据断层的形成机制,断层表达式规则建模方法的基本思想：首先定义算子集合和算子作用的对象类型,然后规定每个算子的目数、结合性、优先级的目数,最后制定表达式的规则。断层表达式规则研究主要思路如下：
（1）断层表达式运算对象的提取
经过对断层形成机制进行分析,运算对象的类型有六面体（G）和断层抽象模型（F）2种子地质体（图1）。子地质体是一个抽象六面体（该六面体每个面可以是平面或曲面,而平面或者曲面在特殊情况下可以是一条线段或曲线段）;断层抽象模型由2个子地质体通过1个共同的断层面组成,而组成断层抽象模型的每个子地质体本身也可以是一个断层抽象模型。
Fig.1 Sub geological unit and abstract fault model

图1 子地质体和断层抽象模型

（2） 断层表达式算子的提取
经过分析断层表达式的算子集合如表1所示。根据算子作用于运算对象属性的类型将算子集合分为以下3类：①几何运算算子：该类算子实现对运算对象的几何属性（即空间位置等）的修改。包括 u 、 d 、 t 、 r 。②拓扑运算算子：该类算子实现对运算对象的拓扑属性（即连接关系等）的修改,包括 p 。③混合运算算子：该类算子实现对几何属性或拓扑属性的修改,包括“=”。
（3） 制定断层表达式的运算规则
通过分析,总结每个算子作用的对象类型及如何作用于运算对象,并运用上下文无关文法（context-free grammar）对表达式规则进行描述。

4 断层表达式规则三维建模的方法

断层表达式规则三维建模方法包括：（1）断层数据规范化,主要用来将地质断层数据进行处理,整理成构建表达式所需要的断层数据表。（2）构建断层表达式生成断层地质体抽象模型,以断层数据规范化步骤所生成二维数据表。首先,生成算子所作用的对象集和算子序列;然后,以算子和运用对象的映射关系生成断层表达式;最后,对断层表达式进行生成断层地质体抽象模型。

4.1 断层数据规范化

断层数据的来源主要分5类[12-14]：地质勘探获取的断层参数数据、钻孔中的断点数据、地质剖面数据、平面地质图、地形图数据。这些数据可共同描述断层的几何形态,但格式不一,为断层构模技术的研究带来不便。因此,将获取的断层相关不同类别数据整理为标准表格形式,是后续断层三维建模及可视化步骤的前提。地质勘探获取的断层参数据主要包括断层的走向、倾角、倾向、断距、类型等。从地质剖面图、钻孔中获取断层的断点数据,从平面地质图、地形图中获取断层线数据。结合专家经验,以及表达式方法断层三维建模及可视化的需要,将数据规范化为相应的二维表（表2）。
Tab.2 Standardization of two-dimensional tables

表2 规范化二维表

表名 字段名
断层描述表 编号、类型、优先级、起始位置、终止位置、倾向、 倾角、位移、旋转扩展、标志位、其他
层线数据表 编号、所属断层
编号、采样点集数据表指针、起始点、终止点、其它
点数据表 编号、所属断层编号、坐标、其它
子地质体数据表 编号、所属断层编号、面集合、其它、编号
采样点集数据表 所属断层线编号、坐标、其它

4.2 构建断层地质体抽象模型

断层地质体抽象模型（即概念模型,因为该步骤没有对地质体进行地层划分及地层面拟合）的构建流程如图2所示。
Fig.2 The process of creating the abstract fault model

图2 断层地质体抽象模型的构建流程

4.2.1 算子的描述
（1）u（Up）
功能：该算子作用于子地质体（G）,实现子地质体沿断层面向上平移的操作。
输入：断层描述表的位移、倾向
构建：位移在倾向的投影
（2）d（Down）
功能：该算子作用于子地质体（G）,实现子地质体沿断层面向下平移的操作。
输入：断层描述表的位移、倾向
构建：位移在倾向的投影
（3）t（Translate）
功能：该算子作用于子地质体（G）,实现子地质体沿断层面向上平移的操作。
输入：断层描述表的位移、倾向
构建：位移在倾向的投影
（4）r（Rotate）
功能：该算子作用于子地质体（G）,实现子地质体沿断层面的旋转轴旋转一定角度的操作。
输入：断层描述表的位移、倾向
（5）p（Topology）
功能：该算子作用于子地质体（G）,将子地质体加入断层抽象模型的拓扑结构中。
输入：子地质体A和子地质体B
（6）=（Assignment）
功能：该算子作用于子地质体（G）及断层抽象模型（F）,实现更新F、G的几何属性（即空间位置）和拓扑属性。
4.2.2 断层表达式规则
运用上下文无关文法（context-free grammar）对断层表达式规则的描述如下：
（1）T（终结符集合）
包括u、d、t、r、p、=、F、G,其中u、d、t、r、p、=分别代表断层表达式的算子,F、G代表算子作用的对象。
（2）N（非终结符集合）
包括Factor、Expr、Item、Term、Right,其中Factor代表断层表达式的基本单元,Expr、Term分别对应对于2个优先级层次,Right对应于右结合性算子。
（3）P（产生式集合）
Right → Term = Right | Term （1）
Term → Item | Term.pItem （2）
Item → F | Factor （3）
Expr → uFactor | dFactor | tFactor | rFactor | Factor （4）
Factor → G ( Expr ) （5）
其中,Right作为断层表达式的开始符号。
4.2.3 生成断层表达式的运算对象
断层表达式的运算对象有子地质体六面体（G）和断层抽象模型（F）2种。若要生成断层抽象模型（F）对象,首先要生成子地质体六面体（G）对象。通过抽象断层线及相应断层的倾角、倾向即可生成抽象断层面,一个抽象断层面将地质体切分为2个子地质体六面体（G）对象,而这2个子地质体就组成了一个断层抽象模型对象（F）。具体过程如下：
（1）提取抽象断层线。该步骤提取断层线数据表每条记录中的起始点和终结点构成抽象断层线(即为仅包含起始点和终止点的断层线)。
（2）断层线相交处理。研究区域可能包含单条或多条断层线。如果包含多条断层线则断层线可能出现相交情形,因而需要进行相交检测并划分子断层线。
（3）构建断层面。结合断层线数据表记录中的所属断层编号字段,提取断层描述表的相应记录中的倾角、倾向以及断点数据表并结合抽象断层线或子断层线构建断层面集合。
（4）断层面切分初始六面体生成运算对象,并建立对象间的拓扑关系。断层面集合将初始六面体地质体划分为子地质体六面体对象（G）集合和断层抽象模型对象（F）集合。步骤如下：
① 以断层描述表中的优先级字段确定切分次序即确定断层面或断层面子集对子地质体的切分次序,并形成有序的面集合 L 。
② 对步骤①中面集合 L 遍历并提取 L 集合中每个断层面元素对应断层的倾角及倾向,对初始六面体或子地质体时行切分。在切分的同时建立子地质体对象间的拓扑关系,并将生成的子地质体对象信息存入子地质体信息表。
③重复步骤②直到所有的 L 集合中所有元素遍历完为止。
下面以存在2条断层 F ( a ) 与 F ( b ) 的地质体 F 为例说明,其可能有以下3种情况：
（1）断层线不相交（图3）
该情况下断层抽象模型对象集合 F 及子地质体对象的集合 G 分别为式（6）-（7）。
F = { F ( a ) , F ( b ) } （6）
G = { G ( a ) , G ( b ) , G ( c ) } （7）
式中： G ( a ) 、 G ( b ) 与 F ( a ) 建立拓扑关系, G ( b ) 、 G ( c ) 与 F ( b ) 建立拓扑关系, F ( a ) 、 F ( b ) 与 F 拓扑关系。
（2）断层线呈T型相交（图4）
该情况下断层抽象模型对象集合 F 及子地质体对象的集合 G 分别为式（8）-（9）。
F = { F ( a ) , F ( b ) } （8）
G = { G ( a ) , G ( b ) , G ( c ) } （9）
式中： G ( a ) 、 F ( b ) 与 F ( a ) 建立拓扑关系, G ( b ) 、 G ( c ) 与 F ( b ) 建立拓扑关系, F ( a ) 、 F ( b ) 与 F 建立拓扑关系。
（3）断层线呈X字相交（图5）
该情况下断层抽象模型对象集合 F 及子地质体对象的集合 G （图6）分别为式（10）-（11）。
F = { F ( a ) , F ( b ) , F ( a - 1 ) , F ( a - 2 ) , F ( b - 1 ) , F ( b - 2 ) } （10）
G = { G ( a ) , G ( b ) , G ( c ) } （11）
式中： G ( a ) 、 G ( b ) 与 F ( a - 1 ) , G ( c ) 、 G ( d ) 与 F ( a - 2 ) , F ( a - 1 ) 、 F ( a - 2 ) 与F（a）建立拓扑关系; G ( a ) 、 G ( c ) 与 F ( b - 1 ) , G ( b ) 、 G ( d ) 与 F ( b - 2 ) , F ( b - 1 ) 、 F ( b - 2 ) 与 F ( b ) 建立拓扑关系, F ( a ) 、 F ( b ) 与 F 建立拓扑关系。
Fig.3 Indication of no intersection

图3 不相交

Fig.4 T type intersection

图4 T型相交

Fig.5 X type intersection

图5 X型相交

Fig.6 The object set

图6 对象集合

4.2.4 生成断层表达式的算子集合
遍历断层数据表的每条记录,根据每条记录中的断层类型、旋转扩展、位移字段构造几何运算算子集合。其中,断层类型及旋转扩展字段决定所生成的算子类型,位移字段则决定算子数值属性。每条记录可能生成一个或多个算子,具体如表3所示。
Tab.3 Operator set of fault expression

表3 断层表达式方法算子集合

断层类型 算子集合(旋转扩展) 算子集合
正断层 d、r d
逆断层 u、r u
平移断层 t、r t
正-平移断层 d、t、r d、t
逆-平移断层 u、t、r u、t
在每次遍历断层数据表的条记录生成算子的同时,结合子地质体数据表中所属断层编号字段,即可将几何运算算子与子地质体运算对象建立映射关系。
4.2.5 建立算子集合与运算对象的映射关系
首先根据上述所建立的2类对象集合间及集合内部对象的拓扑关系生成拓扑运算及混合运算算子,并与建立的几何运算算子与子地质体对象的映射关系相结合,确定3种类型的算子作用的运算 对象。
Fig.7 Structure ofa F(a)

图7 断层F(a)结构图

Fig.8 Structure of F(b)

图8 断层F(b)结构图

以图4为例进行说明。假设断层 F ( a ) 、 F ( b ) 分别为正断层和逆断层,且2条断层的旋转扩展均为空。断层F(a)由断盘A、B构成（图7）,断层 F ( b ) 由断盘C、D构成（图8）。由断层 F ( a ) 、 F ( b ) 的类型可知,断盘B沿F(a)的断层面下降,断盘C沿 F ( b ) 的断层面上升。而断盘B即为 F ( b ) ,断盘C即为 G ( b ) 。由表2可构造出算子 d 、 u 。将算子与上述生成运算对象集合 F = { F ( a ) , F ( b ) } , G = { G ( a ) , G ( b ) , G ( c ) } 建立映射关系如式（12）所示。
d → F ( b ) u → G ( b ) （12）
根据图4情况所构建的拓扑关系,建立如式（13）所示的映射关系。
F ( b ) = dF ( b ) G ( b ) = uG ( b ) F = F ( a ) F ( a ) = G ( a ) pF ( b ) F ( b ) = G ( c ) pG ( b ) （13）
最终算子集合与对象集合建立式（14）所示的映射关系。
d → F ( b ) u → G ( b ) p → G ( a ) , F ( b ) p → G ( b ) , G ( c ) （14）
4.2.6 构建断层地质体抽象模型
根据断层数据表中的每条记录中的优先级（优先级决定的构建次序）字段将所建立的算子与映射关系集合按照断层表达式的文法规则进行重组,生成最终的断层表达式。以图4所代表的地质体 F 为例。假设 F ( a ) 优先级低于 F ( b ) ,图4所代表的断层表达式生成过程如式（15）-（20）所示。
F ( b ) = dF ( b ) （15）
G ( b ) = uG ( b ) （16）
F = F ( a ) （17）
F ( a ) = G ( a ) pF ( b ) （18）
F ( b ) = G ( c ) pG ( b ) （19）
最终断层表达式如式（20）所示。
F = G ( a ) pd ( G ( c ) puG ( b ) ) （20）
对断层表达式 F 进行计算即可生成断层地质体抽象模型,后续可在断层地质体抽象模型的基础上进行三维建模及可视化。
4.2.7 算法描述
构建断层表达式算法描述如下：
输入：断层描述表（F-tab）、断层线数据表（FL-tab）,子地质体数据表（G-tab）、采样点集数据表（Spt-tab）。
算法
{
构建初始六面体G;
遍历FL-tab,构建断层线集合Line,并记录数量n;
If（ n >1）
{
依据断层线所属断层优先级对Line排序;
While（ Line非空）
{
取Line中两条断层线L1、L2;
If（L1与L2相交）//这里的相交指交点不能同时位于L1与L2的端点
{
计算L1与L2交点P;
If（P非L1、L2端点）
{
求交生成子断层线L1-1、L1-2、L2-1、L2-2;
删除L1、L2,并在Line相应位置插入L1-1、L1-2、 L2-1、L2-2;
以P为依据划分Spt-tab;
将L1-1、L1-2、L2-1、L2-2插入FL-tab;
置L1、L2相应于F-tab中记录标志位无效;
}
Else
{ //假设P位于L1端点
求交生成子断层线L2-1、L2-2;
删除L2,并在Line相应位置插入L2-1、L2-2;
以P为依据划分Spt-tab;
将L2-1、L2-2插入FL-tab;
置L2相应于F-tab中标志位无效;
}
}
}
}
While（遍历FL-tab未完毕且相应于F-tab记录中标志位有效）
{
提取倾向、倾角,结合断层线构建断层面S;
将S插入断层面集合Face;
}
While（Face非空）
{ //G（a）为上盘,G（b）为下盘
取断层面S切分G,生成G（a）、G（b）;
将G（a）、G（b）插入G-tab;
建立G（a）、G（b）拓扑关系（即生成算子p）,生成F;
删除S;
}
While（遍历F-tab,标志位有效）
{ // 旋转扩展包括轴位置、轴向量、旋转角度
提取类型、扩展、旋转扩展字段;
以表2为依据生成算子集合O;
遍历G-tab,提取F-tab当前记录对应运算对象G（a）、G（b）;
建立O与G（a）、G（b）映射关系;
根据断层表达式方法规则及O与G（a）、G（b）映射关系构建断层表达式Exp;
}
计算Exp生成断层地质体抽象模型。
}

5 地质断层三维建模系统

本文以具有3条断层的达连河地质体为例,运用表达式规则实现了断层三维建模系统。该断层三维建模系统包括2个基本模块。
（1）构建抽象断层地质体模块
该模块完成原始断层数据的规范化、构建断层表达式并对表达式计算生成断层地质体抽象模型。断层地质体抽象模型的构建是一个迭代的过程,如果断层抽象模型不符合实际勘探数据,则可以对标准表格数据修改,进行迭代处理。
利用该模块对达连河断层地质体原始勘探数据规范为数据表,数据表部分字段部分数据示意如表4-8所示。
Tab.4 Data indication of the fault description table

表4 断层描述表数据示意

编号 类型 优先级 起始位置 终止位置 倾向 倾角 位移 旋转扩展 标志位
fault-A 正断层 3 528 822.6059,5 113 857.152,-75.26 528 498.1747,5 114 089.661,-4.5523 2.2528,
2.086,
-9.897		 72.44 -0.1877,
7.633,
-11.7466		 见表5 1
fault-B 正断层 2 528 700,5 113 772.154,-76.62 528 073.5928,5 113 713.364,-102.62 50.6899,18.393,
-11.71		 82.34 -6.144,
-5.716,
-20.00		 null 1
fault-C 正断层 1 526 804.3254,5 113 154.351,-103.89 526 993.2164,5 114 005.912,6.8673 7.1764,
-9.099,
-20.24		 68.23 17.7343,
-22.972,
-23.68		 null 1
Tab.5 Data indication of the rotation expansion of fault-A

表5 断层fault-A旋转扩展数据示意

轴位置 轴向量 旋转角度
528 498.1747,5 114 089.661,-4.5523 0.5882,0.6657,0.4592 22.34
Tab.6 Data indication of the fault line table

表6 断层线数据表数据示意

编号 所属断层编号 采样点集数据表指针 起始点 终止点
fl-01 falut-A sample-pt 528 824.8587,5 113 859.238,
-85.157		 528 498.275,5 114 096.104,
-14.469
fl-02 falut-B sample-pt 528 750.6899,5 113 790.547,
-88.330		 528 073.5928,5 113 713.364,
-113.330
fl-03 fault-C sample-pt 526 811.5018,5 113 144.252,
-124.130		 526 993.2164,5 114 005.912,
-5.627
Tab.7 Data indication of the fault point table

表7 断点数据表数据示意

编号 所属断层编号 坐标
pt-fl-01-001 fault-01 526 804.3254,5 113 153.351,-103.89
pt-fl-01-002 fault-01 526 868.7078,5 113 518.247, -58.93
pt-fl-01-003 fault-01 526 928.5775,5 113 855.739, -27.66
pt-fl-01-004 fault-01 526 947.4465,5 113 926.165, -20.85
pt-fl-01-005 fault-01 526 993.2164,5 114 005.912, 6.86
… … …
Tab.8 Data indication of the sample points table

表8 采样点集数据表（sample-pt）数据示意

编号 所属断层线编号 坐标
1 spt-fl-01-001 fl-01 526 951.535,5 113 925.42,-33.35
2 spt-fl-01-002 fl-01 526 935.531,5 113 849.70, -43.89
3 spt-fl-01-003 fl-01 526 875.784,5 113 510.98, -76.99
4 spt-fl-01-004 fl-02 528 737.342,5 113 784.59, -88.33
5 spt-fl-01-005 fl-02 528 700.000,5 113 772.15, -87.33
… … …
根据上述数据表构建断层表达式（式（21））。
F = G ( a ) pd ( G ( b ) pdG ( d ) p ( rdG ( c ) ) ) （21）
对上述断层表达（式（21））进行计算,即可得到达连河断层地质体抽象模型。
（2）三维建模及可视化模块
该模块完成对断层地质体抽象模型的可视化,包括拟合断层面、整理地层数据、划分地层数据为不同的区块、拟合地层层面、生成断层地质体块体模型、断层块体模型的三维可视化。
在拟合断层面及地层面时首先运用Kriging空间插值算法[6-8]对断层面控制点进行加密处理,然后运用Delaunay三角剖分算法[15-16]进行拟合。在进行地层层面拟合时,将地层控制点数据与对应的断点数据相结合,从而避免了地层层面与断层面的交割处理,降低了复杂度。运用文献[17]-[19]提出的三维地层建模方法结合每个地层顶、底板生成断层地质体块体模型。
本断层三维建模系统以Visual Studio 2012、Qt作为开发工具,结合Open GL图形库,在Window 7系统下对达连河地质体进行可视化工作。结合达连河研究区域的钻孔数据、断点数据、剖面数据、地层数据等构建出达连河地质体的可视化模型。图9为模型全貌,图10为仅包含断层地层的可视化模型。
Fig.9 The 3D model of Dalian river geologic body

图9 达连河地质体三维模型

Fig.10 The 3D layer model including the faults of DaLian river

图10 达连河断层地层三维可视化模型

从图10可看出,该地质体包含3条断层A、B、C。基中A断距较大,且断盘绕一端发生大角度的旋转（图11）,而断层B、C断距相对较小（图12）。经证明,运用表达式方法构建的达连河地区的三维可视化模型准确地展示了该地区的地层与构造展布特征,并在经过对所建模型数据的提取与分析,结果显示,地层面控制点与断层面控制点数据与真实情况相比,误差在可接受的范围内（误差≤1.5 m）。证明断层表达式建模技术对于地质断层三维建模及可视化应用的正确性。
Fig.11 The close view of fault A

图11 断层A近展示意图

Fig.12 The close view of fault B and C

图12 断层B、C近展示意图

经分析可得出断层表达式建模技术相对于引言所述方法优点如下：
（1）与整体建模技术相比,断层表达式建模技术能够较准确地模拟断距较大的断层（图11）,克服了整体建模技术适用于断距较小的断层局限性问题。
（2）与现有局部建模技术相比,断层表达式建模技术能够较准确地模拟断距较小的断层（图12）,克服了局部建模技术适用于断距较大断层的局限性问题;同时,断层表达式方法仅在对原始数据规范化过程中要求建模人员具有较高的地质学背景,后续步骤基于规范化的二维数据表自动进行断层表达式的构建、计算得到断层抽象模型以及基于断层抽象断层模型进行三建模及可视化,因而大大降低了局部法的复杂度,易于推广。
（3）与统一建模技术相比,断层表达式建模技术延用局部法建模技术将地质体分为不同区块,然后分别对各个区块进行三维建模的思路,避免了对地层断裂处的插值,从而解决了统一建模技术所面临的插值困境问题。另外,在进行地层层面拟合时,将地层数据与相应断点数据相结合,避免了进行地层层面与断层面的交割处理。在降低了建模复杂度的同时,达到较准确地模拟地层层面的目的。

6 结论

本文提出了运用表达式方法对地质断层进行建模。首先,以断层抽象模型（F）和子地质体（G）作为运算对象,以u、d、t、r、p、=作为算子,制定断层表达式的运算规则;然后,通过对原始数据进行规范化形成4个二维数据表,从二维数据表中提取算子和运算对象,形成断层表达式,并对断层表达式进行运算生成断层地质体抽象模型,并以此模型为基础进行三维建模及可视化。实验证明了该方法的有效性。
下一步工作将继续深入研究如何提高断层表达式方法的运算效率,以及如何增加表达式中运算对象和算子的信息,以便为地质体的储量分析等计算提供服务。

The authors have declared that no competing interests exist.

References

Publishing order | Descend order by publishing year | Descend order by cited within
[1]
李忠权,刘顺.构造地质学[M].北京:地质出版社,2010.

[ Li Z Q, Liu S.Structural geology[M]. Beijing: Geological Publishing, 2010. ]

[2]
武强,徐华.虚拟地质建模及可视化[M].北京:科学出版社,2011.

[ Wu Q, Xu H.3D Geological modeling and virtual visualization[M]. Beijing: China Science Publishing, 2010. ]

[3]
朱良峰,潘信,吴信才,等.地质断层三维可视化模型的构建与实现技术[J].软件学报,2008(8):2004-2015.

[ Zhu L F, Pan X, Wu X C, et al.Construction method and actualizing techniques of 3D visual model for geological faults[J]. Journal of Software, 2008,8:2004-1015. ]

[4]
蒋晓路,陈浩凡,鲁才,等.基于曲面约束的地质层面重构方法研究[J].计算机应用研究,2014,31(6):1894-1897.提出了基于曲面约束的地质层面 重构方法,其基本思路是先构建断面形成约束曲面,在地质层面和断层拓扑关系约束条件下,形成基于断层曲面的空间约束的层面插值问题,并在层面穿越约束断面 的部分进行矢量裁剪。通过仿真,该方法有效解决了复杂地质构造的层面的矢量重构问题;并通过实际的工区数据测试,有效支撑了复杂地质构造下的构造成图和构 造建模应用。由此,该方法回避了现有方法存在的层面区域分割所需边界信息难以获取的问题,并且实现了复杂地质构造的层面矢量重构。

DOI

[ Jiang X L, Chen H F, Lu C, et al.Reconstruction algorithm of stratum surface based on constraints of fault surfaces[J]. Application Research of Computers, 2014,31(6):1894-1897. ]

[5]
杨洋,潘懋,吴耕宇,等.一种新的轮廓线三维地质表面重建方法[J].地球信息科学学报,2015,17(3):253-259.<p>在三维地质建模中,由于地质剖面与医学CT有本质上的相似性,因此,轮廓线的三维表面重建技术有广泛的应用。目前,较为成熟的算法并没有考虑地质数据的特殊性,由于地质体形态复杂多变,地质剖面稀疏、数据来源多样等特点,将普适的算法应用到地质建模中,只是解决了重建曲面的合理性,并没有关注曲面的几何质量,导致曲面无法达到地质建模对于模型质量的要求。本文针对地质数据的特点,提出了一种新的轮廓线表面重构方法。在不改变轮廓线连接的合理性与正确性的情况下,通过生成过渡剖面轮廓线,克服传统方法生成曲面过于粗糙的缺点,以提高建模质量。在Creatar 三维地质基础平台上,对多种情况的剖面数据实验结果表明,在数据稀疏的情况下达到了良好的建模效果。该方法的扩展模型满足了地质剖面轮廓线建模中复杂地质情况,以及交叉剖面数据的曲面重建,且算法复杂度低,实现简单,有较强的实用性。</p>

DOI

[ Yang Y, Pan M, Wu G Y, et al.High quality geological surface reconstruction from planar contours[J]. Journal of Geo-information Science, 2015,17(3):253-259. ]

[6]
王英博,王栋,李钟学,等.基于SELA-Kriging算法的三维地质建模[J].系统工程理论与实践,2014,34(11):2913-2928.Kriging方法可对离散数据点进行内插和外推,近年来在地质、测绘各方面得到了普遍应用,但Kriging方法中变异函数理论模型的选择和实验变异函数参数的设定大部分依赖于地质工作人员的专业技能,具有一定的客观性. 引入混合蛙跳算法,对Kriging实验变异函数进行优化,提出SFLA-Kriging算法并结合某露天矿钻孔数据进行三维建模. 实验结果标明,SFLA-Kriging算法与传统Kriging方法相比,克服了一般方法的局限性,误差降低27.89%、三维地质模型精度更高.

[ Wang Y B, Wang D, Li Z X, et al.Research on 3D geological modeling based on shuffled frog leaping algorithm and Kriging[J]. Systems Engineering-Theory&Practices, 2014,34(11):2913-2928. ]

[7]
张芳,王军辉.普通Krging方法建立宏观地层有验证与评价[J].计算机工程与应用,2010,46(20):234-236.空间插值方法对空间变异性和空间相关性的反映程度及其准确性与精度直接影响到三维地质模型的真实性。普通Kriging方法应用于数字地层中，其变异函数模型和参数选择对结果影响很大，利用最值、方差、相关系数、平均误差以及误差分布统计为主要比较指标，采用交叉验证方法结合误差分析，对球状模型、指数模型、线性模型进行比较研究，初步得出线性有块金模型优于其他变异函数模型的结论，对数字地层的工程应用及其真实性评价有一定指导意义。

DOI

[ Zhang F, Wang J H.Variograms model analysis of Kriging in digital stratum[J]. Computer Engineering and Applications, 2010,46(20):234-236. ]

[8]
王靖波,潘懋,张绪定.基于Kriging方法的空间散乱点插值[J].计算机辅助设计与图形学报,1999,1(6):525-529.给出的一种距离加权插值方法借鉴地质统计学中的Ｋｒｉｇｉｎｇ方法，在空间数据场满足给定的统计分布特征的前提假设下进行插值。通过对帽状曲面函数的插值结果的观察，可以看出该方法可以获得良好的插值效果。

DOI

[ Wang J B, Pan M, Zhang X D.The Kringing interpolation method for scattered data points[J]. Journal of Computer Aided Design and Computer Graphics, 1999,1(6):525-529. ]

[9]
Robert W Sebesta.程序设计语言原理[M].张勤,王方矩,译.北京:机械工业出版社,2007:82-87.

[ Robert W Sebesta.Concepts of programing language[M]. Beijing: China Mechine Press, 2007:82-87. ]

[10]
Alfred V Aho, Monica S Lam, Ravi Sethi, et al.编译原理[M].北京:机械工业出版社,2008:125-131.

[ Alfred V Aho, Monica S Lam, Ravi Sethi, et al.Complies: principles, techniques and tools[M]. Beijing: China Machine Press, 2008:125-131. ]

[11]
朱志澄,曾佐勋,樊光明.构造地质学[M].北京:机械工业出版社,2008.

[ Zu Z C, Zheng Z X, Fan G M.Structural geology[M]. Beijing: China Machine Press, 2008. ]

[12]
朱良峰,潘信.地质断层三维构模技术研究[J].岩土力学,2008,29(1):274-275.对包含断层的复杂地质体三维构模技术展开讨论并提出相应的解决方案。分析了实际勘察所获取的 断层数据类型,设计了标准的断层数据库格式;系统地评价了两类基本的断层构模方法（即基于地层恢复的断层构模技术和基于分区插值的断层构模技术）,分析了 各自的特点和不足,在此基础上提出了断层与地层的统一构模技术;该技术既能应用于同沉积断层,又能适应断层终止于地层内部的情况,还能处理具有多值面的逆 断层构造的拟合问题,在理论上能够构造任意复杂的断层系统。通过北京奥运公园场址的断层构模实例验证了该方法的有效性。

DOI

[ Zhu L F, Pan X.Study on computer modeling of geological faults in 3D[M]. Rock and Soil Mechanics, 2008,29(1):274-275. ]

[13]
Wu Q, Xu H.An effective method for 3D geological modeling with multi source data integration[J]. Computers & Geosciences, 2005,31(1):35-43.The existing 3D geological modeling systems rely heavily on large numbers of borehole and cross-section data. However, it is well known that the available geological data are generally sparse and undersampled. In this paper, we propose a stepwise refinement method for 3D modeling with multi-source data integration. The method can naturally simulate geological structures no matter whether the available geological data are sufficient or not. By stepwise refinement on multiple data, the method increases the accuracy of 3D models gradually and effectively. In addition, the mechanisms used in the method for organizing and manipulating information can have an equally important impact upon geologists鈥 thought, the interpretation of geological data, and 3D modeling methodology. A concrete example of using the method to Huai Bei fault and fold belt shows that the method can be applied to broad and complex geological areas.

DOI

[14]
刘喜. 三维地质建模关键算法研究[D].北京:中国科学技术大学,2011.

[ Liu Xi.The research on key algorithm of three dimensional geological modeling[D]. Beijing: University of Science and Technology of China, 2011. ]

[15]
余杰,吕品,郑昌文.Delaunay 三角网构建方法比较研究[J].中国图形图象学报,2010,15(8):1158-1167.Delaunay三角网构建是3维场景可视化领域的一个热点也是难点问题。归纳总结了现有Delaunay三角网构建研究中的3类方法——逐点插入法、三角网生长法和分治法,以及在各自原理框架下的不同实现算法;比较分析了3种不同方法的优缺点和各自代表性算法的时间复杂度,并详细讨论了Delaunay三角网构建方法在大规模场景渲染和地形可视化领域中未来3个研究方向:混合算法研究、算法支撑技术研究和分布式并行算法研究。

[ Yu J, Lu P, Zheng C W.A comparative research on mehthods of delaunay triangulation[J]. Journal of Image and Graphics, 2010,15(8):1158-1167. ]

[16]
周雪梅,黎应飞.基于Bowyer-Watson三角网生成算法的研究[J].计算机工程与应用,2013,49(6):198-200.不规则三角网（Triangulated Irregular Network，TIN）在表示地形的形态方面具有较好的表现，其生成算法一直备受关注。讨论了三角网的数据结构的设计，采用逐点插入算法中的Bowyer-Watson算法思想为研究重点，设计并实现了该算法，对算法实验过程中可能出现的交叉现象进行分析，给出算法的改进。该改进算法已用于地形的可视化建模中，获得了较好的效果，对于三角剖分的相关研究具有一定的价值。

[ Zhou X M, Li Y F.Algorithm research to generate triangulation network based on Bowyer-Watson[J]. Computer Engineering and Applications, 2013,49(6):198-200. ]

[17]
Liu X, Dong L F, Lu D T.Three dimensional stratum modeling and a fast algorithm for arbitrary cutting based on half-edge structure[C]. IEEE International Conference on Computational Intelligence and Software Engineering (IEEE CiSE2009), Wuhan, China, 2009:1-4.

[18]
An C R, Liu Z, Bai Y L.Research on the key technologies of three-dimensional stratum modeling[C]. 2010 International Conference on Computer Application and System Modeling (ICCASM 2010), Taiyuan, China, 2010:354-357.

[19]
张煜,白世伟.一种基于三棱柱体元的三维地层建模方法及应用[J].中国图形图像学报,2001,6(3):285-290.传统的基于表面的数字模型在构建地层时会导致地层模型显示和分析的复杂化 ,例如在基于剖面的地层轮廓重建时 ,不同剖面之间轮廓线的分叉问题以及对应和拼接问题等 .为了有效地显示三维层状地层的结构 ,建立了一种基于三棱柱体体元的数字地层模型 ,同时给出了相关的剖分算法 .由于体元的引入 ,致使相应的剖面和切割算法得到了简化 ,并解决了工程地质中常见的层状地层模型任意剖面和任意切割的显示问题 ,也使工程中的多种操作可以在计算机上得到模拟 .同时 ,还列举了基于该模型的相关应用 .

DOI

[ Zhang Y, Bai S W.An approach of 3D stratum modeling based on tri-prism volume elements[J]. Journal of Image and graphics, 2001,6(3):285-290. ]

Options
Outlines

/

〈 〉