阴极保护电位分布数值计算的研究进展.pdfVIP

阴极保护电位分布数值计算的研究进展 枉艳霞·1， (I．中国石油人学(华东)山东东营257061； 张国忠1，李蓉蓉2 2．中钉化销售有限公司西南分公司，广东茂名525000) 摘要： 阐述r国内外酬极保护电位分布数值计算的研究概况，包括数学模型的建市，边界条件的选择以及所采用 的数值计算方法．提¨l了该领域存舟的问题及今后的发展方向。 关键词：阴极保护：数值计算；数学模型；边界条件；极化特件 1前言 在阴极保护工程中，电位是监视和控制阴极保护效果的一项重要指标，因此掌握被保护体表面 的电位分布是非常重要的。在传统的阴极保护工程中，大多采用实际测景或经验估计的方法来获得 被保护体表面fI勺电位分布。然而，对于某些被保护结构，如海底管道、海洋平台、深埋的钢桩或钢 管、大罐的罐底等，实地测量技术难度很大，或需要昂贵的费用，对于待上马的新项目，更不可能 事先在现场测试，而经验公式无法对复杂结构表面的电竹分布进行准确的预测，可见传统的方法已 难以满足越来越高的安仝可靠性和经济性的要求。随着电化学和计算机技术的发展，人们尝试采用 数值计算方法来获取被保护体表面的电位和电流分布状况，这已成为最近十多年阴极保护领域中十 分活跃的一个方面。 2阴极保护体系的数学模型 数值计算的实质就是利用数值方法来求解描述所研究问题特性的数学模型，数学模型包括问题 的描述方程和相应的边界条件。近几十年来，随着对阴极保护理论的深入了解和高性能计算机的普 及，阴极保护设计中数学模型的研究取得某些实际应用，在计算阴极保护构件电位分布和阳极参数 等方面显示出模型研究的重要作用114s1。 2．1描述方程 国内外对阴极保护体系电位分布的数值研究中，大多采用拉普拉斯(Laplace)方程作为电位分 布描述方程[H7】，方程形式为： 胛2击=0 (1) 式中由为电位，k为所研究区域内介质的电导率。 要保证拉普拉斯方程的有效性，需要假设所研究区域内导电介质均匀，系统内没有电流的得失， 不存在源点或汇点，系统状态不随时间发生变化即处于稳态。 当所考察区域内有场源存在时，电位分布方程的形式变为泊松(Poisson)方程： kV 2≯=q。 (2) 对于实际的阴极保护体系，由于阴极垢层的形成，以及环境参数如土壤含水量、温度等周期变 化的影响，阴极保护系统的状态并不是一成不变的，而是随着时间发生着变化。由于缺乏对变化规 律的理解，至今尚未出现描述刚极保护体系电位分布动态变化过程的描述方程。Santiago 161采用动态 极化till线来模拟阴极表面的特性，在足够小的时间段内假定阴极保护系统处于稳态，应用Laplace方 程进行求解．求得不同时间段所对应的电位和电流分布情况，从而得到阴极表面电位和电流分布随 时间的变化规律。刘曼，殷正安等人I”1分别采用瞬时极化曲线和阴极钙镁盐沉积达到稳定时的极化 由线，采用边界元方法对Laplace方程求解得到了流动管内的瞬时电位分布和稳定电位分布。 DeGiorgi口]在模拟船板的|；Ij极保护系统时所采用的方法是模拟某一瞬时的状态，作为一独立的边界元 进行求解。在模型中没有将阳极定义为电源，而是作为电位值恒定的边界条件，这样保证了对于稳 态情况采用Laplace方程的有效性，即一个特定的计算机求解结果给出了在一定的时刻，一定的阳极 电流水平下的稳定状态。 2．2边界条件 拉普拉斯或泊松方程的定解取决于求解区域的几何布局和边界条件。从数学上讲，满足一个偏 微分方程的解可以有很多，还必须有一些特定的边界条件来补充和限制，才能得到定解。 2．2．1边界条件的选择 边界条件一般根据实际问题的特点提出，能否选择适当的边界条件直接影响数值计算结果的准 确性。在阴极保护系统中通常遇到的边界条件有以下几类： (1)边界上的电位已知，即妒b，y)=九，例如，无穷远处边界电位为零，阳极可看作电位恒 定的等位面，由恒电位控制的Jl,DIl电流阴极保护系统中的控制点电位恒定等。 (2)边界上流入的电流密度已知，即a烈墨y)，ab，y)=q。。例如恒电流控制的外加电流阴 极保护系统其绝缘面是电流等于零的特例，阳极的输出电流恒定等。 (3)边界上流入的电流密度与电位的函数关系已知·即a≯G，y)lOn(x，y)=厂船)·例如给出边界 上的极化曲线函数。 其中以x，y)为点G，y)处的电位值，九G，y)为在该点处表面的法线方向，g^为恒定的电流值， 九为恒定的电位值，，◇)为极化电流密度函数。 DeGiorgilll在模拟船侧推进器的阴极保护电位分布时，所定义的计算模型区域边界由阳极、阴极 和绝缘面组成。阳极作为固定值的边界条件，具有恒定的电位值；阴极采用极化函数进行定义，极 化函数通过实验测定，绝缘面电流为零。Douglasl5】对罐底的阴极保护电位分布进行模拟时，阴极