CN114200526A - 电磁勘探方法及装置 - Google Patents
电磁勘探方法及装置 Download PDFInfo
- Publication number
- CN114200526A CN114200526A CN202010909504.0A CN202010909504A CN114200526A CN 114200526 A CN114200526 A CN 114200526A CN 202010909504 A CN202010909504 A CN 202010909504A CN 114200526 A CN114200526 A CN 114200526A
- Authority
- CN
- China
- Prior art keywords
- point
- wave field
- imaging
- measuring point
- electromotive force
- 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
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/08—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation operating with magnetic or electric fields produced or modified by objects or geological structures or by detecting devices
-
- 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
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/30—Assessment of water resources
Abstract
本公开公开了一种电磁勘探方法及装置,属于地球物理勘探技术领域。所述方法包括:获取每个测点对应的波场电动势集合;确定每个测点与成像点之间的距离;根据每个测点与成像点之间的距离及每个测点的波场电动势集合,确定每个测点在成像点的波场值;将各个测点在成像点的波场值进行叠加,得到成像点的总波场值;根据成像点的总波场值,进行电磁勘探。本公开不是基于传统的视电阻率概念而是基于电磁波在地下的传播速度远低于光速的特性设计而成,该方法无需进行方程变换,通过波场值叠加即可计算出每个成像点的总波场值,因而降低了测量过程的复杂度,且能有效降低电磁法的体积效应,提高勘探的效果。
Description
技术领域
本公开涉及地球物理勘探技术领域,特别涉及一种电磁勘探方法及装置。
背景技术
电磁勘探法是指以电磁场为工具,基于电磁场在地下的传播规律而进行勘探的方法,电磁勘探方法包括直流电阻率法、大地电磁测深法、电偶源瞬变电磁测深法、磁偶源瞬变电磁测深法、探地雷达法等多种方法。受制于体积效果,导致电磁法勘探方法的分辨率较低。然而,随着生产实践对精度的要求越来越高,如何提高电磁勘探方法的分辨率变得至关重要。
由于大地中的电磁场满足扩散方程而不是波段方程,因此,相关技术在进行电磁勘探时,需要进行波场变换,将电磁场的扩散方程变换成波段方程,并将波动方程进行逆时偏移成像,进而基于电磁偏移成像结果进行电磁勘探。
然而,相关技术需要进行复杂的方程变换,导致勘探过程复杂度较高。
发明内容
为了解决相关技术的问题,本公开实施例提供了一种电磁勘探方法及装置。所述技术方案如下:
一方面,提供了一种电磁勘探方法,所述方法包括:
获取每个测点对应的波场电动势集合,所述波场电动势集合包括不同时刻在每个测点采集到的多个电动势;
确定每个测点与成像点之间的距离;
根据每个测点与所述成像点之间的距离及每个测点的波场电动势集合,确定每个测点在所述成像点的波场值;
将各个测点在所述成像点的波场值进行叠加,得到所述成像点的总波场值;
根据所述成像点的总波场值,进行电磁勘探。
在本公开的另一个实施例中,所述确定每个测点与成像点之间的距离,包括:
获取电磁信号从每个测点到所述成像点的往返时间;
获取所述电磁信号的传播速度;
根据所述往返时间和所述传播速度,确定每个测点与所述成像点之间的距离。
在本公开的另一个实施例中,所述根据每个测点与所述成像点之间的距离及每个测点的波场电动势集合,确定每个测点在所述成像点的波场值,包括:
从每个测点的电动势集合中,获取目标电动势;
根据目标电动势及每个测点与所述成像点之间的距离,确定每个测点在所述成像点的波场值。
在本公开的另一个实施例中,所述从每个测点的电动势集合中,获取目标电动势,包括:
计算指定时刻与采集时间间隔的比值,得到所述指定时刻对应的采集编号;
基于所述采集编号,从所述电动势集合中,获取所述目标电动势。
在本公开的另一个实施例中,所述根据每个目标电动势及每个测点与所述成像点之间的距离,确定每个测点在所述成像点的波场值,包括:
对于任一测点,如果所述测点的目标电动势为感应电动势,则将所述目标电动势和所述测点与所述成像点之间的距离的乘积,作为所述测点在所述成像点的波场值;
如果所述目标电动势为电场或磁场的电动势,则获取所述测点与所述成像点之间的距离的平方,并将所述目标电动势与所述距离的平方的乘积,作为所述测点在所述成像点的波场值。
在本公开的另一个实施例中,所述根据所述成像点的总波场值,进行电磁勘探,包括:
根据所述成像点的位置坐标和所述成像点的总波场值,将总波场值相同的成像点相连,得到波场值剖面图;
基于波场值剖面图进行电磁勘探。
另一方面,提供了一种电磁勘探装置,所述装置包括:
获取模块,用于获取每个测点对应的波场电动势集合,所述波场电动势集合包括不同时刻在每个测点采集到的多个电动势;
第一确定模块,用于确定每个测点与成像点之间的距离;
第二确定模块,用于根据每个测点与所述成像点之间的距离及每个测点的波场电动势集合,确定每个测点在所述成像点的波场值;
叠加模块,用于将各个测点在所述成像点的波场值进行叠加,得到所述成像点的总波场值;
勘探模块,用于根据所述成像点的总波场值,进行电磁勘探。
在本公开的另一个实施例中,所述第一确定模块,用于获取电磁信号从每个测点到所述成像点的往返时间;获取所述电磁信号的传播速度;根据所述往返时间和所述传播速度,确定每个测点与所述成像点之间的距离。
在本公开的另一个实施例中,所述第二确定模块,用于从每个测点的电动势集合中,获取目标电动势;根据目标电动势及每个测点与所述成像点之间的距离,确定每个测点在所述成像点的波场值。
在本公开的另一个实施例中,所述第二确定模块,用于计算指定时刻与采集时间间隔的比值,得到所述指定时刻对应的采集编号;基于所述采集编号,从所述电动势集合中,获取所述目标电动势。
在本公开的另一个实施例中,所述第二确定模块,用于对于任一测点,如果所述测点的目标电动势为感应电动势,则将所述目标电动势和所述测点与所述成像点之间的距离的乘积,作为所述测点在所述成像点的波场值;如果所述目标电动势为电场或磁场的电动势,则获取所述测点与所述成像点之间的距离的平方,并将所述目标电动势与所述距离的平方的乘积,作为所述测点在所述成像点的波场值。
在本公开的另一个实施例中,所述勘探模块,用于根据所述成像点的位置坐标和所述成像点的总波场值,将总波场值相同的成像点相连,得到波场值剖面图;基于波场值剖面图进行电磁勘探。
本公开实施例提供的技术方案带来的有益效果是:
不是基于传统的视电阻率概念而是基于电磁波在地下的传播速度远低于光速的特性设计而成,该方法无需进行方程变换,通过波场值叠加即可计算出每个成像点的总波场值,因而降低了测量过程的复杂度,且能有效降低电磁法的体积效应,提高勘探的效果。
附图说明
为了更清楚地说明本公开实施例中的技术方案,下面将对实施例描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本公开的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1是本公开实施例提供的一种电磁勘探方法的流程图;
图2是本公开实施例提供的另一种电磁勘探方法的流程图;
图3是本公开实施例提供的观测点与成像点的示意图;
图4是本公开实施例提供的一种感应电动势的曲线图;
图5是本公开实施例提供的一种偏移成像结果图;
图6是本公开实施例提供的一种电磁勘探装置的结构示意图。
具体实施方式
为使本公开的目的、技术方案和优点更加清楚,下面将结合附图对本公开实施方式作进一步地详细描述。
本公开实施例提供了一种电磁勘探方法,参见图1,本公开实施例提供的方法流程包括:
101、获取每个测点对应的波场电动势集合。
其中,波场电动势集合包括不同时刻在每个测点采集到的多个电动势。
102、确定每个测点与成像点之间的距离。
103、根据每个测点与成像点之间的距离及每个测点的波场电动势集合,确定每个测点在成像点的波场值。
104、将各个测点在成像点的波场值进行叠加,得到成像点的总波场值。
105、根据成像点的总波场值,进行电磁勘探。
本公开实施例提供的方法,不是基于传统的视电阻率概念而是基于电磁波在地下的传播速度远低于光速的特性设计而成,该方法无需进行方程变换,通过波场值叠加即可计算出每个成像点的总波场值,因而降低了测量过程的复杂度,且能有效降低电磁法的体积效应,提高勘探的效果。
在本公开的另一个实施例中,确定每个测点与成像点之间的距离,包括:
获取电磁信号从每个测点到成像点的往返时间;
获取电磁信号的传播速度;
根据往返时间和传播速度,确定每个测点与成像点之间的距离。
在本公开的另一个实施例中,根据每个测点与成像点之间的距离及每个测点的波场电动势集合,确定每个测点在成像点的波场值,包括:
从每个测点的电动势集合中,获取目标电动势;
根据目标电动势及每个测点与成像点之间的距离,确定每个测点在成像点的波场值。
在本公开的另一个实施例中,从每个测点的电动势集合中,获取目标电动势,包括:
计算指定时刻与采集时间间隔的比值,得到指定时刻对应的采集编号;
基于采集编号,从电动势集合中,获取目标电动势。
在本公开的另一个实施例中,根据每个目标电动势及每个测点与成像点之间的距离,确定每个测点在成像点的波场值,包括:
对于任一测点,如果测点的目标电动势为感应电动势,则将目标电动势和测点与成像点之间的距离的乘积,作为测点在成像点的波场值;
如果目标电动势为电场或磁场的电动势,则获取测点与成像点之间的距离的平方,并将目标电动势与距离的平方的乘积,作为测点在成像点的波场值。
在本公开的另一个实施例中,根据成像点的总波场值,进行电磁勘探,包括:
根据成像点的位置坐标和成像点的总波场值,将总波场值相同的成像点相连,得到波场值剖面图;
基于波场值剖面图进行电磁勘探。
上述所有可选技术方案,可以采用任意结合形成本公开的可选实施例,在此不再一一赘述。
本公开实施例提供了一种电磁勘探方法,以进行电磁勘探的设备执行本公开实施例为例,参见图2,本公开实施例提供的方法流程包括:
201、获取每个测点的位置坐标及成像点的位置坐标。
其中,测点是指预先设置的对电磁信号进行接收、采集的观测点。测点的数量可根据实际勘探需要进行设置,如果对勘探精度要求较高,则测点的数量可以相对选取较多数量;如果对勘探精度要求较低,则测点的数量可以相对选取较少的数量。
根据测点的大地坐标计算相对坐标。相对坐标的原点设置在测线的第一个测点。在相对坐标系中,测点的横坐标表示测点与原点之间的距离,测点的纵坐标表示测点的高程,也即是测点的海拔高度。测点的横坐标和纵坐标的单位为米。在本公开实施例中,测点的横坐标采用u表示,纵坐标采用v表示,测点的位置坐标表示为(u,v)。对于第i个测点,其坐标表示为(u(i),v(i))。
其中,成像点为预先设置的待确定波场强度的点。成像点的数量也可根据实际勘探需求进行设置,如果对勘探精度要求较高,则成像点的数量可以相对选取较多数量;如果对勘探精度要求较低,则成像点的数量可以相对选取较少数量。
成像点的横坐标和纵坐标表示水平面上横纵方向的坐标。成像点的横坐标用x表示,成像点的纵坐标用y表示。考虑到不同海拔高度的成像点的波场值是不同的,对于海拔高度为k的第i成像点P的坐标表示为(x(i,k),y(i,k))。
图3示出了测点与成像点的示意图,图中1指示的箭头表示测点,图中2所示的所有网格节点表示成像点,由可3可知,测点数量为5个,成像点数量为15个。
202、获取每个测点对应的波场电动势集合。
对于每个测点,每隔预设时间间隔采集一次电动势,得到多个电动势,并将多个电动势组成波场电动势集合。其中,采集时间间隔可用delta T表示,采集到的多个电动势可用y1、y2、y3、…、yn表示,测点对应的波场电动势集合表示为(y1,y2,y3,…,yn)。
图4示出了中心回线瞬变电磁测深观测方式记录的一个测点的感应电动势曲线,从0时刻开始采样,采样的时间步长是0.8微秒,共记录0.024秒的时间,即30000个感应电动势数据。图4中横轴表示时间,时间单位为秒;纵轴表示感应电动势,感应电动势的单位是毫伏。
203、确定每个测点与成像点之间的距离。
本公开实施例在确定每个测点与成像点之间的距离时,可采用如下方法:
2031、获取电磁信号从每个测点到成像点的往返时间。
其中,往返时间是指电磁信号从测点发出后到达成像点,再从成像点回到测点的时间。
2032、获取电磁信号的传播速度。
考虑到不同频率的电磁波在不同电阻率的大地中传播时,电磁波的传播速度是不同的,因此,本公开实施例提供的方法在获取电磁信号的传播速度时,可获取电磁信号的电阻率,并获取大地的电阻率,进而根据电磁信号的电阻率和大地的电阻率,确定出电磁信号的传播速度。
表1示出了不同频率的电磁波在不同电阻率的大地传播时,电磁波的传播速度。
表1
2033、根据往返时间和传播速度,确定每个测点与成像点之间的距离。
对于任一测点,通过将往返时间与传播速度相乘,得到往返路程,并将往返路程除以2,得到该测点与成像点之间的距离。设定往返时间为T,传播速度为V,则往返路程S=TV,测点与成像点之间的距离R=S/2。
204、根据每个测点与成像点之间的距离及每个测点的波场电动势集合,确定每个测点在成像点的波场值。
本公开实施提供的方法在根据每个测点与成像点之间的距离及每个测点的波场电动势集合,确定每个测点在成像点的波场值时,可采用如下方法:
2041、从每个测点的电动势集合中,获取目标电动势。
由于本公开实施例中每个测点的电动势都对应采集时间,因而获取指定时刻的目标电动势时,可计算指定时刻与采集时间间隔的比值,得到指定时刻对应的采集编号,进而基于采集编号,从电动势集合中,获取目标电动势。其中,指定时刻可以为预先设置的勘探时间。
例如,设定指定时刻为T1,采集时间间隔为deltaT,则采集编号index=T/deltaT,基于该采集编号从测点A的电动势集合(y1,y2,y3,…,yn)中找到第i个数据yi,进而确定目标电动势为yi。
2042、根据目标电动势及每个测点与成像点之间的距离,确定每个测点在成像点的波场值。
由于测点的电动势包括感应电动势、波场或磁场电动势等,对于不同的电动势在计算成像点的波场值时,采取的计算方法是不同的,因此,对于任一测点,针对该测点的目标电动势的不同,在确定该测点在成像点的波场值时,包括但不限于如下方式:
第一种方式、如果测点的目标电动势为感应电动势,则将目标电动势和测点与成像点之间的距离的乘积,作为测点在成像点的波场值。
设定目标电动势为F,测点与成像点之间的距离为R,则当目标电动势为感应电动势时,测点在成像点的波场值为FP=F*R。
第二种方式、如果目标电动势为电场或磁场的电动势,则获取测点与成像点之间的距离的平方,并将目标电动势与距离的平方的乘积,作为测点在成像点的波场值。
设定目标电动势为F,测点与成像点之间的距离为R,则当目标电动势为电场或磁场的电动势时,测点在成像点的波场值为FP=F*R2。
205、将各个测点在成像点的波场值进行叠加,得到成像点的总波场值。
当获取到各个测点在成像点的波场值,可将各个测点在成像点的波场值进行叠加,得到该成像点的总波场值。
206、根据成像点的总波场值,进行电磁勘探。
当得到成像点的总波场值,本公开实施例根据成像点的位置坐标和成像点的总波场值,将总波场值相同的成像点相连,得到波场值剖面图,进而基于波场值剖面图进行电磁勘探。
具体实施时,可基于等值线绘制工具将总波场值相同位置坐标不同的成像点相连,得到波场值剖面图,该波场值剖面图可参见图5。其中,等值线绘制工具可以为surfer等。由于不同物质所形成的波场值剖面图是不同的,因此,基于不同物质的特性及所得到的波场值剖面图,可进行金属矿勘探、地下水勘探、油气资源勘探、工程检测等电磁勘探。
本公开实施例提供的方法,不是基于传统的视电阻率概念而是基于电磁波在地下的传播速度远低于光速的特性设计而成,该方法无需进行方程变换,通过波场值叠加即可计算出每个成像点的总波场值,因而降低了测量过程的复杂度,且能有效降低电磁法的体积效应,提高勘探的效果。
参见图6,本公开实施例提供了一种电磁勘探装置,该装置包括:
获取模块601,用于获取每个测点对应的波场电动势集合,波场电动势集合包括不同时刻在每个测点采集到的多个电动势;
第一确定模块602,用于确定每个测点与成像点之间的距离;
第二确定模块603,用于根据每个测点与成像点之间的距离及每个测点的波场电动势集合,确定每个测点在成像点的波场值;
叠加模块604,用于将各个测点在成像点的波场值进行叠加,得到成像点的总波场值;
勘探模块605,用于根据成像点的总波场值,进行电磁勘探。
在本公开的另一个实施例中,第一确定模块602,用于获取电磁信号从每个测点到成像点的往返时间;获取电磁信号的传播速度;根据往返时间和传播速度,确定每个测点与成像点之间的距离。
在本公开的另一个实施例中,第二确定模块603,用于从每个测点的电动势集合中,获取目标电动势;根据目标电动势及每个测点与成像点之间的距离,确定每个测点在成像点的波场值。
在本公开的另一个实施例中,第二确定模块603,用于计算指定时刻与采集时间间隔的比值,得到指定时刻对应的采集编号;基于采集编号,从电动势集合中,获取目标电动势。
在本公开的另一个实施例中,第二确定模块603,用于对于任一测点,如果测点的目标电动势为感应电动势,则将目标电动势和测点与成像点之间的距离的乘积,作为测点在成像点的波场值;如果目标电动势为电场或磁场的电动势,则获取测点与成像点之间的距离的平方,并将目标电动势与距离的平方的乘积,作为测点在成像点的波场值。
在本公开的另一个实施例中,勘探模块605,用于根据成像点的位置坐标和成像点的总波场值,将总波场值相同的成像点相连,得到波场值剖面图;基于波场值剖面图进行电磁勘探。
综上,本公开实施例提供的装置,不是基于传统的视电阻率概念而是基于电磁波在地下的传播速度远低于光速的特性设计而成,该方法无需进行方程变换,通过波场值叠加即可计算出每个成像点的总波场值,因而降低了测量过程的复杂度,且能有效降低电磁法的体积效应,提高勘探的效果。
本领域普通技术人员可以理解实现上述实施例的全部或部分步骤可以通过硬件来完成,也可以通过程序来指令相关的硬件完成,所述的程序可以存储于一种计算机可读存储介质中,上述提到的存储介质可以是只读存储器,磁盘或光盘等。
以上所述仅为本公开的可选实施例,并不用以限制本公开,凡在本公开的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本公开的保护范围之内。
Claims (10)
1.一种电磁勘探方法,其特征在于,所述方法包括:
获取每个测点对应的波场电动势集合,所述波场电动势集合包括不同时刻在每个测点采集到的多个电动势;
确定每个测点与成像点之间的距离;
根据每个测点与所述成像点之间的距离及每个测点的波场电动势集合,确定每个测点在所述成像点的波场值;
将各个测点在所述成像点的波场值进行叠加,得到所述成像点的总波场值;
根据所述成像点的总波场值,进行电磁勘探。
2.根据权利要求1所述的方法,其特征在于,所述确定每个测点与成像点之间的距离,包括:
获取电磁信号从每个测点到所述成像点的往返时间;
获取所述电磁信号的传播速度;
根据所述往返时间和所述传播速度,确定每个测点与所述成像点之间的距离。
3.根据权利要求1所述的方法,其特征在于,所述根据每个测点与所述成像点之间的距离及每个测点的波场电动势集合,确定每个测点在所述成像点的波场值,包括:
从每个测点的电动势集合中,获取目标电动势;
根据目标电动势及每个测点与所述成像点之间的距离,确定每个测点在所述成像点的波场值。
4.根据权利要求3所述的方法,其特征在于,所述从每个测点的电动势集合中,获取目标电动势,包括:
计算指定时刻与采集时间间隔的比值,得到所述指定时刻对应的采集编号;
基于所述采集编号,从所述电动势集合中,获取所述目标电动势。
5.根据权利要求3所述的方法,其特征在于,所述根据每个目标电动势及每个测点与所述成像点之间的距离,确定每个测点在所述成像点的波场值,包括:
对于任一测点,如果所述测点的目标电动势为感应电动势,则将所述目标电动势和所述测点与所述成像点之间的距离的乘积,作为所述测点在所述成像点的波场值;
如果所述目标电动势为电场或磁场的电动势,则获取所述测点与所述成像点之间的距离的平方,并将所述目标电动势与所述距离的平方的乘积,作为所述测点在所述成像点的波场值。
6.根据权利要求1所述的方法,其特征在于,所述根据所述成像点的总波场值,进行电磁勘探,包括:
根据所述成像点的位置坐标和所述成像点的总波场值,将总波场值相同的成像点相连,得到波场值剖面图;
基于波场值剖面图进行电磁勘探。
7.一种电磁勘探装置,其特征在于,所述装置包括:
获取模块,用于获取每个测点对应的波场电动势集合,所述波场电动势集合包括不同时刻在每个测点采集到的多个电动势;
第一确定模块,用于确定每个测点与成像点之间的距离;
第二确定模块,用于根据每个测点与所述成像点之间的距离及每个测点的波场电动势集合,确定每个测点在所述成像点的波场值;
叠加模块,用于将各个测点在所述成像点的波场值进行叠加,得到所述成像点的总波场值;
勘探模块,用于根据所述成像点的总波场值,进行电磁勘探。
8.根据权利要求7所述的装置,其特征在于,所述第一确定模块,用于获取电磁信号从每个测点到所述成像点的往返时间;获取所述电磁信号的传播速度;根据所述往返时间和所述传播速度,确定每个测点与所述成像点之间的距离。
9.根据权利要求7所述的装置,其特征在于,所述第二确定模块,用于从每个测点的电动势集合中,获取目标电动势;根据目标电动势及每个测点与所述成像点之间的距离,确定每个测点在所述成像点的波场值。
10.根据权利要求9所述的装置,其特征在于,所述第二确定模块,用于计算指定时刻与采集时间间隔的比值,得到所述指定时刻对应的采集编号;基于所述采集编号,从所述电动势集合中,获取所述目标电动势。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010909504.0A CN114200526B (zh) | 2020-09-02 | 2020-09-02 | 电磁勘探方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010909504.0A CN114200526B (zh) | 2020-09-02 | 2020-09-02 | 电磁勘探方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114200526A true CN114200526A (zh) | 2022-03-18 |
CN114200526B CN114200526B (zh) | 2023-08-22 |
Family
ID=80644315
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010909504.0A Active CN114200526B (zh) | 2020-09-02 | 2020-09-02 | 电磁勘探方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114200526B (zh) |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB1279105A (en) * | 1969-10-09 | 1972-06-28 | Chevron Res | Vertical resistivity logging |
US7203599B1 (en) * | 2006-01-30 | 2007-04-10 | Kjt Enterprises, Inc. | Method for acquiring transient electromagnetic survey data |
CN102419453A (zh) * | 2011-07-15 | 2012-04-18 | 中国科学院地质与地球物理研究所 | 长导线源瞬变电磁地空探测方法 |
CN103499842A (zh) * | 2013-09-27 | 2014-01-08 | 长江大学 | 一种微米电磁勘探方法 |
US20160313468A1 (en) * | 2015-04-27 | 2016-10-27 | Pgs Geophysical As | Seismic guided inversion of electromagnetic survey data |
CN106547031A (zh) * | 2016-10-09 | 2017-03-29 | 煤炭科学技术研究院有限公司 | 一种基于场变换确定地质电性界面的瞬变电磁探测法 |
CN106610509A (zh) * | 2015-10-26 | 2017-05-03 | 中国石油天然气股份有限公司 | 一种瞬变电磁资料的时域处理方法 |
CN109752762A (zh) * | 2018-12-11 | 2019-05-14 | 中煤科工集团西安研究院有限公司 | 单发多收观测装置瞬变电场数据动校正方法及装置 |
CN111123264A (zh) * | 2019-12-31 | 2020-05-08 | 桂林理工大学 | 一种考虑介质频散和衰减补偿的探地雷达逆时偏移成像方法 |
-
2020
- 2020-09-02 CN CN202010909504.0A patent/CN114200526B/zh active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB1279105A (en) * | 1969-10-09 | 1972-06-28 | Chevron Res | Vertical resistivity logging |
US7203599B1 (en) * | 2006-01-30 | 2007-04-10 | Kjt Enterprises, Inc. | Method for acquiring transient electromagnetic survey data |
CN102419453A (zh) * | 2011-07-15 | 2012-04-18 | 中国科学院地质与地球物理研究所 | 长导线源瞬变电磁地空探测方法 |
CN103499842A (zh) * | 2013-09-27 | 2014-01-08 | 长江大学 | 一种微米电磁勘探方法 |
US20160313468A1 (en) * | 2015-04-27 | 2016-10-27 | Pgs Geophysical As | Seismic guided inversion of electromagnetic survey data |
CN106610509A (zh) * | 2015-10-26 | 2017-05-03 | 中国石油天然气股份有限公司 | 一种瞬变电磁资料的时域处理方法 |
CN106547031A (zh) * | 2016-10-09 | 2017-03-29 | 煤炭科学技术研究院有限公司 | 一种基于场变换确定地质电性界面的瞬变电磁探测法 |
CN109752762A (zh) * | 2018-12-11 | 2019-05-14 | 中煤科工集团西安研究院有限公司 | 单发多收观测装置瞬变电场数据动校正方法及装置 |
CN111123264A (zh) * | 2019-12-31 | 2020-05-08 | 桂林理工大学 | 一种考虑介质频散和衰减补偿的探地雷达逆时偏移成像方法 |
Non-Patent Citations (4)
Title |
---|
LI XIU ET AL.: "Inverse synthetic aperture imaging of the ground-airborne transient electromagnetic method with a galvanic source", 《CHINESE JOURNAL OF GEOPHYSICS》, vol. 58, no. 02, pages 158 - 169 * |
吕国印: "瞬变电磁法二维逆时偏移", 《物探与化探》, vol. 22, no. 02, pages 139 - 142 * |
程久龙 等: "矿井瞬变电磁法波场变换与数据处理方法研究", 《煤炭学报》, vol. 38, no. 09, pages 1646 - 1650 * |
陈本池 等: "瞬变电磁场的波场变换研究", 《物探与化探》, vol. 23, no. 03, pages 195 - 201 * |
Also Published As
Publication number | Publication date |
---|---|
CN114200526B (zh) | 2023-08-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP6945895B2 (ja) | 表面波探査方法および端末デバイス | |
CN102027390B (zh) | 对地震数据的联合内插及消除虚反射 | |
CN102305944B (zh) | 一种水下淤积量探测分析处理系统 | |
EP3078991B1 (en) | Method for swell effect and mis-tie correction in high-resolution seismic data using multi-beam echo sounder data | |
CN107817523B (zh) | 绕射波偏移速度的分析方法及装置 | |
Rentsch et al. | Fast location of seismicity: A migration-type approach with application to hydraulic-fracturing data | |
CN108415096B (zh) | 基于牛顿迭代法的水下重力梯度目标探测方法 | |
WO2007081855A2 (en) | Traveltime calculation in three dimensional transversely isotropic (3d tti) media by the fast marching method | |
CN102298156A (zh) | 用于反虚反射地震数据的方法和装置 | |
Mukoyama | Estimation of ground deformation caused by the earthquake (M7. 2) in Japan, 2008, from the geomorphic image analysis of high resolution LiDAR DEMs | |
Zeng et al. | Turning a telecom fiber‐optic cable into an ultradense seismic array for rapid postearthquake response in an urban area | |
CN105093300B (zh) | 一种地质体边界识别方法及装置 | |
Ji et al. | Observation of higher‐mode surface waves from an active source in the Hutubi Basin, Xinjiang, China | |
CN108761540A (zh) | 一种频率域天然电场三维勘探方法 | |
CN104297753B (zh) | 一种基于自适应缩减算子的导航雷达图像反演海面风向方法 | |
CN109212598A (zh) | 基于直达波反演的三维空间二次定位方法 | |
CN111487621B (zh) | 一种基于雷达图像的海表流场反演方法及电子设备 | |
US20130046472A1 (en) | Method of determining the relative position of two detectors at the bottom of the sea | |
CN114200526B (zh) | 电磁勘探方法及装置 | |
CN108508479B (zh) | 一种空地井立体重磁数据协同目标位置反演方法 | |
CN112067113B (zh) | 一种瑞利波在土壤中传播速度的确定方法及装置 | |
CN113075732B (zh) | 一种消除高分辨率小多道地震地层异常起伏的方法 | |
CN111474580B (zh) | 一种基于炮检距矢量片的方位角道集提取方法和系统 | |
US20090052279A1 (en) | Method and system for determining azimuth parameters for seismic data processing | |
CN113608262A (zh) | 利用平动分量计算旋转分量的地震数据处理方法及装置 |
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 |