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.     

Age and early metamorphic history of the Sanbagawa belt: Lu–Hf and P–T constraints from the Western Iratsu eclogite

Two distinct age estimates for eclogite-facies metamorphism in the Sanbagawa belt have been proposed: (i) c.  120–110 Ma based on a zircon SHRIMP age for the Western Iratsu unit and (ii) c.  88–89 Ma based on a garnet–omphacite Lu–Hf isochron age from the Seba and Kotsu eclogite units. Despite the contrasting estimates of formation ages, petrological studies suggest the formation conditions of the Western Iratsu unit are indistinguishable from those of the other two units—all ∼20 kbar and 600–650 °C. Studies of the associated geological structures suggest the Seba and Western Iratsu units are parts of a larger semi-continuous eclogite unit. A combination of geochronological and petrological studies for the Western Iratsu eclogite offers a resolution to this discrepancy in age estimates. New Lu–Hf dating for the Western Iratsu eclogite yields an age of 115.9 ± 0.5 Ma that is compatible with the zircon SHRIMP age. However, petrological studies show that there was significant garnet growth in the Western Iratsu eclogite before eclogite facies metamorphism, and the early core growth is associated with a strong concentration of Lu. Pre-eclogite facies garnet (Grt1) includes epidote–amphibolite facies parageneses equilibrated at 550–650 °C and ∼10 kbar, and this is overgrown by prograde eclogite facies garnet (Grt2). The Lu–Hf age of c.  116 Ma is strongly skewed to the isotopic composition of Grt1 and is interpreted to reflect the age of the pre-eclogite phase. The considerable time gap ( c.  27 Myr) between the two Lu–Hf ages suggests they may be related to separate tectonic events or distinct phases in the evolution of the Sanbagawa subduction zone.     

Numerical manifold method for dynamic nonlinear analysis of saturated porous media

It is well known that the Babuska–Brezzi stability criterion or the Zienkiewicz–Taylor patch test precludes the use of the finite elements with the same low order of interpolation for displacement and pore pressure in the nearly incompressible and undrained cases, unless some stabilization techniques are introduced for dynamic analysis of saturated porous medium where coupling occurs between the displacement of solid skeleton and pore pressure. The numerical manifold method (NMM), where the interpolation of displacement and pressure can be determined independently in an element for the solution of u–p formulation, is derived based on triangular mesh for the requirement of high accurate calculations from practical applications in the dynamic analysis of saturated porous materials. The matrices of equilibrium equations for the second‐order displacement and the first‐order pressure manifold method are given in detail for program coding. By close comparison with widely used finite element method, the NMM presents good stability for the coupling problems, particularly in the nearly incompressible and undrained cases. Numerical examples are given to illustrate the validity and stability of the manifold element developed. Copyright © 2006 John Wiley & Sons, Ltd.     

流褶层及韧变带研究

流褶层与韧变带是地壳拉伸变形，顺层固态流变作用下的产物。流褶层是以原始层理为变形面或再经递进变形的褶皱变形岩层或岩石共生组合层位。韧变带具明显的层控性，受岩石成分和应变程度控制，不同环境和不同成分岩石的韧变带具有相异的组合型式和变形机制，井具有一定的递变规律。流褶层和韧变带可分属不同层位，但流褶层可实现向韧变带的转化。     

About Darcy's law in non-Galilean frame

This paper is aimed towards investigating the filtration law of an incompressible viscous Newtonian fluid through a rigid non-inertial porous medium (e.g. a porous medium placed in a centrifuge basket). The filtration law is obtained by upscaling the flow equations at the pore scale. The upscaling technique is the homogenization method of multiple scale expansions which rigorously gives the macroscopic behaviour and the effective properties without any prerequisite on the form of the macroscopic equations. The derived filtration law is similar to Darcy's law, but the tensor of permeability presents the following remarkable properties: it depends upon the angular velocity of the porous matrix, it verifies Hall–Onsager's relationship and it is a non-symmetric tensor. We thus deduce that, under rotation, an isotropic porous medium leads to a non-isotropic effective permeability. In this paper, we present the results of numerical simulations of the flow through rotating porous media. This allows us to highlight the deviations of the flow due to Coriolis effects at both the microscopic scale (i.e. the pore scale), and the macroscopic scale (i.e. the sample scale). The above results confirm that for an isotropic medium, phenomenological laws already proposed in the literature fails at reproducing three-dimensional Coriolis effects in all types of pores geometry. We show that Coriolis effects may lead to significant variations of the permeability measured during centrifuge tests when the inverse Ekman number Ek⁻¹ is 𝒪(1). These variations are estimated to be less than 5% if Ek⁻¹<0.2, which is the case of classical geotechnical centrifuge tests. We finally conclude by showing that available experimental data from tests carried out in centrifuges are not sufficient to determining the effective tensor of permeability of rotating porous media. Copyright © 2004 John Wiley & Sons, Ltd.     

Nutrient budgets and trophic state in a hypersaline coastal lagoon: Lagoa de Araruama, Brazil

Lagoa de Araruama in the state of Rio de Janeiro, Brazil, is a hypersaline lagoon with salinity varying spatially from 45 to 56. We collected water samples during monthly cruises throughout the lagoon, and along the streams feeding the system, from April 1991 to March 1992. Nutrients and other water quality parameters exhibited great spatial and temporal variations. Mass balance calculations indicate large amounts of anthropogenic nutrient inputs. The data indicate that the lagoon currently is oligotrophic but is in a state of transition to become a mesotrophic system. Molar dissolved inorganic nitrogen:dissolved inorganic phosphorus (DIN/DIP) varied between 2.2:1 and 659:1 with a volume-weighted average of 22:1. The high DIN/DIP ratio contrasts with that found in nearby lagoons, suggesting that phytoplankton primary production is limited by phosphorus in Lagoa de Araruama. The major loss of DIP is apparently driven by biological assimilation and diagenic reactions in the sediments. Calculations indicate that the lagoon is slightly net autotrophic at +0.9 mol C m⁻² yr⁻¹. This suggests that the biomass of the primary producers is restricted by phosphorus availability. Phosphorus retention in the sediment and the hypersaline state of the lagoon prevent changes in autotrophic communities and the formation of eutrophic conditions.     

太平洋三海区热液烟囱物的氧同位素特征及形成温度

使用氧同位素方法 ,测试了西太平洋马里亚纳岛弧、马里亚纳海槽、冲绳海槽和东太平洋加拉帕戈斯裂谷的海底热液烟囱和硫化物全岩样品。结果表明 :马里亚纳岛弧上的埃斯梅拉尔达破火山口的热液硫化物的形成温度最高 (达 2 75℃ ) ,而成为高温产物的代表 ;冲绳海槽和加拉帕戈斯裂谷烟囱 ,表现为中高温类型 ,其氧同位素温度在 1 5 0℃左右 ;马里亚纳海槽以蛋白石为主要矿物的烟囱物 ,则显示出典型的低温热液类型 ,少量黄铁矿砂试样则表现出中温类型特征。但上述区域大都不同程度地存在有闪锌矿、黄铁矿、黄铜矿等中高温热液矿物 ,说明研究区除马里亚纳海槽明显有中低温热液活动类型外 ,其它区域皆具备中高温热液活动特征     

井中激电地-井方式井旁球体正反演

以激发极化法中的地-井方式中的体极化球形体作为研究对象,采用解析法对地下半空间的电位场进行求解,利用提取出的井中二次场异常电位差进行反演,得出极化球体的球心离井的距离和球心的埋深。通过计算机进行正反演模拟,验证了进行正反演计算的解析式的正确性,并从物理意义上对正反演过程中出现的现象进行分析。最后给出该方法的适用条件以及注意事项。研究对激发极化法地下勘探目标的定量分析具有明确意义。     

Geochemistry of the Mian-Lue ophiolites in the Qinling Mountains, central China: Constraints on the evolution of the Qinling orogenic belt and collision of the North and South China Cratons

The Anzishan ophiolite, a typical ophiolitic block of early Carboniferous age in the Mian-Lue suture zone of the Qinling Mountains, central China, consists of amphibolites/metabasalts, gabbros and gabbroic cumulates. All of these rocks, as well as those in the Hunshuiguan-Zhuangke (HZ) block, have compositions similar to normal MORB and back-arc basin basalts (BABB) with high ε_Nd(t) values, indicating that they were derived from a depleted mantle source. The Mian-Lue suture zone also contains blocks of other lithologies, e.g., rift volcanic rocks in the Heigouxia block and arc volcanic rocks in the Sanchazi block. Although they are in fault contact with each other, the presence of these different blocks in the Mian-Lue suture zone may represent a complete Wilson cycle, from initial rifting to open ocean basin to final subduction and continent-continent collision, during the late Paleozoic-early Triassic. In this region, the North and South China Cratons were separated by Paleo-Tethys at least until the early Carboniferous, and final amalgamation of both cratons along the Qinling orogenic belt took place in the Triassic.     

Secular variations of the upper crust composition: Implication of geochemical data on the Upper Precambrian shales from the Southern Urals western flank and Uchur-Maya region

In the mid-1980s, it was concluded based on geochemical study that Th, Sc, La concentrations and ratios Th/Sc, La/Sc and Eu/Eu* did not wary significantly in the post-Archean time. It was impossible to judge about compositional variations of upper crust during the Riphean and Vendian, because data of that time characterized a limited number of samples from the post-Archean basins of Australia, New Zealand, and Antarctic. Considered in this work are variations of Eu/Eu*, LREE/HREE, Th/Sc, and La/Sc ratios in Upper Precambrian fine-grained siliciclastic rock of the Southern Urals western flank (Bashkirian meganticlinorium) and Uchur-Maya region (Uchur-Maya plate and Yudoma-Maya belt). As is established, only the Eu anomaly in the studied siliciclastic rocks is practically identical to this parameter of the average post-Archean shale. Three other parameters plot on the Riphean-Vendian variation curves with positive and negative excursions of diverse magnitude, which do not coincide always in time. It is assumed that these excursions likely mark stages of local geodynamic activity, destruction of pre-Riphean cratons, and progressing recycling of sedimentary material during the Riphean.