Phast validation of discharge and atmospheric dispersion for pressurised carbon dioxide releases

The consequence modelling package Phast examines the progress of a potential incident from the initial release to the far-field dispersion including the modelling of rainout and subsequent vaporisation. The original Phast discharge and dispersion models allow the released substance to occur only in the vapour and liquid phases. The latest versions of Phast include extended models which also allow for the occurrence of fluid to solid transition for carbon dioxide (CO₂) releases.As part of two projects funded by BP and Shell (made publicly available via CO2PIPETRANS JIP), experimental work on CO₂ releases was carried out at the Spadeadam site (UK) by GL Noble Denton. These experiments included both high-pressure steady-state and time-varying cold releases (liquid storage) and high-pressure time-varying supercritical hot releases (vapour storage). The CO₂ was stored in a vessel with attached pipework. At the end of the pipework a nozzle was attached, where the nozzle diameter was varied.This paper discusses the validation of Phast against the above experiments. The flow rate was predicted accurately by the Phast discharge models (within 10%; considered within the accuracy at which the BP experimental data were measured), and the concentrations were found to be predicted accurately (well within a factor of two) by the Phast dispersion model (UDM). This validation was carried out with no fitting whatsoever of the Phast extended discharge and dispersion models.     

Modelling of discharge and atmospheric dispersion for carbon dioxide releases

This paper discusses the modelling of the discharge and subsequent atmospheric dispersion for carbon dioxide releases using extensions of models in the consequence modelling package Phast. Phast examines the progress of a potential incident from the initial release to the far-field dispersion including the modelling of rainout and subsequent vaporisation. The original Phast discharge and dispersion models allow the released chemical to occur only in the vapour and liquid phases. As part of the current work these models have been extended to also allow for the occurrence of liquid to solid transition or vapour to solid transition. This applies both for the post-expansion state in the discharge model, as well as for the thermodynamic calculations by the dispersion model. Solid property calculations have been added where necessary. The above extensions are generally valid for fluid releases including CO₂. Using the extended dispersion formulation, a sensitivity study has been carried out for mixing of solid CO₂ with air, and it is demonstrated that solid effects may significantly affect the predicted concentrations.     

Modelling three-phase releases of carbon dioxide from high-pressure pipelines

This paper describes the development and experimental validation of a three-phase flow model for predicting the transient outflow following the failure of pressurised CO₂ pipelines and vessels. The choked flow parameters at the rupture plane, spanning the dense-phase and saturated conditions to below the triple point, are modelled by maximisation of the mass flowrate with respect to pressure and solids mass fraction at the triple point. The pertinent solid/vapour/liquid phase equilibrium data are predicted using an extended Peng–Robinson equation of state.The proposed outflow model is successfully validated against experimental data obtained from high-pressure CO₂ releases performed as part of the FP7 CO2PipeHaz project (www.co2pipehaz.eu).The formation of solid phase CO₂ at the triple point is marked by a stabilisation in pressure as confirmed by both theory and experimental observation. For a fixed diameter hypothetical pipeline at 100 bar and 20 °C, the flow model is used to determine the impact of the pipeline length on the time taken to commence solid CO₂ discharge following its rupture.     

An alternative CFD tool for gas dispersion modelling of heavy gas

The numerical simulation of gas dispersion is of great importance in various areas of engineering such as optimisation, synthesis of chemical process, petroleum industry and process safety. The OpenFOAM (Open Field Operation and Manipulation) code is a free and open source computational fluid dynamics (CFD) program. The current research is focused on the development and customisation of a computational tool for handling gas dispersion of heavy gases, such a LNG and CO₂. The novel CFD tool relies on OpenFOAM framework. The core of the work is based on the OpenFOAM solver rhoReactingBuoyantFoam to handle gas dispersion. A series of CFD simulations has been performed for methane and CO₂. The source term of the former is modelled by HSM (Hybrid Switch Model). The model comprises contribution from HEM (Homogeneous Equilibrium Model) approach, frozen model and non-equilibrium model for CO₂ leak. The novel approach switches between equilibrium and non-equilibrium conditions based on the meta-stable parameter on the grounds of thermodynamics and experimental observations. Good agreement with experimental data is observed. Numerical findings for methane leakage from the proposed CFD tool are compared with experimental data and FLACS. Good agreement is observed.     

The influence of impurities on the transportation safety of an anthropogenic CO2 pipeline

Transportation safety is a key aspect of carbon capture and storage (CCS), which is a major technology used to reduce greenhouse gas emissions. Supercritical CO₂ pipelines have been certified as an optimised choice for CO₂ transportation. The results of this study show that the Peng–Robinson (PR) equation of state is recommended for analysis of the properties of supercritical CO₂. The influence of nonpolar and polar impurities on the two-phase region and the location of the sharp discontinuity in the density are found by analysing the ternary phase equilibrium and physical parameters using the PR equation of state. A transitional area between the supercritical phase and the dense phase, where the density changes abruptly, is defined as the quasi-critical region. This study describes the functional relation between the temperature and the pressure that defines the quasi-critical line by calculating the partial derivative equations and then determines the effect of impurities on the quasi-critical region of transported CO₂. Operational recommendations for pipeline transportation of flue CO₂ are developed using a pipeline operated by Sinopec as an example, demonstrating the influence of impurities in flue CO₂ on saturation pressure for control and prevention of fractures in CO₂ pipelines.     

Validation of discharge and atmospheric dispersion for unpressurised and pressurised carbon dioxide releases

This paper discusses the validation of discharge and subsequent atmospheric dispersion for both unpressurised and pressurised carbon dioxide releases using the consequence modelling package Phast.The paper first summarises the validation of the Phast dispersion model (UDM) for unpressurised releases. This includes heavy gas dispersion from either a ground-level line source (McQuaid wind-tunnel experiments) or an area source (Kit-Fox field experiments). For the McQuaid experiments minor modifications of the UDM were made to support line sources. For the Kit Fox experiments steady-state and 20-s finite-duration releases were simulated for both neutral and stable conditions. Most accurate predictions of the concentrations for finite duration releases were obtained using the UDM Finite Duration Correction method.Using experiments funded by BP and Shell and made available via DNV's CO2PIPETRANS JIP, the paper secondly summarises the validation of the Phast discharge and dispersion models for pressurised CO₂ releases. This modelling accounted for the possible presence of the solid CO₂ phase following expansion to atmospheric pressure. These experiments included both high-pressure steady-state and time-varying cold releases (liquid storage) and high-pressure time-varying supercritical hot releases. Both the flow rate and the concentrations were found to be predicted accurately.The above validation was carried out with no fitting whatsoever of the Phast extended discharge and dispersion models.     

Fracture criterion and control plan on CO2 pipelines: Theory analysis and full-bore rupture (FBR) experimental study

Pressurized pipelines are the most reliable and cost-effective option for the long-distance transportation of CO₂ from an emitter to an onshore storage site. Propagating or unstable factures are considered catastrophic pipeline failures, resulting in a massive escape of inventory within a short period of time. The decompression curve for CO₂ exhibits a large drop in decompression wave speed at the phase transition pressure, leading to a higher driving force for crack propagation. The study of fracture control plans is very important for assessing the possibility of fracture propagation and preventing unstable fracturing along CO₂ pipelines. Three full-bore rupture (FBR) experiments were performed using an industrial-scale (258 m long, 233 mm inner diameter) CO2 pipeline with initial CO2 states of gaseous, dense and supercritical phases, respectively. The relation between the decompression velocity and the pipeline fracture propagation velocity was analyzed during the process of buried CO2 pipeline release. A fracture propagation criterion was established for the buried CO₂ pipeline. For the gaseous CO₂ leakage, the pressure plateau corresponding to the decompression wave velocity only appeared near the closed end of the pipeline. For the dense CO₂ leakage, the pressure plateau corresponding to the decompression wave velocity was observed near the saturation pressure after rapid decompression. For the supercritical CO₂ leakage, the pressure plateau corresponding to the decompression wave velocity was observed in the stage when the supercritical CO₂ transformed into the two phases of gas and liquid. Compared with the gaseous and dense CO₂, for the supercritical CO_2, the initial decompression wave velocity was the smallest, and the requirement of the pipeline safety factor was the highest.     

Thermodynamic method for the prediction of solid CO2 formation from multicomponent mixtures

The increase in GHG concentration has a direct effect on global climate conditions. Among the possible technologies to mitigate GHG emissions, CCS is being accepted to gain emission reduction. Such technology also involves cryogenic CO₂ capture processes based on CO₂ freeze-out or where the formation of solid CO₂ must be avoided. Captured CO₂ is usually transported in pipelines for the reinjection.The risk associated to the release of CO₂ is due to the changing temperatures and pressures the system may experience, which can lead to the deposition of solid CO₂ where it must be avoided. Prolonged exposure to dry ice can cause severe skin damage and its resublimation could pose a danger of hypercapnia. It is, thus, necessary to build up a tool able to predict the conditions in which CO₂ can freeze-out.A thermodynamic methodology based on cubic EoSs has been developed which is able to predict solid–liquid–vapor equilibrium of CO₂ mixtures with n-alkanes or H₂S which are usually found in equipment for acidic gas, mainly natural gas, treatment.The focus is a detailed analysis of the method performances when more than two components are present since, for such a case, literature does not provide significant modeling results.     

Evaluation of multi-phase atmospheric dispersion models for application to Carbon Capture and Storage

A dispersion model validation study is presented for atmospheric releases of dense-phase carbon dioxide (CO₂). Predictions from an integral model and two different Computational Fluid Dynamics (CFD) models are compared to data from field-scale experiments conducted by INERIS, as part of the EU-funded CO₂PipeHaz project.The experiments studied consist of a 2 m³ vessel fitted with a short pipe, from which CO₂ was discharged into the atmosphere through either a 6 mm or 25 mm diameter orifice. Comparisons are made to measured temperatures and concentrations in the multi-phase CO₂ jets.The integral dispersion model tested is DNV Phast and the two CFD models are ANSYS-CFX and a research and development version of FLACS, both of which adopt a Lagrangian particle-tracking approach to simulate the sublimating solid CO₂ particles in the jet. Source conditions for the CFD models are taken from a sophisticated near-field CFD model developed by the University of Leeds that simulates the multi-phase, compressible flow in the expansion region of the CO₂ jet, close to the orifice.Overall, the predicted concentrations from the various models are found to be in reasonable agreement with the measurements, but generally in poorer agreement than has been reported previously for similar dispersion models in other dense-phase CO₂ release experiments. The ANSYS-CFX model is shown to be sensitive to the way in which the source conditions are prescribed, while FLACS shows some sensitivity to the solid CO₂ particle size. Difficulties in interpreting the results from one of the tests, which featured some time-varying phenomena, are also discussed.The study provides useful insight into the coupling of near- and far-field dispersion models, and the strengths and weaknesses of different modelling approaches. These findings contribute to the assessment of potential hazards presented by Carbon Capture and Storage (CCS) infrastructure.     

Release and dispersion behaviour of carbon dioxide released from a small-scale underground pipeline

The development of carbon capture and storage (CCS) brings challenges for safety issues regarding carbon dioxide (CO₂) transmission pipelines. Once a pipeline is punctured or full-bore ruptured, the leaked CO₂ is hazardous to personnel and the environment. Small-scale devices were established with the aim of studying the release and dispersion behaviour of gas and liquid CO₂ from a punctured underground pipeline. A sandbox was built to simulate the underground conditions. The parameters of the sand used in the experiments were tested. CO₂ concentrations on the ground and temperatures around the release orifice in the sand were analysed. The results indicate that in the CO₂ gas release experiments, the CO₂ concentration on the sand surface decreases with increasing horizontal distance in the form of a power function. CO₂ concentrations in upward release are slightly larger than those in horizontal release at the same location but are obviously bigger than values in downward release. The temperature-drop region is much smaller than that in air. A frozen ice ball can be generated near the release orifice during the gas phase of the CO₂-release process. In the liquid phase of CO₂-release experiments, a large amount of dry ice is generated near the release orifice. Dry ice can only be generated in the area close to the release orifice, especially in the near-field area.     

Experimental investigation of the flow characteristics of jettisoning in a CO2 carrier

This experimental study was performed to investigate the flow characteristics in the jettisoning flow line of a liquid CO₂ carrier. When a pressurized liquid CO₂ container loses mechanical integrity, possibly by material or mechanical defects, the liquid inventory should be drained out rapidly for safety reasons using the so-called jettisoning process. In the course of jettisoning, the liquid CO₂ undergoes two phase change stages, from liquid to liquid + vapor and from liquid + vapor to solid + vapor. Consequently, the jettisoning release rate is affected by the characteristics of these phase changes. In this study, liquid CO₂ was discharged through a small tube, representing a jettisoning flow line. The temperature and pressure were measured along the tube, and the locations of the phase changes were inferred from the measured data. The experimental results showed that active nucleation occurred near the tube tip and that the phase change into solid and vapor occurred just after leaving the pipe, irrespective of the tube length in this study.     

A risk assessment methodology for high pressure CO2 pipelines using integral consequence modelling

This paper presents a risk assessment methodology for high pressure CO₂ pipelines developed at the Health and Safety Laboratory (HSL) as part of the EU FP7 CO2Pipehaz project. Until recently, risk assessment of dense phase and supercritical CO₂ pipelines has been problematic because of the lack of suitable source term and integral consequence models that handle the complex behaviour of CO₂ appropriately. The risk assessment presented uses Phast, a commercially available source term and dispersion model that has been recently updated to handle the effects of solid CO₂. A test case pipeline was input to Phast and dispersion footprints to different levels of harm (dangerous toxic load and probit values) were obtained for a set of pipeline specific scenarios. HSL's risk assessment tool QuickRisk was then used to calculate the individual and societal risk surrounding the pipeline. Knowledge gaps that were encountered such as: harm criteria, failure rates and release scenarios were identified and are discussed.     

Simulation study on near-field structure and flow characteristics of high-pressure CO2 released from the pipeline

The accidental release of high-pressure carbon dioxide (CO₂) can cause serious damages to both humans and pipeline equipment. Therefore, it is of great significance to have a deeper understanding about the release characteristics of high-pressure CO₂ for improving the safety level of Carbon Capture and Storage (CCS) technologies. Both industrial-scale and laboratory-scale studies have been carried out to predict the release behaviors. In recent years, computational fluid dynamics (CFD) simulation has become a crucial method to study the instantaneous changes and microscopic details of the fluid behaviors. In this paper, the simulation method was employed to study the near-field structure and flow characteristics of high-pressure CO₂ released from pipelines. The Peng-Robinson Equation of State (EOS) was used to compute the thermodynamic properties of high-pressure CO₂, and SST k-ω model was applied to simulate the structure and physical parameters of the under-expanded jet. In addition, the multi-phase mixture model was introduced to study the phase transition. The non-equilibrium liquid/vapor transition is modeled by introducing ‘source terms’ for mass transfer and latent heat. Compared to the experimental results, the simulation results showed good agreement. Furthermore, the influences of operating conditions, including different stagnation pressure, stagnation temperature, and nozzle size, were analyzed.     

Generalising two-phase homogeneous equilibrium pipeline and jet models to the case of carbon dioxide

The safety of high pressure pipelines is nowadays routinely analysed using simple integral models to predict the nature of a release following a possible breach. Where two-phase flow may result, models making the assumption of homogeneous equilibrium between liquid and gas phases have been used successfully. In the case of a pipeline carrying carbon dioxide, however, solidification of a significant amount of the gas (owing to its high triple point pressure) is a distinct possibility, necessitating consideration of gas, liquid, and solid phases. This paper shows how integral models of two-phase flow can be generalised to accommodate the solid phase, and investigates some of the details of the flow where the triple point pressure occurs and partial solidification begins. It is shown that two-phase homogeneous equilibrium flow models may be generalised straightforwardly to cover the case of carbon dioxide, and some features are seen to emerge independently of more detailed flow modelling considerations.     

Inherently Safer Design for low-temperature process during depressurization: An industrial case study

The exploitation of a low-quality gas field with high CO₂ concentration is more viable through liquid CO₂ produced from cryogenic distillation technology. Despite the bright potential of the technology, there are deficiencies in handling high concentration of CO₂ at low temperature and high pressure during the blowdown condition. This study focuses on the CO₂ blowdown at a cryogenic pilot plant designed to manage high concentrations of CO₂ in the feed gas, high pressures, and low temperatures. A comprehensive design review and risk assessment using Inherent Safer Design (ISD) indexes were carried out in this study. The ISD was performed to identify the current risk level, and the critical parameters that may cause solid CO₂ formation in the piping or equipment as well as to identify mitigation measures to avoid the temperature to drop below the CO₂ freezing point during blowdown. The present findings confirmed that the initial pressure and temperature, as well as CO₂ concentration are key parameters towards significant impact on blowdown conditions. Therefore, the reduction of the feed gas pressure from 80 bar to 70 bars has minimized the Joule Thomson (JT) effect during blowdown and avoided the CO₂ solid formation in the system. Moreover, the relocation of the blowdown valve at the downstream heater resulted in a higher final temperature above the CO₂ freezing point. The ISD indexes confirmed that the cryogenic facilities are inherently safer during blowdown with the mitigation measures adopted.     

The extreme carbon dioxide outburst at the Menzengraben potash mine 7 July 1953

Carbon dioxide is an asphyxiant and an irritant gas. An extreme outburst of carbon dioxide took place 7 July 1953 in a potash mine in the former East Germany. During 25 min, a large amount of CO₂ was blown out of the mine shaft with great force. It was wind still and concentrated CO₂ accumulated in a valley leading to multiple asphyxiation casualties. Based on a review of concentration-response relationships, the location of victims, and other information, it is concluded that concentrations of 10-30% carbon dioxide may have occurred 450 m from the point of release for at least 45 min. It is concluded that 1100-3900 tonnes of CO₂ were blown out of the mine shaft, possibly with intensities around 4 tonnes/s. It is also concluded that the large majority of the gas escaped as a near-vertical high-velocity jet with only little loss of momentum due to impingement. The release was modelled using PHAST. Output from the model is inconsistent with the asphyxiation harm observed. The high-momentum release is predicted to disperse safely and never reach the ground. Carbon dioxide capture and storage (CCS) schemes will involve handling and transportation of unprecedented quantities of CO₂. Case histories to date include sudden releases of CO₂ of up to 50 tonnes only, far too small to provide a suitable empirical perspective on predicted hazard distances for CCS projects. The 1953 outburst contributes to filling this gap.     

Experimental study on CO and CO2 emissions from spontaneous heating of coals at varying temperatures and O2 concentrations

Laboratory experiments were conducted to investigate carbon monoxide (CO) and carbon dioxide (CO₂) emissions from spontaneous heating of three U.S. coal samples in an isothermal oven at temperatures between 50 and 110 °C. The oxygen (O₂) concentration of an oxygen/nitrogen (N₂) mixture flowing through the coal sample was 3, 5, 10, 15, and 21%, respectively. The temperature at the center of the coal sample was continuously monitored, while the CO, CO₂, and O₂ concentrations of the exit gas were continuously measured. The results indicate that the CO and CO₂ concentrations and the CO/CO₂ ratio increased when the initial temperature was increased. As the inlet O₂ concentration increased, the CO and CO₂ concentrations increased, while the CO/CO₂ ratios tended to converge to the same value. The ratio of CO/CO₂ was found to be independent of coal properties, approaching a constant value of 0.2. The maximum CO production rate correlated well with the maximum coal temperature rise. The apparent order of reaction for coal oxidation was estimated to be between 0.52 and 0.72. The experimental results in this study could be used for early detection and evaluation of a spontaneous heating in underground coal mines.     

Experimental research on displacing coal bed methane with supercritical CO2

The high-gas and low-permeability are the common problems of China coal mine, which restrain the mining of coal-seam gas resources safely and efficiently. Hence, to solve the problem of low permeability of coal seam, an experimental system was set up and experimental research was conducted to investigate the effect of the displacement of methane by injecting supercritical CO₂ into coal samples. The experimental results indicated that, the extraction effect of supercritical CO₂ changes the coal’s porosity, and broadens the seepage channel for methane. Thus, the methane could be desorbed effectively from the coal matrix, and flow through more cracks at higher speed.     

Reducing the risk level for pipelines transporting carbon dioxide and hydrogen by means of optimal safety valves spacing

In many countries where electricity generation is based on their natural resources of fossil fuels a need arises to implement new power engineering technologies that allow carbon dioxide capture. Simultaneously, efforts are made to find new energy carriers which, if fired, do not involve carbon dioxide emissions. Hydrogen is one of such fuels with this future potential which is now becoming increasingly popular. Obviously, this means that the two gases mentioned above – carbon dioxide and hydrogen – will be produced in large quantities in future, which in many cases will necessitate their transport over considerable distances. If a pipeline failure occurs, the transport of the gases may pose a serious hazard to people in the immediate vicinity of the leakage site. This paper presents an analysis of the possibility of reducing the level of risk related to pipelines transporting CO₂ and H₂ by means of safety valves. It is shown that for a 50 km long and a 0.4 m diameter pipeline transporting gas with the pressure of 15 MPa the individual risk level can be reduced from 1·10⁻⁴ to 6.5·10⁻⁷ for CO₂ and from 1·10⁻⁶ to 6·10⁻¹⁰ for H₂. The social risk can be diminished in similar proportions.     

Risk assessment in a hypothetical network pipeline in UK transporting carbon dioxide

With the advent of Carbon Capture and Storage technology (CCS) the scale and extent of its handling is set to increase. Carbon dioxide (CO₂) capture plants are expected to be situated near to power plants and other large industrial sources. Afterward CO₂ is to be transported to storage site using one or a combination of transport media: truck, train, ship or pipeline. Transport by pipeline is considered the preferred option for large quantities of CO₂ over long distances. The hazard connected with this kind of transportation can be considered an emerging risk and is the subject of this paper.The paper describes the Quantitative Risk Assessment of a hypothetical network pipeline located in UK, in particular the study of consequences due to a CO₂ release from pipeline.The risk analysis highlighted that some sections of pipeline network cross densely populated areas. For this reason, some changes in the original path of the network have been proposed in order to achieve a significant reduction in the societal risk.