Formation pressure testing at the Mount Elbert Gas Hydrate Stratigraphic Test Well, Alaska North Slope: Operational summary, history matching, and interpretations

In February 2007, the U.S. Department of Energy, BP Exploration (Alaska), and the U.S. Geological Survey, collected open-hole pressure-response data, as well as gas and water sample collection, in a gas hydrate reservoir (the BPXA-DOE-USGS Mount Elbert Gas Hydrate Stratigraphic Test Well) using Schlumberger's Modular Dynamics Formation Tester (MDT) wireline tool. Four such MDT tests, ranging from six to twelve hours duration, and including a series of flow, sampling, and shut-in periods of various durations, were conducted. Locations for the testing were selected based on NMR and other log data to assure sufficient isolation from reservoir boundaries and zones of excess free water. Test stages in which pressure was reduced sufficiently to mobilize free water in the formation (yet not cause gas hydrate dissociation) produced readily interpretable pressure build-up profiles. Build-ups following larger drawdowns consistently showed gas-hydrate dissociation and gas release (as confirmed by optical fluid analyzer data), as well as progressive dampening of reservoir pressure build-up during sequential tests at a given MDT test station.History matches of one multi-stage, 12-h test (the C2 test) were accomplished using five different reservoir simulators: CMG-STARS, HydrateResSim, MH21-HYDRES, STOMP-HYD, and TOUGH + HYDRATE. Simulations utilized detailed information collected across the reservoir either obtained or determined from geophysical well logs, including thickness (11.3 m, 37 ft.), porosity (35%), hydrate saturation (65%), both mobile and immobile water saturations, intrinsic permeability (1000 mD), pore water salinity (5 ppt), and formation temperature (3.3-3.9 °C). This paper will present the approach and preliminary results of the history-matching efforts, including estimates of initial formation permeability and analyses of the various unique features exhibited by the MDT results.     

Use of formation pressure test results over a hydrate interval for long-term production forecasting at the Mount Elbert Gas Hydrate Stratigraphic Test Well, Alaska North Slope: Implications of uncertainties

In the first part of this paper, the pressure and temperature measurements obtained using Schlumberger's Modular Formation Dynamics Tester (MDT) conducted on the C2 interval in the BPXA-DOE-USGS Mount Elbert Gas Hydrate Stratigraphic Test Well (Mount Elbert Well) are history matched, with the following three objectives: (i) to obtain a better understanding of hydrate decomposition and its reformation as conditions cross the p/T stability, (ii) to obtain formation properties (e.g., permeability) that are consistent with the measurements, and (iii) to explore the non-uniqueness in the history match; i.e., to explore the ranges of parameters that allow a reasonable match of the measured quantities. In the second part of this paper, long-term production performance is predicted, and the effect of the uncertain parameters on the predictions is demonstrated. The results are used to demonstrate the range of long-term production that may be expected, when a model is calibrated using the MDT data. Usefulness of short-term tests for long-term forecast prediction is then discussed.     

Mount Elbert Gas Hydrate Stratigraphic Test Well, Alaska North Slope: Overview of scientific and technical program

The Mount Elbert Gas Hydrate Stratigraphic Test Well was drilled within the Alaska North Slope (ANS) Milne Point Unit (MPU) from February 3 to 19, 2007. The well was conducted as part of a Cooperative Research Agreement (CRA) project co-sponsored since 2001 by BP Exploration (Alaska), Inc. (BPXA) and the U.S. Department of Energy (DOE) in collaboration with the U.S. Geological Survey (USGS) to help determine whether ANS gas hydrate can become a technically and commercially viable gas resource. Early in the effort, regional reservoir characterization and reservoir simulation modeling studies indicated that up to 0.34 trillion cubic meters (tcm; 12 trillion cubic feet, tcf) gas may be technically recoverable from 0.92 tcm (33 tcf) gas-in-place within the Eileen gas hydrate accumulation near industry infrastructure within ANS MPU, Prudhoe Bay Unit (PBU), and Kuparuk River Unit (KRU) areas. To further constrain these estimates and to enable the selection of a test site for further data acquisition, the USGS reprocessed and interpreted MPU 3D seismic data provided by BPXA to delineate 14 prospects containing significant highly-saturated gas hydrate-bearing sand reservoirs. The “Mount Elbert” site was selected to drill a stratigraphic test well to acquire a full suite of wireline log, core, and formation pressure test data. Drilling results and data interpretation confirmed pre-drill predictions and thus increased confidence in both the prospect interpretation methods and in the wider ANS gas hydrate resource estimates. The interpreted data from the Mount Elbert well provide insight into and reduce uncertainty of key gas hydrate-bearing reservoir properties, enable further refinement and validation of the numerical simulation of the production potential of both MPU and broader ANS gas hydrate resources, and help determine viability of potential field sites for future extended term production testing. Drilling and data acquisition operations demonstrated that gas hydrate scientific research programs can be safely, effectively, and efficiently conducted within ANS infrastructure. The program success resulted in a technical team recommendation to project management to drill and complete a long-term production test within the area of existing ANS infrastructure. If approved by stakeholders, this long-term test would build on prior arctic research efforts to better constrain the potential gas rates and volumes that could be produced from gas hydrate-bearing sand reservoirs.     

The characteristics of gas hydrates recovered from the Mount Elbert Gas Hydrate Stratigraphic Test Well, Alaska North Slope

Systematic analyses have been carried out on two gas hydrate-bearing sediment core samples, HYPV4, which was preserved by CH₄ gas pressurization, and HYLN7, which was preserved in liquid-nitrogen, recovered from the BPXA-DOE-USGS Mount Elbert Stratigraphic Test Well. Gas hydrate in the studied core samples was found by observation to have developed in sediment pores, and the distribution of hydrate saturation in the cores imply that gas hydrate had experienced stepwise dissociation before it was stabilized by either liquid nitrogen or pressurizing gas. The gas hydrates were determined to be structure Type I hydrate with hydration numbers of approximately 6.1 by instrumentation methods such as powder X-ray diffraction, Raman spectroscopy and solid state ¹³C NMR. The hydrate gas composition was predominantly methane, and isotopic analysis showed that the methane was of thermogenic origin (mean δ¹³C = −48.6‰ and δ_D = −248‰ for sample HYLN7). Isotopic analysis of methane from sample HYPV4 revealed secondary hydrate formation from the pressurizing methane gas during storage.     

“Uniform geothermal gradient” and heat flow in the Qiongdongnan and Pearl River Mouth Basins of the South China Sea

The Pearl River Mouth Basin (PRMB) and Qiongdongnan Basin (QDNB) are oil and gas bearing basins in the northern margin of the South China Sea (SCS). Geothermal survey is an important tool in petroleum exploration. A large data set comprised of 199 thermal conductivities, 40 radioactive heat productions, 543 measured geothermal gradient values, and 224 heat flow values has been obtained from the two basins. However, the measured geothermal gradient data originated from diverse depth range make spatial comparison a challenging task. Taking into account the variation of conductivity and heat production of rocks, we use a “uniform geothermal gradient” to characterize the geothermal gradient distribution of the PRMB and QDNB. Results show that, in the depth interval of 0–5 km, the “uniform geothermal gradient” in the PRMB varies from 17.8 °C/km to 50.2 °C/km, with an average of 32.1 ± 6.0 °C/km. In comparison, the QDNB has an average “uniform geothermal gradient” of 31.9 ± 5.6 °C/km and a range between 19.7 °C/km and 39.5 °C/km. Heat flows in the PRMB and QDNB are 71.3 ± 13.5 mW/m² and 72.9 ± 14.2 mW/m², respectively. The heat flow and geothermal gradient of the PRMB and QDNB tend to increase from the continental shelf to continental slope owing to the lithosphereic/crustal thinning in the Cenozoic.     

Physical properties of sediment from the Mount Elbert Gas Hydrate Stratigraphic Test Well, Alaska North Slope

This study characterizes cored and logged sedimentary strata from the February 2007 BP Exploration Alaska, Department of Energy, U.S. Geological Survey (BPXA-DOE-USGS) Mount Elbert Gas Hydrate Stratigraphic Test Well on the Alaska North Slope (ANS). The physical-properties program analyzed core samples recovered from the well, and in conjunction with downhole geophysical logs, produced an extensive dataset including grain size, water content, porosity, grain density, bulk density, permeability, X-ray diffraction (XRD) mineralogy, nuclear magnetic resonance (NMR), and petrography.This study documents the physical property interrelationships in the well and demonstrates their correlation with the occurrence of gas hydrate. Gas hydrate (GH) occurs in three unconsolidated, coarse silt to fine sand intervals within the Paleocene and Eocene beds of the Sagavanirktok Formation: Unit D-GH (614.4 m-627.9 m); unit C-GH1 (649.8 m-660.8 m); and unit C-GH2 (663.2 m-666.3 m). These intervals are overlain by fine to coarse silt intervals with greater clay content. A deeper interval (unit B) is similar lithologically to the gas-hydrate-bearing strata; however, it is water-saturated and contains no hydrate.In this system it appears that high sediment permeability (k) is critical to the formation of concentrated hydrate deposits. Intervals D-GH and C-GH1 have average “plug” intrinsic permeability to nitrogen values of 1700 mD and 675 mD, respectively. These values are in strong contrast with those of the overlying, gas-hydrate-free sediments, which have k values of 5.7 mD and 49 mD, respectively, and thus would have provided effective seals to trap free gas. The relation between permeability and porosity critically influences the occurrence of GH. For example, an average increase of 4% in porosity increases permeability by an order of magnitude, but the presence of a second fluid (e.g., methane from dissociating gas hydrate) in the reservoir reduces permeability by more than an order of magnitude.     

Alaska North Slope regional gas hydrate production modeling forecasts

A series of gas hydrate development scenarios were created to assess the range of outcomes predicted for the possible development of the “Eileen” gas hydrate accumulation, North Slope, Alaska. Production forecasts for the “reference case” were built using the 2002 Mallik production tests, mechanistic simulation, and geologic studies conducted by the US Geological Survey. Three additional scenarios were considered: A “downside-scenario” which fails to identify viable production, an “upside-scenario” describes results that are better than expected. To capture the full range of possible outcomes and balance the downside case, an “extreme upside scenario” assumes each well is exceptionally productive.Starting with a representative type-well simulation forecasts, field development timing is applied and the sum of individual well forecasts creating the field-wide production forecast. This technique is commonly used to schedule large-scale resource plays where drilling schedules are complex and production forecasts must account for many changing parameters. The complementary forecasts of rig count, capital investment, and cash flow can be used in a pre-appraisal assessment of potential commercial viability.Since no significant gas sales are currently possible on the North Slope of Alaska, typical parameters were used to create downside, reference, and upside case forecasts that predict from 0 to 71 BM³ (2.5 tcf) of gas may be produced in 20 years and nearly 283 BM³ (10 tcf) ultimate recovery after 100 years.Outlining a range of possible outcomes enables decision makers to visualize the pace and milestones that will be required to evaluate gas hydrate resource development in the Eileen accumulation. Critical values of peak production rate, time to meaningful production volumes, and investments required to rule out a downside case are provided. Upside cases identify potential if both depressurization and thermal stimulation yield positive results. An “extreme upside” case captures the full potential of unconstrained development with widely spaced wells. The results of this study indicate that recoverable gas hydrate resources may exist in the Eileen accumulation and that it represents a good opportunity for continued research.     

Formation history and physical properties of sediments from the Mount Elbert Gas Hydrate Stratigraphic Test Well, Alaska North Slope

The synthesis of available geological information and surface temperature evolution in the Alaska North Slope region suggests that: biogenic and deeper thermogenic gases migrated through fault networks and preferentially invaded coarse-grained layers that have relatively high hydraulic conductivity and low gas entry pressures; hydrate started forming before the beginning of the permafrost; eventually, the permafrost deepened and any remaining free water froze so that ice and hydrate may coexist at some elevations. The single tested specimen (depth 620.47-620.62 m) from the D unit consists of uncemented quartzitic fine sand with a high fraction of fines (56% by mass finer than sieve #200). The as-received specimen shows no evidence of gas present. The surface texture of sediment grains is compatible with a fluvial-deltaic sedimentation environment and shows no signs of glacial entrainment. Tests conducted on sediments with and without THF hydrates show that effective stress, porosity, and hydrate saturation are the major controls on the mechanical and geophysical properties. Previously derived relationships between these variables and mechanical/geophysical parameters properly fit the measurements gathered with Mount Elbert specimens at different hydrate saturations and effective stress levels. We show that these measurements can be combined with index properties and empirical geomechanical relationships to estimate engineering design parameters. Volumetric strains measured during hydrate dissociation vanish at 2-4 MPa; therefore, minimal volumetric strains are anticipated during gas production at the Mount Elbert well. However, volume changes could increase if extensive grain crushing takes place during depressurization-driven production strategies, if the sediment has unexpectedly high in situ porosity associated to the formation history, or if fines migration and clogging cause a situation of sustained sand production.     

Pore fluid geochemistry from the Mount Elbert Gas Hydrate Stratigraphic Test Well, Alaska North Slope

The BPXA-DOE-USGS Mount Elbert Gas Hydrate Stratigraphic Test Well was drilled and cored from 606.5 to 760.1 m on the North Slope of Alaska, to evaluate the occurrence, distribution and formation of gas hydrate in sediments below the base of the ice-bearing permafrost. Both the dissolved chloride and the isotopic composition of the water co-vary in the gas hydrate-bearing zones, consistent with gas hydrate dissociation during core recovery, and they provide independent indicators to constrain the zone of gas hydrate occurrence. Analyses of chloride and water isotope data indicate that an observed increase in salinity towards the top of the cored section reflects the presence of residual fluids from ion exclusion during ice formation at the base of the permafrost layer. These salinity changes are the main factor controlling major and minor ion distributions in the Mount Elbert Well. The resulting background chloride can be simulated with a one-dimensional diffusion model, and the results suggest that the ion exclusion at the top of the cored section reflects deepening of the permafrost layer following the last glaciation (∼100 kyr), consistent with published thermal models. Gas hydrate saturation values estimated from dissolved chloride agree with estimates based on logging data when the gas hydrate occupies more than 20% of the pore space; the correlation is less robust at lower saturation values. The highest gas hydrate concentrations at the Mount Elbert Well are clearly associated with coarse-grained sedimentary sections, as expected from theoretical calculations and field observations in marine and other arctic sediment cores.     

Downhole well log and core montages from the Mount Elbert Gas Hydrate Stratigraphic Test Well, Alaska North Slope

The BPXA-DOE-USGS Mount Elbert Gas Hydrate Stratigraphic Test Well was an integral part of an ongoing project to determine the future energy resource potential of gas hydrates on the Alaska North Slope. As part of this effort, the Mount Elbert well included an advanced downhole geophysical logging program. Because gas hydrate is unstable at ground surface pressure and temperature conditions, a major emphasis was placed on the downhole-logging program to determine the occurrence of gas hydrates and the in-situ physical properties of the sediments. In support of this effort, well-log and core data montages have been compiled which include downhole log and core-data obtained from the gas-hydrate-bearing sedimentary section in the Mount Elbert well. Also shown are numerous reservoir parameters, including gas-hydrate saturation and sediment porosity log traces calculated from available downhole well log and core data.     

Pre- and post-drill comparison of the Mount Elbert gas hydrate prospect, Alaska North Slope

In 2006, the United States Geological Survey (USGS) completed a detailed analysis and interpretation of available 2-D and 3-D seismic data, along with seismic modeling and correlation with specially processed downhole well log data for identifying potential gas hydrate accumulations on the North Slope of Alaska. A methodology was developed for identifying sub-permafrost gas hydrate prospects within the gas hydrate stability zone in the Milne Point area. The study revealed a total of 14 gas hydrate prospects in this area.In order to validate the gas hydrate prospecting protocol of the USGS and to acquire critical reservoir data needed to develop a longer-term production testing program, a stratigraphic test well was drilled at the Mount Elbert prospect in the Milne Point area in early 2007. The drilling confirmed the presence of two prominent gas-hydrate-bearing units in the Mount Elbert prospect, and high quality well logs and core data were acquired. The post-drill results indicate pre-drill predictions of the reservoir thickness and the gas-hydrate saturations based on seismic and existing well data were 90% accurate for the upper unit (hydrate unit D) and 70% accurate for the lower unit (hydrate unit C), confirming the validity of the USGS approach to gas hydrate prospecting. The Mount Elbert prospect is the first gas hydrate accumulation on the North Slope of Alaska identified primarily on the basis of seismic attribute analysis and specially processed downhole log data. Post-drill well log data enabled a better constraint of the elastic model and the development of an improved approach to the gas hydrate prospecting using seismic attributes.     

Intercalibration of benthic flux chambers I. Accuracy of flux measurements and influence of chamber hydrodynamics

The hydrodynamic properties and the capability to measure sediment-water solute fluxes, at assumed steady state conditions, were compared for three radically different benthic chamber designs: the “Microcosm”, the “Mississippi” and the “Göteborg” chambers. The hydrodynamic properties were characterized by mounting a PVC bottom in each chamber and measuring mixing time, diffusive boundary layer thickness (DBL thickness) shear velocity (u∗), and total pressure created by the water mixing. The Microcosm had the most even distribution of DBL thickness and u∗, but the highest differential pressure at high water mixing rates. The Mississippi chamber had low differential pressures at high u∗. The Göteborg chamber was in between the two others regarding these properties. DBL thickness and u∗ were found to correlate according to the following empirical formula: DBL=76.18(u∗)−0.933. Multiple flux incubations with replicates of each of the chamber types were carried out on homogenized, macrofauna-free sediments in four tanks. The degree of homogeneity was determined by calculating solute fluxes (of oxygen, silicate, phosphate and ammonium) from porewater profiles and by sampling for porosity, organic carbon and meiofauna. All these results, except meiofauna, indicated that there were no significant horizontal variations within the sediment in any of the parallel incubation experiments. The statistical evaluations also suggested that the occasional variations in meiofauna abundance did not have any influence on the measured solute fluxes. Forty-three microelectrode profiles of oxygen in the DBL and porewater were evaluated with four different procedures to calculate diffusive fluxes. The procedure presented by Berg, Risgaard-Petersen and Rysgaard, 1989 [Limnol. Oceanogr. 43, 1500] was found to be superior because of its ability to fit measured profiles accurately, and because it takes into consideration vertical zonation with different oxygen consumption rates in the sediment. During the flux incubations, the mixing in the chambers was replicated ranging from slow mixing to just noticeable sediment resuspension. In the “hydrodynamic characterizations” these mixing rates corresponded to average DBL thickness from 120 to 550 μm, to u∗ from 0.12 to 0.68 cm/s, and to differential pressures from 0-3 Pa. Although not directly transferable, since the incubations were done on a “real” sediment with a rougher surface while in the characterizations a PVC plate simulated the sediments surface, these data give ideas about the prevailing hydrodynamic condition in the chambers during the incubations. The variations in water mixing did not generate statistically significant differences between the chamber types for any of the measured fluxes of oxygen or nutrients. Consequently it can be concluded that, for these non-permeable sediments and so long as appropriate water mixing (within the ranges given above) is maintained, the type of stirring mechanism and chamber design used were not critical for the magnitude of the measured fluxes. The average measured oxygen flux was 11.2 ± 2.7 (from 40 incubations), while the diffusive flux calculated (from 43 profiles using the Berg et al., 1989 [Limnol. Oceanogr. 43, 1500] procedure) was 11.1 ± 3.0 mmol m⁻² day⁻¹. This strongly suggests that accurate oxygen flux measurements were obtained with the three types of benthic chambers used and that the oxygen uptake is diffusive.     

Gas production from a cold, stratigraphically-bounded gas hydrate deposit at the Mount Elbert Gas Hydrate Stratigraphic Test Well, Alaska North Slope: Implications of uncertainties

As part of an effort to identify suitable targets for a planned long-term field test, we investigate by means of numerical simulation the gas production potential from unit D, a stratigraphically bounded (Class 3) permafrost-associated hydrate occurrence penetrated in the BPXA-DOE-USGS Mount Elbert Gas Hydrate Stratigraphic Test Well on North Slope, Alaska. This shallow, low-pressure deposit has high porosities (? = 0.4), high intrinsic permeabilities (k = 10⁻¹² m²) and high hydrate saturations (S_H = 0.65). It has a low temperature (T = 2.3-2.6 °C) because of its proximity to the overlying permafrost. The simulation results indicate that vertical wells operating at a constant bottomhole pressure would produce at very low rates for a very long period. Horizontal wells increase gas production by almost two orders of magnitude, but production remains low. Sensitivity analysis indicates that the initial deposit temperature is by the far the most important factor determining production performance (and the most effective criterion for target selection) because it controls the sensible heat available to fuel dissociation. Thus, a 1 °C increase in temperature is sufficient to increase the production rate by a factor of almost 8. Production also increases with a decreasing hydrate saturation (because of a larger effective permeability for a given k), and is favored (to a lesser extent) by anisotropy.     

Persistently declining oxygen levels in the interior waters of the eastern subarctic Pacific

Fifty years of measurements at Ocean Station Papa (OSP, 50°N, 145°W) show trends in the interior waters of the subarctic Pacific that are both impacted by short term (few years to bi-decadal) atmospheric or ocean circulation oscillations and by persistent climate trends. Between 1956 and 2006, waters below the ocean mixed layer to a depth of at least 1000 m have been warming and losing oxygen. On density surfaces found in the depth range 100-400 m (σ_θ = 26.3-27.0), the ocean is warming at 0.005-0.012 °C y⁻¹, whereas oxygen is declining at 0.39-0.70 μmol kg⁻¹ y⁻¹ or at an integrated rate of 123 mmol m⁻² y⁻¹ (decrease of 22% over 50 years). During this time, the hypoxic boundary (defined as 60 μmol O₂ kg⁻¹) has shoaled from ∼400 to 300 m. In the Alaska Gyre, the 26.2 isopycnal occasionally ventilates, whereas at OSP 26.0σ_θ has not been seen at the ocean surface since 1971 as the upper ocean continues to stratify. To interpret the 50 year record at OSP, the isopycnal transport of oxygenated waters within the interior of the subarctic Pacific is assessed by using a slightly modified “NO” parameter [Broecker, W., 1974. “NO” a conservative water-mass tracer. Earth and Planetary Science Letters 23, 100-107]. The highest nitrate-oxygen signature in interior waters of the North Pacific is found in the Bering Sea Gyre, Western Subarctic Gyre and East Kamchatka Current region as a consequence of winter mixing to the ∼26.6 isopycnal. By mixing with low NO waters found in the subtropics and Okhotsk Sea, this signature is diluted as waters flow eastward across the Pacific. Evidence of low NO waters flowing north from California is seen along the coasts of British Columbia and SE Alaska. Oxygen in the subsurface waters of the Alaskan Gyre was supplied ∼60% by subarctic and 40% by subtropical waters during WOCE surveys, whereas such estimates are shown to periodically vary by 20% at OSP. Other features discernable in the OSP data include periods of increased ventilation of deeper isopycnals on an ∼18 year cycle and strong, short term (few month) variability caused by passing mesoscale eddies. The potential impacts of declining oxygen on coastal ecosystems are discussed.     

Mud volcanoes and gas hydrates in the Anaximander mountains (Eastern Mediterranean Sea)

Detailed multibeam, sedimentological, and geophysical surveys provide ample new data to confirm that the Anaximander Mountains (Eastern Mediterranean) are an important area for active mud volcanism and gas hydrate formation. More than 3000 km of multibeam track length was acquired during two recent missions and 80 gravity and box cores were recovered. Morphology and backscatter data of the study area have better resolution than previous surveys, and very detailed morphology maps have been made of the known targeted mud volcanoes (Amsterdam, Kazan and Kula), especially the Amsterdam “crater” and the related mud breccia flows. Gas hydrates collected repeatedly from a large area of Amsterdam mud volcano at a sub-bottom depth of around 0.3–1.5 m resemble compacted snow and have a rather flaky form. New gas hydrate sites were found at Amsterdam mud volcano, including the mud flow sloping off to the south. Gas hydrates sampled for the first time at Kazan mud volcano are dispersed throughout the core samples deeper than 0.3 m and display a ‘rice’-like appearance. Relative chronology and AMS dating of interbedded pelagic sediments (Late Holocene hemipelagic, sapropel layer S1 and ash layers) within the mud flows indicate that successive eruptions of Kula mud volcano have a periodicity of about 5–10 kyrs. New mud volcanoes identified on the basis of multibeam backscatter intensity were sampled, documented as active and named “Athina” and “Thessaloniki”. Gas hydrates were sampled also in Thessaloniki mud volcano, the shallowest (1264 m) among all the active Mediterranean sites, at the boundary of the gas hydrate stability zone. Biostratigraphical analyses of mud breccia clasts indicated that the source of the subsurface sedimentary sequences consists of Late Cretaceous limestones, Paleocene siliciclastic rocks, Eocene biogenic limestones and Miocene mudstones. Rough estimations of the total capacity of the Anaximander mud volcanoes in methane gas are 2.56–6.40 km³.     

Investigation of the hydrate reservoir production under different depressurization modes

Abstract

Large reserves of natural gas hydrates exist, and the depressurization method has the greatest potential for gas hydrate reservoir recovery. Currently, the most commonly adopted depressurization simulation method is a constant bottom-hole pressure production scheme. This study proposes a new depressurization mode with decreasing bottom-hole pressure. The production characteristic was numerically investigated using this method. The results show the following: (1) As the depressurization exponent (n) decreases, the development effect improves, and production indexes including cumulative gas production/dissociation and gas-water ratio increase. However, the reservoir energy consumption is higher and the hydrate reformation is more severe. (2) Compared to the proposed depressurization mode, the hydrate production index of the constant bottom-hole pressure production (n?=?0) is better. However, the hydrate reservoir energy consumption is higher and the hydrate reformation is more severe using constant bottom-hole pressure production. (3) To achieve a balance between production and reservoir energy consumption during depressurization production, the bottom-hole pressure should be controlled by selecting a suitable depressurization exponent between n_min and n_max, which can be determined through numerical simulations.     

Particulate organic carbon fluxes to the ocean interior and factors controlling the biological pump: A synthesis of global sediment trap programs since 1983

Particulate organic carbon (POC) is vertically transported to the oceanic interior by aggregates and their ballasts, mainly CaCO₃ and biogenic opal, with a smaller role for lithogenic aerosols through the mesopelagic zone. Diel migrating zooplankton communities effect vertical transport and remineralization of POC in the upper layers of the ocean. Below 1.5 km, the presence of zooplankton is reduced and thus the aggregates travel mainly by gravitational transport. We normalized the fluxes of POC, CaCO₃, and biogenic opal from data published on samples collected at 134 globally distributed, bottom-tethered, time-series sediment trap (TS-trap) stations to annual mole fluxes at the mesopelagic/bathypelagic boundary (m/b) at 2 km and defined them as F_m/bC_org, F_m/bC_inorg, and F_m/bSi_bio. Using this global data set, we investigated (1) the geographic contrasts of POC export at m/b and (2) the supply rate of ∑CO₂ to the world mesopelagic water column. F_m/bC_org varies from 25 (Pacific Warm Pool) to 605 (divergent Arabian Sea) mmolC m⁻² yr⁻¹; F_m/bC_inorg varies from >8 (high latitude Polar Oceans) or 15 (Pacific Warm Pool) to 459 (divergent Arabian Sea) mmolC m⁻² yr⁻¹; and F_m/bSi_bio, the most spatially/temporally variable flux, ranges from 6 (North Atlantic Drift) to 1118 (Pacific Subarctic Gyre) mmolSi m⁻² yr⁻¹. The oceanic region exhibiting the highest POC flux over a significantly large region is the area of the North Pacific Boreal Gyres where the average F_m/bC_org = 213, F_m/bC_inorg = 126, and F_m/bSi_bio = 578 mmol m⁻² yr⁻¹. F_m/bC_org and F_m/bC_inorg are particularly high in large upwelling margins, including the divergent Arabian Sea and off Cape Verde. One of the data sets showing the lowest flux over a significant region/basin is F_m/bC_org = 39, F_m/bC_inorg = 69, and F_m/bSi_bio = 22 mmol m⁻² yr⁻¹ in the North Pacific subtropical/tropical gyres; Pan-Atlantic average fluxes are similar except F_m/bSi_bio fluxes are even lower. Where C_org/C_inorg and Si_bio/C_inorg are <1 defines the “Carbonate Ocean”, and where these ratios are ?1 defines the “Silica Ocean”. The Carbonate Ocean occupies about 80% of the present world pelagic ocean between the two major oceanographic fronts, the North Pacific Polar Front and the Antarctic Polar Front, and the Silica Ocean is found on the polar sides of these fronts. The total global annual fluxes of F_m/bC_org, F_m/bC_inorg, and F_m/bSi_bio at m/b calculated by parameterizations of the export flux data from 134 stations are surprisingly similar; 36.2, 33.8, and 34.6 teramol yr⁻¹ (120, 112, and 114 mmol m⁻² yr⁻¹), respectively, resulting in a near uniform binary ratio between the above three elements of about one. The global ternary % ratios estimated from 152 TS-trap samples of the three elements are 35:32:33. From our global F_m/bC_org and a published model estimate of the global export production, we estimate the regeneration rate of CO₂ through the mesopelagic zone by the biological pump is 441 teramolC yr⁻¹. Based on our global F_m/bC_inorg and recently estimated global primary production of PIC, 36-86 teramolC yr⁻¹ of PIC is assumed to be dissolved within the upper 2 km of the water column.     

Stormwater plume detection by MODIS imagery in the southern California coastal ocean

Stormwater plumes in the southern California coastal ocean were detected by MODIS-Aqua satellite imagery and compared to ship-based data on surface salinity and fecal indicator bacterial (FIB) counts collected during the Bight'03 Regional Water Quality Program surveys in February–March of 2004 and 2005. MODIS imagery was processed using a combined near-infrared/shortwave-infrared (NIR-SWIR) atmospheric correction method, which substantially improved normalized water-leaving radiation (nLw) optical spectra in coastal waters with high turbidity. Plumes were detected using a minimum-distance supervised classification method based on nLw spectra averaged within the training areas, defined as circular zones of 1.5–5.0-km radii around field stations with a surface salinity of S < 32.0 (“plume”) and S > 33.0 (“ocean”). The plume optical signatures (i.e., the nLw differences between “plume” and “ocean”) were most evident during the first 2 days after the rainstorms. To assess the accuracy of plume detection, stations were classified into “plume” and “ocean” using two criteria: (1) “plume” included the stations with salinity below a certain threshold estimated from the maximum accuracy of plume detection; and (2) FIB counts in “plume” exceeded the California State Water Board standards. The salinity threshold between “plume” and “ocean” was estimated as 32.2. The total accuracy of plume detection in terms of surface salinity was not high (68% on average), seemingly because of imperfect correlation between plume salinity and ocean color. The accuracy of plume detection in terms of FIB exceedances was even lower (64% on average), resulting from low correlation between ocean color and bacterial contamination. Nevertheless, satellite imagery was shown to be a useful tool for the estimation of the extent of potentially polluted plumes, which was hardly achievable by direct sampling methods (in particular, because the grids of ship-based stations covered only small parts of the plumes detected via synoptic MODIS imagery). In most southern California coastal areas, the zones of bacterial contamination were much smaller than the areas of turbid plumes; an exception was the plume of the Tijuana River, where the zone of bacterial contamination was comparable with the zone of plume detected by ocean color.     

Phase fractionation and oil-condensate mass balance in the South Marsh Island Block 208–239 area,offshore Louisiana

Phase fractionation can strongly deplete oil of its volatile compounds in a regular and characteristic fashion. This process has affected oils to a remarkably uniform extent throughout the 30 × 15 km South Marsh Island 208–239 and Vermilion 30–31 area (including the Tiger Shoal, Starfak, Mound Point, Lighthouse Point, Amber, Trinity Shoal, and Aquamarine fields) just offshore Louisiana. Fractionation of the original “parent” oil likely occurred in the deep, relatively flat-lying Rob L sand that underlies the area, and produced gas-washed oils (mean API 33°) and gas condensates (mean API 50°) in a volume ratio of 1:3.5. Both fractionated oil and vapor migrated from the fractionation site to shallower reservoirs. However, the estimated ultimate production ratio of gas-washed oil to gas condensate in this group of fields is 1:0.32, about 11 times higher than would be expected on mass balance considerations alone. Thus, there is an apparent deficiency of producible gas condensate relative to the amount of producible oil for the entire study area and for every field in that area. In the case of the Tiger Shoal field, the ratio of industry-estimated ultimately producible oil to gas condensate is 1:1.1. Based on the production data, we conclude that either there is an additional 6.4 × 10⁶ m³ (43 MMbbl) of undiscovered and/or unproduced condensate in the area or that condensate has escaped preferentially in vapor form to the seafloor. The well-studied and nearly depleted Tiger Shoal field provides a good example of how chemical data can be analyzed in a way that contributes insight into the phase fractionation process and the remaining exploration potential of an area.