Poroelastic Backus averaging for anisotropic layered fluid‐ and gas‐saturated sediments

Geophysics ◽  
1997 ◽  
Vol 62 (6) ◽  
pp. 1867-1878 ◽  
Author(s):  
Stephan Gelinsky ◽  
Sergei A. Shapiro

A homogeneous anisotropic effective‐medium model for saturated thinly layered sediments is introduced. It is obtained by averaging over many layers with different poroelastic moduli and different saturating fluids. For a medium consisting of a stack of vertically fractured horizontal layers, this effective medium is orthorhombic. We derive the poroelastic constants that define such media in the long‐wavelength limit as well as the effective large‐scale permeability tensor. The permeability shows strong anisotropy for large porosity fluctuations. We observe pronounced effects that do not exist in purely elastic media. At very low frequencies, seismic waves cause interlayer flow of pore fluid across interfaces from more compliant into stiffer layers. For higher frequencies, the layers behave as if they are sealed, and no fluid flow occurs. The effective‐medium velocities of the quasi‐compressional waves are higher in the no‐flow than in the quasi‐static limit. Both are lower than the high‐frequency, i.e., ray‐theory limit. Partial saturation affects the anisotropy of wave propagation. In the no‐flow limit, gas that is accumulated primarily in the stiffer layers reduces the seismic anisotropy; gas that is trapped mainly in layers with a more compliant frame tends to increase the anisotropy. In the quasi‐static limit, local flow keeps the anisotropy constant independent of partial saturation effects. For dry rock, no‐flow and quasi‐static velocities are the same, and the anisotropy caused by layering is controlled only by fluctuations of the layer shear moduli. If the shear stiffness of all layers is the same and only the compressive stiffness or saturation varies, only the ray‐theory velocity exhibits anisotropy.

2021 ◽  
Vol 4 (1) ◽  
Author(s):  
Soren Wainio-Theberge ◽  
Annemarie Wolff ◽  
Georg Northoff

AbstractSpontaneous neural activity fluctuations have been shown to influence trial-by-trial variation in perceptual, cognitive, and behavioral outcomes. However, the complex electrophysiological mechanisms by which these fluctuations shape stimulus-evoked neural activity remain largely to be explored. Employing a large-scale magnetoencephalographic dataset and an electroencephalographic replication dataset, we investigate the relationship between spontaneous and evoked neural activity across a range of electrophysiological variables. We observe that for high-frequency activity, high pre-stimulus amplitudes lead to greater evoked desynchronization, while for low frequencies, high pre-stimulus amplitudes induce larger degrees of event-related synchronization. We further decompose electrophysiological power into oscillatory and scale-free components, demonstrating different patterns of spontaneous-evoked correlation for each component. Finally, we find correlations between spontaneous and evoked time-domain electrophysiological signals. Overall, we demonstrate that the dynamics of multiple electrophysiological variables exhibit distinct relationships between their spontaneous and evoked activity, a result which carries implications for experimental design and analysis in non-invasive electrophysiology.


2021 ◽  
Vol 2 (1) ◽  
Author(s):  
Cédric P. Legendre ◽  
Li Zhao ◽  
Tai-Lin Tseng

AbstractThe average anisotropy beneath Anatolia is very strong and is well constrained by shear-wave splitting measurements. However, the vertical layering of anisotropy and the contribution of each layer to the overall pattern is still an open question. Here, we construct anisotropic phase-velocity maps of fundamental-mode Rayleigh waves for the Anatolia region using ambient noise seismology and records from several regional seismic stations. We find that the anisotropy patterns in the crust, lithosphere and asthenosphere beneath Anatolia have limited amplitudes and are generally consistent with regional tectonics and mantle processes dominated by the collision between Eurasia and Arabia and the Aegean/Anatolian subduction system. The anisotropy of these layers in the crust and upper mantle are, however, not consistent with the strong average anisotropy measured in this area. We therefore suggest that the main contribution to overall anisotropy likely originates from a deep and highly anisotropic region round the mantle transition zone.


Geophysics ◽  
1994 ◽  
Vol 59 (1) ◽  
pp. 87-92 ◽  
Author(s):  
Gary Mavko ◽  
Diane Jizba

Seismic velocity dispersionin fluid-saturated rocks appears to be dominated by tow mecahnisms: the large scale mechanism modeled by Biot, and the local flow or squirt mecahnism. The tow mechanisms can be distuinguished by the ratio of P-to S-wave dispersions, or more conbeniently, by the ratio of dynamic bulk to shear compliance dispersions derived from the wave velocities. Our formulation suggests that when local flow denominates, the dispersion of the shear compliance will be approximately 4/15 the dispersion of the compressibility. When the Biot mechanism dominates, the constant of proportionality is much smaller. Our examination of ultrasonic velocities from 40 sandstones and granites shows that most, but not all, of the samples were dominated by local flow dispersion, particularly at effective pressures below 40 MPa.


2018 ◽  
Vol 61 ◽  
pp. 1-37 ◽  
Author(s):  
Paola F. Antonietti ◽  
Alberto Ferroni ◽  
Ilario Mazzieri ◽  
Roberto Paolucci ◽  
Alfio Quarteroni ◽  
...  

We present a comprehensive review of Discontinuous Galerkin Spectral Element (DGSE) methods on hybrid hexahedral/tetrahedral grids for the numerical modeling of the ground motion induced by large earthquakes. DGSE methods combine the exibility of discontinuous Galerkin meth-ods to patch together, through a domain decomposition paradigm, Spectral Element blocks where high-order polynomials are used for the space discretization. This approach allows local adaptivity on discretization parameters, thus improving the quality of the solution without affecting the compu-tational costs. The theoretical properties of the semidiscrete formulation are also revised, including well-posedness, stability and error estimates. A discussion on the dissipation, dispersion and stability properties of the fully-discrete (in space and time) formulation is also presented. Here space dis-cretization is obtained based on employing the leap-frog time marching scheme. The capabilities of the present approach are demonstrated through a set of computations of realistic earthquake scenar-ios obtained using the code SPEED (http://speed.mox.polimi.it), an open-source code specifically designed for the numerical modeling of large-scale seismic events jointly developed at Politecnico di Milano by The Laboratory for Modeling and Scientific Computing MOX and by the Department of Civil and Environmental Engineering.


1997 ◽  
Vol 40 (3) ◽  
Author(s):  
L. Margheriti ◽  
C. Nostro ◽  
A. Amato ◽  
M. Cocco

Anisotropy is a common property of the Earth's crust and the upper mantle; it is related to the strain field of the medium and therefore to geodynamics. In this paper we describe the different possible origins of anisotropic behavior of the seismic waves and the seismological techniques used to define anisotropic bodies. In general it is found that the fast polarization direction is parallel to the absolute plate motion in cratonic areas, to the spreading direction near rifts or extensional zones, and to the main structural features in transpressive regimes. The delay times between fast and slow waves reflect the relative strength and penetration at depth of the deformation field. The correspondence between surface structural trends and anisotropy in the upper mantle, found in many regions of the world, strongly suggest that orogenic processes involve not only the shallow crust but the entire lithosphere. Recently in Italy both shear wave splitting analysis and Pn inversion were applied to define the trend of seismic anisotropy. Along the Northern Appeninic arc fast directions follow the strike of the arc (i.e., parallel to the strike of the Miocene-Pleistocene compressional features), whereas in the Tyrrhenian zone fast directions are about E-W SW-NE; parallel to the post-Miocene extension that is thought to have reoriented the mantle minerals fabric in the astenosphere.


2020 ◽  
Author(s):  
Finnigan Illsley-Kemp ◽  
Martha Savage ◽  
Colin Wilson ◽  
S Bannister

© 2019. American Geophysical Union. All Rights Reserved. We use crustal seismic anisotropy measurements in the North Island, New Zealand, to examine structures and stress within the Taupō Volcanic Zone, the Taranaki Volcanic Lineament, the subducting Hikurangi slab, and the Hikurangi forearc. Results in the Taranaki region are consistent with NW-SE oriented extension yet suggest that the Taranaki volcanic lineament may be controlled by a deep-rooted, inherited crustal structure. In the central Taupō Volcanic Zone anisotropy fast orientations are predominantly controlled by continental rifting. However at Taupō and Okataina volcanoes, fast orientations are highly variable and radial to the calderas suggesting the influence of magma reservoirs in the seismogenic crust (≤15 km depth). The subducting Hikurangi slab has a predominant trench-parallel fast orientation, reflecting the pervasive presence of plate-bending faults, yet changing orientations at depths ≥120 km beneath the central North Island may be relics from previous subduction configurations. Finally, results from the southern Hikurangi forearc show that the orientation of stresses there is consistent with those in the underlying subducting slab. In contrast, the northern Hikurangi forearc is pervasively fractured and is undergoing E-W compression, oblique to the stress field in the subducting slab. The north-south variation in fore-arc stress is likely related to differing subduction-interface coupling. Across the varying tectonic regimes of the North Island our study highlights that large-scale tectonic forces tend to dictate the orientation of stress and structures within the crust, although more localized features (plate coupling, magma reservoirs, and inherited crustal structures) can strongly influence surface magmatism and the crustal stress field.


2011 ◽  
Vol 94-96 ◽  
pp. 1941-1945
Author(s):  
Yi Wu ◽  
Chun Yang ◽  
Jian Cai ◽  
Jian Ming Pan

Elasto-plastic analysis of seismic responses of valve hall structures were carried out by using finite element software, and the effect of seismic waves on the seismic responses of the valve hall structures and suspension equipments were studied. Results show that significant torsional responses of the structure can be found under longitudinal and 3D earthquake actions. Under 3D earthquake actions, the seismic responses of the suspension valves are much more significant than those under 1D earthquake actions, the maximum tensile force of the suspenders is about twice of that under 1D action. The seismic responses of the suspension valves under vertical earthquake actions are much stronger than those under horizontal earthquake actions, when suffering strong earthquake actions; the maximum vertical acceleration of the suspension valves is about 4 times of that under horizontal earthquake actions. It is recommended that the effects of 3D earthquake actions on the structure should be considered in seismic response analysis of the valve hall structure.


2021 ◽  
Author(s):  
John Keith Magali ◽  
Thomas Bodin ◽  
Navid Hedjazian ◽  
Yanick Ricard ◽  
Yann Capdeville

<p>Large-scale seismic anisotropy inferred from seismic observations has been loosely interpreted either in terms of intrinsic anisotropy due to Crystallographic Preferred Orientation (CPO) development of mantle minerals or extrinsic anisotropy due to rock-scale Shape Preferred Orientation (SPO). The coexistence of both contributions misconstrues the origins of seismic anisotropy observed in seismic tomography models. It is thus essential to discriminate CPO from SPO in the effective anisotropy of an upscaled/homogenized medium, that is, the best possible elastic model recovered using finite-frequency seismic data assuming perfect data coverage. In this work, we investigate the effects of upscaling an intrinsically-anisotropic and highly-heterogeneous Earth's mantle. The problem is applied to a 2-D marble cake model of the mantle with a binary composition in the presence of CPO obtained from a micro-mechanical model. We compute the long-wavelength effective equivalent of this mantle model using the 3D non-periodic elastic homogenization technique. Our numerical findings predict that overall, upscaling purely intrinsically anisotropic medium amounts to the convection-scale averaging of CPO. As a result, it always underestimates the anisotropy, and may only be overestimated due to the additive extrinsic anisotropy from SPO. Finally, we show analytically (in 1D) and numerically (in 2D) that the full effective radial anisotropy ξ<sup>*</sup> is approximately just the product of the effective intrinsic radial anisotropy ξ<sup>*</sup><sub>CPO</sub> and the extrinsic radial anisotropy ξ<sub>SPO</sub>:</p><p>ξ<sup>* </sup>= ξ<sup>*</sup><sub>CPO </sub>× ξ<sub>SPO</sub></p><p>Based on the above relation, it is imperative to homogenize a texture evolution model first before drawing interpretations from existing anisotropic tomography models. Such a scaling law can therefore be used as a constraint to better estimate the separate contributions of CPO and SPO from the effective anisotropy observed in tomographic models.</p>


Geophysics ◽  
2019 ◽  
Vol 85 (1) ◽  
pp. MR1-MR10 ◽  
Author(s):  
Fuyong Yan ◽  
De-Hua Han ◽  
Tongcheng Han ◽  
Xue-Lian Chen

The layer-induced seismic anisotropy of sedimentary strata is frequency-dependent. At the low-frequency limit, the effective anisotropic properties of the layered media can be estimated by the Backus averaging model. At the high-frequency limit, the apparent anisotropic properties of the layered media can be estimated by ray theory. First, we build a database of laboratory ultrasonic measurement on sedimentary rocks from the literature. The database includes ultrasonic velocity measurements on sandstones and carbonate rocks, and velocity-anisotropy measurements on shales. Then, we simulate the sedimentary strata by randomly selecting a certain number of rock samples and using their laboratory measurement results to parameterize each layer. For each realization of the sedimentary strata, we estimate the effective and apparent seismic anisotropy parameters using the Backus average and ray theory, respectively. We find that, relative to Backus averaging, ray theory usually underestimates the Thomsen parameters [Formula: see text] and [Formula: see text], and overestimates [Formula: see text]. For an effective layered medium consisting of isotropic sedimentary rocks, the differences are significant. These differences decrease when shales with intrinsic seismic anisotropy are included. For the same sedimentary strata, the seismic wave should perceive stronger seismic anisotropy than the ultrasonic wave.


Sign in / Sign up

Export Citation Format

Share Document