Evaluation of CBM reservoir fracture characteristics by azimuthal AVA inversion

Exploitation of coalbed methane (CBM) is directly affected by the strike and density of fractures, particularly because fractures can act as conduits for water and gas (after desorption). To realize the azimuthal amplitude versus angle (AVA) inversion for a CBM reservoir, we treated reservoirs whose fractures are aligned vertically as being equivalent to a horizontal transverse isotropic medium and simplified the P-P wave reflection coefficient formula to an approximately linear equation. Then, we adopted a modified damping least squares algorithm to solve the azimuthal AVA inversion problem. With edited macro-bins, corrected residual moveout, and normalized amplitude, the azimuthal gathers exhibited unambiguous azimuthal AVA characteristics. The #3 CBM reservoir in the Qinshui Basin was selected for a case study. The density and strike of the reservoir fractures were obtained and analyzed, and the reservoir was divided into three classes according to fracture density. Finally, in comparison with the gas production data, we confirmed that the inversion results and reservoir classification were both reasonable and credible.     

∼3,850 Ma tonalites in the Nuuk region,Greenland: geochemistry and their reworking within an Eoarchaean gneiss complex

The Eoarchaean (>3,600 Ma) Itsaq Gneiss Complex of southern West Greenland is dominated by polyphase orthogneisses with a complex Archaean tectonothermal history. Some of the orthogneisses have c. 3,850 Ma zircons, and they vary from rare single phase metatonalites to more common complexly banded migmatites. This is due to heterogeneous strain, in situ anatexis and granitic veining superimposed during younger tectonothermal events. In the single-phase tonalites with c. 3,850 Ma zircon, oscillatory-zoned prismatic zircon is all 3,850 Ma old, but shows patchy ancient loss of radiogenic Pb. SHRIMP spot analyses and laser ablation ICP-MS depth profiling show that thin (usually < 10 μm) younger (3,660–3,590 Ma and Neoarchaean) shells of lower Th/U metamorphic zircon are present on these 3,850 Ma zircons. Several samples with this simple zircon population occur on islands near Akilia. In contrast, migmatites usually contain more complex zircon populations, with often more than one generation of igneous zircon present. Additional zircon dating of banded gneisses across the Complex shows that samples with c. 3,850 Ma igneous zircon are not just a phenomenon restricted to Akilia and adjacent islands. For example, migmatites from Itilleq (c. 65 km from Akilia) contain variable amounts of oscillatory-zoned 3,850 Ma and 3,650 Ma zircon, interpreted, respectively, as the rock age and the time of crustal melting under Eoarchaean metamorphism. With only 110–140 ppm Zr in the tonalites and likely magmatic temperatures of >850°C, zircon solubility–melt composition relationships show that they were only one-third saturated in zircon. Any zircon entrained in the precursor magmas would thus have been highly soluble. Combined with the cathodoluminesence imaging, this demonstrates that the c. 3,850 Ma oscillatory zoned zircon crystallised out of the melt and hence gives a magmatic age. Thus the rare well-preserved tonalites and palaeosome in migmatites testify that c. 3,850 Ma quartzo–feldspathic rocks are a widespread (but probably minor) component in the Itsaq Gneiss Complex. C. 3,850 Ma zircon with negative Eu anomalies (showing growth in felsic systems) also occurs as detrital grains in rare c. 3,800 Ma metaquartzites and as inherited grains in some 3,660 Ma granites (sensu stricto). These demonstrate that still more c. 3,850 Ma rocks were present, but were recycled into Eoarchaean sediments and crustally derived granites. The major and trace element characteristics (e.g. LREE enrichment, HREE depletion, low MgO) of the best-preserved c. 3,850 Ma rocks are typical of Archaean TTG suites, and thus argue for crust formation processes involving important contributions from melting of hydrated mafic crust to the earliest Archaean. Five c. 3,850 Ma tonalites were selected as the best preserved on the basis of field criteria and zircon petrology. Four of these samples have overlapping initial ɛ_Nd (3,850 Ma) values from +2.9 to +3.6± 0.5, with the fourth having a slightly lower value of +0.6. These data provide additional evidence for a markedly LREE-depleted early terrestrial mantle reservoir. The role of c. 3,850 Ma crust should be considered in interpreting isotope signatures of the younger (3,800–3,600 Ma) rocks of the Itsaq Gneiss Complex. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Modeling Seasonal Variations in the Intensity of Internal Gravity Waves in the Lower Thermosphere

Izvestiya, Atmospheric and Oceanic Physics - An updated numerical model of the propagation of a set of spectral harmonics of internal gravity waves (IGWs) in the inhomogeneous atmosphere from the...     

Precambrian mafic magmatism in south Indian granulite terrain

South Indian granulite terrain had witnessed significant part of Precambrian mafic igneous activity in the form of episodic mafic dyke intrusions of the Palaeoproterozoic period. Strike trends of these dykes are not uniform over the region and the dykes are generally fresh, massive, black dolerites except in the Bhavani shear zone bordering the southern fringes of Nilgiri massif. In Agali-Coimbatore area of our study in the western Bhavani shear zone, the dykes appear to be penecontemporaneous with shearing. Isotopic data place age of Agali-Coimbatore dyke intrusions at about 2.1 Ga. The age of these dykes is significant to constrain an early Palaeoproterozoic age for major shearing event in the Bhavani shear zone. Other dyke emplacement ages are placed at about 1.8 Ga and 1.65 Ga based on the Ar/Ar and K-Ar isotopic results of dykes in Dharmapuri and Tiruvannamalai areas. Older ages comparable to those of the Dharwar craton are not known and in this respect future isotopic dating is vital. Geochemically, these dykes are quartz/hypersthene normative subalkalic tholeiites. An attempt is made here to provide insights into the general petrogenetic history of the Precambrian dykes. Compositional trends are explained by the fractional crystallization of ferromagnesian phases and plagioclase control is conspicuous at the advanced stages of fractionation. Geochemical characteristics suggest that the dykes have tapped Fe-rich non-pyrolite mantle sources with LIL and LREE enrichment as in many continental basalts. The data suggest that role of crustal contamination is limited in petrogenesis; crustal signatures are noticed in the more mafic end members formed in early stage of evolution suggesting that contamination was temperature controlled with most primitive high temperature magmas being most vulnerable to the process. Nd-Sr isotopic data, at present restricted to Agali-Coimbatore dykes, suggest that Palaeoproterozoic magmas tapped subcontinental lithosphere that may have stabilized in the Archaean times at about 3 Ga during the major crustal building activity in the shield region. Further work coupled with isotopic and mineral chemistry will improve our knowledge on the petrological evolution of the dyke magmas and mafic magmatism in general.     

Soil moisture dynamics under different land uses on karst hillslope in northwest Guangxi,China

Temporal and spatial dynamics of soil moisture are little known on karst hillslope with shallow soil in subtropical region. The objectives of this paper were (1) to investigate the temporal dynamics of soil moisture at depth of 0–10 cm under different land uses; and (2) to understand the relationship between mean and coefficient of variation (CV) of moisture contents on karst hillslope in northwest Guangxi, China. Soil moisture contents had a moderate variation (CV 17.5–30.3%) over an 8-month period and they had a significant difference among different land uses at the 0.01 level with a decreasing order: native scrubland > abandoned cropland and sloping cropland > economic forestland. There were higher mean and lower CV of moisture contents in rainy season than those in drought season. Mean and CV of moisture contents had a significant negative linear relationship except in abandoned cropland with higher soil and vegetation heterogeneity. This suggested that spatial variability of soil moisture within sampling sites would decrease when soils were wet and increase when soils were arid. Compared with rainy season, more soil samples may be needed and the interval for sampling should be shortened in drought season. Such information provided some insights to better understand the dynamics and variability of soil moisture at a larger scale in karst region of southwest China.     

Experimental modeling of the interaction of subducted carbonates and sulfur with mantle silicates

Experimental studies in the system Fe,Ni–olivine–carbonate–S (P = 6.3 GPa, T = 1050–1550°C, t = 40–60 h) aimed at modeling of the interaction of subducted carbonates and sulfur with rocks of the silicate mantle and at investigation of the likely mechanism of the formation of mantle sulfides were performed. It is shown that an association of olivine + orthopyroxene + magnesite + pyrite coexisting with a sulfur melt/fluid with dissolved Fe, Ni, and O is formed at T ≤ 1250°C. An association of low-Fe olivine, orthopyroxene, and magnesite and two immiscible melts of the carbonate and S–Fe–Ni–O compositions are formed at T ≥ 1350°C. It is shown that the reduced S-bearing fluids may transform silicates and carbonates, extract metals from the solid-phase matrix, and provide conditions for generation of sulfide melts.     

Identifying influence areas with connectivity analysis – application to the local perturbation of heterogeneity distribution for history matching

Numerical representations of a target reservoir can help to assess the potential of different development plans. To be as predictive as possible, these representations or models must reproduce the data (static, dynamic) collected on the field. However, constraining reservoir models to dynamic data – the history-matching process – can be very time consuming. Many uncertain parameters need to be taken into account, such as the spatial distribution of petrophysical properties. This distribution is mostly unknown and usually represented by millions of values populating the reservoir grid. Dedicated parameterization techniques make it possible to investigate many spatial distributions from a small number of parameters. The efficiency of the matching process can be improved from the perturbation of specific regions of the reservoir. Distinct approaches can be considered to define such regions. For instance, one can refer to streamlines. The leading idea is to identify areas that influence the production behavior where the data are poorly reproduced. Here, we propose alternative methods based on connectivity analysis to easily provide approximate influence areas for any fluid-flow simulation. The reservoir is viewed as a set of nodes connected by weighted links that characterize the distance between two nodes. The path between nodes (or grid blocks) with the lowest cumulative weight yields an approximate flow path used to define influence areas. The potential of the approach is demonstrated on the basis of 2D synthetic cases for the joint integration of production and 4D saturation data, considering several formulations for the weights attributed to the links.     

An approximate analytical solution for measuring air permeability of asphalt samples partially saturated with water

In coastal areas, abnormally high pressure may be caused by the tide-induced water table variation under extensive pavements, particularly during rainfall. To simulate the rainfall infiltration effects on the air permeability of asphalt pavements in coastal area, column-shaped asphalt sample was fixed in the upper part of a steel cylinder with its upper surface saturated with ponding water (depth < 5 mm) and open to the atmosphere. The cylinder’s lower part formed an air chamber. The chamber was pressurized and then the air therein was released naturally through the sample. The pressure variation with time in the chamber was recorded for analysis. Based on the Green–Ampt piston model for the surface water infiltration, an approximate analytical solution was derived to describe the pressure–time relationship in the chamber. A new parameter called the escape pressure was introduced to describe the air pressure needed for the chamber air to break through the capillary pressure induced by the ponding water. The analytical solution gave good estimations of both the escape pressures and the harmonic averages of the permeabilities of the wet and dry parts of 14 samples in the sense that excellent fittings were obtained between the observed and predicted air pressures in the air chamber. The estimated escape pressure ranges from 0.0 to 1.74 kPa. The harmonic average of the permeabilities of the wet and dry parts is 5–94% of the dry sample’s permeability.     

Theoretical density distribution of plasma streaming around a cylinder

The distribution of plasma density around a metallic cylinder in a collisionfree plasma is determined by the system of Vlassov equations and the Poisson equation. The plasma may have a streaming velocity perpendicular to the cylinder axis producing a wake behind the cylinder. In the region outside a thin double layer at the cylinder surface the problem allows some simplification. Then numerical solutions become possible-even if the streaming velocity is low. Some numerical solutions are presented. A magnetic field produces an asymmetric shape of the wake, if the gyration radius of the ion becomes comparable to the cylinder radius.