CN110208868A - 一种基于联合剖面的二维电法勘探方法 - Google Patents

一种基于联合剖面的二维电法勘探方法 Download PDF

Info

Publication number
CN110208868A
CN110208868A CN201910386820.1A CN201910386820A CN110208868A CN 110208868 A CN110208868 A CN 110208868A CN 201910386820 A CN201910386820 A CN 201910386820A CN 110208868 A CN110208868 A CN 110208868A
Authority
CN
China
Prior art keywords
point
supply terminals
composite profile
record
data
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.)
Pending
Application number
CN201910386820.1A
Other languages
English (en)
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 CN201910386820.1A priority Critical patent/CN110208868A/zh
Publication of CN110208868A publication Critical patent/CN110208868A/zh
Pending legal-status Critical Current

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/38—Processing data, e.g. for analysis, for interpretation, for correction

Landscapes

  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Environmental & Geological Engineering (AREA)
  • Geology (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

本发明提供一种基于联合剖面的二维电法勘探方法。本发明基于电法勘探中的联合剖面法获取某测线的所有记录点的两个方向的电阻率、极化率等电法数据，对相同记录点、两个方向、相同极距的电阻率、极化率等电法数据分别进行平均值和差值的绝对值等综合处理，主要选择相关电法参数的差值的绝对值的极大值点进行异常体的异常宽度判断，并进行综合解译，从而在不增加野外工作量的前提下，提高联合剖面法的勘探效果。本发明能提高联合剖面法的电法勘探效果和精度，增加对地质异常体的识别度。

Description

一种基于联合剖面的二维电法勘探方法
技术领域
本发明属于电法勘探领域，涉及一种基于联合剖面的二维电法勘探方法。
背景技术
目前电法勘探中的联合剖面对于判断地下地质异常体的产状效果较好，在电法勘探中得到较多的应用，且对于相同测点的相同极距的视电阻率、视极化率等电法参数可以在理论上获得对称四极剖面数据，从而能进一步提高对地下地质异常体的判断精度。但基于目前的联合剖面电法数据处理方式，对于地质异常体的分布范围特别是异常体的宽度判断精度不高，故有必要对联合剖面所获取的视电阻率、视极化率等电法参数做进一步改进，提高对地质异常体的分布范围特别是异常体的宽度的判断精度。
发明内容
本发明的目的是针对现有技术存在的问题，提出一种基于联合剖面的二维电法勘探方法，以便提高联合剖面的电法勘探效果和精度。
为实现以上目的，本发明通过下列技术步骤来实现：
a)按照常规的联合剖面方法获取某条测线的多个记录点的两个不同方向的联合剖面电法数据；
b)对相同记录点的两个不同方向的相同供电极距的视电阻率数据，按照公式计算两个方向的视电阻率数据的平均值，按照公式计算两个方向的视电阻率数据的差值的绝对值；其中Oi为记录点的编号，Ai和Bi为对应记录点Oi的供电点编号，i为自然数；代表Ai供电点供电时，Oi记录点的联合剖面视电阻率；代表Bi供电点供电时，Oi记录点的联合剖面视电阻率；代表Oi记录点的Ai和Bi供电点的联合剖面两个方向的视电阻率平均值；代表Oi记录点的Ai和Bi供电点的联合剖面两个方向的视电阻率差值的绝对值；ABS代表对取绝对值；其中Ai供电点到Oi记录点的距离等于Bi供电点到Oi记录点的距离，且Ai供电点与Bi供电点相对于Oi记录点的方向相反；
c)对相同记录点的两个不同方向的相同供电极距的视极化率数据，按照公式计算两个方向的视极化率数据的平均值，按照公式计算两个方向的视极化率数据的差值的绝对值；其中Oi为记录点的编号，Ai和Bi为对应记录点Oi的供电点编号，i为自然数；代表Ai供电点供电时，Oi记录点的联合剖面视极化率；代表Bi供电点供电时，Oi记录点的联合剖面视极化率；代表Oi记录点的Ai和Bi供电点的联合剖面两个方向的视极化率平均值；代表Oi记录点的Ai和Bi供电点的联合剖面两个方向的视极化率差值的绝对值；ABS代表对取绝对值；其中Ai供电点到Oi记录点的距离等于Bi供电点到Oi记录点的距离，且Ai供电点与Bi供电点相对于Oi记录点的方向相反；
d)按照相同电法参数类型、相同供电极距、相同测线的要求，绘制曲线；把同一条测线所有记录点的相同供电极距的所有电法数据分别绘制曲线，每条曲线所有记录点的的电法参数类型、供电极距、测线相同，纵坐标为电法数据值，横坐标为电法数据的记录点；
e)对同一条测线所有记录点的不同电法曲线进行综合解译，利用曲线主要判断电法异常体的电法特性，利用或曲线的极大值点主要判断地质异常体的分布范围，利用或曲线的异常段数值判断地质异常体的电阻率或极化率特征，从而解决地质异常体的埋藏位置、分布范围等问题。
以上方法步骤中的等联合剖面的原始数据均按照常规的联合剖面方法获得，Ai、Bi供电点分别位于Oi记录点的两侧，且均位于同一条测线上。
附图说明
图1为本发明的方法流程图；
图2为应用本发明的高阻异常体的分布示意图；
图3为采用本发明处理在高阻异常体上联合剖面的与电阻率相关的电法参数曲线对比图(异常体对应图2)；
图4为采用本发明处理在高阻异常体上联合剖面的电阻率平均值和差值绝对值的电法参数曲线对比图(异常体对应图2)；
图5为采用本发明处理在低阻异常体上联合剖面的与电阻率相关的电法参数曲线对比图；
图6为采用本发明处理在低阻异常体上联合剖面的电阻率平均值和差值绝对值的电法参数曲线对比图；
图中RES-A代表Ai供电点供电时记录点上所获取的视电阻率数据(即)所组成的曲线，且同一条曲线中的每个记录点的供电极距相同；RES-B代表Bi供电点供电时记录点上所获取的视电阻率数据(即)所组成的曲线，且同一条曲线中的每个记录点的供电极距相同；RES/2代表Ai和Bi供电点的相同记录点的联合剖面视电阻率平均值(即)所组成的曲线，且同一条曲线中的每个记录点的供电极距相同；RES-ABS代表Ai和Bi供电点的相同记录点的联合剖面视电阻率差值的绝对值(即 )所组成的曲线，且同一条曲线中的每个记录点的供电极距相同。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚，下面将结合图1、图2、图3、图4、图5、图6和实施例对本发明实施方式作进一步详细说明。
实施例1
如图2所示，假设要对20个长度单位的测线开展联合剖面，其在9～12个长度单位位置下面埋藏一个长度为3个长度单位、顶面埋深1个长度单位、厚度为3个长度单位的高阻异常体，联合剖面的供电极距为2个长度单位，测点间距为1个长度单位，MN长度为1个长度单位。
采用图1所示的方法流程，首先按照上述联合剖面的工作参数，在该测线上获取每个记录点的常规联合剖面的两个方向的视电阻率，即和每个参数各获取20个视电阻率数据。其中Oi为记录点的编号，Ai和Bi为对应记录点Oi的供电点编号，i为区间[1,20]之间的自然数。
按照公式计算该测线上相应记录点的联合剖面视电阻率平均值。
按照公式计算该测线上相应记录点的联合剖面视电阻率差值的绝对值。
以记录点为横坐标，以为纵坐标，绘制不同参数的曲线，如图3所示，图3中的RES-A、RES-B、RES/2、RES-ABS曲线分别为的曲线。
根据图3中的RES-A、RES-B曲线，可以初步判断该地质异常体为高阻异常体，RES/2曲线更进一步可以判断该地质异常体为高阻异常体，但根据RES-A、RES-B、RES/2曲线很难判断高阻异常体的宽度信息。综合分析图3中的RES-ABS曲线(即测量值)，可以发现该曲线上在9和12号记录点上出现极大值，故可以简单判断该高阻异常体的平面分布范围为9～12号记录点之间，即高阻异常体位于9～12长度单位范围内。
由于图3中的曲线较多，图4仅绘制RES/2、RES-ABS曲线，从图4中RES/2曲线可以很明显发现地质异常体为高阻异常体，从图4中的RES-ABS曲线可以很明显发现该高阻异常体的分布范围为9～12号记录点之间。且单纯根据RES-ABS曲线在9～12号记录点变化情况也能对地质异常体的高阻异常进行判断，从9～12号记录点之间的曲线数值明显大于0，从而可以对异常体给出高阻异常的判断。该判断结果与实际地质异常体的情况相符。
若联合剖面开展了多个供电极距的测量，则对其他供电极距的数据也采用上述的方式进行处理，从而进一步提高勘探效果。
实施例2
假设要对20个长度单位的测线开展联合剖面，其在8～13个长度单位位置下面埋藏一个长度为5个长度单位、顶面埋深1个长度单位、厚度为3个长度单位的低阻异常体，联合剖面的供电极距为2个长度单位，测点间距为1个长度单位，MN长度为1个长度单位。
采用图1所示的方法流程，首先按照上述联合剖面的工作参数，在该测线上获取每个记录点的常规联合剖面的两个方向的视电阻率，即和每个参数各获取20个视电阻率数据。其中Oi为记录点的编号，Ai和Bi为对应记录点Oi的供电点编号，i为区间[1,20]之间的自然数。
按照公式计算该测线上相应记录点的联合剖面视电阻率平均值。
按照公式计算该测线上相应记录点的联合剖面视电阻率差值的绝对值。
以记录点为横坐标，以为纵坐标，绘制不同参数的曲线，如图5所示，图5中的RES-A、RES-B、RES/2、RES-ABS曲线分别为的曲线。
根据图5中的RES-A、RES-B曲线，可以初步判断该地质异常体为低阻异常体，RES/2曲线更进一步可以判断该地质异常体为低阻异常体，但根据RES-A、RES-B、RES/2曲线很难判断低阻异常体的宽度信息。综合分析图5中的RES-ABS曲线(即测量值)，可以发现该曲线上在8和13号记录点上出现极大值，故可以简单判断该低阻异常体的平面分布范围为8～13号记录点之间，即低阻异常体位于8～13长度单位范围内。
由于图5中的曲线较多，图6仅绘制RES/2、RES-ABS曲线，从图6中RES/2曲线可以很明显发现地质异常体为低阻异常体，从图6中的RES-ABS曲线可以很明显发现该低阻异常体的分布范围为8～13号记录点之间。且单纯根据RES-ABS曲线在8～13号记录点变化情况也能对地质异常体的低阻异常进行判断，从8～13号记录点之间的曲线数值明显接近于0，从而可以对异常体给出低阻异常的判断。该判断结果与实际地质异常体的情况相符。
若联合剖面开展了多个供电极距的测量，则对其他供电极距的数据也采用上述的方式进行处理，从而进一步提高勘探效果。
以上所述仅为本发明的较佳实施例，并不用以限制本发明，凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明内。

Claims (1)

1.一种基于联合剖面的二维电法勘探方法，其特征在于：该方法基于电法勘探中的联合剖面方法获取的记录点的两个不同方向的相同供电极距的电法数据进行后期数据处理，并对其开展综合解译；该方法的步骤为：
a)按照常规的联合剖面方法获取某条测线的多个记录点的两个不同方向的联合剖面电法数据；
b)对相同记录点的两个不同方向的相同供电极距的视电阻率数据，按照公式计算两个方向的视电阻率数据的平均值，按照公式计算两个方向的视电阻率数据的差值的绝对值；其中Oi为记录点的编号，Ai和Bi为对应记录点Oi的供电点编号，i为自然数；代表Ai供电点供电时，Oi记录点的联合剖面视电阻率；代表Bi供电点供电时，Oi记录点的联合剖面视电阻率；代表Oi记录点的Ai和Bi供电点的联合剖面两个方向的视电阻率平均值；代表Oi记录点的Ai和Bi供电点的联合剖面两个方向的视电阻率差值的绝对值；ABS代表对取绝对值；其中Ai供电点到Oi记录点的距离等于Bi供电点到Oi记录点的距离；
c)对相同记录点的两个不同方向的相同供电极距的视极化率数据，按照公式计算两个方向的视极化率数据的平均值，按照公式计算两个方向的视极化率数据的差值的绝对值；其中Oi为记录点的编号，Ai和Bi为对应记录点Oi的供电点编号，i为自然数；代表Ai供电点供电时，Oi记录点的联合剖面视极化率；代表Bi供电点供电时，Oi记录点的联合剖面视极化率；代表Oi记录点的Ai和Bi供电点的联合剖面两个方向的视极化率平均值；代表Oi记录点的Ai和Bi供电点的联合剖面两个方向的视极化率差值的绝对值；ABS代表对取绝对值；其中Ai供电点到Oi记录点的距离等于Bi供电点到Oi记录点的距离；
d)按照相同电法参数类型、相同供电极距、相同测线的要求，绘制曲线；把同一条测线所有记录点的相同供电极距的所有电法数据分别绘制曲线，每条曲线所有记录点的的电法参数类型、供电极距、测线相同，纵坐标为电法数据值，横坐标为电法数据的记录点；
e)对同一条测线所有记录点的不同电法曲线进行综合解译，利用曲线主要判断电法异常体的电法特性，利用或曲线的极大值点主要判断地质异常体的分布范围，利用或曲线的异常段数值判断地质异常体的电阻率或极化率特征，从而解决地质异常体的埋藏位置、分布范围等问题。
CN201910386820.1A 2019-05-10 2019-05-10 一种基于联合剖面的二维电法勘探方法 Pending CN110208868A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910386820.1A CN110208868A (zh) 2019-05-10 2019-05-10 一种基于联合剖面的二维电法勘探方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910386820.1A CN110208868A (zh) 2019-05-10 2019-05-10 一种基于联合剖面的二维电法勘探方法

Publications (1)

Publication Number Publication Date
CN110208868A true CN110208868A (zh) 2019-09-06

Family

ID=67787060

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910386820.1A Pending CN110208868A (zh) 2019-05-10 2019-05-10 一种基于联合剖面的二维电法勘探方法

Country Status (1)

Country Link
CN (1) CN110208868A (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111208570A (zh) * 2020-01-20 2020-05-29 中南大学 一种单个电阻率异常体的频率域传导电法的工作参数优化方法
CN111221043A (zh) * 2020-01-20 2020-06-02 中南大学 一种被动源电场法的工作参数优化方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104267437A (zh) * 2014-08-29 2015-01-07 桂林理工大学 利用联合剖面视电阻率曲线定位人工湿地堵塞区域的方法
JP2015028458A (ja) * 2013-07-30 2015-02-12 韓国地質資源研究院Ｋｏｒｅａ Ｉｎｓｔｉｔｕｔｅ Ｏｆ Ｇｅｏｓｃｉｅｎｃｅ Ａｎｄ Ｍｉｎｅｒａｌ Ｒｅｓｏｕｒｃｅｓ（Ｋｉｇａｍ） 電気比抵抗探査資料獲得方法
CN109356585A (zh) * 2018-11-16 2019-02-19 青海省第三地质矿产勘查院 一种判断金属硫化物矿体产状的方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2015028458A (ja) * 2013-07-30 2015-02-12 韓国地質資源研究院Ｋｏｒｅａ Ｉｎｓｔｉｔｕｔｅ Ｏｆ Ｇｅｏｓｃｉｅｎｃｅ Ａｎｄ Ｍｉｎｅｒａｌ Ｒｅｓｏｕｒｃｅｓ（Ｋｉｇａｍ） 電気比抵抗探査資料獲得方法
CN104267437A (zh) * 2014-08-29 2015-01-07 桂林理工大学 利用联合剖面视电阻率曲线定位人工湿地堵塞区域的方法
CN109356585A (zh) * 2018-11-16 2019-02-19 青海省第三地质矿产勘查院 一种判断金属硫化物矿体产状的方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
L.MARESCOT ET AL.: "Assessing reliability of 2D resistivity imaging in mountain permafrost studies using the depth of investigation index method", 《NEAR SURFACE GEOPHYSICS》 *
傅良魁: "聚焦法与二极法及三极联合剖面法的联系和比较", 《桂林冶金地质学院学报》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111208570A (zh) * 2020-01-20 2020-05-29 中南大学 一种单个电阻率异常体的频率域传导电法的工作参数优化方法
CN111221043A (zh) * 2020-01-20 2020-06-02 中南大学 一种被动源电场法的工作参数优化方法

Similar Documents

Publication Publication Date Title
CN110208867A (zh) 一种基于联合剖面的三维电法勘探方法
Karous et al. Linear filtering of VLF dip‐angle measurements
CN108873083B (zh) 一种人工场源频率域电磁视电阻率测量方法
Edwards The magnetometric resistivity method and its application to the mapping of a fault
Pérez-Flores et al. Imaging low-frequency and dc electromagnetic fields using a simple linear approximation
CN110208868A (zh) 一种基于联合剖面的二维电法勘探方法
CN102798898A (zh) 大地电磁场非线性共轭梯度三维反演方法
CN104265279B (zh) 断层条件下随钻测井曲线预测方法
CN109100808A (zh) 一种多线源瞬变电磁横磁极化场探测方法
CN110007357A (zh) 一种航空tem和航空mt联合反演方法
CN108873073B (zh) 一种基于网络并行电法的三维跨孔电阻率层析成像方法
CN106019386A (zh) 一种回线源内比值视电阻率测定方法
CN109613066A (zh) 一种利用直流充电法测定地下水流速和优势流向的方法
Li et al. Aspects of charge accumulation in dc resistivity experiments 1
CN104318003A (zh) 基于tdcm模型的变电站地表电位计算及选址检测方法
CN112630842B (zh) 一种应用于直流电阻率法勘探的偶极差分装置及方法
CN108535783B (zh) 一种对电阻率断面进行层位校正的方法及装置
CN105700026B (zh) 一种地下目标体导电指数提取方法与装置
Nair et al. Experimental studies on the electromagnetic response of tilted conducting half-planes to a horizontal-loop prospecting system
CN109001824B (zh) 一种校正二极排列非无穷远点影响的方法
CN105137495A (zh) 一种油气检测的方法及系统
CN110109184B (zh) 一种基于多日变点的被动场源类三维电场勘探方法
CN103485772B (zh) 测井设备、方法及装置
CN111812729A (zh) 用于优化电法勘探反演的多维地质背景模型构建方法
CN107203007B (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
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20190906

WD01 Invention patent application deemed withdrawn after publication