CN113917541B - 基于直流电法勘探获取地下介质电各向异性的方法和装置 - Google Patents

基于直流电法勘探获取地下介质电各向异性的方法和装置 Download PDF

Info

Publication number
CN113917541B
CN113917541B CN202111157020.6A CN202111157020A CN113917541B CN 113917541 B CN113917541 B CN 113917541B CN 202111157020 A CN202111157020 A CN 202111157020A CN 113917541 B CN113917541 B CN 113917541B
Authority
CN
China
Prior art keywords
apparent resistivity
background
resistivity
potential difference
normalized
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
CN202111157020.6A
Other languages
English (en)
Other versions
CN113917541A (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.)
Hunan University of Science and Technology
Original Assignee
Hunan University of Science and Technology
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 Hunan University of Science and Technology filed Critical Hunan University of Science and Technology
Priority to CN202111157020.6A priority Critical patent/CN113917541B/zh
Publication of CN113917541A publication Critical patent/CN113917541A/zh
Application granted granted Critical
Publication of CN113917541B publication Critical patent/CN113917541B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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/02—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation operating with propagation of electric current
    • G01V3/04—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation operating with propagation of electric current using dc
    • 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/02—Investigating or analysing materials by the use of electric, electrochemical, or magnetic means by investigating impedance
    • G01N27/04—Investigating or analysing materials by the use of electric, electrochemical, or magnetic means by investigating impedance by investigating resistance
    • G01N27/041—Investigating or analysing materials by the use of electric, electrochemical, or magnetic means by investigating impedance by investigating resistance of a solid body
    • 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

Landscapes

  • Life Sciences & Earth Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Chemical & Material Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Chemical Kinetics & Catalysis (AREA)
  • Electrochemistry (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Geology (AREA)
  • Geophysics (AREA)
  • Environmental & Geological Engineering (AREA)
  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Health & Medical Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Geophysics And Detection Of Objects (AREA)
  • Investigating Or Analyzing Materials By The Use Of Electric Means (AREA)

Abstract

本申请涉及一种基于直流电法勘探获取地下介质电各向异性的方法和装置。所述方法包括：通过在测区任意位置布设发射场源后，采用相互正交的电极在测区内进行全区观测，得到不同测点之间相互正交的两个方向的电位差，根据电位差计算视电阻率，计算归一化视电阻率，根据选取的背景电阻率模型，计算背景视电阻率，计算归一化背景视电阻率，根据实际归一化电阻率和归一化背景视电阻率，计算电阻率各向异性系数。本发明场源布设灵活简单。采用电阻率各向异性系数可有效压制视电阻率零线的影响，很好地反应出地下各向异性结构体在地面投影以及形状，还可避免求解张量电阻率，效率高。采用归一化视电阻率可有效抵消场源附近视电阻率的误差，提升精度。

Description

基于直流电法勘探获取地下介质电各向异性的方法和装置
技术领域
本申请涉及地球物理勘探领域，特别是涉及一种基于直流电法勘探获取地下介质电各向异性的方法、装置、计算机设备和存储介质。
背景技术
直流电法作为一种相对成熟的地球物理方法，广泛应用于矿产勘查，水文环境检测，工程勘探等领域。地球表层岩石由于组分不均匀、地质作用等形成复杂的构造和断层，具有各项异性的特征，可通过利用直流电法勘探在地表或井中观测获取地下介质的电各向异性特征。目前直流电法的主要勘探方法为Bibby(1986)提出的P2不变量视电阻率张量测量。该测量方法为采用相互正交的偶极源AB和CB在测区外发射，并在测区内采用相互正交的电极测量电场另外需计算对应的各向同性背景模型的电流密度/>最后利用表达式/>计算出张量电阻率/>并用来表征地下各向异性介质的视电阻率。很显然，该方法需要在测区外要求布设两个偶极场源，另外还必须求解张量电阻率，最终通过计算才得出P2来表示地下介质电阻率各向异性特征，这样一来野外场源布设复杂，处理过程复杂，工作效率较低，特别是数据量大时，处理起来相当麻烦。
发明内容
基于此，有必要针对上述技术问题，提供一种能够实现全区观测，场源布设简单、数据处理简单可靠的基于直流电法勘探获取地下介质电各向异性的方法、装置、计算机设备和存储介质。
一种基于直流电法勘探获取地下介质电各向异性的方法，所述方法包括：
在测区任意位置布设发射场源后，采用相互正交的电极在测区内进行全区观测，得到不同测点之间的第一方向电位差和与第二方向电位差；所述发射场源为点源或偶极源；所述第一方向与所述第二方向相互正交；
根据所述第一方向电位差和所述第二方向电位差，得到第一方向视电阻率和第二方向视电阻率，根据所述第一方向视电阻率和所述第二方向视电阻率，得到归一化视电阻率；
根据预设的地下介质背景电阻率模型，计算所述不同测点之间的背景第一方向电位差和背景第二方向电位差；
根据所述背景第一方向电位差和所述背景第二方向电位差，得到背景第一方向视电阻率和背景第二方向视电阻率，根据所述背景第一方向视电阻率和所述背景第二方向视电阻率，得到归一化背景视电阻率；
根据所述归一化视电阻率和所述归一化背景视电阻率，得到测区地下介质的电阻率各向异性系数。
在其中一个实施例中，还包括：所述发射场源为点源或偶极源；当所述发射场源为偶极源时，能够沿任意方向布设。
在其中一个实施例中，还包括：所述相互正交的电极可沿任意方向布设。
在其中一个实施例中，还包括：根据所述第一方向电位差和所述第二方向电位差，得到第一方向视电阻率和第二方向视电阻率，根据所述背景第一方向电位差和所述背景第二方向电位差，得到背景第一方向视电阻率和背景第二方向视电阻率的计算公式为：
其中，ρ为视电阻率，可分别为第一方向视电阻率ρx，第二方向视电阻率ρy，背景第一方向视电阻率背景第二方向视电阻率/>ΔU为电位差，可分别对应为背景第一方向视电阻率ΔUx，第二方向电位差ΔUy，背景第一方向电位差/>背景第二方向电位差K为装置系数，I为电流强度。
在其中一个实施例中，还包括：计算得到归一化视电阻率和归一化背景视电阻率的公式为：
其中，ρxy为归一化视电阻率，为归一化背景视电阻率。
在其中一个实施例中，还包括：根据所述归一化视电阻率和所述归一化背景视电阻率，得到测区地下介质的电阻率各向异性系数，包括：
根据所述归一化视电阻率和所述归一化背景视电阻率，得到测区地下介质的电阻率各向异性系数为：
其中，gxy表示地下介质的电阻率各向异性系数。
在其中一个实施例中，还包括：所述发射场源能够布设于地表或井中。
一种基于直流电法勘探获取地下介质电各向异性的装置，所述装置包括：
电位差观测模块，用于在测区任意位置布设发射场源后，采用相互正交的电极在测区内进行全区观测，得到不同测点之间的第一方向电位差和与第二方向电位差；所述发射场源为点源或偶极源；所述第一方向与所述第二方向相互正交；
归一化视电阻率确定模块，用于根据所述第一方向电位差和所述第二方向电位差，得到第一方向视电阻率和第二方向视电阻率，根据所述第一方向视电阻率和所述第二方向视电阻率，得到归一化视电阻率；
背景电位差计算模块，用于根据预设的地下介质背景电阻率模型，计算所述不同测点之间的背景第一方向电位差和背景第二方向电位差；
归一化背景视电阻率确定模块，用于根据所述背景第一方向电位差和所述背景第二方向电位差，得到背景第一方向视电阻率和背景第二方向视电阻率，根据所述背景第一方向视电阻率和所述背景第二方向视电阻率，得到归一化背景视电阻率；
电阻率各向异性系数确定模块，用于根据所述归一化视电阻率和所述归一化背景视电阻率，得到测区地下介质的电阻率各向异性系数。
一种计算机设备，包括存储器和处理器，所述存储器存储有计算机程序，所述处理器执行所述计算机程序时实现以下步骤：
在测区任意位置布设发射场源后，采用相互正交的电极在测区内进行全区观测，得到不同测点之间的第一方向电位差和与第二方向电位差；所述发射场源为点源或偶极源；所述第一方向与所述第二方向相互正交；
根据所述第一方向电位差和所述第二方向电位差，得到第一方向视电阻率和第二方向视电阻率，根据所述第一方向视电阻率和所述第二方向视电阻率，得到归一化视电阻率；
根据预设的地下介质背景电阻率模型，计算所述不同测点之间的背景第一方向电位差和背景第二方向电位差；
根据所述背景第一方向电位差和所述背景第二方向电位差，得到背景第一方向视电阻率和背景第二方向视电阻率，根据所述背景第一方向视电阻率和所述背景第二方向视电阻率，得到归一化背景视电阻率；
根据所述归一化视电阻率和所述归一化背景视电阻率，得到测区地下介质的电阻率各向异性系数。
一种计算机可读存储介质，其上存储有计算机程序，所述计算机程序被处理器执行时实现以下步骤：
在测区任意位置布设发射场源后，采用相互正交的电极在测区内进行全区观测，得到不同测点之间的第一方向电位差和与第二方向电位差；所述发射场源为点源或偶极源；所述第一方向与所述第二方向相互正交；
根据所述第一方向电位差和所述第二方向电位差，得到第一方向视电阻率和第二方向视电阻率，根据所述第一方向视电阻率和所述第二方向视电阻率，得到归一化视电阻率；
根据预设的地下介质背景电阻率模型，计算所述不同测点之间的背景第一方向电位差和背景第二方向电位差；
根据所述背景第一方向电位差和所述背景第二方向电位差，得到背景第一方向视电阻率和背景第二方向视电阻率，根据所述背景第一方向视电阻率和所述背景第二方向视电阻率，得到归一化背景视电阻率；
根据所述归一化视电阻率和所述归一化背景视电阻率，得到测区地下介质的电阻率各向异性系数。
上述基于直流电法勘探获取地下介质电各向异性的方法、装置、计算机设备和存储介质，通过在测区任意位置布设发射场源后，采用相互正交的电极在测区内进行全区观测，得到不同测点之间相互正交的两个方向的电位差，根据电位差计算视电阻率，之后计算归一化视电阻率，根据选取的背景电阻率模型，计算背景视电阻率，之后计算归一化背景视电阻率，根据实际归一化电阻率和归一化背景视电阻率，计算出电阻率各向异性系数。本发明布设的场源只需用点源或偶极源，布设灵活而简单。采用电阻率各向异性系数表征地下电各向异性的特性，可以有效压制视电阻率零线的影响，能很好地反应出地下各向异性结构体在地面投影以及形状。另外，可以避免求解张量电阻率的繁琐计算，提高技术效率。采用归一化视电阻率可以有效抵消场源附近视电阻率的误差，提升计算精度。
附图说明
图1为一个实施例中基于直流电法勘探获取地下介质电各向异性的方法的流程示意图；
图2为一个实施例中装置示意图；
图3为一个具体实施例中模型示意图；
图4为一个具体实施例中沿x方向测量得到的视电阻率ρx平面图；
图5为一个具体实施例中沿y方向测量得到的视电阻率ρy平面图；
图6为一个具体实施例中实例归一化电阻率各向异性系数gxy平面图；
图7为一个实施例中基于直流电法勘探获取地下介质电各向异性的装置的结构框图；
图8为一个实施例中计算机设备的内部结构图。
具体实施方式
为了使本申请的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本申请进行进一步详细说明。应当理解，此处描述的具体实施例仅仅用以解释本申请，并不用于限定本申请。
本申请提供的基于直流电法勘探获取地下介质电各向异性的方法，可以应用于如下应用环境中。其中，终端执行一种基于直流电法勘探获取地下介质电各向异性的方法，通过在测区任意位置布设发射场源后，采用相互正交的电极在测区内进行全区观测，得到不同测点之间相互正交的两个方向的电位差，根据电位差计算视电阻率，之后计算归一化视电阻率，根据选取的背景电阻率模型，计算背景视电阻率，之后计算归一化背景视电阻率，根据实际归一化电阻率和归一化背景视电阻率，计算出电阻率各向异性系数。其中，终端可以但不限于是各种个人计算机、笔记本电脑、平板电脑和便携式设备。
在一个实施例中，如图1所示，提供了一种基于直流电法勘探获取地下介质电各向异性的方法，包括以下步骤：
步骤102，在测区任意位置布设发射场源后，采用相互正交的电极在测区内进行全区观测，得到不同测点之间的第一方向电位差和与第二方向电位差。
第一方向与第二方向相互正交。发射场源为点源或偶极源，布设灵活而简单。
步骤104，根据第一方向电位差和第二方向电位差，得到第一方向视电阻率和第二方向视电阻率，根据第一方向视电阻率和第二方向视电阻率，得到归一化视电阻率。
采用归一化视电阻率可以有效抵消场源附近视电阻率的误差，提升计算精度。
步骤106，根据预设的地下介质背景电阻率模型，计算不同测点之间的背景第一方向电位差和背景第二方向电位差。
步骤108，根据背景第一方向电位差和背景第二方向电位差，得到背景第一方向视电阻率和背景第二方向视电阻率，根据背景第一方向视电阻率和背景第二方向视电阻率，得到归一化背景视电阻率。
步骤110，根据归一化视电阻率和归一化背景视电阻率，得到测区地下介质的电阻率各向异性系数。
采用电阻率各向异性系数表征地下电各向异性的特性，可以有效压制视电阻率零线的影响，能很好地反应出地下各向异性结构体在地面投影以及形状。另外，可以避免求解张量电阻率的繁琐计算，提高技术效率。
上述基于直流电法勘探获取地下介质电各向异性的方法中，通过在测区任意位置布设发射场源后，采用相互正交的电极在测区内进行全区观测，得到不同测点之间相互正交的两个方向的电位差，根据电位差计算视电阻率，之后计算归一化视电阻率，根据选取的背景电阻率模型，计算背景视电阻率，之后计算归一化背景视电阻率，根据实际归一化电阻率和归一化背景视电阻率，计算出电阻率各向异性系数。本发明布设的场源只需用点源或偶极源，布设灵活而简单。采用电阻率各向异性系数表征地下电各向异性的特性，可以有效压制视电阻率零线的影响，能很好地反应出地下各向异性结构体在地面投影以及形状。另外，可以避免求解张量电阻率的繁琐计算，提高技术效率。采用归一化视电阻率可以有效抵消场源附近视电阻率的误差，提升计算精度。
在其中一个实施例中，还包括：发射场源为点源或偶极源；当发射场源为偶极源时，能够沿任意方向布设。布设灵活而简单。
在其中一个实施例中，还包括：相互正交的电极可沿任意方向布设。布设灵活而简单。
在其中一个实施例中，还包括：根据第一方向电位差和第二方向电位差，得到第一方向视电阻率和第二方向视电阻率，根据背景第一方向电位差和背景第二方向电位差，得到背景第一方向视电阻率和背景第二方向视电阻率的计算公式为：
其中，ρ为视电阻率，可分别为第一方向视电阻率ρx，第二方向视电阻率ρy，背景第一方向视电阻率背景第二方向视电阻率/>ΔU为电位差，可分别对应为背景第一方向视电阻率ΔUx，第二方向电位差ΔUy，背景第一方向电位差/>背景第二方向电位差K为装置系数，I为电流强度。
在其中一个实施例中，还包括：计算得到归一化视电阻率和归一化背景视电阻率的公式为：
其中，ρxy为归一化视电阻率，为归一化背景视电阻率。
在其中一个实施例中，还包括：根据归一化视电阻率和归一化背景视电阻率，得到测区地下介质的电阻率各向异性系数，包括：
根据归一化视电阻率和归一化背景视电阻率，得到测区地下介质的电阻率各向异性系数为：
其中，gxy表示地下介质的电阻率各向异性系数。
在其中一个实施例中，还包括：发射场源能够布设于地表或井中。
应该理解的是，虽然图1的流程图中的各个步骤按照箭头的指示依次显示，但是这些步骤并不是必然按照箭头指示的顺序依次执行。除非本文中有明确的说明，这些步骤的执行并没有严格的顺序限制，这些步骤可以以其它的顺序执行。而且，图1中的至少一部分步骤可以包括多个子步骤或者多个阶段，这些子步骤或者阶段并不必然是在同一时刻执行完成，而是可以在不同的时刻执行，这些子步骤或者阶段的执行顺序也不必然是依次进行，而是可以与其它步骤或者其它步骤的子步骤或者阶段的至少一部分轮流或者交替地执行。
如图2为一个实施例中本发明的装置示意图，其中，M1和N1分别为用于测量第一方向电位差的两个电极，M2和N2分别为用于测量第二方向电位差的两个电极。M1N1的距离和M2N2的距离相等。
在一个具体实施例中，如图3所示为各向同性半空间中存在各向异性异常体模型示意图，其中异常体为棱柱体，它的中心点在地面上的投影为坐标系原点，棱柱大小为20m×20m×10m，地面埋深10m，三个主轴电阻率为ρ1/ρ2/ρ3＝10/5/10Ω·m，对应的欧拉角为α/β/χ＝30°/45°/60°，背景电阻率为ρ0＝100Ω·m。使用点源作为场源，位于坐标为(200m，200m，0)的位置。接收电极正交并分别测量x方向和y方向的电位(差)。图4为沿x方向测量得到的视电阻率ρx平面图。图5为本发明实例沿y方向测量得到的视电阻率ρy平面图。图6为归一化电阻率各向异性系数gxy平面图。图4、5、6中黑框表示异常体在地面上的投影位置和大小。由图4、5、6可知，使用归一化电阻率各向异性系数能更好地反应出各向异性异常体在地表投影的位置和形状。
在一个实施例中，如图7所示，提供了一种基于直流电法勘探获取地下介质电各向异性的装置，包括：电位差观测模块702、归一化视电阻率确定模块704、背景电位差计算模块706、归一化背景视电阻率确定模块708和电阻率各向异性系数确定模块710，其中：
电位差观测模块702，用于在测区任意位置布设发射场源后，采用相互正交的电极在测区内进行全区观测，得到不同测点之间的第一方向电位差和与第二方向电位差；第一方向与第二方向相互正交；
归一化视电阻率确定模块704，用于根据第一方向电位差和第二方向电位差，得到第一方向视电阻率和第二方向视电阻率，根据第一方向视电阻率和第二方向视电阻率，得到归一化视电阻率；
背景电位差计算模块706，用于根据预设的地下介质背景电阻率模型，计算不同测点之间的背景第一方向电位差和背景第二方向电位差；
归一化背景视电阻率确定模块708，用于根据背景第一方向电位差和背景第二方向电位差，得到背景第一方向视电阻率和背景第二方向视电阻率，根据背景第一方向视电阻率和背景第二方向视电阻率，得到归一化背景视电阻率；
电阻率各向异性系数确定模块710，用于根据归一化视电阻率和归一化背景视电阻率，得到测区地下介质的电阻率各向异性系数。
电阻率各向异性系数确定模块710还用于根据归一化视电阻率和归一化背景视电阻率，得到测区地下介质的电阻率各向异性系数为：
其中，gxy表示地下介质的电阻率各向异性系数。
关于基于直流电法勘探获取地下介质电各向异性的装置的具体限定可以参见上文中对于基于直流电法勘探获取地下介质电各向异性的方法的限定，在此不再赘述。上述基于直流电法勘探获取地下介质电各向异性的装置中的各个模块可全部或部分通过软件、硬件及其组合来实现。上述各模块可以硬件形式内嵌于或独立于计算机设备中的处理器中，也可以以软件形式存储于计算机设备中的存储器中，以便于处理器调用执行以上各个模块对应的操作。
在一个实施例中，提供了一种计算机设备，该计算机设备可以是终端，其内部结构图可以如图8所示。该计算机设备包括通过***总线连接的处理器、存储器、网络接口、显示屏和输入装置。其中，该计算机设备的处理器用于提供计算和控制能力。该计算机设备的存储器包括非易失性存储介质、内存储器。该非易失性存储介质存储有操作***和计算机程序。该内存储器为非易失性存储介质中的操作***和计算机程序的运行提供环境。该计算机设备的网络接口用于与外部的终端通过网络连接通信。该计算机程序被处理器执行时以实现一种基于直流电法勘探获取地下介质电各向异性的方法。该计算机设备的显示屏可以是液晶显示屏或者电子墨水显示屏，该计算机设备的输入装置可以是显示屏上覆盖的触摸层，也可以是计算机设备外壳上设置的按键、轨迹球或触控板，还可以是外接的键盘、触控板或鼠标等。
本领域技术人员可以理解，图8中示出的结构，仅仅是与本申请方案相关的部分结构的框图，并不构成对本申请方案所应用于其上的计算机设备的限定，具体的计算机设备可以包括比图中所示更多或更少的部件，或者组合某些部件，或者具有不同的部件布置。
在一个实施例中，提供了一种计算机设备，包括存储器和处理器，该存储器存储有计算机程序，该处理器执行计算机程序时实现上述方法实施例中的步骤。
在一个实施例中，提供了一种计算机可读存储介质，其上存储有计算机程序，计算机程序被处理器执行时实现上述方法实施例中的步骤。
本领域普通技术人员可以理解实现上述实施例方法中的全部或部分流程，是可以通过计算机程序来指令相关的硬件来完成，所述的计算机程序可存储于一非易失性计算机可读取存储介质中，该计算机程序在执行时，可包括如上述各方法的实施例的流程。其中，本申请所提供的各实施例中所使用的对存储器、存储、数据库或其它介质的任何引用，均可包括非易失性和/或易失性存储器。非易失性存储器可包括只读存储器(ROM)、可编程ROM(PROM)、电可编程ROM(EPROM)、电可擦除可编程ROM(EEPROM)或闪存。易失性存储器可包括随机存取存储器(RAM)或者外部高速缓冲存储器。作为说明而非局限，RAM以多种形式可得，诸如静态RAM(SRAM)、动态RAM(DRAM)、同步DRAM(SDRAM)、双数据率SDRAM(DDRSDRAM)、增强型SDRAM(ESDRAM)、同步链路(Synchlink)DRAM(SLDRAM)、存储器总线(Rambus)直接RAM(RDRAM)、直接存储器总线动态RAM(DRDRAM)、以及存储器总线动态RAM(RDRAM)等。
以上实施例的各技术特征可以进行任意的组合，为使描述简洁，未对上述实施例中的各个技术特征所有可能的组合都进行描述，然而，只要这些技术特征的组合不存在矛盾，都应当认为是本说明书记载的范围。
以上所述实施例仅表达了本申请的几种实施方式，其描述较为具体和详细，但并不能因此而理解为对发明专利范围的限制。应当指出的是，对于本领域的普通技术人员来说，在不脱离本申请构思的前提下，还可以做出若干变形和改进，这些都属于本申请的保护范围。因此，本申请专利的保护范围应以所附权利要求为准。

Claims (7)

1.一种基于直流电法勘探获取地下介质电各向异性的方法，其特征在于，所述方法包括：
在测区任意位置布设发射场源后，采用相互正交的电极在测区内进行全区观测，得到不同测点之间的第一方向电位差和与第二方向电位差；所述发射场源为点源或偶极源；所述第一方向与所述第二方向相互正交；
根据预设的地下介质背景电阻率模型，计算所述不同测点之间的背景第一方向电位差和背景第二方向电位差；
根据所述第一方向电位差和所述第二方向电位差，得到第一方向视电阻率和第二方向视电阻率，根据所述背景第一方向电位差和所述背景第二方向电位差，得到背景第一方向视电阻率和背景第二方向视电阻率；其中，得到第一方向视电阻率、第二方向视电阻率、背景第一方向视电阻率和背景第二方向视电阻率的计算公式为：
其中，ρ为视电阻率，可分别为第一方向视电阻率ρx，第二方向视电阻率ρy，背景第一方向视电阻率背景第二方向视电阻率/>ΔU为电位差，可分别对应为第一方向电位差ΔUx，第二方向电位差ΔUy，背景第一方向电位差/>背景第二方向电位差/>K为装置系数，I为电流强度；
根据所述第一方向视电阻率和所述第二方向视电阻率，得到归一化视电阻率，根据所述背景第一方向视电阻率和所述背景第二方向视电阻率，得到归一化背景视电阻率；其中，得到归一化视电阻率和归一化背景视电阻率的公式为：
其中，ρxy为归一化视电阻率，为归一化背景视电阻率；
根据所述归一化视电阻率和所述归一化背景视电阻率，得到测区地下介质的电阻率各向异性系数为：
其中，gxy表示地下介质的电阻率各向异性系数。
2.根据权利要求1所述的方法，其特征在于，当所述发射场源为偶极源时，能够沿任意方向布设。
3.根据权利要求1所述的方法，其特征在于，所述相互正交的电极可沿任意方向布设。
4.根据权利要求1至3任意一项所述的方法，其特征在于，所述发射场源能够布设于地表或井中。
5.一种基于直流电法勘探获取地下介质电各向异性的装置，其特征在于，所述装置包括：
电位差观测模块，用于在测区任意位置布设发射场源后，采用相互正交的电极在测区内进行全区观测，得到不同测点之间的第一方向电位差和与第二方向电位差；所述发射场源为点源或偶极源；所述第一方向与所述第二方向相互正交；
背景电位差计算模块，用于根据预设的地下介质背景电阻率模型，计算所述不同测点之间的背景第一方向电位差和背景第二方向电位差；
归一化视电阻率和归一化背景视电阻率确定模块，用于根据所述第一方向电位差和所述第二方向电位差，得到第一方向视电阻率和第二方向视电阻率，根据所述背景第一方向电位差和所述背景第二方向电位差，得到背景第一方向视电阻率和背景第二方向视电阻率；其中，得到第一方向视电阻率、第二方向视电阻率、背景第一方向视电阻率和背景第二方向视电阻率的计算公式为：
其中，ρ为视电阻率，可分别为第一方向视电阻率ρx，第二方向视电阻率ρy，背景第一方向视电阻率背景第二方向视电阻率/>ΔU为电位差，可分别对应为第一方向电位差ΔUx，第二方向电位差ΔUy，背景第一方向电位差/>背景第二方向电位差/>K为装置系数，I为电流强度；根据所述第一方向视电阻率和所述第二方向视电阻率，得到归一化视电阻率，根据所述背景第一方向视电阻率和所述背景第二方向视电阻率，得到归一化背景视电阻率；其中，得到归一化视电阻率和归一化背景视电阻率的公式为：
其中，ρxy为归一化视电阻率，为归一化背景视电阻率；
电阻率各向异性系数确定模块，用于根据所述归一化视电阻率和所述归一化背景视电阻率，得到测区地下介质的电阻率各向异性系数为：
其中，gxy表示地下介质的电阻率各向异性系数。
6.一种计算机设备，包括存储器和处理器，所述存储器存储有计算机程序，其特征在于，所述处理器执行所述计算机程序时实现权利要求1至4中任一项所述方法的步骤。
7.一种计算机可读存储介质，其上存储有计算机程序，其特征在于，所述计算机程序被处理器执行时实现权利要求1至4中任一项所述的方法的步骤。
CN202111157020.6A 2021-09-30 2021-09-30 基于直流电法勘探获取地下介质电各向异性的方法和装置 Active CN113917541B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111157020.6A CN113917541B (zh) 2021-09-30 2021-09-30 基于直流电法勘探获取地下介质电各向异性的方法和装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111157020.6A CN113917541B (zh) 2021-09-30 2021-09-30 基于直流电法勘探获取地下介质电各向异性的方法和装置

Publications (2)

Publication Number Publication Date
CN113917541A CN113917541A (zh) 2022-01-11
CN113917541B true CN113917541B (zh) 2024-02-23

Family

ID=79237506

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111157020.6A Active CN113917541B (zh) 2021-09-30 2021-09-30 基于直流电法勘探获取地下介质电各向异性的方法和装置

Country Status (1)

Country Link
CN (1) CN113917541B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116088059B (zh) * 2023-04-10 2023-07-18 山东省煤田地质规划勘察研究院 基于双源错频发射的人工源频率域电磁勘探方法及***

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101194262A (zh) * 2005-06-09 2008-06-04 埃克森美孚上游研究公司 在海洋电磁勘测中判断地球垂直电各向异性的方法
CN107742015A (zh) * 2017-09-30 2018-02-27 中南大学 基于任意偶极‑偶极装置的直流激电法三维数值模拟方法
CN110187394A (zh) * 2019-06-20 2019-08-30 甘肃省地震局（中国地震局兰州地震研究所） 双场源电磁测深法获取地层电阻率各向异性的方法及装置
CN111551997A (zh) * 2020-03-12 2020-08-18 上海环联生态科技有限公司 一种隐伏断裂层的勘察***及勘察方法
CN112709567A (zh) * 2019-10-25 2021-04-27 斯伦贝谢技术有限公司 用电磁测量值确定各向异性的地下性质
CN113238285A (zh) * 2021-05-08 2021-08-10 桂林理工大学 用于地球物理充电法勘探的电阻率计算方法及***、终端

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6765386B2 (en) * 2002-04-10 2004-07-20 Halliburton Energy Services, Inc. Galvanic method of measuring electrical anisotropy
US8581594B2 (en) * 2009-12-30 2013-11-12 Schlumberger Technology Corporation Microresistivity anisotropy logging tool employing a monopole current injection electrode

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101194262A (zh) * 2005-06-09 2008-06-04 埃克森美孚上游研究公司 在海洋电磁勘测中判断地球垂直电各向异性的方法
CN107742015A (zh) * 2017-09-30 2018-02-27 中南大学 基于任意偶极‑偶极装置的直流激电法三维数值模拟方法
CN110187394A (zh) * 2019-06-20 2019-08-30 甘肃省地震局（中国地震局兰州地震研究所） 双场源电磁测深法获取地层电阻率各向异性的方法及装置
CN112709567A (zh) * 2019-10-25 2021-04-27 斯伦贝谢技术有限公司 用电磁测量值确定各向异性的地下性质
CN111551997A (zh) * 2020-03-12 2020-08-18 上海环联生态科技有限公司 一种隐伏断裂层的勘察***及勘察方法
CN113238285A (zh) * 2021-05-08 2021-08-10 桂林理工大学 用于地球物理充电法勘探的电阻率计算方法及***、终端

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
普通电阻率测井复杂井眼环境校正模板计算正演方法研究;徐宏宇;《中国优秀硕士学位论文全文数据库基础科学辑》(2018年第05期);全文 *

Also Published As

Publication number Publication date
CN113917541A (zh) 2022-01-11

Similar Documents

Publication Publication Date Title
Wang et al. Three-dimensional DC anisotropic resistivity modelling using finite elements on unstructured grids
Lei et al. A numerical study of stress variability in heterogeneous fractured rocks
CN107742015B (zh) 基于任意偶极-偶极装置的直流激电法三维数值模拟方法
Mulchrone Application of Delaunay triangulation to the nearest neighbour method of strain analysis
Liu et al. Rock-soil slope stability analysis by two-phase random media and finite elements
Wang et al. CPT-based subsurface soil classification and zonation in a 2D vertical cross-section using Bayesian compressive sampling
CN113917541B (zh) 基于直流电法勘探获取地下介质电各向异性的方法和装置
Li et al. Inverse calculation of in situ stress in rock mass using the surrogate-model accelerated random search algorithm
Nevitt et al. Evaluation of transtension and transpression within contractional fault steps: Comparing kinematic and mechanical models to field data
Zhang et al. Three-dimensional magnetotelluric modeling using the finite element model reduction algorithm
Myrontsov The method to research equivalent solutions zones for inverse problem of well logging electrometry
Bai et al. 3D edge-based and nodal finite element modeling of magnetotelluric in general anisotropic media
Pang et al. 3D multi-scale resistivity inversion method applied in the tunnel face to borehole observations for tunnel-ahead prospecting
Roudsari et al. Forward modelling and inversion of self-potential anomalies caused by 2D inclined sheets
Tang et al. 2.5-D DC resistivity modeling considering flexibility and accuracy
Chen et al. Probabilistic back analysis for geotechnical engineering based on Bayesian and support vector machine
CN113076678B (zh) 一种频率域二度体重力异常快速数值模拟方法和装置
CN113221409B (zh) 一种有限元和边界元耦合的声波二维数值模拟方法和装置
CN110766794B (zh) 断层几何结构获取方法、装置和计算机设备
Xie et al. The behavior of vertically non-homogeneous elastic solids under internal rectangular loads
CN113268702A (zh) 一种频率域磁梯度张量变换方法、装置和计算机设备
Malecki et al. A new method for absolute underground positioning based on transient electromagnetics
CN113251989A (zh) 边坡变形的监测方法、装置及终端
Deng et al. Forward and inversion approach for direct current resistivity based on an unstructured mesh and its application to tunnel engineering
Prihantoro et al. Numerical solution of 3-D magnetotelluric using vector finite element method

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