MATLAB在海底管线阴极保护电场计算中的应用.pdfVIP

维普资讯 MATLAB在海底管线阴极保护电场计算中的应用 胡 舸，向 斌，张胜涛 (重庆大学 化学化工学院，重庆 400044) 摘要 ：海洋钢结构 阴极保护电场的数值模拟仿真和计算是 当前海洋腐蚀 、阴极保护研究和工 程实践中一个重要的研究方向。研究海底管线的腐蚀状态就必须准确测定管线表面的电位 分布，而实际测量 中深海管线表面电位分布难以直接测量。构建 了海水 中钢管线的阴极保 护数学模型，采用以边界元法为基础 的数值计算方法 ，利用 Matlab开发 了均匀电解质 中电 位分布 的计算软件 ，实测 中通过对腐蚀 电场海水区域环境 电场 的测量，就可实现不接触测量 管线表面电位 的 目的。最后通过海洋实测验证 了计算结果的可靠性 ，证 明了基于Matlab的 数值仿真计算引入 阴极保护设计 中是行之有效的。 关键词 ：海底管线 ；Matlab；边界元 ；数值计算 中图分类号 ：TQ015．9 文献标识码 ：A 文章编号：1000—3096(2007)12—0034—04 在海洋环境下传输油气 的金属管线腐蚀远 比陆 位，非常适合边界元计算程序的设计 ，可以大大提高 地上要严重 ，而且一旦发生腐蚀破坏 ，维修极为困难 ， 编程效率。 后果不堪设想 。海底管线腐蚀 电场 电位分布 ，对准确 1 阴极保护 电场计算模型 获取海底管线的腐蚀状态和牺牲 阳极是否处于正常 的工作状态信息都有重要 的意义。但是实际测量 中 1．1 原理 被保护管线表面及腐蚀介质 中电位难 以一一测量 ， 海底管线阴极保护系统中牺牲 阳极与海水和金 导致监测数据不够全面 ，很难真实地反 映管线 的实 属管线构成 回路 ，当海底管线上的保护层出现老化破 际状态 。本研究的 目的就是针对牺牲 阳极保护下的 损或失效时，牺牲阳极通过管线周 围海水向管体表面 海底管线，以数值计算方法为基础，开发海底管线腐 提供保护电流。从而 ，在管线周 围的海泥及海水环境 蚀 电场分布 电位 的计算软件 ，通过对腐蚀 电场海水 中形成非均匀 电场 。通常，海底管线阴极保护电场可 区域 内电位差测量 ，可反演计算管线表面 电位 ，从而 由3个要素来描述 ：(1)管线周 围环境 中的电场状 获得管线腐蚀情况 。 态 ；(2)管线表面上保护 电位分布 ；(3)发 自牺牲 阳 在数值计算领域 中，最初得 以成功应用 的有限 差分法n及有 限元法[2]，为实现阴极保护的现代设 极的保护 电流强度 。 计方法提供了理论上 的可能途径 。但是，由于有限元 根据腐蚀 电场的边界条件 ，应用边界元理论 ，可 法 的固有缺陷，很难在对应于海洋结构物等大型复 计算 出电场 内部和管线表面 电位分布。如果对不 同 杂空间问题以及无限域 问题上得到有效应用 。而后 海水环境 中的极化特性进行系列数值模拟仿真计算 ， 来发展 的边界元法能够精确地对应均匀电解质 中的 就可 以建立一个完备的 “环境电场——管线电位分布 阴极保护问题口]，分块边界元技术能够 圆满解决非 数据库”。因此 ，在海底管线腐蚀检测 中只需测量海 均匀电解质 中的阴极保护 问题[3“]，极大地推动了以 底管线附近环境中的电位 ，测量结果经反演计算就可 数值模拟计算为基础 的现代阴极保护设计方法的发 以获得管线表面电位分布和保护 电流强度 ，从而知道 展 ～ 。 边界元方法在计算过程 中有大量 的矩阵运算 ， 若采用Fortran、C 等其他高级语言开