Review and classification of structural controls on fluid flow during regional metamorphism

Mechanisms for kilometre-scale, open-system fluid flow during regional metamorphism remain problematic. Debate also continues over the degree of fluid flow channellization during regional metamorphism, and the mechanisms for pervasive fluid flow at depth. The requirements for pervasive long-distance fluid flow are an interconnected porosity and a large regional gradient in fluid pressure and hydraulic head (thermally or structurally controlled) that dominates over local perturbations in hydraulic head due to deformation. In contrast, dynamic or transient porosity interconnection and fluid flow accompanying deformation of heterogeneous rock suites should result in moderately to strongly channellized flow at a range of scales, of which there are many examples in the literature. Classification of fluid flow types based on scale and degree of equilibration between fluid and rock, wallrock permeability, and mode of fluid transport contributes to an understanding of key factors that control fluid flow. Closed-system fluid behaviour, with restricted fluid flow in microcracks or cracks and limited fluid–rock interaction, occurs over a range of strains and crustal depths, but requires low permeabilities and/or small fluid fluxes. Long-distance, open-system fluid flow in channels is favoured in heterogeneous rocks at high strains, moderate (but variable) permeabilities, and moderate to high fluid fluxes. Long-distance, broad, pervasive fluid flow during regional metamorphism requires that the rocks are not accumulating high strains and have high permeabilities, low permeability contrasts, and high fluid fluxes. The ideal situation for such fluid flow is in situations where the rocks are undergoing stress relaxation immediately after a major deformation phase. In the mid-crust, fairly specific conditions are thus required for pervasive fluid flow. During active orogenesis, structurally controlled fluid flow (with focused open systems surrounding regions of closed-system behaviour) predominates in most, but not all, regional metamorphic situations, at a range of scales.     

Dissolution, solution transfer, diffusion versus fluid flow and volume loss during deformation/metamorphism

ABSTRACT Dissolution and solution transfer during deformation/metamorphism are controlled by the partitioning of deformation into progressive shearing and shortening components. Progressive shearing is readily accommodated by slip on the planar crystal structure of phyllosilicates and graphite without accumulating dislocation density gradients across grain boundaries. Progressive shortening is accommodated by the cores of most other minerals (including sulphides). These minerals develop strain, and hence dislocation density gradients, on their rims due to progressive shearing along grain boundaries. These gradients are particularly large when the mineral abuts phyllosilicate or graphite. The resulting chemical potential gradients between the core and rim drive dissolution, causing removal of the highly strained grain margins. Removal of dissolved material by solution transfer is aided by the geometry of shearing of phyllosilicates and graphite around other grains in an active anastomosing foliation. Interlayers and interfaces on boundaries lying at a low angle to the direction of shearing, and oriented relative to the sense of shear such that they can open, gape by small amounts. Water present in these interlayer spaces becomes destructured, considerably enhancing diffusion rates along the foliation. Penetrative volume loss, especially in deforming/metamorphosing pelitic rocks, is large at all metamorphic grades, increasing and becoming more penetrative with depth to at least the transition into granulite and eclogite facies. Transference of material by fluid flow from deep to high levels in the earth's crust is precluded because thousands to tens of thousands of rock volumes of fluid are required, necessitating continual recirculation of fluid from shallow to deep crustal levels in one large or several small sets of cells, unless some extremely large-scale form of fluid channelling is possible. Reassessment of diffusion mechanisms, and hence rates, during deformation and pervasive foliation generation in large volumes of rock where fluid channeling cannot provide enough fluid, indicates that diffusion can proceed with sufficient rapidity that massive recirculation of fluid is no longer required. The amount of fluid can be reduced sufficiently to allow large volume losses by a one-way flow of fluid to the earth's surface, in deforming/metamorphosing environments where the fluid pressure equals or exceeds the hydrostatic pressure. Deformation partitioning-controlled dissolution progressively changes the bulk chemistry of a rock containing phyllosilicates or graphite during deformation/metamorphism because matrix minerals, other than phyllosilicates and graphite, are preferentially removed. The large size of porphyroblasts, if present, tends to preserve them from dissolution. Hence, the bulk chemistry operative during subsequent porphyroblast growth can have changed considerably from that operative when the first porphyroblasts grew, in rocks in which bedding is still well preserved.     

The effect of heterogeneous stress and strain on metamorphic fluid flow,Mary Kathleen,Australia, and a model for large-scale fluid circulation *

In the Proterozoic Mary Kathleen Fold Belt, northern Australia, infiltration of large volumes of externally derived fluid occurred synchronously with regional amphibolite-facies metamorphism and deformation. This paper develops a model of structurally controlled fluid migration by comparing the distribution of fossil fluid pathways with the inferred stress and strain patterns during the deformation. Intense fluid flow was localized within strong, relatively brittle meta-intrusive bodies, and in discrete, veined, brecciated and altered zones around their margins. In metasediments folded in a ductile manner outside these areas, fluid infiltration was negligible. The direct correlation between structural styles and the magnitude of veining and metasomatism suggests control of permeability enhancement, and hence fluid flow, by deformation. Finite difference modelling of a strong body in a weaker matrix has been used to evaluate the variation of stresses during the deformation, from which it is clear that stress and strain heterogeneities have systematically influenced the development and maintenance of metamorphic fluid pathways. Particular regions in which mean stress may be significantly lower than the average lithostatic pressures include the ‘strain shadow’zones adjacent to the strong bodies, other dilatant zones around the bodies, and the bodies themselves. This geometry is favourable not only for localized brittle deformation under amphilobite facies conditions, but also for focused fluid flow in the low mean stress regions, as evidenced by the abundance of veins. Fluid access through these metamorphic aquifers occurred during tensile failure episodes, with particularly large dilations and decimetre-scale veining in areas of strain incompatibility. It appears likely that fluid circulated many times through the Fold Belt, with flow concentrated in the metamorphic aquifers. A model is developed that explains both the structurally focused fluid flow and the postulated multi-pass recirculation by dilatancy pumping, the ‘pump engines’comprising the low mean stress zones.     

区域变质作用中的流体

区域变质条件下流体的流动有 4种标志 :( 1)细脉 ;( 2 )岩石学 ;( 3 )稳定同位素 ;( 4 )常量元素的交代作用。不同级别的区域变质作用中 ,流体影响着岩石的变质反应和变形 ;在高级变质的情况下甚至有熔体出现。在超高压变质条件下 ,流体量比地壳范围区域变质要少得多 ,从大别山超高压变质带的资料可知 ,流体的演化有明显的阶段性 ,局部曾发现熔融包裹体。水流体的介入 ,引起岩石的退变质和元素地球化学变异 ,是超高压变质岩抬升、进入中下地壳的产物。新近的实验岩石学成果说明 ,多硅白云母、角闪石等含羟基的矿物 ,在俯冲达 10 0km以下依然稳定 ,而一些花岗岩体系在超高压的条件下产生的超临界流体 ,乃是花岗岩、片麻岩只能部分保留超高压矿物组合的原因。     

Fluid flow and mass transport at the Valentine wollastonite deposit, Adirondack Mountains, New York State

The Valentine wollastonite skarn in the north-west Adirondack Mountains, New York, is a seven million ton deposit which resulted from channellized infiltration of H₂O-rich, silica-bearing fluids. The wollastonite formed by reaction of these fluids with non-siliceous calcite marble. The skarn formed at the contact of the syenitic Diana Complex and was subsequently overprinted by Grenville-age granulite facies metamorphism and retrograde hydrothermal alteration during uplift. Calcite marbles adjacent to the deposit have generally high δ¹⁸O values (c. 21‰), typical of Grenville marbles which have not exchanged extensively with externally derived fluids. Carbon isotopic fractiona-tions between coexisting calcite and graphite in the marbles indicate equilibration at 675d? C, consistent with the conditions of regional metamorphism. Oxygen isotopic ratios from wollastonite skarn are lower than in the marbles and show a 14‰ variation (-1‰ to 13‰). Some isotopic heterogeneity is preserved from skarn formation, and some represents localized exchange with low-δ¹⁸O retrograde fluids. Detailed millimetre- to centimetre-scale isotopic profiles taken across skarn/marble contacts reveal steep δ¹⁸O gradients in the skarn, with values increasing towards the marble. The gradients reflect isotopic evolution of the fluid as it reacted with high δ¹⁸O calcite to form wollastonite. Calcite in the marble preserves high δ¹⁸O values to within <5 mm of the skarn contact. The preservation of high δ¹⁸O values in marbles at skarn contacts and the disequilibrium fractionation between wollastonite skarn and calcite marble across these contacts indicate that the marbles were not infiltrated with significant quantities of the fluid. Thus, the marbles were relatively impermeable during both the skarn formation and retrograde alteration. Skarn formation may have been episodic and fluid flow was either chaotic or dominantly parallel to lithological contacts. Although these steep isotope gradients resemble fluid infiltration fronts, they actually represent the sides of the major flow system. Because chromatographic infiltration models of mass transport require the assumption of pervasive fluid flow through a permeable rock, such models are not applicable to this hydrothermal system and, by extension, to many other metamorphic systems where low-permeability rocks restrict fluid migration pathways. Minimum time-integrated fluid fluxes have been calculated at the Valentine deposit using oxygen isotopic mass balance, reaction progress of fluid buffering reactions, and silica mass balance. All three approaches show that large volumes of fluid were necessary to produce the skarn, but silica mass balance calculations yield the largest minimum flux and are hence the most realistic.     

Up-temperature flow of surface-derived fluids in the mid-crust: the role of pre-orogenic burial of hydrated fault rocks

The Walter‐Outalpa shear zone in the southern Curnamona Province of NE South Australia is an example of a shear zone that has undergone intensely focused fluid flow and alteration at mid‐crustal depths. Results from this study have demonstrated that the intense deformation and ductile shear zone reactivation, at amphibolite facies conditions of 534 ± 20 °C and 500 ± 82 MPa, that overprint the Proterozoic Willyama Supergroup occurred during the Delamerian Orogeny (c. 500 Ma) (EPMA monazite ages of 501 ± 16 and 491 ± 19 Ma). This is in contrast to the general belief that the majority of basement deformation and alteration in the southern Curnamona Province occurred during the waning stages of the Olarian Orogeny (c. 1610–1580 Ma). These shear zones contain hydrous mineral assemblages that cut wall rocks that have experienced amphibolite facies metamorphism during the Olarian Orogeny. The shear zone rock volumes have much lower δ¹⁸O values (as low as 1‰) than their unsheared counterparts (7–9‰), and calculated fluid δ¹⁸O values (5–8‰) consistent with a surface‐derived fluid source. Hydrous minerals show a decrease in δD_(H2O) from ?14 to ?22‰, for minerals outside the shear zones, to ?28 to ?40‰, for minerals within the shear zones consistent with a contribution from a meteoric source. It is unclear how near‐surface fluids initially under hydrostatic pressure penetrate into the middle crust where fluid pressures approach lithostatic, and where fluid flow is expected to be dominantly upward because of pressure gradients. We propose a mechanism whereby faulting during basin formation associated with the Adelaidean Rift Complex (c. 700 Ma) created broad hydrous zones containing mineral assemblages in equilibrium with surface waters. These panels of fault rock were subsequently buried to depths where the onset of metamorphism begins to dehydrate the fault rock volumes evolving a low δ¹⁸O fluid that is channelled through shear zones related to Delamerian Orogenic activity.     

The stable isotope signature of kilometre-scale fracturedominated metamorphic fluid pathways, Mary Kathleen, Australia

Abstract Large calcite veins and pods in the Proterozoic Corella Formation of the Mount Isa Inlier provide evidence for kilometre-scale fluid transport during amphibolite facies metamorphism. These 10- to 100-m-scale podiform veins and their surrounding alteration zones have similar oxygen and carbon isotopic ratios throughout the 200 × 10-km Mary Kathleen Fold Belt, despite the isotopic heterogeneity of the surrounding wallrocks. The fluids that formed the pods and veins were not in isotopic equilibrium with the immediately adjacent rocks. The pods have δ¹³C_calcite values of –2 to –7% and δ¹⁸O_calcite values of 10.5 to 12.5%. Away from the pods, metadolerite wallrocks have δ¹⁸O_whole-rock values of 3.5 to 7%. and unaltered banded calc-silicate and marble wallrocks have δ¹³C_calcite of –1.6 to –0.6%, and δ¹⁸O_calcite of 18 to 21%. In the alteration zones adjacent to the pods, the δ¹⁸O values of both metadolerite and calc-silicate rocks approach those of the pods. Large calcite pods hosted entirely in calc-silicates show little difference in isotopic composition from pods hosted entirely in metadolerite. Thus, 100- to 500-m-scale isotopic exchange with the surrounding metadolerites and calc-silicates does not explain the observation that the δ¹⁸O values of the pods are intermediate between these two rock types. Pods hosted in felsic metavolcanics and metasiltstones are also isotopically indistinguishable from those hosted in the dominant metadolerites and calc-silicates. These data suggest the veins are the product of infiltration of isotopically homogeneous fluids that were not derived from within the Corella Formation at the presently exposed crustal level, although some of the spread in the data may be due to a relatively small contribution from devolatilization reactions in the calc-silicates, or thermal fluctuations attending deformation and metamorphism. The overall L-shaped trend of the data on plots of δ¹³C vs. δ¹⁸O is most consistent with mixing of large volumes of externally derived fluids with small volumes of locally derived fluid produced by devolatilization of calc-silicate rocks. Localization of the vein systems in dilatant sites around metadolerite/calc-silicate boundaries indicates a strong structural control on fluid flow, and the stable isotope data suggest fluid migration must have occurred at scales greater than at least 1 km. The ultimate source for the external fluid is uncertain, but is probably fluid released from crystallizing melts derived from the lower crust or upper mantle. Intrusion of magmas below the exposed crustal level would also explain the high geothermal gradient calculated for the regional metamorphism.     

变质流体研究新进展

变质流体是变质过程的主要动力学因素之一。目前变质流体研究主要集中在下部地壳麻粒岩相变质流体，俯冲带高压－超高压变质流体和接触变质流体等方面。研究的主要问题是流体流动机制和元素迁移，流体－岩石相互作用和流体来源。下部地壳麻粒岩相变质流体以ＣＯ２为主，具有较低的ａＨ２Ｏ。δ１３Ｃ研究表明大约２／３ＣＯ２是深成的。富ＣＯ２流体流动是紫苏花岗岩形成和热扰动的原因之一，也是麻粒岩形成和大离子亲石元素亏损的主要因素。俯冲带是高压、超高压变质作用发生和流体活动最活跃的场所。流体富含Ｈ２Ｏ、ＣＨ４和ＣＯ２，可以诱导部分熔融反应和岛弧岩浆作用。高压变质条件下的矿物稳定性也与流体有关。同位素研究表明，在超高压变质期间没有化学上完全相同的流体大规模循环。流体－熔体系统模式能更有效地解释下插板片的元素再循环。接触变质流体研究主要集中在含有易于发生流体－岩石反应的不纯碳酸盐岩地区。硅灰石带中流体／岩石比率高达４０∶１，表明接触变质岩石中有大量流体存在。接触变质过程流体成分有较大差异，主要取决于流体来源、原岩性质和侵入体特征。流体流动和循环模式受控于构造变形，岩浆作用和变质过程的动力学条件及流体成分。     

Fracturing,fluid flow and shear zone development: Relationships between chemical and mechanical processes in Proterozoic mafic dykes from southwestern Montana,USA

Fluids can play an important role in the localization of deformation in the deep crust, yet the specific mechanisms active during the complex interactions between metasomatism, metamorphism and deformation remain elusive. Precambrian metagabbronorite dykes in southwest Montana contain fractures filled with Hbl±Grt and discrete cm‐scale shear zones with well‐preserved strain gradients. This system offers an ideal opportunity to constrain the chemical and mechanical processes that facilitated strain localization. An early M₁ assemblage of Grt₁+Cpx₁+Pl₁+Qz developed at conditions of 0.51–0.85 GPa and 500–700°C and is preserved largely as a static replacement of relict igneous phases (Opx, Pgt, Pl) in coronitic textures. An M₂ assemblage characterized by Grt₂+Pl₂±Cpx₂+Hbl+Scp+Qz developed at 0.86–1.00 GPa and 660–730°C coincided with fluid flow and deformation associated with shear zone development. Microstructural observations in marginal protomylonite/mylonite and laminated ultramylonite suggest a shear zone evolution that involved (1) nucleation from pre‐existing fractures that were sites for major fluid infiltration, (2) initial widening coincident with grain‐size reduction by microfracturing, dislocation creep, and synkinematic metamorphic reaction by solution transfer, and (3) a switch in the dominant deformation mechanisms active in the ultramylonite from grain‐size insensitive mechanisms to grain‐size sensitive granular flow accommodated by fluid‐assisted diffusion. Throughout this evolution, the effective bulk compositions of the rock volumes responding to metamorphism changed through a combination of mechanical and metasomatic processes.     

Contact mechanism of a rock fracture subjected to normal loading and its impact on fast closure behavior during initial stage of fluid flow experiment

Fast closure of rock fractures has been commonly observed in the initial stage of fluid flow experiments at environmental temperatures under low or moderate normal stresses. To fully understand the mechanisms that drive this fast closure, the evolution of local stresses acting on contacting asperities on the fracture surfaces prior to fluid flow tests needs to be evaluated. In this study, we modeled numerically the asperity deformation and failure processes during initial normal loading, by adopting both elastic and elastic–plastic deformation models for the asperities on a real rock fracture with measured surface topography data, and estimated their impact on initial conditions for fluid flow test performed under laboratory conditions. Compared with the previous models that simulate the normal contact of a fracture as the approach of two rigid surfaces without deformations, our models of deformable asperities yielded smaller contact areas and higher stresses on contacting asperities at a given normal stress or normal displacement. The results show that the calculated local stresses were concentrated on the contacts of a few major asperities, resulting in crushing of asperity tips. With these higher contact stresses, however, the predicted closure rates by pressure solution are still several orders of magnitude lower than that of the experimental measurements at the initial stage of fluid flow test. This indicates that single pressure solution may not likely to be the principal compaction mechanism for this fast closure, and that the damages on contacting asperities that occur during the initial normal loading stage may play an important role. Copyright © 2015 John Wiley & Sons, Ltd.     

Contrasting mechanisms of fluid flow through adjacent stratigraphic units during regional metamorphism,south-central Maine,USA

The flow pattern of reactive metamorphic fluid through six outcrops of micaceous, carbonate-bearing sandstones from the Vassalboro Formation was determined by calculating and mapping fluid-rock ratios for numerous samples within each outcrop. The ratio of maximum to minimum measured fluid/rock varied by factors of only 1.3-22.9 in each outcrop. Fluid flow was pervasive at metamorphic grades ranging from the biotite through the sillimanite zones. Average fluid-rock ratio for the outcrops increases with increasing grade of metamorphism from 0.4 in the biotite zone to 1.4 in the sillimanite zone.The flow pattern of reactive fluid through impure sandstones of the Vassalboro Formation was different at low and medium grades from fluid flow through the limestone member of the adjacent Waterville Formation. In the biotite and garnet zones, fluid flow through the Waterville Formation was channelized with channelways corresponding to individual lithologic layers that acted as metamorphic aquifers. Fluid-rock ratios recorded by the aquifers are greater than those recorded by the intervening beds by factors of up to 50–60. At the highest grades of metamorphism (sillimanite zone), however, flow through the Waterville Formation was as pervasive as through the Vassalboro Formation.The Waterville and Vassalboro Formations experienced the same metamorphic event. The difference in pattern of fluid flow through the two formations therefore reflects the important control that lithology exerts on the permeability of rocks during metamorphism. Micaceous, carbonate-bearing sandstones evidently were more permeable than argillaceous carbonate rocks. The greater permeability of the sandstones may result from a greater concentration of grain boundaries between unlike minerals in the rocks.     

造山过程中的流体—岩石相互作用和质量传输的评述

在汇聚板块边缘,滑脱面的滑动和沉积岩系的逆冲堆叠,导致孔隙水被压实排出,流体润滑了脱面,从而引起增生楔的生长。流体对流体制是地壳深熔的先驱事件。碱性花岗岩在次固相下有丰富的岩浆水的出溶,并且促进了碱性长石的微组构重组织。上部地壳浅表流体的循环主要受岩浆侵位驱动。韧性剪切带Ｔｉ、Ｆｅ、Ｍｇ残留富集,Ｓｉ、Ｃａ、Ｓｒ带出,流体不混溶和相分离是Ａｕ沉淀的重要机制。断层带尤其深部断层带具有高的流体／岩石比     

剪切带研究的某些进展

简要介绍了一般剪切带的成分变异、体积变化、流体流动的研究进展、着重讨论了榴辉岩相剪切带流体流动与变形变质作用的关系，认为剪切带应变不均匀性、变质作用不均一性和流体存在状态及流体超压是榴辉岩形成和保存的3个主要制约因素。     

Hydrothermal alteration, fluid flow and volume change in shear zones: the layered mafic–ultramafic Kettara intrusion (Jebilet Massif, Variscan belt, Morocco)

During emplacement and cooling, the layered mafic–ultramafic Kettara intrusion (Jebilet, Morocco) underwent coeval effects of deformation and pervasive fluid infiltration at the scale of the intrusion. In the zones not affected by deformation, primary minerals (olivine, plagioclase, clinopyroxene) were partially or totally altered into Ca‐amphibole, Mg‐chlorite and CaAl‐silicates. In the zones of active deformation (centimetre‐scale shear zones), focused fluid flow transformed the metacumulates (peridotites and leucogabbros) into ultramylonites where insoluble primary minerals (ilmenite, spinel and apatite) persist in a Ca‐amphibole‐rich matrix. Mass‐balance calculations indicate that shearing was accompanied by up to 200% volume gain; the ultramylonites being enriched in Si, Ca, Mg, and Fe, and depleted in Na and K. The gains in Ca and Mg and losses in Na and K are consistent with fluid flow in the direction of increasing temperature. When the intrusion had cooled to temperatures prevailing in the country rock (lower greenschist facies), deformation was still active along the shear zones. Intense intragranular fracturing in the shear zone walls and subsequent fluid infiltration allowed shear zones to thicken to metre‐scale shear zones with time. The inner parts of the shear zones were transformed into chlorite‐rich ultramylonites. In the shear zone walls, muscovite crystallized at the expense of Ca–Al silicates, while calcite and quartz were deposited in ‘en echelon’ veins. Mass‐balance calculations indicate that formation of the chlorite‐rich shear zones was accompanied by up to 60% volume loss near the centre of the shear zones; the ultramylonites being enriched in Fe and depleted in Si, Ca, Mg, Na and K while the shear zones walls are enriched in K and depleted in Ca and Si. The alteration observed in, and adjacent to the chlorite shear zones is consistent with an upward migrating regional fluid which flows laterally into the shear zone walls. Isotopic (Sr, O) signatures inferred for the fluid indicate it was deeply equilibrated with host lithologies.     

Fluid flow,metasomatism and amphibole deformation in an imbricated ophiolite,North Cascades,Washington

Metasomatic tremolite-rich mylonites are widespread in imbricate thrust slices of ultramafic rocks of the ophiolitic Ingalls Complex in Washington State. Protoliths for these amphibolite-facies mylonites were peridotite and serpentinite. Abundant syntectonic tremolite veins in the ultramafites record narrowly channelized flow of infiltrating fluids, whereas metasomatic mylonite zones record more pervasive flow. Fluids were probably released mainly by prograde devolatization reactions within serpentinite and mafic ophiolitic rocks that experienced earlier hydrothermal metamorphism.Olivine apparently deformed by dislocation creep in the mylonites. In the tremolite-rich rocks, locally preserved amphibole porphyroclasts deformed mainly by microfracturing. Acicular tremolites, which dominate the mylonites, form syntectonic overgrowths on porphyroclasts and probably record diffusive mass transfer which may have accompanied cataclasis. Acicular tremolites subsequently were folded and define both post-crystalline crenulations and polygonal arcs.Fluid flow, deformation and metamorphism were apparently complexly interrelated in the imbricate zone. Thrusts juxtaposed contrasting rock types that were sources and sinks for fluids, and shear zones focused fluid flow. Metamorphism probably facilitated deformation through the release of fluids during dehydration reactions. High fluid pressure may have led to hydraulic fracturing and may have controlled strain softening in the tremolitic mylonite zones as it favored microcracking and diffusive mass transfer over dislocation creep. Infiltrating metasomatic fluids probably play an important role in the evolution of shear zones in many ultramafic bodies during medium-grade metamorphism.     

Compaction-driven fluid flow in viscoelastic rock

Abstract

Compaction driven fluid flow is inherently unstable such that an obstruction to upward fluid flow (i.e. a shock) may induce fluid-filled waves of porosity, propagated by dilational deformation due to an effective pressure gradient within the wave. Viscous porosity waves have attracted attention as a mechanism for melt transport, but are also a mechanism for both the transport and trapping of fluids released by diagenetic and metamorphic reactions. We introduce a mathematical formulation applicable to compaction driven flow for the entire range of rheological behaviors realized in the lithosphere. We then examine three first-order factors that influence the character of fluid flow: (1) thermally activated creep, (2) dependence of bulk viscosity on porosity, and (3) fluid flow in the limit of zero initial connected porosity. For normal geothermal gradients, thermally activated creep stabilizes horizontal waves, a geometry that was thought to be unstable on the basis of constant viscosity models. Implications of this stabilization are that: (1) the vertical length scale for compaction driven flow is generally constrained by the activation energy for viscous deformation rather than the viscous compaction length, and (2) lateral fluid flow in viscous regimes may occur on greater length scales than anticipated from earlier estimates of compaction length scales. In viscous rock, inverted geothermal gradients stabilize vertically elongated waves or vertical channels. Decreasing temperature toward the earth’s surface can induce an abrupt transition from viscous to elastic deformation-propagated fluid flow. Below the transition, fluid flow is accomplished by short wavelength, large amplitude waves; above the transition flow is by high velocity, low amplitude surges. The resulting transient flow patterns vary strongly in space and time. Solitary porosity waves may nucleate in viscous, viscoplastic, and viscoelastic rheologies. The amplitude of these waves is effectively unlimited for physically realistic models with dependence of bulk viscosity on porosity. In the limit of zero initial connected porosity, arguably the only model relevant for melt extraction, travelling waves are only possible in a viscoelastic matrix. Such waves are truly self-propagating in that the fluid and the wave phase velocities are identical; thus, if no chemical processes occur during propagation, the waves have the capacity to transmit geochemical signatures indefinitely. In addition to solitary waves, we find that periodic solutions to the compaction equations are common though previously unrecognized. The transition between the solutions depends on the pore volume carried by the wave and the Darcyian velocity of the background fluid flux. Periodic solutions are possible for all velocities, whereas solitary solutions require large volumes and low velocities. © Elsevier, Paris     

The role of the fluid phase during regional metamorphism and deformation

Evidence from rock microstructures, mass transfer and isotopic exchange indicates that substantial quantities of aqueous fluids are involved in low- and medium-grade regional metamorphism. Similar conclusions are drawn from many retrograde environments, whereas high-grade metamorphic fluids may be melt dominated. The mobile fluids play essential roles in metamorphic reactions, mass transport and deformation processes. These processes are linked by the mechanical consequences of metamorphic fluid pressures (P_f) generally being greater than or equal to the minimum principal compressive stress. Under such conditions metamorphic porosity comprises grain boundary tubules and bubbles together with continuously generated (and healed) microfractures. Deformation results in significant interconnected porosity and hence enhanced permeability. Lithologically and structurally controlled permeability variations may cause effective fluid channelling.
Simple Rayleigh-Darcy modelling of a uniformly permeable, crustal slab shows that convective instability of metamorphic fluid is expected at the permeabilities suggested for the high P_f metamorphic conditions. Complex, large-scale convective cells operating in overpressured, but capped systems may provide a satisfactory explanation for the large fluid/rock ratios and extensive mass transport demonstrated for many low- and medium-grade metamorphic environments. Such large-scale fluid circulation may have important consequences for heat transfer in and the thermal evolution of metamorphic belts.     

Effects of fluid production on fluid flow during regional and contact metamorphism

Numerical and analytical models of fluid flow that account for fluid production during prograde regional and contact metamorphism show that expulsion of metamorphic fluids dominates the convective flux when crustal permeabilities are less than 0.1–100 μD, depending primarily on the rate of fluid production. When this is the case, fluid circulation is limited or prevented, fluid pressures are elevated above hydrostatic values, and flow throughout most of the model is up and away from the region of maximum fluid production. Fluid circulation is predicted to occur where permeability is high, in dry rocks, or after rates of fluid production decrease as peak temperatures are reached. Large changes in the pattern of flow and influx of externally derived fluids may thus occur in metamorphic terranes when dehydration wanes or ceases and cooling begins. Inclusion of an impermeable horizon in the models further inhibits fluid circulation. Earlier, shallow hydrothermal models and interpretations based on the Rayleigh number may be inappropriate for characterizing fluid flow during prograde metamorphism at depth because they do not account for fluid production.     

部分熔融与高级变质岩流变机制——以内蒙古大青山高级变质岩为例

部分熔融作用与高级变质岩变形作用是相互制约,变形作用能够提高岩石部分熔融程度,降低熔融温度。熔体存在影响和制约岩石强度和变形机制。大青山高级岩经历了下部地壳构造层次变质变形和深熔作用改造,形成了复杂构造要素组合。宏观与微观构造特点表明:高级变质岩变形机制主要为熔体增强颗粒边界扩散和颗粒流动,使岩石发生大规模的塑性流动。在宏观上形成了不对称流动组构、熔融线理、岩石和矿物条带、层内底辟褶皱和大型穹窿构造。但是,在微观上矿物颗粒变形不明显,晶内变形组构不发育,表现为三边平衡结构,与静态结晶变质岩结构相似,形成了地壳深部构造层次上变质构造岩-构造片麻岩。     

A fully implicit and consistent finite element framework for modeling reservoir compaction with large deformation and nonlinear flow model. Part I: theory and formulation

Reservoir depletion results in rock failure, wellbore instability, hydrocarbon production loss, oil sand production, and ground surface subsidence. Specifically, the compaction of carbonate reservoirs with soft rocks often induces large plastic deformation due to rock pore collapse. On the other hand, following the compaction of reservoirs and failure of rock formations, the porosity and permeability of formations will, in general, decrease. These bring a challenge for reservoir simulations because of high nonlinearity of coupled geomechanics and fluid flow fields. In this work, we present a fully implicit, fully coupled, and fully consistent finite element formulation for coupled geomechanics and fluid flow problems with finite deformation and nonlinear flow models. The Pelessone smooth cap plasticity model, an important material model to capture rock compaction behavior and a challenging material model for implicit numerical formulations, is incorporated in the proposed formulation. Furthermore, a stress-dependent permeability model is taken into account in the formulation. A co-rotational framework is adopted for finite deformation, and an implicit material integrator for cap plasticity models is consistently derived. Furthermore, the coupled field equations are consistently linearized including nonlinear flow models. The physical theories, nonlinear material and flow models, and numerical formulations are the focus of part I of this work. In part II, we verify the proposed numerical framework and demonstrate the performance of our numerical formulation using several numerical examples including a field reservoir with soft rocks undergoing serious compaction.