第３９卷第６期 ２０１５年１２月 物探与化探 Ｖｏ１．３９，Ｎｏ．６ Ｄｅｃ．，２０１５ ＧＥＯＰＨＹＳＩＣＡＬ＆ＧＥ０ＣＨＥＭＩＣＡＬ ＥＸＰＬＯＲＡＴＩＯＮ ｄｏｉ：１０．１ １７２０／ｗｔｙｈｔ．２０１５．６．１２ 郭博，潘和平，杨怀杰，等．地一一井方式激电确定矿体空间位置［Ｊ］．物探与化探，２０１５，３９（６）：１１６７—１１７５．ｈｔｔｐ：／／ｄｏｉ．ｏｒｇ／１０．１１７２０／ｗｔｙｈｔ．２０１５． ６．１２ Ｇｕｏ Ｂ，Ｐａｎ Ｈ Ｐ，Ｙａｎｇ Ｈ Ｊ，ｅｔ １．Ｔｈｅ ａｐｐｌａｉｃａｔｉｏｎ ｏｆ ｔｈｅ ｓｕｒｆａｃｅ・ｔｏ—ｂｏｒｅｈｏｌｅ ｉｎｄｕｃｅｄ ｐｏｌａｒｉｚａｔｉｏｎ ｍｅｔｈｏｄ ｔｏ ｄｅｔｅｒｍｉｎｉｎｇ ｔｈｅ ｓｐ￣ｉｌ ｌａｏｃａｔｉｏｎ ｏｆ ｔｈｅ ｏｒｅｂｏｄｙ ［Ｊ］．Ｇｅｏｐｈｙｓｉｃａｌ ａｎｄ Ｇｅｏｃｈｅｍｉｃａｌ Ｅｘｐｌｏｒａｔｉｏｎ，２０１５，３９（６）：１１６７—１１７５．ｈｔｔｐ：／／ｄｏｉ．ｏｒｇ／ｌＯ．１１７２０／ｗｔｙｈｔ．２０１５．６．１２ 地一井方式激电确定矿体空间位置 郭博，潘和平，杨怀杰，姚静 （中国地质大学地球物理与空间信息学院，湖北武汉４３００７４） 摘要：为更好地解决井旁矿体的勘探及盲矿的预测问题，探讨了地一井方式激发极化法对于井旁的规则形体（矿 体）的激电影响；采用解析投影法，以球体和板状体两种模型为例，遵循“单井一多井”的激电规律来判别矿体的空 间位置。最后，归纳总结出地一井方式激电在实际工作中的使用条件，能够快速判别井旁盲矿体的空间位置。结论 指出。合理的施工方式以及适量的观测数据对预测矿体的产状具有积极作用。 关键词：地一井方式：井中激发极化法；深部找矿；井旁矿体 中图分类号：Ｐ６３１ 文献标识码：Ａ 文章编号：１０００—８９１８（２０１５）０６—１１６７—０９ 地一井方式激发极化法是一种在实际矿体勘探 井方式激电指的是将供电电极Ａ、Ｂ布置在地面，其 中Ａ极置于距井口一定距离（Ｌ）处或者置于井口（￡ ＝０），Ｂ极置于无穷远，测量电极 、Ⅳ置于井中（图 １）。 中十分流行的物探方法，具有分辨率高、操作方便等 优点．引起很多学者的关注。国内从２０世纪７０年 始。黄智辉等人就对地一井方式方位测量资料 解释问题进行了探讨…：吴至善等人对井中激发极 化法的理论与实践做出了研究Ｌ２ ；吕玉增、张国艳 等人对地井方位激电异常正反演研究做了很大贡 献［３－７ ３；潘和平、周峰等人对井中激电井旁球体的正 反演进行了深入研究 ．拓宽了其在井中物探的 应用范围。国外的相关研究大多是以三维可视化的 角度．运用一些复杂原理例如有限元的方法对地一 井方式激电进行正反演Ｉ１卜坞］，取得了很好的效果。 此次在前人研究的基础上完成了论证工作，用 ｍａｆｌａｂ实现球体和板状体的激电模型计算．达到了 预期的效果。 因浮土巨目灰岩 图花岗岩■●矿体 图１地一井激电工作方式示意 １工作原理及方法技术 １．１工作原理 地一井方式激电利用井的天然优势使得测量端 接近矿体，矿体激电异常更容易被探测到。在测量 操作中，地面上不同方位或不同距离的 极布置会 改变矿体的极化方向和强度。可以把Ａ极依次布 置在东、南、西、北各个方向上，在每个方位做测量． 通过对各个方向上的激电异常曲线形态和强弱的差 井中激发极化法可分为地一井、井一地、井一井三 种方式，其中地一井方式［】 在近年来的金属矿勘查 中取得了较好的效果 ］。 总体而言，井中激电同地面激电原理一样，都是 以岩矿石的激发极化机理为基础，包括（电子导体） 超电压的形成和（离子导体）薄膜极化的产生。地一 收稿日期：２０１５—０３—２４ 异来推断井旁盲矿体的相关地质信息［１ ２ｏ］。 基金项目：国家自然科学基金项目（４１０７４０８６）；国家高技术研究发展计划（８６３计划）课题（２０ｌ４ＡＡ０６Ａ６０８） ６期 郭博等：地一井方式激电确定矿体空间位置 ・１１６９・ 剖面的对应点，并且均为四方位供电点的角平分点， 圆心角均为４５。。 需要说明的是，Ａ 和Ａ 供电点其极化方向在 垂直主剖面方向的分量效果为零。则主剖面方向分 模拟采用偶极子等效电流场＿１ ，表达式为 ＝Ｋ・（ｃｏｓＯ／Ｒ ）， 其中， 为激电系数，Ｒ为球心至观测点ＭＮ的距 离，０为变化角。在ｘＯｚ坐标中对ｚ微分，得到二次 场异常电位差表达式［】 ： 量投影效果相当于井口Ａ 供电［２ ，最终的思路是 将八方位辅助供电点投影到主剖面上。此时要注 意．投影点 和 两点的等效极化倾角ｉ的数学关 系可以由投影点相对于球心的距离和中心埋深ｈ求 出。若供电点固定，则Ａ Ａ 等效投影距离相同， △ ＝Ｋ 等 Ｏ Ｍ—Ｎ。 Ａ Ａ 的距离也相同，因此对应方位的激发极化异 常曲线趋势两两相同。这点从图４中得到了验证。 ｌＯＯ ２００ ３００ —＼ １ ， 磕 、、 ，，， ＼ ，一、 ／ 一 ｏ Ａ”ｌ Ｉ’ｌ ，，， 口’７ ， 、、 ＼／ 三 ，、、 ， 。 ＼． — 图４ Ｌ－－２００ ｍ时的八方位供电点球体二次场电位差曲线 ２．２单井板状体 为“反ｓ形”，而直立板状体和离井倾斜板状体（板 顶近井）除外，曲线呈“正ｓ形”。 ２）反方位（Ｌ＝一２００ ｉｎ）特征曲线异常减弱，这 主剖面供电是关于板状体研究里最基础的切人 点，因此依然遵循球体的研究思路。即：将图２中球 形矿体模型变为板状体。图５为计算得到的Ｌ＝２００ ｍ时得到的不同产状二次场电位差曲线。 １）主方位（ ＝２００ｍ）特征曲线异常明显，幅值 是因为对于板状体的极化强度减弱，曲线异常幅值 减小，对于所有的产状的板状体特征曲线形态都是 上负下正的“正Ｓ形”。可以根据特征曲线的幅值 比反方位要大。板在井轴上的投影都落在主方位 △ 曲线的正异常区域内。当板延伸有限时．主方 位的△ 曲线基本是上正下负的“反Ｓ形”。当板 强弱和符号变化特征来推断矿体的方位。 ３）板状体延伸越大，主方位上的特征曲线正异 常区沿其延伸方向伸展也越大。延伸很大的板状体 曲线形态主要由其近端决定的，所以当板状体向着 无限延伸时，向井倾斜板（板尾近井）的特征曲线仍 ６期 郭博等：地一井方式激电确定矿体空间位置 ３多井△ 曲线特征 本节研究内容是上文传统单井ＩＰ情况下的延 伸。多井是指虚拟井，是为了满足本次研究的正演 模拟实验。 首先，建立多井ＩＰ模型（图６），以矿体为中心， “正井口一矿体一反井口”为主线，采用同一个对应 剖面来进行研究．这样使得所在同一个剖面上的测 量参数有对应性。除此之外，还应注意将环绕矿体 四周的４个（或８个）钻井看作是以矿体中心在地 面的投影点为圆心，以ｄ为矿体距井距离而均匀分 布。为了方便起见，将规则几何体分为球体和板状 体两类．且各个钻井的距井距离ｄ保持不变。 ３．１球状矿体激电特征曲线 环绕其打４个（或８个）井，然后，分别以每个井为 观测所在位置来总结其规律。 无论是球体还是板状体，初步的研究思路是以 关于多井环绕矿体的激电问题，关键是要找准 主剖面．此时需要对每一个井进行不同供电距离多 次测量，进行模型分类，考虑到球体的各项同性，则 激电特征曲线只研究某一个剖面上不同供电距离的 Ｔ、 ／ － ｌ ｌ ｌ Ｉ Ｉ Ｉ ：・ ・ ： Ｉ ● ．・｛． ： ▲ Ｉ 模型。多井球情况下主剖面的程序参数设定见表 １，计算结果见图７。 表１ 多井球主剖面在不同供电距离下的程序参数 、＿＿ ）未知矿体▲ 图６多井情况下矿体的井口布置示意 △ ／Ｖ △ ／Ｖ Ａ ／Ｖ 一ＩＩＩＩ 主方位供电 ——井口供电 反方位供电 一（ａ）Ｌ＞ｄ （ｃ）Ｌ＜ｄ 图７多井情况下（Ｌ＞ｄ）、（Ｌ＝ｄ）、（Ｌ＜ｄｌ时的二次场异常曲线 在Ｌ＞Ｊ情况下，可以看出主方位供电的曲线呈 “反Ｓ形”，井口供电与主方位曲线关于 轴对称， 反方位曲线呈现“正ｓ”形，特征分辨明显。 在 ＝ｄ情况下，主方位供电点Ａ 刚好位于球 体中心的正上方，主方位激电曲线的辨认度较Ｌ＞ｄ Ｌ＜Ｊ情况下的辨别度是三种中最差的．因此实 际工作中最好不要这样布置供电电极的位置。使供 电距离 大于等于矿体距井距离ｄ。 综上，在多井条件下，供电电极距 最好是大于 或者等于矿体的距井距离，方便探测。这样的施工 方式必定要使供电方式多样化和数据多样化．才能 约束测井结果的不确定性。 ３．２板状体激电特征曲线 的情况下降低。当主方位供电点Ａ，位于球体中心 点正上方时，此时的曲线辨别度较好。主方位的红 色曲线由于二次场激电垂直向下。故关于Ｚ＝１００ ｍ对 称；而随着供电点向反方位移动，曲线逐渐呈现“正 在多井情况下，板状体的供电电极布置方位与 ｓ”形，反方位供电比井Ｅｌ供电明显。 球体是一致的（图８），但是考虑到板状体的各相异 物探与化探 ３９卷 Ｊ２ Ｔ、 ／ ・－Ｉ Ｊ４ ｌ．ｒ＂ － ．・｝Ｊ８ ： ：： ｌ ｌ Ｉ Ｉ 二 三 摊二 图８多井情况下板状体的井口布置示意 性（板长、板厚、板倾斜、延伸等的差异），不同方位 的钻井所测得的激发极化特性曲线比球体情况复杂 得多。由于前面一节探讨的供电电极距条件，因此 设定的参数为Ｌ＞ｄ。 为了方便计算该模型。使其在东西方向倾斜 （对于Ｊ１为离井倾斜，对于Ｊ２为向井倾斜），借鉴单 井情况下板状体的激电问题。同样，其中红色曲线 代表主方位供电。绿色曲线代表井口供电。蓝色曲线 代表反方位供电。 ３．２．１正剖面（Ｊ１一中心投影点一Ｊ２） 为了很清楚地进行对比，将与矿体中心点在地 面上的投影连线的对应矿井作为一个剖面（图９）。 需要注意的是板状体的倾斜方向对于Ｊ１是离井方 向．而对于Ｊ２是向井方向，这样就可以和单井情况 下的板状体激发极化特性曲线相结合起来了。表３ 给出了儿、Ｊ２的主剖面供电参数，图１０为２口井的 特征曲线。 中心投影点 图９多井板状体正剖面示意（Ｊ１一中心投影点－－Ｊ２） 表３ Ｊ１Ｊ２的主剖面供电参数 井名 底端井距／ｍ探测深度／ｍ 板体状态 板厚／ｍ Ｊ１ ８０ ５００ 离井 １Ｏ Ｊ２ ６０ ５００ 向井 １Ｏ Ａ ／ｖ △ ，ｖ 图１Ｏ Ｊ１、Ｊ２特性曲线 ３．２．２副剖面（Ｊ３一中心投影点＿Ｊ４） 因为板状体在副剖面上的激电等效几何形状为 直立板状体。故可以运用单井情况下的板状体曲线 进行模拟（图１１，图１２）。Ｊ３和Ｊ４的主剖面供电参 数为：供电距Ｌ＝２００ ｍ，顶端埋深ｈ ＝１５０ ｍ，底端埋 深ｈ，＝２５０ ｍ，顶端井距ｄ ＝７０，底端井距ｄ ＝７０ ｍ， 探测深度 ＝５００ ｍ，板体直立，板厚１０ ｍ。 中心投影点 图ｎ 多井板状体副剖面示意（Ｊ３一中心投影点—Ｊ４） 图１２ Ｊ３和ｊ４（竖直）特征曲线 ３．２．３斜剖面（Ｊ５一中心投影点一Ｊ６、Ｊ７一中心投影 点一Ｊ８） 首先要说明的是。研究Ｊ５、Ｊ６、Ｊ７和Ｊ８井的激 电曲线情况需要考虑到三维板状体，这将会比二维 ６期 郭博等：地一井方式激电确定矿体空间位置 ・１１７３・ 板状体要更加复杂。解决这种三维几何体的模型其 的思维来解决斜剖面的模型问题。需要说明的是， 实最好的方法是运用数值模拟的方法来进行处理． 由于对称性，其对应的２个斜剖面（Ｊ７一中心投影 但是矩阵理论过于复杂，因此寻求另外一种比较简 点一Ｊ８）的情况是一样的，这次一起进行阐述。 便的方法来阐述，即投影法。 由图１４比较可知，斜剖面（Ｊ５一中心投影点一 由于一个三维倾斜板状体在竖直非垂直剖面的 Ｊ６，Ｊ７一中心投影点一Ｊ８）情况下，其板状体模型倾 投影倾向均与原来一致（图１３），因此可以采取定性 向不变，倾角介于直立与倾斜之间。 Ａ ／Ｖ △ ／Ｖ 中心投黟点 Ａ７２ Ｊ７ ｌ ＡＴＩ Ｊ８ ●Ａａ２ Ａ６２ Ｊ６ Ａ５１ Ａ６１ Ｊ５ Ａ５２ Ｉ＿Ｉ 图１３板状体在斜剖面方位的投影模型 图１４ Ｊ６．＿Ｊ７、Ｊ５—Ｊ８特征曲线 ３．２．４小结 定矿体所在方位，最后划定主方位后可采用多种切 １）对于中心矿体为板状体的多井问题，不能等 线法确定矿体的距井距离和中心埋深１２４－２７］。如何 同于球体来解决。因为板状体的各项异性，使得各 准确判定矿体的所在方位将对后续解释工作至关 个方位的井口所测得的激电曲线的影响要素与距井 重要。 距离ｄ、供电点空间分布和板状体的产状有关。 ４．１矿体不同方位曲线特征 ２）从图中的曲线可以很清楚地看出矿体的倾 矿体的激电特征曲线受矿体几何分布的影响。 斜方向，曲线最大极值的绝对值表示板状体的近井 对于各向同性的球形矿体，可以通过最大幅度的 端．并且应该注意特征曲线随着矿井方位的不同而 “反Ｓ”形曲线确定矿体主方位．而对于情况较为复 产生有规律的变化。 杂的各向异性的板状矿体．需要重点讨论。 ３）如果板状体是向上或者向下无限延伸，则可 该板状体走向为正南北方向，倾向正西，其顶、 以根据单井情况找出规律，无限延伸段收敛于零值。 底埋深分别为１５０．，２５０ ｍ，矿体中心距离各个方位 ４不同方位激电曲线特征及适用条件 井为５０ ｍ（图１５）。 从模拟图像可以看出，极大值点和极小值点分 多井情况下判别矿体空间位置的首要任务是确 别反应矿体顶端和底部深度，受到供电电流方位影 △ ／Ｖ －０．２ Ｏ Ｏ．２ Ｏ ｌ０Ｏ ２００ ３００ ４０Ｏ 特征曲发 特征曲线 Ｓ００ 特征曲线 特征曲线 特征曲线 图１５多方位井的特征曲线 物探与化探 ３９卷 响，极大值点深度即为板状体顶端深度，极小值点深 度比真实值略大。特征曲线从Ｊ２开始．逐步过渡到 Ｊ６／Ｊ７，Ｊ３／Ｊ４，Ｊ５／Ｊ８，最后到Ｊ１，其极小值幅值异常逐 渐增大，极值点深度位置略微升高，极大值幅值异常 也逐渐增大。通过这一点可以总结出判断板状体空 间方位的方法。即：从特征曲线的极大值和极小值相 对大小以及异常区域平缓程度来进行： １）Ｊ２井为主剖面板状体向井方位。激电特征曲 线正值区域较为平缓； ２）Ｊ６／Ｊ７井为斜剖面向井方位，激电特征曲线 正值区域平缓，极大值点深度变浅： ３）Ｊ３／Ｊ４井为板状体走向方位。激电特征曲线 极小值绝对值等于极大值； ４）Ｊ５／Ｊ８井为斜剖面离井方位，激电特征曲线 极小值绝对值略小于极大值； ５）Ｊ１井为主剖面板状体离井方位，激电特征曲 线极小值绝对值明显小于极大值。 因此，在地一井方式激电施工时，应该注意找出 极值点特征最大的特征曲线，即可确定主剖面方位。 对于其余方位井特征曲线，可以加以修正判断出矿 体空间位置 ４．２适用条件 １）确定矿体的与钻孔的相对位置，注意虚拟多 钻孔的特征判别。 ２）注意供电电极的施工布局，其不同方位和供 电距离将会直接影响激电特征曲线形态。 ３）确定主剖面后，先进行井口测量，若激电曲 线特征满足预期再进行多方位确认观测，进而对井 旁矿体的埋深、距离、见矿深度及延伸方向进行判 断。此次用球体和板状体两个简单模型进行正演模 拟，可作为很好的正演拟合模型，因此结合实测数据 可以对相应的矿体进行反演解释。 ５结语 １）地一井方式激电在井中物探领域有广阔的应 用前景，能够很好地判断矿体所在方位、球形矿体中 心埋深深度和板状体近井端埋深深度。本研究从简 单的规则模型出发，模拟出了不同供电方位、不同矿 井方位的激发极化特征曲线，对于辅助方位的激电 特征情况，采用几何投影的方法，最后通过不同情况 下的激电特征曲线对矿体的产状进行预判，再结合 实际的地质、钻孔和物探资料。用多种手段联合分析 矿体的赋存情况。 ２）对于井旁的未知矿体，由于其钻孔相对于矿 体的方位待定，因此可以通过研究环绕矿体的多个 虚拟方位来对其特征曲线进行比较分析．从而对其 空间分布进行判断，最终为实际勘探提供理论借鉴 的依据。 ３）通过研究，对前人的有关成果进行了一系列 验证，在地一井方式井中激电的正演模拟计算中，主 剖面（正方位一井口一反方位）上所得的曲线特征 符合预期情况。值得说明的是，球形矿体特征曲线 形态较为明显，而无限延伸板状体的特征曲线因其 一端具有收敛的趋势同样能够容易被甄别，但是本 次研究表明特征曲线只能对于不同倾向的板状体敏 感，而对于同一倾向情况下不同倾角的特征曲线异 常形态不明显。研究还表明，在实际现场勘探中应 尽可能用较多数量且分布均匀的供电点对地下矿体 进行施工。 ４）对于非规则矿体的激电特征曲线情况．建议 采用有限元的思想来处理，作为下一步研究的方向。 参考文献： ［１］ 黄智辉．井中激电地一井方式方位测量资料解释问题的探讨 ［Ｊ］．物探与化探，１９７９，１（３）：２２—３０． ［２］吴至善．井中激发极化法的理论与实践［Ｊ］．长春地质学院学 报，１９７６（４）：８６—１２３． ［３］ 吕玉增．地一井、井一地ＩＰ三维快速正反演研究［Ｄ］．湖南：中 南大学．２００８． ｆ４］ 吕玉增，阮白尧，彭苏萍．地一井方位激电观测异常特征研究 ［Ｊ］．地球物理学进展，２０１２，２７（１）：２０１—２１６． ［５］ 张国艳，吕玉增，刘正，等．地一井断面ＩＰ正演模拟研究［Ｊ］．工 程地球物理学报，２０１３，ｌＯ（１）：４６—５０． ［６］施国良，张国雄．宏观场论［Ｍ］．武汉：中国地质大学出版 社．１９８７． ［７］ 张国艳．地一井激发极化法有限元正演模拟与反演成像［Ｄ］．广 西：桂林理工大学，２０１２． ［８］ 周峰，潘和平，吴国平，等．井中激电地一井方式井旁球体正反 演［Ｊ］．物探与化探，２００８，３２（３）：３２１—３２５． ［９］ 周峰，潘和平，文，等．基于有限差分的井中激发极化三维 正演［Ｊ］．电波科学学报，２０１０，２５（４）：７８５—７９２． ［１０］Ｐａｎ Ｈ Ｐ，Ｍｅｎｇ Ｑ Ｍ．Ｓｕｒｆａｃｅ ｔｏ ｂｏｒｅｈｏｌｅ ｉｎｄｕｃｅｄ ｐｏｌａｒｉｚａｔｉｏｎ ｍｅａｓｕｒｅｍｅｎｔ ｄａｔａ ｉｎｖｅｒｓｉｏｎ ｗｉｔｈ ｏｐｔｉｍｉｚｅｄ ｓｃａｌｅ ｔｒａｎｓｆｏｒｍａｔｉｏｎ ＢＦＧＳ ｍｅｔｈｏｄ［Ｒ］．８ｔｈ Ｉｎｔｅｒｎａｔｉｏｎａｌ Ｃｏｎｆｅｒｅｎｃｅ ａｎ Ｎａｔｕｒａｌ Ｃｏｍ— ｐｕｔ．ｉｏｎ，２０１２． ［１ １］Ｂａｒｂａｒａ Ｉｎａ Ａｎｄｅｒｓｏｎ．Ｍｏｄｅｌｉｎｇ ａｎｄ ｉｎｖｅｒｓｉｏｎ ｍｅｔｈｏｄｓ ｆｏｒ ｔｈｅ ｉｎ－ ｔｅｒｐｒｅｔａｔｉｏｎ ｏｆ ｒｅｓｉｓｔｉｖｉｔｙ ｌｏｇｇｉｎｇ ｔｏｏｌ ｒｅｓｐｏｎｓｅ［Ｄ］．Ｗｅｓｔｅｎｒ Ｃｏｎ— ｎｅｅｔｉｏｎ Ｓｔａｔｅ Ｕｎｉｖｅｒｓｉｔｙ ｇｅｂｏｒｅｎ ｔｅ Ｄａｎｂｕｒｙ，Ｃｏｎｎｅｃｔｉｃｕｔ， ＵＳＡ，２００１． ［１２］Ｄｅｙ Ａ，Ｍｏｒｒｉｓｏｎ Ｈ Ｆ．Ｒｅｓｉｓｔｉｖｉｔｙ ｍｏｄｅｌｉｎｇ ｆｏｒ ａｒｂｉｔｒａｒｉｌｙ ｓｈａｐｅｄ ｔｈｒｅｅ．ｄｉｍｅｎｓｉｏｎａｌ ｓｔｒｙｃｔｕｒｅｓ［Ｊ］，１９７７（１０）． ［１３］Ｍｉｎｓｌｅｙ Ｂ Ｊ．Ｍｏｄｅｌｉｎｇ ａｎｄ ｉｎｖｅｒｓｉｏｎ ｏｆ Ｓｅｌｆ－Ｐｏｔｅｎｔｉａｌ Ｄａｔａ［Ｊ］．Ｂ Ｓ Ａｐｐｌｉｅｄ Ｐｈｙｓｉｃｓ，Ｐｕｒｄｕｅ Ｕｎｉｖｅｒｓｉｔｙ，１９９７． ［１４］ｗａｌｔ Ｊ Ｒ，Ｇｒｕｓｚｋａ Ｔ Ｐ．Ｒｅｓｉｓｔｉｖｉｔｙ ａｎｄ ｐｏｌａｒｉｚａｔｉｏｎ ｒｅｓｐｏｎｓｅ ｆｏｒ ａ ｂｏｒｅｈｏｌｅ ｍｏｄｅｌ［Ｊ］．Ｔｒａｎｓａｃｔｉｏｎｓ ｏｎ ｇｅｏｓｅｉｅｎｅｅ ａｎｄ ｒｅｍｏｔｅ ６期 ｓｅｎｓｉｎｇ，ＧＥ－２５，１９８７（３）． 郭博等：地一井方式激电确定矿体空间位置 版社 １９８３． ・１１７５・ ［１５］Ｅｖｅｒｅｔｔ Ｍ Ｅ，Ｂａｄｅａ Ｅ Ａ，Ｓｈｅｎ Ｌ Ｃ．３－Ｄ Ｆｉｎｉｔｅ ｅｌｅｍｅｎｔ ａｎａｌｙｓｉｓ ｏｆ ［２１］吴冰．地一井激电法探测盲矿体的数字模拟研究［Ｄ］．西安：长 安大学。２００８． ［２２］潘和平，马火林，祁明松，等．井中激发极化法模型实验设计 ［Ｊ］．实验技术与管理，２０１０，２７（１）：３１—３３． ｉｎｄｕｃｔｉｏｎ ｌｏｇｇｉｎｇ ｉｎ ａ ｄｉｐｐｉｎｇ ｆｏｒｍａｔｉｏｎ［Ｊ］．Ｔｒａｎｓａｃｔｉｏｎｓ ｏｎ ｇｅｏ— ｓｃｉｅｎｃｅ ａｎｄ ｒｅｍｏｔｅ ｓｅｎｓｉｎｇ，２００１，３９（１０）． ［１６］Ｓｐｉｔｚｅｒ Ｋ，Ｃｈｏｕｔｅａｕ Ｍ．Ａ ＤＣ ｒｅｓｉｓｔｉｖｉｔｙ ａｎｄ ＩＰ ｂｏｒｅｈｏｌｅ ｓｕｒｖｅｙ ａｔ ｔｈｅ ｃａｓｅ ｂｅｒａｒｄｉ ｇｏｌｄ ｍｉｎｅ ｉｎ ｎｏｒｔｈｗｅｓｔｅｒｎ［Ｊ］．Ｇｅｏｐｈｓｉｃｓ，２００３ （６８）：４５３—４６３． 『１７］Ｂｅｖｃ Ｄ．Ｍｏｒｒｉｓｏｎ Ｈ Ｆ．Ｂｏｒｅｈｏｌｅ—ｔｏ—ｓｕｒｆａｃｅ ｅｌｅｃｔｒｉｃａｌ ｒｅｓｉｓｔｉｖｉｔｙ ｍｏ． ［２３］潘和平，马火林，祁明松．井中物探教学内容探讨与改革［Ｊ］．中 国地质教育，２００９，（２）：８３—８７． ［２４］潘和平．井中激发极化法在矿产资源勘探中的作用［Ｊ］．物探与 化探，２０１３，３７（４）：６２０ ６２６． ｎｉｔｏｒｉｎｇ ｏｆ ａ ｓａｌｔ ｗａｔｅｒ ｉｎｊｅｃｔｉｏｎ ｅｘｐｅｒｉｍｅｎｔ［Ｊ］．Ｇｅｏｐｈｙｓｉｃｓ，１９９１， ５６：７６９—７７７． ［２５］宋振龄．井中激发极化法在金属勘查中的应用［Ｊ］．河南地质， １９９２，１０（１）：６４—６９． ［１８］Ｓｌａｔｅｒ Ｌ，Ｂｉｎｌｅｙ Ａ Ｍ，Ｄａｉｌｙ Ｗ，ｅｔ ａ１．Ｃｒｏｓｓ－ｈｏｌｅ ｅｌｅｃｔｒｉｃａｌ ｉｍａ ｎｇ ｏｆ ａ ｃｏｎｔｒｏｌｌｅｄ ｓａｌｉｎｅ ｔｒａｃｅｒ，ｉｎｊｅｃｔｉｏｎ［Ｊ］．Ｊｏｕｒｎａｌ ｏｆ Ａｐ— ｐｌｉｅｄ Ｇｅｏｐｈｙｓｉｃｓ，２０００，（４４）：８５—１０２． ［２６］王耀升，刘斌．井中激电法在内蒙查干布拉根银铅锌多金属矿 的应用［Ｊ］．矿产勘查，２０１３（４）：４５３—４５７． ［２７］高长荣．井中激电在西霞矿区的应用［Ｊ］．物探与化探，２００７，３１ （Ｓ）：９８—１０１． ［１９］潘和平，马火林，蔡柏林，等．地球物理测井与井中物探［Ｍ］．北 京：科学出版社，２００９． ［２Ｏ］蔡柏林，黄智辉，谷守民．井中激发极化法［Ｍ］．北京：地质出 Ｔｈｅ ａｐｐｌｉｃａｔｉｏｎ ｏｆ ｔｈｅ ｓｕｒｆａｃｅ－ｔｏ－ｂｏｒｅｈｏｌｅ ｉｎｄｕｃｅｄ ｐｏｌａｒｉｚａｔｉｏｎ ｍｅｔｈｏｄ ｔｏ ｄｅｔｅｒｍｉｎｉｎｇ ｔｈｅ ｓｐａｔｉａｌ ｌｏｃａｔｉｏｎ ｏｆ ｔｈｅ ｏｒｅｂｏｄｙ ＧＵＯ Ｂｏ，ＰＡＮ Ｈｅ—Ｐｉｎｇ，ＹＡＮＧ Ｈｕａｉ－Ｊｉｅ，ＹＡＯ Ｊｉｎｇ （Ｉｎｓｔｉｔｕｔｅ ｏｆ Ｇｅｏｐｈｙｓｉｃｓ ａｎｄ Ｇｅｏｍａｔｉｃｓ，Ｃｈｉｎａ Ｕｎｉｖｅｒｓｉｔｙ ｏｆ Ｇｅｏｓｃｉｅｎｃｅｓ，Ｗｕｈａｎ ４３００７４，Ｃｈｉｎａ） Ａｂｓｔｒａｃｔ：Ｔｈｅ ｓｕｒｆａｃｅ—ｔｏ—ｂｏｒｅｈｏｌｅ ＩＰ ｍｅｔｈｏｄ，ａｓ ａｎ ａｃｔｕａｌｌｙ ｗｉｄｅｌｙ ｕｓｅｄ ｏｒｅ ｅｘｐｌｏｒａｔｉｏｎ ｍｅａｎｓ，ｈａｓ ａｒｏｕｓｅｄ ａ ｌｏｔ ｏｆ ｉｎｔｅｒｅｓｔ ａｍｏｎｇ ｇｅｏｌｏｇｉｓｔｓ ｂｏｔｈ ｉｎ Ｃｈｉｎａ ａｎｄ ａｂｒｏａｄ．Ｉｎ ｔｈｉｓ ｐａｐｅｒ．ｔｈｅ ａｕｔｈｏｒｓ ｐｒｏｂｅｄ ｉｎｔｏ ｔｈｅ ＩＰ ｅｆｆｅｃｔｓ ｏｆ ｔｈｅ ｓｕｒｆａｃｅ—ｔｏ—ｗｅｌｌｓ ＩＰ ｍｅｔｈｏｄ ｏｎ ｔｈｅ ｒｅｇｎ— ｌｙ ｏｒｅｂｏａｄｙ ｂｅｓｉｄｅ ｔｈｅ ｗｅｌｌ ＳＯ ａｓ ｔｏ ｓｏｌｖｅ ｔｈｅ ｐｒｏｂｌｅｍ ｏｆ ｐｒｅｄｉｃｔｉｎｇ ｔｈｅ ｅｘｐｌｏｒａｔｉｏｎ ｏｆ ｔｈｅ ｏｒｅｂｏｄｙ ｂｅｓｉｄｅ ｔｈｅ ｗｅｌｌ ａｎｄ ｔｈｅ ｂｌｉｎｄ ｏｒｅｂｏｄｙ． Ｉｎ ｔｈｉｓ ｐａｐｅｒ，ｔｈｅ ａｕｔｈｏｒｓ ｍａｉｎｌｙ ｓｔｕｄｉｅｄ ｔｈｅ ｓｐｈｅｒｉｃａｌ ａｎｄ ｔｈｅ ｐｌａｔｙ ｏｒｅｂｏｄｉｅｓ ｏｆ ｄｉｆｆｅｒｅｎｔ ＩＰ ｍｏｄｅｓ ｏｎ ｔｈｅ ｇｅｎｅｒａｔｅｄ ｃｕｒｖｅ，ａｎｄ ｔｈｅｎ ｅｘｐｌｏｒｅｄ ｍｕｌｔｉ—ｗｅｌｌ ｃａｓｅｓ．Ｆｉｎａｌｌｙ。ｔｈｅ ａｕｔｈｏｒｓ ｄｉｓｃｕｓｓｅｄ ｈｏｗ ｔｏ ｄｅｔｅｒｍｉｎｅ ｔｈｅ ｓｐａｔｉａｌ ｌｏｃａｔｉｏｎ ｏｆ ｔｈｅ ｏｒｅｂｏｄｙ ｂａｓｅｄ ｏｎ ｍｕｌｔｉ—ｗｅｌｌｓ ｌａｗ． Ｔｈｅ ｒｅｓｕｌｔｓ ｏｂｔａｉｎｅｄ ｂｙ ｔｈｅ ａｕｔｈｏｒｓ ａｙｅ ｈｅｌｐｆｕｌ ｔｏ ｔｈｅ ｒｅｌａｔｅｄ ｒｅｓｅａｒｃｈ ｆｉｅｌｄｓ．Ｉｎ ｃｏｎｃｌｕｓｉｏｎ，ｒｅａｓｏｎａｂｌｅ ｃｏｎｓｔｒｕｃｔｉｏｎ ｍｅｔｈｏｄ ａｎｄ ｐｒｏｐｅｒ ｑｕａｎｔｉｔｙ ｏｆ ｏｂｓｅｒｖｅｄ ｄａｔａ ｈａｖｅ ｐｏｓｉｔｉｖｅ ｅｆｆｅｃｔ ｏｎ ｔｈｅ ｐｒｏｇｎｏｓｉｓ ｏｆ ｔｈｅ ｏｒｅｂｏｄｙ． Ｋｅｙ ｗｏｒｄｓ：ｓｕｒｆａｃｅ—ｔｏ—ｂｏｒｅｈｏｌｅ ｓｔｙｌｅ；ｂｏｒｅｈｏｌｅ ｉｎｄｕｃｅｄ ｐｏｌａｒｉｚａｔｉｏｎ；ｄｅｅｐ ｐｒｏｓｐｅｃｔｉｎｇ；ｔｈｅ ｒｕｌｅｓ ｏｆ ｇｅｏｍｅｔｒｙ 作者简介：郭博（１９９０一），男，硕士研究生，研究方向为地球物理测井与井中物探。 通讯作者：潘和平（１９５３－），男，博士，教授，博士生导师，主要研究方向为地球物理测井与井中物探。
