CN104407386B - 一种电磁勘探数据采集和处理方法及系统 - Google Patents
一种电磁勘探数据采集和处理方法及系统 Download PDFInfo
- Publication number
- CN104407386B CN104407386B CN201410737744.1A CN201410737744A CN104407386B CN 104407386 B CN104407386 B CN 104407386B CN 201410737744 A CN201410737744 A CN 201410737744A CN 104407386 B CN104407386 B CN 104407386B
- Authority
- CN
- China
- Prior art keywords
- survey
- data
- line
- electromagnetic
- survey line
- 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
Links
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
本申请提供一种电磁勘探数据采集和处理方法及系统。所述方法包括S1:在第一条测线上设置待测点A1、A2、…、AM;在第二条测线上设置待测点B1、B2、…、BN;S2:间隔选取第一条测线上以及第二条测线上的待测点为实际测点;S3:将所述实际测点按照A1、B2、A3、B4、…A(M‑1)、BN的顺序依次首尾连接,形成链状测线;在链状测线上的每个实际测点布置数据采集装置,获取所述链状测线的第一电磁勘探数据;S4:按照S1‑S3的步骤获取由其他相邻测线形成的链状测线上的第一电磁勘探数据;S5:基于第一电磁勘探数据计算所有测线上全部测点的第二电磁勘探数据。利用本申请中各个实施例,可以提高作业施工效率,节约成本,提高采集数据的分辨率。
Description
技术领域
本申请属于地球物理勘探技术领域,尤其涉及一种电磁勘探数据采集和处理方法及系统。
背景技术
随着人类日益增长的需求和不可再生资源的日益耗尽,当代社会正面临着资源紧缺的世界难题。如何更加有效、快速、二次深度的寻找以及开发现有资源越来越成为人们关注的话题。
目前,电磁法在地球物理勘探,尤其是油气资源、金属矿物资源勘探以及工程地质测量中有着显著而广泛的应用。该方法主要是通过在一定的测区范围内测量电磁场分量(场源可以是天然源也可以是人工源,可以是时间域也可以是频率域测量),通过计算阻抗和相位等,达到探测目标体的地质结构。现有技术中利用电磁勘探方法在野外采集数据时主要是采用多条测线并行排列方式布局测点,然后在所述测点上分别布设沿测线方向的电场分类Ex、垂直测线方向的电场分量Ey、垂直测线方向的磁场分量Hz等数据接收装置,所述数据接收装置将采集接收的电、磁分量通过数据线传输至数据采集站主机。作业人员可以利用采集的数据分析、获取地质构造。
在一次野外数据采集时,所述测线布置的条数通常设置为8条、16条、或者24条等。各条测线通常为平行设置,每条测线上可以等距离设置测点。上述电磁法在实际现场布置测点时通常受勘测区域地形和地质的影响导致测点设置困难,过多的测点导致测试设备相应的增多、设备损耗增加,增加了野外数据采集成本。由于上述常规数据采集通常一次只能采集单条测线上的数据,增加数据采集成本,同时由于电磁法勘探的体积效应,受旁侧地质体的影响等,采集的数据往往会存在异常干扰,降低了采集数据的分辨率。
发明内容
本申请目的在于提供一种电磁勘探数据采集和处理方法及系统,可以降低野外数据采集成本、提高施工效率、提高采集数据的分辨率。
本申请提供的一种电磁勘探数据采集和处理方法及系统是这样实现的:
一种电磁勘探数据采集和处理方法,其特征在于,所述方法包括:
S1:在勘探区域中P条测线的第一条测线上依次设置待测点A1、A2、A3、…、AM,M为所述第一条测线上待测点的个数;在与所述第一条测线相邻的第二条测线上依次设置待测点B1、B2、B3、…、BN,N为所述第二条测线上待测点的个数,M、N为自然数;
S2:选取所述第一条测线上的A1、A3、A5、…、A(M-1)待测点以及选取所述第二条测线上的B2、B4、B6、…、BN待测点为实际测点;
S3:将所述第一条测线与所述第二条测线上的实际测点按照A1、B2、A3、B4、…A(M-1)、BN的顺序依次首尾连接,形成链状测线;
在所述链状测线上的实际测点布置数据采集装置,获取所述链状测线的第一电磁勘探数据;
S4:按照S1-S3的步骤获取由其他相邻测线形成的链状测线上的第一电磁勘探数据;
S5:基于属于同一链状测线上相邻实际测点AK、B(K+1)的第一电磁勘探数据的电分量和磁分量计算得出所述同一链状测线上BK、A(K+1)待测点的第二电磁勘探数据;完成计算所述P条测线上所有待测点上的第二电磁勘探数据;将所述第一电磁勘探数据与所述第二电磁勘探数据作为所述勘探区域的电磁勘探数据,K为链状测线上实际测点的编号。
优选的实施例中,所述方法还包括:
S6:将属于同一链状测线上相邻实际测点的第一电磁勘探数据的平均值作为所述两个相邻实际测点连线中心点的组合电磁勘探数据;
相应的,所述勘探区域的电磁勘探数据还包括所述组合电磁勘探数据。
进一步的,所述方法还包括:
S7:基于所述采集和处理的电磁勘探数据获取所述勘探区域的地质构造结果。
一种电磁勘探测点的布设方法,所述方法包括:
S11:选取勘探区域中P条测线的第i条测线上依次设置待测点A1、A2、A3、…、AM,M为所述第i条测线上待测点的个数;在与所述第i条测线相邻的第i+1条测线上依次设置待测点B1、B2、B3、…、BN,N为所述第i+1条测线上待测点的个数,M、N为自然数,1≤i≤P;
S22:选取所述第i条测线上的A1、A3、A5、…、A(M-1)待测点以及选取所述第i+1条测线上的B2、B4、B6、…、BN待测点为实际测点;
将所述第i条测线与所述第i+1条测线上的实际测点按照A1、B2、A3、B4、…A(M-1)、BN的顺序依次首尾连接,形成链状测线,并在所述链状测线的实际测点上布置数据采集装置;
S33:按照上述S11、S22的步骤形成第i+2与i+3条测线、第i+4与i+5条测线、…、第P-1与P条测线的链状测线的实际测点上布设数据采集装置。
一种电磁勘探数据处理方法,所述方法包括:
S1’:获取通过链状阵列方法采集的第一电磁勘探数据;
S2’:基于所述第一电磁勘探数据中同一链状测线上相邻实际测点AK、B(K+1)的第一电磁勘探数据的电分量和磁分量计算得出勘探区域所述同一链状测线上BK、A(K+1)待测点的第二电磁勘探数据,A为链状测线上的其中一条测线,B为同一所述链状测线上的另一条测线,K为链状测线上实际测点的编号;
S3’:将所述第一电磁勘探数据与所述第二电磁勘探数据作为所述勘探区域的电磁勘探数据。
优选的实施方式,所述方法还包括:
S4’:将属于同一链状测线上相邻实际测点的第一电磁勘探数据的平均值作为所述两个相邻实际测点连线中心点的组合电磁勘探数据;
相应的,所述勘探区域的电磁勘探数据还包括所述组合电磁勘探数据。
一种电磁勘探数据处理系统,所述系统包括:
数据获取模块,用于获取通过链状阵列方法采集的第一电磁勘探数据;
待测点计算模块,用于基于所述第一电磁勘探数据中同一链状测线上相邻实际测点AK、B(K+1)的第一电磁勘探数据的电分量和磁分量计算得出勘探区域所述同一链状测线上BK、A(K+1)待测点的第二电磁勘探数据,A为链状测线上的其中一条测线,B为同一所述链状测线上的另一条测线,K为链状测线上实际测点的编号;
合并模块,用于将述第一电磁勘探数据与所述第二电磁勘探数据合并后作为所述勘探区域的电磁勘探数据。
优选的实施例中,所述系统还包括:
组合数据模块,用于将属于同一链状测线上相邻实际测点的第一电磁勘探数据的平均值作为所述两个相邻实际测点连线中心点的组合电磁勘探数据;
相应的,所述合并模块合并的勘探区域的电磁勘探数据还包括所述组合电磁勘探数据。
本申请提供一种电磁勘探数据采集和处理方法及系统,采用相邻测线上的测点相间隔分布,通过数据线将其连接形成链状测线的数据采集方式。这样可以实现原本一条测线上的数据采集装备获取两条测线上所有测点的数据,可以提高野外数据采集作业施工效率,减少设备损耗,节约数据采集成本。同时,由于采用非直线的链式采集方法,对于在因地形复杂不能布设数据采集装置的测点也可以获取该测点有效的测试数据,使得该方法适用于复杂地形等多种勘探区域。
附图说明
为了更清楚地说明本申请实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本申请中记载的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1是本申请一种电磁勘探数据采集和处理方法一种实施例的方法流程示意图;
图2是本申请勘探区域测线上待测点与实际测点分布示意图;
图3是本申请一种电磁勘探数据采集和处理方法另一个实施例的方法流程示意图;
图4是本申请一个实施例中包括组合电磁勘探数据的测点分布图;
图5是本申请应用实例勘探区域采用传统直线排列采集方法获取的数据采样点的分布图;
图6为勘探区域采用本申请所述方法获取的数据采样点的分布图;
图7是利用常规采集方法所获取的数据的勘探区域电阻率分布对比图;
图8是利用本申请所述方法获取的数据的勘探区域电阻率分布图;
图9为利用本申请所述方法获取的数据的时频3D切片与地震处理结果对比图;
图10是本申请所述一种电磁勘探数据处理系统一种实施例的模块结构示意图;
图11是本申请所述一种电磁勘探数据处理系统另一种实施例的模块结构示意图;
图12是本申请所述一种电磁勘探数据处理系统另一种实施例的模块结构示意图。
具体实施方式
为了使本技术领域的人员更好地理解本申请中的技术方案,下面将结合本申请实施例中的附图,对本申请实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本申请一部分实施例,而不是全部的实施例。基于本申请中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都应当属于本申请保护的范围。
本申请中所述利用电磁法进行勘探的基本原理为:当地下存在导电地质体时,在交变电磁场(一次场)的作用下,导体中将产生涡流(感应电流),涡流又在其周围产生二次磁场(二次场)。二次场的出现使一次场发生畸变。一般说,一次场和二次场迭加后的总场在强度、相位和方向上与一次场不同。基于麦克斯韦方程组的电磁理论基础,通过分析、研究二次场的强度和随时间衰变或研究总场各分量的强度、空间分布和时间特性等,可推断地下导电体的存在以及分析出导电体的属性、特性等。如果地质体具有高导磁性,在人工接地场源一次场作用下,受人工磁化产生二次磁场,同样可以发现异常并推断地下导磁体的存在。电磁法主要用于油气勘探、金属矿勘探、水文地质调查等。
在人工源电磁勘探进行数据采集的过程中,通常需要设置人工接地场源,用于向勘探区域发射不同频率的交变电流,然后在勘探区域预先设计好的用于数据采集的测点上布设Ex、Ey分量的电偶极观测装置和垂直向下方向Hz的磁场观测装置。通常在数据采集时,可以先建立数据采集坐标系,可以用于布置一条或多条测线以及在测线上设置相应的测点,以及方便后续数据提取、分析和处理。通常以平行于测线的方向为Ex方向,垂直于测线的方向为Ey方向,还可以观测垂直磁场Hz。进一步的,通过分析处理采集的勘探数据获知测试区域的地质构造。
本申请提供一种电磁勘探数据采集和处理方法,从对采集数据的测点的布置以及采集数据的进一步处理方法,不仅可以低野外数据采集成本、提高施工效率,还可以获取可靠性更高的电磁勘探数据,提高采集数据的分辨率。图1是本申请所述的一种电磁勘探数据采集和处理方法一种实施例的方法流程示意图,如图1所示,所述方法可以包括:
S1:在勘探区域中P条测线的第一条测线上依次设置待测点A1、A2、A3、…、AM,M为所述第一条测线上待测点的个数;
在与所述第一条测线相邻的第二条测线上依次设置待测点B1、B2、B3、…、BN,N为所述第二条测线上待测点的个数,M、N为自然数。
可以在需要勘探的区域中预先规划设计布置P条测线,所述布置的P条测线可根据预先勘测的地形建立的坐标系和测试精度等需求进行规划,包括测线条数的选择、每条测线之间的距离以及每条测线上待测点间距的设置等。所述勘探区域的P条测线之间通常设置为互相平行,以便于数据采集和分析。当然,本申请并不限定于所述每条测线之间必须平行,具体的测线排向方向的设置可以根据实际勘探区域的地理情况和测试需求进行设置。
本申请中所述的测线条数P,可以根据实际的勘探区域地形及测试需求进行选取,通常对于复杂勘探区域,所述布置的测线条数可以取较大的值,条与条之间的距离以及所述待测点间距等可以设置较小的值。本实施例中所述的测线条数P可以设置为8-24之间,通常的P取偶数。
所述方法中可以选取所述P条测线上的某一测线作为第一条测线,具体的可以选取靠近所述勘探区域边界的测线为第一条测线。在所述第一条测线上依次设置待测点A1、A2、A3、…、AM,M为所述第一条测线上待测点的个数,M为自然数。所述的待测点通常为预先规划的需要采集获取该处电磁勘探数据的测点。所述每一条测线上待测点的个数可以根据测试需求进行设置,为方便数据采集和后期数据分析,通常在所述P条测线的每一条测线上设置的待测点个数相同,本实施例中所述M与所述N可以相等。本实施例中相邻两条测线上的待测点个数相同即可组成完整的链状测线,如果不同的话,在末端稍微改变一下排列布置同样可以组成链状测线。
然后可以选取与所述第一条测线相邻的测线作为第二条测线,按照上述所述的方法在与所述第一条测线相邻的第二条测线上依次设置待测点B1、B2、B3、…、BN,N为所述第二条测线上待测点的个数,N为自然数。
S2:选取所述第一条测线上的A1、A3、…、A(M-1)待测点以及选取所述第二条测线上的B2、B4、…、BN待测点为实际测点;
在本实施例中,可以每间隔一个所述待测点选取所述第一条测线上的A1、A3、A5、…、A(M-1)待测点为实际测点,每间隔一个所述待测点选取所述第二条测线上的B2、B4、B6、…、BN待测点为实际测点。本实施例中所述的电磁勘探数据采集方法在每一条测线上间隔选取待测点作为数据采集的实际测点,可以减少使用数据采集设备的使用数目和设备安装时间,提高作业施工的效率和减少设备损耗。
S3:将所述第一条测线与所述第二条测线上的实际测点按照A1、B2、A3、B4、A5、B6、…、A(M-1)、BN的顺序依次首尾连接,形成链状测线;
在所述链状测线的实际测点上布置数据采集装置,获取所述链状测线的第一电磁勘探数据。
图2是勘探区域测线上待测点与实际测点分布示意图。如图2所示,A1、B2、A3、B4…AM为实际测点,其余为预先设计的待测点。可以将所述第一条测线上的实际测点与所述第二条测线上的在所述第一链状测线的实际测点依次首尾相连接,例如图2中所述可以将A1、B2、A3、B4、…A(M-1)、BN依次相连,形成链状测线的测点阵列数据观测、采集方式。
上述在所述链状测线的实际测点及链状测线上布置数据采集装置具体的可以包括在所述链状测线的实际测点上布设Ex、Ey分量的电偶极观测装置和垂直向下方向Hz的磁场观测装置,获取所述链状测线上的第一电磁勘探数据。上述中所述磁场观测装置的位置可以选取所述链状测线上干扰最小的测点的位置。所述的第一电磁勘探数据可以包括根据在实际测点采集的电分量、磁分量数据计算得出的电磁勘探数据,具体的例如可以根据实际测点A1采集获取的电分量Ex(A-1)、Ey(A-1)和磁分量计算得出实际测点A1的第一电磁勘探数据。
S4:按照S1-S3的步骤获取由其他相邻测线形成的链状测线上的第一电磁勘探数据。
在上述获取所述第一条测线和第二条测线形成的第一链状测线上的第一电磁勘探数据后,可以按照上述S1-S3的步骤获取所述N条测线中第三条测线和第四条测线、第五条测线和第六到测线、…、第N-1条测线和第N条测线的形成的链状测线上的第一电磁勘探数据。
当然,上述所述获取第一勘探数据的步骤也可以先将每两条测线设置成所述链状测线,在一次人工接地场源激发中同时获取所述勘探区域内所有链状测线上实际测点的第一电磁勘探数据。
S5:基于属于同一链状测线上相邻实际测点AK、B(K+1)的第一电磁勘探数据的电分量和磁分量计算得出所述同一链状测线上BK、A(K+1)待测点的第二电磁勘探数据,依次计算完成所述P条测线上所有待测点上的第二电磁勘探数据;
将所述第一电磁勘探数据与所述第二电磁勘探数据作为所述勘探区域的电磁勘探数据,K为链状测线上实际测点的编号。
如图2所示,在测量获取所述第一条测线上的A1、A3、…、A(M-1)实际测点以及所述第二条测线上的B2、B4、…、BN实际测点上的第一电磁勘探数据后,可以根据相邻实际测点的接收采集的数据计算获得测线上其他待测点的电磁勘探数据。例如可以根据实际测点A1的采集获取的电分量Ex(A-1)、实际测点B2的电分量Ey(B-2)和磁分量Hz计算得出待测点A2的第二电磁勘探数据,同时也可以根据实际测点A1的采集获取的电分量Ey(A-1)、实际测点B2的电分量Ex(B-2)和磁分量Hz计算得出待测点B1的第二电磁勘探数据。所述计算得出的第一电磁勘探数据或第二电磁勘探数据可以为矢量数据,可以表示在某一方向上测点的电磁勘探数值。
当然,参照上述方式,可以依次计算完成同一条测线上其他待测点如A4、B3等上的第二电磁勘探数据,以及进一步计算完所述勘探区域内所述N条测线上所有待测点上的第二电磁勘探数据。这样,可以将所述链状测线上实际测点上的第一电磁勘探数据与由所述第一电磁勘探数据计算获取的待测点上的第二电磁勘探数据作为所述勘探区域的电磁勘探数据,最终完成对所述勘探区域内的电磁勘探数据采集。
本实施例所述的一种电磁勘探数据采集方法,采用相邻测线上的测试相间隔分布,通过数据线将其连接形成链状测线的数据采集方式。这样可以实现原本一条测线上的数据采集装备获取两条测线上所有测点的数据,可以提高野外数据采集作业施工效率,减少设备损耗,节约数据采集成本。同时,由于采用非直线的链式采集方法,对于在因地形复杂不能布设数据采集装置的测点也可以获取该测点有效的测试数据,使得该方法适用于复杂地形等多种勘探区域。
优选的实施例中,利用上述所述的一种电磁勘探数据采集和处理方法获取的电磁勘探数据,还可以进一步的根据同一链式测线上的数据计算获得组合电磁勘探数据。相比于现有技术中的单测线的数据采集方式可以获取更多的有效电磁勘探数据。图3是本申请所述的一种电磁勘探数据采集和处理方法的另一种实施例的方法流程示意图,如图3所示,所述方法还可以包括:
S6:将属于同一链状测线上相邻实际测点AK、B(K+1)的第一电磁勘探数据的平均值作为所述两个相邻实际测点连线中心点的组合电磁勘探数据;
相应的,所述勘探区域的电磁勘探数据还可以包括所述组合电磁勘探数据。
图4是本实施例中所述计算得出的包括所述组合电磁勘探数据的测点分布图。如图4中所示,以同一链状测线上相邻的两个实际测点组合成一组新的电磁勘探数据,具体的实现方式可以包括通过将两个所述相邻实际测点的数据平均,如图4中实际测点A1与B2连线的中心点C1处的组合电磁勘探数据的电分量可以为:Ex(C1)={Ex(A-1)、Ex(B-2)}/2,Ey(C1)={Ey(A-1)、Ey(B-2)}/2,所述C1处的组合电磁勘探数据的磁分量可以共用实际测点A1或B2的磁分量。当然,在其他的实施方式中也可以将图4中的C2设置在实际测点B2与A3之间。
本实施例中所述的电磁勘探数据采集和处理方法,基于不同测线上测点的链式排列的地理位置关系,可以计算获取相邻测点的组合电磁勘探数据。该方式获取的组合电磁勘探数据不仅能有效压制单点实际测点的突发、随机干扰和旁侧效应的影响,提高电磁勘探数据的分辨率和可靠性,同时由于可以获取更多有效的采集数据,提高了数据解释精度,为后期数据处理提供了方便。基于此,本申请进一步的实施例中,所述方法还可以包括:
S7:基于所述采集和处理的电磁勘探数据获取所述勘探区域的地质构造结果。
尤其是对3D电磁勘探数据的处理,根据上述的电磁勘探数据分析获取的地质构造结果的效果大大提升。
下面是利用本申请上述实施例所述方法进行的电磁勘探的一个应用实例。
辽河盆地某区油气开发进入多轮次蒸汽吞吐末期,形成蒸汽空腔,储层复杂化,给开采部署带来难度。为研究汽腔发育和扩展状况、油水分布、剩余油气空间分布。采用本申请所述的链状测线方式进行了时频电磁勘探数据采集和采集后的数据处理。
图5为勘探区域采用传统直线排列采集方法获取的数据采样点的分布图。在图5的作业施工中。勘探区域面积为7km2范围,测网线点距75×50m,通过430个观测阵列,共完成测线布设42条,完成坐标点1838个,施工期为45天(2011年02月22日-2011年04月07日)。
图6为勘探区域采用本申请所述方法获取的数据采样点的分布图。本次作业勘探区域中,测网线点距100×150m。通过链状测线数据采集及采集后的数据处理,获得430个观测阵列,完成测线布设21条、坐标点1838个,但同时可以多获取429个组合电磁勘探数据,施工期减少为25天。由图5与图6的比较可以看出,相比以往使用的常规采集方法本申请所述的数据采集和处理方法获取电磁勘探数据的施工和采集效率大大提高。
图7是利用常规采集方法获取的数据的勘探区域电阻率分布图。图8是利用本申请所述方法获取的数据的勘探区域电阻率分布图。由图7、图8的对比中可以看出,相比常规数据采集方法所取得的剖面资料,本申请获取的资料随机干扰及静态位移得到明显改善。同时有效消除了由于干扰造成的虚假异常,所反映的浅层和深层电性特征更为符合地质规律,也更加清晰、可信。
图9为利用本申请所述方法获取的数据的时频3D切片与地震处理结果对比图。由图9可以看出,时频电磁3D切片所反映的电性异常分布规律与综合其他资料分析得出的地震勘探结果特征吻合度很高,说明利用本申请所述的数据采集和处理方法所取得的资料是可靠的,或者是加可靠的。
本申请还提供一种电磁勘探测点的布设方法,所述方法可以包括:
S11:选取勘探区域中P条测线的第i条测线上依次设置待测点A1、A2、A3、…、AM,M为所述第i条测线上待测点的个数;在与所述第i条测线相邻的第i+1条测线上依次设置待测点B1、B2、B3、…、BN,N为所述第i+1条测线上待测点的个数,M、N为自然数,1≤i≤P;
S22:选取所述第i条测线上的A1、A3、A5、…、A(M-1)待测点以及选取所述第二条测线上的B2、B4、B6、…、BN待测点为实际测点;
将所述第i条测线与所述第i+1条测线上的实际测点按照A1、B2、A3、B4、…A(M-1)、BN的顺序依次首尾连接,形成链状测线,并在所述链状测线的实际测点上布置数据采集装置;
S33:按照上述S11、S22的步骤形成第i+2与i+3条测线、第i+4与i+5条测线、…、第P-1与P条测线的链状测线的实际测点上布设数据采集装置。
本申请提供的一种电磁勘探测点的布设方法,在保障采集的电磁数据有效可靠的同时,可以大幅度减少数据采集装置的使用,减少设备的损害,节约施工成本,提高施工效率。
本申请中可以将上述电磁勘探测点的布设方法称为链状阵列方法,进一步的,本申请提供一种电磁勘探数据处理方法,所述方法可以包括:
S1’:获取通过链状阵列方法采集的第一电磁勘探数据;
S2’:基于所述第一电磁勘探数据中同一链状测线上相邻实际测点AK、B(K+1)的第一电磁勘探数据的电分量和磁分量计算得出勘探区域所述同一链状测线上BK、A(K+1)待测点的第二电磁勘探数据,完成计算得出勘探区域中所有测线上的待测点上的第二电磁勘探数据,A为链状测线上的其中一条测线,B为同一所述链状测线上的另一条测线,K为链状测线上实际测点的编号;
S3’:将所述第一电磁勘探数据与所述第二电磁勘探数据作为所述勘探区域的电磁勘探数据。
本申请所述的一种电磁勘探数据处理方法,基于上述各个实施例数据采集方法获取第一电磁勘探数据进一步计算获取勘探区域上所有测线上的待测点的第二电磁勘探数据,可以压制单测点随机干扰,最后将所述第一电磁勘探数据与所述第二电磁勘探数据作为所述勘探区域的电磁勘探数据。
参照本申请其他实施例,所述电磁勘探数据处理方法还可以包括:
S4’:将属于同一链状测线上相邻实际测点的第一电磁勘探数据的平均值作为所述两个相邻实际测点连线中心点的组合电磁勘探数据;
相应的,所述勘探区域的电磁勘探数据还包括所述组合电磁勘探数据。
在计算获取所述勘探区域的电磁勘探数据后,可以利用所述电磁勘探数据进行地质构造分析,因此上述所述方法还可以包括:
S5’:基于所述采集和处理的电磁勘探数据获取所述勘探区域的地质构造结果。
根据本申请所述的电磁勘探数据采集和处理方法,本申请还提供一种电磁勘探数据处理系统,图10是本申请所述一种电磁勘探数据处理系统的模块结构示意图,如图10所示,所述系统可以包括:
数据获取模块101,可以用于获取通过链状阵列方法采集的第一电磁勘探数据;
待测点计算模块102,可以用于基于所述第一电磁勘探数据中同一链状测线上相邻实际测点AK、B(K+1)的第一电磁勘探数据的电分量和磁分量计算得出勘探区域所述同一链状测线上BK、A(K+1)待测点的第二电磁勘探数据,A为链状测线上的其中一条测线,B为同一所述链状测线上的另一条测线,K为链状测线上实际测点的编号;
合并模块103,可以用于将述第一电磁勘探数据与所述第二电磁勘探数据合并后作为所述勘探区域的电磁勘探数据。
图11是本申请所述一种电磁勘探数据处理系统另一种实施例的模块结构示意图,如图11所示,所述一种电磁勘探数据处理系统还可以包括:
组合数据模块104,可以用于将属于同一链状测线上相邻实际测点的第一电磁勘探数据的平均值作为所述两个相邻实际测点连线中心点的组合电磁勘探数据;
相应的,所述合并模块103合并的勘探区域的电磁勘探数据还包括所述组合电磁勘探数据。
优选的实施例中,所述一种电磁勘探数据处理系统还可以包括:
地质分析模块105,可以用于基于所述勘探区域的电磁勘探数据获取所述勘探区域的地质构造结果。
图12是本申请所述一种电磁勘探数据处理系统另一种实施例的模块结构示意图。
本申请提供的一种电磁勘探数据处理系统,可以处理本申请所述电磁勘探数据采集和处理方法所述的获取的电磁勘探数据,进而可以进一步基于处理的电磁勘探数据分析勘探区域的地质构造。本申请提供的一种电磁勘探数据处理系统可以提高利用本申请所述电磁勘探数据采集和处理方法所述的获取的电磁勘探数据的分辨率,压制单测点随机干扰,使地质分析结果更加真实、可靠。
本说明书中的各个实施例采用递进的方式描述,各个实施例之间相同或相似的部分互相参见即可,每个实施例重点说明的都是与其他实施例的不同之处。本申请所述的数据处理系统可用于众多通用或专用的计算机系统环境或配置中。例如:个人计算机、服务器计算机、手持设备或便携式设备、平板型设备、多处理器系统、基于微处理器的系统、可编程的电子设备、网络PC、小型计算机、大型计算机、包括以上任何系统或设备的分布式计算环境等等。
虽然通过实施例描绘了本申请,本领域普通技术人员知道,本申请有许多变形和变化而不脱离本申请的精神,希望所附的权利要求包括这些变形和变化而不脱离本申请的精神。
Claims (9)
1.一种电磁勘探数据采集和处理方法,其特征在于,所述方法包括:
S1:在勘探区域中P条测线的第一条测线上依次设置待测点A1、A2、A3、…、AM,M为所述第一条测线上待测点的个数;在与所述第一条测线相邻的第二条测线上依次设置待测点B1、B2、B3、…、BN,N为所述第二条测线上待测点的个数,M、N为自然数;
S2:选取所述第一条测线上的A1、A3、A5、…、A(M-1)待测点以及选取所述第二条测线上的B2、B4、B6、…、BN待测点为实际测点;
S3:将所述第一条测线与所述第二条测线上的实际测点按照A1、B2、A3、B4、…A(M-1)、BN的顺序依次首尾连接,形成链状测线;
在所述链状测线上的实际测点上布置数据采集装置,获取所述链状测线的第一电磁勘探数据;
S4:按照S1-S3的步骤获取由其他相邻测线形成的链状测线上的第一电磁勘探数据;
S5:基于属于同一链状测线上相邻实际测点AK、B(K+1)的第一电磁勘探数据的电分量和磁分量计算得出所述同一链状测线上BK、A(K+1)待测点的第二电磁勘探数据;完成计算所述P条测线上所有待测点上的第二电磁勘探数据;将所述第一电磁勘探数据与所述第二电磁勘探数据作为所述勘探区域的电磁勘探数据,K为链状测线上实际测点的编号。
2.如权利要求1所述的一种电磁勘探数据采集和处理方法,其特征在于,所述方法还包括:
S6:将属于同一链状测线上相邻实际测点的第一电磁勘探数据的平均值作为两个相邻实际测点连线中心点的组合电磁勘探数据;
相应的,所述勘探区域的电磁勘探数据还包括所述组合电磁勘探数据。
3.如权利要求1或2所述的一种电磁勘探数据采集和处理方法,其特征在于,所述方法还包括:
S7:基于所述采集和处理的电磁勘探数据获取所述勘探区域的地质构造结果。
4.如权利要求1所述的一种电磁勘探数据采集和处理方法,其特征在于,在所述链状测线上的实际测点布置数据采集装置包括;
在所述链状测线的实际测点上布设电偶极观测装置和垂直向下方向的磁场观测装置;
上述磁场观测装置的位置为所述链状测线上干扰最小的测点的位置。
5.如权利要求1或2所述的一种电磁勘探数据采集和处理方法,其特征在于,所述的第一电磁勘探数据包括:根据在实际测点采集的电分量、磁分量数据计算得出的电磁勘探数据。
6.一种电磁勘探测点的布设方法,其特征在于,所述方法包括:
S11:选取勘探区域中P条测线的第i条测线上依次设置待测点A1、A2、A3、…、AM,M为所述第i条测线上待测点的个数;在与所述第i条测线相邻的第i+1条测线上依次设置待测点B1、B2、B3、…、BN,N为所述第i+1条测线上待测点的个数,M、N为自然数,1≤i≤P;
S22:选取所述第i条测线上的A1、A3、A5、…、A(M-1)待测点以及选取所述第i+1条测线上的B2、B4、B6、…、BN待测点为实际测点;
将所述第i条测线与所述第i+1条测线上的实际测点按照A1、B2、A3、B4、…A(M-1)、BN的顺序依次首尾连接,形成链状测线,并在所述链状测线的实际测点上布置数据采集装置;
S33:按照上述S11、S22的步骤形成第i+2与i+3条测线、第i+4与i+5条测线、…、第P-1与P条测线的链状测线的实际测点上布设数据采集装置。
7.一种电磁勘探数据处理系统,其特征在于,所述系统包括:
数据获取模块,用于获取通过链状阵列方法采集的第一电磁勘探数据;所述链状阵列包括根据权利要求1中S1至S4步骤所示的方法得到;
待测点计算模块,用于基于所述第一电磁勘探数据中同一链状测线上相邻实际测点AK、B(K+1)的第一电磁勘探数据的电分量和磁分量计算得出勘探区域所述同一链状测线上BK、A(K+1)待测点的第二电磁勘探数据,A为链状测线上的其中一条测线,B为同一所述链状测线上的另一条测线,K为链状测线上实际测点的编号;
合并模块,用于将述第一电磁勘探数据与所述第二电磁勘探数据合并后作为所述勘探区域的电磁勘探数据。
8.如权利要求7所述的一种电磁勘探数据处理系统,其特征在于,所述系统还包括:
组合数据模块,用于将属于同一链状测线上相邻实际测点的第一电磁勘探数据的平均值作为两个相邻实际测点连线中心点的组合电磁勘探数据;
相应的,所述合并模块合并的勘探区域的电磁勘探数据还包括所述组合电磁勘探数据。
9.如权利要求7或8所述的一种电磁勘探数据处理系统,其特征在于,所述系统还包括:
地质分析模块,用于基于所述勘探区域的电磁勘探数据获取所述勘探区域的地质构造结果。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410737744.1A CN104407386B (zh) | 2014-12-05 | 2014-12-05 | 一种电磁勘探数据采集和处理方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410737744.1A CN104407386B (zh) | 2014-12-05 | 2014-12-05 | 一种电磁勘探数据采集和处理方法及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104407386A CN104407386A (zh) | 2015-03-11 |
CN104407386B true CN104407386B (zh) | 2017-03-08 |
Family
ID=52645030
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410737744.1A Active CN104407386B (zh) | 2014-12-05 | 2014-12-05 | 一种电磁勘探数据采集和处理方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104407386B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105467460B (zh) * | 2015-12-28 | 2018-01-05 | 中国石油天然气集团公司 | 电磁勘探方法和装置 |
CN106814401A (zh) * | 2016-12-29 | 2017-06-09 | 深圳市赛盈地脉技术有限公司 | 地下资源测量方法及装置 |
CN109828239B (zh) * | 2019-02-26 | 2021-08-24 | 广西兴霖矿业投资咨询有限公司 | 一种用于矿产勘探的具有防护功能的矿井地质雷达仪 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102305948A (zh) * | 2011-05-25 | 2012-01-04 | 湖南继善高科技有限公司 | 测量地下电阻率立体变化的三维人工源电磁勘探方法 |
CN103839371A (zh) * | 2014-03-31 | 2014-06-04 | 内蒙古科技大学 | 一种室内无人值守安全巡检系统 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6566881B2 (en) * | 1999-12-01 | 2003-05-20 | Schlumberger Technology Corporation | Shielding method and apparatus using transverse slots |
EA019577B1 (ru) * | 2005-10-14 | 2014-04-30 | Эксонмобил Апстрим Рисерч Компани | Конструирование сигнала для геофизического исследования |
-
2014
- 2014-12-05 CN CN201410737744.1A patent/CN104407386B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102305948A (zh) * | 2011-05-25 | 2012-01-04 | 湖南继善高科技有限公司 | 测量地下电阻率立体变化的三维人工源电磁勘探方法 |
CN103839371A (zh) * | 2014-03-31 | 2014-06-04 | 内蒙古科技大学 | 一种室内无人值守安全巡检系统 |
Also Published As
Publication number | Publication date |
---|---|
CN104407386A (zh) | 2015-03-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Okpoli | Sensitivity and resolution capacity of electrode configurations | |
Wilson et al. | From outcrop to flow simulation: Constructing discrete fracture models from a LIDAR survey | |
Xue et al. | Developments measurements of TEM sounding in China | |
CN105911603A (zh) | 基于天然电场的四维物探方法 | |
CN104237970A (zh) | 地震电磁联合勘探系统及其数据采集装置和数据采集方法 | |
CN104614774B (zh) | 一种瞬变电磁探测方法、装置和系统 | |
CN104597519A (zh) | 一种深部金矿阶梯式找矿方法 | |
CN106896415A (zh) | 一种地‑井电性源短偏移距瞬变电磁探测方法 | |
CN106706715A (zh) | 基于三维高密度电阻率法的污染土检测方法 | |
CN105467460B (zh) | 电磁勘探方法和装置 | |
Legault et al. | Synthetic model testing and distributed acquisition dc resistivity results over an unconformity uranium target from the Athabasca Basin, northern Saskatchewan | |
CN109100808A (zh) | 一种多线源瞬变电磁横磁极化场探测方法 | |
CN101625413A (zh) | 一种油气勘查方法 | |
CN104407386B (zh) | 一种电磁勘探数据采集和处理方法及系统 | |
US20170249405A1 (en) | Method of intergrowing and coexisting four-in-one mineral coordinated prospecting | |
CN206378448U (zh) | 基于综合物探方法的污染土检测测线布置结构 | |
CN108873073B (zh) | 一种基于网络并行电法的三维跨孔电阻率层析成像方法 | |
CN105891896A (zh) | 一种地下空区的特征信息识别与分析方法 | |
CN105445787A (zh) | 一种最优方位子体相干的裂缝预测方法 | |
Peng et al. | Data field application in removing large P-phase arrival picking errors and relocating a mine microseismic event | |
Wang et al. | Dynamic monitoring of coalbed methane reservoirs using Super-Low Frequency electromagnetic prospecting | |
CN102608665A (zh) | 物探时间域瞬变电磁系统测量数据资料精细处理技术 | |
Xu et al. | Groundwater resources survey of Tongchuan city using the audio magnetotelluric method | |
Rees et al. | Magnetotelluric monitoring of coal seam gas depressurization | |
CN106772630A (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 |