CN107742015B - 基于任意偶极-偶极装置的直流激电法三维数值模拟方法 - Google Patents

基于任意偶极-偶极装置的直流激电法三维数值模拟方法 Download PDF

Info

Publication number
CN107742015B
CN107742015B CN201710913597.2A CN201710913597A CN107742015B CN 107742015 B CN107742015 B CN 107742015B CN 201710913597 A CN201710913597 A CN 201710913597A CN 107742015 B CN107742015 B CN 107742015B
Authority
CN
China
Prior art keywords
model
resistivity
electric field
dipole
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
CN201710913597.2A
Other languages
English (en)
Other versions
CN107742015A (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.)
Central South University
Original Assignee
Central South University
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 Central South University filed Critical Central South University
Priority to CN201710913597.2A priority Critical patent/CN107742015B/zh
Publication of CN107742015A publication Critical patent/CN107742015A/zh
Application granted granted Critical
Publication of CN107742015B publication Critical patent/CN107742015B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • 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

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Remote Sensing (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Theoretical Computer Science (AREA)
  • Geophysics (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Electromagnetism (AREA)
  • Geology (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • Environmental & Geological Engineering (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

本发明公开了一种基于任意偶极‑偶极装置的直流激电法三维数值模拟方法，包括采用三维直流激电任意偶极‑偶极装置全区观测方法及测量参数；实测数据归一化背景电阻率模型电场的计算；实测数据电阻率和极化率归一化传输函数的计算；理论模型电阻率和极化率归一化传输函数的计算。本发明方法能实现复杂条件下三维直流激电任意偶极‑偶极装置的全区观测三维数值模拟，而且不受电场零线限制，布极灵活，可以针对重点勘探区域任意加密测点密度，适应复杂环境的勘探；能够获取更为丰富的地下目标体的三维信息，并为各向异性介质的反演成像提供准备条件；而且本发明方法具有形态差异大、极化信息丰富和分辨率高的特点，能有效提升数值模拟的计算精度。

Description

基于任意偶极-偶极装置的直流激电法三维数值模拟方法
技术领域
本发明具体涉及一种基于任意偶极-偶极装置的直流激电法三维数值模拟方 法。
背景技术
地球物理学是通过定量的物理方法(如：地震弹性波、重力、地磁、地电、 地热和放射能等方法)研究地球以及寻找地球内部矿藏资源的一门综合性学科， 研究范围包括地球的地壳、地幔、地核和大气层等。近年来，地球物理学发展 越来越快，也越来越受到人们的重视。
直流激电法是地球物理学中一个非常重要的探测方法。直流激电法根据激 发和和观测条件的不同可以分为地面直流激电法、井中直流激电法和隧道直流 激电法，该方法作为一种重要的电法勘探方法，在矿产资源、水资源探测、考 古、环境和灾害监测、隧道超前预报、路基岩溶及煤田采空区探测等领域得到 了大量的应用，并取得了良好的实际效果。
直流激电法数值模拟方法与直流电阻率法完全相同，极化率数据的反演也 可以在电阻率反演的基础上得到。直流激电法与直流电阻率法一样，其观测方 式与数值模拟具有直接的联系，直流电法供电装置主要包括单极供电装置和偶 极供电装置，偶极供电装置测量数据由于极化信息丰富和分辨率高，获得了更 广泛的应用。但是在三维勘探中，由于受电场零线的影响，偶极-偶极装置很难 实现全区观测。考虑到野外施工环境和三维正反演方法的限制，基于任意布极 的三维测量面临极大的困难，目前主要包括如下两种观测方式：①采用以 E-SCAN为代表的规则三维测网进行电极的布设；②将多条二维测线组成三维测 网。
E-SCAN方法采用的是单极-单极(偶极)装置，将大量的电极连接到一根 多芯电缆上，由计算机对多路传输系统进行控制，采用电极转换的方式确定供 电电极和测量电极。其基本思想是将地面上的电极按照二维网格布置，每一个 网格节点上有一根电极，假设共有n根电极，则完成一次全测量所需的次数为 N＝n(n-1)。具体操作过程为：从m＝1(m为电极编号)开始为供电电极，其余电 极转换为测量电极，接着m＝2转换为供电电极，其余电极转换为测量电极，直 到m＝n时为止。这种测量方法当电极个数n过大时，对于单通道的二维测量仪 进行三维测量来说工作效率较低。为了提高野外施工工作效率，Jonathan(2005)针对已有的E-SCAN三维测量方式做了改进，采用“十字-对角方向测量”，即以 源点为中心以45°为间隔进行旋转，测量旋转方向上的测点，并通过模型和实测 资料进行了验证，这种测量方式与以往测测量方式几乎不降低地下目标体的信 息量。
但是，E-SCAN测网方式存在的缺点包括：
1)由于采用单极-单极或三极装置，很难向地下供大电流，只能够进行浅层 勘探，此外由于“无穷远”电极的存在，野外施工不方便；
2)所有电极都连接到同一根电缆上，造成野外施工不够灵活，布极要求非 常严格，在地形较为复杂的地区很难施工；
3)测量数据可能过于密集，造成施工资源的浪费。
而二维测线组成三维测网的施工方式虽然简单可行，但由于二维测线获取 地下目标体的三维信息有限，难以获得高精度成像结果。
发明内容
本发明的目的在于提供一种能够实现全区观测、适应性好、数据处理简单 可靠且准确的基于任意偶极-偶极装置的直流激电法三维数值模拟方法。
本发明提供的这种基于任意偶极-偶极装置的直流激电法三维数值模拟方 法，包括如下步骤：
S1.采用三维直流激电任意偶极-偶极装置对待测区域进行全区观测，从而 获得待测区域测点的水平电场
Figure BDA00014254050500000320
和
Figure BDA0001425405050000031
S2.根据步骤S1得到的待测区域测点的水平电场
Figure BDA0001425405050000032
和
Figure BDA0001425405050000033
计算得到实测数据电阻率合成电场
Figure BDA0001425405050000034
和激发极化合成电场
Figure BDA0001425405050000035
S3.计算步骤S2得到的实测数据的归一化背景电场
Figure BDA0001425405050000036
并计算实测数据 归一化传输函数
Figure BDA0001425405050000037
和极化率归一化传输函数
Figure BDA0001425405050000038
S4.根据步骤S1中观测系统参数对理论模型进行有限元数值模拟，获得理 论模型测点电阻率合成电场
Figure BDA0001425405050000039
和激发极化合成电场
Figure BDA00014254050500000310
S5.根据步骤S2中选取的背景电阻率模型和步骤S4中模型网格剖分参数， 计算理论模型归一化背景电场
Figure BDA00014254050500000311
S6.计算理论模型电阻率归一化传输函数
Figure BDA00014254050500000312
和极化率归一化传输函数
Figure BDA00014254050500000313
至此，分别获得了实测数据和理论模型归一化传输函数，从而为三维直流 激电归一化传输函数数据的反演提供条件。
步骤S1所述的采用三维直流激电任意偶极-偶极装置对待测区域进行全区 观测，具体为采用如下规则进行全区观测：
R1.发电电极和测量电极分开；
R2.发射电极为偶极装置，在待测区域内可以沿任意方向布设，且可以布设 于地表或井中；
R3.测量电极为偶极装置，测量方式为矢量观测方式，测点在待测区域内可 以任意布设，且在任意小区域内可以对测点进行加密或稀疏；
R4.矢量测量时，先后沿着x方向和y方向观测测点的水平电场
Figure BDA00014254050500000314
和
Figure BDA00014254050500000315
R5.对供电电极进行断电操作，观测激发极化电场
Figure BDA00014254050500000316
和
Figure BDA00014254050500000317
步骤S2所述的实测数据电阻率合成电场
Figure BDA00014254050500000318
和激发极化合成电场
Figure BDA00014254050500000319
具体为采用如下算式进行计算：
Figure BDA0001425405050000041
Figure BDA0001425405050000042
步骤S3所述的计算实测数据的归一化背景电场
Figure BDA0001425405050000043
具体为采用如下步 骤计算：
A.选取计算归一化背景场的电阻率模型；
B.采用网络加密-收缩算法对步骤A选定的电阻率模型进行网格离散，从 而提升归一化背景场的计算精度和效率；
C.根据步骤B对背景场电阻率模型进行网格离散后，按照野外观测系统参 数进行有限单元法数值模拟，计算得到实测数据归一化背景合成电场
Figure BDA0001425405050000044
并 计算实测数据电阻率归一化传输函数
Figure BDA0001425405050000045
和极化率归一化传输函数
Figure BDA0001425405050000046
具体采 用如下算式进行计算：
Figure BDA0001425405050000047
Figure BDA0001425405050000048
步骤S4所述的理论模型有限单元法数值模拟，具体为采用如下步骤进行计 算：
a.设计电阻率模型和极化率模型，并采用结构化网格对模型进行离散；
b.按照步骤a对模型离散后，根据野外观测系统参数采用有限单元法求解 三维直流激电微分方程组，计算出模型中每个节点的电位值；
c.根据步骤b得到的每个测点电位U后，采用下式计算得到测点水平电场 分量
Figure BDA0001425405050000049
和
Figure BDA00014254050500000410
E＝-▽U
式中▽为梯度运算，即电场E与电位U为负梯度关系。
d.根据步骤c计算得到的水平电场分量计算理论模型电阻率合成电场
Figure BDA0001425405050000051
和激发极化合成电场
Figure BDA0001425405050000052
具体为采用如下算式进行计算：
Figure BDA0001425405050000053
Figure BDA0001425405050000054
步骤S5所述的计算理论模型归一化背景电场，具体为采用如下步骤进行计 算：
y1.选取计算理论模型归一化背景场的电阻率模型，该电阻率模型与步骤 S3中实测数据归一化背景场电阻率模型相同；
y2.采用结构化网格对步骤y1中选取的电阻率模型进行网格离散，其网格 离散方式与步骤S4中模型离散方法完全相同；
y3.按照步骤y2对背景电阻率模型离散后，根据野外观测系统参数采用有 限单元法进行数值模拟，计算得到理论模型归一化背景电场
Figure BDA0001425405050000055
合成背景电 场的计算方式与步骤S4中相同。
步骤S6所述的计算理论模型电阻率归一化传输函数和极化率归一化传输函 数，具体为采用如下算式进行计算：
Figure BDA0001425405050000056
Figure BDA0001425405050000057
本发明提供的这种基于直流激电法的三维数值模拟方法，能实现复杂条件 下三维直流激电任意偶极-偶极装置全区观测三维数值模拟，并为该条件下的三 维反演成像提供准备条件；矢量测量方式能够有效地探测出地下目标体的三维 信息，并为以后各向异性介质的研究提供条件；偶极-偶极装置测量数据获取地 下目标体信息更为丰富，分辨率更高；不需要计算视电阻率，不受电场零线的 限制，其布极更为灵活，可以针对重点勘探区域任意加密测点密度，能够适应 复杂环境的勘探要求；在数值模拟中，由于采用了归一化处理，场源附近节点 误差抵消，大大提升了数值模拟的计算精度。
附图说明
图1为本发明方法的方法流程图。
图2为本发明中的任意偶极-偶极观测装置示意图。
图3为本发明实施算例中场源位置示意图。
图4为本发明实施例一中的双层介质隆起构造模型示意图。
图5为本发明实施例一中的不同背景模型正演数值模拟结果示意图。
图6为本发明实施例一中的不同激发方位的Rρ和Rη计算结果示意图。
图7为本发明实施例一中的归一化传输函数分辨率的示意图。
图8为本发明实施例二中的起伏地形模型的yz剖面图。
图9为本发明实施例二中的起伏地形条件下归一化传输函数平面图。
图10为本发明实施例二中的物性结构不一致模型的yz剖面图。
图11为本发明实施例二中的物性结构不一致时归一化传输函数的平面图。
具体实施方式
如图1所示为本发明的方法流程图：本发明提供的这种基于任意偶极-偶极 装置的直流激电法三维数值模拟方法，包括如下步骤：
S1.采用三维直流激电任意偶极-偶极装置对待测区域进行全区观测，从而 获得待测区域测点的水平电场
Figure BDA0001425405050000061
和
Figure BDA0001425405050000062
具体为采用如下 规则进行全区观测：
R1.发电电极和测量电极分开；
R2.发射电极为偶极装置，在待测区域内可以沿任意方向布设，且可以布设 于地表或井中；
R3.测量电极为偶极装置，测量方式为矢量观测方式，测点在待测区域内可 以任意布设，且在任意小区域内可以对测点进行加密或稀疏；
R4.矢量测量时，先后沿着x方向和y方向观测测点的水平电场
Figure BDA0001425405050000071
和
Figure BDA0001425405050000072
R5.对供电电极进行断电操作，观测激发极化电场
Figure BDA0001425405050000073
和
Figure BDA0001425405050000074
S2.根据步骤S1得到的待测区域测点的水平电场
Figure BDA0001425405050000075
和
Figure BDA0001425405050000076
计算得到实测数据电阻率合成电场
Figure BDA0001425405050000077
和激发极化合成电场
Figure BDA0001425405050000078
合 成电场的计算公式如下：
Figure BDA0001425405050000079
Figure BDA00014254050500000710
S3.计算步骤S2得到的实测数据的归一化背景电场
Figure BDA00014254050500000711
并计算实测数据 归一化传输函数
Figure BDA00014254050500000712
和极化率归一化传输函数
Figure BDA00014254050500000713
具体计算步骤如下：
A.选取计算归一化背景场的电阻率模型；当背景模型为简单模型(如均匀水 平半空间模型，水平层状介质模型)时，归一化背景电场可以根据解析表达式或 者数值滤波方法计算得到；当背景模型包含起伏地形或已知模型约束条件存在 时，采用数值计算获得合成电场数据。由于数值模拟受场源奇异性影响，但实 测数据没有这个问题，因此，实测数据归一化背景场值的计算至少需要满足如 下三个方面的条件：测点不受源点奇异性的影响、数值计算精度高(地表测点) 以及计算效率快；
B.采用网络加密-收缩算法对步骤A选定的电阻率模型进行网格离散，从 而提升归一化背景场的计算精度和效率；具体为：采用网格加密-收缩算法对背 景模型进行网格离散，该算法主要分为两个步骤：
(1)在常规结构网格基础上均匀加密计算区域中水平方向网格节点密度， 从而降低场源影响范围，提升背景场计算精度；
(2)在垂直方向采用网格收缩方法降低网格节点的规模数，从而提升背景 场计算效率；
上述的网格加密-收缩算法可以参见论文：Finite element numericalsimulation of 2.5D direct current method based on mesh refinement andrecoarsement.Applied Geophysics,2016,13(2):257-266。
C.根据步骤B对背景场电阻率模型进行网格离散后，按照野外观测系统参 数进行有限单元法数值模拟，计算得到实测数据归一化背景合成电场
Figure BDA0001425405050000081
并 分别计算实测数据电阻率归一化传输函数
Figure BDA0001425405050000082
和极化率归一化传输函数
Figure BDA0001425405050000083
具 体采用如下算式进行计算：
Figure BDA0001425405050000084
Figure BDA0001425405050000085
S4.根据步骤S1中观测系统参数对理论模型进行有限元数值模拟，获得理 论模型测点电阻率合成电场
Figure BDA0001425405050000086
和激发极化合成电场
Figure BDA0001425405050000087
具体为采用如下 步骤进行计算：
a.设计电阻率模型和极化率模型，并采用结构化网格对模型进行离散；
b.按照步骤a对模型离散后，根据野外观测系统参数采用有限单元法求解 三维直流激电微分方程组，计算出模型中每个节点的电位值；
c.根据步骤b得到的每个测点电位U后，采用下式计算得到测点水平电场 分量
Figure BDA0001425405050000088
和
Figure BDA0001425405050000089
E＝-▽U
式中▽为梯度运算，即电场E与电位U为负梯度关系。
d.根据步骤c计算得到的水平电场分量计算理论模型电阻率合成电场
Figure BDA00014254050500000810
和激发极化合成电场
Figure BDA00014254050500000811
具体为采用如下算式进行计算：
Figure BDA00014254050500000812
Figure BDA00014254050500000813
S5.根据步骤S2中选取的背景电阻率模型和步骤S4中模型网格剖分参数， 计算理论模型归一化背景电场
Figure BDA0001425405050000091
具体计算步骤如下：
y1.选取计算理论模型归一化背景场的电阻率模型，该电阻率模型与步骤 S3中实测数据归一化背景场电阻率模型相同；
y2.采用结构化网格对步骤y1中选取的电阻率模型进行网格离散，其网格 离散方式与步骤S4中模型离散方法完全相同；
y3.按照步骤y2对背景电阻率模型离散后，根据野外观测系统参数采用有 限单元法进行数值模拟，计算得到理论模型归一化背景电场
Figure BDA0001425405050000092
合成背景电 场的计算方式与步骤S4中相同。
S6.计算理论模型电阻率归一化传输函数
Figure BDA0001425405050000093
和极化率归一化传输函数
Figure BDA0001425405050000094
至此，分别获得了实测数据和理论模型归一化传输函数，从而为三维直流 激电归一化传输函数数据的反演提供条件。步骤S6中理论模型归一化传输函数 计算方式如下：
Figure BDA0001425405050000095
Figure BDA0001425405050000096
从上式可以看出：电阻率归一化传输函数需要选取背景模型，而单独计算 极化率归一化传输函数时对背景模型和网格剖分不做特殊要求。对比步骤S3中 实测数据归一化传输函数可以看出：电阻率归一化传输函数计算中最大的区别 在于背景电场的不同，实测数据背景电场是精确的，而理论模型背景电场在场 源附近存在数值误差，但因为理论模型在场源附近同样存在数值误差，并且这 种误差具有相关性，因此归一化处理后能消除场源奇异性的影响。
以下结合两个实施例对本发明方法进行进一步说明：
如图3所示为模型算例中标量场源布设位置图，S1-S5为5对偶极子源，其 中S1﹑S2﹑S3和S5布设于地表，S1﹑S2和S3沿y方向布设，S5沿x方向布 设，偶极距为80m；S4布设于井中，位于隆起构造下部，埋深400m，偶极距 为60m。在2km×2km的测区范围内，矢量测量测点可以任意布设，以计算区域 中101×101节点作为测点，测点间距为20m。
实施例一：水平地形下隆起构造模型的三维IP正演
构造三维双层介质隆起构造电阻率模型，第一层介质电阻率为40Ω·m，极 化率为0，厚度为100m；第二层介质电阻率为100Ω·m，极化率为0；隆起构造 距离顶界面为150m，厚度为100m，顶界面沿x方向走向200m，沿y方向走 向360m,构造为两层介质：上层介质电阻率为10Ω·m，极化率为20％，厚度为 40m，下层介质电阻率为100Ω·m，极化率为0。
如图4所示为双层介质隆起构造模型剖面图，其中，图4中(a)为yz剖面 图，图中隆起构造呈梯形，顶宽360m，底宽760m，厚度为100m；图4中(b) 为xz剖面图，图中隆起构造呈矩形，宽为200m，厚度为100m。
背景模型对电阻率归一化传输函数的影响：
为了研究背景模型的选取对电阻率归一化传输函数的影响，选取了如下5 个模型计算归一化背景电场
Figure BDA0001425405050000101
背景模型一，电阻率为50Ω·m的均匀半空间模型；
背景模型二，电阻率为40Ω·m的均匀半空间模型；
背景模型三，算例中的双层介质电阻率模型(不含隆起构造)；
背景模型四，算例中含隆起构造的双层介质电阻率模型，隆起构造中电阻 率均为100Ω·m。
分别计算不同背景模型时电阻率归一化传输函数：
Figure BDA0001425405050000102
其中
Figure BDA0001425405050000111
为隆起构造中含低阻异常体的电阻率模型计算结果；
如图5所示为偶极源位于地表(S1)时正演数值模拟结果，从图中可以看出：
1)本文定义的电阻率归一化传输函数在整个测区形态光滑,不受传统视电阻 率计算中电场零线的影响；
2)即使背景模型与实际模型存在差异时(如
Figure BDA0001425405050000112
中背景模型给的50Ω·m的均 匀半空间模型，而实际模型第一层介质为40Ω·m)，源点附近的计算误差均能在 归一化中被抵消；
3)在
Figure BDA0001425405050000113
中，由于背景模型电阻率与真实模型第一层介质电阻率相同，且偶 极子源远离异常体，因此源点附近区域测点不受地下高阻层的影响，归一化传 输函数数值为1，距离源点较远的测点，受地下高阻层介质和隆起构造的影响， 归一化传输函数数值逐渐增大。对比
Figure BDA0001425405050000114
可以发现，两个模型计算结果整体形态 相似，但由于背景模型较真实模型第一层介质电阻率更高，在源点附近归一化 传输函数数值低于1；
背景模型越接近真实模型，电阻率归一化传输函数中异常体特征越明显， 且异常差异更大(如
Figure BDA0001425405050000115
和
Figure BDA0001425405050000116
)，从而为具有一定地质信息的测区采用约束反演提 供条件。
场源位置对归一化传输函数的影响：
分别对观测系统(图3)中4对偶极子源(S1-S4)进行正演数值模拟，并计算电 阻率归一化传输函数Rρ和极化率归一化传输函数Rη，其中计算Rρ时背景模型为 双层介质模型(上述算例中模型三)。令
Figure BDA0001425405050000117
为隆起构造中含低阻异常体和极化率的 等效电阻率模型计算结果，可得计算表达式:
Figure BDA0001425405050000118
其中k＝1,2,3,4分别对应偶极源S1-S4。
计算结果如图6所示：
1)场源位置发生改变时，由于激发方位和照射角度的不同，归一化传输函 数具有形态差异很大﹑极化信息丰富和分辨率高的特点；
2)电阻率异常和极化率异常同源时，Rρ和Rη具有相似的形态分布图；
极化率归一化传输函数直接反映了地下极化异常体的形态和在地表的投影 位置，因此，通过对大量理论模型的数值结果进行研究,总结出规律性认识后， 可以根据归一化传输函数形态分布图大致推断出地下目标体和极化体的构造特 征。
归一化传输函数分辨率：
如图3所示，隆起构造在x和y方向长度不等，沿x走向为200m，沿y走 向为360m。因此，本文分别对观测系统(图3)中偶极源S1和S5进行正演数值 模拟，通过对比两次正演计算结果分析归一化传输函数的水平分辨率。算例中， S1沿y方向布设，S5沿x方向布设,两对偶极子源在空间上的关系为S5为S1 沿逆时针旋转90°，偶极距均为80m，供电电流强度均为20A。
计算结果如图7所示，图中白色方框对应隆起构造中低阻异常体和高极化 率体在地表的投影位置。电阻率归一化传输函数计算中背景模型为双层隆起构 造模型(背景模型四)，从而使得在计算Rρ和Rη时电阻率异常体和极化率异常体 同源.从图中可以看出：当偶极子源沿y方向布设时(S1)，归一化传输函数在y 方向分辨率较高，y方向主异常形态与异常体在空间位置上一一对应；当偶极子 源沿x方向布设时(S2)，归一化传输函数在x方向分辨率较高，x方向异常形态 与异常体在空间位置上一一对应。综上可以推断出归一化传输函数沿偶极子源 走向方向具有很高的分辨率。在实际勘探过程中,综合地下目标体的构造走向， 采用偶极子源进行多方位和多角度照射，可以大大提升勘探效果。
实施例二：复杂起伏地形条件下三维IP正演
设计如图8所示(模型yz剖面图)的三维凹嵌地形构造电阻率模型，模型电 阻率为100Ω·m，极化率为1％；凹嵌地形距离水平地表30m，地形底端为正方 向，边长为800m；在地形下方置入一电阻率为20Ω·m﹑极化率为20％的低阻 异常体，异常体距离顶界面70m，厚度为40m，沿x和y方向长度为200m。
地形对归一化传输函数的影响：
以背景电阻率为100Ω·m﹑极化率为0％的凹嵌地形模型为背景场计算模 型，计算归一化电阻率传输函数Rρ和极化率归一化传输函数Rη，其中偶极源S1 位于地表(如图3所示)。计算结果如图9所示，对比水平地形算例计算结果可以 发现：起伏地形对Rρ和Rη没有任何影响，从而可以得出结论：归一化传输函数 既可以解决源点奇异性问题又能抵消地形对数值计算结果的影响，以此为基础 的三维正反演方法可以适用于复杂环境勘探中。
物性结构不一致：
对图8所示的极化率模型做如下改造：取模型低阻异常体中极化率等于背 景模型极化率(1％)，沿y轴对称的方向置入一个与低阻异常体大小﹑埋深相同 的高极化体，极化率为20％，电阻率为100Ω·m。模型yz剖面图如图10所示。
采用与上一个算例中相同的计算方法：以背景电阻率为100Ω·m﹑极化率为 0％的凹嵌地形模型为背景场计算模型，计算归一化电阻率传输函数Rρ和极化率 归一化传输函数Rη，其中偶极子源S5位于地表(如图3所示)。计算结果如图11 所示：当物性结构不一致时，Rρ和Rη平面图均能反映出电阻率异常体和极化率 异常体在地表的投影位置；对比图9可以发现：当异常体为低阻高极化时，极 化率归一化传输函数Rη具有更大的异常数值，这与直流电阻率法对低阻体更敏 感的特性相符。

Claims (5)

1.一种基于任意偶极-偶极装置的直流激电法三维数值模拟方法，包括如下步骤：
S1.采用三维直流激电任意偶极-偶极装置对待测区域进行全区观测，从而获得待测区域测点的水平电场
Figure FDA0002757863990000011
和
Figure FDA0002757863990000012
S2.根据步骤S1得到的待测区域测点的水平电场
Figure FDA0002757863990000013
和
Figure FDA0002757863990000014
计算得到实测数据电阻率合成电场
Figure FDA0002757863990000015
和激发极化合成电场
Figure FDA0002757863990000016
S3.计算步骤S2得到的实测数据的归一化背景电场
Figure FDA0002757863990000017
并计算实测数据归一化传输函数
Figure FDA0002757863990000018
和极化率归一化传输函数
Figure FDA0002757863990000019
具体为采用如下步骤计算：
A.选取计算归一化背景场的电阻率模型；
B.采用网络加密-收缩算法对步骤A选定的电阻率模型进行网格离散，从而提升归一化背景场的计算精度和效率；
C.根据步骤B对背景场电阻率模型进行网格离散后，按照野外观测系统参数进行有限单元法数值模拟，计算得到实测数据归一化背景合成电场
Figure FDA00027578639900000110
并计算实测数据电阻率归一化传输函数
Figure FDA00027578639900000111
和极化率归一化传输函数
Figure FDA00027578639900000112
具体采用如下算式进行计算：
Figure FDA00027578639900000113
Figure FDA00027578639900000114
S4.根据步骤S1中观测系统参数对理论模型进行有限元数值模拟，获得理论模型测点电阻率合成电场
Figure FDA00027578639900000115
和激发极化合成电场
Figure FDA00027578639900000116
具体为采用如下步骤进行计算：
a.设计电阻率模型和极化率模型，并采用结构化网格对模型进行离散；
b.按照步骤a对模型离散后，根据野外观测系统参数采用有限单元法求解三维直流激电微分方程组，计算出模型中每个节点的电位值；
c.根据步骤b得到的每个节点的电位值，可以得到每个测点的电位值U，然后采用下式计算得到测点水平电场分量
Figure FDA0002757863990000021
和
Figure FDA0002757863990000022
E＝-▽U
式中▽为梯度运算；
d.根据步骤c计算得到的水平电场分量计算理论模型电阻率合成电场
Figure FDA0002757863990000023
和激发极化合成电场
Figure FDA0002757863990000024
具体为采用如下算式进行计算：
Figure FDA0002757863990000025
Figure FDA0002757863990000026
S5.根据步骤S3中的背景场电阻率模型和步骤S4中的采用结构化网格对模型进行离散的方法，计算理论模型归一化背景电场
Figure FDA0002757863990000027
S6.计算理论模型电阻率归一化传输函数
Figure FDA0002757863990000028
和极化率归一化传输函数
Figure FDA0002757863990000029
从而获得了实测数据和理论模型归一化传输函数，为三维直流激电归一化传输函数数据的反演提供条件。
2.根据权利要求1所述的基于任意偶极-偶极装置的直流激电法三维数值模拟方法，其特征在于步骤S1所述的采用三维直流激电任意偶极-偶极装置对待测区域进行全区观测，具体为采用如下规则进行全区观测：
R1.发电电极和测量电极分开；
R2.发射电极为偶极装置，在待测区域内可以沿任意方向布设，且可以布设于地表或井中；
R3.测量电极为偶极装置，测量方式为矢量观测方式，测点在待测区域内可以任意布设，且在任意小区域内可以对测点进行加密或稀疏；
R4.矢量测量时，先后沿着x方向和y方向观测测点的水平电场
Figure FDA00027578639900000210
和
Figure FDA0002757863990000031
R5.对供电电极进行断电操作，观测激发极化电场
Figure FDA0002757863990000032
和
Figure FDA0002757863990000033
3.根据权利要求2所述的基于任意偶极-偶极装置的直流激电法三维数值模拟方法，其特征在于步骤S2所述的实测数据电阻率合成电场
Figure FDA0002757863990000034
和激发极化合成电场
Figure FDA0002757863990000035
具体为采用如下算式进行计算：
Figure FDA0002757863990000036
Figure FDA0002757863990000037
4.根据权利要求3所述的基于任意偶极-偶极装置的直流激电法三维数值模拟方法，其特征在于步骤S5所述的计算理论模型归一化背景电场，具体为采用如下步骤进行计算：
y1.选取计算理论模型归一化背景场的电阻率模型；
y2.采用结构化网格对步骤y1中选取的电阻率模型进行网格离散；
y3.按照步骤y2对背景电阻率模型离散后，根据野外观测系统参数采用有限单元法进行数值模拟，计算得到理论模型归一化背景电场
Figure FDA0002757863990000038
5.根据权利要求4所述的基于任意偶极-偶极装置的直流激电法三维数值模拟方法，其特征在于步骤S6所述的计算理论模型电阻率归一化传输函数和极化率归一化传输函数，具体为采用如下算式进行计算：
Figure FDA0002757863990000039
Figure FDA00027578639900000310
CN201710913597.2A 2017-09-30 2017-09-30 基于任意偶极-偶极装置的直流激电法三维数值模拟方法 Active CN107742015B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710913597.2A CN107742015B (zh) 2017-09-30 2017-09-30 基于任意偶极-偶极装置的直流激电法三维数值模拟方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710913597.2A CN107742015B (zh) 2017-09-30 2017-09-30 基于任意偶极-偶极装置的直流激电法三维数值模拟方法

Publications (2)

Publication Number Publication Date
CN107742015A CN107742015A (zh) 2018-02-27
CN107742015B true CN107742015B (zh) 2021-04-23

Family

ID=61236401

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710913597.2A Active CN107742015B (zh) 2017-09-30 2017-09-30 基于任意偶极-偶极装置的直流激电法三维数值模拟方法

Country Status (1)

Country Link
CN (1) CN107742015B (zh)

Families Citing this family (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110609331A (zh) * 2019-09-24 2019-12-24 中国地质调查局西安地质调查中心 一种隐伏沉积型富锰矿识别方法
CN110703345B (zh) * 2019-11-05 2022-01-18 山东科技大学 一种环工作面对角偶极装置三维直流电法数据采集方法
CN111722291B (zh) * 2020-06-08 2023-05-05 重庆交通大学 一种以闪电为刺激源的大尺度地质构造三维层析扫描方法
CN111767669A (zh) * 2020-07-08 2020-10-13 湖南省有色地质勘查研究院 一种新式伪随机激电法有限元数值模拟方法及系统
CN112630842B (zh) * 2020-10-14 2024-04-12 中国地质调查局南京地质调查中心(华东地质科技创新中心) 一种应用于直流电阻率法勘探的偶极差分装置及方法
CN112415602B (zh) * 2020-10-15 2021-11-19 山东大学 基于深度分辨率的隧道电阻率超前探测优化方法及系统
CN113031080B (zh) * 2021-04-21 2022-05-13 中南大学 一种采用简易极化常数的新激电方法
CN113204905B (zh) * 2021-05-08 2023-06-20 桂林理工大学 一种接触式激发极化法有限单元数值模拟方法
CN113221411B (zh) * 2021-05-08 2022-12-09 桂林理工大学 针对任意形状有损耗介质充电电位数值模拟方法、系统、终端
CN113238285B (zh) * 2021-05-08 2023-05-09 桂林理工大学 用于地球物理充电法勘探的电阻率计算方法及系统、终端
CN113051779B (zh) * 2021-05-31 2021-08-10 中南大学 一种三维直流电阻率法数值模拟方法
CN113391362B (zh) * 2021-08-13 2021-10-29 成都理工大学 基于廊带数据约束的大地电磁剖面三维结构化反演方法
CN113917541B (zh) * 2021-09-30 2024-02-23 湖南科技大学 基于直流电法勘探获取地下介质电各向异性的方法和装置
CN118068445B (zh) * 2024-04-18 2024-07-09 山东大学 一种时间域激发极化法极化率参数标准化方法及系统

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2010074593A1 (ru) * 2008-12-22 2010-07-01 Rykhlinskaya Ekaterina Nikolaevna Способ морской геоэлектроразведки с фокусировкой электрического тока
CN103955000A (zh) * 2014-04-02 2014-07-30 中国冶金地质总局山东正元地质勘查院 一种三维电测深方法
CN105893678A (zh) * 2016-04-01 2016-08-24 吉林大学 一种时域有限差分的三维感应-极化双场数值模拟方法
CN107191181A (zh) * 2017-05-24 2017-09-22 中国石油大学(华东) 一种基于电磁散射的井周界面探测方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2010074593A1 (ru) * 2008-12-22 2010-07-01 Rykhlinskaya Ekaterina Nikolaevna Способ морской геоэлектроразведки с фокусировкой электрического тока
CN103955000A (zh) * 2014-04-02 2014-07-30 中国冶金地质总局山东正元地质勘查院 一种三维电测深方法
CN105893678A (zh) * 2016-04-01 2016-08-24 吉林大学 一种时域有限差分的三维感应-极化双场数值模拟方法
CN107191181A (zh) * 2017-05-24 2017-09-22 中国石油大学(华东) 一种基于电磁散射的井周界面探测方法

Also Published As

Publication number Publication date
CN107742015A (zh) 2018-02-27

Similar Documents

Publication Publication Date Title
CN107742015B (zh) 基于任意偶极-偶极装置的直流激电法三维数值模拟方法
Zhou Electrical resistivity tomography: a subsurface-imaging technique
CN110058315B (zh) 一种三维各向异性射频大地电磁自适应有限元正演方法
Berge et al. Electrical resistivity tomography investigations of multilayered archaeological settlements: part I–modelling
Yin et al. 3D time-domain airborne EM forward modeling with topography
Papadopoulos et al. Electrical resistivity tomography for the modelling of cultural deposits and geomophological landscapes at Neolithic sites: a case study from Southeastern Hungary
CN103728667A (zh) 一种视三维高密度电法的浅表层地质结构建模方法
Ingeman-Nielsen et al. Effect of electrode shape on grounding resistances—Part 1: The focus-one protocol
Lu et al. A systematic study of topography effect of ERT based on 3-D modeling and inversion
CN111638556A (zh) 基于地空分解策略的大地电磁正演方法及装置、存储介质
Mukanova et al. The boundary element method in geophysical survey
Akça et al. Extraction of structure-based geoelectric models by hybrid genetic algorithms
Aning et al. Electrical resistivity as a geophysical mapping tool; a case study of the new art department, Knust-Ghana
CN112666612B (zh) 基于禁忌搜索的大地电磁二维反演方法
CN107748393B (zh) 一种基于数值模拟的地层倾角对电阻率影响的校正方法
Fargier et al. 3D assessment of an underground mine pillar by combination of photogrammetric and geoelectric methods
JP2005337746A (ja) 電気探査方法
CN107797148B (zh) 一种基于三维地质建模的航磁异常场分离方法及系统
Tang et al. 2.5-D DC resistivity modeling considering flexibility and accuracy
Günther et al. Electrical Resistivity Tomography (ERT) in geophysical applications-state of the art and future challenges
Mauriello et al. A data-adaptive probability-based fast ERT inversion method
Yan et al. Three-dimensional distributed DC/IP method for altered tectonite-type gold ore deposits exploration: a case study of the Jiaojia gold metallogenic belt, Eastern China
Aizebeokhai et al. Anomaly effects of orthogonal paired-arrays for 3D geoelectrical resistivity imaging
Zhang et al. Electrode array and data density effects in 3D induced polarization tomography and applications for mineral exploration
CN114065577A (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