“Basin scale” versus “localized” pore pressure/stress coupling - Implications for trap integrity evaluation

In petroleum industry, the difference between pore pressure (P_p) and minimum horizontal stress S_h (termed the seal or retention capacity) is of major consideration because it is often assumed to represent how close a system is to hydraulic failure and thus the maximum hydrocarbon column height that can be maintained. While S_h and P_p are often considered to be independent parameters, several studies in the last decade have demonstrated that S_h and P_p are in fact coupled. However, the nature of this coupling relationship remains poorly understood. In this paper, we explore the influences of the spatial pore pressure distribution on S_h/P_p coupling and then on failure pressure predictions and trap integrity evaluation. With analytical models, we predict the fluid pressure sustainable within a reservoir before failure of its overpressured shale cover. We verify our analytical predictions with experiments involving analogue materials and fluids. We show that hydraulic fracturing and seal breach occur for fluid pressure greater than it would be expected from conventional retention capacity. This can be explained by the impact of the fluid overpressure field in the overburden and the pressure diffusion around the reservoir on the principal stresses. We calculate that supralithostatic pressure could locally be reached in overpressured covers. We also define the retention capacity of a cover (RC) surrounding a fluid source or reservoir as the difference between the failure pressure and the fluid overpressure prevailing in shale at the same depth. In response to a localized fluid pressure rise, we show that the retention capacity does not only depend on the pore fluid overpressure of the overburden but also on the tensile strength of the cover, its Poisson’s ratio, and the depth and width of the fluid source.     

Comparison of the LEW88516 and ALHA77005 martian meteorites: Similar but distinct

Abstract By mineral and bulk compositions, the Lewis Cliff (LEW) 88516 meteorite is quite similar to the ALHA77005 martian meteorite. These two meteorites are not paired because their mineral compositions are distinct, they were found 500 km apart in ice fields with different sources for meteorites, and their terrestrial residence ages are different. Minerals in LEW88516 include: olivine, pyroxenes (low- and high-Ca), and maskelynite (after plagioclase); and the minor minerals chromite, whitlockite, ilmenite, and pyrrhotite. Mineral grains in LEW88516 range up to a few mm. Texturally, the meteorite is complex, with regions of olivine and chromite poikilitically enclosed in pyroxene, regions of interstitial basaltic texture, and glass-rich (shock) veinlets. Olivine compositions range from Fo₆₄ to Fo₇₀, (avg. Fo₆₇), more ferroan and with more variation than in ALHA77005 (Fo₆₉ to Fo₇₃). Pyroxene compositions fall between En₇₇Wo₄ and En₆₅Wo₁₅ and in clusters near En₆₃Wo₉ and En₅₃Wo₃₃, on average more magnesian and with more variation than in ALHA77005. Shock features in LEW88516 range from weak deformation through complete melting. Bulk chemical analyses by modal recombination of electron microprobe analyses, instrumental neutron activation, and radiochemical neutron activation confirm that LEW88516 is more closely related to ALHA77005 than to other known martian meteorites. Key element abundance ratios are typical of martian meteorites, as is its non-chondritic rare earth pattern. Differences between the chemical compositions of LEW88516 and ALHA77005 are consistent with slight differences in the proportions of their constituent minerals and not from fundamental petrogenetic differences. Noble gas abundances in LEW88516, like those in ALHA77005, show modest excesses of ⁴⁰Ar and ¹²⁹Xe from trapped (shock-implanted) gas. As with other ALHA77005 and the shergottite martian meteorites (except EETA79001), noble gas isotope abundances in LEW88516 are consistent with exposure to cosmic rays for 2.5–3 Ma. The absence of substantial effects of shielding from cosmic rays suggest LEW88516 spent this time as an object no larger than a few cm in diameter.     

Estimating coseismic deformations from near source strong motion records: methods and case studies

Digital strong-motion accelerographs have opened up the possibility of extracting ground motion characteristics at much lower frequencies than was offered by analogue instruments. High-quality digital data obtained close to the faults have tempted several efforts to retrieve permanent ground displacements after an earthquake. Such attempts have been partly successful, and somewhat subjective, the main reason being the presence of baseline offsets in the accelerometric data. We review existing methods for such applications, discuss their limitations and propose a more objective and improved scheme to make baseline adjustments and obtain permanent displacements. The proposed technique is applied to 26 digital recordings from the 29 May 2008 Ölfus Earthquake in South Iceland and 9 recordings from the 1999 Chi-Chi, Taiwan, Earthquake, and the permanent displacements obtained are compared with published results and GPS measurements from nearby stations. Our case studies show that the proposed technique, in addition to being simple and objective, is effective in making adjustments for baseline errors in accelerometric data.     

Computing the Sensitivity Kernels for 2.5-D Seismic Waveform Inversion in Heterogeneous,Anisotropic Media

2.5-D modeling and inversion techniques are much closer to reality than the simple and traditional 2-D seismic wave modeling and inversion. The sensitivity kernels required in full waveform seismic tomographic inversion are the Fréchet derivatives of the displacement vector with respect to the independent anisotropic model parameters of the subsurface. They give the sensitivity of the seismograms to changes in the model parameters. This paper applies two methods, called ‘the perturbation method’ and ‘the matrix method’, to derive the sensitivity kernels for 2.5-D seismic waveform inversion. We show that the two methods yield the same explicit expressions for the Fréchet derivatives using a constant-block model parameterization, and are available for both the line-source (2-D) and the point-source (2.5-D) cases. The method involves two Green’s function vectors and their gradients, as well as the derivatives of the elastic modulus tensor with respect to the independent model parameters. The two Green’s function vectors are the responses of the displacement vector to the two directed unit vectors located at the source and geophone positions, respectively; they can be generally obtained by numerical methods. The gradients of the Green’s function vectors may be approximated in the same manner as the differential computations in the forward modeling. The derivatives of the elastic modulus tensor with respect to the independent model parameters can be obtained analytically, dependent on the class of medium anisotropy. Explicit expressions are given for two special cases—isotropic and tilted transversely isotropic (TTI) media. Numerical examples are given for the latter case, which involves five independent elastic moduli (or Thomsen parameters) plus one angle defining the symmetry axis.     

Soil nutrients in relation to vertical roots distribution in the riparian zone of Three Gorges Reservoir,China

Since the impoundment of the Three Gorges Reservoir (TGR), the riparian zone has been subjected to numerous environmental changes. This study was conducted to recognize the distribution of grass roots and its impacts on soil nutrients in the water level fluctuation zone of TGR. Roots of four predominant herbaceous plants in the study area, specifically, Cynodon dactylon, Hemarthria altissima, Hemarthria compressa, and Paspalum paspaloides, and their corresponding relation with soil nutrient contents were investigated. Root surface area density was determined with WinRHIZO, and the relationships of root distribution with soil depths and soil nutrient contents were studied. The results indicates that most roots are distributed in the top soil layer of 0–10 cm. Estimated root surface area density for the selected grass species ranges from 0.16 to 13.44 cm2/cm3, and decreases exponentially with an increase in soil depth. Soil organic matter and total nitrogen contents are significantly lower on bare control area than the corresponding values on the grasslands. Total nutrient contents on grasslands of C. dactylon and H. compressa are higher than those of other grass areas. Root length density and root surface area density are significantly correlated with soil organic matter and total nitrogen content for the four grasslands. The present results suggests that plant roots have significant effects on the distribution of soil nutrients in soil profiles in the riparian zone along the TGR. Nevertheless, additional investigations are needed to reveal the specific interactions between plant roots distribution, soil nutrients and water level fluctuations.     

Estimation of groundwater evaporation and salt flux from Owens Lake, California, USA

Groundwater evaporation and subsequent precipitation of soluble salts at Owens Lake in eastern California have created one of the single largest sources of airborne dust in the USA, yet the evaporation and salt flux have not been fully quantified. In this study, we compare eddy correlation, microlysimeters and solute profiling methods to determine their validity and sensitivity in playa environments. These techniques are often used to estimate evaporative losses, yet have not been critically compared at one field site to judge their relative effectiveness and accuracy. Results suggest that eddy correlation methods are the most widely applicable for the variety of conditions found on large playa lakes. Chloride profiling is shown to be highly sensitive to thermal and density-driven fluxes in the near surface and, as a result, appears to underestimate yearly groundwater evaporation. Yearly mean groundwater evaporation from the playa surface estimated from the three study areas was found to range from 88 to 104 mm year⁻¹, whereas mean evaporation from the brine-covered areas was 872 mm year⁻¹. Uncertainties on these mean rates were estimated to be ±25%, based on comparisons between eddy correlation and lysimeter estimates. On a yearly basis, evaporation accounts for approximately 47 × 10⁶ m³ of water loss from the playa surface and open-water areas of the lake. Over the playa area, as much as 7.5 × 10⁸ kg (7.5 × 10⁵ t) of salt are annually concentrated by evaporation at or near the playa surface, much of which appears to be lost during dust storms in area.     

Petrogenesis and paleotectonic history of the Wild Bight Group,an Ordovician rifted island arc in central Newfoundland

The Wild Bight Group (WBG) is a sequence of early and middle Ordovician volcanic, subvolcanic and epiclastic rocks, part of the Dunnage Tectonostratigraphic Zone of the Newfoundland Appalachians. A detailed geochemical and Nd-isotopic study of the volcanic and subvolcanic rocks has been carried out to determine the geochemical characteristics of the rocks, interpret their palcotectonic environments and constrain their petrogenetic history. The lower and central stratigraphic levels of the WBG contain mafic volcanic rocks with island-arc geochemical signatures, including LREE-enriched are tholeiites with _Nd(t) =-0.1 to +2.2 (type A-I), LREE-depleted arc tholeiites with _Nd(t) =+5.6 to +7.1 (type A-II) and an unusual suite of strongly incompatible-element depleted tholeiites in which _Nd(t) ranges from-0.9 to +4.6 and is negatively correlated with¹⁴⁷Sm/¹⁴⁴Nd (type A-III). High-silica, low-K rhyolites occur locally in the central part of the stratigraphy, associated with mafic rocks of arc affinity, and have _Nd(t) =+4.7 to +5.4. The upper stratigraphic levels of the WBG dominantly contain rocks with non-arc geochemical signatures, including alkalic basalts with _Nd(t) =+4.6 to +5.5 (type N-I), strongly LREE- and incompatible element-enriched tholeiites that are transitional between alkalic and non-alkalic rocks with _Nd(t) =+4.4 to +7.0 (type N-II) and rocks with flat to slightly LREE-enriched patterns and _Nd(t) =+5.1 to +7.4 (type N-III). Rocks with non-arc and arc signatures are locally interbedded near the stratigraphic type of the WBG. Nd-isotopic data in the type A-I and A-II rocks are generally compatible with mixing/partial melting models involving depleted mantle, variably contaminated by a subducted crustally-derived sediment. The petrogenesis of type A-III rocks must involve source mixing and multi-stage partial melting, but the details are not clear. The geochemistry and Nd isotope data for types N-I, N-II and N-III rocks are compatible with petrogenetic models involving variable partial melting of a source similar to that postulated for modern oceanic island basalts. Comparison of the WBG with modern analogues suggests a 3-stage developmental model: stage 1) island-arc volcanism (eruption of type mafic volcancs); stage 2) arc-rifting (continued eruption of type A-I, A-I, eruption of types A-II and A-III mafic volcanics and high-silica, low-K rhyolites); and stage 3) back-arc basin volcanism (continued minor eruption of type A-I basalts, eruption of types N-I, N-II, N-III basalts). Stages 1 and 2 volcanism involved partial melting of subduction contaminated mantle, while stage 3 volcanism utilized depleted-mantle sources not affected by the subducting slab. This model provides a basis for interpreting coeval sequences in central Newfoundland and a comparative framework for some early Paleozoic oceanic volcanic sequences elsewhere in the Appalachian orogen.     

Petrology and geochemistry of greywackes from Goa-Dharwar sector,western Dharwar Craton: Implications for volcanoclastic origin

Late Archaean Supracrustals of the Goa-Dharwar sector (GDS) are composed of a thick sequence of greywacke sequence with narrow intercalations of quartzite, BIF and carbonates. Mafic volcanics occupy the base of the sequence. The greywackes are predominantly tuffacious containing chlorite-sericite and hornblende. Arkosic variations containing biotite dominate the western part of the sector. Fine-grained variations occur as isolated narrow lenses within other types of greywackes. The conglomeratic greywackes are localized along the western and the eastern margins of the sector. All of the greywackes are all typically immature containing coarser clasts of mostly plagioclase (18–23%) and quartz (32–34%). Lithic fragments of felsic volcanic rocks are common. The matrix is dominated by mafic material. Biotite and amphibole are related to metamorphic recrystallization. Chlorite, sericite, epidote, carbonate and chert are products of the interplay of diagenesis and low-grade metamorphism. Fe-Ti oxide, sphene, apatite and zircon are usual accessories. But for slight enrichment in K₂O, the major element chemistry of the GDS greywackes is similar to the chemistry of Late Archaean greywackes. They also compare in respect of V, Co, Hf contents, K₂O/Na₂O, SiO₂/Al₂O₃, Na₂O/Al₂O₃, Ba/Rb, Th/U, La/Th, Sm/Nd ratios, steep REE patterns with distinct LREE enrichment and HREE depletion. The GDS greywackes however are distinctly enriched in Rb, Ba, Sr, Th, U, Cu, Zr, Ce/Ce* and depleted in Cr, Ni, and Zn. The conglomeritic and biotite bearing verities contain considerable proportions of clasts derived from the basement tonalitic/granitic terrain. The common tuffacious greywackes containing hornblende and biotite-sericite however include only volcanic clasts and bear evidence of derivation from submarine weathering of predominantly felsic volcanics erupted on a large scale to form a magmatic arc in the later stages of geosynclinal deposition. Geochemical data suggest that the GDS greywackes were laid down in progressively changing basin geometry from a passive to active continental margin and island arc setting.     

A literature review of the saturation state of seawater with respect to calcium carbonate and its possible significance for scale formation on OTEC heat exchangers

An investigation has been made of available data on the saturation state of seawater with respect to calcium carbonate and its possible significance for scale formation on Ocean Thermal Energy Conversion (OTEC) heat exchangers. Pertinent oceanographic data is lacking at or near potential OTEC sites for the calculation of the degree of saturation of seawater with respect to calcium carbonate. Consequently, only “extrapolated” saturation values can be used. These indicate that near surface seawater is probably supersaturated, with respect to the calcium carbonate phases calcite and aragonite, at all potential OTEC sites. The deep seawater that would be brought to the surface at the potential Atlantic Ocean sites is also likely to be supersaturated with respect to calcium carbonate. The deep seawater at the potential Pacific Ocean sites may be slightly undersaturated.The fact that OTEC heat exchangers will be operating in seawater, which is supersaturated with respect to calcium carbonate, means that if nucleation of calcite or aragonite occurs on the heat exchanger surfaces, significant growth rates of calcium carbonate scale may be expected. The potential for calcium carbonate nucleation is highest at cathodic metal surface locations, which are produced as the result of aluminum corrosion in seawater. Consequently, corrosion and scale formation may be closely related. What the possible effects of biofouling may be on this process are not known.     

Classifying zones of suitability for manual drilling using textural and hydraulic parameters of shallow aquifers: a case study in northwestern Senegal

A method is proposed that uses analysis of borehole stratigraphic logs for the characterization of shallow aquifers and for the assessment of areas suitable for manual drilling. The model is based on available borehole-log parameters: depth to hard rock, depth to water, thickness of laterite and hydraulic transmissivity of the shallow aquifer. The model is applied to a study area in northwestern Senegal. A dataset of boreholes logs has been processed using a software package (TANGAFRIC) developed during the research. After a manual procedure to assign a standard category describing the lithological characteristics, the next step is the automated extraction of different textural parameters and the estimation of hydraulic conductivity using reference values available in the literature. The hydraulic conductivity values estimated from stratigraphic data have been partially validated, by comparing them with measured values from a series of pumping tests carried out in large-diameter wells. The results show that this method is able to produce a reliable interpretation of the shallow hydrogeological context using information generally available in the region. The research contributes to improving the identification of areas where conditions are suitable for manual drilling. This is achieved by applying the described method, based on a structured and semi-quantitative approach, to classify the zones of suitability for given manual drilling techniques using data available in most African countries. Ultimately, this work will support proposed international programs aimed at promoting low-cost water supply in Africa and enhancing access to safe drinking water for the population.