我国煤矿水文地质条件十分复杂,水害防治形势愈发严峻,水源性隐蔽致灾体的准确探测仍将是今后防治水工作的重点和难点。煤矿井下水害事故主要发生在巷道掘进期间和工作面回采期间。虽然物探方法在以往实践中提前发现灾害性地质异常体,为保障矿井安全生产发挥巨大作用,但也存在对巷道前方探测距离不够长(一般在100m左右)、假异常多,对工作面穿透距离短(一般不超过200m)、远处分辨率不够等缺点。另外,在回采工作面顶底板含水层进行注浆改造时,对小体积的导水性陷落柱、岩溶、断层等水源性隐蔽致灾体的精确探测存在技术缺口。针对当前技术手段的不足,研究一种在煤矿井下水平钻孔中接收感应二次场、地面回线磁源激发一次场的瞬变电磁探测技术(简称井-地瞬变电磁法)。地面发射源可将回线源和电流尽量加大,增加发射磁矩以尽量激发地下水源性隐蔽致灾体;水平钻孔内布置的测点距离异常体更近,最大限度减少二次场的距离损失。理论上,井-地瞬变电磁法对水源性隐蔽致灾体具有更强的分辨能力。在煤矿巷道掘进前方超前钻孔中应用该项技术,能实现利用单个钻孔长距离准确探查掘进前方水源性隐蔽致灾体的目的;在工作面钻孔中应用该项技术,能对钻孔周围漏掉的水源性隐蔽致灾体进行定位,提高探查准确度和效率。论文以载流回线磁源激发场分布和异常体感应涡流场三分量特征为研究对象,以公式推导、三维时域有限差分数值模拟、理论分析等为研究方法,对异常体的井-地瞬变电磁法响应特征、感应涡流场三分量的空间指向性、浮动系数空间交汇算法和最小二乘反演算法进行研究,实现基于井-地瞬变电磁法的异常体三维空间定位技术。论文取得的主要研究成果如下:(1)推导了激发场的解析解公式,获得了激发场的分布特征和回线源最佳尺寸匹配依据。从毕奥-萨伐尔定律出发,对圆形回线源、长方形回线源和正方形回线源分别进行磁感应强度解析解的推导。对各种激发源载入直流电后在空间形成的磁场分布进行数值计算,分析总结其分布特征。依据回线源在地下不同深度产生磁场的强度分布特征,获得最佳回线源尺寸的选择办法。三种激发源的激发场在空间的分布特征基本一致,激发场强度总体呈浅部强、深部弱的分布特征,磁感应强度主要受与回线源的距离控制,在回线附近的磁场表现为最强值。随着与回线源之间距离的增加,磁感应强度变化逐渐减慢。在与回线源一定距离的平面上,磁感应强度分布呈中间强、边缘弱的特征,等值线依据回线源形状表现为圆形或椭圆形。在回线源中间部分区域,激发场主要为垂向磁场且强度在横向上的变化较小,这是地面瞬变电磁法将回线中心一定区域当作均匀场的原因。以实际生产中常用的正方形回线为例,计算8种常见尺寸回线在地下不同深度产生的磁场。大尺寸回线源在浅部产生的激发场并非更强,但随深度变化相对稳定,且深部的激发场相对更强。在某个固定深度的目标体处,过大或过小尺寸的回线源并不能激发出最大强度的一次磁场。计算结果对井-地瞬变电磁法回线源尺寸的选择具有指导意义。(2)获得了煤矿典型水源性隐蔽致灾体的响应特征,和主要参数改变对三分量响应的影响规律。采用三维时域有限差分算法,对构建的含水陷落柱、小煤窑积水采空区、含导水断层和顶板砂岩富水体等煤矿典型水源性隐蔽致灾体的全空间三分量响应进行数值模拟。各模型的异常体均有明显的三分量磁异常响应,异常测道曲线均以“N”或“V”为基本形状,过零点和极值点分别指向异常体中心在测线上的位置。对回线源尺寸、低阻覆盖层、异常体方位、异常体规模和距离这五方面因素的影响进行数值模拟。回线源尺寸的增加能提高总场和异常场的响应强度,但增加幅度逐渐减弱,且异常场面临幅值极限。低阻覆盖层使得异常响应幅值随覆盖层电阻率的降低而减弱。三分量异常曲线形态组合与异常体相对测线的空间方位具有单一的匹配性,可据此对异常体方位进行判断。异常体与测线距离的变化,只改变异常曲线的幅值,不改变曲线的形态。距离测线越远,异常信号越弱。以仪器精度和大地电磁噪声为基本阈值,确定可分辨最小信号强度。从测线深度改变时总场三个分量的强度与最小信号阈值之间的相对关系,探讨了井-地瞬变电磁法的极限探测深度。(3)以水平电流环辐射磁场各分量矢量的空间指向性为基础,开发出适用于水平测线不同Y偏移距的浮动系数空间交汇算法。通过分析水平电流环辐射磁场的三个分量在XY平面、XZ平面上,水平测线处各测点的三分量矢量分布,确定各分量矢量与电流环中心在XY平面、XZ平面上具有明确的指向性。采用时域有限差分算法模拟异常体的异常场,分析感应涡流磁场随时间的分布规律,通过与水平电流环产生磁场的进行对比,认为异常体在外部产生的异常场,可以用位于异常体内部的电流环所辐射的磁场来代替。研究水平测线上磁场三分量矢量对水平电流环中心的交汇特性,开发出适用于水平测线的浮动系数空间交汇算法,并建立基于不同Y偏移距的浮动系数表。在进行XZ平面交汇时,根据系数表对Z分量进行自适应调节,可使XZ平面的交汇结果更准确。分别位于水平测线不同象限四个模型的试算结果,验证了浮动系数空间交汇算法的准确性。采用此算法对不同倾斜角度的近似水平电流环进行空间定位试算,计算结果说明浮动系数空间交汇算法更适用于水平电流环,倾斜角度的增加会导致交汇准确度的降低。(4)以分离的异常场三分量数据为基础,采用基于电流环理论的最小二乘算法,反演异常体的尺寸、空间姿态和中心坐标。以水平测线正常区段采集的数据为背景场,采用多项式拟合算法对测线异常段的背景场进行曲线拟合,获得全测线的背景场。通过总场数据减去背景场的方式得到异常场。借助三维直角坐标系的旋转公式实现任意倾斜角度的电流环正演。视异常体内部的涡流场为一个等效电流环,赋予其中心点坐标、半径、倾斜角度等变量,使用带约束的最小二乘反演算法对各变量进行反演计算,以最小拟合误差为导向对电流环的参数不断迭代计算,获得异常体中心坐标、空间姿态、规模大小等参数。分别以MAXWELL软件对板状体的正演、时域有限差分算法对立方体的数值模拟、以地面铝板为异常体的现场试验,共三种方式取得的数据为基础,对反演算法进行验证。结果显示,对板状体数据反演得到的电流环中心坐标、空间姿态、尺寸均与模型参数吻合;对立方体数据反演得到的电流环中心坐标和空间姿态与模型参数吻合较好,但尺寸结果不稳定,与模型参数存在一定偏差;对现场试验数据反演得到的电流环中心坐标和空间姿态与模型参数吻合较好,但同样存在尺寸结果不稳定,存在一定偏差的现象。