Comparison of surface seismic sources at the CO2SINK site, Ketzin, Germany

In 2004 three seismic surface sources (VIBSIST, accelerated weight drop and MiniVib) were tested in a pilot study at the Ketzin test site, Germany, a study site for geological storage of CO₂ (EU project CO₂SINK). The main objectives of this pilot study were to 1) evaluate the response of the Ketzin site to reflection seismics, especially at the planned injection depth, 2) test different acquisition parameters and 3) use the results to guide the planning of the 3D survey. As part of these objectives, we emphasize the source performance comparison in this study. The sources were tested along two perpendicular lines of 2.4 km length each. Data were acquired by shooting at all stations (source and receiver spacing of 20 m) on both lines, allowing common‐midpoint stacked sections to be produced. The sources' signal characteristics based on signal‐to‐noise ratio, signal penetration and frequency content of raw shot records were analysed and stacked sections were compared. The results show that all three surface sources are suitable for reflection seismic studies down to a depth of about 1 km and provide enough bandwidth for resolving the geological targets at the site, i.e., the Weser and Stuttgart Formations. Near surface conditions, especially a thick weathering layer present in this particular area, strongly influence the data quality, as indicated by the difference in reflectivity and signal‐to‐noise ratio of the two common‐midpoint lines. The stacked sections of the MiniVib source show the highest frequency signals down to about 500 ms traveltime (approximately 500 m depth) but also the shallowest signal penetration depth. The VIBSIST source generates signals with the highest signal‐to‐noise ratio and greatest signal penetration depth of the tested sources. In particular, reflections below 900 ms (approximately 1 km depth) are best imaged by the VIBSIST source. The weight drop performance lies in between these two sources and might be recommended as an appropriate source for a 3D survey at this site because of the shorter production time compared to the VIBSIST and MiniVib sources.     

Petrogenetic grid for ultrahigh-pressure metamorphism in the model system CaO-MgO-SiO2-CO2-H2O

Abstract Petrogenetic grids for ultrahigh-pressure (UHP) metamorphism were calculated at different X_co2 conditions in the model system CaO-MgO-SiO₂-CO₂-H₂O involving coesite (Co), diopside (Di), dolomite (Do), enstatite (En), forsterite (Fo), magnesite (Ms), quartz (Qz), talc (Tc), tremolite (Tr) using a published internally consistent thermodynamic data set. Two P-T grids at X_co2= 0.01 and 0.5 are described. In the calculated P-T grid at X_co2= 0.01, four out of 10 stable invariant points, Co-En-Ms-Tc, Co-Di-En-Tc-Tr, Co-Di-Ms-Tc-Tr and Di-En-Ms-Tc-Tr lie within the stability field of coesite. If the fluid phase has X_co2= 0.5, no invariant point is stable under UHP conditions. Some magnesite-bearing assemblages are stabilized by the following three reactions: Di + Ms = Do + Fo + CO₂, Ms + Tr = Do + Fo + CO₂+ H₂O and Ms + Tc = Fo+ CO₂+ H₂O at X_co2= 0.01 and by reaction Ms + Tc = Fo + CO₂+ H₂O together with these three at X_co2= 0.5. Ten possible UHP assemblages for mafic and ultramafic compositions at very low X_co2 conditions include the following: Co-Do-Ms, Co-Di-Ms, Co-Di-Tc, Di-Ms-Tc, Di-En-Tc-, Di-En-Ms, Co-Di-En, Di-En-Fo, Di-Fo-Ms, Di-Do-Fo. Among them, talc-bearing assemblages are restricted to X_co2 < 0.02 and their high-P limit is 31.7 kb (749°C) at X_co2= 0.01. Dolomite-magnesite-silica assemblages have large P-T stability fields even if X_co2 is as low as 0.1, and could occur in cold subduction zones with very low geothermal gradients. Reported UHP coesite-dolomite assemblage is restricted only to a calc-silicate rock interlayered with marble where X_co2 is relatively higher; no such assemblage appears for mafic and ultramafic rocks with low X_co2 evidenced by the occurrence of diopside (or omphacite) at the expense of dolomite + coesite. The effect of X_co2 on the stability of coesite-dolomite-magnesite, diopside-enstatite-magnesite, diopside-talc assemblages is examined and the occurrence of coesite-dolomite, magnesite-bearing and talc-bearing assemblages in the Dabie UHP rocks are interpreted by employing the calculated P-T grids.     

Application of the continuous wavelet transform on seismic data for mapping of channel deposits and gas detection at the CO2SINK site, Ketzin, Germany

Conventional seismic data are band limited and therefore, provide limited geological information. Every method that can push the limits is desirable for seismic data analysis. Recently, time‐frequency decomposition methods are being used to quickly extract geological information from seismic data and, especially, for revealing frequency‐dependent amplitude anomalies. Higher frequency resolution at lower frequencies and higher temporal resolution at higher frequencies are the objectives for different time‐frequency decomposition methods. Continuous wavelet transform techniques, which are the same as narrow‐band spectral analysis methods, provide frequency spectra with high temporal resolution without the windowing process associated with other techniques. Therefore, this technique can be used for analysing geological information associated with low and high frequencies that normally cannot be observed in conventional seismic data. In particular, the continuous wavelet transform is being used to detect thin sand bodies and also as a direct hydrocarbon indicator. This paper presents an application of the continuous wavelet transform method for the mapping of potential channel deposits, as well as remnant natural gas detection by mapping low‐frequency anomalies associated with the gas. The study was carried out at the experimental CO₂ storage site at Ketzin, Germany (CO₂SINK). Given that reservoir heterogeneity and faulting will have significant impact on the movement and storage of the injected CO₂, our results are encouraging for monitoring the migration of CO₂ at the site. Our study confirms the efficiency of the continuous wavelet transform decomposition method for the detection of frequency‐dependent anomalies that may be due to gas migration during and after the injection phase and in this way, it can be used for real‐time monitoring of the injected CO₂ from both surface and borehole seismics.     

High-temperature elasticity of the fluorite-structure compounds CaF2, SrF2 and BaF2

The elastic moduli of single-crystal CaF₂, SrF₂ and BaF₂ have been determined by the ultrasonic pulse superposition technique as a function of temperature from T = 298 to T = 650°K. These new data are consistent with other data obtained by ultrasonic pulse techniques in the region of room temperature and are superior to previous high-temperature data from resonance experiments. The elastic moduli (c) are represented by quadratic functions in T over the experimental temperature range with the curvature in the same sense for all the moduli. Evaluation of the temperature derivatives of the elastic moduli at constant volume indicates that the dominant temperature effect is extrinsic for (?K_S/?T)_P and intrinsic for (?μ/?T)_P, where K_S and μ are the isotropic bulk and shear moduli, respectively. For the series CaF₂SrF₂BaF₂, |(?c/?T)_p| decreases with increasing molar volume for all moduli; however there are no theoretical or empirical grounds on which to derive a simple relationship between (?c/?T)_P and crystallographic parameters.     

Melting relationships in CaO-CO2 and MgO-CO2 to 36 kilobars with comments on CO2 in the mantle

The melting curves of CaCO₃ and MgCO₃ have been extended to pressures of 36 kb by experiments in piston-cylinder apparatus. At 30 kb, the melting temperatures of calcite and magnesite are 1610°C and 1585°C, respectively. New data for the magnesite dissociation reaction permit the location of an invariant point for the assemblage magnesite + periclase + liquid + vapor near 26 kb-1550°C. New data are also presented for the calcite-aragonite transition at 800°C, 950°C and 1100°C. At pressures above 36–50 kb, calcite and magnesite melt at temperatures lower than the solidus of dry mantle peridotite. Natural and experimental evidence suggests that carbon dioxide in the Earth's mantle could be present in a variety of forms: (a) a free vapor phase, (b) vapor dissolved in silicate magma, (c) crystalline carbonate, (d) carbonatite liquid, (e) carbon-bearing silicate analogs, or (f) carbonato-silicates (such as scapolite, spurrite, tilleyite, and related compounds).     

High-pressure X-ray diffraction studies on β- and γ-Mg2SiO4

Pressure effects on the lattice parameters of β- and γ-Mg₂SiO₄ have been measured at room temperature and at pressures up to 100 kbar using a multi-anvil high-pressure X-ray diffraction apparatus. The volume changes (ΔV/V₀) at 90 kbar are 5.4 · 10^?2 and 4.2 · 10^?2 for β- and γ-Mg₂SiO₄, respectively. Isothermal bulk moduli at zero pressure have been calculated from least-square fits of the data to straight lines. They turn out to be 1.66 ± 0.4 and 2.13 ± 0.1 Mbar for β- and γ-Mg₂SiO₄, respectively. The α → γ transition obeys Wang's linear V_φ?ρ relation but the α → β transition does not.     

Correlations between13C and ΣCO2 in surface waters and atmospheric CO2

¹³C and ΣCO₂ data from the North and South Atlantic, the Antarctic, and the North and South Pacific are given. The δ¹³C of the ΣCO₂ in the deep water (～3000m) decreases from 1.7‰ in the North Atlantic to ?0.10‰ in the North Pacific. This change is attributed to the addition of about 158 μmoles of CO₂ per kg of seawater. The in-situ oxidation of organic matter accounts for 83% of this increase in ΣCO₂, while the remainder is attributed to dissolution of calcium carbonate.The δ¹³C of the dissolved CO₂ in mid-latitude surface water samples is controlled by a quasi-steady-state equilibrium with atmospheric CO₂ at a mean temperature of 16°C. The δ¹³C and ΣCO₂ values of Antarctic surface water samples suggest that these waters are derived from a mixture of North Atlantic deep water and equilibrated surface water.     

Time-lapse processing of 2D seismic profiles with testing of static correction methods at the CO2 injection site Ketzin (Germany)

The Ketzin project provides an experimental pilot test site for the geological storage of CO₂. Seismic monitoring of the Ketzin site comprises 2D and 3D time-lapse experiments with baseline experiments in 2005. The first repeat 2D survey was acquired in 2009 after 22 kt of CO₂ had been injected into the Stuttgart Formation at approximately 630 m depth. Main objectives of the 2D seismic surveys were the imaging of geological structures, detection of injected CO₂, and comparison with the 3D surveys. Time-lapse processing highlighted the importance of detailed static corrections to account for travel time delays, which are attributed to different near-surface velocities during the survey periods. Compensation for these delays has been performed using both pre-stack static corrections and post-stack static corrections. The pre-stack method decomposes the travel time delays of baseline and repeat datasets in a surface consistent manner, while the latter cross-aligns baseline and repeat stacked sections along a reference horizon.Application of the static corrections improves the S/N ratio of the time-lapse sections significantly. Based on our results, it is recommended to apply a combination of both corrections when time-lapse processing faces considerable near-surface velocity changes. Processing of the datasets demonstrates that the decomposed solution of the pre-stack static corrections can be used for interpretation of changes in near-surface velocities. In particular, the long-wavelength part of the solution indicates an increase in soil moisture or a shallower groundwater table in the repeat survey.Comparison with the processing results of 2D and 3D surveys shows that both image the subsurface, but with local variations which are mainly associated to differences in the acquisition geometry and source types used. Interpretation of baseline and repeat stacks shows that no CO₂ related time-lapse signature is observable where the 2D lines allow monitoring of the reservoir. This finding is consistent with the time-lapse results of the 3D surveys, which show an increase in reflection amplitude centered around the injection well. To further investigate any potential CO₂ signature, an amplitude versus offset (AVO) analysis was performed. The time-lapse analysis of the AVO does not indicate the presence of CO₂, as expected, but shows signs of a pressure response in the repeat data.     

Melting of the silica isotypes SiO 2, BeF 2 and GeO 2 at elevated pressures

The melting curves of the structural analogues SiO ₂, BeF ₂ and GeO ₂ have been studied at pressures ?40 kbar in a piston-cylinder apparatus. The initial slopes dT_m/dP of the β-quartz-liquid boundaries for SiO ₂ and BeF ₂ are ～35° while the slope of the rutile-liquid boundary for GeO ₂ is approximately 32°C/kbar. These large values of dT/dP reflect the unusually low entropies of fusion for these compounds in which strong structural similarities exist between the crystalline phases and the melt. Implications for the extended phase diagram of silica are discussed and it is concluded that either: (1) a maximum exists on the coesite melting curve, or (2) estimates of the melting temperature of stishovite need to be revised upwards.     

High-pressure disproportionation of Co2SiO4 spinel and implications for Mg2SiO4 spinel

Co₂SiO₄ spinel has been found to disproportionate into its isochemically mixed oxides with rocksalt and rutile structures at pressures between 170 and 190 kbar and temperatures between 1400 and 1800°C in a diamond-anvil press. The exact disproportionation pressure is not certain due to transient increases in pressure during the local and rapid heating by a continuous YAG laser. The slope of the phase boundary between the spinel phase and the mixed oxides is calculated to be?33 ± 20bar/deg. This negative slope is consistent with the observed anomalously large entropy of CoO (relative to its isostructural oxides) in entropy vs.(MV)^?1/2 systematics, whereM is the formula weight andV the molar volume. The sign of the slope for a phase boundary in the disproportionation of spinel depends on the values of entropy of the rocksalt oxides as well as the inverse character exhibited in the spinel phases. The normal entropy of MgO suggests that the phase boundary for the disproportionation of Mg₂SiO₄ spinel has positive slope.     

High-pressure phases of Co2GeO4, Ni2GeO4, Mn2GeO4 and MnGeO3: implications for the germanate-silicate modeling scheme and the Earth's mantle

In a diamond-anvil press coupled with YAG laser heating, the spinels of Co₂GeO₄ and Ni₂GeO₄ have been found to disproportionate into their isochemical oxide mixtures at about 250 kbar and 1400–1800°C in the same manner as their silicate analogues. At about the same P-T conditions MnGeO₃ transforms to the orthorhombic perovskite structure (space group Pbnm); the lattice parameters at room temperature and 1 bar are a₀ = 5.084 ± 0.002, b₀ = 5.214 ± 0.002, and c₀ = 7.323 ± 0.003Å with Z = 4 for the perovskite phase. The zero-pressure volume change associated with the ilmenite-perovskite phase transition in MnGeO₃ is ?6.6%. Mn₂GeO₄ disproportionates into a mixture of the perovskite phase of MnGeO₃ plus the rocksalt phase of MnO at P = 250kbar and T = 1400–1800°C. The concept of utilizing germanates as high-pressure models for silicates is valid in general. The results of this study support the previous conclusion that the lower mantle comprises predominantly the orthorhombic perovskite phase of ferromagnesian silicate.     

武汉和阿德莱德大气太阴N2潮及M2潮的对比分析

利用阿德莱德（35°S,138°E）和武汉（30.6°N,114.5°E）的流星雷达观测数据首次给出了80~100 km高度上周期为12.66太阳时的大气太阴半日潮汐（N₂潮）的季节、高度、年度变化及其与周期为12.42太阳时的大气太阴半日潮汐（M₂潮）的对比分析.分析结果表明：武汉和阿德莱德的N₂潮和M₂潮均有明显的季节、高度和年度变化.N₂潮与M₂潮的幅度比值大于其引力势之比0.191,在某些年份的不同季节和高度上,N₂潮的幅度甚至大于M₂潮的幅度.大多幅度之比接近或超过N₂潮和M₂潮引力势之比的2倍.中低热层的大气太阴N₂潮汐值得关注.