CN115099508A - Srb与co2耦合作用的页岩气集输管道腐蚀速率预测方法 - Google Patents
Srb与co2耦合作用的页岩气集输管道腐蚀速率预测方法 Download PDFInfo
- Publication number
- CN115099508A CN115099508A CN202210773683.9A CN202210773683A CN115099508A CN 115099508 A CN115099508 A CN 115099508A CN 202210773683 A CN202210773683 A CN 202210773683A CN 115099508 A CN115099508 A CN 115099508A
- 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.)
- Granted
Links
- 238000005260 corrosion Methods 0.000 title claims abstract description 142
- 230000007797 corrosion Effects 0.000 title claims abstract description 139
- 238000000034 method Methods 0.000 title claims abstract description 41
- 230000001808 coupling effect Effects 0.000 title claims abstract description 17
- 238000012546 transfer Methods 0.000 claims abstract description 42
- 230000009467 reduction Effects 0.000 claims abstract description 41
- 230000007246 mechanism Effects 0.000 claims abstract description 23
- 102000018832 Cytochromes Human genes 0.000 claims abstract description 18
- 108010052832 Cytochromes Proteins 0.000 claims abstract description 18
- 230000008569 process Effects 0.000 claims abstract description 12
- 238000006722 reduction reaction Methods 0.000 claims description 54
- XEEYBQQBJWHFJM-UHFFFAOYSA-N iron Substances [Fe] XEEYBQQBJWHFJM-UHFFFAOYSA-N 0.000 claims description 53
- 238000006243 chemical reaction Methods 0.000 claims description 38
- 239000007789 gas Substances 0.000 claims description 30
- 238000009792 diffusion process Methods 0.000 claims description 20
- 230000004913 activation Effects 0.000 claims description 19
- 230000027756 respiratory electron transport chain Effects 0.000 claims description 14
- 229910052742 iron Inorganic materials 0.000 claims description 13
- 238000004364 calculation method Methods 0.000 claims description 12
- 239000012528 membrane Substances 0.000 claims description 11
- 239000002184 metal Substances 0.000 claims description 11
- 229910052751 metal Inorganic materials 0.000 claims description 11
- 238000003487 electrochemical reaction Methods 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
- 238000006703 hydration reaction Methods 0.000 claims description 6
- 230000032258 transport Effects 0.000 claims description 6
- 238000010349 cathodic reaction Methods 0.000 claims description 5
- 210000000805 cytoplasm Anatomy 0.000 claims description 5
- 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
- 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
- 238000004088 simulation Methods 0.000 claims description 3
- DOXNDFPXMBXOKH-UHFFFAOYSA-N [Au+2] Chemical compound [Au+2] DOXNDFPXMBXOKH-UHFFFAOYSA-N 0.000 claims 1
- 230000005611 electricity Effects 0.000 claims 1
- 239000007788 liquid Substances 0.000 claims 1
- 238000011156 evaluation Methods 0.000 abstract description 2
- 239000000243 solution Substances 0.000 description 6
- 241000894006 Bacteria Species 0.000 description 2
- LFQSCWFLJHTTHZ-UHFFFAOYSA-N Ethanol Chemical compound CCO LFQSCWFLJHTTHZ-UHFFFAOYSA-N 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 150000002500 ions Chemical class 0.000 description 2
- 238000010946 mechanistic model Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000003647 oxidation Effects 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
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000003181 biological factor Substances 0.000 description 1
- 239000008364 bulk solution Substances 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
- 230000001086 cytosolic effect Effects 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
- 230000009545 invasion Effects 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000002503 metabolic effect Effects 0.000 description 1
- 244000005700 microbiome Species 0.000 description 1
- 239000002070 nanowire Substances 0.000 description 1
- 238000011002 quantification Methods 0.000 description 1
- 238000011160 research 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
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)。
其中,pH为SBR细胞质中的pH值;
其中,Km为生物膜中SRB的Monod半速度系数;qmax为最大硫酸盐利用率;Xs为点蚀坑内的SRB浓度;Lf为生物膜厚度。
优选的,考虑H2CO3还原反应同时受活化反应和扩散传质控制,对H2CO3的极限扩散电流密度和电荷转移电流密度进行计算,包括如下步骤:
H2CO3还原的总阴极电流密度见式(21):
通过式(10)计算bc:
其中,αc为电子转移系数;
通过式(26)计算[CO2]b:
优选的,考虑H2O还原受活化反应控制,对H2O的电荷转移电流密度进行计算,包括如下步骤:
通过Tafel方程对H2O还原的总电流密度进行计算,计算式见式(27):
优选的,计算得到各电化学反应总的腐蚀电流密度,根据阳极和阴极反应的电荷平衡计算得到阳极腐蚀电位,并利用阳极腐蚀电位计算出不同时间的阳极腐蚀电流密度,将阳极腐蚀电流密度转换为腐蚀速率与腐蚀深度,最终计算出腐蚀速率,包括如下步骤:
通过Tafel方程计算阳极腐蚀电流密度iFe,见式(9):
通过式(12)计算io,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与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):
通过式(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方程计算还原的极限扩散电流密度。
电子从金属到生物膜的转移速率是跨生物膜的活性电荷转移层的电流密度和金属表面电势的函数。
其中,本实施例中pH代表SRB细胞质中的酸碱度,的还原发生在细胞质中,SRB细胞质中的酸碱度不同于溶液环境的酸碱度,因为受细菌代谢活动的影响,细胞质内部会保持最适宜细菌生存的酸碱度。本实施例中假设细胞质中的pH=7。
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):
其中,[CO2]b为溶解的CO2体积浓度,单位为mol/m3;为CO2在水中的扩散率,单位为m2/s;Khyd为CO2水合反应的平衡常数,取2.58×10-3,无量纲;为CO2水合反应的正向反应速率,无量纲。
通过式(26)计算[CO2]b:
5、考虑H2O还原受活化反应控制,通过Tafel方程对H2O的电荷转移电流密度进行计算;由于水分子向金属表面传输不存在扩散限制,因此可以认为H2O的还原速率仅由Tafel方程建立的电荷转移过程控制:
H2O还原的Tafel斜率与H2CO3还原的Tafel斜率相同,bc=120mV。H2O还原与H+还原具有相同的热力学条件。通过Nernst方程计算:
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 (6)
1.一种SRB与CO2耦合作用的页岩气集输管道腐蚀速率预测方法,其特征在于,包括如下步骤:
确定SRB/CO2溶液体系腐蚀过程中的电化学反应,包括CO2腐蚀的反应方程、SRB腐蚀的氧化反应方程、SRB腐蚀的还原反应方程,得出SRB与CO2腐蚀的总反应方程;
考虑H2CO3还原反应同时受活化反应和扩散传质控制,计算H2CO3还原的总电流密度;
考虑H2O还原受活化反应控制,计算H2O还原的总电流密度;
计算得到各电化学反应总的腐蚀电流密度,根据阳极和阴极反应的电荷平衡计算得到阳极腐蚀电位,并利用阳极腐蚀电位计算出不同时间的阳极腐蚀电流密度,将阳极腐蚀电流密度转换为腐蚀速率与腐蚀深度,最终计算出腐蚀速率。
3.根据权利要求2所述SRB与CO2耦合作用的页岩气集输管道腐蚀速率预测方法,其特征在于,结合胞外电子传递机制中的细胞色素传递方法,建立还原的电荷转移电流密度方程及的传质模型,计算还原的极限扩散电流密度,包括如下步骤:
其中,pH为SBR细胞质中的pH值;
其中,Km为生物膜中SRB的Monod半速度系数;qmax为最大硫酸盐利用率;Xs为点蚀坑内的SRB浓度;Lf为生物膜厚度。
4.根据权利要求3所述SRB与CO2耦合作用的页岩气集输管道腐蚀速率预测方法,其特征在于,考虑H2CO3还原反应同时受活化反应和扩散传质控制,对H2CO3的极限扩散电流密度和电荷转移电流密度进行计算,包括如下步骤:
H2CO3还原的总阴极电流密度见式(21):
通过式(10)计算bc:
其中,αc为电子转移系数;
通过式(26)计算[CO2]b:
6.根据权利要求5所述SRB与CO2耦合作用的页岩气集输管道腐蚀速率预测方法,其特征在于,计算得到各电化学反应总的腐蚀电流密度,根据阳极和阴极反应的电荷平衡计算得到阳极腐蚀电位,并利用阳极腐蚀电位计算出不同时间的阳极腐蚀电流密度,将阳极腐蚀电流密度转换为腐蚀速率与腐蚀深度,最终计算出腐蚀速率,包括如下步骤:
通过Tafel方程计算阳极腐蚀电流密度iFe,见式(9):
通过式(12)计算io,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。
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 true CN115099508A (zh) | 2022-09-23 |
CN115099508B 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) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115099151A (zh) * | 2022-07-01 | 2022-09-23 | 西南石油大学 | 基于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 |
---|
UBONG EDUOK 等: "Accelerated corrosion of pipeline steel in the presence of Desulfovibrio desulfuricans biofilm due to carbon source deprivation in CO2 saturated medium", 《MATERIALS SCIENCE AND ENGINEERING: C》, 15 August 2019 (2019-08-15), pages 1 - 8 * |
XINGYU PENG 等: "Modification of D-W model for corrosion rate of shale gas pipeline under the synergistic corrosion of SRB and CO2", 《ANTI-CORROSION METHODS AND MATERIALS》, 2 July 2021 (2021-07-02), pages 150 - 159 * |
周逸轩 等: "基于LHS-Kriging-DW的页岩气集输管道内腐蚀预测模型研究", 《 热加工工艺》, vol. 53, no. 16, 19 September 2023 (2023-09-19), pages 113 - 117 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115099151A (zh) * | 2022-07-01 | 2022-09-23 | 西南石油大学 | 基于gwo-bp修正机理模型的管道腐蚀速率预测方法 |
CN115099151B (zh) * | 2022-07-01 | 2024-06-11 | 西南石油大学 | 基于gwo-bp修正机理模型的管道腐蚀速率预测方法 |
Also Published As
Publication number | Publication date |
---|---|
CN115099508B (zh) | 2024-06-07 |
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 | |
Jeremiasse et al. | Acetate enhances startup of a H2‐producing microbial biocathode | |
Sridhar | Local corrosion chemistry—a review | |
Borole et al. | Proton transfer in microbial electrolysis cells | |
CN115099508A (zh) | Srb与co2耦合作用的页岩气集输管道腐蚀速率预测方法 | |
CN111798930B (zh) | 一种考虑腐蚀产物膜影响的co2腐蚀速率预测方法 | |
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 | |
Heßelmann et al. | Simulation‐based guidance for improving CO 2 \rmCO_2 reduction on silver gas diffusion electrodes | |
Deng et al. | Vertical galvanic corrosion of pipeline steel in simulated marine thermocline | |
Yang et al. | Development of miniature self-powered single-chamber microbial fuel cell and its response mechanism to copper ions in high and trace concentration | |
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 | |
Matsena et al. | Hexavalent chromium-reducing microbial fuel cell modeling using integrated Monod kinetics and Butler-Volmer equation | |
YG et al. | Biofilms and beyond: a comprehensive review of the impact of Sulphate Reducing Bacteria on steel corrosion | |
Zheng et al. | Investigation of cathodic reaction mechanisms of H2S corrosion using a passive SS304 rotating cylinder electrode | |
Suprapto et al. | Comparation of the analytical and experimental models of 304SS corrosion rate in 0.5 m H2SS4 with bee wax propolis extract | |
Peng et al. | Nitrous oxide production in a granule-based partial nitritation reactor: a model-based evaluation | |
Torrents et al. | Design of a microfluidic respirometer for semi-continuous amperometric short time biochemical oxygen demand (BODst) analysis | |
CN115099151B (zh) | 基于gwo-bp修正机理模型的管道腐蚀速率预测方法 | |
Żurek et al. | Morphology Analysis of the Scale Formed on Crofer 22APU Steel in Atmospheres Containing SO 2 | |
VanderSteen et al. | Mathematical modeling of the transport phenomena and the chemical/electrochemical reactions in solid oxide fuel cells: a review | |
Szydlowski et al. | Novel species identification and deep functional annotation of electrogenic biofilms, selectively enriched in a microbial fuel cell array | |
Dwivedi et al. | Tailoring flow straighteners’ architectures through additive manufacturing to enhance power performance and scale-up microbial fuel cell technology | |
Moß et al. | Optimal Geometric Parameters for 3D Electrodes in Bioelectrochemical Systems: A Systematic Approach | |
Sahin et al. | Electrochemical reduction of nitrite to ammonia for use in a bioreactor |
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 |