Bayesian estimation of reservoir properties—effects of uncertainty quantification of 4D seismic data

This paper shows a history matching workflow with both production and 4D seismic data where the uncertainty of seismic data for history matching comes from Bayesian seismic waveform inversion. We use a synthetic model and perform two seismic surveys, one before start of production and the second after 1 year of production. From the first seismic survey, we estimate the contrast in slowness squared (with uncertainty) and use this estimate to generate an initial estimate of porosity and permeability fields. This ensemble is then updated using the second seismic survey (after inversion to contrasts) and production data with an iterative ensemble smoother. The impact on history matching results from using different uncertainty estimates for the seismic data is investigated. From the Bayesian seismic inversion, we get a covariance matrix for the uncertainty and we compare using the full covariance matrix with using only the diagonal. We also compare with using a simplified uncertainty estimate that does not come from the seismic inversion. The results indicate that it is important not to underestimate the noise in seismic data and that having information about the correlation in the error in seismic data can in some cases improve the results.     

Influence of terrestrial vegetation on leaf wax δD of Holocene lake sediments

The hydrogen isotope composition of terrestrial plant leaf wax in sediments is increasingly used as a paleoclimatic indicator. Modern calibration studies suggest that paleoclimatic interpretation of leaf wax δD values requires consideration of the differences in the apparent fractionation of hydrogen isotopes among different groups of plants. However, it is not common that paleoecological data are used to help interpret leaf wax δD profiles. Here we assess the relative importance of factors influencing millennial-scale shifts in δD values of n-alkanoic acids at Steel Lake (Minnesota, USA), an extensively studied site with independent records of vegetation composition, δD of input water to the lake, and evaporation. The δD values of the n-C₂₈ alkanoic acid (δD_C28) vary between −190 and −168‰, and do not correlate with δD of input water or the extent of evaporation. However, δD_C28 is negatively correlated with the δ¹³C values of the n-C₂₈ alkanoic acid (δ¹³C_C28). The correlation, along with pollen assemblage and carbonate δ¹³C records, suggests that Holocene shifts between forest and grassland and/or in the water use efficiency of C₃ plants influenced the stratigraphic variation in leaf wax δD. Thus, paleoecological information, such as that inferred from pollen assemblages and carbon isotopes of plant-derived compounds, may aid paleoclimatic interpretation of leaf wax δD profiles from lake sediments.     

Leaf wax n-alkane δD values of field-grown barley reflect leaf water δD values at the time of leaf formation

Leaf wax n-alkanes from barley (Hordeum vulgare) from a field in Switzerland exhibited changes in δD values on the order of 20‰ over a growing season, while source water (soil water) and leaf water varied by 40‰. Additionally the seasonal variability in δD values of leaf wax n-alkanes of different barley leaves can only be found across different leaf generations (i.e. leaves that were produced at different times during the growing season) while n-alkane δD values did not vary significantly within a leaf generation. Interestingly, δD values of n-alkanes correlated best with the δD values of leaf water at midday of the sampling day but showed no significant correlation with soil water (e.g. precipitation) δD values. These results provide empirical evidence that leaf wax δD values record leaf water enrichment, and therefore integrate the isotopic effects of precipitation and evapotranspiration. Our results show that leaf wax n-alkane δD values from grasses are ‘locked in’ early during leaf development and hence record the environmental drivers of leaf water enrichment, such as vapor pressure deficit (VPD). Our data have important implications for the interpretation of paleorecords of leaf wax δD. We suggest that leaf wax n-alkane δD values from sedimentary records could be used to estimate changes in the degree of leaf water enrichment and hence VPD.     

Urgent Proposal for Re-consideration of the D/C Boundary Stratatotype

Despite the transition from Siphonodella Praesulcata to S·sulcata“demonstrated”in the excellent biostratigraphic-sedimentologic treatment of the La Serre section by FLAJS＆FEIST(1988),we remain unconvinced of the scientific validity of this section as thestratotype for the D/C Boundary in the compromise package proposal formulated atCourtmacsherry in May1988.Therefore we,as concerned conodont workers and titular     

2D thermomechanical modelling of continent–arc–continent collision

Arc–continent collision is a key process of continental growth through accretion of newly grown magmatic arc crust to older continental margin. We present 2D petrological–thermo-mechanical models of arc–continent collision and investigate geodynamic regimes of this process. The model includes spontaneous slab bending, dehydration of subducted crust, aqueous fluid transport, partial melting of the crustal and mantle rocks and magmatic crustal growth stemming from melt extraction processes. Results point to two end-member types of subsequent arc–continent collisional orogens: (I) orogens with remnants of accretion prism, detached fragments of the overriding plate and magmatic rocks formed from molten subducted sediments; and (II) orogens mainly consisting of the closed back-arc basin suture, detached fragments of the overriding plate with leftovers of the accretion prism and quasi insignificant amount of sediment-derived magmatic rocks. Transitional orogens between these two endmembers include both the suture of the collapsed back-arc basin and variable amounts of magmatic production. The orogenic variability mainly reflects the age of the subducting oceanic plate. Older, therefore colder and denser oceanic plates trigger subduction retreat, which in turn triggers necking of the overriding plate and opening of a backarc basin in which new oceanic lithosphere is formed from voluminous decompression melting of the rising hot asthenosphere. In this case, subducted sediments are not heated enough to melt and generate magmatic plumes. On the other hand, young and less dense slabs do not retreat, which hampers opening of a backarc basin in the overriding plate while subducted sediments may reach their melting temperature and develop trans-lithospheric plumes. We have also investigated the influences of convergence rate and volcanic/plutonic rocks' ratio in newly forming lithosphere. The predicted gross-scale orogenic structures find similarities with some natural orogens, in particular with deeply eroded orogens such as the Variscides in the Bohemian Massif.     

Compound specific δD values of long chain n-alkanes derived from terrestrial higher plants are indicative of the δD of meteoric waters: Evidence from surface soils in eastern China

The stable carbon and hydrogen isotope composition of higher plant-derived long chain n-alkanes (δ¹³C_n-alkanes and δD_n-alkanes) from 45 surface soil samples (within well characterized vegetation zones) from eastern China (18°N–50°N) are reported. The weighted average δD_n-alkanes value for n-C₂₇, n-C₂₉ and n-C₃₁ in the samples and the annual average δD of meteoric water recorded at 12 weather stations proximal to the sampling sites show similar spatial variations. The δD of n-alkanes shows a gradual depletion in value with increasing latitude. The results demonstrate that, on a large spatial scale, the δD values of long chain n-alkanes derived from higher plants have the potential to record the δD of meteoric water, although many other factors can also influence the isotope values. There appears to be no apparent relationship between the δD of the n-alkanes extracted from the surface soil and the overlying vegetation type (i.e. forest/grassland or C₃/C₄ composition). Therefore, palaeoenvironmental studies utilizing δD_n-alkanes from higher plant-derived material in geological samples have the potential to provide additional information with regard to the past hydrological cycle.     

Local mineral equilibrium in metamorphism: development of D. Korzhinskii’s ideas (on 110th anniversary of D.S. Korzhinskii)

In 1950, D. Korzhinskii hypothesized that minerals in metamorphic reactions can be in local equilibrium. Proceeding from this idea, G. Fisher and R. Joesten later developed a realistic model of metamorphism kinetics. The model included diffusive mass transport in intergranular fluid controlled by the concentrations of components and chemical potential gradients, local-equilibrium reactions between minerals and pore fluid, mass balance of components, etc. The Korzhinskii–Fisher–Joesten model can provide clues to the duration of metamorphic events as inferred from mineral chemistry, free energy of phases, and rock textures and structures. The respective modeling of metamorphic reactions, with regard to spatial distribution of mineral grains, has implications for the transport and balance of chemical components. The balance of components (except volatiles) is restricted to domains of hundredth fractions of a cubic millimeter to several cubic meters. This may be the minimum size of an elementary domain in which minerals are in local equilibrium.     

3D Geological Modeling of the Xiangshan Uraniferous Volcanic , Jiangxi Province

正Objective Three-dimensional (3D) subsurface geological models have been gaining considerable attention wordwide. 3D geological surveys and modelling in the Xiangshan uraniferous volcanic basin, Le’an county, Jiangxi Province is a pilot project of the China Geological Survey. The Xiangshan volcanic basin is the third largest volcanic-type uranium ore field worldwide. This basin is composed of Lower Cretaceous rhyodacite and porphyritic lava with     

The compressibility of Fe- and Al-bearing phase D to 30 GPa

High-pressure in situ X-ray diffraction experiment of Fe- and Al-bearing phase D (Mg_0.89Fe_0.14Al_0.25Si_1.56H_2.93O₆) has been carried out to 30.5 GPa at room temperature using multianvil apparatus. Fitting a third-order Birch–Murnaghan equation of state to the P–V data yields values of V ₀ = 86.10 ± 0.05 ų; K ₀ = 136.5 ± 3.3 GPa and K′ = 6.32 ± 0.30. If K′ is fixed at 4.0 K ₀ = 157.0 ± 0.7 GPa, which is 6% smaller than Fe–Al free phase D reported previously. Analysis of axial compressibilities reveals that the c-axis is almost twice as compressible (K _c  = 93.6 ± 1.1 GPa) as the a-axis (K _a  = 173.8 ± 2.2 GPa). Above 25 GPa the c/a ratio becomes pressure independent. No compressibility anomalies related to the structural transitions of H-atoms were observed in the pressure range to 30 GPa. The density reduction of hydrated subducting slab would be significant if the modal amount of phase D exceeds 10%.     

Multi-dimensional Seismic Response Analysis of Base-Isolated Frame Structure with 3D Isolator

INTRODUCTIONAs a new and effective method for withstandingthe horizontal component of earthquake motion on abuilding structure, the base isolation technique hascome to the engineering application phase ( Hu,1988). However, for some important buildings andfoundation facilities in an earthquake region, it is im portant to consider the effects of 3D aseismatic analy sis, including the vertical component, due to themulti dimensional characteristics of earthquakes (Suand …     

Heterogeneous response of hydrogeological systems to the Izmit and Düzce (Turkey) earthquakes of 1999

At four sites in Turkey and Armenia the physico-chemical properties of thermal and mineral waters were monitored continuously during the Izmit and Düzce earthquakes that occurred along the North Anatolian fault in August and November 1999. The epicentral distances between the moment magnitude (M_w) 7.6 Izmit earthquake and the monitoring locations were 313, 488, 1,161, and 1,395 km. At the most distant site, the specific electrical conductivity of mineral water from a flowing artesian well dropped co-seismically and postseismically by 7%. No changes were observed at the other sites, although the estimated earthquake strains and peak ground accelerations are much higher. A similar pattern was observed after the Düzce earthquake, which happened three months after the Izmit event. The response of a hydrogeological system seems to depend on the site characteristics rather than on the nature of the earthquake. A hydrogeological model for the sensitive observation site farthest from the Izmit earthquake explains the observations in terms of a changed mixing ratio between two fluid components. Passing seismic waves may trigger a local pore-pressure increase according to the mechanism of advective overpressure. The preconditions for this mechanism, free gas bubbles in the aquifer in combination with a trap for rising bubbles, is probably not fulfilled by the other groundwater systems. Electronic Publication     

Intrusion versus inversion—a 3D density model of the southern rim of the Northwest German Basin

An unsolved problem of regional importance for both the evolution and structure of the Northwest German Basin is the existence or non-existence of the so-called Bramsche Massif. Explaining the nature of this massif and the cause of a related strong, positive Bouguer anomaly (Bramsche Anomaly) is critical. In the study described here, we tested an existing “intrusion model” against a newer “inversion model” in the southern Northwest German Basin. In the intrusion model, the strongly-positive Bouguer anomaly represents the gravity effect of an intrusion at depths between 6 and 10 km. More recent interpretations invoke tectonic inversion rather than intrusion to explain increased burial and the low level of hydrocarbon maturity found in boreholes. We tested these different interpretations by constructing 3D forward density models to 15 km depth. The intrusion model was updated and adjusted to incorporate recent data and we also modelled pre-Zechstein structures using different scenarios. The final model has a very good fit between measured and modelled gravity fields. Based on currently available seismic and structural models, as well as borehole density measurements, we show that the positive Bouguer anomaly cannot be modeled without a high-density, intrusive-like body at depth. However, further in-sight into the crustal structures of the Bramsche region requires more detailed investigations. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

D检验法

本文简短地介绍了用D检验法来检验大容量样本统计数据的正态性.为说明方法的应用效果,举了一个用于化探的实例.     

Effects of aridity and vegetation on plant-wax δD in modern lake sediments

We analyzed the deuterium composition of individual plant-waxes in lake sediments from 28 watersheds that span a range of precipitation D/H, vegetation types and climates. The apparent isotopic fractionation (ε_a) between plant-wax n-alkanes and precipitation differs with watershed ecosystem type and structure, and decreases with increasing regional aridity as measured by enrichment of ²H and ¹⁸O associated with evaporation of lake waters. The most negative ε_a values represent signatures least affected by aridity; these values were −125 ± 5‰ for tropical evergreen and dry forests, −130‰ for a temperate broadleaf forest, −120 ± 9‰ for the high-altitude tropical páramo (herbs, shrubs and grasses), and −98 ± 6‰ for North American montane gymnosperm forests. Minimum ε_a values reflect ecosystem-dependent differences in leaf water enrichment and soil evaporation. Slopes of lipid/lake water isotopic enrichments differ slightly with ecosystem structure (i.e. open shrublands versus forests) and overall are quite small (slopes = 0-2), indicating low sensitivity of lipid δD variations to aridity compared with coexisting lake waters. This finding provides an approach for reconstructing ancient precipitation signatures based on plant-wax δD measurements and independent proxies for lake water changes with regional aridity. To illustrate this approach, we employ paired plant-wax δD and carbonate-δ¹⁸O measurements on lake sediments to estimate the isotopic composition of Miocene precipitation on the Tibetan plateau.     

Strain memory of 2D and 3D rigid inclusion populations in viscous flows — What is clast SPO telling us?

We model the development of shape preferred orientation (SPO) of a large population of two- and three-dimensional (2D and 3D) rigid clasts suspended in a linear viscous matrix deformed by superposed steady and continuously non-steady plane strain flows to investigate the sensitivity of clasts to changing boundary conditions during a single or superposed deformation events. Resultant clast SPOs are compared to one developed by an identical initial population that experienced a steady flow history of constant kinematic vorticity and reached an identical finite strain state, allowing examination of SPO sensitivity to deformation path. Rotation paths of individual triaxial inclusions are complex, even for steady plane strain flow histories. It has been suggested that the 3D nature of the system renders predictions based on 2D models inadequate for applied clast-based kinematic vorticity gauges. We demonstrate that for a large population of clasts, simplification to a 2D model does provide a good approximation to the SPO predicted by full 3D analysis for steady and non-steady plane strain deformation paths. Predictions of shape fabric development from 2D models are not only qualitatively similar to the more complex 3D analysis, but they display the same limitations of techniques based on clast SPO commonly used as a quantitative kinematic vorticity gauge. Our model results from steady, superposed, and non-steady flow histories with a significant pure shearing component at a wide range of finite strain resemble predictions for an identical initial population that experienced a single steady simple shearing deformation. We conclude that individual 2D and 3D clasts respond instantaneously to changes in boundary conditions, however, in aggregate, the SPO of a population of rigid inclusions does not reflect the late-stage kinematics of deformation, nor is it an indicator of the unique ‘mean’ kinematic vorticity experienced by a deformed rock volume.     

Pressure responses of portlandite and H–D isotope effects on pressure-induced phase transitions

The pressure responses of portlandite and the isotope effect on the phase transition were investigated at room temperature from single-crystal Raman and IR spectra and from powder X-ray diffraction using diamond anvil cells under quasi-hydrostatic conditions in a helium pressure-transmitting medium. Phase transformation and subsequent peak broadening (partial amorphization) observed from the Raman and IR spectra of Ca(OH)₂ occurred at lower pressures than those of Ca(OD)₂. In contrast, no isotope effect was found on the volume and axial compressions observed from powder X-ray diffraction patterns. X-ray diffraction lines attributable to the high-pressure phase remained up to 28.5 GPa, suggesting no total amorphization in a helium pressure medium within the examined pressure region. These results suggest that the H–D isotope effect is engendered in the local environment surrounding H(D) atoms. Moreover, the ratio of sample-to-methanol–ethanol pressure medium (i.e., packing density) in the sample chamber had a significant effect on the increase in the half widths of the diffraction lines, even at pressures below the hydrostatic limit of the pressure medium.     

Research on Ecostratigraphy of Changhsingian at Meishan Section D, Changxing, Zhejiang Province

On the basis of researches over several years on biostratigraphy of Changhsingian at Meishan Section D in Changxing County, the abundant materials were accumulated. This paper studies the escostratigraphy and establishes 7 community zones (CZ), in ascending order, Sinoplatysomus-Gernitzina-Clarkina subcarinata, Tapashanites-Nodosaria-Clarkina subcarinata,Crurithyris-Geinitzina-Clarkina subcarinata,Glomospira-Clarkina changxingensis, Colaniella-Clarkina changxingensis,Rotodiscoceras-Palaeofusulina-Hindeodus typicalis,Clarkina meishanensis-Hypophiceras.The habitat type (HT) is diffent from upper shallow sea to the lower part of the lower shallow sea.