您的当前位置:首页正文

基于电法装置电极排列系数K值计算方法的改进

来源:个人技术集锦
2013年3O卷 贵 州 地 质 Vo1.30 No.1(To1.114)2013 ・54・ GUIZHOU GEOLOGY 第1期(总第114期) 基于电法装置电极排列系数K值计算方法的改进 彭洪军,张应文,代 迪,马 强,代璨怡 (贵州省地矿局102地质大队,贵州遵义563003) [摘要]应用常规方法校正K值时计算复杂、计算工作量大。采用电极点坐标参数(x,Y,z), 应用空间解析几何原理并借助Excel(电子表格应用软件)计算平台计算K值,实现了因电极位 移、地形起伏导致K值变化的校正。该方法准确、简单、方便、快捷、可操作性强,对山区的电法 勘探具有较强的实用性。 [关键词]K值;坐标;地形;水平距离;实际距离;电阻率 [中图分类号]P63l_3 4[文献标识码]A[文章编号]1000—5943(2013)01—0054—07 坐标参数值就可求出K值,同时实现电极位移K 1 前言 值的校正,该方法在水平地形、单一地层、综合地 层中都可运用,同时能提高数据处理效率,对山区 近年来电阻率法在山区的电法勘探中得到了 的电法勘探具有实用性。 广泛的应用,然而在山区的电法勘探中由于受客 观条件限制导致电极发生位移,还有在高程落差 2原理及工作方法 极大的山谷、山脊地段,K值变化通常超出允许误 差范围,从而使视电阻率曲线产生畸变,因此有必 要对K值进行校正。 2.1 原理 在电阻率法勘探中,比较常用的方法有二极 在地下介质均匀各向同性的条件下,电阻率 装置、三极装置、联合剖面装置、对称四极装置、偶 计算式Ps=K・a ,(1),(1)式中的K称为电 极装置、中间梯度装置等,这些方法都涉及电极排 极排列系数,K值数学表达式K=27r/(1/AⅣ一1/ 列系数K值。传统的K值计算方法一般基于水 AN一1/BM+1/BN)(2),其值是一个仅与各电极 平距离、对称剖面,计算比较简单,在地势比较单 间空间位置有关的量;在实际工作中,电极形式可 一平坦的地区十分实用。 多种多样,但电极的排列形式和电极距一经确定, 然而笔者在实际工作中发现,多数情况下的 K值便可由公式(2)计算出来;(2)式为K值的计 K值计算看似简单,其实计算起来比较复杂,比如 算通式,在计算联合剖面、三极测深等装置时B 计算中间梯度旁测剖面、非对称四极、电极发生位 极的处理有两种办法:一是将通式中的BM、BN 移等就极为复杂。在山区的电阻率法勘探中地形 作“。。”大处理,二是将c极当成B极输入电极点 起伏对K值影响较大,尤其是高程落差极大的山 坐标值;(2)式中的AM、AN、BM、BN是各电极距 谷、山脊地形比较突出,如果仍按水平距离计算K 的距离,也是本文要求解的变量,关于AM、AN、 值,通常会导致K值变化值超出允许误差。基于 BM、BN距离的求解方法请看图1。 上述原因,为达到K值计算准确、简单、方便、快 假设A、B、M、N各点坐标分别为A(A ,Ay, 捷,实现因电极位移、地形起伏导致K值变化的 Az)、B(曰 , ,Bz)、M(朋 , , )、N(以,Ⅳ , 校正,这里应用电极点坐标参数值借助Excel平 Ⅳ )。以求解AM为例,根据空间解析几何求距 台对K值求解,该方法不仅简单、方便、快捷、准 离原理,利用A、M 、AM组成的平面直角三角形 确、可操作性强,更重要的是具有较强的灵活性, 应用勾股定理求出平面斜距Lxy。 即在任意装置下只需导入A、B、M、N四个电极点 [收稿日期]2012—12—22[改回日期】2013—03—01 [作者简介]彭洪军(1987一),男,山东临沂人,物探助理工程师,主要从事地球物理勘探工作。 第1期 彭洪军,等:基于电法装置电极排列系数K值计算方法的改进 ・55・ H∞ (O 踮 O ∞ O ;O2  O O 图1任意装置K值计算原理示意图 Fig.I K value counting principle of any array Lxr= = (3) 2.2计算方法 再利用A、M、M 组成的垂面直角三角形应 用勾股定理求出空间斜距AM ,其中AM 计算公 2.2.1建立模板 式为: 本文的K值计算是借助Excel平台下的函数 AMs= 巧= (4) 计算功能,即可以在表格相应的单元格内编辑输 同理可求得AN BM BN。计算公式为: 入相应的计算公式进行相关计算,为了能方便使 用、减少重复输入、提高计算效率,因此建立一个 A^ =√(A —Nx) +(Ay—Nr) +(Az—Nz) (5) Excel模板很有必要。 BMs=√( -Mx) +(By一 ) +(Bz-Mz) (6) 我们计算的K值模板见表1(此模板在Excel BNs= ̄/( 一 ) +( 一Ⅳy) +( 一Ⅳz) (7) 中总体由A—AD 30个列和1、2、3、4的行构成, 在不考虑高程差的情况下,计算水平距离 这里由于篇幅有限,在Word文件中横向无法全 AMP、ANP、BMP、BNP公式简化如下: 部展开,故将表格一分为二转换为纵向展示,即将 A =J(Az一』I ) +(Ay— y) (8) A—O共计15列置于表格上半部分、O列以右的 P—AD共计15列置于表格下半部分)。模板内 ANt,=√(A —Nx) +(Ar—Ny) (9) 容包括:电极及记录点编号;A、M、N、B(c)电极 BMt,=√(日x—Mx) +( y— y) (1o) 点坐标;水平电极距AMP、ANP、BMP、BNP;水平距 曰Ⅳ尸=√( 一心) +( y一Ⅳy) (11) 离计算K值(以下简称水平K值);斜电极距 从上述(4)一(11)计算公式中可以看出,参 AMs、ANs、BMs、BNs;斜电极距计算K值(以下简 与计算两点距离的参数仅为两点的坐标值,这样 称实际K值);AK(即实际K值与水平K值之 就建立起了两点的空间距离与两点坐标参数值的 差,以下简称△k);K偏离度。其中第4行为计算 关系。 模板,A、M、N、B(C)点坐标值系任意充填值,只 是为了编辑公式计算需要,无实际意义,从第5行 起以下为计算区域。 ・56・ 贵州地质 2013年30卷 在建立完成模板框架结构之后,可在相应的 单元格中进行输入计算公式,计算公式详见表2。 表2 Excel计算公式 Table 2 Calculation formula of Excel 在单元格内输入Excel表达公式后,点击编 储于其它介质中的各测点坐标数据按模板格式复 辑栏中的“、/”符号,一个公式输入计算完成,所 制到对应的单元格内,数据导人完成。 有公式输入完成后将该文件保存为模板,K值计 2.2.3数据计算 算模板建立完成。 数据计算是在计算模板中利用Excel强大的 2.2.2数据导人 计算功能,仅需要简单操作即可完成计算。具体 数据导人有手工输入和文件复制两种形式, 操作过程为分别点击选择模板行(第4行)对应 手工导入效率低而且容易出错,因此常用的数据 单元格,将鼠标对准右下角黑点使其变成“+”向 导人形式为文件复制。数据导人方法是在模板中 下拖动,即可实现对模板公式的复制和数据计算。 编制完记录点及A、M、N、B(C)电极编号后,将存 第1期 彭洪军,等:基于电法装置电极排列系数K值计算方法的改进 -57・ 线为中间梯度旁测剖面,与主剖面104线平行相距 3应用实例 本应用实例取材于某铅锌矿区,矿区位于贵 州省从江县内,区内地形切割强烈、地形起伏连 绵、坡度较陡、山脊与冲沟交互出现、悬崖与陡坎 400 m,该测线地形呈“V”字型,山顶(548点)最大 高程866.2 m,山谷(566点)最小高程666.3 m,最 大高程差200 m,测点编号545—583,共计39个点; 山脊地形114线为中间梯度主剖面,该测线地形为 山脊地形,山顶(500点)最大高程1225.5m,山谷 遍布,最大相对高程差可达数百米。 (517点)最小高程913.1 m,最大高程差312.4 m, 区内作过中间梯度测量、测深测量,测量结果 测点编号482—520,共计39个点。 都存在地形起伏对K值的影响,也都存在供电电 将该实例数据导人后,应用Excel模板计算, 极位移对K值的影响。这里以中间梯度装置,山 可得到如表1格式中的水平k值、实际k值、△k 谷、山脊地形为例来展示应用电极点坐标参数值计 值和K偏离度值,对计算数据用绘图软件作出的 算K值过程及地形对K值的影响程度。中间梯度 相关参数曲线图,如图2、图3所示(图中增加了 装置参数(1)AB=4 000 m;(2)MN=120 m;(3)点 一条2%水平直线用于对K值偏离度的研判参 距20 m;(4)测网密度200×20 m。山谷地形100 考,当K值误差变化在2%以上时需作校正)。 塞 ~ \— 由翻 多 … 。一 … \。 图2 100线模板计算各参数曲线图(山谷地形) 图3 114线模板计算各参数曲线图(山脊地形) Fig.2 Curve ofparameters calculated byformwork ofline100 Fig.3 Curve ofparamete ̄calculatde byformworkofline104 当电极距呈水平状态与倾斜状态时,由模板 值,与K值变化、P值变化对比结果如下表3、表 计算的电极距不同状态装置系数K值,及其相应 4。进一步可以看出山谷、山脊及斜坡地形,对K 工作电流(ip)、工作电压(vp)下的对应电阻率P 值以及视电阻率变化的影响。 表3 100线电极距呈水平状态与倾斜状态下计算K值及p值变化对比表(山谷地形) Table 3 K value calculated when the electrode spacing of line 100 is horizonal and gradient and comparison of P value ・58・ 贵州地质 2013年3O卷 表4 114线电极距呈水平状态与倾斜状态下计算K值及p值变化对比表【山脊地形) Table 4 K value calculated when the electrode spacing of line 1 14 is horizonal and gradient and comparison of P value 第1期 彭洪军,等:基于电法装置电极排列系数K值计算方法的改进 ・59・ (1)不论是山谷地形、山脊地形,还是在斜 出现数值的变化影响,造成山谷或山脊及两侧电 坡,地形高程差变化,都会使每个测点的斜电极距 阻率数值的数倍增长变化,总是P斜大于P水,但 比水平电极距要大。以AM为例,斜电极距有 是对电阻率曲线的形态变化没有多大影响。可 +( 一 +( 三 ) (4), 见,改进后的K值计算方法,在山谷或山脊地形 对水平电极距有AMp=J(Ax—Mx) +(A 一M ) 的高程差变化加大了对K值的影响(斜坡影响较 (8),而(Az一 ) ≥0,可知AM, A 。 小),但取得的K值是接近客观实际的。 (2)在表3的山谷地形影响下,山谷位于557 本文中的所有坐标值都是工作设计中的,然 —567号点之间,极值段一带装置系数K斜一K水 而在实际工作中,由于客观条件和测量仪器精度 变化值为7 024~13 916,左侧斜坡位于548— 的影响,往往实际定点坐标与设计坐标会有偏差, 557号点变化值为2 364~9 139,右侧斜坡位于 就会导致AM、AN、BM、BN并不是设计的距离,传 567—580号点变化值为1 786~15 723;则对应 统方法是按照工作设计人为输入AM、AN、BM、 的变化电阻率P斜一P水,山谷一带为(1 480~ BN各点的距离带人公式(2)中求得K值,计算误 3 358)Q.m,左侧斜坡为(120~868)Q.1TI,右侧 差很大。而用改进的方法,是根据GPS或全站仪 斜坡为(42—567)Q.m。可见,对各区K斜一K 实测点位三维坐标值计算K值误差较小,更为合 水即使有一定数值的变化,但是在各区分段对应 理。同样因为AB两极移动也会导致K值出现偏 的电阻率P斜一P水变化差异并不大。 差,在此就不多作陈述。 (3)在表4的山脊地形影响下,山脊位于496 本文列举的实例主要说明地形对K值的影 —506号点之间,极值段一带装置系数K斜一K水 响程度,提示物探工作者重视地形对K值的影 变化值为5 754~16 242,左侧斜坡位于485— 响,应该对K值进行校正。如今校正K值的方法 496号点变化值为11~840,右侧斜坡位于506— 却不多,其中大部分又过于繁琐,在山区的电法勘 517号点变化值为4 443~15 213;则对应的变化 探中实用性不强。笔者通过工作实践总结出应用 电阻率P斜一P水,山脊一带为(22"-59)Q.nl,左侧 电极点坐标参数值结合Excel计算K值的方法技 斜坡为(0~14)Q.m,右侧斜坡为(104~172) 术,该方法技术极大地提高了山区K值计算的工 Q.m。可见,对各区K斜一、K水即使有一定数值的 作效率。 变化,但是在各区分段对应的电阻率P斜一P水变 化差异并不大。 4 结束语 (4)表3、表4中的P斜一P水的分区段变化值 和图2、图3电阻率曲线展布,体现K斜一K水虽然 本文的宗旨在体现利用电极点坐标参数值求 ・60・ 贵州地质 2013年3O卷 K值有许多优点,其优点是十分明显的,其一是准 发展的迫切需求产生的。这种方法相对于传统方 确性,各坐标值来源于测区控制测量、测网测量, 法的诸多优势,让笔者有理由相信这种方法在贵 坐标值真实的反应剖面地形,所求得的K值就是 州矿区的广泛应用是科技进步与时代发展的 实际工作的K值;其二是实用性,此方法不仅可 必然。 以降低地形对K值的影响,更重要的是实用于常 规电法勘探的任意装置,即规则装置、非规则装 [参考文献] 置、电极位移等都可根据实际坐标值求得准确的 [1]刘国兴.电法勘探原理与方法[M].北京:地质出版 K值;其三是普遍性,相对于传统方法只适用于水 社,2005. 平、单一地形中,该方法适用于任何复杂地形;其 [2] 贵州省地矿局102地质大队.贵州省从江县那哥铅锌矿西 矿段电法勘探报告[R].2010. 四是简捷性,此方法在室内处理数据时,只需要导 [3]贵州省地矿局102地质大队.贵州省从江县平忙铅锌多金 入各测点坐标值,填入到计算模板中即可求得K 属矿区电法勘探报告[R].2011. 值,处理过程简单快捷、易于操作。 应用电极点坐标参数值求K值是顺应时代 Improvement of Calculation Method Based on Electrode Array Coefifcient K of Two Resistivities PENG Hong—jun,ZHANG Ying—wen,DAI Di,MA Qiang,DAI Can—yi (Geologic Party 106,Guizhou Bureau of eGology&Mineral Exploration and Development,Zu ̄yi 563003,Guizhou,China) [Abstract]The calculation is complex and has abundant quantity when K value is adjusted by conventional method.By USe electrode coordinate parameter,space analytic geometry and Excel computing platform,it makes adjustment of K value which changes for electrode displacement and topographic relief.This method is accurate,simple,eonnivent,fast and easy to do,it’S very useful for electrical prospecting in the mountain area. [Key words] K vlaue;coordinate;topography;surface distance;actual distance;resistivity (上接第77页) A New Geophy sical Technology:Microtremor Survey Technology WANG Hong (Ge ̄cid and Geochemical Exploratoin Party,Guizhou Nonferrous Metals Geological Exploration Bureau,Duyun,Guizhou 558000,Guizhou,China) [Abstract] Microtremor survey technology is a new survey technology which developed on hte basis of tradi— tional mierotremor depth sounding by Doc.XU Pei—fen and others of Institute of Geology and Geophysics of Chinese Academy of Science in recent years,it has been used in many exploration fields.This technology is USe seismic detector to receive the waves from every direction on surface,take the Rayleigh wave dispersion curves by SPAC method and get S—wave speed by inversion.The technology won’t be influenced by magnetic nad noise,the survey depth is deep.Alhtough there is some limiattion now,the advantages shows this teeh— nology will be a promising new method. [Key words] mierotremor survey;Rayleigh wave;inversion;formaiton wave speed structure;survey 

因篇幅问题不能全部显示,请点此查看更多更全内容