CN108257080B - 电磁波ct纵剖面图中不良地质体异常快速平面映射方法 - Google Patents

电磁波ct纵剖面图中不良地质体异常快速平面映射方法 Download PDF

Info

Publication number
CN108257080B
CN108257080B CN201711446640.5A CN201711446640A CN108257080B CN 108257080 B CN108257080 B CN 108257080B CN 201711446640 A CN201711446640 A CN 201711446640A CN 108257080 B CN108257080 B CN 108257080B
Authority
CN
China
Prior art keywords
hole
abnormal
holes
drill
plane
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
CN201711446640.5A
Other languages
English (en)
Other versions
CN108257080A (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.)
Changjiang Geophysical Exploration & Testing Wuhan Co ltd
Original Assignee
Changjiang Geophysical Exploration & Testing Wuhan Co ltd
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 Changjiang Geophysical Exploration & Testing Wuhan Co ltd filed Critical Changjiang Geophysical Exploration & Testing Wuhan Co ltd
Priority to CN201711446640.5A priority Critical patent/CN108257080B/zh
Publication of CN108257080A publication Critical patent/CN108257080A/zh
Application granted granted Critical
Publication of CN108257080B publication Critical patent/CN108257080B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G06T3/06—
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N27/00—Investigating or analysing materials by the use of electric, electrochemical, or magnetic means
    • G01N27/72—Investigating or analysing materials by the use of electric, electrochemical, or magnetic means by investigating magnetic variables
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V3/00—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
    • G01V3/12—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation operating with electromagnetic waves
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T11/00—2D [Two Dimensional] image generation
    • G06T11/003—Reconstruction from projections, e.g. tomography

Abstract

一种电磁波CT剖面图中不良地质体异常快速平面映射方法，该方法首先通过对电磁波CT纵剖面图中的不良地质体异常区域的水平范围进行唯一编码，再获取水平范围的参数值并存储为记录序列，之后逐条解译参数信息记录，在平面图中进行不良地质体异常区域映射。该方法可通过计算机编程实现电磁波CT不良地质体异常在平面图中的快速合成，使得电磁波CT不良地质体异常平面图的制图效率大幅提升，减少了平面图制图过程中的错误率以及数据校核工作，提高了质量控制水平。

Description

电磁波CT纵剖面图中不良地质体异常快速平面映射方法
技术领域
本发明涉及一种电磁波CT纵剖面图中不良地质体异常快速平面映射方法，属于工程地球物理探测信息化领域,该技术具体是将电磁波CT纵剖面图中圈画的不良地质体异常区域映射到平面图中，用于将相邻的多个电磁波CT纵剖面中的异常区域合成为一个平面异常区域，形成工程探测范围内不良地质体的平面位置分布图。
背景技术
在城市轨道交通岩溶探测时，通常会在预开挖的隧道两侧布置多条地质钻孔测线，每条测线按一定间隔布置钻孔，采用钻孔电磁波CT技术在钻孔间进行两两对穿，实现对隧道区间岩溶发育情况的全覆盖探测。城市轨道交通岩溶探测成果包含各类图表，工作量最大最费时的是电磁波CT纵剖面异常平面合成图的制作，原因在于：钻孔电磁波CT获取的是两孔间的纵剖面图，在生成不良地质体异常平面图时，采用人工的方式将纵剖面图中不良地质体异常区域投影到平面图的测线上，再将相关的异常闭合成为一个区域。
人工映射的过程如下(以一条测线为例进行说明)：将电磁波CT纵剖面图不良地质体异常区域复制到平面图中(以钻孔孔口位置为复制基点)，沿平面图中不良地质体异常区域左、右两侧分别绘制铅垂线，并任意绘制一条水平线与两条铅垂线相交，从水平线中点引出一条垂线与平面图上的测线相交，交点A，计算中点与两条铅锤线的距离D1(左)，D2(右)，之后以交点A为中心，沿平面图中测线左侧方向移动距离D1并绘制一个点P1，沿平面图中测线右侧方向移动距离D2并绘制一个点P2，连接点P1、A、P2，形成一条折线，完成一条测线的映射过程(见图8)。依次完成邻近测线的不良地质体异常区域的平面映射，再依据异常的相关关系闭合成平面区域，形成平面异常分布图。
由此可见，人工映射方式较为繁琐，效率低。一方面由于轨道交通岩溶探测钻孔较密，电磁波CT纵剖面数量多，不良地质体异常解译量大，通过人工方式将电磁波CT纵剖面图中的不良地质体异常映射到平面图的工作量巨大；另一方面人工映射容易出现错误，一旦发生错误，也难以找到错误点，效率低。因此，急需一种电磁波CT纵剖面图中不良地质体异常快速平面映射方法。
发明内容
本发明的目的是为了解决电磁波CT纵剖面不良地质体异常区域人工映射到平面图过程操作繁琐、易犯错、效率低的问题，所提出了一种电磁波CT纵剖面图不良地质体异常的快速平面映射方法。为实现上述目的，本发明的技术解决方案是：
一种电磁波CT剖面图中不良地质体异常快速平面映射方法，所述平面映射方法按以下步骤进行：
一、对电磁波CT纵剖面图中的不良地质体异常区域的水平范围进行唯一编码
①对于异常区域位于两个钻孔m、n之间的情况，编码规则为“n孔左x米至n孔左y米”。
②对于异常区域跨越一个钻孔m的情况，编码规则为“m孔左x米至m孔右y米”。
③对于异常区域跨越两个钻孔m、n的情况，编码规则为“m孔左x米至n孔右y米”。
④对于异常跨越多个钻孔m、n、n1……、k的情况，编码规则为“m孔左x米过n孔、n1孔……至k孔右y米”。
二、获取电磁波CT纵剖面图中不良地质体异常区域水平范围参数值，并存储为记录序列
①对于异常位于两个钻孔m、n之间的情况，通过异常区域与钻孔n的空间位置关系，获取电磁波CT纵剖面图中异常区域最左端、最右端与钻孔n的水平距离参数值a、b。
②对于异常跨越一个钻孔m的情况，通过异常和跨越异常的钻孔m的空间位置关系，获取电磁波CT纵剖面图中异常区域最左端、最右端与钻孔m的水平距离参数值a、b。
③对于异常跨越两个钻孔m、n的情况，通过异常和跨越异常的两个钻孔m、n的空间位置关系，获取电磁波CT纵剖面图中异常区域最左端与钻孔m、最右端与钻孔n的水平距离参数值a、b。
④对于异常跨越多个钻孔m、n、n1……k的情况，通过异常和跨越异常的多个钻孔m、n、n1……k中的m和k钻孔的空间位置关系，获取电磁波CT纵剖面图中异常区域最左端与钻孔m、最右端与钻孔k的水平距离参数值a、b。
三、逐条解译不良地质体异常区域水平范围参数信息记录
①对于异常区域位于两个钻孔之间的情况，按照上述编码规则将水平范围参数信息记录进行解译，解译钻孔孔号n，空间关系“孔左”以及距离参数a、b。
②对于异常跨越一个钻孔的情况，按照上述编码规则将水平范围参数信息记录进行解译，解译钻孔孔号m，空间关系“孔左、孔右”，以及参数a、b。
③对于异常跨越两个钻孔的情况，按照上述编码规则将水平范围参数信息记录进行解译，解译钻孔孔号m、n,空间关系“孔左、孔右”，以及参数a、b。
④对于异常跨越多个钻孔之间的情况，按照上述编码规则将水平范围参数信息记录进行解译，解译跨越异常的钻孔孔号m、n、……、k，以及参数a、b。
四、在平面图中进行不良地质体异常区域映射
①对于异常区域位于两个钻孔之间的情况，从钻孔平面坐标记录中检索钻孔孔号n所对应的平面坐标，计算沿平面测线距钻孔孔号n左a、b米的两个点的平面坐标，将其与钻孔n的平面坐标连成折线，完成映射过程。
②对于异常跨越一个钻孔的情况，从钻孔平面坐标记录中检索钻孔孔号m所对应的平面坐标，计算沿平面测线距钻孔孔号m左a、右b米的两个点的平面坐标，利用图形软件按将其与钻孔m、n孔的平面坐标连成折线，完成映射过程。
③对于异常跨越两个钻孔的情况，从钻孔平面坐标记录中检索钻孔孔号m、n所对应的平面坐标，计算沿平面测线距钻孔孔号m孔左a米、钻孔孔号n孔右b米的两个点的平面坐标，利用图形软件将其与钻孔m、n孔的平面坐标连成折线，完成映射过程。
④对于异常跨越多个钻孔之间的情况，从钻孔平面坐标记录中检索钻孔孔号m、n、……、k所对应的平面坐标，计算沿平面测线距钻孔孔号m孔左a米，钻孔孔号k孔右b米的两个点的平面坐标，利用图形软件将其与钻孔m、n、……、k孔的平面坐标连成折线，完成映射过程。
所述方法用于电磁波CT或地震波CT或声波CT剖面图。
所述不良地质体异常包括溶洞、溶蚀、孔洞、断层破碎带。
由于采用了以上技术方案，在将电磁波CT纵剖面图不良地质体异常区域合并成平面异常时具有以下优点：
(1)电磁波CT纵剖面异常区域可以通过计算机编程自动投映射到平面图中，为实现电磁波CT不良地质体异常在平面图中快速合成提供了坚实基础。
(2)基于计算机编程使得电磁波CT不良地质体异常平面图的制图效率大幅提升。以往1000对电磁波CT剖面需要2个人连续工作1个月才能完成，采用该技术后，2个人仅用2天便可完成该项工作，提高了生产效率，也确保了生产项目能够按期完成任务。
(3)电磁波CT纵剖面不良地质体异常在映射到平面图时错误率少，减少了返工次数。
(4)质量控制更容易保证。制图校核人员无须对每一个电磁波CT纵剖面不良地质体异常进行计算复核，只需要检查每个电磁波CT纵剖面不良地质体异常的记录是否都参与了制图，质量控制更容易，也提高了成果质量。
完成映射过程后，可依据测线的空间位置关系，将平面图中多条具有相关关系的折线(电磁波CT纵剖面图中异常区域映射而成)的两端端点逆时针或顺时针连接起来，形成平面图中的不良地质体异常。
附图说明
图1为编码规则示意图。
图2为实施例的LL’纵剖面示意图。
图3为实施例的MM’纵剖面示意图。
图4为实施例的RR’纵剖面示意图。
图5为实施例的LML纵剖面示意图。
图6为实施例的RMR纵剖面示意图。
图7为实施例的平面地质钻孔测线布置示意图。
图8为背景技术的人工映射方式示意图。
具体实施方式
下面结合附图对发明做进一步详细介绍。
一种电磁波CT剖面图中不良地质体异常快速平面映射方法，所述平面映射方法按以下步骤进行：
一、对电磁波CT纵剖面图中的不良地质体异常区域的水平范围进行唯一编码
①如附图1(a)所示，对于异常区域位于两个钻孔m、n之间的情况，编码规则为“n孔左x米至n孔左y米”。
②如附图1(b)所示，对于异常区域跨越一个钻孔m的情况，编码规则为“m孔左x米至m孔右y米”。
③如附图1(c)所示，对于异常区域跨越两个钻孔m、n的情况，编码规则为“m孔左x米至n孔右y米”。
④如附图1(d)所示,对于异常跨越多个钻孔m、n、n1……、k的情况，编码规则为“m孔左x米过n孔、n1孔……至k孔右y米”。
所述不良地质体异常包括溶洞、溶蚀、孔洞、断层破碎带。
二、获取电磁波CT纵剖面图中不良地质体异常区域水平范围参数值，并存储为记录序列
①对于异常位于两个钻孔m、n之间的情况，通过异常区域与钻孔n的空间位置关系，获取电磁波CT纵剖面图中异常区域最左端、最右端与钻孔n的水平距离参数值a、b,如附图1(a)所示。
②对于异常跨越一个钻孔m的情况，通过异常和跨越异常的钻孔m的空间位置关系，获取电磁波CT纵剖面图中异常区域最左端、最右端与钻孔m的水平距离参数值a、b,如附图1(b)所示。
③对于异常跨越两个钻孔m、n的情况，通过异常和跨越异常的两个钻孔m、n的空间位置关系，获取电磁波CT纵剖面图中异常区域最左端与钻孔m、最右端与钻孔n的水平距离参数值a、b,如附图1(c)所示。
④对于异常跨越多个钻孔m、n、n1……k的情况，通过异常和跨越异常的多个钻孔m、n、n1……k中的m和k钻孔的空间位置关系，获取电磁波CT纵剖面图中异常区域最左端与钻孔m、最右端与钻孔k的水平距离参数值a、b,如附图1(d)所示。
三、逐条解译不良地质体异常区域水平范围参数信息记录
①对于异常区域位于两个钻孔之间的情况，按照上述编码规则将水平范围参数信息记录进行解译，解译钻孔孔号n，空间关系“孔左”以及距离参数a、b。
②对于异常跨越一个钻孔的情况，按照上述编码规则将水平范围参数信息记录进行解译，解译钻孔孔号m，空间关系“孔左、孔右”，以及参数a、b。
③对于异常跨越两个钻孔的情况，按照上述编码规则将水平范围参数信息记录进行解译，解译钻孔孔号m、n,空间关系“孔左、孔右”，以及参数a、b。
④对于异常跨越多个钻孔之间的情况，按照上述编码规则将水平范围参数信息记录进行解译，解译跨越异常的钻孔孔号m、n、……、k，以及参数a、b。
四、在平面图中进行不良地质体异常区域映射
①对于异常区域位于两个钻孔之间的情况，从钻孔平面坐标记录中检索钻孔孔号n所对应的平面坐标，计算沿平面测线距钻孔孔号n左a、b米的两个点的平面坐标，将其与钻孔n的平面坐标连成折线，完成映射过程。
②对于异常跨越一个钻孔的情况，从钻孔平面坐标记录中检索钻孔孔号m所对应的平面坐标，计算沿平面测线距钻孔孔号m左a、右b米的两个点的平面坐标，利用图形软件按将其与钻孔m、n孔的平面坐标连成折线，完成映射过程。
③对于异常跨越两个钻孔的情况，从钻孔平面坐标记录中检索钻孔孔号m、n所对应的平面坐标，计算沿平面测线距钻孔孔号m孔左a米、钻孔孔号n孔右b米的两个点的平面坐标，利用图形软件将其与钻孔m、n孔的平面坐标连成折线，完成映射过程。
④对于异常跨越多个钻孔之间的情况，从钻孔平面坐标记录中检索钻孔孔号m、n、……、k所对应的平面坐标，计算沿平面测线距钻孔孔号m孔左a米，钻孔孔号k孔右b米的两个点的平面坐标，利用图形软件将其与钻孔m、n、……、k孔的平面坐标连成折线，完成映射过程。
所述方法用于电磁波CT或地震波CT或声波CT剖面图。
完成映射过程后，依据测线的空间位置关系，将平面图中多条具有相关关系的折线(电磁波CT纵剖面图中异常区域映射而成)的两端端点逆时针或顺时针连接起来，形成平面图中的不良地质体异常。
具体实施例
以武汉轨道交通机场线为例，对本技术方案进一步说明。武汉轨道交通机场线布置了L1L1’(L-K1、L-K2、L-K3……L-Kn)、M1M1’(M-K1、M-K2、M-K3……M-Kn)、R1R1’(R-K1、R-K2、R-K3……R-Kn)、LM(L-K1、M-K1、L-K2、M-K2……L-Kn、M-Kn)、RM(R-K1、M-K1、R-K2、M-K2……R-Kn、M-Kn)五条平面地质钻孔测线，如附图7所示。采用电磁波CT方法对钻孔进行两两对穿，形成LL’(L-K1、L-K2、L-K3……L-Kn)、MM’(M-K1、M-K2、M-K3……M-Kn)、RR’(R-K1、R-K2、R-K3……R-Kn)、LML(L-K1、M-K1、L-K2、M-K2……L-Kn、M-Kn)、RMR(R-K1、M-K1、R-K2、M-K2……R-Kn、M-Kn)五条电磁波CT纵剖面。将LL’、MM’、RR’、LML、RMR五条电磁波CT纵剖面不良地质体异常区域映射到平面图中，并形成平面不良地质体异常区域，具体如下：
(1)对纵剖面图中的每个不良地质体异常区域的水平范围进行唯一编码
①对于L1异常区域(如附图2所示)，其跨越了L-K2、L-K3两个钻孔，记录规则为“L-K2孔左x米至L-K3孔右y米”。
②对于M1异常区域(如附图3所示)，其跨越了M-K2、M-K3两个钻孔，记录规则为“M-K2孔左x米至M-K3孔右y米”。
③对于R1异常区域(如附图4所示)，其位于R-K1、R-K2两个钻孔之间，记录规则为“R-K2孔左x米至R-K2孔左y米”。
④对于LM1异常区域(如附图5所示)，其跨越了L-K2一个钻孔，记录规则为“L-K2孔左x米至L-K2孔右y米”。
⑤对于LM2异常区域(如附图5所示)，其跨越了M-K2一个钻孔，记录规则为“M-K2孔左x米至M-K2孔右y米”。
⑥对于LM3异常区域(如附图5所示)，其跨越了L-K3一个钻孔，记录规则为“L-K3孔左x米至L-K3孔右y米”。
⑦对于LM4异常区域(如附图5所示)，其跨越了M-K3一个钻孔，记录规则为“M-K3孔左x米至M-K3孔右y米”。
⑧对于RM1异常区域(如附图6所示)，其跨越了M-K2、R-K3、M-K3三个钻孔，记录规则为“M-K2孔左x米过R-K3孔至M-K3孔右y米”。
(2)在LL’、MM’、RR’、LML、RMR五条纵剖面图中，结合不良地质体异常区域与钻孔之间的空间位置关系，获取不良地质体异常区域水平范围参数值，并存储为记录序列
①对于L1异常区域(如附图2所示)，计算区域最左端与L-K2的距离为d1,区域最右端与L-K3钻孔的距离为d2，异常区域的水平范围记录为“L-K2孔左d1米至L-K3孔右d2米”。
②对于M1异常区域(如附图3所示)，计算区域最左端与M-K2的距离为d3,区域最右端与M-K3钻孔的距离为d4，异常区域的水平范围记录为“M-K2孔左d3米至M-K3孔右d4米”。
③对于R1异常区域(如附图4所示)，计算区域最左端与R-K2的距离为d5,区域最右端与R-K2钻孔的距离为d6，异常区域的水平范围记录为“R-K2孔左d5米至R-K2孔左d6米”。
④对于LM1异常区域(如附图5所示)，计算区域最左端与L-K2的距离为d7,区域最右端与L-K2钻孔的距离为d8，异常区域的水平范围记录为“L-K2孔左d7米至L-K2孔右d8米”。
⑤对于LM2异常区域(如附图5所示)，计算区域最左端与M-K2的距离为d9,区域最右端与M-K2钻孔的距离为d10，异常区域的水平范围记录为“M-K2孔左d9米至M-K2孔右d10米”。
⑥对于LM3异常区域(如附图5所示)，计算区域最左端与L-K3的距离为d11,区域最右端与L-K3钻孔的距离为d12，异常区域的水平范围记录为“L-K3孔左d11米至L-K3孔右d12米”。
⑦对于LM4异常区域(如附图5所示)，计算区域最左端与M-K3的距离为d13,区域最右端与M-K3钻孔的距离为d14，异常区域的水平范围记录为“M-K3孔左d13米至M-K3孔右d14米”。
⑧对于RM1异常区域(如附图6所示)，计算区域最左端与M-K2的距离为d15,区域最右端与M-K3钻孔的距离为d16，异常区域的水平范围记录为“M-K2孔左d15米过R-K3孔至M-K3孔右d16米”。
(3)按照上述编码规则，逐条解译不良地质体异常区域水平参数信息记录
①对于L1异常区域(如附图2所示)，解译异常区域水平范围记录“L-K2孔左d1米至L-K3孔右d2米”中的孔号“L-K2”、“L-K3”，空间关系“孔左”、“孔右”以及距离参数d1、d2的值。
②对于M1异常区域(如附图3所示)，解译异常区域水平范围记录“M-K2孔左d3米至M-K3孔右d4米”中的孔号“M-K2”、“M-K3”，空间关系“孔左”、“孔右”以及距离参数d3、d4的值。
③对于R1异常区域(如附图4所示)，解译异常区域水平范围记录“R-K2孔左d5米至R-K2孔左d6米”中的孔号“R-K2”，空间关系“孔左”以及距离参数d5、d6的值。
④对于LM1异常区域(如附图5所示)，解译异常区域水平范围记录“L-K2孔左d7米至L-K2孔右d8米”中的孔号“L-K2”，空间关系“孔左”、“孔右”以及距离参数d7、d8的值。
⑤对于LM2异常区域(如附图5所示)，解译异常区域水平范围记录“M-K2孔左d9米至M-K2孔右d10米”中的孔号“M-K2”，空间关系“孔左”、“孔右”以及距离参数d9、d10的值。
⑥对于LM3异常区域(如附图5所示)，解译异常区域水平范围记录“L-K3孔左d11米至L-K3孔右d12米”中的孔号“L-K3”，空间关系“孔左”、“孔右”以及距离参数d11、d12的值。
⑦对于LM4异常区域(如附图5所示)，解译异常区域水平范围记录“M-K3孔左d13米至M-K3孔右d14米”中的孔号“M-K3”，空间关系“孔左”、“孔右”以及距离参数d13、d14的值。
⑧对于RM1异常区域(如附图6所示)，解译异常区域水平范围记录“M-K2孔左d15米过R-K3孔至M-K3孔右d16米”中的孔号“M-K2”、“R-K3”、“M-K3”，空间关系“孔左”、“孔右”以及距离参数d15、d16的值。
(4)将LL’、MM’、RR’、LML、RMR五条纵剖面图中的异常区域映射到L1L1’、M1M1’、R1R1’、LM、RM平面测线上(如附图7所示)。
①对于L1异常区域，沿L1L1’测线L-K2～L-K3段，计算距L-K2孔左d1米lp1点的平面坐标，距L-K3孔右d2米lp2点的平面坐标，连接点lp1、L-K2、L-K3、lp2形成折线，完成L1异常区域的平面映射。
②对于M1异常区域，沿M1M1’测线M-K2～M-K3段，计算距M-K2孔左d3米mp1点的平面坐标，距M-K3孔右d4米mp2点的平面坐标,连接点mp1、M-K2、M-K3、mp2形成折线，完成M1异常区域的平面映射。
③对于R1异常区域，沿R1R1’测线R-K1～R-K2段，计算距R-K2孔左d5米rp2点的平面坐标，距R-K2孔左d6米rp1点的平面坐标,连接点rp1、rp2、R-K2形成折线，完成R1异常区域的平面映射。
④对于LM1异常区域，沿LM测线M-K1～L-K2～M-K2段，计算距L-K2孔左d7米lmp1点的平面坐标，距L-K2孔右d8米lmp2点的平面坐标,连接点lmp1、L-K2、lmp2形成折线，完成LM1异常区域的平面映射。
⑤对于LM2异常区域，沿LM测线L-K2～M-K2～L-K2段，计算距M-K2孔左d9米lmp5点的平面坐标，距M-K2孔右d10米lmp6点的平面坐标,连接点lmp5、M-K2、lmp6形成折线，完成LM2异常区域的平面映射。
⑥对于LM3异常区域，沿LM测线M-K2～L-K3～M-K3段，计算距L-K3孔左d11米lmp3点的平面坐标，距L-K3孔右d12米lmp4点的平面坐标,连接点lmp3、L-K3、lmp4形成折线，完成LM3异常区域的平面映射。
⑦对于LM4异常区域，沿LM测线L-K3～M-K3～L-K4段，计算距M-K3孔左d13米lmp7点的平面坐标，距M-K3孔右d14米lmp8点的平面坐标,连接点lmp7、M-K3、lmp8形成折线，完成LM4异常区域的平面映射。
⑧对于RM1异常区域，沿RM测线R-K2～M-K2～R-K3～M-K3～R-K4段，计算距R-K2孔左d15米rmp1点的平面坐标，距R-K4孔右d16米rmp2点的平面坐标,连接点rmp1、M-K2～R-K3～M-K3、rmp2形成折线，完成RM1异常区域的平面映射。
完成映射过程后，依据测线的空间位置关系，将平面图中多条具有相关关系的折线(剖面图中异常区域映射而成)的两端端点逆时针或顺时针连接起来，形成平面图不良地质体异常(如附图7所示)。
①连接点lp1、lmp1、lmp2、lmp3、lmp4、lp2，形成平面不良地质体异常区域P1001。
②连接点lmp5、lmp6、lmp7、lmp8、mp2、rmp2、R-K3、rmp1、mp1，形成平面不良地质体异常区域P1002。
③连接点rp1、rp2，形成平面不良地质体异常区域P1003。
所述方法除用于电磁波CT或地震波CT或声波CT。
所述不良地质体异常包括溶洞、溶蚀、孔洞、断层破碎带。
由于采用了以上技术方案，在将武汉轨道机场线的电磁波CT纵剖面图不良地质体异常区域合并成平面异常时具有明显优点：
(1)通过计算机编程实现了电磁波CT不良地质体异常在平面图中的快速合成；
(2)基于计算机编程可实现电磁波CT不良地质体异常平面图的制图效率大幅提升。以往1000对电磁波CT剖面需要2个人连续工作1个月才能完成，采用该技术后，2个人仅用2天便可完成该项工作，提高了生产效率，也确保了生产项目能够按期完成任务；
(3)平面图制图过程中错误率少，减少了返工次数；
(4)减少了制图过程中的数据校核工作，提高了质量控制水平。

Claims (3)

1.一种电磁波CT剖面图中不良地质体异常快速平面映射方法，其特征在于，所述平面映射方法按以下步骤进行：
一、对电磁波CT纵剖面图中的不良地质体异常区域的水平范围进行唯一编码
①对于异常区域位于两个钻孔m、n之间的情况，编码规则为“n孔左x米至n孔左y米”；
②对于异常区域跨越一个钻孔m的情况，编码规则为“m孔左x米至m孔右y米”；
③对于异常区域跨越两个钻孔m、n的情况，编码规则为“m孔左x米至n孔右y米”；
④对于异常跨越多个钻孔m、n、n1……、k的情况，编码规则为“m孔左x米过n孔、n1孔……至k孔右y米”；
二、获取电磁波CT纵剖面图中不良地质体异常区域水平范围参数值，并存储为记录序列
①对于异常位于两个钻孔m、n之间的情况，通过异常区域与钻孔n的空间位置关系，获取电磁波CT纵剖面图中异常区域最左端、最右端与钻孔n的水平距离参数值a、b；
②对于异常跨越一个钻孔m的情况，通过异常和跨越异常的钻孔m的空间位置关系，获取电磁波CT纵剖面图中异常区域最左端、最右端与钻孔m的水平距离参数值a、b；
③对于异常跨越两个钻孔m、n的情况，通过异常和跨越异常的两个钻孔m、n的空间位置关系，获取电磁波CT纵剖面图中异常区域最左端与钻孔m、最右端与钻孔n的水平距离参数值a、b；
④对于异常跨越多个钻孔m、n、n1……k的情况，通过异常和跨越异常的多个钻孔m、n、n1……k中的m和k钻孔的空间位置关系，获取电磁波CT纵剖面图中异常区域最左端与钻孔m、最右端与钻孔k的水平距离参数值a、b；
三、逐条解译不良地质体异常区域水平范围参数信息记录
①对于异常区域位于两个钻孔之间的情况，按照上述编码规则将水平范围参数信息记录进行解译，解译钻孔孔号n，空间关系“孔左”以及距离参数a、b；
②对于异常跨越一个钻孔的情况，按照上述编码规则将水平范围参数信息记录进行解译，解译钻孔孔号m，空间关系“孔左、孔右”，以及参数a、b；
③对于异常跨越两个钻孔的情况，按照上述编码规则将水平范围参数信息记录进行解译，解译钻孔孔号m、n,空间关系“孔左、孔右”，以及参数a、b；
④对于异常跨越多个钻孔之间的情况，按照上述编码规则将水平范围参数信息记录进行解译，解译跨越异常的钻孔孔号m、n、……、k，以及参数a、b；
四、在平面图中进行不良地质体异常区域映射
①对于异常区域位于两个钻孔之间的情况，从钻孔平面坐标记录中检索钻孔孔号n所对应的平面坐标，计算沿平面测线距钻孔孔号n左a、b米的两个点的平面坐标，将其与钻孔n的平面坐标连成折线，完成映射过程；
②对于异常跨越一个钻孔的情况，从钻孔平面坐标记录中检索钻孔孔号m所对应的平面坐标，计算沿平面测线距钻孔孔号m左a、右b米的两个点的平面坐标，利用图形软件按将其与钻孔m、n孔的平面坐标连成折线，完成映射过程；
③对于异常跨越两个钻孔的情况，从钻孔平面坐标记录中检索钻孔孔号m、n所对应的平面坐标，计算沿平面测线距钻孔孔号m孔左a米、钻孔孔号n孔右b米的两个点的平面坐标，利用图形软件将其与钻孔m、n孔的平面坐标连成折线，完成映射过程；
④对于异常跨越多个钻孔之间的情况，从钻孔平面坐标记录中检索钻孔孔号m、n、……、k所对应的平面坐标，计算沿平面测线距钻孔孔号m孔左a米，钻孔孔号k孔右b米的两个点的平面坐标，利用图形软件将其与钻孔m、n、……、k孔的平面坐标连成折线，完成映射过程。
2.如权利要求1所述的一种电磁波CT剖面图中不良地质体异常快速平面映射方法，其特征在于：所述方法用于电磁波CT或地震波CT或声波CT剖面图。
3.如权利要求1所述的一种电磁波CT剖面图中不良地质体异常快速平面映射方法，其特征在于：所述不良地质体异常包括溶洞、溶蚀、孔洞、断层破碎带。
CN201711446640.5A 2017-12-27 2017-12-27 电磁波ct纵剖面图中不良地质体异常快速平面映射方法 Active CN108257080B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201711446640.5A CN108257080B (zh) 2017-12-27 2017-12-27 电磁波ct纵剖面图中不良地质体异常快速平面映射方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201711446640.5A CN108257080B (zh) 2017-12-27 2017-12-27 电磁波ct纵剖面图中不良地质体异常快速平面映射方法

Publications (2)

Publication Number Publication Date
CN108257080A CN108257080A (zh) 2018-07-06
CN108257080B true CN108257080B (zh) 2021-02-23

Family

ID=62724403

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201711446640.5A Active CN108257080B (zh) 2017-12-27 2017-12-27 电磁波ct纵剖面图中不良地质体异常快速平面映射方法

Country Status (1)

Country Link
CN (1) CN108257080B (zh)

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101066211A (zh) * 2006-05-05 2007-11-07 通用电气公司 用于在超声波系统中显示信息的用户界面及方法
US7853045B2 (en) * 2007-10-31 2010-12-14 Saudi Arabian Oil Company Geostatistical analysis and classification of core data
CN103038670A (zh) * 2010-05-28 2013-04-10 阿克斯有限责任公司 对地球物理数据进行处理
CN104299269A (zh) * 2014-09-30 2015-01-21 长江岩土工程总公司（武汉） 一种三维地质模型的地质剖面输出方法
CN104481534A (zh) * 2014-11-06 2015-04-01 中国矿业大学 一种采煤机滚筒自动调高系统
CN205825951U (zh) * 2016-05-26 2016-12-21 武汉固德超前高新科技研发有限公司 用于勘探溶洞内部结构的探测装置
CN106709988A (zh) * 2015-11-16 2017-05-24 天津市勘察院 一种工程地质剖面图构建方法

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101066211A (zh) * 2006-05-05 2007-11-07 通用电气公司 用于在超声波系统中显示信息的用户界面及方法
US7853045B2 (en) * 2007-10-31 2010-12-14 Saudi Arabian Oil Company Geostatistical analysis and classification of core data
CN103038670A (zh) * 2010-05-28 2013-04-10 阿克斯有限责任公司 对地球物理数据进行处理
CN104299269A (zh) * 2014-09-30 2015-01-21 长江岩土工程总公司（武汉） 一种三维地质模型的地质剖面输出方法
CN104481534A (zh) * 2014-11-06 2015-04-01 中国矿业大学 一种采煤机滚筒自动调高系统
CN106709988A (zh) * 2015-11-16 2017-05-24 天津市勘察院 一种工程地质剖面图构建方法
CN205825951U (zh) * 2016-05-26 2016-12-21 武汉固德超前高新科技研发有限公司 用于勘探溶洞内部结构的探测装置

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
CT探测覆岩破坏高度的试验研究;文学宽;《煤炭学报》;19980630;第23卷(第3期);第300-304页 *
Multiscale Hessian fracture filtering for the enhancement and segmentation of narrow fractures in 3D image data;Maarten Voorn,et al.;《Computers & Geosciences》;20130326;第57卷;第44-53页 *
井间地震CT技术及其在铁路岩溶勘察中的应用;张连伟;《中国优秀硕士学位论文全文数据库 工程科技II辑》;20150515(第5期);C033-6 *
采矿工程剖面图计算机生成方法研究;王琳 等;《煤矿机械》;20171130;第38卷(第11期);第159-161页 *
钻孔歪斜在地质部面上的投影方法;赵莉萍 等;《煤炭科技》;19991231(第2期);第48-49页 *

Also Published As

Publication number Publication date
CN108257080A (zh) 2018-07-06

Similar Documents

Publication Publication Date Title
WO2023185735A1 (zh) 全时空采掘过程的矿山应力场孪生建模同化系统及方法
CN103514630B (zh) 断层构造三维建模方法
Ponting Corner point geometry in reservoir simulation
CN102622526B (zh) 数字矿山中的巷道开掘搜索预测方法
CN114663627B (zh) 一种基于三维点云数据库的矿山数字模型建立方法
US10641923B2 (en) Method for characterizing and exploiting a subterranean formation comprising a network of fractures
CN103184887B (zh) 井下巷道钻孔施工及地质信息反演方法
Du et al. A workflow for integrated Barnett Shale gas reservoir modeling and simulation
Maruvanchery et al. Early construction cost and time risk assessment and evaluation of large-scale underground cavern construction projects in Singapore
CN111814298B (zh) 近水平钻孔轨迹约束下采煤工作面煤层迭代建模方法
WO2021072921A1 (zh) 基于bim+gis融合技术的灾害发生追溯方法
CN105064990A (zh) 一种依靠水平井三维可视化地层对比的构造建模方法
CN105321204A (zh) 一种三维地质栅格体模型构建方法
CN110188426A (zh) 一种隧道地质条件协同预测方法
CN106709988A (zh) 一种工程地质剖面图构建方法
CN108257080B (zh) 电磁波ct纵剖面图中不良地质体异常快速平面映射方法
CN108256130B (zh) 一种地质钻孔空间轨迹还原方法
CN104462649B (zh) 一种矿体块段模型储量的自动更新方法
CN113298935B (zh) 一种基于点云成像技术确定巷道岩体结构特征参数的方法
CN110687591B (zh) 基于先验数据的波形匹配确定煤层及围岩物性参数的方法
CN111414705B (zh) 一种丛式井井眼轨迹防碰设计的三维空间可穿越评估方法
CN112862964B (zh) 一种顺煤层巷道掘进地质剖面动态更新方法
Zhifa et al. Three-dimensional back-analysis of displacements in exploration adits—principles and application
CN113989446A (zh) 提高地质三维模型地表精度的地质界线野外复核方法
RU2810922C1 (ru) Способ управления и контроля рудником на основе системы обеспечения геологических съемок и построения интеллектуальной платформы и устройство для его осуществления

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