CN101281557A - 一种应用于复杂地层数字制图的筛分方法 - Google Patents
一种应用于复杂地层数字制图的筛分方法 Download PDFInfo
- Publication number
- CN101281557A CN101281557A CNA2008100281512A CN200810028151A CN101281557A CN 101281557 A CN101281557 A CN 101281557A CN A2008100281512 A CNA2008100281512 A CN A2008100281512A CN 200810028151 A CN200810028151 A CN 200810028151A CN 101281557 A CN101281557 A CN 101281557A
- Authority
- CN
- China
- Prior art keywords
- layer
- boring
- stratum
- line
- type
- 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.)
- Granted
Links
Images
Landscapes
- Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
Abstract
本发明提供一种应用于复杂地层数字制图的筛分方法,通过步骤实现:(1)进行复杂地层分类,以钻孔地层为基础,将复杂地质情况在制图方法上分为基本情形,各种复杂地质情况都包括在其中,或由其组合衍生;(2)对钻孔及每个钻孔的地层进行绘制剖面图所需的信息预处理;(3)对复杂地层进行地层连线,构造不同的地层区域。本发明的优点是能快速有效地实现地层剖面自动化绘制,对正常地质状况和缺失、倒转、尖灭等复杂地质状况均有较好的实现效果。
Description
技术领域
本发明属于利用计算机程序绘制复杂地质情况下地层剖面的数字制图领域,特别是涉及一种应用于复杂地层数字制图的筛分方法。
技术背景
目前通过数字制图方法可绘制正常地质情况和地层缺失等部分复杂地质情况的地层剖面图,但复杂地质情况下地层剖面的数字制图尚缺乏统一有效方法,不能实现所有缺失、倒转、尖灭等地质情况的剖面图数字绘制。
发明内容
本发明的目的在于提供一种复杂地质情况下有效绘制地层剖面的数字制图方法,以弥补现有数字制图技术在地层剖面图绘制中的不足。
为了实现上述发明目的,采用的技术方案如下:
一种应用于复杂地层数字制图的筛分方法,通过如下步骤实现:
(1)进行复杂地层分类,以钻孔地层为基础,将复杂地质情况在制图方法上分为常用的几种基本情形,各种复杂地质情况都包括在其中,或由其组合衍生;
(2)对钻孔及每个钻孔的地层进行绘制剖面图所需的信息预处理;
(3)对复杂地层进行地层连线,构造不同的地层区域。
上述技术方案中,所述步骤(1)将复杂地质情况分为六类基本情形,具体如下:
假定有两钻孔D1与D2,且存在编号分别为1、2、3、4的三种地层:
类型A:D1中层3、2倒转,且D2中存在层2;
类型B:D1中层3、2倒转,但D2中无此两层;
类型C:D1中层3、2倒转,且D2中倒转情况相同;
类型D:D1中层3、2倒转,且在此倒转之上还有层2,D2中缺失层3,形成夹层;
类型E:D1中层2、1倒转,并与下层的层1构成夹层之势,D2中存在层2;
类型F:D1中层3、2倒转,但D2中层2、3依正常层序出现。
所述步骤(2)的信息预处理具体包括如下:
(21)设置一个算法数组,并确定剖面图最左侧的钻孔为算法数组起始钻孔,将其它钻孔按与起始钻孔的距离从小到大排序;
(22)将地层主层号与亚层号相结合作为地层层序,即层序=主层号*10+亚层号,并将地层信息按层深从小到大排序。
所述步骤(3)的地层连线具体实现如下:
(31)从左到右,从第一个未处理的钻孔起选定相邻的两个钻LK1、LK2;
(32)分别选定当前钻孔LK1、LK2中未处理的层为当前层;
(33)比较两个当前层层序,选择层序较小的层记为层A,另一个钻孔的当前层记为层B;
(34)判断层A与层B层序是否相同,不同转步骤(35),相同则转步骤(36);
(35)判断层A与层B+1是否相同,不同则将层A与层B的起点连线造三角形区域,并将层A+1和层B作为新当前层;若层A与层B+1相同则将层A的起点与层B连线造三角形区域,并将层A与层B+1作为新当前层,转步骤(37);
(36)判断层A与层A+2是否相同,如相同则将层A和层A+2与层B连线造倒“V”型区域,同时层A+1夹于倒“V”型区域中造三角形区域,并将层A+3和B+1作为新当前层;如不同则将层A和层B连线造四方形区域,并将层A+1和B+1作为新当前层。转步骤(37);
(37)转步骤(33),直到处理完当前钻孔所有的地层,若某钻孔已到达孔底而另一钻孔仍有未连线地层,则直接利用此孔底与剩余地层构造三角形区域。
本发明还包括步骤(4)如下:
(4)地层连线完成后,在造好的区域内填充相应的地层图例,并对所有数据进行整合即完成剖面图的绘制。
本发明的优点是能快速有效地实现地层剖面自动化绘制,对正常地质状况和缺失、倒转、尖灭等复杂地质状况均有较好的实现效果。
附图说明
附图1为复杂地层六种倒转情况。
附图2为连线造区方案。
附图3为复杂地层连线流程图。
附图4为实施例连线例子过程图。
具体实施方式
下面结合附图对本发明做进一步的说明。
本发明包括如下步骤:
1、进行复杂地层分类,以钻孔地层为基础,将复杂地质情况在制图方法上分为基本情形,各种复杂地质情况都包括在其中,或由其组合衍生;
本发明以钻孔地层为基础,将复杂地质情况在制图方法上分为附图1所示六类基本情形,其它各种复杂地质情况可由其组合衍生。假定有两钻孔D1与D2,且存在编号分别为1、2、3、4的三种地层:
(1)类型A:如图1(a),D1中层3、2倒转,且D2中存在层2;
(2)类型B:如图1(b),D1中层3、2倒转,但D2中无此两层;
(3)类型C:如图1(c),D1中层3、2倒转,且D2中倒转情况相同;
(4)类型D:如图1(d),D1中层3、2倒转,且在此倒转之上还有层2,D2中缺失层3,形成夹层;
(5)类型E:如图1(e),D1中层2、1倒转,并与下层的层1构成夹层之势,D2中存在层2;
(6)类型F:如图1(f),D1中层3、2倒转,但D2中层2、3依正常层序出现。
(2)对钻孔及每个钻孔的地层进行绘制剖面图所需的信息预处理;
(3)对复杂地层进行地层连线,构造不同的地层区域。
2、信息预处理。对绘制剖面图所需的钻孔及每个钻孔的地层进行信息预处理:
(1)确定剖面图最左侧的钻孔为算法数组起始钻孔,将其它钻孔按与起始钻孔的距离从小到大排序;
(2)地层信息一般包含层深、主层号、次层号三个组成部分,主层号表示大的地质层,次层号表示大的地质层下细分的子层。主层号、次层号的编号都是从1开始,次层号为个位数。层序是反映地层信息的序列,为了层序比较的方便,将地层主层号与亚层号相结合作为地层层序,即层序=主层号*10+亚层号,并将地层信息按层深从小到大排序。
3、复杂地层连线有效筛分,如附图3所示,在对复杂地层进行上述分类与预处理的基础上,利用下述有效筛选算法进行地层连线:
(1)从左到右,从第一个未处理的钻孔起选定相邻的两个钻LK1、LK2;
(2)分别选定当前钻孔LK1、LK2中未处理的层为当前层;
(3)比较两个当前层层序,选择层序较小的层记为层A,另一个钻孔的当前层记为层B;
(4)判断层A与层B层序是否相同,不同转步骤(5),同则转步骤(6);
(5)判断层A与层B+1是否相同,不同则如图2(a)连线造区,并将层A+1和层B作为新当前层;若层A与层B+1相同则进行如图2(b)连线造区,并将层A与层B+1作为新当前层,转步骤(7);
(6)判断层A与层A+2是否相同,如相同则如图2(c)连线造区,并将层A+3和B+1作为新当前层;如不同则如图2(d)连线造区,并将层A+1和B+1作为新当前层。转步骤(7);
(7)转步骤(3),直到处理完当前钻孔所有的地层。若某钻孔已到达孔底而另一钻孔仍有未连线地层,则直接利用此孔底与剩余地层构造三角区。
实施例一:
本实施例利用图4所示两钻孔进行地层连线和填充,钻孔LK1的地层顺序为13242,钻孔LK2的地层顺序为12324。利用本筛选方法连线造区的过程为:
(1)当前层A为LK1-1,B为LK2-1,二者相同;转步骤(6),A+2序号为2与A不同,因此按照图2(d)连线造区a,并更新A为LK1-2,B为LK2-2;
(2)A与B不同,转步骤(5);A与B+1的序号都为3,因此按照图2(b)连线造区b,并更新A为LK1-2,B为LK2-3;
(3)A与B相同,转步骤(6);A+2序号为4与A不同,因此按照图2(d)连线造区c,并更新A为LK1-3,B为LK2-4;
(4)A与B相同,转步骤(6);A+2序号为2与A相同,因此按照图2(c)连线造区d和e。此时A+3已为钻孔孔底,可直接与钻孔LK2剩余层构造三角区。此示例中只构造一个三角区外可。
实施例二:
本实施例将复杂地层数字制图的筛分方法应用在东改工程地下环境信息系统(UEEIS)和广东梅河高速公路工程高速公路空间地理信息可视化管理系统(HSIVS)的研制中。应用过程如下:
1、根据工程单位提供的钻孔柱状图,建立钻孔信息数据库;
2、根据用户在系统中指定的连线钻孔,从钻孔信息数据库中,提取相应的钻孔坐标数据和土层信息,并对各钻孔的地层信息进行预处理:
(1)确定剖面图最左侧的钻孔为算法数组起始钻孔,将其它钻孔按与起始钻孔的距离从小到大排序;
(2)将地层主层号与亚层号相结合作为地层层序,即层序=主层号*10+亚层号,并将地层信息按层深从小到大排序。
3、在上述钻孔信息预处理的基础上,采用地层连线有效筛选算法进行连线,具体步骤如下:
(1)从左到右,从第一个未处理的钻孔起选定相邻的两个钻孔ZK1、ZK2;
(2)分别选定当前钻孔ZK1、ZK2中未处理的层为当前层;
(3)比较两个当前层层序,选择层序较小的层记为层A,另一个钻孔的当前层记为层B;
(4)判断层A与层B层序是否相同,不同转步骤(5),同则转步骤(6);
(5)判断层A与层B+1是否相同,不同则如图2(a)连线造区,并将层A+1和层B作为新当前层;若层A与层B+1相同则进行如图2(b)连线造区,并将层A与层B+1作为新当前层,转步骤(7);
(6)判断层A与层A+2是否相同,如相同则如图2(c)连线造区,并将层A+3和B+1作为新当前层;如不同则如图2(d)连线造区,并将层A+1和B+1作为新当前层。转步骤(7);
(7)转步骤(3),直到处理完当前钻孔所有的地层。若某钻孔已到达孔底而另一钻孔仍有未连线地层,则直接利用此孔底与剩余地层构造三角区;
(8)按步骤(2)-(6)处理完所有钻孔。
4、地层连线完成后,在造好的区域内,根据各区域地层层号填充相应的地层图例,并结合钻孔坐标数据,在适当位置绘制标尺、文字等附属信息,通过对所有图形数据的叠加整合,生成可供用户浏览、编辑的地层剖面图。
Claims (5)
1、一种应用于复杂地层数字制图的筛分方法,其特征在于通过如下步骤实现:
(1)进行复杂地层分类,以钻孔地层为基础,将复杂地质情况在制图方法上分为基本情形,各种复杂地质情况都包括在其中,或由其组合衍生;
(2)对钻孔及每个钻孔的地层进行绘制剖面图所需的信息预处理;
(3)对复杂地层进行地层连线,构造不同的地层区域。
2、根据权利要求1所述的应用于复杂地层数字制图的筛分方法,其特征在于所述步骤(1)将复杂地质情况分为六类基本情形,具体如下:
假定有两钻孔D1与D2,且存在编号分别为1、2、3、4的三种地层:
类型A:D1中层3、2倒转,且D2中存在层2;
类型B:D1中层3、2倒转,但D2中无此两层;
类型C:D1中层3、2倒转,且D2中倒转情况相同;
类型D:D1中层3、2倒转,且在此倒转之上还有层2,D2中缺失层3,形成夹层;
类型E:D1中层2、1倒转,并与下层的层1构成夹层之势,D2中存在层2;
类型F:D1中层3、2倒转,但D2中层2、3依正常层序出现。
3、根据权利要求1或2所述的应用于复杂地层数字制图的筛分方法,其特征在于所述步骤(2)的信息预处理具体包括如下:
(21)设置一个算法数组,并确定剖面图最左侧的钻孔为算法数组起始钻孔,将其它钻孔按与起始钻孔的距离从小到大排序;
(22)将地层主层号与亚层号相结合作为地层层序,即层序=主层号*10+亚层号,并将地层信息按层深从小到大排序。
4、根据权利要求3所述的应用于复杂地层数字制图的筛分方法,其特征在于所述步骤(3)的地层连线具体实现如下:
(31)从左到右,从第一个未处理的钻孔起选定相邻的两个钻LK1、LK2;
(32)分别选定当前钻孔LK1、LK2中未处理的层为当前层;
(33)比较两个当前层层序,选择层序较小的层记为层A,另一个钻孔的当前层记为层B;
(34)判断层A与层B层序是否相同,不同转步骤(35),相同则转步骤(36);
(35)判断层A与层B+1是否相同,不同则将层A与层B的起点连线造三角形区域,并将层A+1和层B作为新当前层;若层A与层B+1相同则将层A的起点与层B连线造三角形区域,并将层A与层B+1作为新当前层,转步骤(37);
(36)判断层A与层A+2是否相同,如相同则将层A和层A+2与层B连线造倒“V”型区域,同时层A+1夹于倒“V”型区域中造三角形区域,并将层A+3和B+1作为新当前层;如不同则将层A和层B连线造四方形区域,并将层A+1和B+1作为新当前层。转步骤(37);
(37)转步骤(33),直到处理完当前钻孔所有的地层,若某钻孔已到达孔底而另一钻孔仍有未连线地层,则直接利用此孔底与剩余地层构造三角形区域。
5、根据权利要求4所述的应用于复杂地层数字制图的筛分方法,其特征在于还包括步骤(4)如下:
(4)地层连线完成后,在造好的区域内填充相应的地层图例,并对所有数据进行整合即完成剖面图的绘制。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2008100281512A CN101281557B (zh) | 2008-05-16 | 2008-05-16 | 一种应用于复杂地层数字制图的筛分方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2008100281512A CN101281557B (zh) | 2008-05-16 | 2008-05-16 | 一种应用于复杂地层数字制图的筛分方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101281557A true CN101281557A (zh) | 2008-10-08 |
CN101281557B CN101281557B (zh) | 2011-07-06 |
Family
ID=40014026
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2008100281512A Expired - Fee Related CN101281557B (zh) | 2008-05-16 | 2008-05-16 | 一种应用于复杂地层数字制图的筛分方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101281557B (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102354322A (zh) * | 2011-09-19 | 2012-02-15 | 北京城建勘测设计研究院有限责任公司 | 关于工程勘察剖面图内业整理的方法及装置 |
CN102646141A (zh) * | 2012-02-17 | 2012-08-22 | 上海同岩土木工程科技有限公司 | 一种非等深钻孔的地质剖面图自动成图方法 |
CN106504328A (zh) * | 2016-10-27 | 2017-03-15 | 电子科技大学 | 一种基于稀疏点云曲面重构的复杂地质构造建模方法 |
CN106709988A (zh) * | 2015-11-16 | 2017-05-24 | 天津市勘察院 | 一种工程地质剖面图构建方法 |
CN112908145A (zh) * | 2019-11-19 | 2021-06-04 | 核工业二一六大队 | 一种砂岩型铀矿砂体等厚度图成图方法 |
CN114048542A (zh) * | 2022-01-14 | 2022-02-15 | 四川省交通勘察设计研究院有限公司 | 一种公路工程数据空间的可视化方法及系统 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1151021A (zh) * | 1995-11-22 | 1997-06-04 | 中国海洋石油总公司海洋石油勘探开发研究中心 | 确定、展示和使用地下岩石弹性模量和密度相对变化的方法 |
-
2008
- 2008-05-16 CN CN2008100281512A patent/CN101281557B/zh not_active Expired - Fee Related
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102354322A (zh) * | 2011-09-19 | 2012-02-15 | 北京城建勘测设计研究院有限责任公司 | 关于工程勘察剖面图内业整理的方法及装置 |
CN102646141A (zh) * | 2012-02-17 | 2012-08-22 | 上海同岩土木工程科技有限公司 | 一种非等深钻孔的地质剖面图自动成图方法 |
CN102646141B (zh) * | 2012-02-17 | 2015-06-24 | 上海同岩土木工程科技有限公司 | 一种非等深钻孔的地质剖面图自动成图方法 |
CN106709988A (zh) * | 2015-11-16 | 2017-05-24 | 天津市勘察院 | 一种工程地质剖面图构建方法 |
CN106504328A (zh) * | 2016-10-27 | 2017-03-15 | 电子科技大学 | 一种基于稀疏点云曲面重构的复杂地质构造建模方法 |
CN112908145A (zh) * | 2019-11-19 | 2021-06-04 | 核工业二一六大队 | 一种砂岩型铀矿砂体等厚度图成图方法 |
CN114048542A (zh) * | 2022-01-14 | 2022-02-15 | 四川省交通勘察设计研究院有限公司 | 一种公路工程数据空间的可视化方法及系统 |
CN114048542B (zh) * | 2022-01-14 | 2022-07-15 | 四川省交通勘察设计研究院有限公司 | 一种公路工程数据空间的可视化方法及系统 |
Also Published As
Publication number | Publication date |
---|---|
CN101281557B (zh) | 2011-07-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101281557B (zh) | 一种应用于复杂地层数字制图的筛分方法 | |
CN103513293B (zh) | 隧道地质综合超前预报专家系统及其实现方法 | |
CN102646141B (zh) | 一种非等深钻孔的地质剖面图自动成图方法 | |
CN106547941A (zh) | 一种基于revit体量模型构建BIM工程勘察信息系统的方法 | |
CN104732359A (zh) | 油田地理信息与勘探开发协同工作平台系统 | |
CN104834518A (zh) | 油田勘探开发工作平台系统 | |
CN101617101A (zh) | 井和泄流位置的自动化油田开发方案 | |
CN106023313A (zh) | 一种三维地质数字勘察报告生成方法 | |
CN110211231B (zh) | 一种三维地质灾害信息模型建模方法 | |
Abrar et al. | AHP based analysis of groundwater potential in the western escarpment of the Ethiopian rift valley | |
CN109725347A (zh) | 一种层间氧化带砂岩型铀矿三维地质体模型构建方法 | |
CN106709988A (zh) | 一种工程地质剖面图构建方法 | |
CN109615696A (zh) | 一种工程地质信息三维可视化处理方法 | |
CN103377196B (zh) | 一种基于地理信息系统的岩相古地理资料处理方法 | |
CN108256259B (zh) | 基于模糊层次分析确定煤层底板注浆靶区的方法及装置 | |
CN102053981A (zh) | 一种顾及语义的地上下集成模型的多尺度表达方法 | |
CN104008147B (zh) | 三维地质结构模型多体索引构建方法 | |
Al-Khazraji et al. | Development of heterogeneous immature Brownfield with Waterdrive using dynamic opportunity index: a case study from Iraqi oilfields | |
CN107194966B (zh) | 一种二维地质断面图地质界线及其属性信息的辨识方法 | |
CN107013207B (zh) | 一种获取复杂断块油藏砂体地下井位信息的方法 | |
CN113253893B (zh) | 一种地质信息显示管理方法、装置、设备及可读存储介质 | |
CN116612248B (zh) | 一种基于标志层和控制点的地层空间扩展与建模方法 | |
Kayode et al. | Automatic Placement of Infill Wells During Simulation Run-Time | |
Bialowas et al. | Greater Sierra Jean Marie Tight Gas Carbonate: Multidisciplinary Approach Drives Decade of Development | |
CN118113798A (zh) | 煤田数字化地质模型中的标准地层定义方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
C17 | Cessation of patent right | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20110706 Termination date: 20140516 |