工程地质三维模型精细构建方法研究.pdf
《工程地质三维模型精细构建方法研究.pdf》由会员分享,可在线阅读,更多相关《工程地质三维模型精细构建方法研究.pdf(10页珍藏版)》请在咨信网上搜索。
1、上海国土资源 Shanghai Land&Resources34 2023Vol.44.3 doi:10.3969/j.issn.2095-1329.2023.03.006工程地质三维模型精细构建方法研究王 军1,2(1.上海市地质调查研究院,上海 200072;2.上海市地质资料信息专业技术服务平台,上海 200072)摘 要:上海三维城市地质信息系统有很好的地质建模基础,但是由于地质数据及其应用本身所具有的复杂性和不确定性,建立的三维地质模型与实际的地质条件之间还有不小的差距。本文结合上海地质条件和钻孔数据情况,对当前的钻孔自动建模方法流程进行了分解、对实现过程缺陷进行了剖析,从工程地质条
2、件分析、数据预处理方法、建模算法优化三个层面开展了建模流程的精细优化,并选择典型的中山公园场地进行建模验证和精度分析,通过将建模结果与人工生成的成果图进行精度评估,取得了较好的建模效果。关键词:工程地质;三维模型;精细建模中图分类号:P642.4;TP311.521 文献标志码:A 文章编号:2095-1329(2023)03-0034-09三维地质建模是一门综合运用现代空间信息理论来研究地质体及其内部物理、化学属性的信息处理、数据组织、空间建模与数字表达,并运用科学计算可视化技术对其进行真三维再现与交互的科学与技术。由于地质数据及其应用本身所具有的复杂性、不确定性、信息不完备等特点,使得三维
3、地质建模成为制约三维 GIS 研究及其应用的主要难点之一。而钻孔资料因其直观、准确、详细的特性在三维地层模拟中具有至关重要的意义,根据钻孔数据直接构建三维地层实体模型也一直是国内外研究的热点。多年来,上海开展的三维城市地质调查、地质资料信息服务集群化与产业化等工作,积累了大量的地质资料,尤其工程地质勘察的钻孔资料,入库钻孔数多达 84 万个,数据基础非常扎实,并且上海现有三维信息系统中有专门的“工程地质钻孔地质建模”功能模块,专业人员在软件中指定任意的空间范围,即可自动建立对应区域的三维模型,对建立的模型还可以进行任意切割、基坑开挖、隧道模拟等,实时建模,方便快捷。但是在实际使用过程中发现,该
4、功能只适应较简单的地质条件,稍微复杂的地质条件,生成的模型与专业人员设想的地层分布形态相差较远。而对于整个建模过程,专业人员又无法直接参与进去,虽然程序更新过多个版本,但该功能模块未见实质性改进。随着城市地下空间资源开发利用工作的深入开展,业务工作对三维地质模型提出了更高的精度要求。目前“三维城市地质信息系统”中三维地质建模结果的准确性,已成为限制三维地质模型进一步服务于城市地下空间开发利用的瓶颈。1 国内外研究现状自 从 20 世 纪 90 年 代 初 期 加 拿 大 的 Simon W.Houlding 提出三维地质模拟概念以来,地质信息的三维可视化逐渐受到普遍的重视,一些国家相继推出三维
5、地质工程项目,以推动三维地质模拟技术的发展。如:澳大利亚于 1990 年实施了目的为获取深部地球信息的三维研究计划;加拿大于1989年提出“垂直地质填图”的要求;英国地质调查局从 20 世纪 90 年代初开始实施“伦敦计算机化地下与地表地质”(LOCUS)项目等。相关领域的地质工作者也一直密切关注地质体三维建模技术的研究进展。1992 年,国际勘探地球物理学家协会和欧洲勘探地球物理学家协会成立了 SEG/EAEG 3D 建模委员会,开展了 3D 建模工程(SEM)。19961999 年分别在英国Leeds、新西兰 Otago、英国 Bristot、美国 Fredericsburg举行了四次地质
6、计算机会议,内容包括地质建模、模拟收稿日期:2023-04-25修回日期:2023-07-31作者简介:王军,1984 年生,男,硕士,工程师,主要从事自然资源信息化研究。电子邮箱:基金项目:上海市科委研发公共服务平台资助项目上海国土资源Shanghai Land&Resources 2023Vol.44.3 35和可视化。1997 年在巴赛罗那召开的国际数学地质会议上,Graeme B.C.等强调地质材料三维重建、建模及可视化的重要性。1997 年 9 月,一个由青年科学家组成的国际组织在德国召开了地球科学信息可视化研讨会,涉及的主要问题就是以地理为参照的空间和时间数据可视化建模。国内在三维
7、地质建模领域也逐步开始研究。中国地质调查局杨东来组织专家进行了三维地质建模方法的研究,并出版专著;中国地质大学刘修国教授带领的团队结合城市地质项目中的地下三维地层结构情况进行了研究,开发出 MapGIS-TDE 三维地质建模平台;北京航空航天大学杨钦教授带领的网格天地(GridWorld)团队在深部地质结构建模方面取得了丰硕的研究成果;北京大学地球科学学院潘懋教授带领的超维创想(Creatar)团队在基于交叉剖面的三维地质建模方面取得了一定的进展。2 钻孔自动建模流程分析作为一直以来的研究热点,三维地质建模方法很多。本研究是在上海大量工程地质勘察钻孔数据的基础上进行的,本文中后续谈到建模方法都
8、是特指“基于标准地层的简单地质体全自动建模”。2.1 建模基本假定“基于标准地层的简单地质体全自动建模”方法适用于“简单层状地质体”,即连续层状地质体,是指未经变形或轻微变形的、形态相对简单、连续的层状地质实体,它没有被断层或节理切割,层与层之间为整合接触,在每一层内没有垂直方向属性的变化,它们具有水平或倾斜的状态,形态较为规则。根据大量的实例分析,归纳出简单层状地质体具有以下基本的几何特征:简单层状地质体可以按照一定的标准划分成多个地层;各个地层在研究区域内是连续分布的(如果因地层尖灭出现不连续,则可按照一定的规则作特殊处理);各个地层由上、下两界面封闭而成,上、下界面在水平面上的投影完全重
9、合;在研究区域内,任一平面点 P(x,y)在上层界面或下层界面上对应的高程值是唯一的,即具有单值性;在研究区域内,可按照一定的规则或标准整理出一套标准的地层分层顺序。2.2 建模流程分解整个自动建模,采用“主 TIN”(Primary TIN)的思想,即所有层面的三角网在平面上的投影都是同一模板,如此可沿地层深度自上而下推延至建模区域的全部地层。基于各个地层的建模数据,通过插值算法,结合一些尖灭、缺失信息,即可构建出各个地层面。地层面之间若有交错,则可根据交错关系,进行相交处理,将所有交点更新至“主TIN”模板中,通过如图 1 所示的一系列处理步骤来完成地层实体构建。2.3 实现过程缺陷分析通
10、过钻孔全自动建立三维模型的方法,不仅大大降低了建模的繁琐程度、节省了建模时间,又兼顾了模型的准确性,可大大提高建模效率。但经实际验证,这一建模方法在具体实现时有如下几点缺陷:(1)地层层序约束过于严格要求在建模区内有一套完整的地层层序表,并要求钻孔分层数据严格按照这一地层层序解释,不能出现地层倒转、互层、透镜体等情况,否则程序无法正常完成建模过程。(2)地层级别处理简化在地质专业人员的常识中,地层按照年代成因是存在大层、亚层、次亚层等级别关系的,如亚层的地层分界面是不会穿越大层的。而建模算法中虽然也有对地层级别的处理,但仅限于处理钻孔地层是否需要合并,建模过程中层与层之间都是对等的。(3)地层
11、缺失区边界计算粗放在早期版本的程序中,出现地层尖灭后,程序是根据图 1 根据钻孔数据构建三维地层模型流程图Fig.1 Flowchart of constructing 3D stratigraphic model based on borehole data上海国土资源 Shanghai Land&Resources36 2023Vol.44.3 用户指定的尖灭位置(如全局统一 1/2 位置等),但是这种尖灭方式与孔的相对位置直接相关,与实际情况不符;后来的版本也尝试了多种办法自动计算尖灭位置,但始终效果不理想(如图 2 a)。(4)地层“0 厚度”的假定不符合地质思维在上一步的缺失区边界计
12、算,实际上是根据各个层位的高程点调整出来的,为了在程序实现过程中严格满足地层层序,引入了“0 厚度”地层(如图 2 b),虽然建模结果看不到该地层,但该假定会较大程度的影响建模精度。例如工程地质层中的层中出现 a 透镜体,原系统中建模如要体现 a,需要在场地分层中将分为1、a、2,钻孔中的层也需要进行人为拆分,以保证钻孔地层的顺序满足标准地层层序。(5)对钻孔最后一层的处理错误钻孔分层数据中,最后一层虽然也指定了层底深度,但实际并不代表该层的实际底界位置,仅表示本层到该深层位置还没钻穿,实际底界位置一定大于等于该深度。3 工程地质三维模型精细构建根据前面对现有三维系统中建模算法的分解可知,建模
13、效果不理想,很大一部分原因在于程序在实现建模算法过程中,单纯从数学方法、DEM 层位叠加来考虑,尤其是对数据处理的过程,少有考虑地质现象本身沉积演化规律导致。本研究是在现有三维系统中“基于标准地层的简单地质体全自动建模”的基础上,结合上海工程地质条件和钻孔数据特点,针对原有建模缺陷进行优化。优化的基本原则是:以钻孔资料作为构建地层模型的源数据,并尽可能多的将地质沉积演化思维融合到自动建模过程中,然后进行相应的建模工作,提升所建模型的准确性。3.1 数据预处理优化建模过程是以用户数据为基础的,建模前对数据的预处理是将用户数据转换为建模使用数据的过程,数据预处理的方式方法、精细度会直接影响到最终的
14、建模效果。3.1.1 标准层整理标准地层是一个自上而下的一个有严格层次关系的地层级别树。钻孔、地质剖面或者地层顶底板等值线数据中的地层可以引用标准地层中的地层信息。对于场地范围,作为简单层状体地质体,地层的上下关系、层次关系会严格遵循标准地层的整个要求,根据钻孔数据自动建模方法正是以此为基础进行的。本文在对标准地层的处理方法上,数据模型主要有两点变化:(1)强化地层级别概念参照“三维系统”中剖面图自动成图的原理,将原大层、亚层、次亚层的层次关系作为强约束关系来限定建模过程,后续的建模过程严格遵循“先连大层再连小层”的原则,保证层位的逻辑关系合理。(2)增加对透镜体的支持若研究区范围内部分地层内
15、含类透镜体对象,程序支持在标准地层中将该层进行标记,不作为亚层处理,待其它层面都完成后再单独处理,避免了原“三维系统”中需要地质专业人员将夹层人为分为上下两层的困惑。3.1.2 钻孔分层细分钻孔分层数据是工程地质钻孔自动建模最主要的数据源,对该部分数据的精细处理是保证建模精度非常重要的环节。结合整个建模过程,本研究主要对钻孔分层出现(或缺失)进行了如下细分:(1)钻孔中完整出现的地层钻孔中完整出现的地层是指在整个钻孔分层数据中,该层有出现,且顶板、底板高程确定的地层。(2)钻孔中缺失的地层钻孔中缺失的地层指按照场地标准地层层序,该层位下部的地层已经出现,该地层还未出现的情况,如钻孔中出现了、,
16、按照标准地层层序,本该在上部出现而未出现,即判为缺失。(3)钻孔中未钻穿的地层钻孔中未钻穿的地层是指在钻孔最后一层,虽然该层在数据库中也有明确的层底深度数据,但实际并不代表该层的实际底界位置,仅表示本层到该深度还未钻穿,实际底界位置一定大于等于该深度。(4)钻孔中未揭露的地层aa112a212b图 2 地层尖灭及透镜体处理Fig.2 Formation pinch-out and lens processing上海国土资源Shanghai Land&Resources 2023Vol.44.3 37钻孔中未揭露的地层指按照场地标准地层层序,钻孔最后一层下部的地层均为未揭露地层。3.2 建模算法
17、优化在建模数据处理优化的基础上,对原建模算法各个环节进行了针对性的优化,力争做到在不影响建模效率情况下更准确的展现地下地质体的形态。3.2.1 地层分级约束在原建模流程中,整个过程是按照标准地层层序,从上至下,逐层构建地质面,最终完成整个地质体的构建。本研究中,因引入“先连大层再连小层”的地层级别思维,将原建模过程的大循环改为递归,具体步骤如下(图3 所示):(1)初始建模时,首先将所有钻孔的地层合并到大层,并按照标准地层层序,循环构建大层的各个底部地层面,并完成后续的地层面高层调整,上下面求交得到相应地层缺失区边界;(2)后续按照标准地层层序,确立当前建模的地层级别和地层层号,将所有钻孔的分
18、层经过整理,取到与当前级别、当前层号匹配的逻辑钻孔数据(每次传入的并不是整个钻孔,而是钻孔的某一段分层数据,并且分层是根据当前地层级别合并过的),循环构建完成该级别的所有地层面,并且构建的地层面高程也受父地层的顶板、底板高程来约束;(3)递归完成子层的地层面构建,直到达到建模需要的地层级别或没有子层为止。3.2.2 沉积过程反演在整个地质建模过程中,对于相邻钻孔之间地层对比一致的情况,直接连接地层底部即可得到地层面,而最重要、最需要精细处理的部分,即为相邻钻孔地层对比不一致的情况,通过对尖灭点的计算得到地层缺失区边界。在前文建模流程分析部分,对现有程序模块的计算过程进行了详细分解,并总结了数据
19、处理过程中的缺陷。结合实际问题,在本次程序编写过程中也着重对该部分进行了反复尝试和验证,具体理解和实践过程如下:(1)上海地区工程地质条件总体变化不大,工程建设影响范围内(埋深 100 m 以浅)主要为晚更新世以来海陆交替相、以海相为主的松散堆积物;(2)沉积地层是在长久的沉积作用下,碎屑沉积物不断堆积形成的,地层缺失部分区域则主要是由于不断的海进海退、河流侵蚀(淤积)等地质过程的作用导致;(3)某地层缺失区的边界计算基本原理:相邻的两个钻孔,一个钻孔存在该地层,另一个钻孔缺失该地层,那么该地层的缺失区边界一定在两个钻孔之间的某个位置;(4)地层尖灭点为位置,不仅与两个钻孔的相对位置相关,还与
20、两个钻孔中上下地层底板高程、地层厚度等相关,不能简单的按照钻孔间距直接计算;通过与地质专业人员的交流,将出现地层局部缺失的原因概化为如下两种情况:(1)沉积地层被河流侵蚀:沉积地层在海退暴露地表后,受到地表形成的河流影响,河流所在区域的地层被不断侵蚀,导致厚度变薄甚至完全缺失;(2)侵蚀区填充地层被侵蚀:暴露在地表的地层因地质过程作用,再一次回到沉积作用为主的环境中,地层被侵蚀区域被新的碎屑沉积物填充,形成新的地层。以上两种情况只是现象说明,实际的地质过程一定是极为复杂的,海进海退及沉积间断交替进行的。综合前文分析情况,通过不断的总结尝试,程序中最终也按以上思维方式来计算地层尖灭点的具体位置,
21、具体如图 4 所示。图 4 中层为典型沉积地层被河流侵蚀(层被侵蚀后,层沉积填充),在程序处理过程中,首先将缺失层的 B 钻孔标记为不参与插值,通过其他钻孔层的底板高程,插值得出 B 钻孔所在位置的高程值,可视为被侵蚀前的当前位置层的层底高程,将这三点连起来与层界线的交点即为地层的尖灭点。这种计算方法看似简单,实则与以往直接平推或者 1/2 尖灭有本质区别,121312213424ABCD121312213424ABCD图 3 原不分级循环(A、C)与分级别递归(B、D)的 效果对比Fig.3 Comparison of the effects of the original ungraded
22、 loop(A,C)and hierarchical recursion(B,D)上海国土资源 Shanghai Land&Resources38 2023Vol.44.3 该计算方法计算出的尖灭点位置与相邻钻孔的相对位置、相邻钻孔中地层厚度以及全区地层的整体走势都相关。图 5 中层为典型侵蚀区填充地层被侵蚀(层被侵蚀后,层沉积填充后被侵蚀,在钻孔B处层完全消失),在程序处理过程中,首先将缺失层的 B 钻孔标记为不参与插值,通过其他钻孔层的底板高程,插值得出 B钻孔所在位置的厚度值,可视为被侵蚀前当前位置层的厚度。该厚度越大,则说明上部地层沉积过程中对层的侵蚀强度越大。本研究中根据该厚度值反算
23、出的高程值,作为当前位置层的层底高程,将这三点连起来与层界线的交点即为地层的尖灭点。以上两种情况都是通过厚度反演的逻辑得到虚拟控制点,然后通过二次插值的方式与上下地层面求交计算出实际尖灭位置的,一致的逻辑也助于程序的具体实现。3.2.3 地层最小厚度根据与地质专业人员交流,在实际根据原始数据对钻孔进行分层的过程中,待划分的目标土层厚度小于 0.5 m时,一般是忽略掉该层,与上(下)层进行合并处理的。在本研究编制的程序中,也引入了该判断逻辑,可避免模型出现一些超薄地层导致建模结果偏离实际。3.2.4 未钻穿地层的处理对于钻孔分层数据中最后一层的层底深度数据,实际底界位置一定大于等于该深度。在实际
24、建模过程中,如图 6,首先将未钻穿层的钻孔 B 标记不参与插值,通过其他钻孔层的底板高程,插值得出 B 钻孔所在位置的深度值,若得到的高程值大于 B 孔层深度,则取该深度值为层的层底深度,否则取孔本身的层深度。3.2.5 未揭露地层的处理对于某钻孔分层数据中未揭露某层,处理方式与未钻穿地层相似,但计算出的高程受上部地层的底板高程约束。如图 6 中,B 孔中未揭露的层,计算层层位是B 孔不参与计算,但计算结果层的底界线不得超过层底的计算值。3.3 插值算法选择由于钻孔之间的距离稀疏程度、方向、数据值存在差异,钻孔以外未知的地质特性需要插值和推断,散乱点插值在地学领域有着广泛的应用前景,常用的空间
25、插值方法有以下几种:(1)距离倒数加权法距离倒数加权法用距离的倒数作为权值,算法简单,ABCABCABC图 5 侵蚀区填充层被侵蚀的尖灭点计算方法Fig.5 Calculation method of pinch-out point of eroded filling layer in erosion areaABCABC图 6 对钻孔未钻穿层的处理Fig.6 Treatment of undrilled layersABC ABC图 4 沉积地层被河流侵蚀的尖灭点计算方法Fig.4 Calculation method of pinch-out point of sedimentary st
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 工程地质 三维 模型 精细 构建 方法 研究
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【自信****多点】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【自信****多点】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。