瞬变电磁超前探水数值模拟与应用.pdfVIP

· 80· 路基工程 2009年第4期 (总第 145期) 瞬变电磁超前探水数值模拟与应用 谭代明 漆泰岳 (西南交通大学土木工程学院 四川成都 610031) 摘 要 在数值模拟 中通过改变隧道施工掌子面前方低阻体的位置 ，获得瞬变磁场在相应条件下 的衰变值，将实测值与数值模拟结果进行曲线拟合，通过比较拟合结果参数中的衰减因子值，找到与 实测曲线最接近的模拟结果曲线。工程应用实例表明，数值模拟计算结果与反演软件结果的结合，可 以去除反演结果中的伪解。 关键词 隧道 瞬变电磁法 数值模拟 反演 探水 目前在我国应用于隧道超前地质预报的物探方法 的值 。 中以微震波反射法 (aSP)和探地雷达 (GPR)较为 (4)式即为电流源激磁的瞬变磁场有限元方程。 常见。TSP对岩石变化、断层带和不良地质夹层等的 通过求解该方程，就可得出不同时刻的矢量磁位A的 探测预报有相当高的准确性，对不良地质体中是否存 值 ，再由磁感应强度 日与A的旋度关系式B= ×A， 在突泥突水等问题尚有较大差距…。GPR发射高频电 就可得到瞬变磁场的变化规律。 磁波，对浅层 目标的探测比较有效 。 瞬变电磁法勘探是近年来国内外发展得较快 、探 2 计算模型及参数说明 水效果较好的一种物探方法，但对于它的理论研究却 2．1 计算模型 滞后于工程的需要 ，应用于实际工程中的探测数据的 实际的瞬变电磁法问题都是三维的，对三维的计 反演仍基于一维层状理论，这会造成结果的多解性。 算需要时间太长，而采用小模型时得到的结果误差较 本文将反演结果与数值模拟的正演结果相结合 ，有针 大。在岩溶隧道中，正对掌子面前方的充水溶洞最易 对性地调整数值模拟计算中的地电参数，能有效地消 引起突水、突泥灾害。应用于隧道瞬变电磁法激励源 除反演结果中的多解性。 一 般采用边长为2～3m的多匝发射线圈，线圈平面 平行掌子面，形成纵向磁偶极源 (环向电流)，通过 1 二维有限元理论 磁探头或小半径接收线圈得到空间某点处瞬变磁场。 瞬变电磁法脉冲频谱能量集中于 r,／～n×100Hz 将计算模型简化为二维轴对称问题，从而能在个人计 的范围内，在忽略位移电流条件下，可以得到二维轴 算机上实现。这种简化显然比一维层状模型更符合隧 对称扩散方程 道探测实际。 ( )+0I1 )=一+了or 在ANSYS中将场域分为有限区和无限区两部份， 有限区采用 PLANE53电磁单元 ，并取AZ 自由度，无 (1) 限区采用 INF110远场单元。 式中 为激发源电流密度；x／为磁导率；or为电 2．2 参数说明 导率。 工程探测中激励 电流的接通与关断在 ANSYS中 取权函数等于形状函数 ，利用伽辽金加权余量 通过阶跃加载的方式来实现。计算模型及参数如图 法在求解场域 范围内对 (1)式进行加权积分得 1，其中P。表示空气 ，P表示高阻围岩，P表示低阻 ＼(’ Adn+＼J 一＼专 =0 岩溶水体。计算模型如图1所示 (图中未表示远场 区)。