Frequency‐domain waveform modelling and inversion for coupled media using a symmetric impedance matrix

To simulate the seismic signals that are obtained in a marine environment, a coupled system of both acoustic and elastic wave equations is solved. The acoustic wave equation for the fluid region simulates the pressure field while minimizing the number of degrees of freedom of the impedance matrix, and the elastic wave equation for the solid region simulates several elastic events, such as shear waves and surface waves. Moreover, by combining this coupled approach with the waveform inversion technique, the elastic properties of the earth can be inverted using the pressure data obtained from the acoustic region. However, in contrast to the pure acoustic and elastic cases, the complex impedance matrix for the coupled media does not have a symmetric form because of the boundary (continuity) condition at the interface between the acoustic and elastic elements. In this study, we propose a manipulation scheme that makes the complex impedance matrix for acoustic–elastic coupled media to take a symmetric form. Using the proposed symmetric matrix, forward and backward wavefields are identical to those generated by the conventional approach; thus, we do not lose any accuracy in the waveform inversion results. However, to solve the modified symmetric matrix, LDLT factorization is used instead of LU factorization for a matrix of the same size; this method can mitigate issues related to severe memory insufficiency and long computation times, particularly for large‐scale problems.     

Fluid injection monitoring using electrical resistivity tomography — five years of injection at Ketzin,Germany

Between the years 2008 and 2013, approximately 67 kilotons of CO₂ have been injected at the Ketzin site, Germany. As part of the geophysical monitoring programme, time‐lapse electrical resistivity tomography has been applied using crosshole and surface‐downhole measurements of electrical resistivity tomography. The data collection of electrical resistivity tomography is partly based on electrodes that are permanently installed in three wells at the site (one injection well and two observation wells). Both types of ERT measurements consistently show the build‐up of a CO₂‐related resistivity signature near the injection point. Based on the imaged resistivity changes and a petrophysical model, CO₂ saturation levels are estimated. These CO₂ saturations are interpreted in conjunction with CO₂ saturations inferred from neutron‐gamma loggings. Apart from the CO₂–brine substitution response in the observed resistivity changes, significant imprints from the dynamic behaviour of the CO₂ in the reservoir are observed.     

Characterisation and multifaceted anisotropy assessment of Corvio sandstone for geological CO2 storage studies

We present a comprehensive characterisation of the physical, mineralogical, geomechanical, geophysical, and hydrodynamic properties of Corvio sandstone. This information, together with a detailed assessment of anisotropy, is needed to establish Corvio sandstone as a useful laboratory rock‐testing standard for well‐constrained studies of thermo–hydro–mechanical–chemical coupled phenomena associated with CO₂ storage practices and for geological reservoir studies in general. More than 200 core plugs of Corvio sandstone (38.1 and 50 mm diameters, 2:1 length‐to‐diameter ratio) were used in this characterisation study, with a rock porosity of 21.7 ± 1.2%, dry density 2036 ± 32 kg m^?3, and unconfined compressive and tensile strengths of 41 ± 3.28 and 2.3 ± 0.14 MPa, respectively. Geomechanical tests show that the rock behaves elastically between ～10 and ～18 MPa under unconfined conditions with associated Young's modulus and Poisson's ratio of 11.8 ± 2.8 GPa and 0.34 ± 0.01 GPa, respectively. Permeability abruptly decreases with confining pressure up to ～10 MPa and then stabilises at ～1 mD. Ultrasonic P‐ and S‐wave velocities vary from about 2.8–3.8 km s^?1 and 1.5–2.4 km s^?1, respectively, over confining and differential pressures between 0.1 and 35 MPa, allowing derivation of associated dynamic elastic moduli. Anisotropy was investigated using oriented core plugs for electrical resistivity, elastic wave velocity and attenuation, permeability, and tracer injection tests. Corvio sandstone shows weak transverse isotropy (symmetry axis normal to bedding) of <10% for velocity and <20% for attenuation.     

Microseismic events enhancement and detection in sensor arrays using autocorrelation‐based filtering

Passive microseismic data are commonly buried in noise, which presents a significant challenge for signal detection and recovery. For recordings from a surface sensor array where each trace contains a time‐delayed arrival from the event, we propose an autocorrelation‐based stacking method that designs a denoising filter from all the traces, as well as a multi‐channel detection scheme. This approach circumvents the issue of time aligning the traces prior to stacking because every trace's autocorrelation is centred at zero in the lag domain. The effect of white noise is concentrated near zero lag; thus, the filter design requires a predictable adjustment of the zero‐lag value. Truncation of the autocorrelation is employed to smooth the impulse response of the denoising filter. In order to extend the applicability of the algorithm, we also propose a noise prewhitening scheme that addresses cases with coloured noise. The simplicity and robustness of this method are validated with synthetic and real seismic traces.     

Uncertainty analysis and probabilistic segmentation of electrical resistivity images: the 2D inverse problem

In this paper, we present the uncertainty analysis of the 2D electrical tomography inverse problem using model reduction and performing the sampling via an explorative member of the Particle Swarm Optimization family, called the Regressive‐Regressive Particle Swarm Optimization. The procedure begins with a local inversion to find a good resistivity model located in the nonlinear equivalence region of the set of plausible solutions. The dimension of this geophysical model is then reduced using spectral decomposition, and the uncertainty space is explored via Particle Swarm Optimization. Using this approach, we show that it is possible to sample the uncertainty space of the electrical tomography inverse problem. We illustrate this methodology with the application to a synthetic and a real dataset coming from a karstic geological set‐up. By computing the uncertainty of the inverse solution, it is possible to perform the segmentation of the resistivity images issued from inversion. This segmentation is based on the set of equivalent models that have been sampled, and makes it possible to answer geophysical questions in a probabilistic way, performing risk analysis.     

Impacts of Seasonal Pumping on Stream‐Aquifer Interactions in Miryang,Korea

Large agricultural fields in South Korea are located mostly on alluvial plains, where a significant amount of groundwater is used for heating of water‐curtain insulated greenhouses. Such greenhouses are commonly used for crop cultivation during the winter dry season from November to March. After use the groundwater is discharged directly into streams, causing groundwater depletion. A hydrogeological study was carried out in a typical agricultural area of this type, located on an alluvial aquifer near the Nakdong River. Groundwater levels, chemical characteristics, and temperatures from 68 observation wells were analyzed to determine the impacts of seasonal groundwater pumping on the groundwater system and stream‐aquifer interactions. Our results show that the groundwater system has not yet reached a state of dynamic equilibrium. Decades of excessive seasonal pumping have caused a gradual decline of groundwater levels, leading to groundwater depletion, especially in areas further from the river. Seasonal pumping has also significantly affected groundwater quality in the aquifer near the river. Groundwater temperature is decreasing (in this case a disadvantage), and saline groundwater is being diluted by induced recharge. The results of this study provide a basic outline for effective integrated water management that is widely applicable in South Korea.     

Geochemistry,and zircon U–Pb and molybdenite Re–Os geochronology of Jilongshan Cu–Au deposit,southeastern Hubei Province,China

The Jilongshan skarn Cu–Au deposit is located at the Jiurui ore cluster region in the southwestern part of the Middle–Lower Yangtze River valley metallogenic belt. The region is characterized by NW‐, NNW‐ and EW‐trending faults and the mineralization occurs at the contact of lower Triassic carbonate rocks and Jurassic granodiorite porphyry intrusions. The intrusives are characterized by SiO₂, K₂O, and Na₂O concentrations ranging from 61.66 to 67.8 wt.%, 3.29 to 5.65 wt.%, and 2.83 to 3.9 wt.%, respectively. Their A/CNK (A/CNK = n(Al₂O₃)/[n(CaO) + n(Na₂O) + n(K₂O)]) ratio, δEu, and δCe vary from 0.77 to 1.17, 0.86 to 1, and 0.88 to 0.96, respectively. The rocks show enrichment in light rare earth elements ((La/Yb)_N = 7.61–12.94) and large ion lithophile elements (LILE), and depletion in high field strength elements (HFSE), such as Zr, Ti. They also display a peraluminous, high‐K calc‐alkaline signature typical of intrusives associated with skarn and porphyry Cu–Au–Mo polymetallic deposits. Laser ablation inductively coupled plasma spectrometry (LA‐ICP‐MS) zircon U–Pb age indicates that the granodiorite porphyry formed at 151.75 ± 0.70 Ma. A few inherited zircons with older ages (677 ± 10 Ma, 848 ± 11 Ma, 2645 ± 38 Ma, and 3411 ± 36 Ma) suggest the existence of an Archaean basement beneath the Middle–Lower Yangtze River region. The temperature of crystallization of the porphyry estimated from zircon thermometer ranges from 744.3 °C to 751.5 °C, and 634.04 °C to 823.8 °C. Molybdenite Re–Os dating shows that the Jilongshan deposit formed at 150.79 ± 0.82 Ma. The metallogeny and magmatism are correlated to mantle–crust interaction, associated with the subduction of the Pacific Plate from the east. Copyright © 2013 John Wiley & Sons, Ltd.