CN115099508B - Srb与co2耦合作用的页岩气集输管道腐蚀速率预测方法 - Google Patents
Srb与co2耦合作用的页岩气集输管道腐蚀速率预测方法 Download PDFInfo
- Publication number
- CN115099508B CN115099508B CN202210773683.9A CN202210773683A CN115099508B CN 115099508 B CN115099508 B CN 115099508B CN 202210773683 A CN202210773683 A CN 202210773683A CN 115099508 B CN115099508 B CN 115099508B
- Authority
- CN
- China
- Prior art keywords
- corrosion
- current density
- formula
- equation
- srb
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000005260 corrosion Methods 0.000 title claims abstract description 139
- 230000007797 corrosion Effects 0.000 title claims abstract description 139
- 238000000034 method Methods 0.000 title claims abstract description 40
- 230000001808 coupling effect Effects 0.000 title claims abstract description 13
- 230000009467 reduction Effects 0.000 claims abstract description 51
- 230000007246 mechanism Effects 0.000 claims abstract description 29
- 102000018832 Cytochromes Human genes 0.000 claims abstract description 19
- 108010052832 Cytochromes Proteins 0.000 claims abstract description 19
- 230000008569 process Effects 0.000 claims abstract description 13
- 238000006722 reduction reaction Methods 0.000 claims description 61
- XEEYBQQBJWHFJM-UHFFFAOYSA-N iron Substances [Fe] XEEYBQQBJWHFJM-UHFFFAOYSA-N 0.000 claims description 50
- 238000012546 transfer Methods 0.000 claims description 43
- 238000006243 chemical reaction Methods 0.000 claims description 37
- 238000004364 calculation method Methods 0.000 claims description 29
- 239000007789 gas Substances 0.000 claims description 25
- 238000009792 diffusion process Methods 0.000 claims description 22
- 230000004913 activation Effects 0.000 claims description 19
- 230000027756 respiratory electron transport chain Effects 0.000 claims description 15
- 229910052742 iron Inorganic materials 0.000 claims description 12
- 239000002184 metal Substances 0.000 claims description 12
- 229910052751 metal Inorganic materials 0.000 claims description 12
- 238000003487 electrochemical reaction Methods 0.000 claims description 11
- 239000012528 membrane Substances 0.000 claims description 9
- QAOWNCQODCNURD-UHFFFAOYSA-L Sulfate Chemical compound [O-]S([O-])(=O)=O QAOWNCQODCNURD-UHFFFAOYSA-L 0.000 claims description 8
- 238000007254 oxidation reaction Methods 0.000 claims description 7
- 210000000805 cytoplasm Anatomy 0.000 claims description 6
- 238000006703 hydration reaction Methods 0.000 claims description 6
- 239000001257 hydrogen Substances 0.000 claims description 5
- 229910052739 hydrogen Inorganic materials 0.000 claims description 5
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 5
- UFHFLCQGNIYNRP-UHFFFAOYSA-N Hydrogen Chemical compound [H][H] UFHFLCQGNIYNRP-UHFFFAOYSA-N 0.000 claims description 4
- 230000032258 transport Effects 0.000 claims description 4
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 claims description 3
- 230000009056 active transport Effects 0.000 claims description 3
- 150000001875 compounds Chemical class 0.000 claims description 3
- XLYOFNOQVPJJNP-ZSJDYOACSA-N heavy water Substances [2H]O[2H] XLYOFNOQVPJJNP-ZSJDYOACSA-N 0.000 claims description 3
- 229910000015 iron(II) carbonate Inorganic materials 0.000 claims description 3
- 230000005540 biological transmission Effects 0.000 abstract description 6
- 238000011156 evaluation Methods 0.000 abstract description 2
- 239000000243 solution Substances 0.000 description 6
- 230000001580 bacterial effect Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 150000002500 ions Chemical class 0.000 description 2
- 230000003647 oxidation Effects 0.000 description 2
- 238000011002 quantification Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 108091006149 Electron carriers Proteins 0.000 description 1
- 241000364021 Tulsa Species 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000003181 biological factor Substances 0.000 description 1
- 239000008364 bulk solution Substances 0.000 description 1
- 238000010349 cathodic reaction Methods 0.000 description 1
- 210000004027 cell Anatomy 0.000 description 1
- 210000002421 cell wall Anatomy 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000010351 charge transfer process Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000005553 drilling Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000006056 electrooxidation reaction Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 150000002431 hydrogen Chemical class 0.000 description 1
- 230000000977 initiatory effect Effects 0.000 description 1
- 230000009545 invasion Effects 0.000 description 1
- 230000002503 metabolic effect Effects 0.000 description 1
- 244000005700 microbiome Species 0.000 description 1
- 238000005065 mining Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000002070 nanowire Substances 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 239000000758 substrate Substances 0.000 description 1
- 230000004083 survival effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N17/00—Investigating resistance of materials to the weather, to corrosion, or to light
- G01N17/02—Electrochemical measuring systems for weathering, corrosion or corrosion-protection measurement
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
- G16C20/00—Chemoinformatics, i.e. ICT specially adapted for the handling of physicochemical or structural data of chemical particles, elements, compounds or mixtures
- G16C20/10—Analysis or design of chemical reactions, syntheses or processes
Landscapes
- Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- Economics (AREA)
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Human Resources & Organizations (AREA)
- Strategic Management (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Marketing (AREA)
- General Health & Medical Sciences (AREA)
- Analytical Chemistry (AREA)
- Tourism & Hospitality (AREA)
- General Business, Economics & Management (AREA)
- Chemical Kinetics & Catalysis (AREA)
- Environmental Sciences (AREA)
- Operations Research (AREA)
- Crystallography & Structural Chemistry (AREA)
- Entrepreneurship & Innovation (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Computing Systems (AREA)
- Biodiversity & Conservation Biology (AREA)
- Ecology (AREA)
- Environmental & Geological Engineering (AREA)
- Quality & Reliability (AREA)
- Game Theory and Decision Science (AREA)
- Biochemistry (AREA)
- Development Economics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Public Health (AREA)
- Water Supply & Treatment (AREA)
- Primary Health Care (AREA)
- Testing Resistance To Weather, Investigating Materials By Mechanical Methods (AREA)
Abstract
本发明公开了一种SRB与CO2耦合作用的页岩气集输管道腐蚀速率预测方法,在综合SRB和CO2腐蚀因素的基础上,考虑以细胞色素传递机制量化SRB腐蚀电化学过程中的还原引起的腐蚀电流密度,建立一个SRB/CO2共存环境的页岩气集输管道内腐蚀机理模型,以期充分预测SRB/CO2共存环境下的管道内腐蚀速率,由于在获得页岩气集输管道腐蚀速率的过程中考虑了SRB和CO2耦合作用,以及考虑细胞色素传递机制的影响,可以使评价结果更加科学、准确、合理。
Description
技术领域
本发明涉及管道腐蚀速率预测技术领域,特别是一种SRB与CO2耦合作用的页岩气集输管道腐蚀速率预测方法。
背景技术
页岩气集输管道内腐蚀是影响在役管道剩余寿命以及页岩气集输流程的重要原因之一。页岩气因为其特殊的开采方式,压裂液,钻井泥浆和蓄水都可能将有害的微生物如SRB引入页岩气藏。同时,页岩气集输环境中又含有CO2、Cl-等腐蚀性介质,加剧页岩气集输管道的电化学腐蚀,造成管道内腐蚀甚至穿孔,引发安全事故和财产损失。因此,在内检测技术受限的情况下,通过形成一套内腐蚀速率预测技术体系,也对保障页岩气集输安全运行具有重要的研究与工程意义。
引起页岩气集输管道内腐蚀的原因比较复杂,最主要的两个因素为SRB和CO2,国内外学者对于CO2腐蚀速率预测模型的研究颇多,多用常规的数学建模方法,如常见的经验模型有NorsokM-506模型、Ohio模型以及Corpos模型,常见半经验模型的有De Waard-Milliams模型、Cassandra模型、ECE模型、Predict模型等,常见的机理模型有Nesic模型、Tulsa模型等。相较之下SRB腐蚀速率预测模型的研究很少,基本为机理模型,且SRB机理模型大多是基于模拟腐蚀金属基底上的生物膜模型,并包括表征系统的各种参数,例如生物因素、输运过程、金属表面的化学反应等建立的,而关于同时考虑SRB和CO2的内腐蚀速率预测的相关内容较少见到。胞外电子传递机制有3种,分别为可溶性氧化还原介体的电子穿梭、细胞色素传递和纳米导线传递,而目前几乎未见将某种EET机制引入SRB机理模型的研究。
发明内容
为解决现有技术中存在的问题,本发明提供了一种SRB与CO2耦合作用的页岩气集输管道腐蚀速率预测方法,在综合SRB和CO2腐蚀因素的基础上,考虑以细胞色素传递机制量化SRB腐蚀电化学过程中的还原引起的腐蚀电流密度,建立一个SRB/CO2共存环境的页岩气集输管道内腐蚀机理模型,以期充分预测SRB/CO2共存环境下的管道内腐蚀速率。
本发明采用的技术方案是:
一种SRB与CO2耦合作用的页岩气集输管道腐蚀速率预测方法,包括如下步骤:
确定SRB/CO2溶液体系腐蚀过程中的电化学反应,包括CO2腐蚀的反应方程、SRB腐蚀的氧化反应方程、SRB腐蚀的还原反应方程,得出SRB与CO2腐蚀的总反应方程;
结合胞外电子传递机制中的细胞色素传递方法,建立还原的电荷转移电流密度方程及/>的传质模型,计算/>还原的总电流密度;
考虑H2CO3还原反应同时受活化反应和扩散传质控制,计算H2CO3还原的总电流密度;
考虑H2O还原受活化反应控制,计算H2O还原的总电流密度;
计算得到各电化学反应总的腐蚀电流密度,根据阳极和阴极反应的电荷平衡计算得到阳极腐蚀电位,并利用阳极腐蚀电位计算出不同时间的阳极腐蚀电流密度,将阳极腐蚀电流密度转换为腐蚀速率与腐蚀深度,最终计算出腐蚀速率。
优选的,CO2腐蚀的反应方程分别见式(1)和式(2):
2H2O+2e-→H2+2OH- (2)
SRB腐蚀的氧化反应方程分别见式(3)和(4):
CH3CHOHCOO-+H2O→CH3COO-+CO2+4H++4e- (3)
4Fe→4Fe2++8e- (4)
SRB腐蚀的还原反应方程见式(5):
SO4 2-+9H++8e-→HS-+4H2O (5)
SRB与CO2腐蚀的总反应方程分别见式(6)、(7)和(8):
Fe2++HS-→FeS+H+ (6)
Fe+CO2+H2O→FeCO3+H2 (8)
优选的,结合胞外电子传递机制中的细胞色素传递方法,建立还原的电荷转移电流密度方程及/>的传质模型,计算/>还原的极限扩散电流密度,包括如下步骤:
根据式(14)的总电阻方程计算硫酸盐还原产生的总电流密度式(14)如下所示:
其中,为/>还原的总电流密度;/>为/>还原的极限扩散电流密度;/>为/>还原的电荷转移电流密度;
电子从金属到生物膜的转移速率是跨生物膜的活性电荷转移层的电流密度和金属表面电势的函数,则通过Tafel方程建立的表达式即为式(15):
其中,rc为电子从金属到生物膜的跨电位降的传输速率;F为Faraday常数;为/>还原的交换电流密度;α为电子转移系数;E为阳极腐蚀电位;/>为/>还原平衡电位;R为通用气体常数;T为实际工况下的绝对温度;
使用Nerst方程计算见式(16):
其中,pH为SBR细胞质中的pH值;
最终,如下式(17):
其中,为细胞色素与生物膜之间电子转移的速率;LB为主动传输层长度;CZ为细胞色素浓度;
由Nernst-Planck方程计算见式(18):
其中,ne为参与反应的电子数;为生物膜中/>的扩散率;/>为生物膜内的/>浓度;
通过修正后的Fick第二定律可模拟在生物膜内部的分布,见式(19):
其中,为生物膜中SRB对/>消耗的速率,通过Monod方程计算,见式(20):
其中,Km为生物膜中SRB的Monod半速度系数;qmax为最大硫酸盐利用率;Xs为点蚀坑内的SRB浓度;Lf为生物膜厚度。
优选的,考虑H2CO3还原反应同时受活化反应和扩散传质控制,对H2CO3的极限扩散电流密度和电荷转移电流密度进行计算,包括如下步骤:
H2CO3还原的总阴极电流密度见式(21):
其中,为H2CO3还原的总电流密度;/>为H2CO3还原的电荷转移电流密度;/>为H2CO3还原的极限扩散电流密度;
通过Tafel方程建立计算方程见式(22):
其中,为H2CO3还原的交换电流密度;bc为Tafel斜率;/>为H2CO3还原的平衡电位;
通过式(10)计算bc:
其中,αc为电子转移系数;
通过式(23)计算
其中,为氢气的分压;
通过式(24)计算
其中,为H2CO3还原的参考交换电流密度;/>为H2CO3还原的活化焓;Tref为参考温度;
通过Vetter方程建立的计算方式见式(25):
其中,[CO2]b为溶液中溶解的CO2体积浓度;为CO2在水中的扩散率;Khyd为CO2水合反应的平衡常数;/>为CO2水合反应的正向反应速率;
通过式(26)计算[CO2]b:
其中,为CO2的亨利常数;/>为CO2分压。
优选的,考虑H2O还原受活化反应控制,对H2O的电荷转移电流密度进行计算,包括如下步骤:
通过Tafel方程对H2O还原的总电流密度进行计算,计算式见式(27):
其中,为H2O还原的总电流密度;/>为H2O还原的交换电流密度;/>为H2O还原的平衡电位;H2O还原的Tafel斜率与H2CO3还原的Tafel斜率相同;
通过式(28)计算
通过Arrhenius方程计算计算式见式(29):
其中,为H2O还原的参考交换电流密度;/>为H2O还原的活化焓。
优选的,计算得到各电化学反应总的腐蚀电流密度,根据阳极和阴极反应的电荷平衡计算得到阳极腐蚀电位,并利用阳极腐蚀电位计算出不同时间的阳极腐蚀电流密度,将阳极腐蚀电流密度转换为腐蚀速率与腐蚀深度,最终计算出腐蚀速率,包括如下步骤:
通过Tafel方程计算阳极腐蚀电流密度iFe,见式(9):
其中,iFe为阳极腐蚀电流密度;io,Fe为阳极交换电流密度;为阳极平衡电位;
通过式(12)计算io,Fe:
其中,为Fe的参考交换电流密度;ΔH(Fe)为Fe的活化焓;
根据Faraday第二定律将阳极腐蚀电流密度转换为腐蚀速率Vcorr,见式(13)所示,此处的腐蚀速率Vcorr相当于一维点蚀模型中的点蚀速率;
其中,MFe为铁的原子量;ρFe为铁的密度;
通过电荷平衡方程计算E,见式(30):
其中,na和nc分别为阳极和阴极反应的总数;ia和ic分别为阳极电荷和阴极电荷;展开后即为式(31):
通过有限差分法对式(19)求解,得到不同时间下生物膜上浓度,再将T、/>pH、SRB数量分别代入包含该参数的各公式中,将式(15)和(18)代入式(14)得到式(32),将式(22)和(25)代入式(21)得到式(33),式(32)和式(33)如下所示:
此时式(9)、式(27)、式(32)和式(33)中的未知数仅为阳极腐蚀电位E;将式(9)、式(27)、式(32)和式(33)通过电荷平衡方程式(31)联立,求解阳极出腐蚀电位E,再以阳极腐蚀电位E作为式(9)的输入,即可求出不同时间的阳极腐蚀电流密度iFe,最后将阳极腐蚀电流密度iFe代入式(13)即可求解腐蚀速率Vcorr。
本发明的有益效果是:
考虑了以细胞色素传递机制量化SRB腐蚀电化学过程中的还原引起的腐蚀电流密度;建立了SRB/CO2共存环境下页岩气集输管道内腐蚀机理模型;由于在获得页岩气集输管道腐蚀速率的过程中考虑了SRB和CO2耦合作用,以及考虑细胞色素传递机制的影响,可以使评价结果更加科学、准确、合理。
具体实施方式
下面结合实施例对本发明的实施例进行详细说明。
实施例
本实施例提供一种SRB与CO2耦合作用的页岩气集输管道腐蚀速率预测方法,该方法考虑细胞色素传递机制的SRB与CO2耦合作用,以预测页岩气集输管道的内腐蚀速率趋势,保障页岩气集输管道安全运行
1、首先确定SRB/CO2溶液体系腐蚀过程中主要的电化学反应。
SRB和CO2共存时诱发的管道腐蚀机理比较复杂。首先分析CO2的腐蚀机理,传统上认为CO2腐蚀的阴极反应涉及析氢反应、H2CO3的直接还原反应以及H2O还原反应,根据现场水质分析发现不涉及析氢反应。
CO2腐蚀的反应方程分别见式(1)和式(2):
2H2O+2e-→H2+2OH- (2)
SRB腐蚀的阴极反应,目前国内外最先进的SRB腐蚀机理是顾停月等人发现的BCSR机理,由于铁氧化发生在SRB外部,而硫酸盐还原发生在SRB内部,因此需要EET的方法使电子穿过SRB的细胞壁,考虑以EET机制中的细胞色素传递作为电子载体的方式,对SO4 2-还原的阴极电流密度进行建模,由此量化SO4 2-还原中因EET产生的腐蚀电流密度。
SRB腐蚀的反应方程如下:
氧化反应方程分别见式(3)和(4):
CH3CHOHCOO-+H2O→CH3COO-+CO2+4H++4e- (3)
4Fe→4Fe2++8e- (4)
还原反应方程见式(5):
SO4 2-+9H++8e-→HS-+4H2O (5)
因此,SRB与CO2腐蚀的总反应方程分别见式(6)、(7)和(8):
Fe2++HS-→FeS+H+ (6)
Fe+CO2+H2O→FeCO3+H2(8)
2、通过Tafel方程计算阳极铁氧化的电化学反应的腐蚀电流密度;阳极表面的电流密度可由Tafel方程计算:见式(9):
其中,iFe为阳极腐蚀电流密度,单位为A/m2;io,Fe为阳极交换电流密度,单位为A/m2;bc为Tafel斜率,无量纲;E为阳极腐蚀电位,单位为V;为阳极平衡电位,单位为V。
通过式(10)计算bc:
其中,αc为电子转移系数,取0.5,无量纲;F为Faraday常数,取96485.34C/mol;R为通用气体常数,取8.31J/(mol·K);T为实际工况下的绝对温度,单位为K。
如果实验条件偏离标准电极电位定义的条件,可使用Nerst方程计算平衡电位,见式(11):
对于交换电流密度,Pour-Ghaz等提出了一种简单的方法揭示温度对交换电流密度的影响,见式(12):
其中,为25℃下Fe的参考交换电流密度,取0.33A/m2;ΔH(Fe)为Fe的活化焓,取37.5KJ/mol;Tref为参考温度,代表的是不同实验研究下所取得不同参考温度值,本实施例中Tref统一取同一值为298K。
根据Faraday第二定律将阳极腐蚀电流密度转换为腐蚀速率Vcorr,见式(13)所示,此处的腐蚀速率Vcorr相当于一维点蚀模型中的点蚀速率。
其中,腐蚀速率Vcorr的单位为mm/a;ne为参与反应的电子数,取2;ρFe为铁的密度,取7870kg/cm3;MFe为铁的原子量,取0.056kg/mol。
3、结合胞外电子传递机制中的细胞色素传递方法,通过Tafel方程建立还原的电荷转移电流密度,利用Monod方程修正Fick第二定律建立/>的传质模型,再根据Nernst-Planck方程计算/>还原的极限扩散电流密度。
硫酸盐还原产生的总电流密度根据式(14)总电阻方程计算:
其中,为/>还原的总电流密度,单位为A/m2;/>为/>还原的极限扩散电流密度,单位为A/m2;/>为/>还原的电荷转移电流密度,单位为A/m2。
电子从金属到生物膜的转移速率是跨生物膜的活性电荷转移层的电流密度和金属表面电势的函数。
则通过Tafel方程建立的还原的电荷转移电流密度的表达式即为式(15):
其中,rc为电子从金属到生物膜的跨电位降的传输速率,单位为s-1;为还原的交换电流密度;α为电子转移系数,取0.5;/>为/>还原平衡电位,单位为V。
使用Nerst方程计算见式(16):
其中,本实施例中pH代表SRB细胞质中的酸碱度,的还原发生在细胞质中,SRB细胞质中的酸碱度不同于溶液环境的酸碱度,因为受细菌代谢活动的影响,细胞质内部会保持最适宜细菌生存的酸碱度。本实施例中假设细胞质中的pH=7。
最终,如下式(17):
其中,为细胞色素与生物膜之间电子转移的速率,取0.08s-1;LB为主动传输层长度,取10-3m;CZ为细胞色素浓度,取3.6×10-3mol/L。
由Nernst-Planck方程计算见式(18):
其中,为生物膜中/>的扩散率,单位为m2/s;/>为生物膜内的/>浓度,单位为mol/m3。
Fick第二定律是描述离子在扩散机理作用下离子浓度随时间的变化,能很好地反映生物膜中离子侵入的过程。通过修正后的Fick第二定律可模拟在生物膜内部的分布,见式(19):
其中,为生物膜中SRB对/>消耗的速率,单位为mol/(m3·s)。
Fick第二定律初始与边界条件见表所示。
表1 Fick第二定律初始与边界条件
对于腐蚀深坑,坑底和坑壁可能存在SRB生物膜,而坑的大部分空间填充有主体溶液。从溶液进入生物膜,到达坑底接受电子,导致MIC发生。在扩散过程中,生物膜中大量固着细胞消耗了硫酸盐,/>应表现为负值,该项通过Monod方程计算,见式(20):
其中,Km为生物膜中SRB的Monod半速度系数,取10.69×10-5mg/m3;qmax为最大硫酸盐利用率,取1.08mmol/(mg·s-1);Xs为点蚀坑内的SRB浓度,假设点蚀坑中的SRB浓度为溶液中SRB浓度的1%,单位为个/mL;Lf为生物膜厚度,取239μm。
4、考虑H2CO3还原反应同时受活化反应和扩散传质控制,通过Tafel方程和Vetter方程对H2CO3的极限扩散电流密度和电荷转移电流密度进行计算。
H2CO3还原的总阴极电流密度见式(21):
其中,为H2CO3还原的总电流密度,单位为A/m2;/>为H2CO3还原的电荷转移电流密度,单位为A/m2;/>为H2CO3还原的极限扩散电流密度,单位为A/m2。
通过Tafel方程建立计算方程见式(22):
其中,为H2CO3还原的交换电流密度,单位为A/m2;/>为H2CO3还原的平衡电位,单位为V。
根据研究,H2CO3还原与H+还原具有相同的热力学条件。通过式(23)计算
其中,为氢气的分压,在本实施例中为简化研究过程,/>取值为0。
通过式(24)计算
其中,为H2CO3还原的参考交换电流密度,根据/>的研究,在本实施例中,H2CO3浓度为1×10-4mol/L的参考条件下,/>为0.018A/m2;/>为H2CO3还原的活化焓,取50KJ/mol。
通过Vetter方程建立的计算方式见式(25):
其中,[CO2]b为溶解的CO2体积浓度,单位为mol/m3;为CO2在水中的扩散率,单位为m2/s;Khyd为CO2水合反应的平衡常数,取2.58×10-3,无量纲;/>为CO2水合反应的正向反应速率,无量纲。
通过式(26)计算[CO2]b:
其中,为CO2的亨利常数,取0.1MPa·m3/mol;/>为CO2分压,单位为MPa。
5、考虑H2O还原受活化反应控制,通过Tafel方程对H2O的电荷转移电流密度进行计算;由于水分子向金属表面传输不存在扩散限制,因此可以认为H2O的还原速率仅由Tafel方程建立的电荷转移过程控制:
其中,为H2O还原的总电流密度,单位为A/m2;/>为H2O还原的交换电流密度,单位为A/m2;/>为H2O还原的平衡电位,单位为V。
H2O还原的Tafel斜率与H2CO3还原的Tafel斜率相同,bc=120mV。H2O还原与H+还原具有相同的热力学条件。通过Nernst方程计算:
通过Arrhenius方程计算计算式见式(29):/>
其中,在为H2O还原的参考交换电流密度,在Tref环境下的实验中确定 为H2O还原的活化焓,/>
6、根据步骤2~步骤5计算可得到各电化学反应总的腐蚀电流密度,利用阳极和阴极反应的电荷平衡计算得到腐蚀电位,进一步通过利用Faraday第二定律求得腐蚀速率。
在阴极和阳极反应建模完成之后,通过电荷平衡方程计算腐蚀电位,见式(30):
其中,na和nc分别为阳极和阴极反应的总数;ia和ic分别为阳极电荷和阴极电荷;展开后即为式(31):
通过有限差分法对式(19)求解,得到不同时间下生物膜上浓度,再将T、pH、SRB数量分别代入包含该参数的各公式中,将式(15)和(18)代入式(14)得到式(32),将式(22)和(25)代入式(21)得到式(33),式(32)和式(33)如下所示:
此时式(9)、式(27)、式(32)和式(33)中的未知数仅为阳极腐蚀电位E;将式(9)、式(27)、式(32)和式(33)通过电荷平衡方程式(31)联立,求解出阳极腐蚀电位E,再以阳极腐蚀电位E作为式(9)的输入,即可求出不同时间的阳极腐蚀电流密度iFe,最后将阳极腐蚀电流密度iFe代入式(13)即可求解腐蚀速率Vcorr。
通过本实施例提供的SRB与CO2耦合作用的页岩气集输管道腐蚀速率预测方法,对多组样本进行腐蚀速率预测,其预测结果如表2所示。
表2管道腐蚀数据及机理模型预测结果
/>
由表2可知,本实施例的机理模型预测的点蚀速率与现场实测点蚀速率的大部分相对误差在10%以内,只有个别相对误差在10%~15%,总体相对误差较小,验证了本实施例的SRB/CO2共存环境的页岩气集输管道内腐蚀机理模型的有效性。
以上所述实施例仅表达了本发明的具体实施方式,其描述较为具体和详细,但并不能因此而理解为对本发明专利范围的限制。应当指出的是,对于本领域的普通技术人员来说,在不脱离本发明构思的前提下,还可以做出若干变形和改进,这些都属于本发明的保护范围。
Claims (2)
1.一种SRB与CO2耦合作用的页岩气集输管道腐蚀速率预测方法,其特征在于,包括如下步骤:
确定SRB/CO2溶液体系腐蚀过程中的电化学反应,包括CO2腐蚀的反应方程、SRB腐蚀的氧化反应方程、SRB腐蚀的还原反应方程,得出SRB与CO2腐蚀的总反应方程;
结合胞外电子传递机制中的细胞色素传递方法,建立还原的电荷转移电流密度方程及/>的传质模型,计算/>还原的总电流密度;
考虑H2CO3还原反应同时受活化反应和扩散传质控制,计算H2CO3还原的总电流密度;
考虑H2O还原受活化反应控制,计算H2O还原的总电流密度;
计算得到各电化学反应总的腐蚀电流密度,根据阳极和阴极反应的电荷平衡计算得到阳极腐蚀电位,并利用阳极腐蚀电位计算出不同时间的阳极腐蚀电流密度,将阳极腐蚀电流密度转换为腐蚀速率与腐蚀深度,最终计算出腐蚀速率;
结合胞外电子传递机制中的细胞色素传递方法,建立还原的电荷转移电流密度方程及/>的传质模型,计算/>还原的极限扩散电流密度,包括如下步骤:
根据式(14)的总电阻方程计算硫酸盐还原产生的总电流密度式(14)如下所示:
其中,为/>还原的总电流密度;/>为/>还原的极限扩散电流密度;为/>还原的电荷转移电流密度;
电子从金属到生物膜的转移速率是跨生物膜的活性电荷转移层的电流密度和金属表面电势的函数,则通过Tafel方程建立的表达式即为式(15):
其中,rc为电子从金属到生物膜的跨电位降的传输速率;F为Faraday常数;为还原的交换电流密度;α为电子转移系数;E为阳极腐蚀电位;/>为/>还原平衡电位;R为通用气体常数;T为实际工况下的绝对温度;
使用Nerst方程计算见式(16):
其中,pH为SBR细胞质中的pH值;
最终,如下式(17):
其中,为细胞色素与生物膜之间电子转移的速率;LB为主动传输层长度;CZ为细胞色素浓度;
由Nernst-Planck方程计算见式(18):
其中,ne为参与反应的电子数;为生物膜中/>的扩散率;/>为生物膜内的浓度;
通过修正后的Fick第二定律可模拟在生物膜内部的分布,见式(19):
其中,为生物膜中SRB对/>消耗的速率,通过Monod方程计算,见式(20):
其中,Km为生物膜中SRB的Monod半速度系数;qmax为最大硫酸盐利用率;Xs为点蚀坑内的SRB浓度;Lf为生物膜厚度;
考虑H2CO3还原反应同时受活化反应和扩散传质控制,对H2CO3的极限扩散电流密度和电荷转移电流密度进行计算,包括如下步骤:
H2CO3还原的总阴极电流密度见式(21):
其中,为H2CO3还原的总电流密度;/>为H2CO3还原的电荷转移电流密度;为H2CO3还原的极限扩散电流密度;
通过Tafel方程建立计算方程见式(22):
其中,为H2CO3还原的交换电流密度;bc为Tafel斜率;/>为H2CO3还原的平衡电位;
通过式(10)计算bc:
其中,αc为电子转移系数;
通过式(23)计算
其中,为氢气的分压;
通过式(24)计算
其中,为H2CO3还原的参考交换电流密度;/>为H2CO3还原的活化焓;Tref为参考温度;
通过Vetter方程建立的计算方式见式(25):
其中,[CO2]b为溶液中溶解的CO2体积浓度;为CO2在水中的扩散率;Khyd为CO2水合反应的平衡常数;/>为CO2水合反应的正向反应速率;
通过式(26)计算[CO2]b:
其中,为CO2的亨利常数;/>为CO2分压;
考虑H2O还原受活化反应控制,对H2O的电荷转移电流密度进行计算,包括如下步骤:
通过Tafel方程对H2O还原的总电流密度进行计算,计算式见式(27):
其中,为H2O还原的总电流密度;/>为H2O还原的交换电流密度;/>为H2O还原的平衡电位;H2O还原的Tafel斜率与H2CO3还原的Tafel斜率相同;
通过式(28)计算
通过Arrhenius方程计算计算式见式(29):
其中,为H2O还原的参考交换电流密度;/>为H2O还原的活化焓;
计算得到各电化学反应总的腐蚀电流密度,根据阳极和阴极反应的电荷平衡计算得到阳极腐蚀电位,并利用阳极腐蚀电位计算出不同时间的阳极腐蚀电流密度,将阳极腐蚀电流密度转换为腐蚀速率与腐蚀深度,最终计算出腐蚀速率,包括如下步骤:
通过Tafel方程计算阳极腐蚀电流密度iFe,见式(9):
其中,iFe为阳极腐蚀电流密度;io,Fe为阳极交换电流密度;为阳极平衡电位;
通过式(12)计算io,Fe:
其中,为Fe的参考交换电流密度;ΔH(Fe)为Fe的活化焓;
根据Faraday第二定律将阳极腐蚀电流密度转换为腐蚀速率Vcorr,见式(13)所示,此处的腐蚀速率Vcorr相当于一维点蚀模型中的点蚀速率;
其中,MFe为铁的原子量;ρFe为铁的密度;
通过电荷平衡方程计算E,见式(30):
其中,na和nc分别为阳极和阴极反应的总数;ia和ic分别为阳极电荷和阴极电荷;展开后即为式(31):
通过有限差分法对式(19)求解,得到不同时间下生物膜上浓度,再将T、/>pH、SRB数量分别代入包含该参数的各公式中,将式(15)和(18)代入式(14)得到式(32),将式(22)和(25)代入式(21)得到式(33),式(32)和式(33)如下所示:
此时式(9)、式(27)、式(32)和式(33)中的未知数仅为阳极腐蚀电位E;将式(9)、式(27)、式(32)和式(33)通过电荷平衡方程式(31)联立,求解出阳极腐蚀电位E,再以阳极腐蚀电位E作为式(9)的输入,即可求出不同时间的阳极腐蚀电流密度iFe,最后将阳极腐蚀电流密度iFe代入式(13)即可求解腐蚀速率Vcorr。
2.根据权利要求1所述SRB与CO2耦合作用的页岩气集输管道腐蚀速率预测方法,其特征在于,CO2腐蚀的反应方程分别见式(1)和式(2):
2H2O+2e-→H2+2OH- (2)
SRB腐蚀的氧化反应方程分别见式(3)和(4):
CH3CHOHCOO-+H2O→CH3COO-+CO2+4H++4e- (3)
4Fe→4Fe2++8e- (4)
SRB腐蚀的还原反应方程见式(5):
SO4 2-+9H++8e-→HS-+4H2O (5)
SRB与CO2腐蚀的总反应方程分别见式(6)、(7)和(8):
Fe2++HS-→FeS+H+ (6)
Fe+CO2+H2O→FeCO3+H2 (8)。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210773683.9A CN115099508B (zh) | 2022-07-01 | 2022-07-01 | Srb与co2耦合作用的页岩气集输管道腐蚀速率预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210773683.9A CN115099508B (zh) | 2022-07-01 | 2022-07-01 | Srb与co2耦合作用的页岩气集输管道腐蚀速率预测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115099508A CN115099508A (zh) | 2022-09-23 |
CN115099508B true CN115099508B (zh) | 2024-06-07 |
Family
ID=83294101
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210773683.9A Active CN115099508B (zh) | 2022-07-01 | 2022-07-01 | Srb与co2耦合作用的页岩气集输管道腐蚀速率预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115099508B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115099151B (zh) * | 2022-07-01 | 2024-06-11 | 西南石油大学 | 基于gwo-bp修正机理模型的管道腐蚀速率预测方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102713280A (zh) * | 2009-08-27 | 2012-10-03 | 麦卡利斯特技术有限责任公司 | 通过可再生能的一体化全面生产实现可持续经济发展的系统和方法 |
WO2012149487A1 (en) * | 2011-04-27 | 2012-11-01 | Ohio University | Methods and devices for the detection of biofilms |
CN109668820A (zh) * | 2019-01-22 | 2019-04-23 | 北京市燃气集团有限责任公司 | 埋地管道外腐蚀风险综合评判方法 |
CN111487184A (zh) * | 2020-04-11 | 2020-08-04 | 华东交通大学 | 一种模拟地铁列车杂散电流泄漏及腐蚀的实验平台 |
CN111798930A (zh) * | 2020-07-17 | 2020-10-20 | 西南石油大学 | 一种考虑腐蚀产物膜影响的co2腐蚀速率预测方法 |
CN216160421U (zh) * | 2021-08-18 | 2022-04-01 | 中国科学院金属研究所 | 一种输水管道内微生物膜监测探头 |
-
2022
- 2022-07-01 CN CN202210773683.9A patent/CN115099508B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102713280A (zh) * | 2009-08-27 | 2012-10-03 | 麦卡利斯特技术有限责任公司 | 通过可再生能的一体化全面生产实现可持续经济发展的系统和方法 |
WO2012149487A1 (en) * | 2011-04-27 | 2012-11-01 | Ohio University | Methods and devices for the detection of biofilms |
CN109668820A (zh) * | 2019-01-22 | 2019-04-23 | 北京市燃气集团有限责任公司 | 埋地管道外腐蚀风险综合评判方法 |
CN111487184A (zh) * | 2020-04-11 | 2020-08-04 | 华东交通大学 | 一种模拟地铁列车杂散电流泄漏及腐蚀的实验平台 |
CN111798930A (zh) * | 2020-07-17 | 2020-10-20 | 西南石油大学 | 一种考虑腐蚀产物膜影响的co2腐蚀速率预测方法 |
CN216160421U (zh) * | 2021-08-18 | 2022-04-01 | 中国科学院金属研究所 | 一种输水管道内微生物膜监测探头 |
Non-Patent Citations (3)
Title |
---|
Accelerated corrosion of pipeline steel in the presence of Desulfovibrio desulfuricans biofilm due to carbon source deprivation in CO2 saturated medium;Ubong Eduok 等;《Materials Science and Engineering: C》;20190815;第1-8页 * |
Modification of D-W model for corrosion rate of shale gas pipeline under the synergistic corrosion of SRB and CO2;Xingyu Peng 等;《Anti-Corrosion Methods and Materials》;20210702;第150-159页 * |
基于LHS-Kriging-DW的页岩气集输管道内腐蚀预测模型研究;周逸轩 等;《 热加工工艺》;20230919;第53卷(第16期);第113-117页 * |
Also Published As
Publication number | Publication date |
---|---|
CN115099508A (zh) | 2022-09-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Marciales et al. | Mechanistic microbiologically influenced corrosion modeling—A review | |
Lekbach et al. | Microbial corrosion of metals: The corrosion microbiome | |
Hall et al. | An evaluation of corrosion processes affecting copper-coated nuclear waste containers in a deep geological repository | |
Jeremiasse et al. | Acetate enhances startup of a H2‐producing microbial biocathode | |
Breck | Redox levels in the sea | |
CN115099508B (zh) | Srb与co2耦合作用的页岩气集输管道腐蚀速率预测方法 | |
Avci et al. | A new analysis of two phase flow on hydrogen production from water electrolysis | |
Orcutt et al. | Constraints on mechanisms and rates of anaerobic oxidation of methane by microbial consortia: process-based modeling of ANME-2 archaea and sulfate reducing bacteria interactions | |
Borole et al. | Proton transfer in microbial electrolysis cells | |
Downing et al. | Effect of bulk liquid BOD concentration on activity and microbial community structure of a nitrifying, membrane-aerated biofilm | |
Deng et al. | Vertical galvanic corrosion of pipeline steel in simulated marine thermocline | |
Ma et al. | Biomachining properties of various metals by microorganisms | |
Tognia et al. | Multi-physics modelling and simulation approach with experimental validation for electrocatalytic filtration process: Part A. Optimization of electrodeposition of metal oxides on carbon membranes | |
Lan et al. | High current density with spatial distribution of Geobacter in anodic biofilm of the microbial electrolysis desalination and chemical-production cell with enlarged volumetric anode | |
Peng et al. | Nitrous oxide production in a granule-based partial nitritation reactor: a model-based evaluation | |
Landolt et al. | A review of materials and corrosion issues regarding canisters for disposal of spent fuel and high-level waste in Opalinus Clay | |
CN115099151B (zh) | 基于gwo-bp修正机理模型的管道腐蚀速率预测方法 | |
Żurek et al. | Morphology Analysis of the Scale Formed on Crofer 22APU Steel in Atmospheres Containing SO 2 | |
Liu et al. | Experimental study on corrosion resistance of coiled tubing welds in high temperature and pressure environment | |
Prakash et al. | Sediment diagenesis module for CE-QUAL-W2 Part 2: numerical formulation | |
YG et al. | Biofilms and beyond: a comprehensive review of the impact of Sulphate Reducing Bacteria on steel corrosion | |
Moß et al. | Optimal Geometric Parameters for 3D Electrodes in Bioelectrochemical Systems: A Systematic Approach | |
Zhongying et al. | Investigations on corrosion behaviors of super-strength sucker rod FG20 steel in high SO42− environment | |
León‐Fernandez et al. | Modelling of a bioelectrochemical system for metal‐polluted wastewater treatment and sequential metal recovery | |
Prakash et al. | The oil sands pit lake model–sediment diagenesis module |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant |