Three dimensional simulation of fluid flow in X‐ray CT images of porous media

A numerical scheme is developed in order to simulate fluid flow in three dimensional (3‐D) microstructures. The governing equations for steady incompressible flow are solved using the semi‐implicit method for pressure‐linked equations (SIMPLE) finite difference scheme within a non‐staggered grid system that represents the 3‐D microstructure. This system allows solving the governing equations using only one computational cell. The numerical scheme is verified through simulating fluid flow in idealized 3‐D microstructures with known closed form solutions for permeability. The numerical factors affecting the solution in terms of convergence and accuracy are also discussed. These factors include the resolution of the analysed microstructure and the truncation criterion. Fluid flow in 2‐D X‐ray computed tomography (CT) images of real porous media microstructure is also simulated using this numerical model. These real microstructures include field cores of asphalt mixes, laboratory linear kneading compactor (LKC) specimens, and laboratory Superpave gyratory compactor (SGC) specimens. The numerical results for the permeability of the real microstructures are compared with the results from closed form solutions. Copyright © 2004 John Wiley & Sons, Ltd.     

A continuum‐discrete hydromechanical analysis of granular deposit liquefaction

A coupled continuum‐discrete hydromechanical model was employed to analyse the liquefaction of a saturated loose deposit of cohesionless particles when subjected to a dynamic base excitation. The pore fluid flow was idealized using averaged Navier–Stokes equations and the discrete element method was employed to model the solid phase particles. A well established semi‐empirical relationship was utilized to quantify the fluid–particle interactions. The conducted simulations revealed a number of salient micro‐mechanical mechanisms and response patterns associated with the deposit liquefaction. Space and time variation of porosity was a major factor which affected the coupled response of the solid and fluid phases. Pore fluid flow was within Darcy's regime. The predicted response exhibited macroscopic patterns consistent with experimental results and case histories of the liquefaction of granular soil deposits. Copyright © 2004 John Wiley & Sons, Ltd.     

Effect of metamorphism and partial melting of host rocks on zircons

Abstract Zircons have been studied from different layers of migmatites (from Arvika, western Sweden and Nelaug, southern Norway) and from a paragneiss (from Arvika) associated with one of the migmatites. The main purpose of the investigation is to establish whether or not information about zircons can help in the elucidation of the parentage and rock-forming processes of migmatites.
The elongation ratio of zircons from all layers is small and characteristic of sedimentary zircons. Further, the absence of characteristic colours and the growth trends of the zircons (indicated by the reduced major axes) observed in the various samples both support a sedimentary parentage for these rocks. The zircons of all layers exhibit secondary growth (overgrowth, outgrowth and multiple growth) due to metamorphism. Compared with the zircons from the paragneiss, those of the migmatite layers are more clouded and less rounded, some of them becoming opaque or even skeletal; this is especially true of the zircons from the leucosomes. These observations indicate an alteration of the original sedimentary zircons in the migmatite, especially in the leucosomes, in response to the migmatization process, previously interpreted as partial melting.     

差分法周跳探测与修复方法改进

在目前常用的周跳探测与修复方法基础上 ,提出了首先将观测资料按照观测历元不连续分成若干小弧段 ,然后利用差分法进行周跳探测 ,根据差分后周跳放大的特性判断周跳和野值 ,并确定其位置利用宽带组合和电离层组合的方法解算周跳大小。通过实例验证了其有效性。     

Microboudin structures of piemontite along the Sambagawa metamorphic belt, Japan: implications for lateral variation of differential stress during regional metamorphism

Evidence is presented of a lateral variation in differential stress during metamorphism along a regional metamorphic belt on the basis of the proportion of microboudinaged piemontite grains (p) in a quartz matrix in metacherts. It is proposed that p is a practical indicator of relative differential stress. Analysis of 123 metacherts from the 800 km long Sambagawa metamorphic belt, Japan, reveals that p‐values range from < 0.01 to 0.7 in this region. Most samples from Wakayama in the mid‐belt area have p‐values of 0.4–0.6, whereas those from western Shikoku have p‐values of < 0.1. This difference cannot be explained by variations in metamorphic temperature, and is instead attributed to a regional, lateral variation in differential stress during metamorphism.     

Petrogenesis of Adakitic Intrusive Rocksin Xu-Huai Area, Eastern China

1IntroductionAdakite has specially sense in the resuming structure background of rock formed,because of its distinctive geochemistry feature and implying partic-ular petrogenesis and ore genesis mechanism (De-fant and Drummond,1990),as a result,it arouses numerous researchers widely concern. After adakite conception was drawn by Chinese (Wang etal.,2000), it was widely concerned by domestic re-searchers (Wang etal.,2000,2001a,2001b;Xu etal.,2000;Pan etal.,2001;Qian,2001;Xu etal.,2001;X…     

Seismic rays in an ellipsoidal earth model: theoretical analysis and estimation of deviations due to the ellipticity

Summary. A first-order form of the Euler's equations for rays in an ellipsoidal model of the Earth is obtained. The conditions affecting the velocity law for a monotonic increase, with respect to the arc length, in the angular distance to the epicentre, and in the angle of incidence, are the same in the ellipsoidal and spherical models. It is therefore possible to trace rays and to compute travel times directly in an ellipsoidal earth as in the spherical model. Thus comparison with the rays of the same coordinates in a spherical earth provides an estimate of the various deviations of these rays due to the Earth's flattening, and the corresponding travel-time differences, for mantle P -waves and for shallow earthquakes. All these deviations are functions both of the latitude and of the epicentral distance. The difference in the distance to the Earth's centre at points with the same geocentric latitude on rays in the ellipsoidal and in the spherical model may reach several kilometres. Directly related to the deformation of the isovelocity surfaces, this difference is the only cause of significant perturbation in travel times. Other differences, such as that corresponding to the ray torsion, are of the first order in ellipticity, and may exceed 1 km. They induce only small differences in travel time, less than 0.01s. Thus, we show that the ellipticity correction obtained by Jeffreys (1935) and Bullen (1937) by a perturbational method can be recovered by a direct evaluation of the travel times in an ellipsoidal model of the Earth. Moreover, as stated by Dziewonski & Gilbert (1976), we verify the non-dependence of this correction on the choice of the velocity law.     

Crustal and upper mantle structures of the brazilian coast

The Rayleigh wave phase and group velocities in the period range of 24–39 sec, obtained from two earthquakes which occurred in northeastern brazil and which were recorded by the Brazilian seismological station RDJ (Rio de Janeiro), have been used to study crustal and upper mantle structures of the Brazilian coastal region. Three crustal and upper mantle models have been tried out to explain crustal and upper mantle structures of the region. The upper crust has not been resolved, due basically to the narrow period range of the phase and group velocities data. The phase velocity inversions have exhibited good resolutions for both lower crust and upper mantle, with shear wave velocities characteristic of these regions. The group velocity data inversions for these models have showed good results only for the lower crust. The shear wave velocities of the lower crust (3.86 and 3.89 km/sec), obtained with phase velocity inversions, are similar to that (=3.89 km/sec) found byHwang (1985) to the eastern South American region, while group velocity inversions have presented shear velocity (=3.75 km/sec) similar to that (=3.78 km/sec) found byLazcano (1972) to the Brazilian shield. It was not possible to define sharply the crust-mantle transition, but an analysis of the phase and group velocity inversions results has indicated that the total thickness of the crust should be between 30 and 39 km. The crustal and upper mantle model, obtained with phase velocity inversion, can be used as a preliminary model for the Brazilian coast.     

Partial crustal melting beneath the Betic Cordillera (SE Spain): The case study of Mar Menor volcanic suite

The Neogene Volcanic Province (NVP) within the Betic Cordillera (SE Spain) consists of three main metapelitic enclave suites (from SW to NE: El Hoyazo, Mazarrón and Mar Menor). Since the NVP represents a singular place in the world where crustal enclaves were immediately quenched after melting, their microstructures provide a “photograph” of the conditions at depth just after the moment of the melting.

The thermobarometric information provided by the different microstructural assemblages has been integrated with the geophysical and geodynamical published data into a model of the petrologic evolution of the Mar Menor enclaves. They were equilibrated at 2–3 kbar, 850–900 °C, and followed a sequence of heating melt producing reactions. A local cooling event evidenced by minor melt crystallization preceded the eruption.

The lower crustal studies presented in this work contribute to the knowledge of: (i) the partial melting event beneath the Mar Menor volcanic suite through a petrologic detailed study of the enclaves; (ii) how the microstructures of fast cooled anatectic rocks play an important role in tracing the magma evolution in a chamber up to the eruption, and how they can be used as pseudothermobarometers; (iii) the past and current evolution of the Alborán Domain (Betic Cordillera) and Mediterranean Sea, and how the base of a metapelitic crust has melted within an active geodynamic setting.     

Partial Melting and Counterclockwise P T Path of Subducted Oceanic Crust (Sierra del Convento Melange, Cuba)

Partial melting of subducted oceanic crust has been identifiedin the Sierra del Convento mélange (Cuba). This serpentinite-matrixmélange contains blocks of mid-ocean ridge basalt (MORB)-derivedplagioclase-lacking epidote ± garnet amphibolite intimatelyassociated with peraluminous trondhjemitic–tonalitic rocks.Field relations, major element bulk-rock compositions, mineralassemblages, peak metamorphic conditions (c. 750°C, 14–16kbar), experimental evidence, and theoretical phase relationssupport formation of the trondhjemitic–tonalitic rocksby wet melting of subducted amphibolites. Phase relations andmass-balance calculations indicate eutectic- and peritectic-likemelting reactions characterized by large stoichiometric coefficientsof reactant plagioclase and suggest that this phase was completelyconsumed upon melting. The magmatic assemblages of the trondhjemitic–tonaliticmelts, consisting of plagioclase, quartz, epidote, ±paragonite, ± pargasite, and ± kyanite, crystallizedat depth (14–15 kbar). The peraluminous composition ofthe melts is consistent with experimental evidence, explainsthe presence of magmatic paragonite and (relict) kyanite, andplaces important constraints on the interpretation of slab-derivedmagmatic rocks. Calculated P–T conditions indicate counterclockwiseP–T paths during exhumation, when retrograde blueschist-faciesoverprints, composed of combinations of omphacite, glaucophane,actinolite, tremolite, paragonite, lawsonite, albite, (clino)zoisite,chlorite, pumpellyite and phengite, were formed in the amphibolitesand trondhjemites. Partial melting of subducted oceanic crustin eastern Cuba is unique in the Caribbean realm and has importantconsequences for the plate-tectonic interpretation of the region,as it supports a scenario of onset of subduction of a youngoceanic lithosphere during the early Cretaceous (c. 120 Ma).The counterclockwise P–T paths were caused by ensuingexhumation during continued subduction. KEY WORDS: amphibolite; Cuba; exhumation; partial melting; trondhjemite; subduction