Computational fluid dynamics (CFD) based approach to consequence assessment of accidental release of hydrocarbon during storage and transportation

This thesis investigated the risk of accidental release of hydrocarbons during transportation and storage. Transportation of hydrocarbons from an offshore platform to processing units through subsea pipelines involves risk of release due to pipeline leakage resulting from corrosion, plastic deformation caused by seabed shakedown or damaged by contact with drifting iceberg. The environmental impacts of hydrocarbon dispersion can be severe. Overall safety and economic concerns of pipeline leakage at subsea environment are immense. A large leak can be detected by employing conventional technology such as, radar, intelligent pigging or chemical tracer but in a remote location like subsea or arctic, a small chronic leak may be undetected for a period of time. In case of storage, an accidental release of hydrocarbon from the storage tank could lead pool fire; further it could escalate to domino effects. This chain of accidents may lead to extremely severe consequences. Analyzing past accident scenarios it is observed that more than half of the industrial domino accidents involved fire as a primary event, and some other factors for instance, wind speed and direction, fuel type and engulfment of the compound. In this thesis, a computational fluid dynamics (CFD) approach is taken to model the subsea pipeline leak and the pool fire from a storage tank. A commercial software package ANSYS FLUENT Workbench 15 is used to model the subsea pipeline leakage. The CFD simulation results of four different types of fluids showed that the static pressure and pressure gradient along the axial length of the pipeline have a sharp signature variation near the leak orifice at steady state condition. Transient simulation is performed to obtain the acoustic signature of the pipe near leak orifice. The power spectral density (PSD) of acoustic signal is strong near the leak orifice and it dissipates as the distance and orientation from the leak orifice increase. The high-pressure fluid flow generates more noise than the low-pressure fluid flow. In order to model the pool fire from the storage tank, ANSYS CFX Workbench 14 is used. The CFD results show that the wind speed has significant contribution on the behavior of pool fire and its domino effects. The radiation contours are also obtained from CFD post processing, which can be applied for risk analysis. The outcome of this study will be helpful for better understanding of the domino effects of pool fire in complex geometrical settings of process industries. The attempt to reduce and prevent risks is discussed based on the results obtained from the numerical simulations of the numerical models.

Efficient Computation of Electromagnetic Waves in Hydrocarbon Exploration Using the Improved Numerical Mode Matching (NMM) Method

In this study, we developed and improved the numerical mode matching (NMM) method which has previously been shown to be a fast and robust semi-analytical solver to investigate the propagation of electromagnetic (EM) waves in an isotropic layered medium. The applicable models, such as cylindrical waveguide, optical fiber, and borehole with earth geological formation, are generally modeled as an axisymmetric structure which is an orthogonal-plano-cylindrically layered (OPCL) medium consisting of materials stratified planarly and layered concentrically in the orthogonal directions.

In this report, several important improvements have been made to extend applications of this efficient solver to the anisotropic OCPL medium. The formulas for anisotropic media with three different diagonal elements in the cylindrical coordinate system are deduced to expand its application to more general materials. The perfectly matched layer (PML) is incorporated along the radial direction as an absorbing boundary condition (ABC) to make the NMM method more accurate and efficient for wave diffusion problems in unbounded media and applicable to scattering problems with lossless media. We manipulate the weak form of Maxwell's equations and impose the correct boundary conditions at the cylindrical axis to solve the singularity problem which is ignored by all previous researchers. The spectral element method (SEM) is introduced to more efficiently compute the eigenmodes of higher accuracy with less unknowns, achieving a faster mode matching procedure between different horizontal layers. We also prove the relationship of the field between opposite mode indices for different types of excitations, which can reduce the computational time by half. The formulas for computing EM fields excited by an electric or magnetic dipole located at any position with an arbitrary orientation are deduced. And the excitation are generalized to line and surface current sources which can extend the application of NMM to the simulations of controlled source electromagnetic techniques. Numerical simulations have demonstrated the efficiency and accuracy of this method.

Finally, the improved numerical mode matching (NMM) method is introduced to efficiently compute the electromagnetic response of the induction tool from orthogonal transverse hydraulic fractures in open or cased boreholes in hydrocarbon exploration. The hydraulic fracture is modeled as a slim circular disk which is symmetric with respect to the borehole axis and filled with electrically conductive or magnetic proppant. The NMM solver is first validated by comparing the normalized secondary field with experimental measurements and a commercial software. Then we analyze quantitatively the induction response sensitivity of the fracture with different parameters, such as length, conductivity and permeability of the filled proppant, to evaluate the effectiveness of the induction logging tool for fracture detection and mapping. Casings with different thicknesses, conductivities and permeabilities are modeled together with the fractures in boreholes to investigate their effects for fracture detection. It reveals that the normalized secondary field will not be weakened at low frequencies, ensuring the induction tool is still applicable for fracture detection, though the attenuation of electromagnetic field through the casing is significant. A hybrid approach combining the NMM method and BCGS-FFT solver based integral equation has been proposed to efficiently simulate the open or cased borehole with tilted fractures which is a non-axisymmetric model.

(Table 1) Summary of low-molecular-weight hydrocarbon concentrations, organic carbon contents, Rock-Eval pyrolysis and lithology at DSDP Leg 79 Holes

A series of core samples taken during Cruise 79 of Glomar Challenger, drilling offshore Morocco (Mazagan Plateau), is analyzed for their low-molecular-weight hydrocarbon contents. Fifty-four samples from DSDP Holes 544A, 545, 547A, and 547B, deep frozen on board immediately after recovery, are studied by a hydrogen-stripping/thermovaporization technique combined with capillary gas chromatography. Thirty-eight compounds in the C2-C8 molecular range, including saturated, olefinic, and aromatic hydrocarbons, are identified. Because of large differences in organic carbon contents, the total C2-C8 hydrocarbon concentrations vary from about 20 to 1500 ng/g dry sediment weight in the whole sample series. Organic-carbon normalized values are about 3.2 x 10**4 ng/g Corg for Lithologic Subunits IIIA and IIIB at Site 545 (Cenomanian to Aptian) and 1.0 x 10**5 ng/g Corg for Unit V at Site 547 (Cenomanian to Albian) reflecting the slightly more advanced maturity stage at the latter site. Values exceeding 10**5 ng/g Corg (Site 545) and 2 x 10**5 ng/g Corg (Site 547) are associated with samples that are very lean in organic carbon and are generally rich in carbonate. These samples are enriched by small amounts of gaseous hydrocarbons. A detailed study of individual hydrocarbon concentrations, plotted against depth, reveal additional indications for migration phenomena. At Site 547, for instance, the most mobile hydrocarbons studied (e.g., ethane) appear to migrate by diffusion or a related process from more than 700 m depth toward the surface.

(Table 1) Low molecular weight hydrocarbon concentrations and organic carbon contents at DSDP Hole 71-511

C2-C8 hydrocarbons (36 compounds identified) from 56 shipboard sealed, deep-frozen core samples of DSDP Leg 71, Site 511, Falkland Plateau, South Atlantic, were analyzed by a combined hydrogen stripping-thermovaporization method. Concentrations, which represent hydrocarbons dissolved in the pore water and adsorbed to the mineral surfaces of the sediment, vary from 24 ng/g of dry weight sediment in Lithologic Unit 4 to 17,400 ng/g in Lithologic Unit 6 ("black shale" unit). Likewise, the organic carbon normalized C2-C8 hydrocarbon concentrations range from 104 to 3.5 x 105 ng/g Corg. The latter value is more than one order of magnitude lower than expected for petroleum source beds in the main phase of oil generation. The low maturity at 600 meters depth is further supported by light hydrocarbon concentration ratios. The change of the kerogen type from Lithologic Unit 5 (Type III) to 6 (Type II) is evidenced by changes in the C6 and C7 hydrocarbon composition. Redistribution phenomena are observed close to the Tertiary-Cretaceous unconformity and at the contact between the "black shale" unit and the overlying Cretaceous chalks and claystones. Otherwise, the low molecular weight hydrocarbons in Hole 511 are formed in situ and remain at their place of formation. The core samples turned out to be contaminated by large quantities of acetone, which is routinely used as a solvent during sampling procedures onboard Glomar Challenger.

(Table 1) Hydrocarbon gas composition of sediment and gas hydrate collected in video-guided grabs from Chapopote

In the Campeche Knolls, in the southern Gulf of Mexico, lava-like flows of solidified asphalt cover more than 1 square kilometer of the rim of a dissected salt dome at a depth of 3000 meters below sea level. Chemosynthetic tubeworms and bivalves colonize the sea floor near the asphalt, which chilled and contracted after discharge. The site also includes oil seeps, gas hydrate deposits, locally anoxic sediments, and slabs of authigenic carbonate. Asphalt volcanism creates a habitat for chemosynthetic life that may be widespread at great depth in the Gulf of Mexico.

Hydrocarbon analyses of submarine mud volcanoes (MV) in the Kumano forearc basin

Twelve submarine mud volcanoes (MV) in the Kumano forearc basin within the Nankai Trough subduction zone were investigated for hydrocarbon origins and fluid dynamics. Gas hydrates diagnostic for methane concentrations exceeding solubilities were recovered from MVs 2, 4, 5, and 10. Molecular ratios (C1/C2<250) and stable carbon isotopic compositions (d13C-CH4 >-40 per mil V-PDB) indicate that hydrate-bound hydrocarbons (HCs) at MVs 2, 4, and 10 are derived from thermal cracking of organic matter. Considering thermal gradients at the nearby IODP Sites C0009 and C0002, the likely formation depth of such HCs ranges between 2300 and 4300 m below seafloor (mbsf). With respect to basin sediment thickness and the minimum distance to the top of the plate boundary thrust we propose that the majority of HCs fueling the MVs is derived from sediments of the Cretaceous to Tertiary Shimanto belt below Pliocene/Pleistocene to recent basin sediments. Considering their sizes and appearances hydrates are suggested to be relicts of higher MV activity in the past, although the sporadic presence of vesicomyid clams at MV 2 showed that fluid migration is sufficient to nourish chemosynthesis-based organisms in places. Distributions of dissolved methane at MVs 3, 4, 5, and 8 pointed at fluid supply through one or few MV conduits and effective methane oxidation in the immediate subsurface. The aged nature of the hydrates suggests that the major portion of methane immediately below the top of the methane-containing sediment interval is fueled by current hydrate dissolution rather than active migration from greater depth.

Analyses of hydrocarbon gases and sulphate in sediments from the King George Basin, Antarctica

Thermogenic hydrocarbons, formed by the thermal alteration of organic matter, are encountered in several piston core stations in the King George Basin, Anatarctica. These hemipelagic sediments are being deposited in an area of active hydrothermalism, associated with the back-arc spreading in the Bransfield Strait. The lateral extent of sediments infiltrated by the hydrothermally influenced interstitial fluids is characterized by basalt diapiric intrusions and is delineated by an acoustically turbid zone in the sediments of the eastern part of the basin. Iron-sulphide-bearing veins and fractures cut across the sediment in several cores; they appear to be conduits for flow of hydrothermally altered fluids. These zones have the highest C2+ and ethene contents. The thermogenic hydrocarbons have molecular C1/(C2 + C3) ratios typically < 50 and delta13CH4 values between -38? and -48?, indicating an organic source which has undergone strong thermal stress. Several sediment cores also have mixed gas signatures, which indicate the presence of substantial amounts of bacterial gas, predominantly methane. Hydrocarbon generation in the King George Basin is thought to be a local phenomenon, resulting from submarine volcanism with temperatures in the range 70-150°C. There are no apparent seepages of hydrocarbons into the water column, and it is not believed that significant accumulation of thermogenic hydrocarbons reside in the basin.

Seasonal variation of n-alkanes and polycyclic aromatic hydrocarbon concentrations in PMsub(10) samples collected at urban sites of São Paulo state, Brazil

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Evaluation and comparison of durability of polycyclic aromatic hydrocarbon compounds (PAHs) and its relation to the lipid content in smoked fishes, Silver carp (Hypophthalmichthys molitrix), Rutilus frisii kutum and Herring (Alosa caspia)

Cold smoking method is one of the commonest ways for fish smoking. It is done by the smoke that is the result of burning hard and soft woods is smoking rooms. Smoke includes a number of chemical constructs and its main part is poly aromatic hydrocarbons. More than one hundred kinds of these constructs are recognized in smoke that is produced from saturated hydrocarbons resulted from the solution of the woods Ligno cellulose in high temperature and lack of oxygen conditions. The high poisoning potentials and carcinogenic features sixteen constructs among them are proved and observed on humans. In this research, the PAH compounds were identified and observed in a three month period after smoking and during storing among three types of smoked fishes Silver carp and Caspian sea Sefid and herring. They are the most produced and consumed smoked fish in Iran. To find the relationship between the concentrations of PAH constructs and the amount of lipid in fish, first, the amount of lipid were determined separately in the skin and flesh of 30 samples of each type. The method used was Bligh and Dyer (1959). PAH compounds derivation were made for all skin and flesh samples smoked fish using organic solvents with Soxeleh and the derived samples were injected to gas chromatography (GC) by Hamilton injectors for determining their components quality and their quantity. The height of the used column was 25 meters and its diameter was 0.32 mm with the silica filler, nitrogen gas as carrier and flame ionization detector (FID) that are special for these constructs. For data analysis, Statistical tests were used by computer soft ware identified that the difference in the amount of lipid within the flesh and skin of each species and also among each other is significant. The largest amount was in Herrings flesh and skin, 18.74% in skin and 14.47% in flesh. The least amount in the skin 4.19% and the flesh 3.10% of Sefid. The amount in Silver carp was 13.28%in skin and 8.16% in flesh. The examination of the PAH compounds in smoked fish showed that is carcinogenic compounds; exist in these in these fish with different quantities in each. It seems that its amount is directly related to the amount of their lipid. The amount is different in flesh and skin. One of the most important reasons is the direct content of smoke and the concentration of lipid in tissues of all three types. The maintenance of the smoked fish for three months showed that most of PAH compounds were solved and their density decreased. The changes in density within time in different in each type and in flesh and skin. The amount of their receiving in human through the consumption of the smoked fish depends on the resulted density, the way and the amount of consumption and now we can determine and execute standards for the maximum dosage per day and per month regarding effective factors.

A truncation in the Aryl Hydrocarbon Receptor of the CRL:WI(Han) rat does not affect the developmental toxicity of TCDD

The Aryl Hydrocarbon Receptor (AhR) is required for the toxicity of TCDD, and so the AhR of CRL:WI and CRL:WI(Han) rats was characterised. Western blot showed AhR proteins of ~110 and ~97 kDa in individual rats from both strains. The AhR cDNA from a CRL:WI(Han) rat with the ~110kDa protein revealed a sequence that was identical to that of the CRL:WI and SD rat. However, cloning of the AhR from a rat with the ~97kDa protein revealed a point mutation, and five variants encoding two C-terminally truncated variants of the AhR protein, arising from a point mutation in the intron/exon junction and consequent differential splicing. These C-terminally truncated variants were expressed and shown to give rise to a protein of ~97kDa; the recombinant AhR bound TCDD with an affinity that was not statistically different from the full-length protein. A single-nucleotide polymorphism (SNP) assay was developed, and showed that both alleles were represented in a Hardy-Weinberg equilibrium in samples of CRL:WI and CRL:WI(Han) populations; both alleles are abundant. Rats from two studies of TCDD developmental toxicity were genotyped, and the association with toxicity investigated using statistical analysis. There was no plausible evidence that the AhR allele had a significant effect on the toxic endpoints examined. These data show that the two AhR alleles are common in two strains of Wistar rat, and that the AhR alleles had no effect on TCDD-induced developmental toxicity in two independent studies.

Recombinant expression of Aryl Hydrocarbon Receptor for quantitative ligand-binding analysis

Recombinant expression of the Aryl Hydrocarbon Receptor (AhR) yields small amounts of ligand- binding competent AhR. Therefore, Spodoptera frugiperda (Sf9) cells and baculovirus have been evaluated for high level and functional expression of AhR. Rat and human AhR were expressed as soluble protein in significant amounts. Expression of ligand-binding competent AhR was sensitive to the protein concentration of Sf9 extract, and co-expression of the chaperone p23 failed to affect the yield of functional ligand-binding AhR. The expression system yielded high levels of functional protein, with the ligand-binding capacity (Bmax) typically 20- fold higher than that obtained with rat liver cytosol. Quantitative estimates of the ligand-binding affinity of human and rat AhR were obtained; the Kd for recombinant rat AhR was indistinguishable from that of native rat AhR, thereby validating the expression system as a faithful model for native AhR. The human AhR bound TCDD with significantly lower affinity than the rat AhR. These findings demonstrate high-level expression of ligand-binding competent AhR, and sufficient AhR for quantitative analysis of ligand-binding.

2-(4-Amino-3-methylphenyl)-5-fluorobenzothiazole is a ligand and shows species-specific partial agonism of the Aryl Hydrocarbon Receptor

2-(4-Amino-3-methylphenyl)-5-fluorobenzothiazole (5F 203) and related compounds are a series of anti-cancer candidate pharmaceuticals (Table 1.), that have been shown to activate the AhR. We show that these compounds are high affinity ligands for the rat AhR, but a quantitative assay for their ability to induce CYP1A1 RNA in H4IIEC3 cells, a measure of activation of the AhR, showed a poor relationship between affinity for the AhR and ability to induce CYP1A1 RNA. 5F 203, an agonist with low potency, was able to antagonise the induction of CYP1A1 RNA by TCDD, while IH 445, a potent agonist, did not antagonise the induction of CYP1A1 RNA by TCDD, and Schild analysis confirmed 5F 203 to be a potent antagonist of the induction of CYP1A1 RNA by TCDD in H4IIEC3 cells. In contrast, several benzothiazoles show potent induction of CYP1A1 RNA in human MCF-7 cells, and 5F 203 is unable to detectably antagonise the induction of CYP1A1 RNA in MCF-7 cells, showing a species difference in antagonism. Evaluation of the antiproliferative activity of benzothiazoles showed that the ability to agonise the AhR correlated with growth inhibition both in H4IIEC3 cells for a variety of benzothiazoles, and between H4IIEC3 and MCF-7 cells for 5F 203, suggesting an important role of agonism of the AhR in the anti-proliferative activity of benzothiazoles.