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.     

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.     

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.     

Mars: Experimental study of albedo changes caused by dust fallout

An experimental study was undertaken to determine how the spectral and photometric properties of representative Martian areas are affected by fallout of atmospheric dust suspended during dust forms. A laboratory apparatus was used to simulate the uniform fallout and deposition of particles 1 to 5 μm in diameter. Spectral measurements from 0.4- to 1.2-μm wavelengths and photometric measurements at several wavelengths were made for a number of Mars-analog materials before and after deposition of 6 × 10^?5 to 1.5 × 10^?3 g/cm² of simulated fallout. These results indicate that the spectral and photometric properties of Martian regions can be affected significantly even by minute amounts of fallout. For instance, the reflectance at 0.56 μm of an average dark area will increase by 35% after deposition of only 9 × 10^?5 g/cm², and by 70% after deposition of 1.5 × 10^?4 g/cm². Thus the fallout from even one dust storm season (～2 × 10^?3 g/cm²) is sufficient to change significantly the spectral and photometric characteristics of the substrate material, if the fallout were ubiquitous over the surface and if no competing processes of dust removal from surface grains occured.     

Local dynamical equivalence of certain food webs

An important challenge in theoretical ecology is to find good coarse-grained representations of complex food webs. Here, we use the approach of generalized modeling to show that it may be possible to formulate a coarse-graining algorithm that conserves the local dynamics of the model exactly. We show examples of food webs with a different number of species that have exactly identical local bifurcation diagrams. Based on these observations, we formulate a conjecture governing which populations of complex food webs can be grouped together into a single variable without changing the local dynamics. As an illustration, we use this conjecture to show that chaotic regions generically exist in the parameter space of a class of food webs with more than three trophic levels. While our conjecture is at present only applicable to relatively special cases, we believe that its applicability could be greatly extended if a more sophisticated mapping of parameters were used in the model reduction.     

On the application of the mode-acceleration method to structural engineering problems

Mode-superposition has been extensively used in computing the dynamic response of complex structures. Two versions of mode-superposition, namely the mode-displacement method and the mode-acceleration method, have been employed. The present paper summarizes the results of a systematic study comparing the accuracy of the mode-displacement and mode-acceleration methods when applied to structures with various levels of damping or various excitation frequencies. The paper also discusses several details concerning the implementation of the mode-acceleration method.     

Response of the overturning circulation to high-latitude fresh-water perturbations in the North Atlantic

Studies have suggested that sea-ice cover east and west of Greenland fluctuates out-of phase as a part of the Atlantic decadal climate variability, and greater changes are possible under global warming conditions. In this study, the response of the Atlantic meridional overturning circulation (MOC) to the distribution of surface fresh-water flux is explored using a global isopycnal ocean model. An Arctic ice related fresh-water flux of 0.1 Sv entering the Nordic Seas is shown to reduce the maximum overturning by 1 to 2 Sv (10⁶ m³ s^–1). A further decrease of 3 to 5 Sv in the MOC is observed when the fresh-water flux is shifted from the Fram Strait to the southern Baffin Bay area. Surprisingly, the salinity in much of the upper Nordic Seas actually increases when the Arctic fresh-water source is the strongest there, as a result of enhanced global overturning. It reflects the great influence of Labrador Sea convection on this models MOC. By applying a weaker surface fresh-water transport perturbation (0.02 Sv) on the Baffin Bay area and therefore perturbing the Labrador Sea Water (LSW) formation, we have also investigated the interaction between the overflows across the Greenland–Scotland Ridge and the LSW and find that, with the same surface forcing conditions in the Nordic Seas, volume transport of the overflows weakens when the LSW formation intensifies.     

Oscillations of galactic cosmic rays and solar indices before the arrival of relativistic solar protons

Using modern wavelet analysis techniques, we have made an attempt to search for oscillations of intensity of galactic cosmic rays (GCR), sunspot numbers (SS) and magnitudes of coronal index (CI) implying that the time evolution of those oscillations may serve as a precursor of Ground Level Enhancements (GLEs) of solar cosmic rays (SCR). From total number of 70 GLEs registered in 1942–2006, the four large events — 23 February 1956, 14 July 2000, 28 October 2003, and 20 January 2005 — have been chosen for our study. By the results of our analysis, it was shown that a frequency of oscillations of GCR decreases as time approaches to the event day. We have also studied a behaviour of common periodicities of GCR and SCR within the time interval of individual GLE. The oscillations of GLE occurrence rate (OR) at different stages of the solar activity (SA) cycle is of special interest. We have found some common periodicities of SS and CI in the range of short (2.8, 5.2, 27 and 60 days), medium (0.3, 0.5, 0.7, 1.3, 1.8 and 3.2 years) and long (4.6 and 11.0 years) periods. Short and medium periodicities, in general, are rather concentrated around the maxima of solar cycles and display the complex phase relations. When comparing these results with the behaviour of OR oscillations we found that the period of 11 years is dominating (controlling); it is continuous over the entire time interval of 1942–2006, and during all this time it displays high synchronization and clear linear ratios between the phases of oscillations of η, SS and CI. It implies that SCR generation is not isolated stochastic phenomena characteristic exclusively for chromospheric and/or coronal structures. In fact, this process may have global features and involve large regions in the Sun’s atmosphere.