堤坝探测中电阻率成像3D有限单元法模拟.pdfVIP

1 ·802· 中国地球物理201 (21)工程地球物理 堤坝探测中的电阻率成像3D有限单元法模拟 3DFEMSimulationof onEmbankment ResistivityImaging Detecting 王锐冷元宝 周扬 Rui Zhou Wang LengYuangbaoYang 1．黄河水利科学研究院 郑州450003；2．黄河流域水资源保护局郑州450003 在堤坝安全管理和堤坝隐患探测工程中广泛的采用了电阻率检测和成像技术。尤其是在近几年，电阻 率成像在堤坝内部侵蚀探测中的作用越来越大。这是因为这个方法效率高，且结果包含了大量的安全信息。 经过多年的研究和发展，硬件设施已经达到成熟。电压和电流的精确测量已经接近了极限，这虽然是获得 好的成像结果的主要因素之一。但是，关键还是反演的理论和技术问题。 电阻率成像和其他的电阻率方法广泛的应用于探测堤坝中的潜在问题。堤坝有着复杂的形状且有多种材料 建成。有限单元法(FEM)是解决像堤坝这类的复杂目标体的正演问题的有效的方法。本文介绍了怎样用ANSYS 建立3D堤坝模型以及模拟内部侵蚀的探测。我们模拟了通过温纳阵列探测堤坝中一个渗漏道的过程，并列出 了结果和APDL源程序。ANSYS和APDL是一个解决类似于堤坝等复杂目标体的止演问题的有效方法。 和快速的工具。APDL能够自动的进行多阶模拟，或者有更少的介入和干扰。传统的转化和反演的理论是 基于水平、半无限、各向同性的导电体模型。这样的模型非常简单，降低了解释的精度。本文研发的有限 单元法软件类似于ANSYS的专业软件，使用该软件，我们可以建立3D堤坝模型进行转换研究。我们不 仅可以比较不同类型阵列的有效性，而且还能够讨论堤坝形状的影响。这正是·般的2D有限单元和有限 差分计算分析方法所无法完成的。 ANSYS有限元网格划分直接影响着后续数值计算分析结果的精确性。而在计算资源相对有限的条件 下，网格划分的优化能够在尽量保证计算精度的前提下，减少计算时间，提高工作效率。网格划分涉及单 元的形状及其拓扑类型、单元类型、网格生成器的选择、网格的密度、单元的编号以及几何体素。在计算 数据变化梯度较大的部位(如电流集中处)，为了较好地反映数据变化规律，需要采用比较密集的网格。而 在计算数据变化梯度较小的部位，为减小模型规模，则应划分相对稀疏的网格。在此次计算中，电极排列 附近是电流的加载点和电压的取值点，需要加密划分网格，较远的区域网格划分较为稀疏。 我们设计了包括不同隐患类型的堤防模型(球体、圆柱体等)进行有限单元法模拟计算。基本堤防计算模 型横截面为梯形，上底宽5m，高5m，斜边坡度为1：3，长度40m，为模拟堤防接地条件，将堤防模型设置在 长40m，宽50m，高5m的台上。电极排列在堤顶中线居中位置布设，平行于堤防走向。在基本模型的基础上， 通过增加不同的异常体模拟堤防常见的隐患，包括层状异常体、球状异常体、管状异常体、直立板状异常体。 几何模型建立后，首先进行网格剖分，为优化网格划分参数，开展了不同网格划分条件F有限元计算 的误差分析工作，以在保证计算精度的前提下，尽可能地减少网格的节点数。我们通过10组不同的网格 剖分进行计算，并与理论值进行了对比。对比结果表明，10组不同剖分方案中，相对误差最大的为2．42％， 最小的为0．02％。最终选择了相对误差最小的划分方案，该方案以测线为中心，从内到外分3层进行划分， 内层与中层为长方体区域，外层为堤防与地台表面。在内层区域沿测线方向增加了分割线，以更好的控制 节点位置。划分采用4面体单元solid23l，对堤防中的异常体及其附近区域进行加密划分，其单元尺寸控 制在与电极附近单元尺寸相当的水平。 为模拟不同电极距、不同装置对模拟隐患的探测效果，按照电法装置的实际工作情况，在划分好网格 的有限元模型上，依次在模拟电极的位置加载电流，在测量电极的位置提取电压数据。通过计算得到视电 阻率值后，结合装置系数等参数，组成电测剖面，绘制视电阻率成果图，并进行反演。从成果町以看出， 由于堤防边界的