Velocity variations in carbonate rocks due to dual porosity and wave-induced fluid flow

Stiffness variations in carbonates may be described as resulting from different concentrations of flat compliant pores or cracks, which can have a significant effect on the effective stiffness and acoustic properties (e.g., velocities and attenuations) of dry as well as saturated carbonates, although they carry extremely little porosity. As shown in this paper, the effects of dual porosity and wave-induced fluid flow or pore pressure communication may also play a significant role. On the basis of a previously published T-matrix approach to model the effective viscoelastic properties of cracked porous media, we illustrate the (frequency-dependent) effects of wave-induced fluid flow (mainly squirt flow) or pore pressure communication for a model structure consisting of a mixture of fluid-saturated porous grains and fluid-saturated cavities (vugs, etc.) that are embedded in a solid matrix associated with carbonates. We assume that the pores within the porous grains are decoupled from the pores in the solid matrix (and possibly saturated with different fluids) but that each pore system at the micro and/or mesoscale may or may not be connected. For each of four different connectivity models, we present numerical results for four different cases of microstructure (that emphasize the importance of cracks and flat compliant pores). Our numerical results indicate that the velocity and attenuation spectra of carbonates vary significantly, even when the crack density and all other volume concentrations are constant.     

Palaeo-hydrogeological control on groundwater As levels in Red River delta,Vietnam

To study the geological control on groundwater As concentrations in Red River delta, depth-specific groundwater sampling and geophysical logging in 11 monitoring wells was conducted along a 45 km transect across the southern and central part of the delta, and the literature on the Red River delta’s Quaternary geological development was reviewed. The water samples (n = 30) were analyzed for As, major ions, Fe²⁺, H₂S, NH₄, CH₄, δ¹⁸O and δD, and the geophysical log suite included natural gamma-ray, formation and fluid electrical conductivity. The SW part of the transect intersects deposits of grey estuarine clays and deltaic sands in a 15–20 km wide and 50–60 m deep Holocene incised valley. The NE part of the transect consists of 60–120 m of Pleistocene yellowish alluvial deposits underneath 10–30 m of estuarine clay overlain by a 10–20 m veneer of Holocene sediments. The distribution of δ¹⁸O-values (range −12.2‰ to −6.3‰) and hydraulic head in the sample wells indicate that the estuarine clay units divide the flow system into an upper Holocene aquifer and a lower Pleistocene aquifer. The groundwater samples were all anoxic, and contained Fe²⁺ (0.03–2.0 mM), Mn (0.7–320 μM), SO₄ (<2.1 μM–0.75 mM), H₂S (<0.1–7.0 μM), NH₄ (0.03–4.4 mM), and CH₄ (0.08–14.5 mM). Generally, higher concentrations of NH₄ and CH₄ and low concentrations of SO₄ were found in the SW part of the transect, dominated by Holocene deposits, while the opposite was the case for the NE part of the transect. The distribution of the groundwater As concentration (<0.013–11.7 μM; median 0.12 μM (9 μg/L)) is related to the distribution of NH₄, CH₄ and SO₄. Low concentrations of As (?0.32 μM) were found in the Pleistocene aquifer, while the highest As concentrations were found in the Holocene aquifer. PHREEQC-2 speciation calculations indicated that Fe²⁺ and H₂S concentrations are controlled by equilibrium for disordered mackinawite and precipitation of siderite. An elevated groundwater salinity (Cl range 0.19–65.1 mM) was observed in both aquifers, and dominated in the deep aquifer. A negative correlation between aqueous As and an estimate of reduced SO₄ was observed, indicating that Fe sulphide precipitation poses a secondary control on the groundwater As concentration.     

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.     

Transport Through the Contraction Area in the Little Belt

The Great Belt, the Øresund and the Little Belt connect the central Baltic Sea and the Kattegat. A fixed station was moored in the contraction area in the Little Belt during the period 18–28 July 1995, measuring temperature, salinity and current in two levels, while discharge was measured by the RVDana. The composite Froude number calculated at the fixed station shows that the two layer flow through this area was most often supercritical. The discharges were satisfactorily related to the currents measured at the fixed station, and time-series of transports through the Little Belt were established. When compared to the transports through the Øresund the water transport ratio (Øresund:Little Belt) was found to be 4·4, while the salt transport ratio was found to be 3·0. The resistance of the Little Belt, when considering the differences in sea level from Gedser to Hornbæk, was 1839×10⁻¹² s² m⁻⁵. On the basis of water level and surface salinity measurements made during the period 1931–76, a net discharge of 2300 m³ s⁻¹and a net salt transport of 36 tonnes s⁻¹through the Little Belt from the central Baltic Sea were found.     

Estimation of winter leaf area index and sky view fraction for snow modelling in boreal coniferous forests: consequences on snow mass and energy balance

Leaf area index (LAI) and canopy coverage are important parameters when modelling snow process in coniferous forests, controlling interception and transmitting radiation. Estimates of LAI and sky view factor show large variability depending on the estimation method used, and it is not clear how this is reflected in the calculated snow processes beneath the canopy. In this study, the winter LAI and sky view fraction were estimated using different optical and biomass‐based approximations in several boreal coniferous forest stands in Fennoscandia with different stand density, age and site latitude. The biomass‐based estimate of LAI derived from forest inventory data was close to the values derived from the optical measurements at most sites, suggesting that forest inventory data can be used as input to snow hydrological modelling. Heterogeneity of tree species and site fertility, as well as edge effects between different forest compartments, caused differences in the LAI estimates at some sites. A snow energy and mass balance model (SNOWPACK) was applied to detect how the differences in the estimated values of the winter LAI and sky view fraction were reflected in simulated snow processes. In the simulations, an increase in LAI and a decrease in sky view fraction changed the snow surface energy balance by decreasing shortwave radiation input and increasing longwave radiation input. Changes in LAI and sky view fraction affected directly snow accumulation through altered throughfall fraction and indirectly snowmelt through the changed surface energy balance. Changes in LAI and sky view fraction had a greater impact on mean incoming radiation beneath the canopy than on other energy fluxes. Snowmelt was affected more than snow accumulation. The effect of canopy parameters on evaporation loss from intercepted snow was comparable with the effect of variation in governing meteorological variables such as precipitation intensity and air temperature. Copyright © 2012 John Wiley & Sons, Ltd.     

Homotopy scattering series for seismic forward modelling with variable density and velocity

We have derived a convergent scattering series solution for the frequency-domain wave equation in acoustic media with variable density and velocity. The convergent scattering series solution is based on the homotopy analysis of a vectorial integral equation of the Lippmann–Schwinger type. By using the Green's function and partial integration, we have derived the vectorial integral equation of the Lippmann–Schwinger type that involves the pressure gradient field as well as the pressure field from the wave equation. The vectorial Lippmann–Schwinger equation can in principle be solved via matrix inversion, but the computational cost of matrix inversion scales like $N^3$, where $N$ is the number of grid blocks. The computational cost can be significantly reduced if one solves the vectorial Lippmann–Schwinger equation iteratively. A simple iterative solution is the Born series, but it is only convergent when the scattering potential is sufficiently small. In this study, we have used the so-called homotopy analysis method to derive an iterative solution for the vectorial Lippmann–Schwinger equation which can be made convergent even in strongly scattering media. The computational cost of our convergent scattering series scales as $N^2$. Our algorithm, which is based on the homotopy analysis method, involves a convergence control operator that we select using hierarchical matrices. We use a three-layer model and a resampled version of the SEG/EAGE salt model to show the performance of the developed convergent scattering series.     

Immiscible silicate liquid partition coefficients: implications for crystal-melt element partitioning and basalt petrogenesis

This study investigates partitioning of elements between immiscible aluminosilicate and borosilicate liquids using three synthetic mixtures doped with 32 trace elements. In order to get a good spatial separation of immiscible liquids, we employed a high-temperature centrifuge. Experiments were performed at 1,050–1,150°C, 1 atm, in sealed Fe and Pt containers. Quenched products were analysed by electron microprobe and LA ICP-MS. Nernst partition coefficients (D’s) between the Fe-rich and Si-rich aluminosilicate immiscible liquids are the highest for Zn (3.3) and Fe (2.6) and the lowest for Rb and K (0.4–0.5). The plots of D values against ionic potential Z/r in all the compositions show a convex upward trend, which is typical also for element partitioning between immiscible silicate and salt melts. The results bear upon the speciation and structural position of elements in multicomponent silicate liquids. The ferrobasalt–rhyolite liquid immiscibility is observed in evolved basaltic magmas, and may play an important role in large gabbroic intrusions, such as Skaergaard, and during the generation of unusual lavas, such as ferropicrites.