Influence of a low-viscosity layer between rigid inclusion and viscous matrix on inclusion rotation and matrix flow: A numerical study

We have used 2-D finite element modelling to investigate the influence of a permanent low-viscosity layer between matrix and inclusion on matrix flow and inclusion rotation under viscous simple shear flow. Rigid inclusions of different shape (circle, square, ellipse, lozenge, rectangle and skewed rectangles) and aspect ratio (R) were used. The calculated matrix flow pattern is neither bow tie nor eye-shaped. It is a new flow pattern that we call cat eyes-shaped, which is characterized by: (i) straight streamlines that slightly bend inwards at the inclusion's crests; (ii) elongate eye-shaped streamlines on each side of the inclusion; (iii) stagnation points in the centre of the eyes; (iv) absence of closed streamlines surrounding the inclusion; (v) changes in flow configuration with inclusion orientation; the lines of flow reversal bend and tilt, closed streamline circuits may disappear, and streamlines may bend outwards at the inclusion's crests.Concerning inclusion rotation, the numerical results show that: (i) a low-viscosity layer (LVL) makes inclusions with R = 1 rotate synthetically, but the rotation rate depends upon shape (circle or square) and orientation. Therefore, shape matters in the slipping mode. (ii) All studied shapes with R > 1 rotate antithetically when starting with the greatest principal axis (e₁) parallel to the shear direction ( = 0°); (iii) rotation is limited because there is a stable equilibrium orientation (_se) for all studied shapes with R > 1. (iii) There is also an unstable equilibrium orientation (_ue), and both _se and _ue depend upon inclusion's R and shape.The present numerical results closely agree with previous results of analogue experiments with a permanent low viscosity interface. Only minor deviations related with small shape differences were detected.     

Behaviour of an isolated rimmed elliptical inclusion in 2D slow incompressible viscous flow

For 2D linear viscous flow, it is shown that the rates of rotation and stretch of an isolated elliptical inclusion with a coaxial elliptical rim are fully determined by two corresponding scalar values. For power-law viscosity, effective viscosity ratios of the inclusion and rim to the matrix depend on orientation and the system is more complex but, in practice, the simplification with two scalar values still provides a good approximation. Finite-element modelling (FEM) is used to determine the two characteristic values across a wide parameter space for the linear viscous case, with a viscosity ratio (relative to the matrix) of the inclusion from 10⁶ to 1, of the rim from 10⁻⁶ to 1, axial ratios from 1.00025 to 20, and rim thicknesses relative to the inclusion axes of 5 to 20%. Results are presented in a multi-dimensional data table, allowing continuous interpolation over the investigated parameter range. Based on these instantaneous rates, the shape fabric of a population of inclusions is forward modelled using an initial value Ordinary Differential Equation (ODE) approach, with the simplifying but unrealistic assumption that the rim remains elliptical in shape and coaxial with respect to the inclusion. However, comparison with accurate large strain numerical experiments demonstrates that this simplified model gives qualitatively robust predictions and, for a range of investigated examples, also remarkably good quantitative estimates for shear strains up to at least γ = 5. A statistical approach, allowing random variation in the initial orientation, axial ratio and rim viscosity, can reproduce the characteristic shape preferred orientation (SPO) of natural porphyroclast populations. However, vorticity analysis based on the SPO or the interpreted stable orientation of inclusions is not practical. Varying parameters, such as inclusion and rim viscosity, rim thickness, and power law-exponents for non-linear viscosity, can reproduce the range of naturally observed behaviour (e.g., back-rotation, effectively stable orientations at back-rotated angles, a cut-off axial ratio separating rotating from stable inclusions) even for constant simple shear and these features are not uniquely characteristic of the vorticity of the background flow.     

Analogue modelling of the influence of shape and particle/matrix interface lubrication on the rotational behaviour of rigid particles in simple shear

The rotational behaviour of a rigid particle embedded in a linear viscous matrix undergoing cylindrical simple shear (Couette) flow was studied in 2D rock-analogue experiments. The influence of particle shape (elliptical vs. monoclinic), aspect ratio and the nature of the matrix/particle interface (lubricated vs. unlubricated) was investigated. Both matrix (PDMS) and lubricant (liquid soap) were linear viscous, with a viscosity ratio of ca. 10⁴. Without lubricant, the rotational behaviour of all particles closely approximates the Jeffery theory. Lubricated monoclinic particles with the long diagonal initially parallel to the shear direction show back rotation and approach a stable position. Lubricated elliptical particles initially parallel to the shear direction also show back rotation but only transiently stabilize. Weak planar zones in the matrix adjacent to unlubricated elliptical particles do not induce backward rotation. In general for elliptical particles, rotation rate as a function of orientation depends on axial ratio and thickness of the lubricant mantle. For thick mantles (initially >10% of the volume of the particle), rotation rates are faster than Jeffery theory. For very thin mantles they are markedly slower compared with thick mantles, particularly when the long axis is nearly parallel to the shear direction. Rotation rates are never strictly zero, so true stabilization does not occur. However, for more elongate particles (axial RATIO=6) rotation rates are so slow that a very strong shape preferred orientation would develop in a lubricated elliptical particle population. In experiments, the volume of lubricant is constant and the thickness adjacent to the long side of the particle progressively decreases with increasing strain. In natural examples of porphyroclast systems, the weak mantle continually develops by recrystallization and/or cataclasis of the rigid clast core and a steady state between production and thinning could be attained, potentially leading to true stabilization for particles with a high axial ratio.     

2D rotation of rigid inclusions in confined bulk simple shear flow: a numerical study

We have used incompressible Navier–Stokes in 2D finite element modelling to investigate rigid inclusion rotation under confined bulk simple shear flow. Confinement is defined as the ratio (S) between the channel width (H) and the inclusion's least axis (e₂)(S=H/e₂). The numerical results show that (i) inclusion rotation is strongly influenced by S and, when the confinement is effective, aspect ratio (R) and shape also play an important role. (ii) Back rotation is limited because inclusions reach a stable equilibrium orientation (_se). (iii) There is also an unstable equilibrium orientation (_ue), which defines an antithetic rotation field with _se, and both _se and _ue depend on S, and inclusion R and shape.     

2-D rotation behavior of a rigid ellipse in confined viscous simple shear: numerical experiments using FEM

The rotation behavior of rigid elliptical inclusions adherent to the viscous matrix in simple shear flow is investigated using a 2-D finite element numerical model. Several simulations were performed using different ratios (S) between shear zone width and inclusion's least principal axis. A computational strategy was devised to calculate pressure and viscous forces exerted on the inclusion and deduce its angular velocity. For large S values, results agree remarkably well with theoretical predictions, while for small S values results deviate significantly from theory but are in agreement with previous analogue experiments. The numerical model provided detailed and coherent information about the physical parameters involved in the process (e.g., pressure, strain rate and vorticity distributions within the model).     

Pulsating strains

Pulsating strains are modelled analytically using a continuum model with a well defined heterogeneous rheology. The mechanical continuum comprises single circular cylindrical inclusions of isotropic homogeneous viscosity hosted in a less viscous matrix. The analysis concentrates on the behaviour of the competent inclusion and demonstrates that the progressive deformation of the inclusion is entirely controlled by its viscosity ratio to the host and the orientation of the far field stress. The rate of deformation is determined by the stress magnitude and the viscosity of the inclusion, but does not affect the deformation path followed. An extreme case is provided by pure shear deformations where, by definition, the shear rate vanishes from the deformation tensor. Vorticity will then be zero, and strains will develop coaxially without pulsation.     

Characterisation of deformation and flow mechanics around porphyroclasts in a calcite marble ultramylonite by means of EBSD analysis

This paper studies the flow heterogeneity around porphyroclasts associated with greenschist facies deformation of a calcite marble shear zone. Microstructural data from electron backscatter diffraction analyses (EBSD) are used to constrain the flow mechanics of this dominantly non-coaxial type of deformation. The microstructure of the undisturbed ultramylonite (grain-size range 5–100 μm, mean 40 μm) is interpreted to represent steady-state (time-independent) flow conditions with flow planes parallel to the shear zone boundary. Single calcite porphyroclasts (grain-size 1–3 mm) caused flow perturbation in the fine-grained marble ultramylonite. It is the shape, in particular, of these rigid porphyroclasts that controls their rotational behaviour during deformation and, therefore, the development of specific flow fabrics. The flow planes around elongated-rhomboidal, stable porphyroclasts change the orientation to become roughly parallel to the porphyroclast margin, whereas the geometry of flow planes around nearly equant, rotating porphyroclasts describes a δ-type flow pattern. We infer that to some extent decoupling at the clast–matrix interface has occurred to guarantee a stable orientation of elongated porphyroclasts, but was not sufficient to reduce the rotation rate of equant clasts to zero. According to the flow deflection, the general crystallographic preferred orientation (CPO) with its single c-axis maximum perpendicular to the flow plane is rotated about an axis which is (sub)parallel to the kinematic rotation axis of the shear zone. Ultramylonite microstructures, CPOs and misorientation data are best explained by the dual operation of grain-size-insensitive (dislocation creep with recovery and recrystallization) and grain-size-sensitive (diffusion creep) mechanisms. The limited grain-size reduction around porphyroclasts suggests that the grain-size-insensitive mechanisms controlled rheology.     

Rotated staurolite porphyroblasts in the Littleton Schist at Bolton, Connecticut, USA

Staurolite porphyroblasts, 1.5–8cm in length and 0.3–2cm in width, in the Littleton Schist at Bolton, Connecticut, contain curved quartz inclusion trails which document synkinematic rotations of at least 135°. The orientations of long axes of these staurolite crystals define a weak preferred orientation in a plane approximately parallel to the external foliation. Serial sections of four differently orientated crystals and U-stage measurements of the orientations of their inclusion trails demonstrate that the inflection hinge line and the statistical 'symmetry axis' characterizing the foliation within a porphyroblast are unrelated to the orientations of external crenulations and are, in all cases, parallel to the long axis of the porphyroblast. The cumulative rotation reflected in the curvature of the inclusion trails is a maximum in a c -axis section through the initial core of a crystal. The amount of rotation about the c -axis decreases linearly along the length of the crystal away from the nucleation site.
The sense and amount of rotation recorded by a porphyroblast is related to its orientation. A tightly constrained transition from clockwise to anticlockwise rotation defines a slip direction that coincides with the preferred orientation of the staurolite c -axes. The total rotation reflected by the inclusion trails increases as a function of the angle between the c -axes of the staurolite crystals and the slip direction.
Initially random staurolite porphyroblasts rotated during growth, as a consequence of laminar shear in the surrounding viscous matrix. This interpretation is quantitatively consistent with: the staurolite preferred orientation; its coincidence with the apparent slip direction; the correlation between both the sense and the amount of rotation and the orientation of the long axis of the porphyroblast; and the twisted conical shape of the family of surfaces defined by the inclusion trails.     

Controls on the geometry of tails around rigid circular inclusions: insights from analogue modelling in simple shear

We used analogue modelling to investigate the factors controlling tail geometry in porphyroclast systems. Results show that: (1) σ inclusions can develop in both slipping and non-slipping modes, but δ-inclusions only form in the latter. (2) σ inclusions develop when the mantle production rate is constant and the mantle is transected by the separatrix. δ inclusions form when the mantle is initially outside the separatrix and later comes inside this line. (3) In the slipping mode, the wedge-shaped tail of σ-inclusions always has a straight external surface parallel to the shear plane, whereas in the non-slipping mode the external surface is curved inwards (external embayments). (4) Together with earlier theoretical results, σ- and δ inclusions always show stair-stepping of tails when embedded in a viscous matrix under homogeneous simple shear deformation. (5) Maximum stair-stepping occurs in the slipping mode and is at least equal to the inclusion diameter. If our models bear significant similarity to nature, then (i) the straight or curved character of σ-inclusions could mean that they had, respectively, a slipping or non-slipping interface with the surrounding recrystallized matrix, and (ii) δ-inclusions may result from shear deformation under retrogressive metamorphic conditions in thrust systems.     

Quartz deformation mechanisms during Barrovian metamorphism: Implications from crystallographic orientation of different generations of quartz in pelites

The behaviour of quartz during metamorphism is studied based on two case studies from the Barrovian terrains of Sulitjelma in arctic Scandinavia and Loch Tay in the Central Highlands Dalradian of Scotland. Both terrains preserve evidence for metamorphism in pelites involving nucleation and growth of garnet at different times in the deformation history. Data are presented on the size, shape and crystallographic orientation of quartz preserved as inclusions in garnet and as grains in the surrounding matrix. While quartz-grains remain small and dispersed between mica grains, deformation appears to be dominated by grain-boundary sliding accommodated by dissolution–precipitation. At amphibolite facies, textural coarsening occurs by dissolution of small quartz grains and growth of larger quartz grains, coupled with segregation of quartz from mica. As a result, quartz deforms by dislocation creep, developing crystallographic preferred orientations (CPO) consistent with both coaxial and non-coaxial strain. Quartz CPOs with <0001> axes lying parallel to foliation and stretching direction are commonly developed, and best explained by mechanical rotation of inequant (detrital?) quartz grains. There is no evidence for selective entrapment of quartz inclusions in garnet on the basis of quartz crystallographic orientation.     

流体作用下石榴石溶解-沉淀蠕变过程——以红河-哀牢山剪切带内石榴夕线片麻岩研究为例

石榴夕线片麻岩是中、下地壳主要组成岩石之一,岩石内石榴石和夕线石的结晶学优选方位会显著影响地壳深部流变性质,因此探讨特征变质矿物的变形机制和主要受控因素对构造带深部演化过程有深远意义。本文选取红河-哀牢山韧性剪切带内石榴夕线片麻岩为研究对象,通过定向切片内显微构造、电子探针、X-ray成分扫描、电子背散射衍射(EBSD)和相平衡模拟综合研究,揭示出石榴石在溶解沉淀反应过程中具有明显的粒度敏感性,不同粒径石榴石表现出截然不同的长宽比、成分环带、包裹体排列方式和压力影发育情况。石榴石表面流体活动明显截切早期生长环带。EBSD分析揭示石榴石破碎颗粒以绕〈112〉轴机械旋转为主,溶解过程主要集中于颗粒表面和裂隙内高曲率位置。夕线石的EBSD结果表明基质内夕线石以绕〈010〉轴旋转为主,而流体作用明显区域夕线石则以(100)[001]滑移为主。岩石相平衡模拟限定岩石变质峰期P-T条件达高压麻粒岩相,退变过程中同剪切变形导致大量流体渗入而形成降温降压退变轨迹,由~9.5kbar、760℃演化至~6.0kbar、500~600℃,并在粗粒石榴石内保存早期进变质环带,剪切抬升过程中石榴石内普遍发育垂直剪切方向的裂隙,并在流体作用下进一步改造其形态。此研究揭示红河-哀牢山剪切带内除前人报道的石榴石高温韧性变形外,还存在大量中-上地壳层次同变质反应下的溶解-沉淀蠕变作用。因此,石榴石变质-变形的综合研究有助于揭示变质杂岩带挤压-剪切-伸展多阶段构造演化过程。     

Pulsating oblate and prolate three-dimensional strains

The progressive ductile deformation of competent spherical inclusions is modeled analytically. Results of this study may help to understand better the limitations connected to geological field methods using competent inclusions for strain analyses. Parameters studied and quantified here are the strain magnitude, the progressive change in inclusion shape, the orientation of the finite strain axes, the frequency of pulsation, and the coupling between the strain ellipticity and viscosity contrast. Competent inclusions develop pulsating apparent strains if the host material is subjected to a component of simple shear and provided time or strain rate is sufficient to complete the strain cycle. The disparity between the strain magnitude inferred from competent viscous inclusions and that undergone by the host rock, increases for larger viscosity between them. The pulsation of the inclusion may suggest zero strain after a strain cycle has been completed, even though strain in the host rock is extremely large. The inclusion will develop pulsating oblate strains if a shortening rate is superposed normal to the plane of pulsation. Conversely, pulsating prolate strains occur if an extension rate is superposed instead of shortening. Stretching lineations outlined by deformed competent inclusions within shear zones beneath collapsing nappe sheets may even point perpendicular to the direction of nappe transport. This finding offers an explanation for the occurrence of mutually perpendicular pebble elongations in nearby locations within the Bygdin conglomerate beneath the Jotun nappe, Norwegian Caledonides.     

Modelling the internal stress field in argillaceous rocks under humidification/desiccation

This paper deals with analytical and numerical modelling of the internal stress generated in argillaceous rocks during humidification/desiccation processes, which is an essential issue for damage study. This local stress field arises from two mechanisms: (i) complex interactions between free swelling/shrinking clay matrix and non‐strained inclusions of carbonate and quartz and (ii) a self‐restraint effect induced by the moisture gradient during the transient moisture exchange process. The inclusion–matrix interaction is investigated in different cases. Firstly, the analytical solution of the stress around a cylindrical inclusion embedded in an infinite swelling matrix is derived: The inclusion would suffer tension (compression) under humidification (desiccation), and the resulting cracking patterns are discussed. Then, the problem of two inclusions with different distances in an infinite swelling matrix is considered, and it is shown that the local stress around an inclusion will be perturbed and amplified by neighbouring inclusions. Finally, an inclusion outcropping at the free surface of a swelling matrix is modelled as to investigate the effect of free surface: The inclusion–matrix interface undergoes shear stresses of which the maximum is found at the free surface. In addition to the inclusion–matrix interaction, the self‐restraint effect is investigated: The induced stress is maximal at the beginning of humidification/desiccation processes and vanishes gradually with time. The quantity of the self‐restraint stress is strongly controlled by the hydric loading rate. Copyright © 2014 John Wiley & Sons, Ltd.     

Microlite orientation in obsidian flow measured by synchrotron X-ray diffraction

Clinopyroxene and plagioclase (andesine) microlites in an obsidian flow from Glass Mountain (NE California, USA) display strong alignment. Synchrotron X-ray diffraction, coupled with Rietveld analysis, was used to quantify crystallographic-preferred orientation (CPO). Clinopyroxene, with a rod-shaped morphology, shows a strong alignment of [001] in the flow direction and (010) aligned parallel to the inferred flow plane. Andesine, with a platy morphology, displays an alignment of (010) platelets in the flow plane. Some pole densities exceed 90 multiples of random distribution. Applying a model of rigid ellipsoidal inclusions in a viscous matrix, the local pure shear strains are between 2 and 3.     

Numerical modeling of the motion of rigid ellipsoidal objects in slow viscous flows: A new approach

A simple algorithm for modeling the rotation of rigid ellipsoidal objects in viscous flows based on Jeffery's (1922, Proceedings of the Royal Society of London A102, 161–179) theory is presented and is implemented in a fully graphic mathematics application Mathcad^® (http://www.mathsoft.com). The orientation of ellipsoidal objects is specified in terms of polar coordinate angles that can be easily converted to the trend and plunge angles of the three principal axes rather than the Euler angles. With the Mathcad worksheets presented in the supplementary data associated with this paper, modeling the rotation paths of individual rigid objects, the development of inclusion trail geometry within syn-kinematic porphyroblasts, and the development of preferred orientation and shape fabrics for a population of rigid objects becomes as easy a task as using a spreadsheet. The shape and preferred orientation fabrics for a population of rigid objects can be presented in both a three-dimensional form and a two-dimensional form, allowing easy comparison between field data and model predictions. The modeler can customize the type and format of the output to best fit the purpose of the investigation and to facilitate the comparison of model predictions with geological observations. Application examples are presented for various types of modeling involving rigid objects.     

Porphyroblast rotation: eppur si muove*?

In a number of recent papers, the theory has been postulated that porphyroblasts as a rule do not rotate with respect to geographical coordinates, and can be used to determine the original orientation of older foliations. Complex inclusion patterns in spiral garnets have even been used to advocate a new model of orogenesis, involving several alternating phases of horizontal shortening and extension. Critical assessment of the assumptions and data used to support the theory of irrotational porphyroblasts reveals numerous flaws. Millipede structures, used as proof for flow partitioning, can also form by other flow geometries. Evidence quoted to support irrotational behaviour of porphyroblasts is unsound. Porphyroblasts do occur in sets with a preferred orientation of the internal foliation trace, but these cannot be shown to represent original orientations. Microstructures which resemble truncation planes in spiral garnets are used as evidence that these structures developed by several phases of deformation and as proof for periodic extension and horizontal shortening in orogenesis. They can, however, also be explained by intermittent growth of a rotating porphyroblast during a single phase of deformation. Finally, porphyroblast sets in which orientation is a function of aspect ratio indicate that porphyroblast rotation with respect to kinematic axes does occur in at least some situations.     

Rheological and kinematical responses to flow of two-phase rocks

Numerical simulations have been performed to investigate the strain-dependent behaviour of rheological and kinematical responses to flow of two-phase rocks using the commercial finite-difference program FLAC2D. It was assumed that the two phases have Maxwell rheology. Plane strain and velocity boundary condition, which produces a simple shear deformation, were also assumed. Two types of geometries were considered: strong phase supported (SPS) and weak phase supported (WPS). We calculated strain-dependent variations of effective viscosity and partitioning of strain rate, vorticity and kinematic vorticity number during deformation in both SPS and WPS structure models.The results show that the strain-dependent behaviour is largely influenced by the geometry of the composite. SPS models show both strain hardening and strain softening during the simulations, with strain hardening preceding strain softening. A critical shear strain is necessary to begin the strain softening behaviour. Strain hardening and strain softening are accompanied by a reduction and an increase of the partition of strain rate into the weak phase, respectively. On the other hand, WPS models show only weak strain hardening and strain softening, being the strain-dependent behaviour close to a steady state flow. In addition, the following results are obtained on vorticity and kinematic vorticity number; (1) in both SPS and WPS models the partition of vorticity into weak phase increases with progressive shear strain, i.e. the strong phase becomes less rotational, (2) in SPS models weak inclusions changes from sub-simple shear to super-simple shear with progressive strain, whereas the strong matrix changes from super-simple shear to sub-simple shear, (3) in WPS models the strong inclusions with high viscosity contrasts are less rotational but can be in super-simple shear condition to high strains.The observed strain-dependent behaviours have been compared with previous proposed analytical models. The degree of agreement is variable. Balshin and Ryshkewitch–Duckworth models are only applicable to SPS models. Ji-generalized mixture rule model is applicable to both models.The results suggest that polyphase rocks with SPS structure during ductile shear deformation respond as strain softening materials, after an initial strain hardening stage that may drive to the strain localization into the material.     

Fluid entrapment and reequilibration during subduction and exhumation: A case study from the high-grade Nestos shear zone, Central Rhodope, Greece

A fluid inclusion study on metamorphic minerals of successive growth stages was performed on highly deformed paragneisses from the Nestos Shear Zone at Xanthi (Central Rhodope), in which microdiamonds provide unequivocal evidence for ultrahigh-pressure (UHP) metamorphism. The correlation of fluid inclusion density isochores and fluid inclusion reequilibration textures with geothermobarometric data and the relative chronology of micro- and macro-scale deformation stages allow a better understanding of both the fluid and metamorphic evolution along the P–T–d path. Textural evidence for subduction towards the NE is recorded by the orientation of intragranular NE-oriented fluid inclusion planes and the presence of single, annular fluid inclusion decrepitation textures. These textures occur within quartz “foam” structures enclosed in an earlier generation of garnets with prolate geometries and rarely within recrystallized matrix quartz, and reequilibrated both in composition and density during later stages of exhumation. No fluid inclusions pertaining to the postulated ultrahigh-pressure stage for microdiamond-bearing garnet–kyanite–gneisses have yet been found. The prolate shape of garnets developed during the earliest stages of exhumation that is recorded structurally by (L  S) tectonites, which subsequently accommodated progressive ductile SW shearing and folding up to shallow crustal levels. The majority of matrix kyanite and a later generation of garnet were formed during SW-directed shear under plane-strain conditions. Fluid inclusions entrapped in quartz during this stage of deformation underwent density loss and transformed to almost pure CO₂ inclusions by preferential loss of H₂O. Those inclusions armoured within garnet retained their primary 3-phase H₂O–CO₂ compositions. Reequilibration of fluid inclusions in quartz aggregates is most likely the result of recrystallization along with stress-induced, preferential H₂O leakage along dislocations and planar lattice defects which results in the predominance of CO₂ inclusions with supercritical densities. Carbonic fluid inclusions from adjacent kyanite–corundum-bearing pegmatoids and, the presence of shear-plane-parallel fluid inclusion planes within late quartz boudin structures consisting of pure CO₂-fluid inclusions with negative crystal shapes, bear witness of the latest stage of deformation by NE-directed extensional shear.This study shows that the textures of early fluid inclusions that formed already during the prograde metamorphic path can be preserved and used to derive information about the kinematics of subduction that is difficult to obtain from other sources. The textures of early inclusions, together with later generations of unaltered primary and secondary inclusions in metamorphic index minerals that can be linked to specific deformation stages and even P–T conditions, are a welcome supplement for the reconstruction of a rather detailed P–T–d path.     

Petroleum inclusions in sedimentary basins: systematics, analytical methods and applications

Fluid dynamics in sedimentary basins is of tremendous interest, both from a scientific and an economic point of view. Integration between fluid inclusion and present-day fluid data provides the time aspect necessary for reconstruction of fluid flow paths, and can be used for mapping fluid dynamics both on a regional basin scale or on the more local scale of petroleum reservoirs. This paper presents a review of analytical and modelling methods for petroleum in fluid inclusions. Essentially, four types of data for petroleum inclusions can be established through analysis and modelling: (1) textures (which give indirect information of the time aspect), (2) fluid composition, (3) fluid properties and (4) pressure–temperature of trapping. During the last decade, development of analytical methods for determination of inclusion fluid compositions has taken place. Traditional correlation studies for characterisation of maturity and source facies may now therefore include fluid inclusion data. The development has also been directed towards an improved understanding of the physical properties of petroleum inclusions. Although these methods are in their early stages concerning precision and accuracy, data necessary for fluid flow modelling, such as fluid densities, viscosities, pressure and temperature, can be estimated.