CN1013001B - 电磁阵列断面勘测方法 - Google Patents

电磁阵列断面勘测方法

Info

Publication number
CN1013001B
CN1013001B CN86103551A CN86103551A CN1013001B CN 1013001 B CN1013001 B CN 1013001B CN 86103551 A CN86103551 A CN 86103551A CN 86103551 A CN86103551 A CN 86103551A CN 1013001 B CN1013001 B CN 1013001B
Authority
CN
China
Prior art keywords
magnetic field
geophysics
electromagnetism
investigation method
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.)
Expired
Application number
CN86103551A
Other languages
English (en)
Other versions
CN86103551A (zh
Inventor
弗朗西斯·X·伯斯蒂克
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.)
University of Texas System
Original Assignee
University of Texas System
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 University of Texas System filed Critical University of Texas System
Publication of CN86103551A publication Critical patent/CN86103551A/zh
Publication of CN1013001B publication Critical patent/CN1013001B/zh
Expired legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V3/00Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
    • G01V3/08Electric 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
    • G01V3/087Electric 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 the earth magnetic field being modified by the objects or geological structures
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V3/00Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
    • G01V3/08Electric 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

Landscapes

  • Engineering & Computer Science (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • Environmental & Geological Engineering (AREA)
  • Geology (AREA)
  • General Physics & Mathematics (AREA)
  • Electromagnetism (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Geophysics (AREA)
  • Geochemistry & Mineralogy (AREA)
  • Geophysics And Detection Of Objects (AREA)
  • Road Signs Or Road Markings (AREA)
  • Measuring Magnetic Variables (AREA)
  • Investigating Or Analysing Materials By Optical Means (AREA)

Abstract

用于地球物理勘探的电磁勘测方法,其中在勘测区中的一点测量沿两个非平行方向的地磁场变化。同时在沿勘测线的若干点测量平行于勘测线的大地电场变化。把这些测得的变化变换到频域,并随后计算出垂直于测得的电场方向的磁场水平分量,计算出勘测线上各测量点的作为频率的函数的阻抗,并利用就多个预定频率的阻抗加权平均值,计算出地下导电率分布。该阻抗加权值是利用对应于对电场进行低通滤波的零相位有限长度权函数计算出来的。

Description

本发明涉及利用测量自然产生或人工感生的地磁场和感生地电(大地)场(总称大地电磁场)的地球物理勘探。
流过地壳的大地电流流量取决于特定点处的地壳构造的导电率或电阻率。如果能测量并绘制这种导电率或电阻率,就能够获得关于该构造的信息,特别是关于碳氢化合物、矿物或地热资源的信息。这在能够采用地球物理勘测的地震方法的地区是特别有用的,例如在沉积岩为厚火山岩层所覆盖的地区。
大地电流的大小、方向和极性是不断变化的,并构成具有各种频率的分量的复杂频谱。
第一种实用的电磁勘测方法是由卡格马亚德(Cagnird)在美国专利第2,677,801号中描述的。它包括测量并记录大地场的一个水平分量在一段时间中的变化,并同时测量并记录地磁场的正交分量的变化。随后用付立叶分析法把这些测量结果变换成频率分量。电场频率分量与磁场频率分量的比值是一种波阻抗,它是频率的函数。因为电磁波透入大地的透入深度与该波的频率和该地的导电电率成反比,故可用该波阻抗来估算沿通过地面的垂直方向的导电率分布。卡格尼亚德用一种数学模型来做这种估算,在他的这种模型中,导电率仅随深度变化,即所谓一维(1-D)模型。
这种方法随后受到了另一些人的发展,他们采用了导电率随水平坐标和深度变化的模型,即二维(2-D)模型。在此模型中,导电率沿其不变的坐标被称作“走向”(strike)。在这种模型中, 根据大地电磁场的极性,有两种要考虑的情况。在这两种情况中,电场分别平行和垂直于走向发生极化。理论研究表明,可结合为一维而发展的技术,用E平行阻抗函数来相当准确地估算该函数计算地点正下方的导电率。然而,E垂直阻抗函数有非常困难的特性,并且不能直接对其进行变换,以给出足够准确的导电率估算。解决该问题普通的方法是利用在一个以上地点得到的曲线,并且迭代法修正该模型中的导电率分布,直到获得与观察结果符合得最好的理论E垂直曲线为止。该方法有许多缺点,特别是它需要大量的计算,并且无法保证所得解的正确性。
当勘测区导电率确实仅沿一或两个方向变化时,这些方法是有效的。然而不幸的是,这种两维的变化是少见的。当对具有三维导电率分布的构造进行大地电磁测量时,就会碰到下列问题。首先,尽管存在若干确实一套坐标轴的特定方法,但要象在二维情况下那样去确定一套原则坐标轴一般是不可能的。第二,不论用什么坐标系,要想把电场分成两种具有不同性质的不同情况(象二维情况下那种E平行和E垂直)是不可能的。所希望的E平行情况的性质首先消失了。而两种情况都有与E垂直相类似的性质。虽然为解决这个问题而发展了许多特定的方法,但它们的成功是有效的。第三个问题是显著增加的复杂性使确保对所得数据的正确理解的难度大大增大了。
减少第三个问题的一种方法是进行更多的测量。在经典的大地电磁方法中,磁场测量通常比电场测量更困难,使得大量测量既费时又费钱。然而,一个观测到的事实是磁水平分量随距离的变化通常比电场水平分量的要慢得多。这导致了大地电磁一大地方法,它包括在勘测区域上分布的有限个地点测量磁场的两个正交分量和在更多的地点测量电场的两个正交分量。随后,用与大地电磁模型中所描述的相类似 的技术。对这些数据进行处理。然而,当导电率是三维分布时,其结果仍是不十分可靠的。在美国专利第4,286,218中,描述了这种技术的一种类型,它采用了与地震勘探的“滚动”(roll-along)方法相似的方法。沿勘测线的电场重迭测量改善了信噪比,而在该线的两端测量磁场,而用插入法导出电场测量点的中间磁场。
所有这些传统方法的主要缺点在于当地下构造的电和磁特性沿所有三个方向变化时,它们给出的结果非常不可靠,以前克服这种不可靠性的努力包括大量计算,而这就所得结果来说是难以证明的。
传统大地电磁勘测方法的另一缺点是在每个勘测点都要测量两个最好是正交方向的电场。这使该方法很难用于海上勘测。在美国专利第4,210,869号中描述了一种海上大地电磁勘测方法,但此方法依赖于测量海底上的电极的位置和取向的声学方法。这些测量的准确性是有限的,从而使得对勘测结果的分析变得不可靠了。
传统大地电磁方法的再一缺点是不论勘测区域多小,都要在勘测区域内的数个点进行测量。另外,为区分视在E平行和E垂直分量,必须在各个E测量点确定垂直分量。如上所述,如未进行这种区别,就很可能会对所得数据作出错误的分析。
先有技术的上述和其他缺点可通过提供一种电磁地球物理勘测方法而得到克服,在勘测点的大地构造的导电率具有三维变化时,该方法能用较简单地计算方法给出可靠的结果。
根据本发明的一个方面,该电磁地球物理勘测方法包括:在勘测区域中的至少一个地点测量地磁场沿至少两个非平行方向的变化,并同时测量该磁场的变化;在沿勘测线的多个勘测点测量平行于该勘测线的地电场变化,随后把测出的数据转换成频率分量。随后的步骤 是:从两个非平行方向的地磁场测量结果,计算出作为频率的函数的、与各地点的测出电场方向相正交的地磁场水平分量;计算出各地点作为频率的函数的阻抗,该阻抗为在各点测得的电场与垂直于该测得电场方向的地磁场水平分量之间的比值。随后,对多个预定的频率,以适当的方法计算出该阻抗的加权平均值,使计入各加权平均值的阻抗数目随频率的减小而增加,以使该数目与该频率的电磁波透入大地的有效透入深度大致成正比。最后,由这些阻抗加权平均值,计算出作为频率的函数的该勘测线下的大地导电率的分布。
在一个实施方案中,测量大地磁场的两个方向是相互垂直的。
在另一个实施方案中,测量平行于勘测线的地电场变化的步骤包括测量与大地电接触并沿勘测线以大致相等的间距设置的一或多对电极之间的电位差,并用该电极间的距离除所测得的该电位差。
在又一种实施方案中,所有电极间的电位差是同时测量的。
在一种供选用的实施方案中,同时测量一组相邻电极的成员之间的电位差(该组电极由少于全部电极的电极组成)并顺序地对各组进行测量。
在另一种可供选用的实施方案中,一组中的至少一个电极同时也是相邻组中的成员。
在最佳实施方案中,阻抗的加权平均值是经过下列步骤计算出来的;利用零相位有限长权函数,对每一频率,通过选择一适当宽度确定该权函数的宽度;利用选出的函数获得加权平均阻抗;利用如此得到的加权平均值计算近似透入深度;把计算出的透入深度同期望透入深度相比较,利用计算和期望深度之差选择更合适的宽度,以迭代方法重复该过程,直至透入深度的计算达到预定的精确度。
根据本发明的另一个方面,该电磁地球物理勘测方法包括:在勘 测区域中的至少一个点测量在至少两个非平行水平方向上的地磁场变化;在测量该地磁场变化的同时,测量与大地电接触的相邻电极之间的电位差,这些电极沿勘测线以大致相等的间隔设置,且同时进行该电位差测量,并将测得的电位差变化变换成频率的函数。随后,对于多个预定频率,以适当方法计算测得电位差的加权平均值,以使计入各加权平均值的电位差的数目随频率的降低而增大,从而使该数目大致正比于该频率的电磁波进入大地的有效透入深度。由沿两个非平行方向的该磁场测量结果,计算出垂直于勘测线方向的地磁场水平分量,从电位差的加权平均值同该磁场水平垂直分量的比值计算出沿勘测线的加权阻抗,并从该加权阻抗计算出该勘测线下方作为深度的函数的大地导电率分布。
根据本发明的再一个方面,该电磁地球物理勘测方法还包括在该勘测线的多个点测量地磁场的垂直分量,并利用这些测得垂直磁场分量连同计算出的勘测线下导电率分布来确定勘测线下的地球物理构造的步骤。
本发明的一个优点是提供了当勘测线下方的大地地下构造的导电率分布在三个方向上都发生变化时估算该导电率分布的一种新颖的改进方法。本发明的另一优点是在该勘测方法中电和磁场测量装置的设置很简便。再一个优点是本发明估算复杂导电率分布所需的计量比先有方法简便。
通过下面结合附图对最佳实施方案所作的详细说明,本发明的上述和其它目的及优点将变得更为明显。
在附图中:
图1是表示勘测场的示意图。
在图1中显示了根据本发明的勘测场布置。在勘测区中选定用于 磁场测量的适当地点,它最好远离局部的非自然电磁干扰源,如发电机、高压输电线,无线电发射机及繁忙的道路。设置传统形式的磁场检测器10和12,以测量两个非平行方向上的水平地磁场分量的变化。在最佳实施方案中,这些检测器彼此成直角,其中一个检测器平行于选定的勘测线,可用无线电遥测装置(如图1所示)或用直接电缆传输或通过电话线把来自检测器10和12的读数实时地传送到处理站14。
在一种供选择的替代方案中,该数据被记录在诸如位于记录场地的磁盘、电子计算机存贮器或磁带的存贮介质中。可把该存贮数据延时传送到处理站,也可将包含该数据的存贮介质本身运送到处理站。这样能建立独立的、自动运转的记录站。
虽然单一磁测量站在提供一致的基准信号方面是较为可取的,但如果情况需要,也可采用一个以上的磁场站配置。如果随着勘测的进展,需要移动磁场站,可以在新的位置设置第二个磁场站,它具有与第一磁场站类似地对准的两个检测器,并且进行一适当时期的同时测量,以将用第二站得到的数据平移到与第一站所获得的数据相同的基准系中,这是内行人所熟知的步骤。
随后确定勘测线16。该线最好是直的,但该方法也可用于弯曲的勘测线。沿该勘测线设置与大地接触的电极18,它们最好彼此相距同样的距离。用导线把相邻的电极连到电位差检测器20。用无线电、电缆线电话把电位差检测器20的续数传送到处理站14,也可在记录场地将其记录在大容量存贮介质上,以便随后送到处理站,就象磁场读数那样。
E(电)场分量是通过用电极间距L除测得的两电极间电位差而得到的。这种E场测量结果实际上是沿连接两个电极的导线的E场分 量的加权平均值。因此,该设置起波数空间中的空间低通滤波器的作用,它具有矩形或矩形波串函数的因有加权因子。检测器在波数空间中的响应正比于正弦(Six/x)函数,并有稍小于L的倒数的截止波数。因此,内行人可通过选择适当的L值来提供所需的勘测分辨率。
沿勘测线16进行一系列的E场分量测量。如图1所示,这些检测器是头尾相接地设置的。这种设置沿勘测路径对连续的E场分量进行空间低通滤波,并随后以固定的空间间隔L对这种经过滤波的分量进行取样。这种取样是与众所周知的取样定理为防止取样数据的混淆所作的要求相符合的。在要取得导电率数据的点两侧的测量跨距应大致等于所要勘测的最大深度。E场测量可在整个勘测线上同时进行,也可一次只对一组相邻电极进行,然后再对另一组进行,直到顺序地覆盖了整条勘测线。当勘测较深的深度因而所需的测量跨距很长,因而难以同时在整条勘测线上同时进行测量时,后一种方法是较好的。当进行顺序测量时,可在相邻组的成员间造成小的重迭,以提供有助于减小噪声的冗余,但这样做会降低勘测速度,所以最好只在噪声成为严重问题的情况下采用此方法。
当勘测线是曲线时,该方法必须进行修正。首先,对阻抗施行具有加权和形式的低通滤波,该阻抗描述了与一组相邻路径段中的各段相切的电场分量同与为该组路径段确定的有效方向垂直的磁场分量之间的关系。简单的余弦钟形权函数也被用于弯曲勘测线。
位于路径组中心点a的勘测线段的有效方向是通过确定该组的该有效方向同任意基准方向(如北)之间的夹角φ0(a)而确定的。该角的确定关系为:
Figure 86103551_IMG2
其中,Wi(a)是余弦钟形权,N1(a)和N2(a)是要包括余弦钟形函数长度区域内所有非零权的和的下和上界,而φi是勘测线段同基准方向间的夹角。
随后用由余弦钟形函数加权的线段阻抗和来计算近似透入深度。并且勘测点a附近的勘测线段,由下列确定关系计算有效长度LE(a):
Figure 86103551_IMG3
其中,L是各路径段的长度(假设所有路径段的长度都是相等的),而Wmax是包括在该和中的Ni权中的最大者。对弯曲勘测路径的情况,要求有效长度LE(a),而不是余弦钟形函数的长度,为透入深度的一个给定倍数。
所要求的LE(a)值是用迭代法得到的。该方法开始于为余弦钟形函数假设适当的长度。随后计算出权Wi(a)并利用它得到角φ0(a)和有效长度LE(a)。从线段阻抗的加权和计算出近似透入深度。比较有效长度和透入深度,并利用计算出的和给定倍数之差来为余弦钟形函数确定新的长度。以迭代方式重复该过程,直至达到适当的收敛程度。对勘测线16的每一待测位置和频率,都采用该方法。
在最佳实施方案中,随后计算出阻抗,它是E场测量结果的频率分量与同时测得的沿E场正交方向的磁场分量的比值。这种表为Zxy的阻抗是水平电场同磁场空间分量间的二阶阻抗张量关系的元素。Zxy由下列方程确定:
Ex=ZxxHx+ZxyHy
其中Hx是平行于E场测量结果的磁场测量结果的变换,而Hy是与E场垂直的磁场测量结果的变换。可由记录数据和大地电磁先有技术领域的内行人所熟知任一种技术给出的限定方程来估算阻抗Zxy
在本发明的一个供选择的实施方案中,也可沿与勘测线成直角的方向进行电场测量。在此情况下,可由测得数据的频率分量估算表为Zxy并由方程
Ey=ZyxHx+ZyyHy
限定的另一阻抗张量元素。
作为本发明的一部分,发明人发现除在大波数值外,函数Zxy 在三维解析中的表现同E平行阻抗函数在传统两维大地电磁分析中的表现相似。在本发明中,通过用相对波数空间而言的低通滤波器来衰减在大波数值处由Zxy反映的电场相关性,就可在三维大地导电率变化的情况下获得E平行函数的所需特性。
另外,还发现低通滤波操作的关键因素在于截止波数的值。已经确定出,对于增大的波数值,在等于电磁波透入大地的透入深度的波长值处,Zxy同二维E平行函数开始有显著的不同,因此该波长值的倒数就是一个合适的截止值,尽管在本发明的范围内也可采用其他的值。
先有技术中众所周知的是,在波数和空间域中,与低通滤波相对应的线性运算是与核权函数的卷积。该核函数的空间长度为截止波数的倒数,因此该核的长度随波数的减小而增大。为在所有波数都能适当地进行加权运算,在要获得导电率探测的各点就必须有足够的E场测量跨距,以便能适应Exy阻抗的空间低通滤波所要求的最长核函数。在最佳实施方案中,可用简单的余弦钟形函数作加权核函数。
在最佳实施方案中在沿勘测线16的任一位置用于任一频率的余弦钟形函数的准确长度是用迭代法确定的。在每一要勘测的任置和频率,对余弦钟形函数的长度作出适当的初始假定。随后用该函数对阻抗进行加权,并用该阻抗的和计算近似透入深度。该余弦函数的长度被要求为这个透入深度的指定倍数。该最佳实施方案采用等于一的乘数。对余弦钟形函数的长度和透入深度进行比较,并用计算出的和指定的乘数间的差为余弦钟形函数确定新的长度。重复这一过程,直至计算出的倍数相对指定的倍乘数达到适当的收敛程度。在实践中发现,该方法是行之有效的,并且收敛很快。
当在整个勘测路径上同时进行所有的E场测量时,可采用一种供 选择的方法来得到加权阻抗。该加权方法仍采用具有适当长度的余弦钟形函数,并对转换到频域的E场测量结果进行,并随后把这些加权测量结果用于阻抗张量中,以获得加权Zxy值。
一旦该阻抗在所有感兴趣的频率都得到了空间加权,就对该函数进行反演,以获得导电率沿勘测线随深度的变化。该实施方案中所用的这种反演方法就是那些在先有技术的一维分析中所用的方法。如果采用更复杂的方法所获得的精确度增加足以抵挡它们所要求的计算量增加时,也可采用这种复杂方法。
进行各种计算的方法并非本发明的组成部分,但对内行人明显的是借助程序数字电子计算机可以方便地进行这种计算。进行付立叶变换,余弦钟形加权函数和一维反演过程的程序都是容易得到的。
在根据本发明方法进行的勘测中,在任何点的估算导电率是勘测线正下方及其横向的地下导电率加权平均值。当大体积的导电部分位于勘测线侧方但在该截面中仍显而易见时,就称该截面受到了“侧击”(Side-swiped)。根据已公开的数学研究结果,可用勘测点磁场的垂直分量来判断侧击。因此,在本发明一个供选择的实施方案中,在沿勘测路径的点处测量该磁场的垂直分量,并根据已公开的分析方法,用这些测量结果确定用上述方法确定的截面是否含有侧击。
当碰到或预见到磁场的非常不均匀性时,或在有不寻常噪声的环境下,或是要求高精确度的情况下,可进行沿两个非平行水平方向的附加磁场测量。这些测量可沿勘测线进行,但最好是在遍布勘测区的各点进行。这两个方向最好是正交的,其中一个方向最好平行于该勘测线。
如果在N个独立点同时或不同时测量磁场的x和y分量,则假设 对电场阵列的每一间距要估算的阻抗Zax和Zay(该间距就是两相邻电极间的距离)满足下列关系:
E a =Z ax ( Σ i=1 N U i H xi )+Zay( Σ i=1 N V i H yi )
其中Ui和Vi是适当的加权因子,它们所用内行人所熟知的各种方法确定。
把基准地点的x和y场分量表为Rx和Ry,这些场是电场、磁场或它的组合,则其与第i地点的x和y磁场分量Hxi和Hyi的关系由下列式给出:
Hxi=TxxiRx+TxyiRy
Hyi=TyxiRx+TyyiRy
转换系数Txxi、Txyi、Tyxi和Tyyi可用大地电磁勘测中的任何已知传统基准技术来确定。
将这些方程式代入关于Ea的表达式,得出:
E a =Z ax Σ i=1 N (V i T xxi +V i T yxi )R x
+Zay(ViTxyi+ViTyyi)Ry
从这个式子,可以用传统装置测得的Ea、Rx和Ry值得到Zax和Zay。随后用同上述与核权函数的卷积相等价的方法,对如此为电场阵列的各个间距确定的阻抗进行滤波。虽然附加磁测量会使测量费用更高,但以此方式计算出的导电率分布的精确度应高于那些仅包含单个地点磁测量的导电率分布。
本发明的上述实施方案涉及到利用自然或人工感生的交变电流的勘测方法。本发明还可用于采用直流、受控激励源的电勘测。在这种勘测中,信号进入大地的透入深度同源与电场测量点之间的距离有关。用直流受控源在大地中激发出电场(与此有关的技术对内行人是熟知的)。沿与该源相距一段距离的勘测线,利用上述技术,即测量阵列中各对相邻电极间的电位差并将该测得电位差除以电极间距离,测量该激发电场。随后,用传统处理技术,把测得电场数据转换成视在电阻,它是测得电场同受控源输出的电流的比值。对多个相距源的预定距离用适当方法获得该视在电阻的加权平均值,该方法等价于沿勘测线对电场采用利用核权函数的低通滤波技术,如上所述。在此情况下,该核函数的空间长度同截止波数成倒数,该截止波数同源和电场测量阵列间的距离成反比。
本发明的一个优点在于它可用结构比先有大地电磁勘测方法中的更简单的电场检测器。这在海上勘测及非常复杂地形的勘测中是特别有益的。在同样环境下可获得好处的另一优点是能以更简单的方式测量磁场,只需要一个磁场检测器地点,在进行靠近海岸线的海上勘测时,可将该地点选在了海岸上;或将该点取在海上平台附近,另外,与传统大地电磁方法不同,该勘测不需要进行磁场垂直分量的测量,虽然为了增加与勘测线侧方的构造有关的信息,最好进行这种测量。
本发明的另一优点是它易于适用在较宽范围的勘测应用中。在不 以任何方式限制此处公开并提出权利要求的本发明范围的情况下,可将其用于矿物和碳氢化合物勘探、地下和公共管道的位置探测、放射性或毒性物质从存贮地的泄漏的检测,地下水源探测、考古研究、以及体内组织、骨骼和器官导电率的检测。
对于内行人来说,其他优点和修正都是显而易见的。因而就其广义方面而言本发明并不限于所显示和描述的说明例子、代表装置或特定细节。相应地,在不背离所公布的一般发明概念的情况下,是可以在细节上有所不同的。

Claims (13)

1、一种电磁地球物理勘测方法,包括:
在勘测区内的至少一个点测量沿至少两个水平非平行方向的大地磁场变化;
在测量该磁场变化的同时,在多个适当分布的点测量平行勘测区中的勘测线的大地电场变化,以对沿该勘测线的电场空间变化进行充分取样,作法是同时测量与大地电接触的一对或多对相邻电极间的电位差变化并用测得的各对相邻电极间电位差除以该对电极间的距离,所述电极对沿勘测线基本上等间距设置。
2、如权利要求1的电磁地球物理勘测方法,其中测量大地磁场的两个方向是彼此垂直的。
3、如权利要求1的电磁地球物理勘测方法,其中测量大地磁场变化的一个方向大致与勘测线平行。
4、权利要求1的电磁地球物理勘测方法,其中勘测线大致是直的。
5、如权利要求1的电磁地球物理勘测方法,其中同时测量由少于所有电极的电极构成的一组相邻电极中的各对相邻电极间的电位差,并顺序地对各组进行测量。
6、如权利要求5的电磁地球物理勘测方法,其中一个组中的至少一个电极对也是相邻组的成员。
7、一种电磁地球物理勘测方法,包括:
在勘测区中的至少一点测量沿至少两个水平非平行方向的大地磁场变化;
在测量磁场变化的同时,在基本上等间距分布的多个点同时测量平行和垂直于勘测区中的勘测线的大地电场变化,以对沿该勘测线的被测电场分量的空间变化进行充分取样;
其中,测量平行于勘测线的大地电场变化的步骤包括测量与大地电接触的一对或多对相邻电极间的电位差变化,并用各相邻电极间的所测得的电位差除以该对电极间的距离;
其中测量垂直于勘测线的大地电场变化的步骤包括测量与大地电接触的位于与垂直于勘测线的一条线上的一对或多对相邻电极间的电位差变化并用各对相邻电极间的所测得电位差除以该对电极间的距离。
8、如权利要求7的电磁地球物理勘测方法,其中测量大地磁场的两个方向是彼此垂直的。
9、权利要求7的电磁地球物理勘测方法,其中测量大地磁场变化的一个方向大致同勘测线平行。
10、如权利要求7的电磁地球物理勘测方法,其中在由少于沿勘测线的所有点的构成的一组相邻点同时测量大地电场变化,并顺序地对各组进行测量。
11、如权利要求1的电磁地球物理勘测方法,其中一组中的至少一点也是相邻组中的成员。
12、如权利要求1或7的电磁地球物理勘测方法,进一步包括在沿勘测线的各个点同时测量大地磁场在垂直方向的变化的步骤。
13、一种电磁地球物理勘测方法包括在勘测区中的第一基准点测量沿至少两个水平非平行方向的大地磁场变化并同时在至少沿一个平行于该勘测区中的勘测线的多个点测量大地电场变化,所述点相互隔开,以沿该勘测线对被测电场分量的空间变化进行充分取样,在这种勘测方法中,为测量大地磁场变化确定第二基准点的方法,包括:
在该第二基准点测量沿至少两个水平方向的大地磁场变化,该水平方向大致平行于在第一基准点进行大地磁场测量的方向;
在足以使由第二基准点获得的数据同由第一基准点获得的数据相关的时间间隔内,同时继续在第一基准点测量大地磁场的变化。
CN86103551A 1986-04-25 1986-05-26 电磁阵列断面勘测方法 Expired CN1013001B (zh)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
US854,147 1986-04-25
US06/854,147 US4757262A (en) 1984-07-31 1986-04-25 Method for geophysical exploration using electromagnetic array

Publications (2)

Publication Number Publication Date
CN86103551A CN86103551A (zh) 1987-11-04
CN1013001B true CN1013001B (zh) 1991-06-26

Family

ID=25317853

Family Applications (1)

Application Number Title Priority Date Filing Date
CN86103551A Expired CN1013001B (zh) 1986-04-25 1986-05-26 电磁阵列断面勘测方法

Country Status (7)

Country Link
US (1) US4757262A (zh)
EP (1) EP0242424B1 (zh)
JP (1) JPS62257083A (zh)
CN (1) CN1013001B (zh)
AT (1) ATE75055T1 (zh)
AU (3) AU593800B2 (zh)
DE (1) DE3684911D1 (zh)

Families Citing this family (29)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4835473A (en) * 1984-07-31 1989-05-30 Board Of Regents, The University Of Texas System Method using a d.c. source for determining resistivity distribution of the earth as a function of depth
JPS63269083A (ja) * 1987-04-27 1988-11-07 Maakurando Kk 電磁式地下探査方法
JPS63269084A (ja) * 1987-04-27 1988-11-07 Maakurando Kk 電磁式地下探査装置
JPS63293486A (ja) * 1987-05-27 1988-11-30 Maakurando Kk 電磁探査方法
JPS6449995A (en) * 1987-08-20 1989-02-27 Maakurando Denshi Kogyo Kk Electromagnetic surveying method
JPH0641983B2 (ja) * 1988-04-18 1994-06-01 動力炉・核燃料開発事業団 商用周波数信号を利用した地下探査法及び装置
JPH0641985B2 (ja) * 1988-04-18 1994-06-01 動力炉・核燃料開発事業団 地下比抵抗探査におけるデータ処理方法および装置
US4862089A (en) * 1988-09-23 1989-08-29 Amoco Corporation Method of magnetotelluric exploration using a zigzag array
US5043667A (en) * 1989-04-21 1991-08-27 Amoco Corporation Method of magnetotelluric exploration using areal arrays
US5745384A (en) * 1995-07-27 1998-04-28 Lucent Technologies, Inc. System and method for detecting a signal in a noisy environment
US5770945A (en) * 1996-06-26 1998-06-23 The Regents Of The University Of California Seafloor magnetotelluric system and method for oil exploration
CA2441226C (en) * 2000-07-26 2013-03-12 Northern Digital, Inc. Method for determining the position of a sensor element
US6950747B2 (en) * 2003-01-30 2005-09-27 Kent Byerly Methods of processing magnetotelluric signals
CN100483154C (zh) * 2006-08-10 2009-04-29 中国地质大学(武汉) 感应式地磁场测量实验仪
US8055446B2 (en) * 2006-10-11 2011-11-08 Byerly Kent A Methods of processing magnetotelluric signals
CN101556340B (zh) * 2008-04-10 2011-08-03 中国石油集团东方地球物理勘探有限责任公司 三维小面元大地电磁连续阵列数据采集方法
US9335432B2 (en) 2010-08-30 2016-05-10 King Abdulaziz City For Science And Technology Semi-permeable terrain geophysical data acquisition
GB201016556D0 (en) * 2010-10-01 2010-11-17 Senergy Technology Ltd Methods for providing for correcting data and associated apparatus
CN103389514B (zh) * 2013-07-18 2016-04-20 中国科学院地质与地球物理研究所 一种大地电磁去噪装置及大地电磁去噪方法
CA2919205C (en) * 2013-07-24 2021-05-11 Schlumberger Canada Limited Method to predict local geomagnetic disturbance field and its practical application
US9651700B2 (en) * 2013-08-29 2017-05-16 Saudi Arabian Oil Company Mapping resistivity distribution within the earth
CN106321085A (zh) * 2015-07-02 2017-01-11 中石化石油工程技术服务有限公司 激发极化电位测井仪归一化方法
CN107229077B (zh) * 2017-05-26 2019-03-26 湖南普奇地质勘探设备研究院(普通合伙) 物探仪及其剖面图成像方法
CN108508281B (zh) * 2018-03-19 2020-04-17 江苏伏波海洋探测科技有限公司 舰船静电场的基于点电源频域法的深度换算方法
CN109342978B (zh) * 2018-11-06 2021-03-30 中国石油天然气集团有限公司 大地电磁各向异性采集系统、方法及设备
CN109884714B (zh) * 2019-03-05 2021-11-16 中国地质科学院地球物理地球化学勘查研究所 一种可控源电磁测量方法、装置及其存储介质
CN112162323B (zh) * 2020-09-28 2024-03-08 长春工程学院 基于人工场源频率域电法的地下不良地质勘测方法及系统
CN113281401B (zh) * 2021-04-29 2022-08-19 同济大学 无砟轨道隐蔽病害的检测方法、系统与装置
KR102501035B1 (ko) * 2022-05-17 2023-02-21 (주)미도지리정보 도로 매립형 이중구조 수치지도용 수치측량 기준점 장치

Family Cites Families (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US2677801A (en) * 1950-10-06 1954-05-04 Centre Nat Rech Scient Geophysical exploration method
US3188558A (en) * 1961-12-06 1965-06-08 California Research Corp Telluric sounding method in which recorded variations in the telluric field with time are converted into amplitude-versus-frequency spectra
JPS52801B2 (zh) * 1974-01-25 1977-01-11
FR2272442B1 (zh) * 1974-05-21 1977-03-11 Aquitaine Petrole
FR2307279A1 (fr) * 1975-04-07 1976-11-05 Procedes Electro Magnetiques Procede et dispositif pour l'etude du sous-sol et de sa nature electro-magnetique
US4041372A (en) * 1975-09-08 1977-08-09 Continental Oil Company Apparatus for multi-channel induced polarization surveying
US4210869A (en) * 1978-02-06 1980-07-01 Mobil Oil Corporation Rotating telluric field measurements
US4286218A (en) * 1979-08-06 1981-08-25 Mobil Oil Corporation Multiple site magnetotelluric measurements
US4349781A (en) * 1980-01-07 1982-09-14 The Regents Of The University Of California Superconducting gradiometer-magnetometer array for magnetotelluric logging
US4473800A (en) * 1980-11-14 1984-09-25 Mobil Oil Corporation Method for detecting magnetotelluric subsurface anomalies
US4449099A (en) * 1981-08-18 1984-05-15 Mobil Oil Corporation Recording decimated magnetotelluric measurements with coherence
US4467283A (en) * 1981-09-23 1984-08-21 Geo-Optics, Ltd. Electrical geophysical exploration system with addressable current probes
JPS60133387A (ja) * 1983-12-19 1985-07-16 Chinetsu Gijutsu Kaihatsu Kk 地下探査方法及びその探査装置
US4591791A (en) * 1984-07-31 1986-05-27 Board Of Regents, University Of Texas System Electromagnetic array profiling survey method

Also Published As

Publication number Publication date
AU4547289A (en) 1990-03-08
EP0242424A3 (en) 1988-12-21
EP0242424B1 (en) 1992-04-15
AU593800B2 (en) 1990-02-22
DE3684911D1 (de) 1992-05-21
AU5733486A (en) 1987-11-19
US4757262A (en) 1988-07-12
CN86103551A (zh) 1987-11-04
AU4547189A (en) 1990-03-08
ATE75055T1 (de) 1992-05-15
AU607646B2 (en) 1991-03-07
EP0242424A2 (en) 1987-10-28
JPS62257083A (ja) 1987-11-09
AU607647B2 (en) 1991-03-07

Similar Documents

Publication Publication Date Title
CN1013001B (zh) 电磁阵列断面勘测方法
Okpoli Sensitivity and resolution capacity of electrode configurations
Legchenko et al. Inversion of surface NMR data
Tsourlos Modelling, interpretation and inversion of multielectrode resistivity survey data.
Papadopoulos et al. Two‐dimensional and three‐dimensional resistivity imaging in archaeological site investigation
CN1168997C (zh) 远程储层电阻率映像
CA2358124C (en) Water saturation and sand fraction determination from borehole resistivity imaging tool, transverse induction logging and a tensorial water saturation model
Akingboye et al. Insight into seismic refraction and electrical resistivity tomography techniques in subsurface investigations
US3256480A (en) Method of detecting geologically anomalous bodies lateral to a well bore by comparing electrical resistivity measurements made using short-spaced and long-spaced electrode systems
CN102707323B (zh) 用于地质勘探的可控源音频磁场测深法
US9891339B2 (en) Method and apparatus for detecting and mapping subsurface anomalies
US4591791A (en) Electromagnetic array profiling survey method
Brunner et al. Investigation of a Tertiary maar structure using three-dimensional resistivity imaging
US11892587B2 (en) Distributed three-dimensional (3D) induced polarization (IP) data acquisition and processing device and method
Van Zijl A deep Schlumberger sounding to investigate the electrical structure of the crust and upper mantle in South Africa
CN1209499A (zh) 阵列电极系成象测井方法
Gupta et al. Straightforward inversion of vertical electrical sounding data
Hautot et al. Groundwater electromagnetic imaging in complex geological and topographical regions: A case study of a tectonic boundary in the French Alps
Domenzain et al. Efficient inversion of 2.5 D electrical resistivity data using the discrete adjoint method
Brass et al. Resistivity profiling with different electrode arrays over a graphite deposit
Drahor et al. An application of the self‐potential (SP) method in archaeogeophysical prospection
CN111239835A (zh) 一种适用于狭长地形的正交传导电测深方法
Nemesi et al. Telluric map of Transdanubia
CN112860820B (zh) 一种通用地球物理数据存储方法及装置
Akingboye et al. PREGLED TEHNIKA SEIZMIČKE REFRAKCIJE I TOMOGRAFIJE ELEKTRIČNOM OTPORNOŠĆU U ISTRAŽIVANJU PODZEMLJA

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C13 Decision
GR02 Examined patent application
C14 Grant of patent or utility model
GR01 Patent grant
C19 Lapse of patent right due to non-payment of the annual fee
CF01 Termination of patent right due to non-payment of annual fee