基于滤波算法的隧道全时域激发极化超前探测正演方法  

Forward method for full time domain induced polarization advance detection of tunnel based on filtering algorithm

作　　者：李志强 褚国安 聂利超[1] 裴文兵 郭一凡 周维 胡军 闫冰 邓朝阳 LI Zhiqiang;CHU Guoan;NIE Lichao;PEI Wenbing;GUO Yifan;ZHOU Wei;HU Jun;YAN Bing;DENG Zhaoyang(Geotechnical&Structural Engineering Research Center of Shandong University,Jinan 250061,China;Yunnan Water Diversion Phase II Project Co.Ltd.,Kunming 650205,China)

机构地区：[1]山东大学岩土与结构工程研究中心,山东济南250061 [2]云南省滇中引水二期工程有限公司,云南昆明650205

出　　处：《中南大学学报（自然科学版）》2024年第2期638-648,共11页Journal of Central South University:Science and Technology

基　　金：国家重点研发专项课题(2021YFC3100805);国家自然科学基金资助项目(52079076);山东省自然科学基金资助项目(ZR2022YQ46,ZR2022QD014);云南省科技厅重大科技专项计划项目(202002AF080003)。

摘　　要：为进一步提高隧道全时域激发极化多参数超前探测的效果,采用滤波算法进行计算时间域等效电阻率,采用有限单元法进行计算,提出一种基于滤波算法的隧道全时域激发极化超前探测正演方法。该方法针对时间域Cole-Cole表达式计算收敛性差的问题,引入滤波算法进行时间域激发极化等效电阻率的计算,提高算法的收敛性;针对正演过程中大型矩阵求解影响计算效率的问题,引入OpenMP并行算法,实现隧道全时域激发极化超前探测正演的加速运算。最后,开展数值算例和物理模型试验。研究结果表明:该正演方法对异常体有敏感响应,且在加快正演速度具有显著优势。通过基于滤波算法的隧道全时域激发极化超前探测获得的观测数据与实际探测数据较为相符,这为基于全时域激发极化隧道超前探测提供了一种可行的正演算法。In order to meet the requirements of high detection efficiency in actual tunnel engineering,the filtering algorithm was used to calculate the equivalent resistivity in the time domain,and the finite element method was used for calculation.A forward modeling method of full-time induced polarization advance tunnel detection based on the filtering algorithm was proposed.To address the problem of poor convergence of Cole-Cole expression in time domain,the filtering algorithm was introduced to calculate the equivalent resistivity in time domain induced polarization,improving the convergence of algorithm.To solve the problem that the large-scale matrix solution affects the computational efficiency in the forward modeling process,the OpenMP parallel algorithm was introduced to realize the accelerated operation of the full-time domain induced polarization advanced detection forward modeling of the tunnel.Finally,numerical examples and physical model tests were carried out.The results show that the forward modeling method has a sensitive response to abnormal bodies.The full-time domain induced polarization forward acceleration algorithm has significant advantages in accelerating forward modeling speed.The forward data obtained by the full-time domain induced polarization forward detection method based on the filtering algorithm is in good agreement with the actual detection data,which lays a foundation for the application of the full-time domain induced polarization inversion method in the actual project of tunnel geological prediction.

关 键 词：隧道超前地质预报 全时域激发极化 滤波算法 快速三维正演模拟 并行算法 

分 类 号：P631[天文地球—地质矿产勘探]