阴极保护数值模拟计算边界条件的确定.pdfVIP

文章编号：1000．8241(2011)07．0504．04 阴极保护数值模拟计算边界条件的确定 郝宏娜’李自力’王太源’ 丁延鹏 衣华磊 (1．中国石油大学(华东)油气储运工程系，山东青岛 266555；2．中海油研究总院，北京 100027) 郝宏娜等．阴极保护数值模拟计算边界条件的确定．油气储运，201l，30(7)：504 507 摘要：数值模拟技术是求解阴极保护电位模型的有效方法，国内外流行的阴极保护数学模型包括分 布型模型和时变型模型两种。在数值模拟过程中，边界条件的确定及其处理方式直接影响计算结果 的准确性。为此，对国内外阴极保护数学模型的边界条件进行了分类，包括：边界电位已知的第一类 边界条件 (Dirichlet边界)，电流密度 已知的第二类边界条件 (Neumann边界)，以及边界电位和电流 密度函数关系已知的第三类边界条件。系统阐述 了阳极和阴极边界条件的确定方法以及非线性边 界条件的处理方式，与Newton—Raphson迭代法相比，分段线性拟合法具有更好的灵活性和实用性， 对于实测极化曲线的处理效果非常明显。指出在 阴极边界条件方面，应加强对极化模型的研究。 关键词：阴极保护；数值模拟；边界条件；极化；非线性 中图分类号：TE88 文献标识码：A DOI：CNKI：13—1093／TE1305．001 在阴极保护系统中，掌握金属表面的阴极保护电 情况下电位分布的控制方程为Laplace方程： 位分布对阴极保护系统的设计、优化、日常管理和维护 V (一 V )= 0 具有重要的意义。近年来，随着电化学理论和计算机 式中：V为Laplace算子， 为介质的电导率，为阴极 技术 的发展，数值模拟计算成为阴极保护电位研究领 保护电位。 域的新热点，该方法对于大型储罐、海底管道等 电位不 当计算范围内的介质均匀， 为常数时，方程简化 易测量设施的阴极保护系统的设计和评价具有重要意 为：V 一0。 义。数值模拟计算的本质是利用数值方法求解计算描 述实际问题的数学模型，其中：数学模型由控制方程和 2边界条件 边界条件组成。在数值模拟过程中，边界条件的确定 及其处理方式直接影响计算结果的准确性。 从数学的角度，满足Laplace方程的解有很多，要 得到定解，必须对其研究区域进行 限制并确定相应的 1控制方程 边界条件，实际问题模型数值计算结果的准确性往往 取决于边界条件的选取是否恰当。 目前，国内外流行的阴极保护数学模型包括分布 2．1边界条件的分类 型模型和时变型模型两种。分布型模型用以研究空间 要计算 1个体系的阴极保护 电位Q，首先需要 上的阴极保护 电位和电流密度的关系，假设环境介质 确定体系的范围 ，通常阴极保护系统由阳极、阴极禾1 均匀，对于不均匀的体系则分解成若干均匀体系进行 绝缘面组成。与之对应的边界条件主要有3类：① 研究，同时，假设整个体系达到稳态，即不考虑阴极极 边界 电位 已知，这属于数理方程的第 一类边界条件 化行为随时间的变化，采用静态场进行描述。时变型 (Dirichlet边界)，例如：阳极工作面的电位 已知，无穷 模型则考虑 了阴极极化行为随时间的变化，可以计算 远处电位为0等；②给出边界电位的导数，即电流密度 阴极保护 电位随时间变化的分布规律，比分布型模型 已知，这属于数理方程的第■类边界条件 (Neumann 研究难度大，在此不做赘述。 边界)，例如：在稳态的阴极保护体系中