1 引言
随着我国经济的迅速发展，水利工程行业充满了发展机遇，同时也面临严峻考验 [1]。在水利工程的设计和施工阶段，避开地下管线是至关重要的一个环节。地下管线（如石油管道、燃气管道、国防通信管线、高压电力管线等）承担着输送石油、天然气、通信、电力等功能，是城市发展的基础 [2]。地下管线的复杂性高、信息化程度低、管理不统一等，难以通过资料搜集和整理获得准确的地下管线信息。特别是施工阶段，挖断管线时有发生。不仅存在重大安全隐患，而且造成经济损失 [3]。如何获得准确地下管线信息，成为了一个重要课题。
随着科学技术的发展，地下管线的探测方法也愈加多样，主要包括直接开井测量和物探方法间接测量。使用物探方法探测地下管线的物性基础在于：目标管线和周围介质存在物性差异，如介电常数、磁导率、电导率、密度等。物探方法如高密度电法、井中磁法、瞬变电磁法等 [4]，受到场地和测区环境的限制，应用场景有限，管线仪法和探地雷达法仍是基础、高效、便捷的地下管线探测方法。文章将从探测原理、探测方法、数值模拟和实际应用几方面进行探讨。
2 管线仪的工作原理和应用
2.1 管线仪的探测原理
常用的管线仪如 RD8100 和 RD8000 等，探测原理为电磁感应，即：管线仪发射机发射电磁信号，在金属管线中会产生感应电流 ( 也称为等效电流 )，感应电流往周围激发出变化电磁场，接收机接收到该电磁信号。根据接收到的电磁信号，就可以判定目标管线的位置和埋深等信息 [2，3]。
2.2 管线仪的探测方法
根据管线仪发射机激发信号的方式不同，可以将常用的探测方法分为 3 种：夹钳法、直连法、感应法 [3]。
（1）夹钳法是指将管线仪配套的环形夹钳闭合地夹在管线上，使得管线穿过夹钳。主要用于通信、电力、路灯等管线的探测，优点在于抗干扰能力强，
缺点在于必须有出露的管线才可以使用。
（2）直连法是指将管线仪接线的负极 ( 黑色 ) 接到接地电极，正极 ( 红色 ) 接到管线上。接地电极在情况允许下，尽可能远且垂直于管线走向，正极连
接管线的位置要除锈、刮开保护壳 ( 刮开油漆或者塑料层 )。主要应用于有示踪线的燃气管、有电位测试桩的石油管或燃气管、铸铁输水管等。优点在于
可以较长距离的追踪管线，缺点是抗干扰能力相对较弱。
（3）感应法是指将发射机放置在目标管线上方，接收机沿着管线走向确定其位置和深度大方法。该方法的优点在于适用性强，缺点在于易受到其他信号的干扰。在实际工作中，以上两种方法不适用才使用该方法。
根据管线仪接收机接收信号确定位置和深度的原理不同，可以将常用探测方法分为以下 3 类：极大值法 ( 峰值法 )、极小值法 ( 谷值法 )、极大 / 极小值
法 ( 峰值 / 谷值法 )。将管线中等效电流简化为无限长直导线，根据右手螺旋定则，管线周围会产生“环状”的磁场，如图 2 所示。
在目标管线正上方的磁场水平分量达到*大，垂直分量*小。测量水平分量对应极大值法，测量垂直分量对应极小值法。RD8100 接收机内部有 2个水平线圈 ( 用于极大值法 ) 和 1 个垂直线圈 ( 用于极小值 )，采用极大值法定位时 2 个水平线圈工作，抗干扰能力较强；采用极小值定位时垂直线圈工作，抗干扰能力相对较弱；采用极大 / 极小值法定位时，结合以上 2 种方法，但隔一段距离需要极大值法或极小值法验证。
2.3 管线仪的数值模拟
等效电流在周围激发磁场的基本解析表达式 [2，5] 如式（1）、（2）：Hx= μIh2π（x2+h2）（1）
Hz= μIx2π（x2+h2）（2）
式中：I 为电流大小，μ 为管线周围介质的磁导率，h 为地面到管线的距离，x 为探测点到目标管线地面投影的距离。以某测区实际探测的情况为例，电流为 10 mA，h 埋深为 0.8 m，周围介质为土壤，磁导率为 0.001。地面的信号 Hx 和 Hz 如图 3 所示，可以看到极大值法和极小值法管线磁场信号变化规
律。极大值法在目标管线正上方达到*大，离目标管线越远，Hx 的值越小，趋近于 0，如图 3 中线 1 所示。极小值法在目标管线正上方为 0 值，距离目标管线为 h 时有*大值，远离目标管线 Hz 的值越小，趋近于零。若考虑场的方向，极小值的规律如图 3线 2 所示；只考虑大小，极小值法的规律如图 3 中线3 所示。图 3 中极小值法线 2 和线 3 在场值为正时，曲线完全重合。
若目标管线附近存在其他干扰源，假设距离目标管线 0.5 m 处有其他信号干扰，等效电流为 9 mA，埋深为 0.8 m，可以对干扰情况下的磁场进行模拟，模拟结果如图 4~ 图 7。
如图 4 所示，干扰源和目标管线中的等效电流一致时，叠加信号强度变大，峰值往干扰源方向偏移；如图 5 所示，等效电流方向相反极大值法探测时，叠加信号强度变小，且出现了两个极大值点，不对应目标管线，在实际中会造成探测结果误差。
如图 6 所示，目标管线和干扰管线等效电流方向一致时，叠加信号的极小值偏向干扰管线方向。如图 7 所示，等效电流相反极小值法探测时，叠加场
中有两个极小值点，与目标管线不对应，在实际探测中会产生误差。
2.4 管线仪实际应用
（1）上节讨论的干扰源与目标管线的等效电流大小近似，产生了如图 5 和图 7 中的误差情况。实际探测中可以改变发射机的发射频率，重复探测来确定管线的位置和埋深。
（2）若测区中无明显管线出露，可以采用感应法：一人手持发射机，另一人手持接收机，间距 15~20 m同时移动在测区内感应，确保无遗漏。
（3）管线探测时，可以在有出露或埋深较浅的位置开挖验证，确保探测准确性。
3 探地雷达的工作原理和应用
3.1 探地雷达的探测原理
基于电磁波遇到不同的介质会发生反射和透射的原理，探地雷达发射天线向地下发射电磁波，遇到管线发生反射和透射，反射回地面的信号被接收天线接收，透射信号向下传播，直到信号能量耗尽 [4,6]。
根据反射信号的走时和强度，经过数据处理可以确定地下管线的材质、管径、埋深等信息。
3.2 探地雷达的探测方法
探地雷达根据工作方式可以分为以下 3 类：A-scan、B-scan、C-scan。其中 A-scan 即为连续采集模式，发射天线和接收天线位置固定，采集数据对应地面某一点的地下情况；B-scan 即为沿着剖面线采集的方式，发射天线和接收天线相对位置固定，同时沿着某一设定好的测线进行探测，采集的数据对应测线下剖面的地下情况，常见的“推车”工作方式即为该工作方式；C-scan 即为在地面布置一定间距的测线，在每条测线上进行 B-scan 测量，采集数据对应地下三维情况。
3.3 探地雷达的数值模拟
探地雷达数值模拟的物理方程是Maxwell方程，对该方程进行微分表示，推导出符合探地雷达模拟的二维横磁波（TM）方程 [7]，如式（3）：
∂Ez∂t = 1ε（∂Ey∂x - ∂Hx∂y -σE）z∂Hx∂t =- 1μ（∂Ez∂y+ρHx）∂Hy∂t = 1μ（∂Ez∂x -ρHy）（3）
式中：Ez 对应电场的 z 分量，Hx、Hy 对应磁场的 x 和 y 分量，ε、μ 对应介电常数和磁导率，σ、ρ 对应电导率和电荷密度。使用有限差分法对式（3）进行离散，就可以进行正演模拟。
GprMax 是探地雷达时间域有限差分正演软件，使用该软件需要对以下参数进行定义：模型大小、网格尺寸、时窗长度、介质参数、子波类型，发射天线和接收天线位置及步进距离，介质在模型中的位置等；吸收边界和并行等参数为默认。
设置模型大小为4 m×3 m，网格尺寸为0.005 m，时窗为 60 ns，子波是主频 200 MHz 的雷克子波，天线步进距离为 0.1m，介质参数 [6] 见表 1 所示。
如图 8 所示，模型顶部介质为空气表示地表，两侧和底部是匹配吸收层 (PML)，用于去除有限差分法人工边界的反射，管道为厚度 5 cm 的铸铁管，管
道内部充满水。
如图 9 和图 10 所示，在 35 ns 可以看到铸铁管对应的反射数据，地面回波占据了大部分能量。去除地面回波后，反射曲线更清晰，便于相位、振幅等的研究。再以 PVC 材质的管线为例进行模拟，同时放置
3 个不同深度的管线，对比深度对模拟数据的影响。如图 11 所示。设置模型大小为3 m×2 m，网格尺寸为0.005 m，时窗为 60 ns，子波是主频 200 MHz 的雷克子波，天线步进距离为 0.05 m，介质参数见表 1 所示。
如图 12 和图 13 中所示，A、B、C 分别对应不同埋深�� PVC 管线的反射曲线，随着深度增加信号强度降低，同时出现了多次反射波。
再以水泥管为例进行模拟，设置模型大小为6 m×3 m，网格尺寸为 0.005 m，时窗为 80 ns，子波是主频 200 MHz 的雷克子波，天线步进距离为 0.1 m，
介质参数如表 1 所示。左侧管中为水，右侧管中为空气。如图 14 所示。
如图 15 和 16 中所示，充满水的管对应反射曲线为 B，空管对应曲线为 A。曲线 B 的强度相对较大，同时稍深位置产生了一个较弱的曲线 C。在两个管之间产生了较为复杂的多次反射曲线 D。
3.4 探地雷达的实际应用
以湛江市饮调水工程某测区内的 DN1200 铸铁给 水 管 为 例，BS-M 探 地 雷 达 发 射 天 线 频 率为 200 MHz，相对介电参数为 7，外业采集数据如图
17，可以看到在深度约为 1.4 m 处有明显的反射曲线。经过现场验证，判断为 DN1200 的铸铁水管。
该工程某测区有 DN315 材质 PVC 的燃气管，使用直连法无信号，推测示踪线被破坏。因此使用探地雷达探测，相对介电常数为 10，采集的雷达数据如图 18。可以看到在深度约为 1.0 m 处有一曲线，沿着管线走向前后探测，确定该信号对应燃气管。
该工程另一测区有 DN1000 材质水泥的排水管，可直接开井探测。为了验证探地雷达信号，发射天线频率为 200 MHz，相对介电参数为 10，外业采集数据如图 19。可以看到在深度约为 0.7 m 处有一曲线，对应水泥管，与实际开井量测深度一致。
4 结束语
文章讨论了水利工程中地下管线探测的管线仪法和探地雷达法的探测原理、探测方法、数值模拟、实际应用等。对管线仪探测极大值法和极小值法产生误差的原因进行讨论，采用 GprMax 软件对不同材质、管径、埋深、管内是否含水情况下的管线进行正演模拟，同时对 3 个实地探地雷达数据进行讨论。
给出以下建议：
（1）地下管线埋设情况多样，管线仪法和探地雷达法是获得地下管线信息的重要手段。管线仪法在金属管线探测方面优势较大，探地雷达法则弥补了
非金属管线探测的不足，结合两者可以更好地服务
地下管线探测。
（2）管线仪探测方法的选择取决于探测管线的材质和埋设情况，不同管线需选用不同的探测方法，具体情况具体分析。
（3）探地雷达是一种高效的探测方法，实测数据由于地面回波、场地不平整、噪声等影响，目标管线的反射曲线会被“压制、扭曲”等，不易识别。通过
GprMax 软件正演模拟探地雷达的电磁波传播，可以研究不同地下情况的探地雷达信号，有利于实际数据的识别。
