基于数值模拟的海底管道外表面阴极保护系统优化设计.docxVIP

基于数值模拟的海底管道外表面阴极保护系统优化设计 0 海底管道外表面阴极保护的组成 原油管道是海洋石油资源开发和生产中最重要的管道之一。该项目的投资成本高，服务时间长。海洋环境腐蚀条件苛刻, 海底管道外表面在下水安装后将自然地呈现腐蚀倾向。资料显示, 在海底管道结构损伤案例中, 由腐蚀问题导致的结构失效约占35%, 可见对海底管道的腐蚀防护至关重要。 海底管道的腐蚀损伤源于两个方面:管道内腐蚀和管道外腐蚀。在实际工程中, 作为海底管道外表面的防腐措施, 多采用防腐涂层结合牺牲阳极阴极保护的复合防腐技术。海底管道外表面阴极保护系统利用电化学手段, 由手镯状阳极向管道外表面提供充分的保护电流, 以使管道外表面上的电位充分负移, 从而达到防腐的目的。海底管道外表面全寿命期的安全性和可靠性对阴极保护系统的设计水平提出了越来越严格的要求。国内的阴极保护系统设计采用的是传统的基于规范的设计方法, 无法实现对海底管道外表面全寿命期防腐系统状态演变的动态描述。目前, 数值模拟计算技术在阴极保护系统设计中的作用已越来越重要。 针对上述问题, 本文将设计可以描述实际工程中海底管道管型的样本, 开展大规模的数值模拟计算, 对输入输出参数进行相关性分析;搭建人工神经网络非线性映射模型, 利用L-M算法和贝叶斯正则化算法进行训练, 获取高精度的神经网络映射关系;输出参数可构成精细化的设计方案, 实现海底管道外表面阴极保护系统的优化设计。 1 海底管道的阴极保护 在海洋钢结构牺牲阳极阴极保护系统设计中, 应用最广泛的规范是阴极保护设计 (DNV-RP-B401) 。针对海底管道外表面阴极保护系统设计, 主要的设计规范包括海底管道牺牲阳极阴极保护系统 (DNV-RP-F103) 和管线传输系统阴极保护第2部分:海底管道 (ISO-15589-2) 。由于B401规范的设计结果偏于保守且应用时间较久, 因而在新建海底管道外表面阴极保护系统设计工程项目中仍有应用。这两部规范的设计流程类似, 均是根据结构湿表面积、初始/平均/最终电流密度需求和牺牲阳极发出电流等输入参数来计算求得牺牲阳极的总块数或单块镯状牺牲阳极的保护长度。对于由普通碳钢、低合金钢构成的海底管道外表面, 规范规定的保护电位数值为-0.80 V (相对于Ag/Ag Cl/海水参比电极) 。对于埋覆状态下输送常温流质的海底管道, 与-0.80 V (相对于Ag/Ag Cl/海水参比电极) 对应的平均保护电流密度取值为0.02 A/m2。 2 计算海底管道外表面尹维护系统的数值模拟计算 2.1 海底管道失稳耦合控制方程 近年来, 计算机分析技术开始逐步应用于海洋工程结构防腐系统设计领域, 显著提高了腐蚀控制水平。本文研究的海底管道外表面阴极保护系统优化设计技术的核心为基于边界元法的数值分析方法, 其通过应用分块边界元技术来圆满解决海底管道所处的海水、海泥双重介质环境下阴极保护系统的数值分析问题。该问题的数值模型定义如下: 式中:q代表电流密度;?为电位;ρ为介质电阻率;下标m, w分别表示处于海泥电解质区域和海水电解质区域。式 (1) 的前2个方程描述的是达到稳定状态时海底管道外表面阴极保护系统伴生电场的域内控制方程;后4个方程代表全部边界条件的选取S=Sc+Sw+Sm+Swm, 分别为海底管道外表面和牺牲阳极表面Sc的电位和电流密度关系满足材料极化曲线。位于足够远处, 海水电解质区域边界Sw与海泥电解质区域边界Sm的电流密度为零;在海水电解质区域Ωw与海泥电解质区域Ωm的交界处Swm, 电位和电流密度数值分别对应相等。 2.2 敏感性分析结果 为了获得满足实际工程要求的优化设计样本数据库, 一方面需要设计具有代表性、完备性特点的样本方案, 实施大规模数值模拟计算, 另一方面, 对节点数为N的模型, 在求解方程组时, 边界元数值模拟计算需要形成维度为O (n2) 的双精度型矩阵, 这对计算机的硬件及求解时间提出了非常高的要求。因此, 在实施大规模样本计算前, 需要进行各输入参数对海底管道外表面阴极保护电场参数的敏感性评估, 合理设计参数点的数量及间隔, 在保证样本完整性的同时, 最大限度地降低计算时间成本。 本文设计了如表1所示的小规模样本, 基于数值分析结果, 对阳极长度Lanode和阳极厚度Tanode进行了敏感性评估。敏感性分析的主要评价指标为牺牲阳极发出电流I。 根据表1中的样本设计情况, 阳极长度、阳极厚度与牺牲阳极发出电流之间的相关性曲线如图1所示。 通过上述数值分析可以得出结论, 阳极长度和阳极厚度与海底管道外表面阴极保护伴生电场有一定的关联性。因此, 在大规模数值模拟计算方案中, 需要体现这两种因素的存在。但由于这种关联性并不强烈, 因而可以选择较大样本间距的参数点作为设计样本。 2