CN111175831B - 一种估算地下水水平流动特性的传导电法勘探方法 - Google Patents

一种估算地下水水平流动特性的传导电法勘探方法 Download PDF

Info

Publication number
CN111175831B
CN111175831B CN202010108096.9A CN202010108096A CN111175831B CN 111175831 B CN111175831 B CN 111175831B CN 202010108096 A CN202010108096 A CN 202010108096A CN 111175831 B CN111175831 B CN 111175831B
Authority
CN
China
Prior art keywords
power supply
horizontal flow
underground water
measuring electrodes
measuring
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 - Fee Related
Application number
CN202010108096.9A
Other languages
English (en)
Other versions
CN111175831A (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 CN202010108096.9A priority Critical patent/CN111175831B/zh
Publication of CN111175831A publication Critical patent/CN111175831A/zh
Application granted granted Critical
Publication of CN111175831B publication Critical patent/CN111175831B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

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/02Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation operating with propagation of electric current
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01DMEASURING NOT SPECIALLY ADAPTED FOR A SPECIFIC VARIABLE; ARRANGEMENTS FOR MEASURING TWO OR MORE VARIABLES NOT COVERED IN A SINGLE OTHER SUBCLASS; TARIFF METERING APPARATUS; MEASURING OR TESTING NOT OTHERWISE PROVIDED FOR
    • G01D21/00Measuring or testing not otherwise provided for
    • G01D21/02Measuring two or more variables by means not covered by a single other subclass
    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02ATECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/30Assessment of water resources

Landscapes

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

Abstract

一种估算地下水水平流动特性的传导电法勘探方法,该方法通过在与地下水相连的部位添加易溶的矿物质(如食盐),并通过至少3个测量电极开展传导电法勘探,探测勘探区在添加易溶物质之前和之后电阻率变化情况,采用相关公式对地下水的水平流向、流速幅值、流量进行估算,从而获得地下水的水平流动特性。若获取多个供电电极的勘探结果,则能获取地下水不同深度内的水平流动特性,即获得水平流动特性的三维特征。该方法具有快速、无损、成本低等优势,适用于对上游有出露部位的中浅层地下水进行水平流向、流速、流量等特性的估算。

Description

一种估算地下水水平流动特性的传导电法勘探方法
技术领域
本发明涉及一种勘查地球物理领域的估算地下水水平流动特性的传导电法,该方法能实现对上游有出露部位的地下水进行水平流向、流速、流量等流动特性的估算。
背景技术
目前区域地下水流场主要是考虑区域地势整体高差,对照相关地质环境图集来判别地下水流向,该方法判断的地下水流向精度较低,一般只作为区域上的地下水流场的一种粗略、概略性的判别方法。
局部的地下水的流向等流动特性主要是通过打观测井,加上一些辅助试验来判别。
如常规化学离子示踪法,该方法基本原理是通过观测井将某种化学物质投入到地下水后,溶解成化学离子。若被检测到的试验中有投放的化学离子,即可证明投放点和接收点的地下水是连通的,再根据各个井监测到时间的差,概略估计地下水流向。该方法需要多个观测点或观测井,观测的辅助设施的成本过高。
如等水位线法三点法,该方法确定地下水流向和流速的步骤为:1)选择不在一条直线上的三个井点分别测定井的高程,分别测出地下水位埋深(井口地面到水面的距离);2)井口地面高程减去地下水位埋深等于地下水面高程,将三井点准确定位在地形底图上,三点连线做一个三角形,在每条边上内插高程数值,相等高程值连成曲线即是等水位线;3)垂直等水位线方向由高水位到低水位即是地下水流向;4)地下水流向确定后,进行抽水试验,确定含水层参数;5)最后,按照达西定律求出地下水流速值。这种方法需要打大口径抽水井,对于含水层埋深大时尤为不利。由于含水层的非均质性等原因,致使相关参数都可能存在较大误差,故求出的地下水流速可靠性差。用这种方法求出的地下水流向、水力坡度、渗透系数和流速值都是调查区一定范围内该参数在各点的加权平均值,各参数代表范围较大。该方法的实施需要多个观测井,且观测井越多,精度越高;故辅助设施的成本也过高,但该方法的估算结果精度较高。
如环境同位素年龄法。它是应用氢或碳等环境同位素年龄测定配合地质、水文地质等因素的分析确定地下水的流向流速。如果在区域性地下水流系统中知道地下水流的总方向,那么沿水流方向所取的各个样品中同位素年龄的差异便可用来计算区域平均地下水流速·用同位素鉴别并圈定出地下水的年龄古老而流速低的地下水系统的区段。该方法精度较高,但同位素年龄测定效率低,难以在现场直接获得地下水的流动特性,另该方法也需要有多个取样点,故辅助设施的费用高。
如自然电场法。该方法是通过测量某个中心不同方位上的自然电位,从而根据自然电位的8字形特征分析地下水的流向。该方法具有成本低、效率高等优势,但无法对地下水的流速和流量等特性进行估算。且自然电场的测量精度较低,抗干扰能力弱。
如充电法。该方法是以揭露地下水的钻孔为中心布置测量电极,通过钻孔投放食盐,并通过钻孔供电,测量等电位点,从而对地下水的流向和流速进行判断。该方法具有精度高的特点,但是该方法在实际勘探中由于需要寻找到等电位点,从而需要进行多个测点的多次勘探以便寻找到真正的等电位点,从而其勘探效率较低,另也导致该方法的环境适应能力较弱,很难满足地形陡峭、植被发育、建筑物多的勘探区的地下水勘探工作。另该方法无法对地下水的流量进行估计。
频域电磁感应电法开展地下水的流动特性估算具有可以通过改变供电信号的频率而实现不同的勘探深度内的估算工作,但电磁感应方法相对于传导电法而言,信号较弱,从而影响其估算精度;故利用传导电法方式估算地下水的流动特性将具有其独特的优势。
发明内容:
本发明是基于现有地下水的流动特性的测量方法所具有的优点和缺点的基础上,特别是基于充电法的优缺点基础上,提出的一种新的估算地下水流动特性的传导型电法勘探方法。该方法能对地下水的水平流向、水平流速、水平流量等水平流动特性进行估算,从而为地下水的水平流动特性测定提供了一种新的选择。该方法相对于充电法而言,具有如下的优点:1)工作效率更高,勘探成本低;该方法无需寻找等电位点,仅布置几个固定测量电极,就能获得地下水的水平流动特性;2)可以获得地下水的水平流动三维特性;该方法通过布置多个供电电极,从而可以获得不同深度范围内的地下水的流动特性,而不仅仅是对地下水的平面流动特性进行估算;3)可以获得地下水的时间特性;该方法通过测量不同时间的地下水流动特性,从而可以测量出地下水随着时间变化的流动特性,即可以获得地下水的四维流动特性;4)可以获得地下水的流量特性;该方法由于供电电极不布置在地下水的出露部位,从而可以获得深度信息,从而可以对地下水的流量特性进行估算。
一种估算地下水水平流动特性的传导电法勘探方法,其具体步骤如下:
a)选择需要估算地下水水平流动特性的勘探区,布置至少3个测量电极(且任意3个测量电极不在一条直线上)和1个供电电极以及1个电法勘探意义上的无穷远供电电极,并获取测量电极和供电电极的水平坐标,测量电极的水平坐标为(CXi,CYi),供电电极的水平坐标为(GX,GY);其中i∈[1,n],n为测量电极的总数量;优选地,所有测量电极以投放易溶矿物质的地下水的出露点的中心部位为中心,等角度、与投放易溶矿物质的地下水的出露点的中心部位相等水平距离的均匀分布。供电电极到投放易溶矿物质的地下水出露点的中心部位S的水平距离大于任一测量电极到投放易溶矿物质的地下水出露点的中心部位S的水平距离;且供电电极不在任意2个测量电极连线的中垂线上。若要获得地下水的水平流动特性的三维特征,则布置多个供电电极,每个供电电极供电时,均按照后续的单个供电电极的所有步骤开展,从而能获得不同深度内的地下水的水平流动特性的三维特征。若要获得地下水的更高精度的流动特性,则可以增加测量电极数量,即测量电极数量越多,则勘探精度越高,但需要注意任意3个测量电极不在一条直线上,以免造成某些测量电极对上的信号过弱,导致测量精度降低的现象。固体的易溶矿物质一般选用食盐,一方面可以减少污染,另一方面减少勘探成本,第三方面食盐容易获取。
b)通过供电电极配合无穷远供电电极供电,测量相邻2个测量电极上的电信号,即采用电法勘探的三极测量装置;供电电极供电时,所有相邻2个测量电极上的电信号同时测量,并换算相应的视电阻率参数,从而获得一系列的视电阻率参数ρi,j,k,t0;其中i和j为2个测量电极的编号,且i和j不相等;k为供电电极的编号;t0为k号供电电极供电时,i和j号测量电极上第一次测量电信号的时刻;i∈[1,n],j∈[1,n],n为测量电极的总数量;所有相邻的测量电极之间的电信号同时测量的主要目的是减少随机干扰,若测量电极数量过多,而电法仪的测量通道数不足以同时测量,则尽量减少滞后时间,以便降低随机干扰的影响程度。若需要获取三维特征,则该步骤中需要对每个供电电极供电时的视电阻率进行测量。
c)第一次测量电信号完毕后,立即通过地下水的出露点,向地下水投放无害的、固体的易溶矿物质(如食盐),并保持在测量过程中地下水中的固体的易溶矿物质总存在固体状态的易溶矿物质;记录投放无害的、固体易溶矿物质的地下水出露点中心部位的坐标(Xs,Ys);S表示该地下水出露点中心部位;之所以要保持在测量过程中的固体易溶矿物质有固体状态存在,是尽量保证地下水中的该种易溶矿物质处于饱和状态,从而尽量保证流动特性的估算精度。
d)投放易溶矿物质后,每隔一段时间,通过供电电极和无穷远供电电极供电,同时测量任意相邻的2个测量电极上的电信号,并换算相应的视电阻率参数,从而获得一系列视电阻率参数ρi,j,k,t;其中i和j为2个测量电极的编号,且i和j不相等;k为供电电极的编号,i∈[1,n],j∈[1,n],n为测量电极的总数量;t为投放易溶矿物质后k号供电电极供电时,i和j号测量电极上测量电信号的时刻;该步骤中的间隔时间需要基于选用的电法仪等设备的情况而定,若电法仪适合快速测量,则可以减少间隔时间,并适当考虑地下水的大体情况,如地下水流速快,则可以减少间隔时间;若流速慢,则可以加大间隔时间。若需要获取三维特征,则该步骤中需要对每个供电电极供电时的视电阻率进行测量。
e)通过公式(1)选择t0和tD时刻的视电阻率数值,其中tD>t0,tD为投放易溶矿物质后,供电电极供电时,测量电极上测量信号的某个时刻,且该时刻的相同测量电极的视电阻率ρi,j,k,tD与t0时刻的视电阻率ρi,j,k,t0不全部相等;根据测量电极和投放无害的易溶矿物质的地下水出露点中心部位S的水平坐标,先求取每个向量
Figure BDA0002389037100000061
再求取所有向量的和,从而估算地下水的水平流向,即地下水的水平流向为向量
Figure BDA0002389037100000071
的方位;
Figure BDA0002389037100000072
其中i和j不相等,i∈[1,n],j∈[1,n],k为供电电极的编号,n为测量电极的总数量;ρi,j,k,t0表示k号供电电极供电时t0时刻在i和j测量电极上测量的视电阻率;ρi,j,k,tD为k号供电电极供电时tD时刻在i和j测量电极上测量的视电阻率;tD的时刻选择比较重要和关键,若投入的易溶矿物质还未通过测量电极时,则此时所有测量电极上所获取的视电阻率相对于t0时刻的视电阻率来说,还未发生变化,即ρi,j,k,tD=ρi,j,k,t0,则公式(1)的计算结果为0,但并不是地下水不流动,而是测量电极还未测量到地下水的变化状态,故tD的时刻需要满足“该时刻的相同测量电极的视电阻率ρi,j,k,tD与t0时刻的视电阻率ρi,j,k,t0不全部相等”的条件,从而保证该时刻已经测量到地下水的流动特性。若需要获取三维特征,则该步骤中需要对每个供电电极供电时的结果进行单独计算,即获取不同供电电极供电时的水平流向结果。也可以单独求取多个不同的tD时刻的水平流向,再求取水平流向平均值,从而获得精度更高的估算结果。
f)通过公式(2)估算地下水水平流速的幅值V;
Figure BDA0002389037100000081
其中i∈[1,n],j∈[1,n],i和j不相等,n为测量电极的总数量;
其中S为投放易溶矿物质的地下水出露点中心部位,Oi,j表示第i号和第j号测量电极连线的中点,SOi,j表示投放易溶矿物质的地下水的出露点中心部位S到第i号和第j号测量电极的中点Oi,j的水平距离,单位为米;tmin,i,j为第i号和第j号测量电极上的视电阻率第一次达到极小值的时刻;t0为未投放易溶矿物质时,第一次测量第i号和第j号测量电极上的视电阻率的时刻,由于第一次测量是同一时刻测量,故公式(2)中该参数为相同值;S-Oi,j表示S与Oi,j之间的连线,θi,j为地下水的水平流向(即公式(1)中向量
Figure BDA0002389037100000082
的方位)与S-Oi,j的锐夹角;max()为求取最大值的函数;把公式(2)中最大值函数中的最大值结果所对应的i和j测量电极编号分别命名为p和q,若公式(2)中最大值函数中的变量存在多个最大值,则任选一个最大值所对应的i和j测量电极编号分别命名为p和q;若需要获取三维特征,则该步骤中需要对每个供电电极供电时的结果进行单独计算,即获取不同供电电极供电时的水平流量幅值的结果。
g)通过公式(3)估算地下水的水平流量W;
Figure BDA0002389037100000091
其中ρp,q,k,t0为k号供电电极供电时t0时刻在p和q测量电极上测量的视电阻率;
Figure BDA0002389037100000093
为k号供电电极供电时tmin,p,q时刻在p和q测量电极上测量的视电阻率;tmin,p,q为第p和第q号测量电极上测量的视电阻率第一次达到极小值时的时刻,p和q为步骤f)所确定的测量电极编号;V为公式(2)的水平流速的幅值计算结果,即地下水的水平流速幅值,单位为米/小时;L为k号供电电极到第p号和第q号测量电极连线的中点Op,q的水平距离,单位为米,即
Figure BDA0002389037100000094
Dp,q为p号和q号测量电极之间的水平距离,单位为米,即
Figure BDA0002389037100000095
S-Op,q表示投放易溶矿物质的地下水的出露点中心部位S与Op,q之间的连线,θp,q为地下水的水平流向(即公式(1)中向量
Figure BDA0002389037100000096
的方位)与S-Op,q的锐夹角。若需要获取三维特征,则该步骤中需要对每个供电电极供电时的结果进行单独计算,即获取不同供电电极供电时(也代表不同深度内)的水平流量结果。
通过以上的数据获取及计算,从而获得地下水的L/3深度范围内水平方向的水平流向、水平流速的幅值、水平流量等流动特性的估算。
如果需要对地下水不同深度内的水平方向的流动特性进行估算,则需要增加供电电极数量,改变供电电极的位置,并按照以上a)~g)步骤进行不同深度内(或说不同供电电极供电时)的地下水流动特性估算。把不同深度内的地下水的流向、流速、流量等流动特性进行综合,从而可以获得地下水的三维的流动特性,可以更加清晰表达地下水的流动特性的空间特征。如把不同时间段的地下水的三维流动特性进行获取,即获取地下水的流动特性的时空特征,则能对地下水的流动特性在时间和三维空间的波动情况进行更好地刻画。
附图说明:
图1为本发明的一种估算地下水水平流动特性的传导电法勘探方法流程图;
图2为本发明的一种估算地下水水平流动特性的传导电法勘探方法的野外布置示意图;
其中图2中的十字形符号及编号A,表示供电电极位置及编号;矩形空心符号及编号∞,表示无穷远供电电极位置及编号;黑色圆形实心点及编号(1/2/3)为测点及编号;黑色圆形空心点及编号(O1,2、O2,3、O1,3)为相邻测量电极的连线的中点及编号;黑色三角形实心点及编号(S)为地下水出露点的中心部位,并作为投放固体的易溶矿物质的部位;
具体实施方式:
以下参照图1、图2,结合具体实施方式对本发明做进一步说明。
按照图1所示的一种估算地下水水平流动特性的传导电法勘探方法流程图中的步骤实施本方法,其具体步骤如下:
a)按照图2所示,选择一个需要估算地下水水平流动特性的勘探区,在勘探区内布置3个测量电极(图2中的1/2/3测量电极)(且3个测量电极不在一条直线上)和1个供电电极(图2中的A供电电极)以及1个电法勘探意义上的无穷远供电电极(如图2中的∞为无穷远供电电极),并获取测量电极和供电电极的坐标,测量电极的坐标为(CXi,CYi)(如图2中的1号测量电极坐标为(CX1=10,CY1=50),2号测量电极坐标为(CX2=70,CY2=10),3号测量电极坐标为(CX3=90,CY3=70)),供电电极的坐标为(GXA,GYA)(如图2中的A供电电极的坐标为(GXA=60,GYA=90));供电电极A到投放易溶矿物质的地下水出露点的中心部位(图2中的S)的距离大于任一测量电极(图2中的1/2/3测量电极)到投放易溶矿物质的地下水出露点的中心部位S的距离;且供电电极(图2中的A供电电极)不在任意2个测量电极(图2中的1/2/3测量电极)连线的中垂线上。
b)通过供电电极(如图2中的A)配合无穷远供电电极(如图2中的∞)供电,测量相邻2个测量电极(如图2中的1和2、2和3、3和1测量电极)上的电信号,即采用电法勘探的三极测量装置;供电电极A供电时,所有相邻2个测量电极(如图2中的1和2、2和3、3和1测量电极)上的电信号同时测量,并换算相应的电阻率参数,从而获得一系列的电阻率参数ρi,j,A,t0;其中i和j为2个测量电极的编号,且i和j不相等;A为供电电极的编号;t0为A号供电电极供电时,i和j号测量电极上第一次测量电信号的时刻;i∈[1,3],j∈[1,3],即有3个测量电极;
c)第一次测量电信号完毕后,立即通过地下水的出露点,向地下水投放食盐,并保持地下水中的食盐处于饱和状态;记录投放食盐的地下水出露点中心部位S的坐标(Xs,Ys)(如图2中的投放点的坐标为(Xs=60,Ys=50));
d)投放食盐后,每隔一段时间(如每隔1个小时测量一次),通过供电电极A和无穷远供电电极∞供电,同时测量任意相邻的2个测量电极(如图2中的1和2、2和3、3和1测量电极)上的电信号,并换算相应的电阻率参数,从而获得一系列电阻率参数ρi,j,A,t;其中i和j为2个测量电极的编号,且i和j不相等;A为供电电极的编号,i∈[1,3],j∈[1,3];t为投放易溶矿物质后A供电电极供电时,i和j号测量电极上测量电信号的时刻;
e)通过公式(1)选择t0和tD时刻的电阻率数值,其中tD>t0,根据测量电极(如图2中的1和2、2和3、3和1测量电极)和投放食盐的地下水出露点中心部位S的水平坐标,先求取每个向量
Figure BDA0002389037100000131
再求取所有向量的和,从而估算地下水的水平流向,即地下水的水平流向为向量
Figure BDA0002389037100000132
的方位;
Figure BDA0002389037100000133
其中i和j不相等,i∈[1,3],j∈[1,3],k∈[1];ρi,j,A,t0表示A供电电极供电时t0时刻在i和j测量电极上测量的电阻率;ρi,j,A,tD为A供电电极供电时tD时刻在i和j测量电极上测量的电阻率;时刻tD在公式(1)中为一固定值,即公式(1)中的tD维持不变;
f)通过公式(2)估算地下水的水平流速幅值V;
Figure BDA0002389037100000135
其中i∈[1,3],j∈[1,3],i和j不相等,3为测量电极的总数量;其中S为投放食盐的地下水出露点中心部位,Oi,j表示第i号和第j号测量电极连线的中点,SOi,j表示投放易溶矿物质的地下水的出露点中心部位S到第i号和第j号测量电极的中点Oi,j的距离,单位为米;tmin,i,j为第i号和第j号测量电极上的电阻率达到极小值的时刻(根据不同时刻相同供电电极供电,相同测量电极对所获取的视电阻率,从而获得极小值的时刻);t0为未投放食盐时,第一次测量第i号和第j号测量电极上的电阻率的时刻,由于第一次测量是同一时刻测量,故公式(2)中该参数为相同值;S-Oi,j表示S与Oi,j之间的连线,θi,j为地下水的水平流向(即公式(1)中向量
Figure BDA0002389037100000141
的方位)与S-Oi,j的锐夹角;max()为求取最大值的函数;把公式(2)中最大值函数中的最大值结果所对应的i和j测量电极编号分别命名为p和q,若公式(2)中最大值函数中的变量存在多个最大值,则任选一个最大值所对应的i和j测量电极编号分别命名为p和q;
g)通过公式(3)估算地下水的水平流量W;
Figure BDA0002389037100000142
其中ρp,q,A,t0为A供电电极供电时t0时刻在p和q测量电极上测量的电阻率;ρp,q,A,tmin,p,q为A供电电极供电时tmin,p,q时刻在p和q测量电极上测量的电阻率;tmin,p,q为第p和第q号测量电极上测量的电阻率达到极小值时的时刻,p和q为步骤f)所确定的测量电极编号;V为公式(2)的水平流速幅值计算结果,即地下水的水平流速幅值,单位为米/小时;L为A供电电极到第p号和第q号测量电极连线的中点Op,q的距离,单位为米,即
Figure BDA0002389037100000151
Dp,q为p号和q号测量电极之间的距离,单位为米,即
Figure BDA0002389037100000152
S-Op,q表示投放易溶矿物质的地下水的出露点中心部位S与Op,q之间的连线,θp,q为地下水的水平流向与S-Op,q的锐夹角。
以上所述仅为本发明的较佳实施例,并不用以限制本发明,凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明内。

Claims (7)

1.一种估算地下水水平流动特性的传导电法勘探方法,其具体步骤如下:
a)选择需要估算地下水水平流动特性的勘探区,布置至少3个测量电极和1个供电电极以及1个电法勘探意义上的无穷远供电电极,且任意3个测量电极不在一条直线上,并获取测量电极和供电电极的水平坐标,测量电极的水平坐标为(CXi,CYi),供电电极的水平坐标为(GX,GY);其中i∈[1,n],n为测量电极的总数量;
b)通过供电电极配合无穷远供电电极供电,测量相邻2个测量电极上的电信号,即采用电法勘探的三极测量装置;供电电极供电时,所有相邻2个测量电极上的电信号同时测量,并换算相应的视电阻率参数,从而获得一系列的视电阻率参数
Figure FDA0002955536010000011
其中i和j为2个测量电极的编号,且i和j不相等;k为供电电极的编号;t0为k号供电电极供电时,i和j号测量电极上第一次测量电信号的时刻;i∈[1,n],j∈[1,n],n为测量电极的总数量;
c)第一次测量电信号完毕后,立即通过地下水的出露点,向地下水投放无害的、固体的易溶矿物质,并保持在测量过程中地下水中的固体的易溶矿物质总存在固体状态的易溶矿物质;记录投放无害的、固体易溶矿物质的地下水出露点中心部位的坐标(Xs,Ys);S表示该地下水出露点中心部位;
d)投放易溶矿物质后,每隔一段时间,通过供电电极和无穷远供电电极供电,同时测量任意相邻的2个测量电极上的电信号,并换算相应的视电阻率参数,从而获得一系列视电阻率参数ρi,j,k,t;其中i和j为2个测量电极的编号,且i和j不相等;
k为供电电极的编号,i∈[1,n],j∈[1,n],n为测量电极的总数量;t为投放易溶矿物质后k号供电电极供电时,i和j号测量电极上测量电信号的时刻;
e)通过公式(1)选择t0和tD时刻的视电阻率数值,其中tD>t0,tD为投放易溶矿物质后,供电电极供电时,测量电极上测量信号的某个时刻,且该时刻的相同测量电极的视电阻率ρi,j,k,tD与t0时刻的视电阻率
Figure FDA0002955536010000021
不全部相等;根据测量电极和投放无害的易溶矿物质的地下水出露点中心部位S的水平坐标,先求取每个向量
Figure FDA0002955536010000022
再求取所有向量的和,从而估算地下水的水平流向,即地下水的水平流向为向量
Figure FDA0002955536010000023
的方位;
Figure FDA0002955536010000024
其中i和j不相等,i∈[1,n],j∈[1,n],k为供电电极的编号,n为测量电极的总数量;
Figure FDA0002955536010000025
表示k号供电电极供电时t0时刻在i和j测量电极上测量的视电阻率;ρi,j,k,tD为k号供电电极供电时tD时刻在i和j测量电极上测量的视电阻率;
f)通过公式(2)估算地下水水平流速的幅值V;
Figure FDA0002955536010000031
其中i∈[1,n],j∈[1,n],i和j不相等,n为测量电极的总数量;其中S为投放易溶矿物质的地下水出露点中心部位,Oi,j表示第i号和第j号测量电极连线的中点,SOi,j表示投放易溶矿物质的地下水的出露点中心部位S到第i号和第j号测量电极的中点Oi,j的水平距离,单位为米;tmin,i,j为第i号和第j号测量电极上的视电阻率第一次达到极小值的时刻;t0为未投放易溶矿物质时,第一次测量第i号和第j号测量电极上的视电阻率的时刻,由于第一次测量是同一时刻测量,故公式(2)中该参数为相同值;S-Oi,j表示S与Oi,j之间的连线,θi,j为地下水的水平流向与S-Oi,j的锐夹角;max()为求取最大值的函数;把公式(2)中最大值函数中的最大值结果所对应的i和j测量电极编号分别命名为p和q,若公式(2)中最大值函数中的变量存在多个最大值,则任选一个最大值所对应的i和j测量电极编号分别命名为p和q;
g)通过公式(3)估算地下水的水平流量W;
Figure FDA0002955536010000041
其中
Figure FDA0002955536010000042
为k号供电电极供电时t0时刻在p和q测量电极上测量的视电阻率;
Figure FDA0002955536010000043
为k号供电电极供电时tmin,p,q时刻在p和q测量电极上测量的视电阻率;tmin,p,q为第p和第q号测量电极上测量的视电阻率第一次达到极小值时的时刻,p和q为步骤f)所确定的测量电极编号;V为公式(2)的水平流速的幅值计算结果,即地下水的水平流速幅值,单位为米/小时;L为k号供电电极到第p号和第q号测量电极连线的中点Op,q的水平距离,单位为米,即
Figure FDA0002955536010000044
Dp,q为p号和q号测量电极之间的水平距离,单位为米,即
Figure FDA0002955536010000045
Figure FDA0002955536010000046
S-Op,q表示投放易溶矿物质的地下水的出露点中心部位S与Op,q之间的连线,θp,q为地下水的水平流向与S-Op,q的锐夹角。
2.如权利要求1所述的一种估算地下水水平流动特性的传导电法勘探方法,其特征在于:优选地,所有测量电极以投放易溶矿物质的地下水的出露点的中心部位S为中心,等角度、与投放易溶矿物质的地下水的出露点的中心部位S相等的水平距离的均匀分布。
3.如权利要求1所述的一种估算地下水水平流动特性的传导电法勘探方法,其特征在于:供电电极到投放易溶矿物质的地下水出露点的中心部位S的水平距离大于任一测量电极到投放易溶矿物质的地下水出露点的中心部位S的水平距离;且供电电极不在任意2个测量电极连线的中垂线上。
4.如权利要求1所述的一种估算地下水水平流动特性的传导电法勘探方法,其特征在于:若获取不同的供电电极供电情况下,相同测量电极上的电信号,根据公式(1)、(2)和(3)获取地下水相应的水平流向、水平流速幅值、水平流量的多个估算结果,则对地下水的不同深度的水平流动特性进行分析和判断。
5.如权利要求1所述的一种估算地下水水平流动特性的传导电法勘探方法,其特征在于:无害的、固体的易溶矿物质优选用食盐。
6.如权利要求1所述的一种估算地下水水平流动特性的传导电法勘探方法,其特征在于:地下水是指未投放无害的、固体的易溶矿物质时的水中所包含的易溶矿物质未达到饱和状态的水。
7.如权利要求1所述的一种估算地下水水平流动特性的传导电法勘探方法,其特征在于:利用公式(1)计算地下水的水平流向时,可以单独求取多个不同的tD时刻的水平流向,再求取平均值,获得精度更高的估算结果。
CN202010108096.9A 2020-02-21 2020-02-21 一种估算地下水水平流动特性的传导电法勘探方法 Expired - Fee Related CN111175831B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010108096.9A CN111175831B (zh) 2020-02-21 2020-02-21 一种估算地下水水平流动特性的传导电法勘探方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010108096.9A CN111175831B (zh) 2020-02-21 2020-02-21 一种估算地下水水平流动特性的传导电法勘探方法

Publications (2)

Publication Number Publication Date
CN111175831A CN111175831A (zh) 2020-05-19
CN111175831B true CN111175831B (zh) 2021-04-23

Family

ID=70658405

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010108096.9A Expired - Fee Related CN111175831B (zh) 2020-02-21 2020-02-21 一种估算地下水水平流动特性的传导电法勘探方法

Country Status (1)

Country Link
CN (1) CN111175831B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112433263B (zh) * 2020-12-01 2021-09-28 中南大学 一种快捷高效的地下水富集区直接探测方法
CN114859074B (zh) * 2022-04-25 2024-04-16 无锡水文工程地质勘察院有限责任公司 一种地下水流向测量方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101702034A (zh) * 2009-10-28 2010-05-05 湖南金宏源电子科技有限公司 基于地球电磁感应原理的地下水源探测方法及探测仪
CN102444111A (zh) * 2011-10-19 2012-05-09 邓业灿 地下溶洞造影探测法
WO2015030993A2 (en) * 2013-08-29 2015-03-05 Saudi Arabian Oil Company Mapping resistivity distribution within the earth
CN106032754A (zh) * 2016-05-18 2016-10-19 陕西煤业化工技术研究院有限责任公司 一种基于地下水流速流向测定的煤矿开采防治水方法
CN110146939A (zh) * 2019-06-24 2019-08-20 辽宁省水利水电勘测设计研究院有限责任公司(原名称为辽宁省水利水电勘测设计研究院) 一种深部地下水流速测定方法

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7386402B2 (en) * 2006-06-22 2008-06-10 Bryant Consultants, Inc. Remotely reconfigurable system for mapping structure subsurface geological anomalies
WO2017127108A1 (en) * 2016-01-22 2017-07-27 Halliburton Energy Services, Inc. Determining downhole wettability

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101702034A (zh) * 2009-10-28 2010-05-05 湖南金宏源电子科技有限公司 基于地球电磁感应原理的地下水源探测方法及探测仪
CN102444111A (zh) * 2011-10-19 2012-05-09 邓业灿 地下溶洞造影探测法
WO2015030993A2 (en) * 2013-08-29 2015-03-05 Saudi Arabian Oil Company Mapping resistivity distribution within the earth
CN106032754A (zh) * 2016-05-18 2016-10-19 陕西煤业化工技术研究院有限责任公司 一种基于地下水流速流向测定的煤矿开采防治水方法
CN110146939A (zh) * 2019-06-24 2019-08-20 辽宁省水利水电勘测设计研究院有限责任公司(原名称为辽宁省水利水电勘测设计研究院) 一种深部地下水流速测定方法

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
An improved method for shallow water flow velocity measurement with practical electrolyte inputs;Tingwu Lei et al.;《Journal of Hydrology》;20101231;第45-56页 *
An Inexpensive Multilevel Array of Sensors for Direct Ground Water Velocity Measurement;J.F. Devlin et al.;《Ground Water Monitoring & Remediation》;20090930;第29卷(第2期);第73-77页 *
Field comparison of the p oint velocity probe with other groundwater velocity measurement methods;W. Labaky et al.;《WATER RESOURCES RESEARCH》;20090314;第45卷;第1-9页 *
充电法和高密度电法在典型岩溶区勘查中的应用分析;陈松等;《工程地球物理学报》;20170131;第14卷(第1期);第6-12页 *
地下水渗流速度测量方法研究进展;戴传山等;《河北工业大学学报》;20181031;第47卷(第5期);第106-111页 *

Also Published As

Publication number Publication date
CN111175831A (zh) 2020-05-19

Similar Documents

Publication Publication Date Title
CA2921822C (en) Borehole electric field survey with improved discrimination of subsurface features
CA2183131C (en) Apparatus and method for measuring the resistivity of underground formations
US8030934B2 (en) Method for hydrocarbon reservoir mapping and apparatus for use when performing the method
US4446434A (en) Hydrocarbon prospecting method with changing of electrode spacing for the indirect detection of hydrocarbon reservoirs
CN101382599B (zh) 一种确定储层孔隙各向异性的瞬变电磁方法
CN101263404A (zh) 高分辨率的电阻率地层成像仪
EA012880B1 (ru) Способ получения данных о геофизических свойствах коллектора
US20230333276A1 (en) Method and system for exploring hidden karst pipelines
CN111175831B (zh) 一种估算地下水水平流动特性的传导电法勘探方法
Cardimona zywvutsrqponmlihgfed
CN103821495A (zh) 测井方法
CN111175832B (zh) 一种估算地下水三维流动特性的频域电磁感应勘探方法
Barrash et al. Boise Hydrogeophysical Research Site (BHRS): objectives, design, initial geostatistical results
EP0039959A2 (en) Hydrocarbon prospecting method and apparatus for the indirect detection of hydrocarbon reservoirs
CN206439047U (zh) 极板内部发射回流的电成像测量装置
Zonge et al. Effect of electrode contact resistance on electric field measurements
Bibby Electrical resistivity mapping in the central volcanic region of New Zealand
Acworth The electrical image method compared with resistivity sounding and electromagnetic profiling for investigation in areas of complex geology: A case study from groundwater investigation in a weathered crystalline rock environment
JP2790597B2 (ja) 比抵抗利用の地下流れ特性測定方法及び装置
Kirsch et al. Geoelectrical methods
CN104991281B (zh) 一种煤层埋藏深度的探测方法和装置
CN110454155B (zh) 一种应用磁化率测井方法确定第四纪地层年代的方法
Domenzain et al. 3D DC inversion, visualization, and processing of dense time-lapse data in fine domains applied to remediation monitoring
Santoso et al. Presumption of Ground Water Depth Using the Schlumberger Configuration Geoelectrical Method
Abbey et al. Determination of ground water potential using electrical resistivity 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
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20210423

CF01 Termination of patent right due to non-payment of annual fee