Vježbe: What Is A Fluid Inclusion?
Vježbe: What Is A Fluid Inclusion?
Vježbe: What Is A Fluid Inclusion?
Literatura
What is a Fluid Inclusion?
Bakker, R.J., 2003. Package FLUIDS 1. Computer programs for analysis of fluid
inclusion data and for modelling bulk fluid properties. Chemical Geology, 194, 3–
23. Cavity in a mineral that may contain 1 or more phases
vapor (V) - H2O, CO2, CH4, N2, H2S
http://fluids.unileoben.ac.at/Home.html liquid (L)- H2O, CO2, petroleum
Henry Clifton Sorby
(1826-1908) solid (S) - NaCl, KCl, hematite, anhydrite, muscovite,
English microscopist and magnetite, carbonates
Brown P.E., 1989. FLINCOR; a microcomputer program for the reduction and geologist
investigation of fluid-inclusion data. American Mineralogist, 74/11-12, 1390-1393. 1858 On the Microscopical
Structure of Crystals
Roedder, E., 1984. Fluid inclusions. Mineralogical Society America, Review in (Quart. Journ. Geol. Soc.)
Mineralogy 12, Washington, 644 pp.
V
Shepherd, T.J., Rankin, A.H., Alderton, D.H.M., 1985. A practical guide to fluid
inclusion studies. Blackie and Son Ltd, Glasgow, 239 pp. L
The liquid of the inclusion is normally an aqueous solutions with dissolved ions of Na+, Cl-, Ca2+,
Mg2+, SO42-, HCO32-, CO32-
The concentration of the salts ranges from <1 wt. % to >50 wt. %
1
Top 10 minerals Abundance of FI-s in Average No of FI-s occupying
1. Quartz single crystal if the total FI-s size 0.1% volume
Occurrence and distribution
2. Fluorite population occupy 0.1% 1 mm 1
3. Halite
volume 100 µm 103
4. Calcite
Earth Crust – magmatic, metamorphic and sedimentary rocks, 5. Apatite Size and abundance 10 µm 106
ore deposits, fault zones 6. Dolomite 1 µm 109
7. Sphalerite
8. Barite
9. Topaz
10. Cassiterite
Size of fluid
inclusions in minerals
What can I study? Which type of material do I use? What can I get from this method?
Ore deposits (hydrothermal, porphyry, Large, well developed crystals in
skarn) vughs and druses, all phases of vein
developing, alterations
1. Composition of fluid form which mineral precipitated
Igneous rocks Quartz (remnant of magmatic fluid or
late stage hydrothermal circulation) 2. Changes of fluid composition during precipitation (mixing, dilution,
boiling, cooling)
Apatite from carbonatite
Phenocrysts from volcanic 3. Minimum temperature and pressure at the time of precipitation
Pegmatites Quartz, beril, tourmaline 4. True temperature and pressure applying pressure correction (i.e.
independent geothermometer or boiling fluid)
5. Depth of formation (i.e. overlying deposits)
Metamorphic rocks Quartz from veins, pods, segregations
2
What type of samples do I need? How to prepare doubly polished wafers?
P = Primary
PS = PseudoSecondary
S = Secondary
inclusions
Roedder, 1984,
American Mineralogist
3
Primary FI-s classification
Mechanisms of trapping of primary FI-s
(a) Diagnostic criteria for classifying fluid inclusions as primary (after Roedder, 1979)
(b) Different occurrences of primary fluid inclusions in relation to growth zoning (compilation) A dendritic growth
(from: Van den Kerkhof & Hein, 2001, Lithos 55, 27-47) B partial disolution
C between growth spirals
D sub-parallel block growth
E fracture during growth
F foreign object
a) main distinction is
made between transgranular,
intergranular, and intragranular
inclusions
4
Working example 1: FI-s content - classification
(from: Sheperd, 1985)
Determine relative age
(Hansteen, 1988)
Single-phase Two-phase
Liquid or Liquid-rich
growth zoning
CH4, CO2
Working example 2:
Classify the inclusions
Basic Assumptions
HOWEVER
15 µm 10 µm 10 µm 1. Bulk leakage
2. Leakage through diffusion
3. Stretching
4. Re-equilibration
5. Necking-down
6. Migration
5
Recognition
Necking-down
1. Bulk leakage – soft, easy cleaved material – visual estimation on constancy of L to V
ratio, reproducibility of thermometric results
2. Leakage through diffusion – proved as extremely low (H+ found by Raman)
3. Stretching – change of volume without changing of composition
4. Re-equilibration – irregularly shaped FI-s tend to change shape and morphology
during time into negative crystal of spherical – insoluble minerals as quartz -
isochemical at constant volume
5. Necking-down – if original FI is large, flat and irregular will split into many small but
more regular FI-s. If the process occur during homogenous stage – good, preserved
L to V ratio, if not L to V ratio disturbed – erroneous values time
6. Migration – in a thermal gradient – important for water soluble materials
Working example 3:
What happened? Estimation of volume fractions relative to FI-s size – cylindrical FI-s
Why?
How do you know?
6
Estimation of volume fractions of vapor-rich
inclusions – here comes the problem
Estimation of volume fractions of spherical inclusions
Cylindrical
Trapping mode
Spherical
7
Naziv faze Sastav Stupanj Kristalni Habitus
Accidental trapping of solids refleksije sustav
Tekući ugljik- CO2 (l) 1,195
dioksid
Led H2O 1,31 Heksagonski Zaobljen,anizotropan,
izgleda izotropno
Hidrohalit NaCl×2H2O 1,41 Monoklinski Sitna zrnca
svjetlucavog izgleda
Silvin KCl 1,49 Kubični Kockice
In most cases we can assue that density of the vapor phase is zero, thus:
1 0 15 25 wt % NaCl equ.
ρTOT = ρLF
Technical specifications
Degree of fill (F)
8
Eutectic properties of
SALT SOLUTIONS Diagnostic features
(Hein, 1990)
H2O system
vapor
liquid
cooling
Two phase L+ V
phase boundaries
Tripple point
Critical point
9
H2O system H2O system
Principle of fluid inclusion geothermometry PT diagram Consider an inclusion trapped at a given
for pure water temperature and pressure (Tt, Pt)
1.0 1.0
Liquid Liquid
water Isochore Isochore
(g/cc) Pt (g/cc)
Pressure KBar
Pressure KBar
Critical Critical
point point
0.5 0.5
0
0. 6
0. 6
0.80
0.80
0.95
0.95
e
ling
Dry steam curv
Boi curve -Vapour L= V Vapour
0 L= V 0
50 150 350 Temp oC 50 150 Tt 350 Temp oC
On cooling, the inclusion follows an isochoric PT path Beyond this point the inclusion cools along the L=V
until it meets the L=V curve curve and a vapour bubble nucleates
1.0 1.0
Liquid Liquid
Isochore Isochore
Pt (g/cc) Pt (g/cc)
Pressure KBar
Pressure KBar
Critical Critical
point point
0.5 0.5
0
0
0. 6
0. 6
0.80
0.80
0.95
0.95
e e
curv curv
L= V Vapour L= V Vapour
0 0
50 150 Tt 350 Temp oC 50 150 Tt 350 Temp oC
10
H2O system H2O system
Continued cooling results in further shrinkage of liquid On heating along the V/L curve, the liquid expands and
and growth of the vapour bubble the bubble shrinks
1.0 1.0
Liquid Liquid
Isochore Isochore
Pt (g/cc) Pt (g/cc)
Pressure KBar
Pressure KBar
Critical Critical
point point
0.5 0.5
0
0. 6
0. 6
0.80
0.80
0.95
0.95
e e
curv curv
L= V Vapour L= V Vapour
0 0
50 150 Tt 350 Temp oC 50 150 Tt 350 Temp oC
1.0 1.0
Liquid Liquid
Isochore Isochore
Pt (g/cc) Pt (g/cc)
Pressure KBar
Pressure KBar
Critical Critical
point point
ore
0.5 0.5
Isoch
0
0
0. 6
0. 6
0.80
0.80
0.95
0.95
e e
curv curv
L= V Vapour L= V Vapour
0 0
50 150 Th Tt 350 Temp oC 50 150 Th Tt 350 Temp oC
11
H2O system H2O system
With continued heating the inclusion follows the original If Pt is known, or estimated, the trapping temperature
isochore (Tt) can be determined
1.0 1.0
Liquid Liquid
Isochore Isochore
Pt (g/cc) Pt (g/cc)
Pressure KBar
Pressure KBar
Critical Critical
point point
0.5 0.5
0
0. 6
0. 6
0.80
0.80
0.95
0.95
e e
curv curv
L= V Vapour L= V Vapour
0 0
50 150 Th Tt 350 Temp oC 50 150 Th Tt 350 Temp oC
The difference between Th and Tt is known as the Principle of fluid inclusion geothermometry based on
Pressure Correction PVT diagram for pure water
1.0 1.0
Liquid Liquid
Isochore Isochore
Pt (g/cc) Pt (g/cc)
Pressure KBar
Pressure KBar
Critical Critical
point point
0.5 0.5
0
0
0. 6
0. 6
0.80
0.80
0.95
0.95
e e
curv curv
L= V Vapour L= V Vapour
0 0
50 150 Th Tt 350 Temp oC 50 150 Th Tt 350 Temp oC
12
H2O system
SC
F
H2O isochores modified from Fischer (1976).
Enlargement of the low pressure region, in the box, appears to the right. 3D phase diagram –
two component system Liquid
So
l A
+L
Sol B + L
Sol A + sol B
So
lA
+
V
r
Va po
po
r va
A +V
com So
lB lB
pos + So
itio V
n
Heating
Cooling
13
350oC 350oC
25oC 25oC
0oC 0oC
-100oC -100oC
350oC 350oC
25oC 25oC
0oC 0oC
-100oC -100oC
14
350oC 350oC
25oC 25oC
0oC 0oC
-100oC -100oC
350oC 350oC
Temperature of Homogenisation Th
25oC 25oC
0oC 0oC
-100oC -100oC
15
350oC 350oC
25oC 25oC
0oC 0oC
-100oC -100oC
25oC
0oC Heating
Cooling
-100oC
16
350oC 350oC
25oC 25oC
0oC 0oC
Eutectic properties of
SALT SOLUTIONS
350oC
25oC
0oC
17
350oC 350oC
25oC 25oC
0oC 0oC
-100oC -100oC
350oC 350oC
25oC 25oC
0oC 0oC
-100oC -100oC
18
350oC 350oC
25oC 25oC
0oC 0oC
25 25
Temperature oC
0 0.1oC 0 0.1oC
NaCl.2H2O NaCl.2H2O
Ice + L+ V +L+ V Ice + L+ V +L+ V
-20.8oC -20.8oC
-25 -25
-50 -50
0 10 20 30 0 10 20 30
Weight % NaCl Weight % NaCl
19
NaCl-H2O system NaCl-H2O system
Rapid cooling below the eutectic temperature (Te) On heating first melting (Tfm) occurs at -20.8
is usually needed before the inclusion freezes (Te), evident by “unlocking” of the vapour bubble
25 25
Temperature oC
0 0.1oC 0 0.1oC
NaCl.2H2O NaCl.2H2O
Ice + L+ V +L+ V Tfm Ice + L+ V +L+ V
-20.8oC -20.8oC
-25 -25
-50 -50
0 10 20 30 0 10 20 30
Weight % NaCl Weight % NaCl
25 25
Temperature oC
0 0.1oC 0 0.1oC
Tm(ice) Tm(ice)
NaCl.2H2O NaCl.2H2O
Tfm Ice + L+ V +L+ V Tfm Ice + L+ V +L+ V
-20.8oC -20.8oC
-25 -25
-50 -50
0 10 20 30 0 10 20 30
Weight % NaCl Weight % NaCl
20
Sample calculations from fluid inclusion observations
Temperature of ice melting directly determine the salinity
For diferent composition of inclusions (Te) we apply different 300 bars fluid pressure (hydrostatic)
phase diagrame to calculate salinities Fluid T = Rock T
Average fluid density = 1.0 g/cc
=> Depth = 3 km
Temperature = 300°C
=> Thermal gradient = 100°C /km
21
Ternary water salt systems
22