２０１３年３Ｏ卷 贵 州 地 质 Ｖｏ１．３０ Ｎｏ．１（Ｔｏ１．１１４）２０１３ ・５４・ ＧＵＩＺＨＯＵ ＧＥＯＬＯＧＹ 第１期（总第１１４期） 基于电法装置电极排列系数Ｋ值计算方法的改进 彭洪军，张应文，代 迪，马 强，代璨怡 （贵州省地矿局１０２地质大队，贵州遵义５６３００３） ［摘要］应用常规方法校正Ｋ值时计算复杂、计算工作量大。采用电极点坐标参数（ｘ，Ｙ，ｚ）， 应用空间解析几何原理并借助Ｅｘｃｅｌ（电子表格应用软件）计算平台计算Ｋ值，实现了因电极位 移、地形起伏导致Ｋ值变化的校正。该方法准确、简单、方便、快捷、可操作性强，对山区的电法 勘探具有较强的实用性。 ［关键词］Ｋ值；坐标；地形；水平距离；实际距离；电阻率 ［中图分类号］Ｐ６３ｌ＿３ ４［文献标识码］Ａ［文章编号］１０００—５９４３（２０１３）０１—００５４—０７ 坐标参数值就可求出Ｋ值，同时实现电极位移Ｋ １ 前言 值的校正，该方法在水平地形、单一地层、综合地 层中都可运用，同时能提高数据处理效率，对山区 近年来电阻率法在山区的电法勘探中得到了 的电法勘探具有实用性。 广泛的应用，然而在山区的电法勘探中由于受客 观条件限制导致电极发生位移，还有在高程落差 ２原理及工作方法 极大的山谷、山脊地段，Ｋ值变化通常超出允许误 差范围，从而使视电阻率曲线产生畸变，因此有必 要对Ｋ值进行校正。 ２．１ 原理 在电阻率法勘探中，比较常用的方法有二极 在地下介质均匀各向同性的条件下，电阻率 装置、三极装置、联合剖面装置、对称四极装置、偶 计算式Ｐｓ＝Ｋ・ａ ，（１），（１）式中的Ｋ称为电 极装置、中间梯度装置等，这些方法都涉及电极排 极排列系数，Ｋ值数学表达式Ｋ＝２７ｒ／（１／ＡⅣ一１／ 列系数Ｋ值。传统的Ｋ值计算方法一般基于水 ＡＮ一１／ＢＭ＋１／ＢＮ）（２），其值是一个仅与各电极 平距离、对称剖面，计算比较简单，在地势比较单 间空间位置有关的量；在实际工作中，电极形式可 一平坦的地区十分实用。 多种多样，但电极的排列形式和电极距一经确定， 然而笔者在实际工作中发现，多数情况下的 Ｋ值便可由公式（２）计算出来；（２）式为Ｋ值的计 Ｋ值计算看似简单，其实计算起来比较复杂，比如 算通式，在计算联合剖面、三极测深等装置时Ｂ 计算中间梯度旁测剖面、非对称四极、电极发生位 极的处理有两种办法：一是将通式中的ＢＭ、ＢＮ 移等就极为复杂。在山区的电阻率法勘探中地形 作“。。”大处理，二是将ｃ极当成Ｂ极输入电极点 起伏对Ｋ值影响较大，尤其是高程落差极大的山 坐标值；（２）式中的ＡＭ、ＡＮ、ＢＭ、ＢＮ是各电极距 谷、山脊地形比较突出，如果仍按水平距离计算Ｋ 的距离，也是本文要求解的变量，关于ＡＭ、ＡＮ、 值，通常会导致Ｋ值变化值超出允许误差。基于 ＢＭ、ＢＮ距离的求解方法请看图１。 上述原因，为达到Ｋ值计算准确、简单、方便、快 假设Ａ、Ｂ、Ｍ、Ｎ各点坐标分别为Ａ（Ａ ，Ａｙ， 捷，实现因电极位移、地形起伏导致Ｋ值变化的 Ａｚ）、Ｂ（曰 ， ，Ｂｚ）、Ｍ（朋 ， ， ）、Ｎ（以，Ⅳ ， 校正，这里应用电极点坐标参数值借助Ｅｘｃｅｌ平 Ⅳ ）。以求解ＡＭ为例，根据空间解析几何求距 台对Ｋ值求解，该方法不仅简单、方便、快捷、准 离原理，利用Ａ、Ｍ 、ＡＭ组成的平面直角三角形 确、可操作性强，更重要的是具有较强的灵活性， 应用勾股定理求出平面斜距Ｌｘｙ。 即在任意装置下只需导入Ａ、Ｂ、Ｍ、Ｎ四个电极点 ［收稿日期］２０１２—１２—２２［改回日期】２０１３—０３—０１ ［作者简介］彭洪军（１９８７一），男，山东临沂人，物探助理工程师，主要从事地球物理勘探工作。 第１期 彭洪军，等：基于电法装置电极排列系数Ｋ值计算方法的改进 ・５５・ Ｈ∞ （Ｏ 踮 Ｏ ∞ Ｏ ；Ｏ２ Ｏ Ｏ 图１任意装置Ｋ值计算原理示意图 Ｆｉｇ．Ｉ Ｋ ｖａｌｕｅ ｃｏｕｎｔｉｎｇ ｐｒｉｎｃｉｐｌｅ ｏｆ ａｎｙ ａｒｒａｙ Ｌｘｒ＝ ＝ （３） ２．２计算方法 再利用Ａ、Ｍ、Ｍ 组成的垂面直角三角形应 用勾股定理求出空间斜距ＡＭ ，其中ＡＭ 计算公 ２．２．１建立模板 式为： 本文的Ｋ值计算是借助Ｅｘｃｅｌ平台下的函数 ＡＭｓ＝ 巧＝ （４） 计算功能，即可以在表格相应的单元格内编辑输 同理可求得ＡＮ ＢＭ ＢＮ。计算公式为： 入相应的计算公式进行相关计算，为了能方便使 用、减少重复输入、提高计算效率，因此建立一个 Ａ＾ ＝√（Ａ —Ｎｘ） ＋（Ａｙ—Ｎｒ） ＋（Ａｚ—Ｎｚ） （５） Ｅｘｃｅｌ模板很有必要。 ＢＭｓ＝√（ －Ｍｘ） ＋（Ｂｙ一 ） ＋（Ｂｚ－Ｍｚ） （６） 我们计算的Ｋ值模板见表１（此模板在Ｅｘｃｅｌ ＢＮｓ＝￣／（ 一 ） ＋（ 一Ⅳｙ） ＋（ 一Ⅳｚ） （７） 中总体由Ａ—ＡＤ ３０个列和１、２、３、４的行构成， 在不考虑高程差的情况下，计算水平距离 这里由于篇幅有限，在Ｗｏｒｄ文件中横向无法全 ＡＭＰ、ＡＮＰ、ＢＭＰ、ＢＮＰ公式简化如下： 部展开，故将表格一分为二转换为纵向展示，即将 Ａ ＝Ｊ（Ａｚ一』Ｉ ） ＋（Ａｙ— ｙ） （８） Ａ—Ｏ共计１５列置于表格上半部分、Ｏ列以右的 Ｐ—ＡＤ共计１５列置于表格下半部分）。模板内 ＡＮｔ，＝√（Ａ —Ｎｘ） ＋（Ａｒ—Ｎｙ） （９） 容包括：电极及记录点编号；Ａ、Ｍ、Ｎ、Ｂ（ｃ）电极 ＢＭｔ，＝√（日ｘ—Ｍｘ） ＋（ ｙ— ｙ） （１ｏ） 点坐标；水平电极距ＡＭＰ、ＡＮＰ、ＢＭＰ、ＢＮＰ；水平距 曰Ⅳ尸＝√（ 一心） ＋（ ｙ一Ⅳｙ） （１１） 离计算Ｋ值（以下简称水平Ｋ值）；斜电极距 从上述（４）一（１１）计算公式中可以看出，参 ＡＭｓ、ＡＮｓ、ＢＭｓ、ＢＮｓ；斜电极距计算Ｋ值（以下简 与计算两点距离的参数仅为两点的坐标值，这样 称实际Ｋ值）；ＡＫ（即实际Ｋ值与水平Ｋ值之 就建立起了两点的空间距离与两点坐标参数值的 差，以下简称△ｋ）；Ｋ偏离度。其中第４行为计算 关系。 模板，Ａ、Ｍ、Ｎ、Ｂ（Ｃ）点坐标值系任意充填值，只 是为了编辑公式计算需要，无实际意义，从第５行 起以下为计算区域。 ・５６・ 贵州地质 ２０１３年３０卷 在建立完成模板框架结构之后，可在相应的 单元格中进行输入计算公式，计算公式详见表２。 表２ Ｅｘｃｅｌ计算公式 Ｔａｂｌｅ ２ Ｃａｌｃｕｌａｔｉｏｎ ｆｏｒｍｕｌａ ｏｆ Ｅｘｃｅｌ 在单元格内输入Ｅｘｃｅｌ表达公式后，点击编 储于其它介质中的各测点坐标数据按模板格式复 辑栏中的“、／”符号，一个公式输入计算完成，所 制到对应的单元格内，数据导人完成。 有公式输入完成后将该文件保存为模板，Ｋ值计 ２．２．３数据计算 算模板建立完成。 数据计算是在计算模板中利用Ｅｘｃｅｌ强大的 ２．２．２数据导人 计算功能，仅需要简单操作即可完成计算。具体 数据导人有手工输入和文件复制两种形式， 操作过程为分别点击选择模板行（第４行）对应 手工导入效率低而且容易出错，因此常用的数据 单元格，将鼠标对准右下角黑点使其变成“＋”向 导人形式为文件复制。数据导人方法是在模板中 下拖动，即可实现对模板公式的复制和数据计算。 编制完记录点及Ａ、Ｍ、Ｎ、Ｂ（Ｃ）电极编号后，将存 第１期 彭洪军，等：基于电法装置电极排列系数Ｋ值计算方法的改进 －５７・ 线为中间梯度旁测剖面，与主剖面１０４线平行相距 ３应用实例 本应用实例取材于某铅锌矿区，矿区位于贵 州省从江县内，区内地形切割强烈、地形起伏连 绵、坡度较陡、山脊与冲沟交互出现、悬崖与陡坎 ４００ ｍ，该测线地形呈“Ｖ”字型，山顶（５４８点）最大 高程８６６．２ ｍ，山谷（５６６点）最小高程６６６．３ ｍ，最 大高程差２００ ｍ，测点编号５４５—５８３，共计３９个点； 山脊地形１１４线为中间梯度主剖面，该测线地形为 山脊地形，山顶（５００点）最大高程１２２５．５ｍ，山谷 遍布，最大相对高程差可达数百米。 （５１７点）最小高程９１３．１ ｍ，最大高程差３１２．４ ｍ， 区内作过中间梯度测量、测深测量，测量结果 测点编号４８２—５２０，共计３９个点。 都存在地形起伏对Ｋ值的影响，也都存在供电电 将该实例数据导人后，应用Ｅｘｃｅｌ模板计算， 极位移对Ｋ值的影响。这里以中间梯度装置，山 可得到如表１格式中的水平ｋ值、实际ｋ值、△ｋ 谷、山脊地形为例来展示应用电极点坐标参数值计 值和Ｋ偏离度值，对计算数据用绘图软件作出的 算Ｋ值过程及地形对Ｋ值的影响程度。中间梯度 相关参数曲线图，如图２、图３所示（图中增加了 装置参数（１）ＡＢ＝４ ０００ ｍ；（２）ＭＮ＝１２０ ｍ；（３）点 一条２％水平直线用于对Ｋ值偏离度的研判参 距２０ ｍ；（４）测网密度２００×２０ ｍ。山谷地形１００ 考，当Ｋ值误差变化在２％以上时需作校正）。 塞 ～ ＼— 由翻 多 … 。一 … ＼。 图２ １００线模板计算各参数曲线图（山谷地形） 图３ １１４线模板计算各参数曲线图（山脊地形） Ｆｉｇ．２ Ｃｕｒｖｅ ｏｆｐａｒａｍｅｔｅｒｓ ｃａｌｃｕｌａｔｅｄ ｂｙｆｏｒｍｗｏｒｋ ｏｆｌｉｎｅ１００ Ｆｉｇ．３ Ｃｕｒｖｅ ｏｆｐａｒａｍｅｔｅ￣ｃａｌｃｕｌａｔｄｅ ｂｙｆｏｒｍｗｏｒｋｏｆｌｉｎｅ１０４ 当电极距呈水平状态与倾斜状态时，由模板 值，与Ｋ值变化、Ｐ值变化对比结果如下表３、表 计算的电极距不同状态装置系数Ｋ值，及其相应 ４。进一步可以看出山谷、山脊及斜坡地形，对Ｋ 工作电流（ｉｐ）、工作电压（ｖｐ）下的对应电阻率Ｐ 值以及视电阻率变化的影响。 表３ １００线电极距呈水平状态与倾斜状态下计算Ｋ值及ｐ值变化对比表（山谷地形） Ｔａｂｌｅ ３ Ｋ ｖａｌｕｅ ｃａｌｃｕｌａｔｅｄ ｗｈｅｎ ｔｈｅ ｅｌｅｃｔｒｏｄｅ ｓｐａｃｉｎｇ ｏｆ ｌｉｎｅ １００ ｉｓ ｈｏｒｉｚｏｎａｌ ａｎｄ ｇｒａｄｉｅｎｔ ａｎｄ ｃｏｍｐａｒｉｓｏｎ ｏｆ Ｐ ｖａｌｕｅ ・５８・ 贵州地质 ２０１３年３Ｏ卷 表４ １１４线电极距呈水平状态与倾斜状态下计算Ｋ值及ｐ值变化对比表【山脊地形） Ｔａｂｌｅ ４ Ｋ ｖａｌｕｅ ｃａｌｃｕｌａｔｅｄ ｗｈｅｎ ｔｈｅ ｅｌｅｃｔｒｏｄｅ ｓｐａｃｉｎｇ ｏｆ ｌｉｎｅ １ １４ ｉｓ ｈｏｒｉｚｏｎａｌ ａｎｄ ｇｒａｄｉｅｎｔ ａｎｄ ｃｏｍｐａｒｉｓｏｎ ｏｆ Ｐ ｖａｌｕｅ 第１期 彭洪军，等：基于电法装置电极排列系数Ｋ值计算方法的改进 ・５９・ （１）不论是山谷地形、山脊地形，还是在斜 出现数值的变化影响，造成山谷或山脊及两侧电 坡，地形高程差变化，都会使每个测点的斜电极距 阻率数值的数倍增长变化，总是Ｐ斜大于Ｐ水，但 比水平电极距要大。以ＡＭ为例，斜电极距有 是对电阻率曲线的形态变化没有多大影响。可 ＋（ 一 ＋（ 三 ） （４）， 见，改进后的Ｋ值计算方法，在山谷或山脊地形 对水平电极距有ＡＭｐ＝Ｊ（Ａｘ—Ｍｘ） ＋（Ａ 一Ｍ ） 的高程差变化加大了对Ｋ值的影响（斜坡影响较 （８），而（Ａｚ一 ） ≥０，可知ＡＭ， Ａ 。 小），但取得的Ｋ值是接近客观实际的。 （２）在表３的山谷地形影响下，山谷位于５５７ 本文中的所有坐标值都是工作设计中的，然 —５６７号点之间，极值段一带装置系数Ｋ斜一Ｋ水 而在实际工作中，由于客观条件和测量仪器精度 变化值为７ ０２４～１３ ９１６，左侧斜坡位于５４８— 的影响，往往实际定点坐标与设计坐标会有偏差， ５５７号点变化值为２ ３６４～９ １３９，右侧斜坡位于 就会导致ＡＭ、ＡＮ、ＢＭ、ＢＮ并不是设计的距离，传 ５６７—５８０号点变化值为１ ７８６～１５ ７２３；则对应 统方法是按照工作设计人为输入ＡＭ、ＡＮ、ＢＭ、 的变化电阻率Ｐ斜一Ｐ水，山谷一带为（１ ４８０～ ＢＮ各点的距离带人公式（２）中求得Ｋ值，计算误 ３ ３５８）Ｑ．ｍ，左侧斜坡为（１２０～８６８）Ｑ．１ＴＩ，右侧 差很大。而用改进的方法，是根据ＧＰＳ或全站仪 斜坡为（４２—５６７）Ｑ．ｍ。可见，对各区Ｋ斜一Ｋ 实测点位三维坐标值计算Ｋ值误差较小，更为合 水即使有一定数值的变化，但是在各区分段对应 理。同样因为ＡＢ两极移动也会导致Ｋ值出现偏 的电阻率Ｐ斜一Ｐ水变化差异并不大。 差，在此就不多作陈述。 （３）在表４的山脊地形影响下，山脊位于４９６ 本文列举的实例主要说明地形对Ｋ值的影 —５０６号点之间，极值段一带装置系数Ｋ斜一Ｋ水 响程度，提示物探工作者重视地形对Ｋ值的影 变化值为５ ７５４～１６ ２４２，左侧斜坡位于４８５— 响，应该对Ｋ值进行校正。如今校正Ｋ值的方法 ４９６号点变化值为１１～８４０，右侧斜坡位于５０６— 却不多，其中大部分又过于繁琐，在山区的电法勘 ５１７号点变化值为４ ４４３～１５ ２１３；则对应的变化 探中实用性不强。笔者通过工作实践总结出应用 电阻率Ｐ斜一Ｐ水，山脊一带为（２２＂－５９）Ｑ．ｎｌ，左侧 电极点坐标参数值结合Ｅｘｃｅｌ计算Ｋ值的方法技 斜坡为（０～１４）Ｑ．ｍ，右侧斜坡为（１０４～１７２） 术，该方法技术极大地提高了山区Ｋ值计算的工 Ｑ．ｍ。可见，对各区Ｋ斜一、Ｋ水即使有一定数值的 作效率。 变化，但是在各区分段对应的电阻率Ｐ斜一Ｐ水变 化差异并不大。 ４ 结束语 （４）表３、表４中的Ｐ斜一Ｐ水的分区段变化值 和图２、图３电阻率曲线展布，体现Ｋ斜一Ｋ水虽然 本文的宗旨在体现利用电极点坐标参数值求 ・６０・ 贵州地质 ２０１３年３Ｏ卷 Ｋ值有许多优点，其优点是十分明显的，其一是准 发展的迫切需求产生的。这种方法相对于传统方 确性，各坐标值来源于测区控制测量、测网测量， 法的诸多优势，让笔者有理由相信这种方法在贵 坐标值真实的反应剖面地形，所求得的Ｋ值就是 州矿区的广泛应用是科技进步与时代发展的 实际工作的Ｋ值；其二是实用性，此方法不仅可 必然。 以降低地形对Ｋ值的影响，更重要的是实用于常 规电法勘探的任意装置，即规则装置、非规则装 ［参考文献］ 置、电极位移等都可根据实际坐标值求得准确的 ［１］刘国兴．电法勘探原理与方法［Ｍ］．北京：地质出版 Ｋ值；其三是普遍性，相对于传统方法只适用于水 社，２００５． 平、单一地形中，该方法适用于任何复杂地形；其 ［２］ 贵州省地矿局１０２地质大队．贵州省从江县那哥铅锌矿西 矿段电法勘探报告［Ｒ］．２０１０． 四是简捷性，此方法在室内处理数据时，只需要导 ［３］贵州省地矿局１０２地质大队．贵州省从江县平忙铅锌多金 入各测点坐标值，填入到计算模板中即可求得Ｋ 属矿区电法勘探报告［Ｒ］．２０１１． 值，处理过程简单快捷、易于操作。 应用电极点坐标参数值求Ｋ值是顺应时代 Ｉｍｐｒｏｖｅｍｅｎｔ ｏｆ Ｃａｌｃｕｌａｔｉｏｎ Ｍｅｔｈｏｄ Ｂａｓｅｄ ｏｎ Ｅｌｅｃｔｒｏｄｅ Ａｒｒａｙ Ｃｏｅｆｉｆｃｉｅｎｔ Ｋ ｏｆ Ｔｗｏ Ｒｅｓｉｓｔｉｖｉｔｉｅｓ ＰＥＮＧ Ｈｏｎｇ—ｊｕｎ，ＺＨＡＮＧ Ｙｉｎｇ—ｗｅｎ，ＤＡＩ Ｄｉ，ＭＡ Ｑｉａｎｇ，ＤＡＩ Ｃａｎ—ｙｉ （Ｇｅｏｌｏｇｉｃ Ｐａｒｔｙ １０６，Ｇｕｉｚｈｏｕ Ｂｕｒｅａｕ ｏｆ ｅＧｏｌｏｇｙ＆Ｍｉｎｅｒａｌ Ｅｘｐｌｏｒａｔｉｏｎ ａｎｄ Ｄｅｖｅｌｏｐｍｅｎｔ，Ｚｕ￣ｙｉ ５６３００３，Ｇｕｉｚｈｏｕ，Ｃｈｉｎａ） ［Ａｂｓｔｒａｃｔ］Ｔｈｅ ｃａｌｃｕｌａｔｉｏｎ ｉｓ ｃｏｍｐｌｅｘ ａｎｄ ｈａｓ ａｂｕｎｄａｎｔ ｑｕａｎｔｉｔｙ ｗｈｅｎ Ｋ ｖａｌｕｅ ｉｓ ａｄｊｕｓｔｅｄ ｂｙ ｃｏｎｖｅｎｔｉｏｎａｌ ｍｅｔｈｏｄ．Ｂｙ ＵＳｅ ｅｌｅｃｔｒｏｄｅ ｃｏｏｒｄｉｎａｔｅ ｐａｒａｍｅｔｅｒ，ｓｐａｃｅ ａｎａｌｙｔｉｃ ｇｅｏｍｅｔｒｙ ａｎｄ Ｅｘｃｅｌ ｃｏｍｐｕｔｉｎｇ ｐｌａｔｆｏｒｍ，ｉｔ ｍａｋｅｓ ａｄｊｕｓｔｍｅｎｔ ｏｆ Ｋ ｖａｌｕｅ ｗｈｉｃｈ ｃｈａｎｇｅｓ ｆｏｒ ｅｌｅｃｔｒｏｄｅ ｄｉｓｐｌａｃｅｍｅｎｔ ａｎｄ ｔｏｐｏｇｒａｐｈｉｃ ｒｅｌｉｅｆ．Ｔｈｉｓ ｍｅｔｈｏｄ ｉｓ ａｃｃｕｒａｔｅ，ｓｉｍｐｌｅ，ｅｏｎｎｉｖｅｎｔ，ｆａｓｔ ａｎｄ ｅａｓｙ ｔｏ ｄｏ，ｉｔ’Ｓ ｖｅｒｙ ｕｓｅｆｕｌ ｆｏｒ ｅｌｅｃｔｒｉｃａｌ ｐｒｏｓｐｅｃｔｉｎｇ ｉｎ ｔｈｅ ｍｏｕｎｔａｉｎ ａｒｅａ． ［Ｋｅｙ ｗｏｒｄｓ］ Ｋ ｖｌａｕｅ；ｃｏｏｒｄｉｎａｔｅ；ｔｏｐｏｇｒａｐｈｙ；ｓｕｒｆａｃｅ ｄｉｓｔａｎｃｅ；ａｃｔｕａｌ ｄｉｓｔａｎｃｅ；ｒｅｓｉｓｔｉｖｉｔｙ （上接第７７页） Ａ Ｎｅｗ Ｇｅｏｐｈｙ ｓｉｃａｌ Ｔｅｃｈｎｏｌｏｇｙ：Ｍｉｃｒｏｔｒｅｍｏｒ Ｓｕｒｖｅｙ Ｔｅｃｈｎｏｌｏｇｙ ＷＡＮＧ Ｈｏｎｇ （Ｇｅ￣ｃｉｄ ａｎｄ Ｇｅｏｃｈｅｍｉｃａｌ Ｅｘｐｌｏｒａｔｏｉｎ Ｐａｒｔｙ，Ｇｕｉｚｈｏｕ Ｎｏｎｆｅｒｒｏｕｓ Ｍｅｔａｌｓ Ｇｅｏｌｏｇｉｃａｌ Ｅｘｐｌｏｒａｔｉｏｎ Ｂｕｒｅａｕ，Ｄｕｙｕｎ，Ｇｕｉｚｈｏｕ ５５８０００，Ｇｕｉｚｈｏｕ，Ｃｈｉｎａ） ［Ａｂｓｔｒａｃｔ］ Ｍｉｃｒｏｔｒｅｍｏｒ ｓｕｒｖｅｙ ｔｅｃｈｎｏｌｏｇｙ ｉｓ ａ ｎｅｗ ｓｕｒｖｅｙ ｔｅｃｈｎｏｌｏｇｙ ｗｈｉｃｈ ｄｅｖｅｌｏｐｅｄ ｏｎ ｈｔｅ ｂａｓｉｓ ｏｆ ｔｒａｄｉ— ｔｉｏｎａｌ ｍｉｅｒｏｔｒｅｍｏｒ ｄｅｐｔｈ ｓｏｕｎｄｉｎｇ ｂｙ Ｄｏｃ．ＸＵ Ｐｅｉ—ｆｅｎ ａｎｄ ｏｔｈｅｒｓ ｏｆ Ｉｎｓｔｉｔｕｔｅ ｏｆ Ｇｅｏｌｏｇｙ ａｎｄ Ｇｅｏｐｈｙｓｉｃｓ ｏｆ Ｃｈｉｎｅｓｅ Ａｃａｄｅｍｙ ｏｆ Ｓｃｉｅｎｃｅ ｉｎ ｒｅｃｅｎｔ ｙｅａｒｓ，ｉｔ ｈａｓ ｂｅｅｎ ｕｓｅｄ ｉｎ ｍａｎｙ ｅｘｐｌｏｒａｔｉｏｎ ｆｉｅｌｄｓ．Ｔｈｉｓ ｔｅｃｈｎｏｌｏｇｙ ｉｓ ＵＳｅ ｓｅｉｓｍｉｃ ｄｅｔｅｃｔｏｒ ｔｏ ｒｅｃｅｉｖｅ ｔｈｅ ｗａｖｅｓ ｆｒｏｍ ｅｖｅｒｙ ｄｉｒｅｃｔｉｏｎ ｏｎ ｓｕｒｆａｃｅ，ｔａｋｅ ｔｈｅ Ｒａｙｌｅｉｇｈ ｗａｖｅ ｄｉｓｐｅｒｓｉｏｎ ｃｕｒｖｅｓ ｂｙ ＳＰＡＣ ｍｅｔｈｏｄ ａｎｄ ｇｅｔ Ｓ—ｗａｖｅ ｓｐｅｅｄ ｂｙ ｉｎｖｅｒｓｉｏｎ．Ｔｈｅ ｔｅｃｈｎｏｌｏｇｙ ｗｏｎ’ｔ ｂｅ ｉｎｆｌｕｅｎｃｅｄ ｂｙ ｍａｇｎｅｔｉｃ ｎａｄ ｎｏｉｓｅ，ｔｈｅ ｓｕｒｖｅｙ ｄｅｐｔｈ ｉｓ ｄｅｅｐ．Ａｌｈｔｏｕｇｈ ｔｈｅｒｅ ｉｓ ｓｏｍｅ ｌｉｍｉａｔｔｉｏｎ ｎｏｗ，ｔｈｅ ａｄｖａｎｔａｇｅｓ ｓｈｏｗｓ ｔｈｉｓ ｔｅｅｈ— ｎｏｌｏｇｙ ｗｉｌｌ ｂｅ ａ ｐｒｏｍｉｓｉｎｇ ｎｅｗ ｍｅｔｈｏｄ． ［Ｋｅｙ ｗｏｒｄｓ］ ｍｉｅｒｏｔｒｅｍｏｒ ｓｕｒｖｅｙ；Ｒａｙｌｅｉｇｈ ｗａｖｅ；ｉｎｖｅｒｓｉｏｎ；ｆｏｒｍａｉｔｏｎ ｗａｖｅ ｓｐｅｅｄ ｓｔｒｕｃｔｕｒｅ；ｓｕｒｖｅｙ