Quantitative comparison between simulations of seismic wave propagation in heterogeneous poro-elastic media and equivalent visco-elastic solids for marine-type environments

There is increasing evidence to suggest that the presence of mesoscopic heterogeneities constitutes an important seismic attenuation mechanism in porous rocks. As a consequence, centimetre-scale perturbations of the rock physical properties should be taken into account for seismic modelling whenever detailed and accurate responses of specific target structures are desired, which is, however, computationally prohibitive. A convenient way to circumvent this problem is to use an upscaling procedure to replace each of the heterogeneous porous media composing the geological model by corresponding equivalent visco-elastic solids and to solve the visco-elastic equations of motion for the inferred equivalent model. While the overall qualitative validity of this procedure is well established, there are as of yet no quantitative analyses regarding the equivalence of the seismograms resulting from the original poro-elastic and the corresponding upscaled visco-elastic models. To address this issue, we compare poro-elastic and visco-elastic solutions for a range of marine-type models of increasing complexity. We found that despite the identical dispersion and attenuation behaviour of the heterogeneous poro-elastic and the equivalent visco-elastic media, the seismograms may differ substantially due to diverging boundary conditions, where there exist additional options for the poro-elastic case. In particular, we observe that at the fluid/porous-solid interface, the poro- and visco-elastic seismograms agree for closed-pore boundary conditions, but differ significantly for open-pore boundary conditions. This is an important result which has potentially far-reaching implications for wave-equation-based algorithms in exploration geophysics involving fluid/porous-solid interfaces, such as, for example, wavefield decomposition.

P-wave velocities of samples from the Penninic of the western Alps along NFP20-west seismic-reflection profiles

Using the transit pulse method, we have determined compressional wave velocities of rocks from various geological units belonging to the Penninic zone along the NFP20-West profiles of the Swiss western Alps. The velocities have been measured at confining pressures up to 400 MPa, along three orthogonal axes defined by the macrostructure of the rocks. The samples analysed show a degree of metamorphism ranging from greenschist to eclogite facies. This collection includes schists, dolomites, gneisses and ophiolitic rocks. The mean velocities range from 5.9 km/s for a quartzitic calcschist to 7.9 km/s for an eclogitic metagabbro. The velocity anisotropy is as high as 20 %. The range of acoustic impedance is wide, from 15 to 27 10(6) kg/m2s. From these measurements, normal incident reflection coefficients for likely rock assemblages within and between geological units were estimated in order to interpret zone of the strong reflections recorded along the seismic profiles. Reflection coefficients as high as 0.17 could be determined.

Optical evidence for the proximity to a spin-density-wave metallic state in Na0.7CoO2

We present the optical properties of Na0.7CoO2 single crystals, measured over a broad spectral range as a function of temperature (T). The capability to cover the energy range from the far-infrared up to the ultraviolet allows us to perform reliable Kramers-Kronig transformation, in order to obtain the absorption spectrum (i.e., the complex optical conductivity). To the complex optical conductivity we apply the generalized Drude model, extracting the frequency dependence of the scattering rate (Gamma) and effective mass (m*) of the itinerant charge carriers. We find that Gamma(omega) at low temperatures and for similar to omega. This suggests that Na0.7CoO2 is at the verge of a spin-density-wave metallic phase.

Power spectral changes of the superimposed M wave during isometric voluntary contractions of increasing strength.

INTRODUCTION: We examined the power spectral changes of the compound muscle action potential (M wave) evoked during isometric contractions of increasing strength. METHODS: Surface electromyography (sEMG) of the vastus lateralis and medialis was recorded from 20 volunteers who performed 4-s step-wise isometric contractions of different intensities. A maximal M wave was elicited by a single stimulus to the femoral nerve and superimposed on the voluntary contractions. The spectral characteristics (Fmean and Fmedian) of sEMG and M-wave signals were calculated. RESULTS: M-wave spectral indicators increased systematically with contraction intensity up to 60% MVC and then leveled off at higher forces. Over the 10-60% MVC range, the increase in spectral indicators was 3 times higher for M waves (36%) than for sEMG (12%). CONCLUSIONS: The consistent increase in M-wave spectral characteristics with force is due to the fact that the number of motor units recruited by the superimposed supramaximal stimulus is approximately stable.

Two-dimensional probabilistic inversion of plane-wave electromagnetic data: Methodology, model constraints and joint inversion with electrical resistivity data

Probabilistic inversion methods based on Markov chain Monte Carlo (MCMC) simulation are well suited to quantify parameter and model uncertainty of nonlinear inverse problems. Yet, application of such methods to CPU-intensive forward models can be a daunting task, particularly if the parameter space is high dimensional. Here, we present a 2-D pixel-based MCMC inversion of plane-wave electromagnetic (EM) data. Using synthetic data, we investigate how model parameter uncertainty depends on model structure constraints using different norms of the likelihood function and the model constraints, and study the added benefits of joint inversion of EM and electrical resistivity tomography (ERT) data. Our results demonstrate that model structure constraints are necessary to stabilize the MCMC inversion results of a highly discretized model. These constraints decrease model parameter uncertainty and facilitate model interpretation. A drawback is that these constraints may lead to posterior distributions that do not fully include the true underlying model, because some of its features exhibit a low sensitivity to the EM data, and hence are difficult to resolve. This problem can be partly mitigated if the plane-wave EM data is augmented with ERT observations. The hierarchical Bayesian inverse formulation introduced and used herein is able to successfully recover the probabilistic properties of the measurement data errors and a model regularization weight. Application of the proposed inversion methodology to field data from an aquifer demonstrates that the posterior mean model realization is very similar to that derived from a deterministic inversion with similar model constraints.