Seismic
anisotropy
Anisotropy () is the property of a material which allows it to change or assume different properties in different directions, as opposed to isotropy. It can be defined as a difference, when measured along different axes, in a material's physic ...
is a term used in
seismology
Seismology (; from Ancient Greek σεισμός (''seismós'') meaning "earthquake" and -λογία (''-logía'') meaning "study of") is the scientific study of earthquakes and the propagation of elastic waves through the Earth or through other ...
to describe the directional dependence of the velocity of
seismic wave
A seismic wave is a wave of acoustic energy that travels through the Earth. It can result from an earthquake, volcanic eruption, magma movement, a large landslide, and a large man-made explosion that produces low-frequency acoustic energy ...
s in a medium (
rock
Rock most often refers to:
* Rock (geology), a naturally occurring solid aggregate of minerals or mineraloids
* Rock music, a genre of popular music
Rock or Rocks may also refer to:
Places United Kingdom
* Rock, Caerphilly, a location in Wales ...
) within the
Earth
Earth is the third planet from the Sun and the only astronomical object known to harbor life. While large volumes of water can be found throughout the Solar System, only Earth sustains liquid surface water. About 71% of Earth's surfa ...
.
Description
A
material
Material is a substance or mixture of substances that constitutes an object. Materials can be pure or impure, living or non-living matter. Materials can be classified on the basis of their physical and chemical properties, or on their geolo ...
is said to be anisotropic if the value of one or more of its properties varies with direction. Anisotropy differs from the property called heterogeneity in that anisotropy is the variation in values with direction at a point while heterogeneity is the variation in values between two or more points.
Seismic Anisotropy can be defined as the dependence of seismic velocity on direction or upon angle. General anisotropy is described by a 4th order elasticity tensor with 21 independent elements. However, in practice observational studies are unable to distinguish all 21 elements, and anisotropy is usually simplified. In the simplest form, there are two main types of anisotropy, both of them are called
transverse isotropy (it is called transverse isotropy because there is isotropy in either the horizontal or vertical plane) or polar anisotropy. The difference between them is in their axis of symmetry, which is an axis of rotational invariance such that if we rotate the formation about the axis, the material is still indistinguishable from what it was before. The symmetry axis is usually associated with regional stress or gravity.
* TIV- transverse isotropy with a vertical axis of symmetry, this is also called VTI (vertical transverse isotropy). This kind of anisotropy is associated with
layering
Layering has evolved as a common means of vegetative propagation of numerous species in natural environments. Layering is also utilized by horticulturists to propagate desirable plants.
Natural layering typically occurs when a branch touches ...
and
shale and is found where
gravity
In physics, gravity () is a fundamental interaction which causes mutual attraction between all things with mass or energy. Gravity is, by far, the weakest of the four fundamental interactions, approximately 1038 times weaker than the stro ...
is the dominant factor.
* TIH- transverse isotropy with a horizontal axis of symmetry, this is also called HTI (horizontal transverse isotropy). This kind of anisotropy is associated with cracks and fractures and is found where regional stress is the dominant factor.
The transverse anisotropic matrix has the same form as the isotropic matrix, except that it has five non-zero values distributed among 12 non-zero elements.
Transverse isotropy is sometimes called transverse anisotropy or anisotropy with hexagonal symmetry. In many cases the axis of symmetry will be neither horizontal nor vertical, in which case it is often called "tilted".
History of the recognition of anisotropy
Anisotropy dates back to the 19th century following the theory of Elastic wave propagation.
Green
Green is the color between cyan and yellow on the visible spectrum. It is evoked by light which has a dominant wavelength of roughly 495570 nm. In subtractive color systems, used in painting and color printing, it is created by a combi ...
(1838) and
Lord Kelvin
William Thomson, 1st Baron Kelvin, (26 June 182417 December 1907) was a British mathematician, mathematical physicist and engineer born in Belfast. Professor of Natural Philosophy at the University of Glasgow for 53 years, he did important ...
(1856) took anisotropy into account in their articles on wave propagation. Anisotropy entered
seismology
Seismology (; from Ancient Greek σεισμός (''seismós'') meaning "earthquake" and -λογία (''-logía'') meaning "study of") is the scientific study of earthquakes and the propagation of elastic waves through the Earth or through other ...
in the late 19th century and was introduced by
Maurice Rudzki. From 1898 till his death in 1916,
Rudzki attempted to advance the theory of anisotropy, he attempted to determine the wavefront of a transversely isotropic medium (TI) in 1898 and in 1912 and 1913 he wrote on surface waves in transversely isotropic half space and on Fermat's principle in anisotropic media respectively.
With all these, the advancement of anisotropy was still slow and in the first 30 years (1920-1950) of exploration seismology only a few papers were written on the subject. More work was done by several scientists such as
Helbig (1956) who observed while doing seismic work on
Devonian schist
Schist ( ) is a medium-grained metamorphic rock showing pronounced schistosity. This means that the rock is composed of mineral grains easily seen with a low-power hand lens, oriented in such a way that the rock is easily split into thin flakes ...
s that velocities along the foliation were about 20% higher than those across the foliation. However the appreciation of anisotropy increased with the proposition of a new model for the generation of anisotropy in an originally isotropic background and a new exploration concept by Crampin (1987). One of the main points by Crampin was that the polarization of three component shear waves carries unique information about the internal structure of the rock through which they pass, and that
shear wave splitting
Shear wave splitting, also called seismic birefringence, is the phenomenon that occurs when a Polarization (waves), polarized shear wave enters an anisotropic medium (Fig. 1). The incident shear wave splits into two polarized shear waves (Fig. 2). ...
may contain information about the distribution of
crack orientations.
With these new developments and the acquisition of better and new types of data such as three component 3D
seismic data, which clearly show the effects of shear wave splitting, and wide Azimuth 3D data which show the effects of Azimuthal anisotropy, and the availability of more powerful computers, anisotropy began to have great impact in
exploration
Exploration refers to the historical practice of discovering remote lands. It is studied by geographers and historians.
Two major eras of exploration occurred in human history: one of convergence, and one of divergence. The first, covering most ...
seismology in the past three decades.
Concept of seismic anisotropy
Since the understanding of seismic anisotropy is closely tied to the shear wave splitting, this section begins with a discussion of shear wave splitting.
Shear waves have been observed to split into two or more fixed polarizations which can propagate in the particular ray direction when entering an anisotropic medium. These split phases propagate with different polarizations and velocities. Crampin (1984) amongst others gives evidence that many rocks are anisotropic for shear wave propagation. In addition, shear wave splitting is almost routinely observed in three-component
VSPs. Such shear wave splitting can be directly analyzed only on three component geophones recording either in the subsurface, or within the effective shear window at the free surface if there are no near surface low-velocity layers. Observation of these shear waves show that measuring the orientation and polarization of the first arrival and the delay between these split shear waves reveal the orientation of cracks and the crack density . This is particularly important in reservoir characterization.
In a linearly elastic material, which can be described by
Hooke's law
In physics, Hooke's law is an empirical law which states that the force () needed to extend or compress a spring by some distance () scales linearly with respect to that distance—that is, where is a constant factor characteristic of ...
as one in which each component of stress is dependent on every component of strain, the following relationship exists:
:
where ''σ'' is the stress, ''C'' is the
elastic moduli
An elastic modulus (also known as modulus of elasticity) is the unit of measurement of an object's or substance's resistance to being deformed elastically (i.e., non-permanently) when a stress is applied to it. The elastic modulus of an object is ...
or stiffness constant, and ''e'' is the strain.
The elastic modulus matrix for an anisotropic case is
:
The above is the elastic modulus for a vertical transverse isotropic medium (VTI), which is the usual case. The elastic modulus for a horizontal transverse isotropic medium (HTI) is;
:
For an anisotropic medium, the directional dependence of the three phase velocities can be written by applying the elastic moduli in the wave equation is; The direction dependent
wave speeds for
elastic wave
Linear elasticity is a mathematical model of how solid objects deform and become internally stressed due to prescribed loading conditions. It is a simplification of the more general nonlinear theory of elasticity and a branch of continuum mec ...
s through the material can be found by using the
Christoffel equation and are given by
:
where
is the angle between the axis of symmetry and the wave propagation direction,
is mass density and the
are elements of the
elastic stiffness matrix. The Thomsen parameters are used to simplify these expressions and make them easier to understand.
Seismic anisotropy has been observed to be weak, and Thomsen (1986)
rewrote the velocities above in terms of their deviation from the vertical velocities as follows;
:
where
:
are the P and S wave velocities in the direction of the axis of symmetry (
) (in geophysics, this is usually, but not always, the vertical direction). Note that
may be further linearized, but this does not lead to further simplification.
The approximate expressions for the wave velocities are simple enough to be physically interpreted, and sufficiently accurate for most geophysical applications. These expressions are also useful in some contexts where the anisotropy is not weak.
The Thomsen parameters are anisotropic and are three non-dimensional combinations which reduce to zero in isotropic cases, and are defined as
:
Origin of anisotropy
Anisotropy has been reported to occur in the Earth's three main layers; the
crust,
mantle and the
core
Core or cores may refer to:
Science and technology
* Core (anatomy), everything except the appendages
* Core (manufacturing), used in casting and molding
* Core (optical fiber), the signal-carrying portion of an optical fiber
* Core, the centra ...
.
The origin of seismic anisotropy is non-unique, a range of phenomena may cause Earth materials to display seismic anisotropy. The anisotropy may be strongly dependent on wavelength if it is due to the average properties of aligned or partially aligned heterogeneity. A solid has intrinsic anisotropy when it is homogeneously and sinuously anisotropic down to the smallest particle size, which may be due to crystalline anisotropy. Relevant crystallographic anisotropy can be found in the
upper mantle
The upper mantle of Earth is a very thick layer of rock inside the planet, which begins just beneath the crust (at about under the oceans and about under the continents) and ends at the top of the lower mantle at . Temperatures range from appr ...
. When an otherwise isotropic rock contains a distribution of dry or liquid-filled cracks which have preferred orientation it is named crack induced anisotropy. The presence of aligned cracks, open or filled with some different material, is an important mechanism at shallow depth, in the crust. It is well known that the small-scale, or microstructural, factors include (e.g. Kern & Wenk 1985; Mainprice et al. 2003): (1) crystal lattice preferred orientation (LPO) of constituent mineral phases; (2) variations in spatial distribution of grains and minerals; (3) grain morphology and (4) aligned fractures, cracks and pores, and the nature of their infilling material (e.g. clays, hydrocarbons, water, etc.). Because of the overall microstructural control on seismic anisotropy, it follows that anisotropy can be diagnostic for specific rock types. Here, we consider whether seismic anisotropy can be used as an indicator of specific
sedimentary
Sedimentary rocks are types of rock that are formed by the accumulation or deposition of mineral or organic particles at Earth's surface, followed by cementation. Sedimentation is the collective name for processes that cause these particles ...
lithologies within the Earth's crust.
In sedimentary rocks, anisotropy develops during and after deposition. For anisotropy to develop, there needs to be some degree of homogeneity or uniformity from point to point in the deposited clastics. During deposition, anisotropy is caused by the periodic layering associated with changes in sediment type which produces materials of different grain size, and also by the directionality of the transporting medium which tends to order the grains under gravity by grain sorting. Fracturing and some
diagenetic
Diagenesis () is the process that describes physical and chemical changes in sediments first caused by water-rock interactions, microbial activity, and compaction after their deposition. Increased pressure and temperature only start to play a ...
processes such as
compaction and dewatering of
clays
Clay is a type of fine-grained natural soil material containing clay minerals (hydrous aluminium phyllosilicates, e.g. kaolin, Al2 Si2 O5( OH)4).
Clays develop plasticity when wet, due to a molecular film of water surrounding the clay par ...
, and alteration etc. are post depositional processes that can cause anisotropy .
The importance of anisotropy in hydrocarbon exploration and production
In the past two decades, the seismic anisotropy has dramatically been gaining attention from academic and industry, due to advances in anisotropy parameter estimation, the transition from post stack imaging to pre stack depth migration, and the wider offset and azimuthal coverage of 3D surveys. Currently, many seismic processing and inversion methods utilize anisotropic models, thus providing a significant enhancement over the seismic imaging quality and resolution. The integration of anisotropy velocity model with seismic imaging has reduced uncertainty on internal and bounding-
fault positions, thus greatly reduce the risk of investment decision based heavily on seismic interpretation.
In addition, the establishment of correlation between anisotropy parameters, fracture orientation, and density, lead to practical reservoir characterization techniques. The acquisition of such information, fracture spatial distribution and density, the drainage area of each producing well can be dramatically increased if taking the fractures into account during the drilling decision process. The increased drainage area per well will result in fewer wells, greatly reducing the drilling cost of exploration and production (E&P) projects.
The application of the anisotropy in petroleum exploration and production
Among several applications of seismic anisotropy, the following are the most important: anisotropic parameter estimation, prestack depth anisotropy migration, and fracture characterization based on anisotropy velocity models.
Anisotropy parameter estimation
The anisotropy parameter is most fundamental to all other anisotropy application in E&P area. In the early days of seismic petroleum exploration, the geophysicists were already aware of the anisotropy-induced distortion in
P-wave imaging (the major of petroleum exploration seismic surveys). Although the anisotropy-induced distortion is less significant since the poststack processing of narrow-azimuth data is not sensitive to velocity. The advancement of seismic anisotropy is largely contributed by the Thomsen's work on anisotropy notation and also by the discovery of the P-wave time-process parameter
. These fundamental works enable to parametrize the transverse isotropic (TI) models with only three parameters, while there are five full independent stiff tensor element in transverse isotropic (VTI or HTI) models. This simplification made the measurement of seismic anisotropy a plausible approach.
Most anisotropy parameter estimation work is based on shale and
silt
Silt is granular material of a size between sand and clay and composed mostly of broken grains of quartz. Silt may occur as a soil (often mixed with sand or clay) or as sediment mixed in suspension with water. Silt usually has a floury feel ...
s, which may be due to the fact that shale and silts are the most abundant sedimentary rocks in the Earth's crust. Also in the context of
petroleum geology, organic shale is the
source rock
In petroleum geology, source rock is rock which has generated hydrocarbons or which could generate hydrocarbons. Source rocks are one of the necessary elements of a working petroleum system. They are organic-rich sediments that may have been depo ...
as well as seal rocks that trap oil and gas. In seismic exploration, shales represent the majority of the wave propagation medium overlying the petroleum
reservoir
A reservoir (; from French ''réservoir'' ) is an enlarged lake behind a dam. Such a dam may be either artificial, built to store fresh water or it may be a natural formation.
Reservoirs can be created in a number of ways, including contro ...
. In conclusion, seismic properties of shale are important for both exploration and reservoir management.
Seismic velocity anisotropy in shale can be estimated from several methods, including deviated-well sonic logs, walkway VSP, and core measurement. These methods have their own advantages and disadvantages: the walkway VSP method suffers from scaling issues, and core measure is impractical for shale, since shale is hard to be cored during drilling.
Walkway VSP
The Walkway VSP array several seismic surface sources at different offset from the well. Meanwhile, a vertical receiver array with constant interval between receivers is mounted in a vertical well. The sound arrival times between multiple surface sources and receivers at multiple depths are recorded during measurement. These arrival times are used to derive the anisotropy parameter based on the following equations
:
Where
is the arrival time from source with
offset,
is the arrival time of zero offset,
is NMO velocity,
is Thompson anisotropy parameter.
The layout of surface sources and receivers positions is shown in the following diagram.
Core measurement
Another technique used to estimate the anisotropy parameter is directly measure them from the core which is extracted through a special hollow drill bit during drill process. Since coring a sample will generate large extra cost, only limited number of core samples can be obtained for each well. Thus the anisotropy parameter obtained through core measurement technique only represent the anisotropy property of rock near the borehole at just several specific depth, rending this technique often provides little help on the field seismic survey application. The measurements on each shale plug require at least one week.
From the context of this article, wave propagation in a vertically transverse medium can be described with five elastic constants, and ratios among these parameters define the rock anisotropy . This anisotropy parameter can be obtained in the laboratory by measuring the velocity travel speed with transducer ultrasonic systems at variable saturation and pressure conditions. Usually, three directions of wave propagation on core samples are the minimum requirement to estimate the five elastic coefficients of the stiffness tensor. Each direction in core plug measurement yields three velocities (one P and two S).
The variation of wave propagation direction can be achieved by either cutting three samples at 0°, 45° and 90° from the cores or by using one core plug with transducers attached at these three angles. Since most shales are very friable and fissured, it is often difficult to cut shale core plug. Its edges break off easily. Thus the cutting sample method can only be used for hard, competent rocks. The cutting position of samples can be explained by the following diagram.
Another way to get the wave propagation velocity at three directions is to arrange the ultrasonic transducer onto several specific location of the core sampler. This method avoids the difficulties encounter during the cutting of shale core sample. It also reduces the time of measurement by two thirds since three pairs of ultrasonic transducer work at the same time. The following diagram gives us a clear image on the arrangement of the transducers.
Once the velocities at three directions are measured by one of the above two methods, the five independent elastic constants are given by the following equations:
:
The P-wave anisotropy of a VTI medium can be described by using Thomsen's parameters
. The
quantifies the velocity difference for wave propagation along and perpendicular to the symmetry axis, while
controls the P-wave propagation for angles near the symmetry axis.
Deviated well sonic log
The last technique can be used to measure the seismic anisotropy is related to the sonic logging information of a deviated well. In a deviated well, the wave propagation velocity is higher than the wave propagation velocity in a vertical well at the same depth. This difference in velocity between deviated well and vertical well reflects the anisotropy parameters of the rocks near the borehole. The detail of this technique will be shown on an example of this report.
Anisotropic prestack depth migration
In the situation of complex geology, e.g. faulting, folding, fracturing, salt bodies, and unconformities, pre-stack migration (PreSM) is used due to better resolution under such complex geology. In PreSM, all traces are migrated before being moved to zero-offset. As a result, much more information is used, which results in a much better image, along with the fact that PreSM honours velocity changes more accurately than post-stack migration. The PreSM is extremely sensitive to the accuracy of the velocity field. Thus the inadequacy of isotropic velocity models is not suitable for the pre stack depth migration. P-wave anisotropic prestack depth migration (APSDM) can produce a seismic image that is very accurate in depth and space. As a result, unlike isotropic PSDM, it is consistent with well data and provides an ideal input for reservoir characterization studies. However, this accuracy can only be achieved if correct anisotropy parameters are used. These parameters cannot be estimated from seismic data alone. They can only be determined with confidence through analysis of a variety of geoscientific material – borehole data and geological history.
During recent years, the industry has started to see the practical use of anisotropy in seismic imaging. We show case studies that illustrate this integration of the geosciences. We show that much better accuracy is being achieved. The logical conclusion is that, this integrated approach should extend the use of anisotropic depth imaging from complex geology only, to routine application on all reservoirs.
Fracture characterization
After considering applications of anisotropy that improved seismic imaging, two approaches for exploiting anisotropy for the analysis of fractures in the formation are worthy of discussing. Ones uses azimuthal variations in the amplitude versus offset (AVO) signature when the wave is reflected from the top or base of an anisotropic material, and a second exploits the polarizing effect that the fractures have on a transmitted shear-wave. In both cases, the individual fractures are below the resolving power of the seismic signal and it is the cumulative effect of the fracturing that is recorded. Based on the idea behind them, both approaches can be divided into two steps. The first step is to get the anisotropy parameters from seismic signals, and the second steps is to retreat the information of fractures from anisotropy parameters based on the fracture induce anisotropy model.
Fractures-azimuthal variations
Aligned subseismic-scale fracturing can produce seismic anisotropy (i.e., seismic velocity varies with direction) and leads to measurable directional differences in traveltimes and reflectivity.
If the fractures are vertically aligned, they will produce azimuthal anisotropy (the simplest case being horizontal transverse isotropy, or HTI) such that reflectivity of an interface depends on azimuth as well as offset. If either of the media bounding the interface is azimuthally anisotropic, the AVO will have an azimuthal dependence . The P-P wave reflection coefficient have the following relation with the azimuthal if anisotropy exist in the layers:
:
Where
is the azimuth from data acquisition grid, the terms
are coefficients describing anisotropy parameter.
Fractures- shear-wave splitting
The behavior of shear waves as they pass through anisotropic media has been recognized for many years, with laboratory and field observations demonstrating how the shear wave splits into two polarized components with their planes aligned parallel and perpendicular to the anisotropy. For a fractured medium, the faster shear wave is generally aligned with the strike direction and the time delay between the split shear waves related to the fracture density and path length traveled. For layered medium, the shear wave polarized parallel to the layering arrives first.
Examples of the application of anisotropy
Example of anisotropy in petroleum E&P
Two examples will be discussed in there to show the anisotropy application in Petroleum E&P area. The first related to anisotropy parameter estimation via deviated well sonic logging tool. And the second example reflects the image quality improvement by PreStack Depth Migration technology.
Example of deviated well sonic logging
In this case, the sonic velocity in a deviated well is obtained by dipole sonic logging tool . The formation is mostly composed of shale. In order to use the TI model, several assumptions are made:
* Rock should be in normally pressured regime.
* Rock should have similar burial history.
Satisfying the above conditions, the following equation hold for a TI model:
:
Where
is the deviated angle of the well, and
,
are anisotropy parameter.
The following plot shows typical velocity distribution vs density in a deviated well. The color of each data point represents the frequency of this data point. The red color means a high frequency while the blue color represents a low frequency. The black line shows a typical velocity trend without the effect of anisotropy. Since the existence of anisotropy effect, the sound velocity is higher than the trend line.
From the well logging data, the velocity vs
plot can be drawn. On the basis of this plot, a no liner regression will give us an estimate of
and
. The following plot show the non-linear regression and its result.
Put the estimated
and
into the following equation, the correct
can be obtained.
:
By doing the above correction calculation, the corrected
is plot vs density in the following plot. As be seen in the plot, most of the data point falls on the trend line. It validate the correctness of the estimate of anisotropy parameter.
Example of prestack depth migration Imaging
In this case, the operator conducted several seismic surveys on a gas field in the north sea over the period of 1993-1998 . The early survey does not take anisotropy into account, while the later survey employs the PreStack Depth Migration imaging. This PSDM was done on a commercial seismic package developed by Total. The following two plots clearly reveal the resolution improvement of the PSDM method. The top plot is a convention 3D survey without anisotropy effect. The bottom plot used PSDM method. As can be seen in the bottom plot, more small structure features are revealed due to the reduce of error and improved resolution.
Limitations of seismic anisotropy
Seismic anisotropy relies on shear waves, shear waves carry rich information which can sometimes impede its utilization. Shear waves survey for anisotropy requires multi component (usually 3 component) geophones which are oriented at angles, these are more expensive than the widely used vertical oriented single component geophones. However, while expensive 3 component seismometers are much more powerful in their ability to collect valuable information about the Earth that vertical component seismometers simply cannot. While seismic waves do attenuate, large earthquakes (moment magnitude > 5) have the ability to produce observable shear waves. The second law of thermodynamics ensures a higher attenuation of shear wave reflected energy, this tends to impede the utilization of shear wave information for smaller earthquakes.
Crustal anisotropy
In the Earth's crust, anisotropy may be caused by preferentially aligned joints or microcracks, by layered bedding in sedimentary formations, or by highly foliated metamorphic rocks. Crustal anisotropy resulting from aligned cracks can be used to determine the state of stress in the crust, since in many cases, cracks are preferentially aligned with their flat faces oriented in the direction of minimum compressive stress. In active tectonic areas, such as near faults and volcanoes, anisotropy can be used to look for changes in preferred orientation of cracks that may indicate a rotation of the stress field.
Both seismic
P-waves
A P wave (primary wave or pressure wave) is one of the two main types of elastic body waves, called seismic waves in seismology. P waves travel faster than other seismic waves and hence are the first signal from an earthquake to arrive at any ...
and
S-waves
__NOTOC__
In seismology and other areas involving elastic waves, S waves, secondary waves, or shear waves (sometimes called elastic S waves) are a type of elastic wave and are one of the two main types of elastic body waves, so named because th ...
may exhibit anisotropy. For both, the anisotropy may appear as a (continuous) dependence of velocity upon the direction of propagation. For S-waves, it may also appear as a (discrete) dependence of velocity upon the direction of polarization. For a given direction of propagation in any homogeneous medium, only two polarization directions are allowed, with other polarizations decomposing trigonometrically into these two. Hence, shear waves naturally "split" into separate arrivals with these two polarizations; in optics this is called
birefringence.
Crustal anisotropy is very important in the production of oil reservoirs, as the seismically fast directions can indicate preferred directions of fluid flow.
In crustal geophysics, the anisotropy is usually weak; this enables a simplification of the expressions for seismic velocities and reflectivities, as functions of propagation (and polarization) direction. In the simplest geophysically plausible case, that of
polar anisotropy, the analysis is most conveniently done in terms of
Thomsen Parameters.
Mantle anisotropy
In the mantle, anisotropy is normally associated with crystals (mainly
olivine
The mineral olivine () is a magnesium iron silicate with the chemical formula . It is a type of nesosilicate or orthosilicate. The primary component of the Earth's upper mantle, it is a common mineral in Earth's subsurface, but weathers quickl ...
) aligned with the mantle flow direction called lattice preferred orientation (LPO). Due to their elongate crystalline structure, olivine crystals tend to align with the flow due to
mantle convection
Mantle convection is the very slow creeping motion of Earth's solid silicate mantle as convection currents carrying heat from the interior to the planet's surface.
The Earth's surface lithosphere rides atop the asthenosphere and the two for ...
or small scale convection. Anisotropy has long been used to argue whether plate tectonics is driven from below by mantle convection or from above by the plates, i.e. slab pull and ridge push.
The favored methods for detecting seismic anisotropy are
shear wave splitting
Shear wave splitting, also called seismic birefringence, is the phenomenon that occurs when a Polarization (waves), polarized shear wave enters an anisotropic medium (Fig. 1). The incident shear wave splits into two polarized shear waves (Fig. 2). ...
,
seismic tomography of surface waves and body waves, and converted-wave scattering in the context of a
receiver function. In shear-wave splitting, the S wave splits into two orthogonal polarizations, corresponding to the fastest and slowest wavespeeds in that medium for that propagation direction. The period range for mantle splitting studies is typically 5-25-sec. In seismic tomography, one must have a spatial distribution of seismic sources (earthquakes or man-made blasts) to generate waves at multiple wave-propagation azimuths through a 3-D medium. For receiver functions, the P-to-S converted wave displays harmonic variation with earthquake back azimuth when the material at depth is anisotopic. This method allows determination of layers of anisotropic material at depth beneath a station.
In the transition zone,
wadsleyite
Wadsleyite is an orthorhombic mineral with the formula β-(Mg,Fe)2SiO4. It was first found in nature in the Peace River meteorite from Alberta, Canada. It is formed by a phase transformation from olivine (α-(Mg,Fe)2SiO4) under increasing p ...
and/or
ringwoodite
Ringwoodite is a high-pressure phase of Mg2SiO4 (magnesium silicate) formed at high temperatures and pressures of the Earth's mantle between depth. It may also contain iron and hydrogen. It is polymorphous with the olivine phase forsterite (a ...
could be aligned in LPO. Below the
transition zone, the three main minerals,
periclase
Periclase is a magnesium mineral that occurs naturally in contact metamorphic rocks and is a major component of most basic refractory bricks. It is a cubic form of magnesium oxide ( Mg O). In nature it usually forms a solid solution with wüstit ...
,
silicate perovskite
Silicate perovskite is either (the magnesium end-member is called bridgmanite) or ( calcium silicate known as davemaoite) when arranged in a perovskite structure. Silicate perovskites are not stable at Earth's surface, and mainly exist in the l ...
(
bridgmanite
Silicate perovskite is either (the magnesium end-member is called bridgmanite) or (calcium silicate known as davemaoite) when arranged in a perovskite structure. Silicate perovskites are not stable at Earth's surface, and mainly exist in the l ...
), and
post-perovskite Post-perovskite (pPv) is a high-pressure phase of magnesium silicate (MgSiO3). It is composed of the prime oxide constituents of the Earth's rocky mantle (MgO and SiO2), and its pressure and temperature for stability imply that it is likely to occur ...
are all
anisotropic and could be generating anisotropy observed in the
D" region (a couple hundred kilometer thick layer about the core-mantle boundary).
References
{{Reflist
Sources
* Helbig, K., Thomsen, L., 75-plus years of anisotropy in exploration and reservoir seismics: A historical review of concepts and methods: Geophysics. VOL. 70, No. 6 (November–December 2005): p. 9ND–23ND http://www.geo.arizona.edu/geo5xx/geo596f/Readings/Helbig%20and%20Thomsen,%202005,%20historical%20review%20anisotropy%201.pdf
* Crampin, S., 1984, Evaluation of anisotropy by shear wave splitting: Applied Seismic Anisotropy: Theory, Background, and Field Studies, Geophysics Reprint series, 20, 23–33.
* Ikelle, L.T., Amundsen, L., Introduction to Petroleum Seismology, Investigations in Geophysics series No.12.
* Thomsen, L., 1986, Weak elastic anisotropy: Applied Seismic Anisotropy: Theory, Background, and Field Studies, Geophysics Reprint series, 20, 34–46
* Anderson et al., Oilfield Anisotropy: Its Origins and Electrical Characteristics: Oil field review, 48–56. https://www.slb.com/~/media/Files/resources/oilfield_review/ors94/1094/p48_56.pdf
* Thomsen, L., : Geophysics, 51, 1954–1966, Weak elastic anisotropy.
* Tsvankin, I., : Geophysics, 62, 1292-1309.1997, Anisotropic parameters and P-wave velocity for orthorhombic media.
* Tsvankin, I., Seismic signatures and analysis of reflection data in anisotropic media: Elsevier Science Publ, 2001,.
* Stephen A. H. and J-Michael K. GEOPHYSICS, VOL. 68, NO. 4, P1150–1160. Fracture characterization at Valhall: Application of P-wave amplitude variation with offset and azimuth (AVOA) analysis to a 3D ocean-bottom data set
* Tushar P. and Robert V. SPE 146668. Improved Reservoir Characterization through Estimation of Velocity Anisotropy in Shales.
* Jeffrey S., Rob R., Jean A., et al. www.cgg.com/technicalDocuments/cggv_0000000409.pdf Reducing Structural Uncertainties Through Anisotropic Prestack Depth Imaging: Examples from the Elgin/Franklin/Glenelg HP/HT Fields Area, Central North Sea
* Helbig, K., 1984, Shear waves – what they are and how they are and how they can be used: Applied Seismic Anisotropy: Theory, Background, and Field Studies, Geophysics Reprint series, 20, 5–22.
External links
* http://www1.gly.bris.ac.uk/~wookey/MMA/index.htm
* https://web.archive.org/web/20050909171919/http://geophysics.asu.edu/anisotropy/
* http://www.geo.arizona.edu/geo5xx/geo596f/Readings/Helbig%20and%20Thomsen,%202005,%20historical%20review%20anisotropy%201.pdf
* https://www.slb.com/~/media/Files/resources/oilfield_review/ors94/1094/p48_56.pdf
Elasticity (physics)
Petroleum geology
Geophysics