Effects of the geometry of two‐dimensional fractures on their hydraulic aperture and on the validity of the local cubic law

Flow through rough fractures is investigated numerically in order to assess the validity of the local cubic law for different fracture geometries. Two‐dimensional channels with sinusoidal walls having different geometrical properties defined by the aperture, the amplitude, and the wavelength of the walls' corrugations, the corrugations asymmetry, and the phase shift between the two walls are considered to represent different fracture geometries. First, it is analytically shown that the hydraulic aperture clearly deviates from the mean aperture when the walls' roughness, the phase shift, and/or the asymmetry between the fracture walls are relatively high. The continuity and the Navier–Stokes equations are then solved by means of the finite element method and the numerical solutions compared to the theoretical predictions of the local cubic law. Reynolds numbers ranging from 0.066 to 66.66 are investigated so as to focus more particularly on the effect of flow inertial effects on the validity of the local cubic law. For low Reynolds number, typically less than 15, the local cubic law properly describes the fracture flow, especially when the fracture walls have small corrugation amplitudes. For Reynolds numbers higher than 15, the local cubic law is valid under the conditions that the fracture presents a low aspect ratio, small corrugation amplitudes, and a moderate phase lag between its walls.     

Solution for a plane strain rough‐walled hydraulic fracture driven by turbulent fluid through impermeable rock

The impact of turbulent flow on plane strain fluid‐driven crack propagation is an important but still poorly understood consideration in hydraulic fracture modeling. The changes that hydraulic fracturing has experienced over the past decade, especially in the area of fracturing fluids, have played a major role in the transition of the typical fluid regime from laminar to turbulent flow. Motivated by the increasing preponderance of high‐rate, water‐driven hydraulic fractures with high Reynolds number, we present a semianalytical solution for the propagation of a plane strain hydraulic fracture driven by a turbulent fluid in an impermeable formation. The formulation uses a power law relationship between the Darcy‐Weisbach friction factor and the scale of the fracture roughness, where one specific manifestation of this generalized friction factor is the classical Gauckler‐Manning‐Strickler approximation for turbulent flow in a rough‐walled channel. Conservation of mass, elasticity, and crack propagation are also solved simultaneously. We obtain a semianalytical solution using an orthogonal polynomial series. An approximate closed‐form solution is enabled by a choice of orthogonal polynomials embedding the near‐tip asymptotic behavior and thus giving very rapid convergence; a precise solution is obtained with 2 terms of the series. By comparison with numerical simulations, we show that the transition region between the laminar and turbulent regimes can be relatively small so that full solutions can often be well approximated by either a fully laminar or fully turbulent solution.     

Micromechanical approach to effective viscoelastic properties of micro‐fractured geomaterials

The aim of this paper is to formulate a micromechanics‐based approach to non‐aging viscoelastic behavior of materials with randomly distributed micro‐fractures. Unlike cracks, fractures are discontinuities that are able to transfer stresses and can therefore be regarded from a mechanical viewpoint as interfaces endowed with a specific behavior under normal and shear loading. Making use of the elastic‐viscoelastic correspondence principle together with a Mori‐Tanka homogenization scheme, the effective viscoelastic behavior is assessed from properties of the material constituents and damage parameters related to density and size of fractures. It is notably shown that the homogenized behavior thus formulated can be described in most cases by means of a generalized Maxwell rheological model. For practical implementation in structural analyses, an approximate model for the isotropic homogenized fractured medium is formulated within the class of Burger models. Although the approximation is basically developed for short‐term and long‐term behaviors, numerical applications indicate that the approximate Burger model accurately reproduce the homogenized viscoelastic behavior also in the transient conditions.     

Late Quaternary submarine bedforms and ice‐sheet flow in Gerlache Strait and on the adjacent continental shelf,Antarctic Peninsula

Geophysical data from Gerlache Strait, Croker Passage, Bismarck Strait and the adjacent continental shelf reveal streamlined subglacial bedforms that were produced at the bed of the Antarctic Peninsula Ice Sheet (APIS) during the last glaciation. The spatial arrangement and orientation of these bedforms record the former drainage pattern and flow dynamics of an APIS outlet up‐flow, and feeding into, a palaeo‐ice stream in the Western Bransfield Basin. Evidence suggests that together, they represent a single ice‐flow system that drained the APIS during the last glaciation. The ice‐sheet outlet flowed north/northeastwards through Gerlache Strait and Croker Passage and converged with a second, more easterly ice‐flow tributary on the middle shelf to form the main palaeo‐ice stream. The dominance of drumlins with low elongation ratios suggests that ice‐sheet outlet draining through Gerlache Strait was comparatively slower than the main palaeo‐ice stream in the Western Bransfield Basin, although the low elongation ratios may also partly reflect the lack of sediment. Progressive elongation of drumlins further down‐flow indicates that the ice sheet accelerated through Croker Passage and the western tributary trough, and fed into the main zone of streaming flow in the Western Bransfield Basin. Topography would have exerted a strong control on the development of the palaeo‐ice stream system but subglacial geology may also have been significant given the transition from crystalline bedrock to sedimentary strata on the inner–mid‐shelf. In the broader context, the APIS was drained by a number of major fast‐flowing outlets through cross‐shelf troughs to the outer continental shelf during the last glaciation. Copyright © 2004 John Wiley & Sons, Ltd.     

Evolution of 2:1 layered silicates in low-grade metamorphosed Liassic shales of Central Switzerland

White mica from the Liassic black shales and slates in Central Switzerland was analysed by transmission electron microscopy (TEM) and electron microprobe to determine its textural and compositional evolution during very low-grade prograde metamorphism. Samples were studied from the diagenetic zone, anchizone and epizone (T ≈100°–450 °C). Phyllosilicate minerals analysed include illite/smectite (I/S), phengite, muscovite, brammallite, paragonite, margarite and glauconite. Textural evolution primarily is towards larger, more defect-free grains with compositions that approach those of their respective end-members. The smectite-to-illite transformation reduced the amounts of the exchange components SiK_?1Al_?1, MgSiAl_?2, and Fe³⁺Al_?1. These trends continue to a lesser degree in the anchizone and epizone. Correlations between the proportion of smectite in I/S and the composition of I/S indicate that smectite layers may contain a high layer charge. Illite in I/S bears a compositional resemblance to macrocrystalline phengite in some samples, but is different in others. Paragonite first appears in the upper diagenetic zone or lower anchizone as an interlayer-deficient brammallite, and it may be mixed with muscovite on the nanometre scale. Owing to the small calculated structure factor for paragonite-muscovite superstructures, conventional X-ray powder diffraction cannot distinguish between mixed-layer structures and a homogeneous compositionally intermediate solid solutions. However, indirect TEM evidence shows that irregularly shaped domains of Na- and K-rich mica exist below 10 nm. Subsequent coarsening of domains at higher grades produced discrete paragonite grains at the margins of muscovite crystals or in laths parallel to the basal plane of the host muscovite. Margarite appears in the epizone and follows a textural evolution similar to paragonite in that mixtures of margarite, paragonite, and muscovite may initially occur on the nanometre scale. However, no evidence of interlayer-poor margarite has been found.     

High-pressure eclogites in northern Jiangsu – southern Shandong province, eastern China

Abstract Eclogites are distributed for more than 500 km along a major tectonic boundary between the Sino-Korean and Yangtze cratons in central and eastern China. These eclogites usually have high-P assemblages including omphacite + kyanite and/or coesite (or its pseudomorph), and form a high-P eclogite terrane. They occur as isolated lenses or blocks 10 cm to 300 m long in gneisses (Type I), serpentinized garnet peridotites (Type II) and marbles (Type III). Type I eclogites were formed by prograde metamorphism, and their primary metamorphic mineral assemblage consists mainly of garnet [pyrope (Prp) = 15–40 mol%], omphacite [jadeite (Jd) = 34–64 mol%], pargasitic amphibole, kyanite, phengitic muscovite, zoisite, an SiO₂ phase, apatite, rutile and zircon. Type II eclogites characteristically contain no SiO₂ phase, and are divided into prograde eclogites and mantle-derived eclogites. The prograde eclogites of Type II are petrographically similar to Type I eclogites. The mantle-derived eclogites have high MgO/(FeO + Fe₂O₃) and Cr₂O₃ compositions in bulk rock and minerals, and consist mainly of pyrope-rich garnet (Prp = 48–60 mol%), sodic augite (Jd = 10–27 mol%) and rutile. Type III eclogites have an unusual mineral assemblage of grossular-rich (Grs = 57 mol%) garnet + omphacite (Jd = 30–34 mol%) + pargasite + rutile. Pargasitic and taramitic amphiboles, calcic plagioclase (An₆₈), epidote, zoisite, K-feldspar and paragonite occur as inclusions in garnet and omphacite in the prograde eclogites. This suggests that the prograde eclogites were formed by recrystallization of epidote amphibolite and/or amphibolite facies rocks with near-isothermal compression reflecting crustal thickening during continent–continent collision of late Proterozoic age. Equilibrium conditions of the prograde eclogites range from P > 26 kbar and T= 500–750°C in the western part to P > 28 kbar and T= 810–880°C in the eastern part of the high-P eclogite terrane. The prograde eclogites in the eastern part are considered to have been derived from a deeper position than those in the western part. Subsequent reactions, manifested by (1) narrow rims of sodic plagioclase or paragonite on kyanite and (2) symplectites between omphacite and quartz are interpreted as an effect of near-isothermal decompression during the retrograde stage. The conditions at which symplectites re-equilibrated tend to increase from west (P < 10 kbar and T < 580°C) to east (P > 9 kbar and T > 680°C). Equilibrium temperatures of Type II mantle-derived eclogites and Type III eclogite are 730–750°C and 680°C, respectively.     

Late Quaternary landscapes in Central Australia: sedimentary history and palaeoecology of Puritjarra rock shelter

Puritjarra rock shelter provides a long record of late Quaternary vegetation in the Australian arid zone. Analysis of the sedimentary history of this rock shelter is combined with reanalysis of charcoal and phytolith records to provide a first‐order picture of changing landscapes in western Central Australia. These show a landscape responding to increasing aridity from 45 ka with deflation of clay‐rich red palaeosols (<45 ka) and sharp declines in grassland and other vegetation at 40–36 ka, and at the beginning of the Last Glacial Maximum (LGM) (24 ka). Vegetation in the catchment of the rock shelter recovered after 15 ka with expansion of both acacia woodland and spinifex grasslands, registering stronger summer rainfall in the interior of the continent. By 8.3 ka re‐vegetation of local palaeosols and dunes had choked off sediment supply to the rock shelter and the character of the sediments changed abruptly. Poaceae values peaked at 5.8 ka, suggesting the early–mid Holocene climatic optimum in Central Australia is bracketed between 8.3 and 5.8 ka. Local vegetation was disrupted in the late Holocene with a sharp decline in Poaceae at 3.8 ka, coinciding with an abrupt intensification of ENSO. Local grasslands recovered over the next two millennia and by 1.5 ka the modern vegetation appears to have become established. Copyright © 2009 John Wiley & Sons, Ltd.     

江西浒坑钨矿含矿石英脉的地质特征及成矿构造演化

位于赣西武功山地区的江西浒坑钨矿是一个大型石英脉型钨矿床,与赣南典型的“五层楼”模式不同,矿脉主要集中在岩体内接触带。含矿石英脉除了普通的块状构造以外,还普遍存在条带状构造。详细的野外编录和镜下研究表明含矿块状石英脉在形成后遭受不同程度的韧性剪切,从而形成不同成分的分异条带。剪切带以出现新生的面理和线理（拉伸线理）为特征。岩石学研究表明条带状矿石属于糜棱岩,具有丰富的塑性变形显微构造特征。通过对含矿石英脉的构造解析,认为与岩浆岩有关的原生节理系统基本上是同时形成的,沿节理充填的石英脉经历了长期的、多期次的构造运动,从而形成了复杂的石英脉和矿石类型。建立了原生节理形成→石英充填→韧性剪切发生→二期石英充填→脆性碎裂的矿区构造演化序列。韧性剪切作用的动力学背景可能是受后期岩浆侵入的影响,最后的脆性破裂应与区域上存在南北向的压应力有关。