CN105572753A - 一种瞬变电磁探测数据等值线成图与高阻区域识别方法 - Google Patents
一种瞬变电磁探测数据等值线成图与高阻区域识别方法 Download PDFInfo
- Publication number
- CN105572753A CN105572753A CN201510916764.XA CN201510916764A CN105572753A CN 105572753 A CN105572753 A CN 105572753A CN 201510916764 A CN201510916764 A CN 201510916764A CN 105572753 A CN105572753 A CN 105572753A
- Authority
- CN
- China
- Prior art keywords
- isoline
- color
- apparent resistivity
- region
- new
- 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
Classifications
-
- 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/38—Processing data, e.g. for analysis, for interpretation, for correction
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Image Analysis (AREA)
- Developing Agents For Electrophotography (AREA)
Abstract
本发明实施例公开了一种瞬变电磁探测数据等值线成图与高阻区域识别方法,包括,把不同测点不同时间段的感应电动势数据反演为空间离散的视电阻率数据，再利用不规则三角网绘制等值线，建立单等值线集数组、以及与所述单等值线集数组相对应的存储颜色值的颜色数组，最后通过视电阻率比较和颜色查找方法进行颜色填充，提高了瞬变电磁法探测数据绘制为等值线图的精度与速度；在填色过程中根据给定的高阻区域视电阻率值下限以及相对应的指定颜色值，凡大于或等于该视电阻率值下限值的区域都置为所述指定颜色值，本发明实施例对于高阻区域的识别给出了简便方法，对于有些高阻区，如空洞区域等能够根据其视电阻率范围进行快速识别。
Description
技术领域
本发明涉及物探成图技术领域，特别是涉及一种基于不规则三角网的瞬变电磁探测数据等值线成图与高阻区域识别方法。
背景技术
瞬变电磁法通过不接地回线或接地线源向空间发送一次脉冲磁场，在一次脉冲磁场间歇期间，采用不接地线圈接收感应二次磁场，并在线圈中形成感应电动势。一个测点上有不同时间间隔的离散感应电动势数据，通过各个时间段的感应电动势随时间变化规律，可得到该测点不同深度的视电阻率值。收集一条测线上多个测点的视电阻率数据，可以得到二维剖面空间分布的离散视电阻率数据。
等值线图能用图形的形式来反映离散数据，在GIS、气象学、物探等学科中有着广泛的应用。在物探领域，等值线图能够直观地显示测得的地质构造的某一参量(如视电阻率)值。对瞬变电磁法探测得到的二维剖面离散视电阻率进行等值线图的绘制，可以更直观地反映地质结构特征。等值线的绘制越精确，越能保证数据值分布情况以图形方式反映的准确性。等值线围成区域颜色表示区域的数据视电阻率范围，对区域颜色必须进行准确填充。
等值线的绘制与填充通常采用规则网格法，但很多时候，需要在一个由不规则测量点组成的区域内进行等值线的绘制。在这种情况下，如采用规则矩形网格方式绘制等值线，首先要用空间插值法将离散点数据内插成格点资料，但经过网格化后，等值线在经过测点时往往与实测值不符，其误差大小与网格间距关系密切，为了接近实际值，网格间距应尽可能小，但这势必降低运算速度，而且也不能解决不规则边界的问题。而且，如果进行颜色填充，必然造成等值线围成区域颜色反映的物理数据也存在一定误差。采用TIN(TriangularIrregularNetwork，不规则三角网)方式，无需对离散数据网格化，更能忠实于原始离散数据的分布。基于不规则三角网的属性值趋势判断方法是根据围成区域等值线的视电阻率差来判断区域颜色的，但是对于没有显著趋势的由一条单等值线围成的区域颜色判断往往不够准确。
发明内容
本发明实施例中提供了一种瞬变电磁探测数据等值线成图与高阻区域识别方法，在保证瞬变电磁法探测数据绘制为等值线图与高阻区域识别速度的同时，提高待填充区域颜色判断的精度。
为了解决上述技术问题，本发明实施例公开了如下技术方案：
一种基于不规则三角网的瞬变电磁探测数据等值线成图与高阻区域识别方法，包括：
用瞬变电磁仪获取待绘制区域若干测点对应多个测道的感应电压数据；
将所述感应电压数据传输给数据处理单元，所述数据处理单元将所述感应电压数据反演成离散的视电阻率数据；
以测点号为空间横坐标、同时以测道号为视深坐标确定所述视电阻率数据的空间分布，生成离散数据点集；
根据所述离散数据点集生成不规则三角网TIN和有不同视电阻率的单等值线集数组Isoline[n]；
根据所述不规则三角网TIN，获得TIN'的边界点集ST；
根据所述单等值线集数组，获得等值线端点集SL；
将所述等值线端点集SL插入所述边界点集ST中，合并形成新的点集S；
对所述点集S中的边界点和等值线端点进行空间上的排序，形成有序点集S’；
根据所述单等值线集数组Isoline[n]，对应建立存储颜色值的颜色数组Color[n+1]，给出指定高阻区域视电阻率最小值D以及与所述视电阻率最小值D对应的颜色ColorS；
根据所述有序点集S’中的所有等值线端点，搜索待填充开区域Ri，并通过视电阻率值比较和颜色查找法确定所述填充开区域Ri的填充色，其中，当所述填充开区域Ri的视电阻率大于或等于D时，则所述填充开区域Ri的填充色为所述颜色ColorS；
判断所述有序点集S’中是否存在封闭等值线；
如果所述有序点集S’中存在封闭等值线，则通过视电阻率值比较和颜色查找法确定各所述封闭等值线所围成封闭区域的填充色，其中，当所述封闭区域的视电阻率大于或等于D时，则所述封闭区域的填充色为所述颜色ColorS；
如果所述有序点集S’中不存在封闭等值线，则成图结束。
优选地，所述根据所述有序点集S’中的所有等值线线端点，搜索待填充开区域Ri，并通过视电阻率值比较和颜色查找法确定所述填充开区域Ri的填充色，其中，当所述填充开区域Ri的视电阻率大于或等于D时，则所述填充开区域Ri的填充色为所述颜色ColorS，包括：
依次序获取所述有序点集S’中新等值线端点；
判断所获取的新等值线端点是否为已获取的重复等值线端点；
如果所获取的新等值线端点为已获取的重复等值线端点，则结束开区域填充；
如果所获取的新等值线端点不是已获取的重复等值线端点，则以所述新等值线端点为起点搜索一个待填充开区域Ri；
获取在所述有序点集S’中、与所述新等值线相邻且排在所述新等值线端点前面的一个有序点，其中，所述有序点为等值线端点或边界点；
确定所述新等值线端点所在单等值线集的索引index；
判断所述新等值线端点的视电阻率是否大于所述有序点的视电阻率；
如果所述新等值线端点的视电阻率大于所述有序点的视电阻率，则用index查找所述颜色数组来确定所述开区域Ri的颜色；
如果所述新等值线端点的视电阻率小于或等于所述有序点的视电阻率，则用index+1查找所述颜色数组来确定所述开区域Ri的颜色；
判断所述开区域Ri对应的视电阻率是否大于或等于D；
如果所述开区域Ri对应的视电阻率大于或等于D，则所述填充开区域Ri的填充色为所述颜色ColorS。
优选地，所述如果所述有序点集S’中存在封闭等值线，则通过视电阻率值比较和颜色查找法确定各所述封闭等值线所围成封闭区域的填充色，其中，当所述封闭区域的视电阻率大于或等于D时，则所述封闭区域的填充色为所述颜色ColorS，包括：
把所有封闭等值线按照在所述空间横坐标中的位置从小到大依次排序，并存入队列Q；
依次序获取所述队列Q中新封闭等值线；
判断所获取的新封闭等值线是否为已获取的重复封闭等值线；
如果所获取的新封闭等值线为已获取的重复封闭等值线，则结束封闭等值线填充；
如果所获取的新封闭等值线不是已获取的重复封闭等值线，则取所述新封闭等值线的线头点所在TIN线段的两个端点T1和T2；
取T1和T2中处于所述新封闭等值线所围成区域内部的点T，确定所述线头点所在单等值线集的索引CIndex；
判断所述线头点的视电阻率是否大于所述点T的视电阻率；
如果所述线头点的视电阻率大于所述点T的视电阻率，则用CIndex查找所述颜色数组来确定所述新封闭等值线所围成区域的颜色；
如果所述线头点的视电阻率小于或等于所述点T的视电阻率，则用CIndex+1查找所述颜色数组来确定所述新封闭等值线所围成区域的颜色；
判断所述新封闭等值线所围成区域对应的视电阻率是否大于或等于D；
如果所述新封闭等值线所围成区域对应的视电阻率大于或等于D，则所述新封闭等值线所围成区域的填充色为所述颜色ColorS。
优选地，所述对所述点集S中的边界点和等值线端点进行空间上的排序，包括：
对所述点集S中的边界点和等值线端点进行空间上的逆时针排序。
优选地，所述根据所述单等值线集数组Isoline[n]，对应建立存储颜色值的颜色数组Color[n+1]，所述单等值线集数组Isoline[n]和所述颜色数组Color[n+1]的对应关系包括：
所述单等值线集数组中一条单等值线Isoline[i]两侧的区域颜色为Color[i]和Color[i+1]。
优选地，所述依次序获取所述有序点集S’中新等值线端点，包括依逆时针次序获取所述有序点集S’中新等值线端点。
优选地，所述成图结束之后，还包括把成图与高阻区识别结果显示在终端屏幕和/或保存到数据库中。
由以上技术方案可见，本发明实施例提供的一种瞬变电磁探测数据等值线成图与高阻区域识别方法,把不同测点不同时间段的感应电动势数据反演为空间离散视电阻率数据，再利用不规则三角网绘制等值线，建立单等值线集数组、以及与所述单等值线集数组相对应的存储颜色值的颜色数组，通过视电阻率比较与颜色查找方法来填充颜色，提高了瞬变电磁法探测数据绘制为等值线图的精度与速度；在填色过程中根据给定高阻区域视电阻率值下限以及相对应的指定颜色值，凡大于等于该值的区域都置为所述指定颜色值，本发明实施例对于高阻区域的识别给出了简便方法，对于有些高阻区，如空洞区域等能够根据其视电阻率范围进行快速识别。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，对于本领域普通技术人员而言，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1为本发明实施例提供的一种瞬变电磁探测数据等值线成图与高阻区域识别方法的流程示意图；
图2为本发明实施例提供的一种等值线开区域颜色填充的流程示意图；
图3为本发明实施例提供的一种等值线封闭区域颜色填充的流程示意图；
图4为本发明实施例提供的一种瞬变电磁探测系统的基本结构示意图。
具体实施方式
为了使本技术领域的人员更好地理解本发明中的技术方案，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都应当属于本发明保护的范围。
参见图1，为本发明实施例提供的一种瞬变电磁探测数据等值线成图与高阻区域识别方法的流程示意图，具体包括：
步骤S101：获取待绘制区域若干测点对应多个测道的感应电压数据。
具体的，可以用瞬变电磁仪探测待绘制区域中若干测点对应多个测道的感应电压数据。
步骤S102：将所述感应电压数据传输给数据处理单元，所述数据处理单元将所述感应电压数据反演成离散的视电阻率数据。
具体的，将所述感应电压数据传输给数据处理单元，通过所述数据处理单元中预先加载好的反演软件通过一定的算法将所述感应电压数据反演成离散的视电阻率数据。
步骤S103：以测点号为空间横坐标、同时以测道号为视深坐标确定所述视电阻率数据的空间分布，生成离散数据点集。
步骤S104：根据所述离散数据点集生成不规则三角网TIN和有不同视电阻率的单等值线集数组Isoline[n]。
具体的，将所述离散数据点集根据三角网生长算法生成不规则三角网TIN；根据所述离散数据点集生成有n种不同视电阻率的单等值线集数组Isoline[n]，所述单等值线集数组Isoline[n]中的所有单等值线集按照视电阻率的大小升序排列，其中Isoline[i](i>＝0,i<n)包含一个具有相同视电阻率的单等值线集SingleLines[mi]，mi为第i个单等值线集的单等值线条数。
步骤S105：根据所述不规则三角网TIN，获得TIN'的边界点集ST。
步骤S106：根据所述单等值线集数组，获得等值线端点集SL。
步骤S107：将所述等值线端点集SL插入所述边界点集ST中，合并形成新的点集S。
步骤S108：对所述点集S中的边界点和等值线端点进行空间上的排序，形成有序点集S’。
具体的，对所述点集S中的边界点和等值线端点在步骤S103中坐标系中的位置进行空间上的逆时针排序，形成有序点集S’，当然也可以是空间上的顺时针排序。
步骤S109：根据所述单等值线集数组Isoline[n]，对应建立存储颜色值的颜色数组Color[n+1]，给出指定高阻区域视电阻率最小值D以及与所述视电阻率最小值D对应的颜色ColorS。
其中，根据所述单等值线集数组Isoline[n]，对应建立存储颜色值的颜色数组Color[n+1]，所述单等值线集数组Isoline[n]与所述颜色数组Color[n+1]的对应关系表述为Isoline[i]中一条单等值线两侧的区域颜色为Color[i]和Color[i+1]。
步骤S110：根据所述有序点集S’中的所有等值线端点，搜索待填充开区域Ri，并通过视电阻率值比较和颜色查找法确定所述填充开区域Ri的填充色，其中，当所述填充开区域Ri的视电阻率大于或等于D时，则所述填充开区域Ri的填充色为所述颜色ColorS。
如图2所示，具体包括如下步骤：
步骤S201：依次序获取所述有序点集S’中新等值线端点。
具体的，依照逆时针次序获取所述有序点集S’中新等值线端点，在本实施例中，获取的所述新等值线端点，在本实施例例中所述新等值线端点命名为等值线端点A，当然，还可以依照顺时针次序获取所述有序点集S’中新等值线端点。
步骤S202：判断所获取的新等值线端点是否为已获取的重复等值线端点。
如果所获取的所述等值线端点A之前已获取的重复等值线端点，即两次获取到同一个等值线端点，则执行步骤S203、即结束开区域填充。
如果所获取的所述等值线端点A不是已获取的重复等值线端点，则执行步骤S204、即以所述等值线端点A为起点搜索一个待填充开区域Ri。
步骤S205：获取在所述有序点集S’中、与所述新等值线相邻且排在所述新等值线端点前面的一个有序点A'，其中，所述有序点为等值线端点或边界点。
步骤S206：确定所述等值线端点A所在单等值线集的索引index。
具体为，确定所述等值线端点A在所述单等值线集数组中具体位于哪一单等值线中。
步骤S207：判断所述等值线端点A的视电阻率是否大于所述有序点A'的视电阻率。
如果所述等值线端点A的视电阻率大于所述有序点A'的视电阻率，则执行步骤S208：用index查找所述颜色数组来确定所述开区域Ri的颜色。
如果所述新等值线端点的视电阻率小于或等于所述有序点的视电阻率，则则执行步骤S209：用index+1查找所述颜色数组来确定所述开区域Ri的颜色。
通过步骤S205至步骤S209中的视电阻率比较法，确定开区域对应的电阻率，然后通过步骤S109中建立的视电阻率和颜色数组的对应关系，确定开区域的填充色，提高了绘图速度，同时因为采用与所获取的等值线端点邻近点的视电阻率比较、以及基于不规则三角网的采样方法，提高了开区域的绘图精度。
步骤S210：判断所述开区域Ri对应的视电阻率是否大于或等于D。
如果所述开区域Ri对应的视电阻率大于或等于指定高阻区域视电阻率最小值D，则执行步骤S211，将所述填充开区域Ri的填充色定位为所述颜色ColorS，这样便可以对开区域中的高阻区置为指定的颜色。
将所有开区域填充完毕后，再进行封闭等值线围成区域的填充。
步骤S111：判断所述有序点集S’中是否存在封闭等值线。
如果所述有序点集S’中不存在封闭等值线，则执行步骤S113，成图结束。
其中，成图结束后，还可以成图识别结果显示在终端屏幕上和/或保存到数据库中。
如果所述有序点集S’中存在封闭等值线，则执行步骤S112：通过视电阻率值比较和颜色查找法确定各所述封闭等值线所围成封闭区域的填充色，其中，当所述封闭区域的视电阻率大于或等于D时，则所述封闭区域的填充色为所述颜色ColorS。
如图3所示，具体可以包括如下步骤：
步骤S301：把所有封闭等值线按照在所述空间横坐标中的位置从小到大依次排序，并存入队列Q。
步骤S302：依次序获取所述队列Q中新封闭等值线。
步骤S303：判断所获取的新封闭等值线是否为已获取的重复封闭等值线。
如果所获取的新封闭等值线为已获取的重复封闭等值线，则执行步骤S304：结束封闭等值线填充，即成图结束。
如果所获取的新封闭等值线不是已获取的重复封闭等值线，则执行步骤S305：取所述新封闭等值线的线头点所在TIN线段的两个端点T1和T2，其中，所述线头点为在所述新封闭等值线中选取的任一点。
步骤S306：取T1和T2中处于所述新封闭等值线所围成区域内部的点、并命名为点T，确定所述线头点所在单等值线集的索引CIndex。
其中，所述确定所述线头点所在单等值线集的索引CIndex，具体为确定所述线头点在所述单等值线集数组中具体位于哪一封闭单等值线中。
步骤S307：判断所述线头点的视电阻率是否大于所述点T的视电阻率。
如果所述线头点的视电阻率大于所述点T的视电阻率，则执行步骤S308：用CIndex查找所述颜色数组来确定所述新封闭等值线所围成区域的颜色。
如果所述线头点的视电阻率小于或等于所述点T的视电阻率，则则执行步骤S309：用CIndex+1查找所述颜色数组来确定所述新封闭等值线所围成区域的颜色。
通过步骤S306至步骤S309中的视电阻率比较法，确定封闭等值线对应的电阻率，然后通过步骤S109中建立的视电阻率和颜色数组的对应关系，确定封闭等值线所围成区域的填充色，提高了绘图速度；同时因为步骤S305和S306中基于不规则三角网的视电阻率比较点的选取方法，提高了绘图精度。
步骤S310：判断所述新封闭等值线所围成区域对应的视电阻率是否大于或等于D。
如果所述新封闭等值线所围成区域对应的视电阻率大于或等于D，则执行步骤S311：所述新封闭等值线所围成区域的填充色为所述颜色ColorS这样便可以对封闭等值线所围封闭区域中的高阻区置为指定的颜色。
为对本发明实施例提供的成图和高阻区识别方法作进一步说明，下面将结合具体实施例说明瞬变电磁探测数据等值线成图与高阻区域识别方法的过程。
如图4所示，瞬变电磁仪1以蓄电池3供电、电压为12V，发射与接收线圈2的匝数均为1匝，形状均为正方形。其边长均为20m。
瞬变电磁仪1将采集的感应电压数据发送给计算机4；数据库5与所述计算机4相连接，用于存储成图和高阻区识别结果。
用图4中的瞬变电磁探测系统对某一指定区域进行探测，离散视电阻率数据的最大视电阻率maxZ预设为616、最小视电阻率minZ预设为340，等值线视电阻率个数n预设为16，颜色表Color[n+1]有17种不同的颜色，高阻区域视电阻率最小值D设为550，指定高阻区域颜色为一个颜色，然后依据上述成图与高阻区识别方法，便可以获取识别出高阻区的等值线填充图。
本发明实施例提供的瞬变电磁探测数据等值线成图方法，把不同测点不同时间段的感应电动势数据反演为空间离散视电阻率数据，再利用不规则三角网绘制等值线，通过视电阻率比较与颜色查找方法填充颜色，提高了瞬变电磁法探测数据绘制为等值线图的精度与速度。
需要说明的是，在本文中，诸如“第一”和“第二”等之类的关系术语仅仅用来将一个实体或者操作与另一个实体或操作区分开来，而不一定要求或者暗示这些实体或操作之间存在任何这种实际的关系或者顺序。而且，术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、物品或者设备所固有的要素。在没有更多限制的情况下，由语句“包括一个……”限定的要素，并不排除在包括所述要素的过程、方法、物品或者设备中还存在另外的相同要素。
以上所述仅是本发明的具体实施方式，使本领域技术人员能够理解或实现本发明。对这些实施例的多种修改对本领域的技术人员来说将是显而易见的，本文中所定义的一般原理可以在不脱离本发明的精神或范围的情况下，在其它实施例中实现。因此，本发明将不会被限制于本文所示的这些实施例，而是要符合与本文所公开的原理和新颖特点相一致的最宽的范围。
Claims (7)
1.一种瞬变电磁探测数据等值线成图与高阻区域识别方法，其特征在于，包括：
获取待绘制区域若干测点对应多个测道的感应电压数据；
将所述感应电压数据传输给数据处理单元，所述数据处理单元将所述感应电压数据反演成离散的视电阻率数据；
以测点号为空间横坐标、同时以测道号为视深坐标确定所述视电阻率数据的空间分布，生成离散数据点集；
根据所述离散数据点集生成不规则三角网TIN和有不同视电阻率的单等值线集数组Isoline[n]；
根据所述不规则三角网TIN，获得TIN'的边界点集ST；
根据所述单等值线集数组，获得等值线端点集SL；
将所述等值线端点集SL插入所述边界点集ST中，合并形成新的点集S；
对所述点集S中的边界点和等值线端点进行空间上的排序，形成有序点集S’；
根据所述单等值线集数组Isoline[n]，对应建立存储颜色值的颜色数组Color[n+1]，给出指定高阻区域视电阻率最小值D以及与所述视电阻率最小值D对应的颜色ColorS；
根据所述有序点集S’中的所有等值线端点，搜索待填充开区域Ri，并通过视电阻率值比较和颜色查找法确定所述填充开区域Ri的填充色，其中，当所述填充开区域Ri的视电阻率大于或等于D时，则所述填充开区域Ri的填充色为所述颜色ColorS；
判断所述有序点集S’中是否存在封闭等值线；
如果所述有序点集S’中存在封闭等值线，则通过视电阻率值比较和颜色查找法确定各所述封闭等值线所围成封闭区域的填充色，其中，当所述封闭区域的视电阻率大于或等于D时，则所述封闭区域的填充色为所述颜色ColorS；
如果所述有序点集S’中不存在封闭等值线，则成图结束。
2.根据权利要求1所述的瞬变电磁探测数据等值线成图与高阻区域识别方法，其特征在于，所述根据所述有序点集S’中的所有等值线线端点，搜索待填充开区域Ri，并通过视电阻率值比较和颜色查找法确定所述填充开区域Ri的填充色，其中，当所述填充开区域Ri的视电阻率大于或等于D时，则所述填充开区域Ri的填充色为所述颜色ColorS，包括：
依次序获取所述有序点集S’中新等值线端点；
判断所获取的新等值线端点是否为已获取的重复等值线端点；
如果所获取的新等值线端点为已获取的重复等值线端点，则结束开区域填充；
如果所获取的新等值线端点不是已获取的重复等值线端点，则以所述新等值线端点为起点搜索一个待填充开区域Ri；
获取在所述有序点集S’中、与所述新等值线相邻且排在所述新等值线端点前面的一个有序点，其中，所述有序点为等值线端点或边界点；
确定所述新等值线端点所在单等值线集的索引index；
判断所述新等值线端点的视电阻率是否大于所述有序点的视电阻率；
如果所述新等值线端点的视电阻率大于所述有序点的视电阻率，则用index查找所述颜色数组来确定所述开区域Ri的颜色；
如果所述新等值线端点的视电阻率小于或等于所述有序点的视电阻率，则用index+1查找所述颜色数组来确定所述开区域Ri的颜色；
判断所述开区域Ri对应的视电阻率是否大于或等于D；
如果所述开区域Ri对应的视电阻率大于或等于D，则所述开区域Ri的填充色为所述颜色ColorS。
3.根据权利要求1所述的瞬变电磁探测数据等值线成图与高阻区域识别方法，其特征在于，所述如果所述有序点集S’中存在封闭等值线，则通过视电阻率值比较和颜色查找法确定各所述封闭等值线所围成封闭区域的填充色，其中，当所述封闭区域的视电阻率大于或等于D时，则所述封闭区域的填充色为所述颜色ColorS，包括：
把所有封闭等值线按照在所述空间横坐标中的位置从小到大依次排序，并存入队列Q；
依次序获取所述队列Q中新封闭等值线；
判断所获取的新封闭等值线是否为已获取的重复封闭等值线；
如果所获取的新封闭等值线为已获取的重复封闭等值线，则结束封闭等值线填充；
如果所获取的新封闭等值线不是已获取的重复封闭等值线，则取所述新封闭等值线的线头点所在TIN线段的两个端点T1和T2；
取T1和T2中处于所述新封闭等值线所围成区域内部的点T，确定所述线头点所在单等值线集的索引CIndex；
判断所述线头点的视电阻率是否大于所述点T的视电阻率；
如果所述线头点的视电阻率大于所述点T的视电阻率，则用CIndex查找所述颜色数组来确定所述新封闭等值线所围成区域的颜色；
如果所述线头点的视电阻率小于或等于所述点T的视电阻率，则用CIndex+1查找所述颜色数组来确定所述新封闭等值线所围成区域的颜色；
判断所述新封闭等值线所围成区域对应的视电阻率是否大于或等于D；
如果所述新封闭等值线所围成区域对应的视电阻率大于或等于D，则所述新封闭等值线所围成区域的填充色为所述颜色ColorS。
4.根据权利要求1所述的瞬变电磁探测数据等值线成图与高阻区域识别方法，其特征在于，所述对所述点集S中的边界点和等值线端点进行空间上的排序，包括：
对所述点集S中的边界点和等值线端点进行空间上的逆时针排序。
5.根据权利要求1所述的瞬变电磁探测数据等值线成图与高阻区域识别方法，其特征在于，所述根据所述单等值线集数组Isoline[n]，对应建立存储颜色值的颜色数组Color[n+1]，其中：所述单等值线集数组Isoline[n]和所述颜色数组Color[n+1]的对应关系包括、所述单等值线集数组中一条单等值线Isoline[i]两侧的区域颜色为Color[i]和Color[i+1]。
6.根据权利要求2所述的瞬变电磁探测数据等值线成图与高阻区域识别方法，其特征在于，所述依次序获取所述有序点集S’中新等值线端点，包括：
依逆时针次序获取所述有序点集S’中新等值线端点。
7.根据权利要求1-6任一所述的瞬变电磁探测数据等值线成图与高阻区域识别方法，其特征在于，所述成图结束之后，还包括：
把成图与高阻区识别结果显示在终端屏幕和/或保存到数据库中。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510916764.XA CN105572753B (zh) | 2015-12-11 | 2015-12-11 | 一种瞬变电磁探测数据等值线成图与高阻区域识别方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510916764.XA CN105572753B (zh) | 2015-12-11 | 2015-12-11 | 一种瞬变电磁探测数据等值线成图与高阻区域识别方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105572753A true CN105572753A (zh) | 2016-05-11 |
| CN105572753B CN105572753B (zh) | 2017-12-05 |
Family
ID=55883094
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510916764.XA Active CN105572753B (zh) | 2015-12-11 | 2015-12-11 | 一种瞬变电磁探测数据等值线成图与高阻区域识别方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105572753B (zh) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108415081A (zh) * | 2018-02-27 | 2018-08-17 | 山东大学 | 一种瞬变电磁探测陆地日遗化武的方法 |
| CN112859182A (zh) * | 2021-01-19 | 2021-05-28 | 中煤科工集团西安研究院有限公司 | 一种矿井音频电透视反演成像假异常压制方法 |
| CN117092415A (zh) * | 2023-10-18 | 2023-11-21 | 深圳市城市公共安全技术研究院有限公司 | 一种区域电磁环境监测方法、装置、设备及介质 |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101738638A (zh) * | 2008-11-24 | 2010-06-16 | 中国石油集团东方地球物理勘探有限责任公司 | 石油地震构造图等值线的填充方法 |
| CN102314711A (zh) * | 2010-07-01 | 2012-01-11 | 中国地质科学院矿产资源研究所 | 矿产资源评价信息的三维可视化方法及其装置 |
| CN102759755A (zh) * | 2011-04-29 | 2012-10-31 | 黄桂芝 | 采用旋转tin网和非剖面法直接制作平、立面图的地质勘查方法 |
| JP5243476B2 (ja) * | 2010-04-13 | 2013-07-24 | 大日本コンサルタント株式会社 | 空中電磁探査法による地盤の解析方法 |
| CN104166158A (zh) * | 2013-05-17 | 2014-11-26 | 中国石油天然气集团公司 | 一种地震数据非均匀间距等值线交互追踪方法 |
-
2015
- 2015-12-11 CN CN201510916764.XA patent/CN105572753B/zh active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101738638A (zh) * | 2008-11-24 | 2010-06-16 | 中国石油集团东方地球物理勘探有限责任公司 | 石油地震构造图等值线的填充方法 |
| JP5243476B2 (ja) * | 2010-04-13 | 2013-07-24 | 大日本コンサルタント株式会社 | 空中電磁探査法による地盤の解析方法 |
| CN102314711A (zh) * | 2010-07-01 | 2012-01-11 | 中国地质科学院矿产资源研究所 | 矿产资源评价信息的三维可视化方法及其装置 |
| CN102759755A (zh) * | 2011-04-29 | 2012-10-31 | 黄桂芝 | 采用旋转tin网和非剖面法直接制作平、立面图的地质勘查方法 |
| CN104166158A (zh) * | 2013-05-17 | 2014-11-26 | 中国石油天然气集团公司 | 一种地震数据非均匀间距等值线交互追踪方法 |
Non-Patent Citations (2)
| Title |
|---|
| 付晓东 等: ""不规则区域内基于D-TIN的等值线生成算法"", 《昆明理工大学学报(理工版)》 * |
| 何俊 等: ""基于Delaunay三角网的二维电磁态势可视化方法"", 《系统工程理论与实践》 * |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108415081A (zh) * | 2018-02-27 | 2018-08-17 | 山东大学 | 一种瞬变电磁探测陆地日遗化武的方法 |
| CN112859182A (zh) * | 2021-01-19 | 2021-05-28 | 中煤科工集团西安研究院有限公司 | 一种矿井音频电透视反演成像假异常压制方法 |
| CN112859182B (zh) * | 2021-01-19 | 2022-11-18 | 中煤科工集团西安研究院有限公司 | 一种矿井音频电透视反演成像假异常压制方法 |
| CN117092415A (zh) * | 2023-10-18 | 2023-11-21 | 深圳市城市公共安全技术研究院有限公司 | 一种区域电磁环境监测方法、装置、设备及介质 |
| CN117092415B (zh) * | 2023-10-18 | 2024-01-19 | 深圳市城市公共安全技术研究院有限公司 | 一种区域电磁环境监测方法、装置、设备及介质 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105572753B (zh) | 2017-12-05 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Oldenburg et al. | Geophysical inversion for mineral exploration: A decade of progress in theory and practice | |
| CN105572753A (zh) | 一种瞬变电磁探测数据等值线成图与高阻区域识别方法 | |
| CN107247926B (zh) | 一种人体检测方法及装置 | |
| CN104504754B (zh) | 一种油气储层多点统计建模的方法及装置 | |
| CN104809756B (zh) | 基于X-ray CT图像的沥青混合料空隙空间结构重构方法 | |
| Sun et al. | Underground power cable detection and inspection technology based on magnetic field sensing at ground surface level | |
| CN106702995A (zh) | 基于bim的岩土工程监测模型构建方法 | |
| CN104252549A (zh) | 一种基于克里金插值的分析布井方法 | |
| CN112231874A (zh) | 地下管线模型的建立方法、装置、计算机设备和存储介质 | |
| CN108873073B (zh) | 一种基于网络并行电法的三维跨孔电阻率层析成像方法 | |
| CN103868958A (zh) | 一种电阻抗断层成像植物根系构型原位观测方法 | |
| CN106019386A (zh) | 一种回线源内比值视电阻率测定方法 | |
| CN106772620B (zh) | 十字形三维电法勘探方法 | |
| CN107703542B (zh) | 一种逆断层地区的地震层位的确定方法和装置 | |
| CN108415081B (zh) | 一种瞬变电磁探测陆地日遗化武的方法 | |
| Mulkal et al. | Inverse distance weight spatial interpolation for topographic surface 3D modelling | |
| CN113406707A (zh) | 一种大地电磁多尺度、多时段探测方法 | |
| CN108875184A (zh) | 基于数字露头模型的页岩有机碳含量预测方法及装置 | |
| CN109884714A (zh) | 一种可控源电磁测量方法、装置及其存储介质 | |
| CN115238566A (zh) | 一种基于频域电磁响应的人工空洞识别方法 | |
| CN106772625B (zh) | 一字形三维电法勘探方法 | |
| CN107371171A (zh) | 一种区域聚合的方法及装置 | |
| CN107591080B (zh) | 一种道路压盖关系的绘制方法及装置 | |
| Bobrov et al. | Physical and mathematical modeling of transient electromagnetic soundings over salt-dome structures | |
| Cascone et al. | ACLAS: A new automatic method of defining potential field lineaments using coherency analysis |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20191115 Address after: Yunda economic and Technological Development Zone in Yunnan province Kunming city 650217 West Road No. 105 Patentee after: Electric Power Research Institute of Yunnan Power System Ltd Address before: Yunda economic and Technological Development Zone in Yunnan province Kunming city 650217 West Road No. 105 Co-patentee before: KUNMING LIGONG GUANGZHI DETECTION TECHNOLOGY CO., LTD. Patentee before: Electric Power Research Institute of Yunnan Power System Ltd |
|
| TR01 | Transfer of patent right |