SlideShare a Scribd company logo
1 of 41
Thomsen Parameters and
their constraints for TI media
Kamal Aghazade
Part1:
Thomsen parameters
Part 2:
 physical constraints on Thomsen parameter ( )𝛿
Part 1:
Introducing Thomsen parameters to seismic
anisotropy
3
Seismic isotropy vs anisotropy
Isotropy comes from the Greek words isos (equal)
and tropos (way) and means uniform in all
directions. Isotropic materials like glass exhibit the
same material properties in all directions.
Physical properties are
not direction dependent.
Anisotropy: physical properties
are direction dependent.
Thomsen,
2014
Seismic anisotropy:
The dependent of seismic velocity
upon angle.
4
Beyond the simple definition of seismic anisotropy
a) isotropic, b) anisotropic media
Shear wave splitting in anisotropic media
https://en.wikipedia.org/wiki/Shear_wave_splitting
5
Beyond the simple definition of seismic anisotropy
isotropic wave equation simulation
6
anisotropic wave equation simulation
Lets go further: mathematical description Elastic media
Linear Elastic material and Hook’s law
7
seismic anisotropy media
(Transverse Isotropy)
Common coring configuration for anisotropy
measurements. Meléndez-Martínez (2014).
Rüger, 1997
http://www.glossary.oilfield.slb.com
8
Stiffness matrix in Transverse Isotropy 9
Christoffel equation
(eigenvalue
problem)
Directional dependence of three phase velocities
(Daley and Horn, 1977)
Algebraically complex 
10Transverse Isotropy
A simpler and meaningful structure for anisotropic analysis (Tsvankin, 2012).
Continue on next page 
11
The problems mentioned above are already solved by Thomsen (1986)
12
Weak elastic anisotropy and “anisotropies”
We can recast the equation for P-SV and SH waves velocity by
using notation involving:
1 – two elastic moduli (Vertical P- and S-waves velocity )
+
2- Three measures of anisotropy (We call them anisotropies).
131 – Simplifying Velocity equations for
types of waves
2 – anisotropies are non-dimensional.
So, one may speak of X percent
anisotropy.
3 – reduce to zero in the degenerate
case of isotropy. So, material with small
value of anisotropy may be denoted:
“Weakly anisotropic”
14
the algebraic complexity of new equations impedes a clear
understanding of their physical content. Progress may be
made. however, by observing that most
rocks are only weakly anisotropic. even though many of
their constituent minerals are highly anisotropic (Thomsen,
1986).
Thomsen’s revolution
14
Weak anisotropy
Based on laboratory data Thomsen (1986) showed
that most of evaluated rocks have anisotropy in the
weak to moderate (i.e. less than 0.2 ).
We use Taylor series to expand equations for
anisotropies and just retain only linear terms.
Now, we can talk about anisotropies
based on angles. 
15
Some notes about Thomesn parameters 16
By measurements at 0, 45 and 90 degrees we
have:
Lets consider error propagation in 𝛿 17
For more information about previous topics see:
Thomsen, L., [2014], Understanding Seismic Anisotropy
in Exploration and Exploitation, Second edition, SEG.
Tsvankin, I.,[2012] Seismic Signatures and Analysis of
Reflection Data in Anisotropic Media
,Third Edition, Society of exploration Geophysicists.
And related references.
Leon Thomsen
Ilya Tsvankin
Part 2:
physical constraints on 𝐶13 & 𝛿
for
transversely isotropic hydrocarbon source rocks
Yan et.al (2015)
20
21
22
Coordinate system for
following relations
Young’s modulus and Poisson’s ration for TI media (see King,1964) 23
Schema of deformation of vertical plug
(left) and horizontal plug (right)
of organic shale under axial
compressional testing. Dark grey
represents plugs
before deformation, and light grey
represents plugs after deformation
24
1- measurements
uncertainty
2 – material should not
be classified as a TI
medium
If a TI medium is
infinitely stronger in
the horizontal direction
compared with vertical
direction.
𝜐 𝐻𝑉 → 1
25
26
Setting constraint on 𝐶13
27On behavior of 𝛿
We have a constraint on 𝐶13
Substituting in 𝛿
Rewriting the
relation based on
other Thomsen’s
parameters
28
Where
29
Constraints on anellipticity parameter
“anellipticity”
parameter (η) that describes the degree of deviation from
elliptic anisotropy.
The anellipticity parameter η is important for anisotropic
seismic data processing because it determines the relation
between the normal moveout velocity and the horizontal
velocity (Tsvankin 2012).
30
Lower and upper bounds for 𝜂
Laboratory data and the constraints 31
several data points have negative values. The
corresponding 𝛿 values are above the high bound,
and they tend to have higher values of δ.
𝜐 𝐻𝐻
here are quite a few points with 𝜐 𝐻𝐻 > 𝜐 𝐻𝑉 . The
corresponding values are lower than the low bound,
and they tend to have lower values of δ .
𝐶13
About two-thirds of the data points lie in the center area,
where believed that all the hydrocarbon source rocks with
TI anisotropy should lie within.
There is some uncertainty in the middle part of the figure.
there are more data points lying below the low bound than
above the high bound.
Uncertainty in labratory velocity anisotropy measurements
Laboratory velocity anisotropy measurement on TI media
requires at least five velocity component measurements,
among which one velocity measurement must be made in
oblique direction.
32
Case I: Negative angle error
make about 20% of the data points lie
below the low bound
2 𝑜
Case II: Negative 5 𝑜
angle error
make about 62% of the data points
lie below the low bound
Case III: Positive 5 𝑜 angle error
make
about less than 8% of the data points lie above
the high bound
1) If the phase velocity in 45 𝑜
is
underestimated by 1%, 22% of the data
points move below the low bound.
2) If the phase velocity in 45 𝑜
is
overestimated by 1%, only one data
point moves above the high bound
33
difference between group and phase velocities
ray tracing of ultrasonic velocity measurement on the 45ºplug (left),
transmission time versus angle (right)
34
I. if the transducer is not wide enough (or
the sample is too long), the first
arriving energy might be missed by the
receiving transducer and the phase
velocity tends to be underestimated
Phase to group correction effect 35
Applications
The effect of the other Thomsen parameters on δ .
36
37
Histogram of from
laboratory velocity anisotropy measurement.
𝐶11−2𝐶66
here are only 2 data points with
One data point is due to data entry error and the other data
point is due to signals of substandard quality.
𝐶11−2𝐶66 < 0
37
The trends of the approximated bounds comply well with the laboratory
measured data if data points lying outside of the δ bounds are not displayed.
𝛿+
𝛿−
Ignore data
points outside
the bounds
38
39
Conclusions
 The physical constraints on the Thomsen parameter δ can help us understand the relation
between δ and the other Thomsen parameters.
 Generally, δ increases with ε and decreases with increasing γ . Variation of β0/α0 of the
hydrocarbon source rocks in a certain area is usually small so that δ is less sensitive to β0/α0.
 δ can be approximately predicted by the other Thomsen parameters .
 Using these constraints, there exist significant uncertainties in laboratory velocity anisotropy
measurement.
on Thomesn's strange anisotropy parameter

More Related Content

What's hot

Seismic data processing 15, kirchhof migration
Seismic data processing 15, kirchhof migrationSeismic data processing 15, kirchhof migration
Seismic data processing 15, kirchhof migrationAmin khalil
 
Seismic attribute analysis using complex trace analysis
Seismic attribute analysis using complex trace analysisSeismic attribute analysis using complex trace analysis
Seismic attribute analysis using complex trace analysisSomak Hajra
 
Basics of seismic interpretation
Basics of seismic interpretationBasics of seismic interpretation
Basics of seismic interpretationAmir I. Abdelaziz
 
Simple seismic processing workflow
Simple seismic processing workflowSimple seismic processing workflow
Simple seismic processing workflowAli M. Abdelsamad
 
Evaluation of low resistivity zones
Evaluation of low resistivity zonesEvaluation of low resistivity zones
Evaluation of low resistivity zonesShahnawaz Mustafa
 
inverse theory and inversion of seismic
inverse theory and inversion of seismic inverse theory and inversion of seismic
inverse theory and inversion of seismic Abdullah Abderahman
 
Using 3-D Seismic Attributes in Reservoir Characterization
Using 3-D Seismic Attributes in Reservoir CharacterizationUsing 3-D Seismic Attributes in Reservoir Characterization
Using 3-D Seismic Attributes in Reservoir Characterizationguest05b785
 
Structur Alanalysis
Structur AlanalysisStructur Alanalysis
Structur AlanalysisShah Naseer
 
Seismic interpretation - Fluvial Deltaic System
Seismic interpretation - Fluvial Deltaic SystemSeismic interpretation - Fluvial Deltaic System
Seismic interpretation - Fluvial Deltaic SystemAndi Anriansyah
 
2 d and 3d land seismic data acquisition and seismic data processing
2 d and 3d land seismic data acquisition and seismic data processing2 d and 3d land seismic data acquisition and seismic data processing
2 d and 3d land seismic data acquisition and seismic data processingAli Mahroug
 
Interpretation 23.12.13
Interpretation 23.12.13Interpretation 23.12.13
Interpretation 23.12.13Shashwat Sinha
 
Ultra sonic caliper log
Ultra sonic caliper logUltra sonic caliper log
Ultra sonic caliper logShah Naseer
 
ambiguity in geophysical interpretation_dileep p allavarapu
ambiguity in geophysical interpretation_dileep p allavarapuambiguity in geophysical interpretation_dileep p allavarapu
ambiguity in geophysical interpretation_dileep p allavarapuknigh7
 
Introduction to some basics
Introduction to some basicsIntroduction to some basics
Introduction to some basicsmohull
 
Geophysics-Gravity Method
Geophysics-Gravity Method Geophysics-Gravity Method
Geophysics-Gravity Method Shimul Al-razi
 
Introduction to seismic interpretation
Introduction to seismic interpretationIntroduction to seismic interpretation
Introduction to seismic interpretationAmir I. Abdelaziz
 
Role of Seismic Attributes in Petroleum Exploration_30May22.pptx
Role of Seismic Attributes in Petroleum Exploration_30May22.pptxRole of Seismic Attributes in Petroleum Exploration_30May22.pptx
Role of Seismic Attributes in Petroleum Exploration_30May22.pptxNagaLakshmiVasa
 

What's hot (20)

Seismic data processing 15, kirchhof migration
Seismic data processing 15, kirchhof migrationSeismic data processing 15, kirchhof migration
Seismic data processing 15, kirchhof migration
 
Seismic attribute analysis using complex trace analysis
Seismic attribute analysis using complex trace analysisSeismic attribute analysis using complex trace analysis
Seismic attribute analysis using complex trace analysis
 
Basics of seismic interpretation
Basics of seismic interpretationBasics of seismic interpretation
Basics of seismic interpretation
 
Principles of seismic data interpretation m.m.badawy
Principles of seismic data interpretation   m.m.badawyPrinciples of seismic data interpretation   m.m.badawy
Principles of seismic data interpretation m.m.badawy
 
Simple seismic processing workflow
Simple seismic processing workflowSimple seismic processing workflow
Simple seismic processing workflow
 
Evaluation of low resistivity zones
Evaluation of low resistivity zonesEvaluation of low resistivity zones
Evaluation of low resistivity zones
 
inverse theory and inversion of seismic
inverse theory and inversion of seismic inverse theory and inversion of seismic
inverse theory and inversion of seismic
 
Principles of seismic data processing m.m.badawy
Principles of seismic data processing m.m.badawyPrinciples of seismic data processing m.m.badawy
Principles of seismic data processing m.m.badawy
 
Using 3-D Seismic Attributes in Reservoir Characterization
Using 3-D Seismic Attributes in Reservoir CharacterizationUsing 3-D Seismic Attributes in Reservoir Characterization
Using 3-D Seismic Attributes in Reservoir Characterization
 
Structur Alanalysis
Structur AlanalysisStructur Alanalysis
Structur Alanalysis
 
Seismic interpretation - Fluvial Deltaic System
Seismic interpretation - Fluvial Deltaic SystemSeismic interpretation - Fluvial Deltaic System
Seismic interpretation - Fluvial Deltaic System
 
2 d and 3d land seismic data acquisition and seismic data processing
2 d and 3d land seismic data acquisition and seismic data processing2 d and 3d land seismic data acquisition and seismic data processing
2 d and 3d land seismic data acquisition and seismic data processing
 
Rock Physics: Definitions
Rock Physics: DefinitionsRock Physics: Definitions
Rock Physics: Definitions
 
Interpretation 23.12.13
Interpretation 23.12.13Interpretation 23.12.13
Interpretation 23.12.13
 
Ultra sonic caliper log
Ultra sonic caliper logUltra sonic caliper log
Ultra sonic caliper log
 
ambiguity in geophysical interpretation_dileep p allavarapu
ambiguity in geophysical interpretation_dileep p allavarapuambiguity in geophysical interpretation_dileep p allavarapu
ambiguity in geophysical interpretation_dileep p allavarapu
 
Introduction to some basics
Introduction to some basicsIntroduction to some basics
Introduction to some basics
 
Geophysics-Gravity Method
Geophysics-Gravity Method Geophysics-Gravity Method
Geophysics-Gravity Method
 
Introduction to seismic interpretation
Introduction to seismic interpretationIntroduction to seismic interpretation
Introduction to seismic interpretation
 
Role of Seismic Attributes in Petroleum Exploration_30May22.pptx
Role of Seismic Attributes in Petroleum Exploration_30May22.pptxRole of Seismic Attributes in Petroleum Exploration_30May22.pptx
Role of Seismic Attributes in Petroleum Exploration_30May22.pptx
 

Similar to on Thomesn's strange anisotropy parameter

Probing spin dynamics from the Mott insulating to the superfluid regime in a ...
Probing spin dynamics from the Mott insulating to the superfluid regime in a ...Probing spin dynamics from the Mott insulating to the superfluid regime in a ...
Probing spin dynamics from the Mott insulating to the superfluid regime in a ...Arijit Sharma
 
20-12-11 PhD Day Turin
20-12-11 PhD Day Turin20-12-11 PhD Day Turin
20-12-11 PhD Day TurinFabioDiRienzo
 
The Analytical/Numerical Relativity Interface behind Gravitational Waves: Lec...
The Analytical/Numerical Relativity Interface behind Gravitational Waves: Lec...The Analytical/Numerical Relativity Interface behind Gravitational Waves: Lec...
The Analytical/Numerical Relativity Interface behind Gravitational Waves: Lec...Lake Como School of Advanced Studies
 
Mate 280 characterization of powders and porous materials
Mate 280 characterization of powders and porous materialsMate 280 characterization of powders and porous materials
Mate 280 characterization of powders and porous materialsSami Ali
 
DISPERSION OF AEROSOLS IN ATMOSPHERIC FLUID FLOW
DISPERSION OF AEROSOLS IN ATMOSPHERIC FLUID FLOWDISPERSION OF AEROSOLS IN ATMOSPHERIC FLUID FLOW
DISPERSION OF AEROSOLS IN ATMOSPHERIC FLUID FLOWijscmcj
 
Goodness Dispersion Curves for Ultrasonic Guided Wave based SHM
Goodness Dispersion Curves for Ultrasonic Guided Wave based SHMGoodness Dispersion Curves for Ultrasonic Guided Wave based SHM
Goodness Dispersion Curves for Ultrasonic Guided Wave based SHMInnerspec Technologies
 
Basic Concepts of UV & IR Spectroscopy
Basic Concepts of UV & IR SpectroscopyBasic Concepts of UV & IR Spectroscopy
Basic Concepts of UV & IR SpectroscopyDrBasavarajaiahSm
 
The potential of terahertz imaging for cancer diagnosis
The potential of terahertz imaging for cancer diagnosisThe potential of terahertz imaging for cancer diagnosis
The potential of terahertz imaging for cancer diagnosisZahid Qaisar
 
20150930 Yokohama Protostellar discs
20150930 Yokohama Protostellar discs20150930 Yokohama Protostellar discs
20150930 Yokohama Protostellar discsGareth Murphy
 
MRI PHYSICS PART 3 Susceptibility-weighted images BY GKM .pptx
MRI PHYSICS PART 3 Susceptibility-weighted images BY GKM .pptxMRI PHYSICS PART 3 Susceptibility-weighted images BY GKM .pptx
MRI PHYSICS PART 3 Susceptibility-weighted images BY GKM .pptxGulshan Verma
 
23125135 the-basics-of-nmr
23125135 the-basics-of-nmr23125135 the-basics-of-nmr
23125135 the-basics-of-nmrdharma281276
 

Similar to on Thomesn's strange anisotropy parameter (20)

Zhang weglein-2008
Zhang weglein-2008Zhang weglein-2008
Zhang weglein-2008
 
Probing spin dynamics from the Mott insulating to the superfluid regime in a ...
Probing spin dynamics from the Mott insulating to the superfluid regime in a ...Probing spin dynamics from the Mott insulating to the superfluid regime in a ...
Probing spin dynamics from the Mott insulating to the superfluid regime in a ...
 
Swi seminar final
Swi seminar finalSwi seminar final
Swi seminar final
 
20-12-11 PhD Day Turin
20-12-11 PhD Day Turin20-12-11 PhD Day Turin
20-12-11 PhD Day Turin
 
The Analytical/Numerical Relativity Interface behind Gravitational Waves: Lec...
The Analytical/Numerical Relativity Interface behind Gravitational Waves: Lec...The Analytical/Numerical Relativity Interface behind Gravitational Waves: Lec...
The Analytical/Numerical Relativity Interface behind Gravitational Waves: Lec...
 
Mate 280 characterization of powders and porous materials
Mate 280 characterization of powders and porous materialsMate 280 characterization of powders and porous materials
Mate 280 characterization of powders and porous materials
 
Quantum theory research overview
Quantum theory research overview Quantum theory research overview
Quantum theory research overview
 
Surface Wave Tomography
Surface Wave TomographySurface Wave Tomography
Surface Wave Tomography
 
Surface Wave Tomography
Surface Wave TomographySurface Wave Tomography
Surface Wave Tomography
 
DISPERSION OF AEROSOLS IN ATMOSPHERIC FLUID FLOW
DISPERSION OF AEROSOLS IN ATMOSPHERIC FLUID FLOWDISPERSION OF AEROSOLS IN ATMOSPHERIC FLUID FLOW
DISPERSION OF AEROSOLS IN ATMOSPHERIC FLUID FLOW
 
NMRpresentationCombined
NMRpresentationCombinedNMRpresentationCombined
NMRpresentationCombined
 
Goodness Dispersion Curves for Ultrasonic Guided Wave based SHM
Goodness Dispersion Curves for Ultrasonic Guided Wave based SHMGoodness Dispersion Curves for Ultrasonic Guided Wave based SHM
Goodness Dispersion Curves for Ultrasonic Guided Wave based SHM
 
Aadrsh kumar tiwari bbau
Aadrsh kumar tiwari bbauAadrsh kumar tiwari bbau
Aadrsh kumar tiwari bbau
 
REU_paper
REU_paperREU_paper
REU_paper
 
Basic Concepts of UV & IR Spectroscopy
Basic Concepts of UV & IR SpectroscopyBasic Concepts of UV & IR Spectroscopy
Basic Concepts of UV & IR Spectroscopy
 
The potential of terahertz imaging for cancer diagnosis
The potential of terahertz imaging for cancer diagnosisThe potential of terahertz imaging for cancer diagnosis
The potential of terahertz imaging for cancer diagnosis
 
20150930 Yokohama Protostellar discs
20150930 Yokohama Protostellar discs20150930 Yokohama Protostellar discs
20150930 Yokohama Protostellar discs
 
MRI PHYSICS PART 3 Susceptibility-weighted images BY GKM .pptx
MRI PHYSICS PART 3 Susceptibility-weighted images BY GKM .pptxMRI PHYSICS PART 3 Susceptibility-weighted images BY GKM .pptx
MRI PHYSICS PART 3 Susceptibility-weighted images BY GKM .pptx
 
23125135 the-basics-of-nmr
23125135 the-basics-of-nmr23125135 the-basics-of-nmr
23125135 the-basics-of-nmr
 
Spectrophotometry in Nutshell
Spectrophotometry in NutshellSpectrophotometry in Nutshell
Spectrophotometry in Nutshell
 

More from Inistute of Geophysics, Tehran university , Tehran/ iran (10)

Presentation_Wind_meeting.pptx
Presentation_Wind_meeting.pptxPresentation_Wind_meeting.pptx
Presentation_Wind_meeting.pptx
 
EAGE_prsentation_Anderson.pptx
EAGE_prsentation_Anderson.pptxEAGE_prsentation_Anderson.pptx
EAGE_prsentation_Anderson.pptx
 
Frequency-domain Finite-difference modelling by plane wave interpolation
Frequency-domain Finite-difference modelling by plane wave interpolationFrequency-domain Finite-difference modelling by plane wave interpolation
Frequency-domain Finite-difference modelling by plane wave interpolation
 
Colour role in seismic interpretation
Colour role in seismic interpretationColour role in seismic interpretation
Colour role in seismic interpretation
 
Seismic Migration
Seismic MigrationSeismic Migration
Seismic Migration
 
Least Squares method
Least Squares methodLeast Squares method
Least Squares method
 
Warrp sesimic_ aghazde
Warrp sesimic_ aghazdeWarrp sesimic_ aghazde
Warrp sesimic_ aghazde
 
Seismic velocity analysis _ aghazade
Seismic velocity analysis _ aghazadeSeismic velocity analysis _ aghazade
Seismic velocity analysis _ aghazade
 
Levenberg - Marquardt (LM) algorithm_ aghazade
Levenberg - Marquardt (LM) algorithm_ aghazadeLevenberg - Marquardt (LM) algorithm_ aghazade
Levenberg - Marquardt (LM) algorithm_ aghazade
 
Huber(FWI) _aghazade -- amini
Huber(FWI) _aghazade -- aminiHuber(FWI) _aghazade -- amini
Huber(FWI) _aghazade -- amini
 

Recently uploaded

GBSN - Microbiology (Unit 2)
GBSN - Microbiology (Unit 2)GBSN - Microbiology (Unit 2)
GBSN - Microbiology (Unit 2)Areesha Ahmad
 
Presentation Vikram Lander by Vedansh Gupta.pptx
Presentation Vikram Lander by Vedansh Gupta.pptxPresentation Vikram Lander by Vedansh Gupta.pptx
Presentation Vikram Lander by Vedansh Gupta.pptxgindu3009
 
Pests of cotton_Sucking_Pests_Dr.UPR.pdf
Pests of cotton_Sucking_Pests_Dr.UPR.pdfPests of cotton_Sucking_Pests_Dr.UPR.pdf
Pests of cotton_Sucking_Pests_Dr.UPR.pdfPirithiRaju
 
Chemistry 4th semester series (krishna).pdf
Chemistry 4th semester series (krishna).pdfChemistry 4th semester series (krishna).pdf
Chemistry 4th semester series (krishna).pdfSumit Kumar yadav
 
Formation of low mass protostars and their circumstellar disks
Formation of low mass protostars and their circumstellar disksFormation of low mass protostars and their circumstellar disks
Formation of low mass protostars and their circumstellar disksSérgio Sacani
 
GBSN - Microbiology (Unit 1)
GBSN - Microbiology (Unit 1)GBSN - Microbiology (Unit 1)
GBSN - Microbiology (Unit 1)Areesha Ahmad
 
PossibleEoarcheanRecordsoftheGeomagneticFieldPreservedintheIsuaSupracrustalBe...
PossibleEoarcheanRecordsoftheGeomagneticFieldPreservedintheIsuaSupracrustalBe...PossibleEoarcheanRecordsoftheGeomagneticFieldPreservedintheIsuaSupracrustalBe...
PossibleEoarcheanRecordsoftheGeomagneticFieldPreservedintheIsuaSupracrustalBe...Sérgio Sacani
 
Botany krishna series 2nd semester Only Mcq type questions
Botany krishna series 2nd semester Only Mcq type questionsBotany krishna series 2nd semester Only Mcq type questions
Botany krishna series 2nd semester Only Mcq type questionsSumit Kumar yadav
 
Disentangling the origin of chemical differences using GHOST
Disentangling the origin of chemical differences using GHOSTDisentangling the origin of chemical differences using GHOST
Disentangling the origin of chemical differences using GHOSTSérgio Sacani
 
Botany 4th semester series (krishna).pdf
Botany 4th semester series (krishna).pdfBotany 4th semester series (krishna).pdf
Botany 4th semester series (krishna).pdfSumit Kumar yadav
 
9654467111 Call Girls In Raj Nagar Delhi Short 1500 Night 6000
9654467111 Call Girls In Raj Nagar Delhi Short 1500 Night 60009654467111 Call Girls In Raj Nagar Delhi Short 1500 Night 6000
9654467111 Call Girls In Raj Nagar Delhi Short 1500 Night 6000Sapana Sha
 
DIFFERENCE IN BACK CROSS AND TEST CROSS
DIFFERENCE IN  BACK CROSS AND TEST CROSSDIFFERENCE IN  BACK CROSS AND TEST CROSS
DIFFERENCE IN BACK CROSS AND TEST CROSSLeenakshiTyagi
 
Biopesticide (2).pptx .This slides helps to know the different types of biop...
Biopesticide (2).pptx  .This slides helps to know the different types of biop...Biopesticide (2).pptx  .This slides helps to know the different types of biop...
Biopesticide (2).pptx .This slides helps to know the different types of biop...RohitNehra6
 
Traditional Agroforestry System in India- Shifting Cultivation, Taungya, Home...
Traditional Agroforestry System in India- Shifting Cultivation, Taungya, Home...Traditional Agroforestry System in India- Shifting Cultivation, Taungya, Home...
Traditional Agroforestry System in India- Shifting Cultivation, Taungya, Home...jana861314
 
TEST BANK For Radiologic Science for Technologists, 12th Edition by Stewart C...
TEST BANK For Radiologic Science for Technologists, 12th Edition by Stewart C...TEST BANK For Radiologic Science for Technologists, 12th Edition by Stewart C...
TEST BANK For Radiologic Science for Technologists, 12th Edition by Stewart C...ssifa0344
 
GFP in rDNA Technology (Biotechnology).pptx
GFP in rDNA Technology (Biotechnology).pptxGFP in rDNA Technology (Biotechnology).pptx
GFP in rDNA Technology (Biotechnology).pptxAleenaTreesaSaji
 
Stunning ➥8448380779▻ Call Girls In Panchshil Enclave Delhi NCR
Stunning ➥8448380779▻ Call Girls In Panchshil Enclave Delhi NCRStunning ➥8448380779▻ Call Girls In Panchshil Enclave Delhi NCR
Stunning ➥8448380779▻ Call Girls In Panchshil Enclave Delhi NCRDelhi Call girls
 
Chromatin Structure | EUCHROMATIN | HETEROCHROMATIN
Chromatin Structure | EUCHROMATIN | HETEROCHROMATINChromatin Structure | EUCHROMATIN | HETEROCHROMATIN
Chromatin Structure | EUCHROMATIN | HETEROCHROMATINsankalpkumarsahoo174
 

Recently uploaded (20)

GBSN - Microbiology (Unit 2)
GBSN - Microbiology (Unit 2)GBSN - Microbiology (Unit 2)
GBSN - Microbiology (Unit 2)
 
Presentation Vikram Lander by Vedansh Gupta.pptx
Presentation Vikram Lander by Vedansh Gupta.pptxPresentation Vikram Lander by Vedansh Gupta.pptx
Presentation Vikram Lander by Vedansh Gupta.pptx
 
Pests of cotton_Sucking_Pests_Dr.UPR.pdf
Pests of cotton_Sucking_Pests_Dr.UPR.pdfPests of cotton_Sucking_Pests_Dr.UPR.pdf
Pests of cotton_Sucking_Pests_Dr.UPR.pdf
 
Chemistry 4th semester series (krishna).pdf
Chemistry 4th semester series (krishna).pdfChemistry 4th semester series (krishna).pdf
Chemistry 4th semester series (krishna).pdf
 
Formation of low mass protostars and their circumstellar disks
Formation of low mass protostars and their circumstellar disksFormation of low mass protostars and their circumstellar disks
Formation of low mass protostars and their circumstellar disks
 
GBSN - Microbiology (Unit 1)
GBSN - Microbiology (Unit 1)GBSN - Microbiology (Unit 1)
GBSN - Microbiology (Unit 1)
 
PossibleEoarcheanRecordsoftheGeomagneticFieldPreservedintheIsuaSupracrustalBe...
PossibleEoarcheanRecordsoftheGeomagneticFieldPreservedintheIsuaSupracrustalBe...PossibleEoarcheanRecordsoftheGeomagneticFieldPreservedintheIsuaSupracrustalBe...
PossibleEoarcheanRecordsoftheGeomagneticFieldPreservedintheIsuaSupracrustalBe...
 
CELL -Structural and Functional unit of life.pdf
CELL -Structural and Functional unit of life.pdfCELL -Structural and Functional unit of life.pdf
CELL -Structural and Functional unit of life.pdf
 
Botany krishna series 2nd semester Only Mcq type questions
Botany krishna series 2nd semester Only Mcq type questionsBotany krishna series 2nd semester Only Mcq type questions
Botany krishna series 2nd semester Only Mcq type questions
 
Disentangling the origin of chemical differences using GHOST
Disentangling the origin of chemical differences using GHOSTDisentangling the origin of chemical differences using GHOST
Disentangling the origin of chemical differences using GHOST
 
Botany 4th semester series (krishna).pdf
Botany 4th semester series (krishna).pdfBotany 4th semester series (krishna).pdf
Botany 4th semester series (krishna).pdf
 
The Philosophy of Science
The Philosophy of ScienceThe Philosophy of Science
The Philosophy of Science
 
9654467111 Call Girls In Raj Nagar Delhi Short 1500 Night 6000
9654467111 Call Girls In Raj Nagar Delhi Short 1500 Night 60009654467111 Call Girls In Raj Nagar Delhi Short 1500 Night 6000
9654467111 Call Girls In Raj Nagar Delhi Short 1500 Night 6000
 
DIFFERENCE IN BACK CROSS AND TEST CROSS
DIFFERENCE IN  BACK CROSS AND TEST CROSSDIFFERENCE IN  BACK CROSS AND TEST CROSS
DIFFERENCE IN BACK CROSS AND TEST CROSS
 
Biopesticide (2).pptx .This slides helps to know the different types of biop...
Biopesticide (2).pptx  .This slides helps to know the different types of biop...Biopesticide (2).pptx  .This slides helps to know the different types of biop...
Biopesticide (2).pptx .This slides helps to know the different types of biop...
 
Traditional Agroforestry System in India- Shifting Cultivation, Taungya, Home...
Traditional Agroforestry System in India- Shifting Cultivation, Taungya, Home...Traditional Agroforestry System in India- Shifting Cultivation, Taungya, Home...
Traditional Agroforestry System in India- Shifting Cultivation, Taungya, Home...
 
TEST BANK For Radiologic Science for Technologists, 12th Edition by Stewart C...
TEST BANK For Radiologic Science for Technologists, 12th Edition by Stewart C...TEST BANK For Radiologic Science for Technologists, 12th Edition by Stewart C...
TEST BANK For Radiologic Science for Technologists, 12th Edition by Stewart C...
 
GFP in rDNA Technology (Biotechnology).pptx
GFP in rDNA Technology (Biotechnology).pptxGFP in rDNA Technology (Biotechnology).pptx
GFP in rDNA Technology (Biotechnology).pptx
 
Stunning ➥8448380779▻ Call Girls In Panchshil Enclave Delhi NCR
Stunning ➥8448380779▻ Call Girls In Panchshil Enclave Delhi NCRStunning ➥8448380779▻ Call Girls In Panchshil Enclave Delhi NCR
Stunning ➥8448380779▻ Call Girls In Panchshil Enclave Delhi NCR
 
Chromatin Structure | EUCHROMATIN | HETEROCHROMATIN
Chromatin Structure | EUCHROMATIN | HETEROCHROMATINChromatin Structure | EUCHROMATIN | HETEROCHROMATIN
Chromatin Structure | EUCHROMATIN | HETEROCHROMATIN
 

on Thomesn's strange anisotropy parameter

  • 1. Thomsen Parameters and their constraints for TI media Kamal Aghazade
  • 2. Part1: Thomsen parameters Part 2:  physical constraints on Thomsen parameter ( )𝛿
  • 3. Part 1: Introducing Thomsen parameters to seismic anisotropy 3
  • 4. Seismic isotropy vs anisotropy Isotropy comes from the Greek words isos (equal) and tropos (way) and means uniform in all directions. Isotropic materials like glass exhibit the same material properties in all directions. Physical properties are not direction dependent. Anisotropy: physical properties are direction dependent. Thomsen, 2014 Seismic anisotropy: The dependent of seismic velocity upon angle. 4
  • 5. Beyond the simple definition of seismic anisotropy a) isotropic, b) anisotropic media Shear wave splitting in anisotropic media https://en.wikipedia.org/wiki/Shear_wave_splitting 5
  • 6. Beyond the simple definition of seismic anisotropy isotropic wave equation simulation 6 anisotropic wave equation simulation
  • 7. Lets go further: mathematical description Elastic media Linear Elastic material and Hook’s law 7
  • 8. seismic anisotropy media (Transverse Isotropy) Common coring configuration for anisotropy measurements. Meléndez-Martínez (2014). Rüger, 1997 http://www.glossary.oilfield.slb.com 8
  • 9. Stiffness matrix in Transverse Isotropy 9
  • 10. Christoffel equation (eigenvalue problem) Directional dependence of three phase velocities (Daley and Horn, 1977) Algebraically complex  10Transverse Isotropy
  • 11. A simpler and meaningful structure for anisotropic analysis (Tsvankin, 2012). Continue on next page  11
  • 12. The problems mentioned above are already solved by Thomsen (1986) 12
  • 13. Weak elastic anisotropy and “anisotropies” We can recast the equation for P-SV and SH waves velocity by using notation involving: 1 – two elastic moduli (Vertical P- and S-waves velocity ) + 2- Three measures of anisotropy (We call them anisotropies). 131 – Simplifying Velocity equations for types of waves 2 – anisotropies are non-dimensional. So, one may speak of X percent anisotropy. 3 – reduce to zero in the degenerate case of isotropy. So, material with small value of anisotropy may be denoted: “Weakly anisotropic”
  • 14. 14 the algebraic complexity of new equations impedes a clear understanding of their physical content. Progress may be made. however, by observing that most rocks are only weakly anisotropic. even though many of their constituent minerals are highly anisotropic (Thomsen, 1986). Thomsen’s revolution 14
  • 15. Weak anisotropy Based on laboratory data Thomsen (1986) showed that most of evaluated rocks have anisotropy in the weak to moderate (i.e. less than 0.2 ). We use Taylor series to expand equations for anisotropies and just retain only linear terms. Now, we can talk about anisotropies based on angles.  15
  • 16. Some notes about Thomesn parameters 16
  • 17. By measurements at 0, 45 and 90 degrees we have: Lets consider error propagation in 𝛿 17
  • 18. For more information about previous topics see: Thomsen, L., [2014], Understanding Seismic Anisotropy in Exploration and Exploitation, Second edition, SEG. Tsvankin, I.,[2012] Seismic Signatures and Analysis of Reflection Data in Anisotropic Media ,Third Edition, Society of exploration Geophysicists. And related references. Leon Thomsen Ilya Tsvankin
  • 19. Part 2: physical constraints on 𝐶13 & 𝛿 for transversely isotropic hydrocarbon source rocks Yan et.al (2015)
  • 20. 20
  • 21. 21
  • 22. 22
  • 23. Coordinate system for following relations Young’s modulus and Poisson’s ration for TI media (see King,1964) 23
  • 24. Schema of deformation of vertical plug (left) and horizontal plug (right) of organic shale under axial compressional testing. Dark grey represents plugs before deformation, and light grey represents plugs after deformation 24
  • 25. 1- measurements uncertainty 2 – material should not be classified as a TI medium If a TI medium is infinitely stronger in the horizontal direction compared with vertical direction. 𝜐 𝐻𝑉 → 1 25
  • 28. We have a constraint on 𝐶13 Substituting in 𝛿 Rewriting the relation based on other Thomsen’s parameters 28
  • 30. Constraints on anellipticity parameter “anellipticity” parameter (η) that describes the degree of deviation from elliptic anisotropy. The anellipticity parameter η is important for anisotropic seismic data processing because it determines the relation between the normal moveout velocity and the horizontal velocity (Tsvankin 2012). 30 Lower and upper bounds for 𝜂
  • 31. Laboratory data and the constraints 31 several data points have negative values. The corresponding 𝛿 values are above the high bound, and they tend to have higher values of δ. 𝜐 𝐻𝐻 here are quite a few points with 𝜐 𝐻𝐻 > 𝜐 𝐻𝑉 . The corresponding values are lower than the low bound, and they tend to have lower values of δ . 𝐶13 About two-thirds of the data points lie in the center area, where believed that all the hydrocarbon source rocks with TI anisotropy should lie within. There is some uncertainty in the middle part of the figure. there are more data points lying below the low bound than above the high bound.
  • 32. Uncertainty in labratory velocity anisotropy measurements Laboratory velocity anisotropy measurement on TI media requires at least five velocity component measurements, among which one velocity measurement must be made in oblique direction. 32
  • 33. Case I: Negative angle error make about 20% of the data points lie below the low bound 2 𝑜 Case II: Negative 5 𝑜 angle error make about 62% of the data points lie below the low bound Case III: Positive 5 𝑜 angle error make about less than 8% of the data points lie above the high bound 1) If the phase velocity in 45 𝑜 is underestimated by 1%, 22% of the data points move below the low bound. 2) If the phase velocity in 45 𝑜 is overestimated by 1%, only one data point moves above the high bound 33
  • 34. difference between group and phase velocities ray tracing of ultrasonic velocity measurement on the 45ºplug (left), transmission time versus angle (right) 34 I. if the transducer is not wide enough (or the sample is too long), the first arriving energy might be missed by the receiving transducer and the phase velocity tends to be underestimated
  • 35. Phase to group correction effect 35
  • 36. Applications The effect of the other Thomsen parameters on δ . 36
  • 37. 37 Histogram of from laboratory velocity anisotropy measurement. 𝐶11−2𝐶66 here are only 2 data points with One data point is due to data entry error and the other data point is due to signals of substandard quality. 𝐶11−2𝐶66 < 0 37
  • 38. The trends of the approximated bounds comply well with the laboratory measured data if data points lying outside of the δ bounds are not displayed. 𝛿+ 𝛿− Ignore data points outside the bounds 38
  • 39. 39
  • 40. Conclusions  The physical constraints on the Thomsen parameter δ can help us understand the relation between δ and the other Thomsen parameters.  Generally, δ increases with ε and decreases with increasing γ . Variation of β0/α0 of the hydrocarbon source rocks in a certain area is usually small so that δ is less sensitive to β0/α0.  δ can be approximately predicted by the other Thomsen parameters .  Using these constraints, there exist significant uncertainties in laboratory velocity anisotropy measurement.