Indonesian Journal of Science & Technology 3 . 64-74 Yuli Panca Asmara. Tedi Kurniawan. Corrosion Prediction for Corrosion Rate of Carbon. | 64 Indonesian Journal of Science & Technology Journal homepage: http://ejournal. edu/index. php/ijost/ Corrosion Prediction for Corrosion Rate of Carbon Steel in Oil and Gas Environment: A Review Yuli Panca Asmara*. Tedi Kurniawan Faculty of Mechanical Engineering. Universiti Malaysia Pahang, 26600 Pekan. Pahang. Malaysia Correspondence: E-mail: ypanca@hotmail. Phone: 6094246239. Fax: 609424622 ABSTRACT Corrosion predictions are essential for corrosion and material It is used to prepare pre-Front End Engineering Design . re-FEED). FEED guides to select appropriate materials, planning test schedule, work over management, and estimate future repair for cost analyses. Corrosion predictions also calculate life of pipeline and equipment systems during operational As oil and gas environments are corrosive for carbon steel, it is important to account the corrosion rate of carbon steels in those environmental conditions. There are many existing corrosion predictions software, which are available in oil and gas industries. However, corrosion predictions only can be used for particular ranges of environmental conditions because they use different input parameters. To select the most applicable of corrosion predictions software, engineers have to understand theoretical background and fundamental concept of the software. This paper reviews the applications of existing corrosion prediction software in calculating corrosion rate of carbon steel in oil and gas environmental systems. The concept philosophy of software is discussed. Parameters used and range of conditions are also studied. From the results of studies, there are limitations and beneficial impacts in using corrosion software. Engineers should understand the fundamental theories of the corrosion mechanism. Knowing limitations of the models, the appropriate model can be correctly selected and interpretation of corrosion rate will close to the real data conditions. ARTICLE INFO Article History: Submitted/ Received 01 Dec 2017 First Revised 05 Jan 2018 Accepted 15 Feb 2018 Available online 09 Mar 2018 Publication Date 01 Apr 2018 ____________________ Keyword: Corrosion predictions models. Oil and gas environments. Carbon steel. A 2018 Tim Pengembang Jurnal UPI INTRODUCTION Corrosion rate in oil and gas industries involve complex factors (Kurniawan, et al. Those models basically are focused on CO2 environments (Asmara & Ismail, 2011. The additional factors such as H2S. H2S scaling effects, oil wetting effects, multiphase flow, precipitation of corrosion product films, inhibitor are, sometimes, also accounted in predictions (Asmara, et al. , 2011. 65 | Indonesian Journal of Science & Technology. Volume 3 Issue 1. April 2018 Hal 64-74 Because concepts and methods to predict corrosion rate are calculated in different ways, it is frequently found that there are differences of the corrosion data calculated by different models under nearly identical conditions. Some models provide over-estimate of corrosion rates than others(Leong, et al. Corrosion models are formulated using various scenarios among theories and experiments (Alia, et al. , 2. Those models use parameters and formulas from literatures, experimental data and their own experiences to calculate corrosion rate data (Asmara, et al. The types of corrosion models can be classified into three types. The first is mechanistic models which uses theoretical background and physical formulas to express the mechanisms of corrosion reactions. Secondly is empirical corrosion prediction model (Asmara, & Ismail, 2. It is developed based on best-fit parameter in experimental And the third is semi empirical models (Asmara, et al. , 2. This model is developed using parameters and formula from literatures and based on the researchersAo experiences. The models are also useful to manage corrosion protection strategy for the life of the metallic structural by optimizing the correct material selections. Thus, it is essential to use the appropriate corrosion model which is suitable for certain conditions in order to obtain correct corrosion rate data. So, the most cost-effective is achieved. MECHANISTIC MODELS The main concepts of mechanistic models are using electrochemical reactions and physical changes of mechanism formulas (Asmara & Ismail, 2. They includes state properties, thermodynamics theories of all It includes electrochemical reactions and diffusion process. Mechanistic model states that corrosion process covers mass transfer . , electrochemical reactions which are oxidation/reduction reactions. The model focuses on cathodic and anodic reactions which occur in the system involving several species. The mechanism of anodic dissolution depends on the dissolution rate and on the activity of hydroxide ions. While cathodic processes are related to the reduction of the species involved. Examples of mechanistic corrosion models are models derived by . e Waard & Milliams, 1975. Nesic 2. , etc. Electrochemistry Processes To investigate the corrosion mechanisms, electrochemical processes consider metal surface reactions which are transport process for each species that involves in the The model focuses in cathodic and anodic reactions which happen in the systems. The electrochemical is chemical reactions where electrons are transferred between molecules which is called oxidation/reduction reactions in anodic and cathodic site (Nordsveen, et al. , 2. Anodic Reactions In anodic reactions, there is dissolution of It is around the corrosion potential which can be under activation control or passivation process. The mechanism of activation control was proposed by (Bockris, et al. Anodic charge transfer in carbon steel is expressed as (Wang, 2. E A A Fe F E A E Fe i( F. A i exp E A iFe aOH a1H. Cathodic Reactions In metallic corrosion, cathodic processes are due to the reduction of separate species, . H ,H2O. H2CO3 and HCO3-. In acidic solutions, the reduction of H is the dominant cathodic reaction. There are two possibilities DOI: http://dx. org/10. 17509/ijost. p- ISSN 2528-1410 e- ISSN 2527-8045 Yuli Panca Asmara. Tedi Kurniawan. Corrosion Prediction for Corrosion Rate of Carbon. | 66 reactions in cathodic site, which are diffusion limiting current density and activation current Activations current density are given by expression (Wang, 2. E A A H F E A E Ho E i H ,a A i exp E The exchange current density is given by iHo A iH* aH0. 5 aH2. where aH is activation of hydrogen ions and aH2O is activation of H2O. The limiting current density results from diffusion-limited transport of protons to the metal surface and can be calculated as iH ,lim A K m FaH where km is the mass transfer coefficient. The value of km can be calculated if the flow regime, diffusion coefficient of H ions and solution viscosity are known. Scale formation In the case of corrosion product on the metal surface with film formation, it can be expressed mathematically. These corrosion mechanisms are based on several assumptions which can be described as follows: convective diffusion, molecular diffusion, and diffusion via solid film. Corrosion mechanism which happens in solutions as a combination of mix gases can be expressed from the following equation. This expression is the case of corrosion rate of steel due to mixed species (Wang, 2. CRi A Ai e RTK E A 1 EE cb,i A CRi EE oc A E Au Ao m ,i E E i c s ,i . where Km is the mass transfer coefficient of species i . Cb is the bulk concentration of species i . ol/m. Co is the interfacial concentration of species i at outer scale/solution interface . ol/m. Di is the diffusion coefficient for dissolved species i . Au is the outer scale porosity. A is tortuosity factor. Ci is interfacial concentration of species i. Aos is the thickness of outer film scale. Ahbl is the turbulence boundary layer thickness. Ambl is the mass transfer boundary layer thickness. Af is the film thickness. A is the Arhrhenius constants. Tk is the temperature (Kelvi. , and cs is the surface concentration. CORROSION PREDICTIONS SOFTWARE Many industrial corrosion companies produce corrosion predictions software, which are based on empirical models and combination of empirical and mechanistic approach. Empirical corrosion prediction models are developed based on best-fit parameter in experimental regression so called Semiempirical models. The semi empirical models are developed using parameters and formula from literatures combined with their experimental data. In oil and gas industries, there are many corrosion productions software. These include Norsok. Cassandra. Electrochemical Corrosion Engineering (ECE). Ohio Model etc. They develop models by using their own experiments, field data combined with fundamental formulas and its many subsequent derivatives. All of these were developed based on different systems and assumptions. DOI: http://dx. org/10. 17509/ijost. p- ISSN 2528-1410 e- ISSN 2527-8045 67 | Indonesian Journal of Science & Technology. Volume 3 Issue 1. April 2018 Hal 64-74 NORSOK The NORSOK standard is owned by Norwegian Oil Industry Association and Federation of Norwegian Manufacturing Industries. The program covers only calculation of corrosion rates where CO2 is the corrosive agent. does not include the corrosivity, e. contamination of O2. H2S etc. The model is an empirical model for CO2 at different temperatures, pH. CO2 fugacities, wall shear stresses, and temperatures from 20 to160AC. The model calculates pH and wall shear stress. The effect of acetic acid is not account for in this model, but it is still valid to predict corrosion rate if the concentration of acetic acid is less than 100 ppm. ECE (Electronic Chemical Engineerin. ECE program software calculates corrosion rate based on the modified model by . e Waard & Milliams 1. ECE model includes oil wetting correlation based on field correlation. ECE proposes a corrosion prediction expression using corrosion reactions and mass transfer effects. The mass transfer represents the main part of the dependence on flow velocity and pipe diameter. In ECE, they develop corrosion prediction by involving several variables such as gas fugacity, formation of protective films, effect of ferrous ions, presence of oil, effect of condensing water, and effect of multiple phase. Cassandra (DWM . Cassandra is a model developed based on the experiences of . e Waard & Milliams The input includes pH. CO2 concentration, temperature, and water contaminant. This model does not consider scaling temperature. The user must set an assumption of the scaling temperature. This model has important aspects which influence rate of corrosion, namely corrosion inhibitor availability and corrosion risk categories. The model also accounts for the presence of acetate in water as acetic acid. The major input to the model are: CO2 mole %, temperature, total pressure, liquid velocity and water chemistry. Besides that, the model has secondary input, such as hydraulic diameter and glycol concentration, oil type . rude or condensat. and water type . ondensed water or formation wate. The effect of oil wetting in this model is not included. Parameters Comparison of NORSOK and CASSANDRA Models Both NORSOK and Cassandra have limitations in predict the corrosion rate. Table 1 shows the parameter comparison between NORSOK and Cassandra models. Experimental Corrosion Predictions Models in Oil and Gas Environments The corrosion mechanism of mild steel in the presence of CO2 in various conditions has been a widely reviewed by researchers. Tthere are many experiments and field studies have been conducted. The initial studies in CO2 corrosion was conducted by . e Waard & Milliams 1. that have become a fundamental concept for the further studies on the CO2 corrosion phenomenon. The newest of corrosion mechanism was studied by (Nesic & coworkers 2. who have claimed to be successful in modeling CO2 corrosion rate based on theoretically. Table 2 compiles corrosion predictions formulas based on experimental Effect of Flow in Corrosion Modelling Effect of flow on steel corrosion is a type of corrosion caused by a combination between mechanical and electrochemical effects (Silverman, 1. Mechanical effects due to water motion causes impingement that leads to metal removal and material DOI: http://dx. org/10. 17509/ijost. p- ISSN 2528-1410 e- ISSN 2527-8045 Yuli Panca Asmara. Tedi Kurniawan. Corrosion Prediction for Corrosion Rate of Carbon. | 68 Water that flows to the surface can wear the corrosion product film or create shear stress to the surface. Thus, corrosion will occur faster (Eisenberg et al. , 1. Corrosion rate also can increase due to effects of differences in velocity turbulence across the Parameters that influence flow induced corrosion are hydrodynamic boundary layer and rate of momentum transfer from the bulk to the wall. In this conditiion, corrosion may be controlled by the rate of mass transfer of a reactant or the rate of corrosion products. As calculated by two models, they are confirmed that increasing flow rate corrosion rate will also increase. However, empirical model calculates corrosion rate higher than freecorp model. Empirical model states that increasing flow, corrosion rate continues increase, while freecorp finds that the highest corrosion rate occurs at 1000 rotation speed. FreeCorp indicates that when rotation speeds are more than 1000 rpm, corrosion rate will remain constant which refers to limiting current density (Silverman, 2. Table 1. Comparison Parameters considered in some corrosion predictions software (NORSOK. ECE. FreeCorp, and Cas-sandr. Parameters Temperature (C) Total Pressure . Total mass flow CO2 fugacity . Wall Shear Stress . Glycol concentration Inhibitor efficiency (%) Diameter . Liquid velocity . H2 S Pipeline orientation Polarization graph Formation water/condensed water Types of materials Cost analyses Length of pipelines Oil contents Corrosion prediction software ECE N/A N/A Norsok N/A N/A N/A N/A N/A N/A N/A N/A DOI: http://dx. org/10. 17509/ijost. p- ISSN 2528-1410 e- ISSN 2527-8045 Cassandra N/A N/A N/A N/A N/A N/A N/A N/A N/A N/A N/A N/A N/A N/A N/A FreeCorp N/A N/A N/A N/A N/A En 69 | Indonesian Journal of Science & Technology. Volume 3 Issue 1. April 2018 Hal 64-74 Table 2. Corrosion formulas calculated using some experimental data CR= eOe 9225/T x10OepH x5. For worst case calculation of the corrosion rate (CR) of low alloy steel in slightly sour conditions OIycE yaycI = 31. 15 yaycyceycyceyc yaycycaycycyccyce ( ) 10. 3 ycO 0. 6 ycEyaycC2 ycNyce Oe2671/ycN ya Crude oil type and water cut (W<30%) OIycE yaycI = 31. 15 yaycyceycyceyc yaycycaycycyccyce ( ya ) 10. 3 ycO1. 6 ycEyaycC2 ycNyce Oe2671/ycN Crude oil type and water cut (W>30%) CR is corrosion rate . m/y. , iP/L is the pressure gradient (N/m. W is water cut. PCO2 is partial pressure of carbon dioxide (MP. T is temperature (K). Crcrude is factor to account for crude oil type, based on . < Crcrude <. Crfreq (= 0. 023 (F) 0. is the normalized factor to account for slug frequency . < Crfreq <. , and F is slug frequency to a maximum of 35 slugs/minute. Vcor = 31. 15 Crfreq Croil . 38 F. 3 wc0. 6 PCO20. 8 T exp(-2671/T) where Crfreq is term for slug frequency. Croil is term for crude oil type, wc is water fraction ( %). T is temperature in K. Fr is Froude number in the liquid film defined as: Fr is (Vt Ae V. / . with Vt is translational velocity of the slug . Vf is average velocity of the stratified liquid film . G is acceleration due to gravity. , and hf is effective height of the liquid film. CR = 8. 856x10-3 (O. 48x10-7(O. - 1. H) 4. 934x10-2(CO. (H2S)4. 8231x10-5(CO. (O. 372x10-3(H2S)(O. 113x10-3(O. H) where CR is general corrosion rate, (O. is O2 concentration of gas . , (CO. is CO2 concentration of gas . H2S is H2S concentration of gas . , pH = pH initial OIycE 3. ) 1 ycO ycEyaycC2 ycNyceycuycyOe2671/ycN ya CR is general corrosion rate. FF is the frequency factor, (DP/L)0. 33 is the flow pressure gradient in Pascal/meter, (PCO. 8 is the CO2 partial pressure in MPascal. T is the Kelvin temperature. yaycI = 31yaya ( Effects of Temperature on Corrosion Rate of Carbon Steel In oil and gas environments (CO2 syste. , temperature affects the conditions for formation of the protective carbonate layers and affects corrosion rate in a different manner. At temperatures lower than 60 AC, the solubility of FeCO 3 is high and the precipitation rate is slow. thus protective films will not form until the pH is increased more than solubility product (Nafday. Above 60AC the solubility of FeCO 3 decreases and the protectiveness of the iron carbonate layer increases with temperature. thus, the corrosion rate is reduced. Scaling temperature is the temperature where corrosion rate reaches a maximum value. DOI: http://dx. org/10. 17509/ijost. p- ISSN 2528-1410 e- ISSN 2527-8045 Yuli Panca Asmara. Tedi Kurniawan. Corrosion Prediction for Corrosion Rate of Carbon. | 70 Corrosion rate . FreeCorp Emp. Model Rotation speed (RPM) Figure 1. Effects of flow on corrosion rate as calculated by empirical model and FreeCorp model. Corrosion rate . Figure 2 shows corrosion rate of carbon steel in CO2 environments at 1 bar and pH 5. Models of corrosion prediction provide different values. The most conservative is obtained by Cassandra model. And the lowest corrosion rate is calculated using norsok software. From the figure, it is clearly that all of software shows increasing of corrosion rate when temperature increases. Only Cassandra states that corrosion rate increases until temperature 60 oC. It tends to decrease when temperature exceeds 60 oC. The temperature is called as scaling temperature. From the figure, it also shows that there are obvious differences among the models. It means that understanding effects of temperature on corrosion rate are still under debatable. Every models uses different approach to count corrosion rate. The differences come from different assumption regarding effects of film formation, quality of film, solubility of FeCO 3, saturation of pH, water cut, interaction among species (Silverman, 2. DW93 Model Cassan Temperature . K) Figure 2 Comparison effects of temperature on corrosion rate of carbon steel at pH5 and 1 bar as calculated by some corrosion models. EFFECTS OF pH pH is an important parameter for corrosion process. Normally, corrosion rate will be lower at higher pH,. The pH of a wet gas in CO2 system is often assumed as equal to the saturation of FeCO3 precipitation. Figure 3 shows the predicted corrosion rates as effects of pH which differs depending on the models. DOI: http://dx. org/10. 17509/ijost. p- ISSN 2528-1410 e- ISSN 2527-8045 71 | Indonesian Journal of Science & Technology. Volume 3 Issue 1. April 2018 Hal 64-74 DW93 Cassandra DW95 MODEL Norsok Corrosion rate . Figure 3 Comparison effects of pH on corrosion rate of carbon steel at pH5 and 1 bar as calculated by some corrosion models. Design of Experiment (DOE) and Statistical Modeling for CO2 System Recently, empirical models of corrosion processes have been used to predict corrosion process involving several independent However, most of the empirical models do not predict the corrosion rate in several variables simultaneously (Mokhtar. Using empirical methods, modeling interactions effects among the species and the operational conditions simultaneously, require large number of experiment which is costly and takes time. These limitations can be overcome by using design experiment of response surface methodology (RSM). This is a simple method and believed can represent overall unselected variables. Contour P lot of Cor. m/ . vs Temperature . C), pH Cor. < 3. 0 Ae 3. 5 Ae 4. 0 Ae 4. 5 Ae 5. 0 Ae 5. 5 Ae 6. > 6. Temperat ure . From the figure, corrosion predictions software present large difference when pH is But at pH 6, almost all models give the same values. These are because various assumptions in calculating corrosion rate as effects of pH. The pH in CO2 system is often assumed as equal to the saturation of FeCO3 The pH can also be calculated by involving the concentration of species such as CO2. H2CO3. HCO3-. CO32-. H2S. HS-. S2. H . OH-. H2O. Fe2 . CH3COOH . cetic aci. CH3COO-. Sometimes, pH calculation can be done by incorporating the FeCO3 precipitation kinetics. The pH is also influenced by H ions concentration, temperature, pressure, and ionic strength. Thus pH in calculation can be different among the models software. Dissolved iron bicarbonate as the initial corrosion product will also contribute to increase the pH of solution. pH = 4. Temperature . C) = 56. Cor. = 6. Hold Values HAc . Rot. 3500 Figure 4. Response surface graph for corrosion rate as a function of temperature and pH (HAc at 30 ppm and rotation speed at 3500 rp. The application of RSM allows visualization of the experimental results in a 3-D DOI: http://dx. org/10. 17509/ijost. p- ISSN 2528-1410 e- ISSN 2527-8045 Yuli Panca Asmara. Tedi Kurniawan. Corrosion Prediction for Corrosion Rate of Carbon. | 72 RSM is used to determine optimal levels for variables input. RSM is a sequential procedure for constructing empirical relation for the experimental data. Using response information, the optimum data between factors can be developed and model improvements can be achieved. It has been proven that researchers have used response surface method (RSM) to process data systematically that can allow to apply multiple regression Response surface design methodology is also often used to refine models to obtain an optimum design. RSM is useful to obtain critical points in the experimental The surfaces generated by linear or polynomial models can be used to indicate the direction in which the original design must be started to attain the optimal conditions. For polynomial models, the critical point can be characterized as maximum, minimum, or saddle. Using RSM, it is possible to calculate the coordinates of the critical point through the first derivative of the mathematical function (Asmara, et al. , 2. First derivative equals to zero indicates that critical points is located. They have studied corrosion rate in CO2 environments by using RSM successfully (Mune, et al. , 2. They are all claim that RSM can reduce number of experiments with satisfied results. CONCLUSION There are differences of corrosion rate predicted by various corrosion models. Variations of results are caused by assumptions made by corrosion models. Particular parameters such as H2S concentration, scale effects, effects of wall shear stress and hydrodynamic condition of the solutions conditions will impact on results. ECE and FreeCorp are the models with more parameters inputs. These two models are flexible which can be applied for any environmental conditions. Other model like Cassandra gives more conservative which contribute a greater over design. Empirical methods combined with RSM propose improvements techniques. RSM can overcome limitation of pure empirical methods by simplifying experiments models. Selecting the best corrosion models require further interpretation to explain real conditions, thus mechanistic methods are more realistic. The user should also understand comparison regarding limitation, advantages and scope of assumptions of the models to obtain appropriate data. ACKNOWLEDGMENTS The authors are thankful to Universiti Teknologi Petronas and Universiti Malaysia Pahang for providing grant and facilities for the research. AUTHORSAo NOTE The author. that there is no conflict of interest regarding the publication of this article. Authors confirmed that the data and the paper are free of plagiarism. REFERENCES