一种冲蚀-二氧化碳腐蚀耦合作用下的天然气注采井管柱腐蚀速率的确定方法
【专利摘要】本发明属于油田设备技术领域，具体提供了一种冲蚀?二氧化碳腐蚀耦合作用下的天然气注采井管柱腐蚀速率的确定方法，包括步骤：1)建立管柱内组分传热、流动与扩散模型；2)建立管柱壁面二氧化碳腐蚀模型，并确定二氧化碳腐蚀过程中产物层对腐蚀速率的影响；3)建立介质冲刷对管壁基材、腐蚀层的磨损模型。4)将以上三模型进行耦合，确定天然气注采井管柱的腐蚀速率。本发明通过天然气注采井管柱中的冲蚀?电化学腐蚀耦合作用机理，建立气井管柱腐蚀速率预测模型，能够有效预测天然气注采井管柱的腐蚀速率，为携二氧化碳及液态水的天然气注采井管柱安全性提供技术参考。
【专利说明】
-种冲蚀-二氧化碳腐蚀輔合作用下的天然气注采井管柱腐 蚀速率的确定方法
技术领域
[0001 ]本发明属于油田设备技术领域，设及一种天然气井垂直管柱安全服役评估方法， 尤其是一种冲蚀-二氧化碳腐蚀禪合作用下的天然气注采井管柱腐蚀速率的确定方法。
【背景技术】
[0002] 随着我国对清洁能源需求的增长，天然气的开采量不断增大。然而在天然气的开 采过程中，绝大部分区块含有C〇2、也0等腐蚀介质，并进入集输和储存系统中，运不仅影响天 然气采集输管材和处理设备的使用寿命、对油气工业造成较大的经济损失，也对环境造成 了很大的污染。随着天然气的逐年开发，储层深度亦不断加深，携带的C〇2、也0等腐蚀介质含 量亦随之变化，加之集输和储存过程中的高溫、高压、流速和流态变化的相互作用，管材在 流场作用下的腐蚀问题也愈加突出。因此，研究流场作用下的管内腐蚀预测问题，对于实验 指导在役天然气注采井管柱系统的检测工作量，保障管管柱统的安全运行，具有重大的现 实意义。
[0003] 天然气注采井管柱的内壁腐蚀不仅预测、检测困难，而且会严重影响输送效率、设 备安全和管柱的可靠性。管柱的内腐蚀不仅会造成管柱的维护与更换周期变短，而且腐蚀 产物层会影响管柱的传热效率降低输送介质的流速，使管柱介质的输送动力能源消耗增 大。多相流的复杂流动会加速材料在应力和腐蚀环境协同作用下发生的开裂及断裂失效现 象。天然气注采井管柱中的多相流输送特点是气相为连续相，液相为非连续相，兼有冲刷 腐蚀及二氧化碳腐蚀的特点。本发明是首次将天然气的冲刷腐蚀及二氧化碳腐蚀结合起 来，结合天然气输送过程中的流场数据，建立注采井管柱的腐蚀速率预测模型。
[0004] 鉴于用物理模拟试验，很难模拟深层地层下的天然气井垂直管柱中高溫高压的复 杂流场W及腐蚀过程，因此，用计算机数值模拟技术为主要手段是今后天然气井垂直管柱 腐蚀机理与腐蚀速率预测研究发展方向。
[0005] 现有的冲蚀模型中大多集中在气固/液固两相流冲蚀研究中，W数学的方式描述 金属材料的冲蚀率随祀材属性、颗粒属性、环境因素的变化规律，从颗粒冲击祀材的碰撞过 程入手，通过求解颗粒运动方程和碰撞的能量方程，建立了基于微切削和变形磨损的塑性 材料冲蚀模型，对于电化学腐蚀过程的讨论很少，不适用于不含砂的携液酸性天然气井管 柱环境。
【发明内容】
[0006] 本发明的目的在于提供一种冲蚀-电化学腐蚀禪合作用下的天然气注采井管柱腐 蚀速率的确定方法，通过天然气注采井管柱中的冲蚀-电化学腐蚀禪合作用机理，建立气井 管柱腐蚀速率预测模型，能够有效预测天然气注采井管柱的腐蚀速率。
[0007] 本发明的技术方案是:一种冲蚀-二氧化碳腐蚀禪合作用下的天然气注采井管柱 腐蚀速率的确定方法，包括如下步骤：
[0008] 步骤一:根据流体力学理论与天然气注采设及的运动状态，建立各组分传热、流动 与扩散数学模型，计算管壁处压力、溫度、速度梯度、各组分浓度分布W及剪切力；
[0009] 步骤二:根据二氧化碳腐蚀理论，在对气井管柱腐蚀数据的详细分析的基础上，建 立二氧化碳腐蚀模型，确定二氧化碳腐蚀过程中腐蚀产物层对于该环境下电化学腐蚀速率 的影响，引入腐蚀层厚度对腐蚀速率的影响因子；
[0010] 步骤根据冲刷磨损理论，对天然气井注采数据的详细分析基础上，建立注采过 程的冲蚀模型，W流体能量损失数据为基础，采用能量磨损理论建立冲蚀模型，并确定冲刷 过程对腐蚀产物层及管壁基体冲刷磨损速率系数；
[0011] 步骤四：将W上步骤一至步骤=所述各模型进行禪合计算，计算出管内的溫度、压 力、入口流速对腐蚀速率的影响，确定天然气气井管柱的安全服役参数。
[0012] 上述步骤一中所描述流动、传热模型采用流体力学N-S方程计算;液滴在气流中的 运动引入Wu曳力公式及牛顿第二定律计算;物质在气流中的扩散采用对流扩散模型计算； 通过上述模型的共同计算，可W得到管壁处压力、溫度、速度梯度、各组分浓度分布W及剪 切力。
[0013] 上述步骤二中所描述的二氧化碳腐蚀模型如下：
[0014]
[001引式中：a，b，d腐蚀速率的修正参数，T是溫度，与0:是C02分压，抑是介质值，Ccof是腐 蚀速率修正系数;V。是二氧化碳腐蚀速率;然后引入腐蚀产物层对腐蚀速率的影响因子来 对腐蚀速率进行修正:Vc ' = [^d (1+C) +a ] Vc，式中：a，0是产物层对腐蚀速率的修正因子，hd 是腐蚀产物层厚度，C是腐蚀速率修正系数，V。'是修正后的腐蚀速率。
[0016] 上述二氧化碳腐蚀模型其中的a，b，d是根据实际腐蚀速率代入模型确定的参数； a, 0是通过建立实际腐蚀速率W及腐蚀产物层厚度关系得到。
[0017] 上述步骤=中所述的根据冲刷磨损理论所建立的冲蚀模型如下：
[001 引
[0019]其中曲是磨损速率;a是磨损速率系数;Fw是剪切力，y是介质流体动力粘度。
[0020] 上述冲蚀模型其中的a是根据实际冲蚀速率代入模型确定的参数。
[0021] 上述步骤四中所述禪合计算方法为:通过流动、传热W及传质模型计算出管壁处 压力、溫度、速度梯度各组分浓度分布W及剪切力；将上述参数导入二氧化碳腐蚀模型，计 算出腐蚀速率W及腐蚀产物的生成速率;通过冲蚀模型计算出在流体的作用下腐蚀产物W 及管壁的冲蚀速率;将冲蚀模型计算结果再重新导入二氧化碳腐蚀模型的腐蚀产物层修正 项，算得最终的腐蚀速率。
[0022] 通过导入不同的环境参数，算得对应的管柱腐蚀速率，选择安全的腐蚀速率，确定 步骤四中所述天然气井管柱的安全服役参数。
[0023] 上述步骤一中所述的传热、流动与扩散数学模型的求解采用fluent软件计算；而 二氧化碳腐蚀、冲刷磨损模型采用fluent软件的UDF用户自定义函数进行禪合计算。
[0024] 使用Fluent软件，将天然气井管柱模型导入后将电化学腐蚀-冲刷腐蚀禪合机理 模型写入UDF，将实际环境中的压力、组分浓度、溫度W及流速参数输入软件进行双精度计 算即可得到计算结果。
[0025] 本发明的有益效果:本发明的模型及方法集中讨论了天然气注采井垂直管柱中， W气相为连续相，液相为离散相的天然气传输过程中冲蚀-电化学腐蚀对管壁的禪合作用。 此过程介质中的含液体积分数低于1%，流速低于25m/s。并使用计算机数值模拟技术进行 计算，计算结果与实际腐蚀结果较吻合，能够为携二氧化碳及液态水的天然气井管柱安全 性提供技术参考。
[0026] W下将对本发明做进一步详细说明。
【具体实施方式】
[0027] 实施例1: 一种冲蚀-二氧化碳腐蚀禪合作用下的天然气注采井管柱腐蚀速率的确 定方法，首先基于流体力学理论与天然气注采设及的运动状态，采用理论分析方法建立流 场数学模型;然后根据二氧化碳腐蚀理论，对天然气井腐蚀数据的详细分析的基础上，建立 二氧化碳腐蚀模型；之后根据二氧化碳腐蚀理论，对天然气井腐蚀数据的详细分析的基础 上，建立二氧化碳腐蚀模型；最后将上述模型写入用户自定义函数导入商业软件进行仿真 计算。本发明方法包括如下步骤：
[0028] 步骤一:根据流体力学理论与天然气注采设及的运动状态，采用理论分析方法建 立流场数学模型，流场数学模型包括传热、流动、传质模型。采用建模软件建立模拟区域的 管柱模型，根据流场计算每个单元的节点坐标数据和流场计算的块中屯、网格数据的集合对 应关系，将建立的管柱网格模型转化为流场计算所需的有限元数值模型，计算管壁处压力、 溫度、速度梯度、各组分浓度分布W及剪切力等参数。其中流动、传热模型采用流体力学N-S 方程计算;液滴在气流中的运动引入Wu曳力公式及牛顿第二定律计算;物质在气流中的扩 散采用对流扩散模型计算;通过上述模型的共同计算，可W得到管壁处压力、溫度、速度梯 度、各组分浓度分布W及剪切力。本步骤中的传热、流动与扩散数学模型的求解采用fluent 软件计算;而二氧化碳腐蚀、冲刷磨损模型采用fluent软件的UDF用户自定义函数进行禪合 计算。使用Fluent软件，将天然气井管柱模型导入后将电化学腐蚀-冲刷腐蚀禪合机理模型 写入UDF，将实际环境中的压力、组分浓度、溫度W及流速参数输入软件进行双精度计算即 可得到计算结果。
[0029] 步骤二:根据二氧化碳腐蚀理论，在对天然气井管柱腐蚀数据的详细分析的基础 上，建立二氧化碳腐蚀模型，确定二氧化碳腐蚀过程中腐蚀产物层对于该环境下电化学腐 蚀速率的影响，引入腐蚀层厚度对腐蚀速率的影响因子。
[0030] 其中二氧化碳腐蚀模型如下：
[0031]
[0032] 式中：a，b，d腐蚀速率的修正参数，T是溫度，巧〇,是0)2分压，抑是介质值，Ccof是腐 蚀速率修正系数;V。是二氧化碳腐蚀速率;然后引入腐蚀产物层对腐蚀速率的影响因子来 对腐蚀速率进行修正:Vc ' = [^d (1+C) +a ] Vc，式中：a，0是产物层对腐蚀速率的修正因子，hd 是腐蚀产物层厚度，C是腐蚀速率修正系数，V。'是修正后的腐蚀速率。二氧化碳腐蚀模型中 的a，b，d是根据实际腐蚀速率代入模型确定的参数;a，e是通过建立实际腐蚀速率W及腐蚀 产物层厚度关系得到。
[0033] 步骤=:根据冲刷腐蚀理论，对天然气井注采数据的详细分析基础上，建立注采过 程的冲蚀模型，W流体能量损失数据为基础，采用能量磨损理论建立冲蚀模型，并确定冲刷 过程对腐蚀产物层及管壁基体冲刷磨损速率系数;根据冲刷磨损理论所建立的冲蚀模型如 下：
[0034]
[0035] 其中出是磨损速率;a是磨损速率系数;Fw是剪切力，y是介质流体动力粘度。其中冲 蚀模型其中的a是根据实际冲蚀速率代入模型确定的参数。
[0036] 步骤四：将W上步骤一至步骤=所述各模型进行禪合计算，计算出管内的溫度、 压力、入口流速对腐蚀速率的影响，确定天然气气井管柱的安全服役参数。具体为:将步骤 一建立的管柱模型导入fluent,在fluent软件中设置好流场边界和初始条件，根据步骤一 中的流场模型、边界类型及初始条件在fluent中选择对应的流场模拟计算模型;并结合注 采井的实际情况，输入井的注采数据，在fluent中设置好模型的边界及初始条件，计算模拟 区域的溫度、压力、流速变化;将步骤二与步骤=中的二氧化碳腐蚀模型与冲刷磨损模型进 行禪合，写入UDF。将步骤一中的流场模型导入fluent进行计算，可W得到管柱内流场中的 溫度、压力、流速、二氧化碳分压、携液量等瞬时状态数据，将UDF导入进行计算。最终计算出 管内的溫度、压力、入口流速对腐蚀速率的影响。本步骤中所述禪合计算方法为:通过流动、 传热W及传质模型计算出管壁处压力、溫度、速度梯度各组分浓度分布W及剪切力;将上述 参数导入二氧化碳腐蚀模型，计算出腐蚀速率W及腐蚀产物的生成速率;通过冲蚀模型计 算出在流体的作用下腐蚀产物W及管壁的冲蚀速率;将冲蚀模型计算结果再重新导入二氧 化碳腐蚀模型的腐蚀产物层修正项，算得最终的腐蚀速率。通过导入不同的环境参数，算得 对应的管柱腐蚀速率，选择安全的腐蚀速率，确定本步骤中所述天然气井管柱的安全服役 参数。
[0037] 本发明的具体计算步骤如下：
[0038] 根据流体力学理论与天然气注采设及的运动状态，采用理论分析方法建立流场数 学模型，流场数学模型包括传热、流动、传质模型，组分扩散模型，和流场边界条件。
[0039] (1)建立质量守恒方程(连续性方程）：
[0040]
[0041] 式中，P是流体密度，V是流体速度。[0042] (2)建立动量守恒方程：
[0043] X方向
[0044] y方向
[0045] Z方向：
[0046] 式中，y、u、w是速度分量，P是流体各向同性压强，T是与流体粘性有关的剪切力，f 是体积力。
[0047] (3)建立能量守恒方程：
[004引
[0049] 式中，e是热力学能，q是对系统传输的热流，k是热传导系数。
[0050] (4)建立组分传输扩散模型：
[0051] 通过第巧巾物质的对流扩散方程预估每种物质的质量分数，Yi。其守恒方程如下：
[0化2]
[0053] 其中，P是密度，V是扩散流动速度，Ji是混合物中第i种物质的扩散通量，Ri是化学 反应的净产生速率，Si为离散相及用户定义的源项导致的额外产生速率。i = l，2,3分别代 表邸4、C〇2、出0立种物质，需要求解两个运种形式的方程。
[0054] 针对管内端流流动，W如下形式计算质量扩散：
[0化5]
[0056] 其中，Di,m是扩散系数，Sct是端流施密特数，化是端流粘性力。
[0057] 根据二氧化碳腐蚀理论，对天然气井腐蚀数据的详细分析的基础上，建立二氧化 碳腐蚀模型;结合冲刷腐蚀模型，建立天然气注采井注采过程的冲蚀-电化学腐蚀禪合数学 模型。二氧化碳腐蚀模型表达式为：
[0化引
[0化9] 式中：曰,b, d都是根据实际腐蚀速率确走的参数,T是溫度，巧..是C〇2分压,pH是介 质值，Cwf是腐蚀速率修正系数。V。是二氧化碳腐蚀速率。之后引入腐蚀产物层对腐蚀速率 的影响因子来对腐蚀速率进行修正:￥。' =[护4(1+(3)+0]￥。，式中：〇，0是产物层对腐蚀速率 的修正因子，hd是腐蚀产物层厚度，C是腐蚀速率修正系数，V。'是修正后的腐蚀速率。
[0060] (5)根据冲刷腐蚀理论，对天然气井注采数据的详细分析基础上，建立注采过程的 冲蚀模型，W流体能量损失数据为基础，采用能量磨损理论建立冲蚀模型。该模型表达式 为：
[0061]
[0062] 其中出是磨损速率;a是磨损速率系数，由实际实验数据推出;Fw是剪切力，y是介质 流体动力粘度。
[0063] (6)采用建模软件Gambit建立模拟区域的管柱模型，根据流场计算每个单元的节 点坐标数据和流场计算的块中屯、网格数据的集合对应关系，将Gambit建立的管柱网格模型 转化为流场计算所需的有限元数值模型；
[0064] 将上述模型写入UDF，并将(6)中的建立的管柱模型导入fluent,在fluent软件中 设置好流场边界和初始条件，根据步骤一中的流场模型、边界类型及初始条件在fluent中 选择对应的流场模拟计算模型;并结合注采井的实际情况，输入井的注采数据，在fluent中 设置好模型的边界及初始条件，计算模拟区域的溫度、压力、流速变化;将UDF导入fluent， 计算出注采井管壁腐蚀速率。
[0065] 本发明提供的模型及方法集中讨论了天然气注采井垂直管柱中，W气相为连续 相，液相为离散相的天然气传输过程中冲蚀-电化学腐蚀对管壁的禪合作用。此过程介质中 的含液体积分数低于1%，流速低于25m/s。并使用计算机数值模拟技术进行计算，计算结果 与实际腐蚀结果较吻合，能够为携二氧化碳及液态水的天然气井管柱安全性提供技术参 考。
[0066] 本实施方式中没有详细叙述的部分或计算方法属本行业的公知的常用手段，运里 不一一叙述。W上例举仅仅是对本发明的举例说明，并不构成对本发明的保护范围的限制， 凡是与本发明相同或相似的设计均属于本发明的保护范围之内。
【主权项】
1. 一种冲蚀-二氧化碳腐蚀耦合作用下的天然气注采井管柱腐蚀速率的确定方法，其 特征在于，包括如下步骤： 步骤一:根据流体力学理论与天然气注采涉及的运动状态，建立各组分传热、流动与扩 散数学模型，计算管壁处压力、温度、速度梯度、各组分浓度分布以及剪切力； 步骤二:根据二氧化碳腐蚀理论，在对气井管柱腐蚀数据的详细分析的基础上，建立二 氧化碳腐蚀模型，确定二氧化碳腐蚀过程中腐蚀产物层对于该环境下电化学腐蚀速率的影 响，引入腐蚀层厚度对腐蚀速率的影响因子； 步骤三:根据冲刷磨损理论，对天然气井注采数据的详细分析基础上，建立注采过程的 冲蚀模型，以流体能量损失数据为基础，采用能量磨损理论建立冲蚀模型，并确定冲刷过程 对腐蚀产物层及管壁基体冲刷磨损速率系数； 步骤四：将以上步骤一至步骤三所述各模型进行耦合计算，计算出管内的温度、压力、 入口流速对腐蚀速率的影响，确定天然气气井管柱的安全服役参数。2. 根据权利要求1所述的一种冲蚀-二氧化碳腐蚀耦合作用下的天然气注采井管柱腐 蚀速率的确定方法，其特征在于，步骤一中所描述流动、传热模型采用流体力学N-S方程计 算;液滴在气流中的运动引入Wu曳力公式及牛顿第二定律计算;物质在气流中的扩散采用 对流扩散模型计算;通过上述模型的共同计算，可以得到管壁处压力、温度、速度梯度、各组 分浓度分布以及剪切力。3. 根据权利要求1所述的一种冲蚀-二氧化碳腐蚀耦合作用下的天然气注采井管柱腐 蚀速率的确定方法，其特征在于，步骤二中所描述的二氧化碳腐蚀模型如下：式中：a，b，d腐蚀速率的修正参数，T是温度，A02是⑶2分压，pH是介质值，Cccif是腐蚀速 率修正系数;V。是二氧化碳腐蚀速率;然后引入腐蚀产物层对腐蚀速率的影响因子来对腐 蚀速率进行修正:￥/。=[妒 <1(1+(：)+〇]￥。，式中：〇，0是产物层对腐蚀速率的修正因子，]1(1是腐 蚀产物层厚度，c是腐蚀速率修正系数，V。是修正后的腐蚀速率。4. 根据权利要求3所述的一种冲蚀-二氧化碳腐蚀耦合作用下的天然气注采井管柱腐 蚀速率的确定方法，其特征在于，所述二氧化碳腐蚀模型其中的a，b，d是根据实际腐蚀速率 代入模型确定的参数;α，β是通过建立实际腐蚀速率以及腐蚀产物层厚度关系得到。5. 根据权利要求1所述的一种冲蚀-二氧化碳腐蚀耦合作用下的天然气注采井管柱腐 蚀速率的确定方法，其特征在于，步骤三中所述的根据冲刷磨损理论所建立的冲蚀模型如 下：其中Hf是磨损速率;α是磨损速率系数;Fw是剪切力，μ是介质流体动力粘度。6. 根据权利要求5所述的一种冲蚀-二氧化碳腐蚀耦合作用下的天然气注采井管柱腐 蚀速率的确定方法，其特征在于，所述冲蚀模型其中的α是根据实际冲蚀速率代入模型确定 的参数。7. 根据权利要求1所述的一种冲蚀-二氧化碳腐蚀耦合作用下的天然气注采井管柱腐 蚀速率的确定方法，其特征在于，步骤四中所述親合计算方法为:通过流动、传热以及传质 模型计算出管壁处压力、温度、速度梯度各组分浓度分布以及剪切力;将上述参数导入二氧 化碳腐蚀模型，计算出腐蚀速率以及腐蚀产物的生成速率;通过冲蚀模型计算出在流体的 作用下腐蚀产物以及管壁的冲蚀速率;将冲蚀模型计算结果再重新导入二氧化碳腐蚀模型 的腐蚀产物层修正项，算得最终的腐蚀速率。8. 根据权利要求1所述的一种冲蚀-二氧化碳腐蚀耦合作用下的天然气注采井管柱腐 蚀速率的确定方法，其特征在于，通过导入不同的环境参数，算得对应的管柱腐蚀速率，选 择安全的腐蚀速率，确定步骤四中所述天然气井管柱的安全服役参数。9. 根据权利要求1所述的一种冲蚀-二氧化碳腐蚀耦合作用下的天然气注采井管柱腐 蚀速率的确定方法，其特征在于，步骤一中所述的传热、流动与扩散数学模型的求解采用 fluent软件计算;而二氧化碳腐蚀、冲刷磨损模型采用fluent软件的UDF用户自定义函数进 行親合计算。10. 根据权利要求9所述的一种冲蚀-二氧化碳腐蚀耦合作用下的天然气注采井管柱腐 蚀速率的确定方法，其特征在于，使用Fluent软件，将天然气井管柱模型导入后将电化学腐 蚀-冲刷腐蚀耦合机理模型写入UDF，将实际环境中的压力、组分浓度、温度以及流速参数输 入软件进行双精度计算即可得到计算结果。
【文档编号】G06F17/50GK106021659SQ201610304136
【公开日】2016年10月12日
【申请日】2016年5月10日
【发明人】刘伟, 付钢旦, 杨斌, 李琼玮, 李治, 董晓焕, 张振云, 朱方辉, 李慧, 孙雨来
【申请人】中国石油天然气股份有限公司