高温高压下橄榄石含水实验和软流圈成因

高温高压橄榄石含水实验表明，水进入橄榄石能加速其塑性形变，促进了地幔塑性流动，为上部软流圈成因提供一种新的可能依据。     

Compressibility of the calcium aluminosilicate, CAS, phase to 44 GPa

In situ X-ray diffraction measurements on a calcium aluminosilicate (CAS) phase have been carried out using a laser-heated diamond anvil cell up to a pressure of 44 GPa, employing a synchrotron radiation source. CAS is the major mineral formed from sediments subducted into the Earth's mantle. The sample was heated using a YAG laser after each pressure increment to relax the deviatoric stress in the sample. X-ray diffraction measurements were carried out at T = 300 K using an angle-dispersive technique. The pressure was calculated using an internal platinum metal pressure calibrant. The Birch–Murnaghan equation of state for the CAS phase obtained from the experimental unit cell parameters showed a density of ρ₀ = 3.888 g/cm³ and a bulk modulus of K₀ = 229 ± 9 GPa for K^′₀ = 4.7 ± 0.7. When the first pressure derivative of the bulk modulus was fixed at K^′₀ = 4, then the value of K₀ = 239 ± 2 GPa. From the experimental compressibility, the density of the CAS phase was observed to be lower than the density of co-existing Al-bearing stishovite, calcium perovskite, calcium ferrite-type phases, and (Fe,Al)-bearing Mg-perovskite in subducted sediments in the lower mantle. Therefore, the density of subducted sediments in the lower mantle decreases with increasing mineral proportion of the CAS phase.     

Phase relations in hydrous MORB at 18-28 GPa: implications for heterogeneity of the lower mantle

The phase relations in hydrous and anhydrous mid-ocean ridge basalt were determined at pressures of 18-28 GPa. Liquidus phase relations in hydrous and anhydrous MORB are different. Garnet is the liquidus phase at pressures below 21 GPa, Ca-Al (CAS) phase and stishovite are the liquidus phases at pressures of 22-27 GPa, and stishovite and Ca-perovskite are the liquidus phases above 27 GPa, whereas Ca-perovskite is a liquidus phase of anhydrous MORB at pressures above 23 GPa. Under subsolidus conditions, we have found that in the hydrous MORB system the stability fields of Al-bearing perovskite and Na-Al (NAL) phase might shift to lower pressure by about 1.5 GPa compared to the dry MORB system. This shift could be explained by oxidation of a garnet-bearing assemblage by hydrous fluid and formation of Fe³⁺-bearing aluminous perovskite at lower pressures relative to the anhydrous system and/or differences in water solubility of the phases existing in perovskite-bearing assemblages. Our data indicate that hydrous basaltic crust remains denser than peridotite along the geotherm of a subducting slab, i.e. there is no density crossover between peridotite and basalt. Therefore, in slabs going through the 660 km discontinuity, basalt would gravitationally sink into the lower mantle under relatively hydrous conditions. The delamination of former basaltic crust near the 660 km discontinuity might be possible under relatively dry conditions of subduction. There are no stable highly hydrous phases in MORB above 10 GPa even at lower temperatures corresponding to subducting slabs. Therefore, MORB cannot be an important carrier of water to the deep Earth interior. However, it can be constantly supplied by water-bearing fluid from the underlying peridotite part of the descending slab. Thus, it is plausible that water can control subduction of the oceanic crust into the lower mantle.     

Thermal equation of state of magnesite to 32 GPa and 2073 K

Pressure–volume–temperature relations have been measured to 32 GPa and 2073 K for natural magnesite (Mg_0.975Fe_0.015Mn_0.006Ca_0.004CO₃) using synchrotron X-ray diffraction with a multianvil apparatus at the SPring-8 facility. A least-squares fit of the room-temperature compression data to a third-order Birch–Murnaghan equation of state (EOS) yielded K₀ = 97.1 ± 0.5 GPa and K′ = 5.44 ± 0.07, with fixed V₀ = 279.55 ± 0.02 Å³. Further analysis of the high-temperature compression data yielded the temperature derivative of the bulk modulus (∂K_T/∂T)_P = −0.013 ± 0.001 GPa/K and zero-pressure thermal expansion α = a₀ + a₁T with a₀ = 4.03 (7) × 10⁻⁵ K⁻¹ and a₁ = 0.49 (10) × 10⁻⁸ K⁻². The Anderson–Grüneisen parameter is estimated to be δ_T = 3.3. The analysis of axial compressibility and thermal expansivity indicates that the c-axis is over three times more compressible (K_Tc = 47 ± 1 GPa) than the a-axis (K_Tc = 157 ± 1 GPa), whereas the thermal expansion of the c-axis (a₀ = 6.8 (2) × 10⁻⁵ K⁻¹ and a₁ = 2.2 (4) × 10⁻⁸ K⁻²) is greater than that of the a-axis (a₀ = 2.7 (4) × 10⁻⁵ K⁻¹ and a₁ = −0.2 (2) × 10⁻⁸ K⁻²). The present thermal EOS enables us to accurately calculate the density of magnesite to the deep mantle conditions. Decarbonation of a subducting oceanic crust containing 2 wt.% magnesite would result in a 0.6% density reduction at 30 GPa and 1273 K. Using the new EOS parameters we performed thermodynamic calculations for magnesite decarbonation reactions at pressures to 20 GPa. We also estimated stability of magnesite-bearing assemblages in the lower mantle.     

人工多晶橄榄石样品孔隙压缩高温高压实验研究

人工制备多晶橄榄石样品的高孔隙度是实验关注的问题之一。本文着重研究在给定温度条件下人工多晶橄榄石样品高温高压孔隙压缩问题。文中给出了含水的湿样品和干燥样品孔隙度与实验压力的关系曲线,并提出了衡量样品压实程度的压实系数Ｖｏ／Ｖ的概念。     

高温高压条件下纳米材料强度的测试和分析方法

基于高温高压同步辐射X射线衍射和飞行时间中子衍射实验,我们提出了一种通过衍射研究纳米材料本征力学性能的方法.实验通过高压压缩实验对纳米晶和微米晶Ni的变形行为进行了对比研究;二者在不同压力下的晶胞体积数据表明,纳米晶Ni具有明显的弹性软化现象,体弹模量比微米晶Ni减小约10%.研究结果表明,纳米晶Ni的整体抗压缩性减小是由于该材料纳米晶壳层呈现为张应力状态,具有较低的抗压缩性,这与分子动力学模拟结果一致.基于该方法,通过定量分析衍射峰的宽化和偏移,我们可以获得晶粒之间的高应力聚集区所产生的"微观/局部"屈服,以及整个样品在应力下所产生的"宏观/整体"屈服.我们的应变/应力图形分析方法还可以研究材料的屈服强度、高温高压条件下的晶粒破碎/长大、加工硬化/软化,以及多晶体中的固有残余/表面应变.本文提出的方法可以用于岩石矿物的高温高压流变强度的定量测定,为研究岩石圈以及地幔物质的蠕变对流提供了高科技的研究手段.     

3 GPa固体传压三轴流变仪实验方法介绍

介绍“3 GPa固体传压三轴流变仪”的结构、性能及其实验方法与用途,阐述在此设备上进行的关于深部(上地幔高温高压条件下)岩石流变特性的实验以及岩石蠕变.应力松驰现象的研究.     

Conductivity of NaCl solution at 0.4–5.0 GPa and 25–500 °C

NaCI-H₂O is the most fundamental ternary system in geology. Until now, the measurements of electrical conductivity of NaCl solutions are still little at high pressures (> O.5 GPa) We measured the conductivity of 0.01 m NaCl solution at 0.4–5.0 GPa and 25-500°C. The results are consistent with that of Quist and Marshall (1968) at 0.4 GPa. The conductivity of NaCl solution increases with increasing temperature. The results also show that the conductivity of NaCl solution changes little with increasing pressure below 1.5 GPa and changes rapidly with increasing pressure above 1.5 GPa. The rapid increase of the conductivity of NaCl solution may play an important role in many geological processes (such as the genesis of ore deposits under hydrothermal condition) and other fields. Project supported by the National Natural Science Foundation of China.     

Dehydration melting of solid amphibolite at 2.0 GPa: Effects of time and temperature

The dehydration melting of the natural rock at high pressure is important to investigating the magma formation in the earth’s interior. Since the 1970s, a lot of geological scientists have paid more atten- tion to the dehydration melting of the natural rock[1―5]. Previous experiments of dehydration melting and observations of fieldwork argued that the dehy- dration melting of the rock was probably the most important fashion for the melting of the lower crust rock[6―12]. The genesis of most …     

Volume-temperature relationship for iron at 330 GPa and the Earth's core density deficit

Estimates of core density deficit (cdd) of the Earth's outer core recently reported by Anderson and Isaak [Another look at the core density deficit of Earth's outer core, Phys. Earth Planet Int. 131 (2002) 19-27] are questionable in view of the serious errors in the pressure-volume and bulk modulus data due to an inadequacy in the calibration process used by Mao et al. [Static compression of iron to 300 GPa and Fe_0.8Ni_0.2 alloy to 200 GPa: implications for the core, J. Geophys. Res. 94 (1990) 21737-21742]. The data used by Anderson and Isaak deviate significantly from the corresponding values derived from seismology. In the present study we have used the input data on density, isothermal bulk modulus and its pressure derivative from Stacey and Davis [High pressure equations of state with application to lower mantle and core, Phys. Earth Planet Int. 142 (2004) 137-184] which are consistent with the seismological data. Volumes of hexagonal close-packed iron have been calculated at different temperatures under isobaric conditions at P = 330 GPa, the inner core boundary (ICB) pressure using the relationship between thermal pressure and volume expansion based on the lattice potential theory originally due to Born and Huang [Dynamical Theory of Crystal Lattices, Oxford University Press, Oxford, 1954, p. 50]. The formulation for thermal pressure used by Anderson and Isaak has been modified by taking into account the variations of thermal expansivity α and isothermal bulk modulus K_T with temperature. Values of cdd are then estimated corresponding to different temperatures ranging from 4000 to 8000 K. The results for cdd at different temperatures obtained in the present study are significantly higher than those estimated by Anderson and Isaak suggesting that the cdd for the Earth's outer core is nearly 10%. The effects of nickel when an Fe-Ni alloy replaces Fe are estimated and found to be insignificant.     

Nonequilibrium effects in models of three-phase flow in porous media

In this paper we extend to three-phase flow the nonequilibrium formalism proposed by Barenblatt and co-workers for two-phase porous media flow. The underlying idea is to include nonequilibrium effects by introducing a pair of effective water and gas saturations, which are linked to the actual saturations by a local evolution equation. We illustrate and analyze how nonequilibrium effects lead to qualitative and quantitative differences in the solution of the three-phase flow equations.     

Scale dependent dynamic capillary pressure effect for two-phase flow in porous media

Causes and effects of non-uniqueness in capillary pressure and saturation (P^c–S) relationship in porous media are of considerable concern to researchers of two-phase flow. In particular, a significant amounts of discussion have been generated regarding a parameter termed as dynamic coefficient (τ) which has been proposed for inclusion in the functional dependence of P^c–S relationship to quantify dynamic P^c and its relation with time derivative of saturation. While the dependence of the coefficient on fluid and porous media properties is less controversial, its relation to domain scale appears to be dependent on artefacts of experiments, mathematical models and the intra-domain averaging techniques. In an attempt to establish the reality of the scale dependency of the τ–S relationships, we carry out a series of well-defined laboratory experiments to determine τ–S relationships using three different sizes of cylindrical porous domains of silica sand. In this paper, we present our findings on the scale dependence of τ and its relation to high viscosity ratio (μ_r) silicone oil–water system, where μ_r is defined as the viscosity of non-wetting phase over that of the wetting phase. An order of magnitude increase in the value of τ was observed across various μ_r and domain scales. Also, an order of magnitude increase in τ is observed when τ at the top and the bottom sections in a domain are compared. Viscosity ratio and domain scales are found to have similar effects on the trend in τ–S relationship. We carry out a dimensional analysis of τ which shows how different variables, e.g., dimensionless τ and dimensionless domain volume (scale), may be correlated and provides a means to determine the influences of relevant variables on τ. A scaling relationship for τ was derived from the dimensionless analysis which was then validated against independent literature data. This showed that the τ–S relationships obtained from the literature and the scaling relationship match reasonably well.     

Determination of phase transition pressures of ZnTe under quasihydrostatic conditions

Pressure behavior of ZnTe at room temperature was studied using an X-ray energy dispersive method on a DIA type cubic anvil apparatus (SAM-85) at NSLS-X17B1. By using powdered polyethylene, the sample and NaCl for a pressure scale were held under quasihydrostatic conditions, which were confirmed by X-ray diffraction method. Two high-pressure phase transitions were confirmed using X-ray powder diffraction simultaneously with electrical resistance measurements. The phase transition pressures under quasihydrostatic conditions were determined to be 9.6 GPa, at which the resistance increased, and 12.0 GPa, which was the midpoint of a large resistance decrease. Errors in the pressure determinations were estimated to be less than 0.2 GPa. These pressure values may depend on grain size and anisotropic stress effects on the calibrant. From X-ray observation of ZnTe, the bulk modulus of the zinc blende structure was calculated to beK ₀=51(3) GPa andK ₀ =3.6(0.8), and the first transition at 9.6 GPa was found to have about 9% volume change. It was consistent with an anomaly in the pressure generating curves.     

Closure of the Fe-S-Si liquid miscibility gap at high pressure

A high pressure investigation of melting relationships in the Fe-S-Si system has been conducted in a multi-anvil apparatus from 10 to 27 GPa and up to 2343 K. At 1 atm, the Fe-S-Si ternary system exhibits a vast miscibility gap [Raghavan, V., 1988. Phase diagrams of ternary iron alloys. Part 2: Ternary systems containing iron and sulphur. Indian Institute of Metals, Calcutta]. Quenched samples from experiments conducted at 10 and 12 GPa show an emulsion of immiscible liquids (an Fe-S melt and an Fe-Si melt). The liquid miscibility gap persists to at least 2343 K at 10 GPa. At 15 GPa, only one liquid is quenched, with a fine homogeneous dendritic texture. The results provide a mechanism to incorporate both S and Si as the light elements into the Earth’s core during a moderately high-pressure differentiation, consistent with geochemical models predicting up to 15 wt.% of light elements in the Earth’s core with 2-5 wt.% S and 7-10 wt.% Si. In contrast, for small planets such as Mars and Ganymede, differentiation took place within the pressure range of the miscibility gap. The composition of these cores is likely to be S-rich but Si-poor.     

中地壳断层带内微裂隙愈合与高压流体形成条件的模拟实验研究

中地壳断层带内发现的接近静岩压力的高压流体能够合理解释汶川MS8.0级地震断层的高角度逆冲滑动, 而高压流体的产生与断层带的微裂隙愈合紧密相关.利用熔融盐固体介质三轴高温高压实验系统,我们采用含水和烘干的Carrara大理岩样品开展了微裂隙愈合实验,研究中地壳断层带内高压流体的形成条件.实验分为三类:A类、A+B类和A+B+C类,其中A阶段实验在室温条件下将样品压裂,形成一系列共轭破裂面,B阶段实验在600℃、围压700 MPa和应变速率10-6s-1条件下愈合了A阶段破碎的样品,实验样品从以碎裂变形为主向以韧性变形为主转变,C阶段实验通过快速降低轴压模拟一个扩容过程,再以相同实验条件重新加载样品,通过比较实验样品强度来检验样品的愈合程度.样品显微结构和实验样品强度表明,动态重结晶作用能够愈合微裂隙和孔隙,水能促进矿物的动态重结晶作用,较高的水含量和较大的应变有利于微裂隙和孔隙的愈合,从而有利于高压流体的形成.     

Immiscible two-liquid regions in the Fe-O-S system at high pressure: Implications for planetary cores

We have determined phase relations in the Fe-O and Fe-O-S systems in the range of 15-21 GPa and 1825-2300 °C. Below the liquidus temperatures, solid FeO and metallic liquids are observed in both the Fe-O and the Fe-O-S systems. An immiscible two-liquid region exists in the Fe-O binary system in the pressure range investigated, and the immiscibility gap between Fe-rich metallic liquid and FeO-rich ionic liquid does not greatly change with either pressure or temperature. On the other hand, an immiscible two-liquid region in the Fe-O-S ternary system narrows significantly with increasing pressure at constant temperature and vice versa, and it almost disappears at 21 GPa, and 2300 °C. Immiscible two-liquid regions are thus not expected to exist in the Fe-O-S system in the Earth's core, suggesting that both oxygen and sulfur can be incorporated into the core. Our results are consistent with a geochemical model for the core containing 5.8 wt.% oxygen and 1.9 wt.% sulfur as proposed by McDonough and Sun [McDonough, W.F., Sun, S.-S., 1995. The composition of the Earth. Chem. Geol. 120, 223-253].     

高温高压下上地幔岩石电导率实验研究

为了建立具有普遍适用性的上地幔电性结构,本文利用Kawai-1000t压机和Solartron IS-1260阻抗/增益-相位分析仪,在4.0~14.0 GPa、873~1673 K的条件下,采用交流阻抗谱法（频率范围10^-1~10⁶Hz）测量了不含水的地幔岩电导率.实验结果显示,岩石的电导率随温度升高而大幅度的增大;在较大的温度范围内岩石的导电机制发生了变化,中低温时为小极化子导电,此时激活焓为0.94 eV（±0.13）eV,激活体积为0.11（±0.92）cm³·mol^-1,高温时为和镁空穴相关的离子导电,此时激活焓为1.6~3.17 eV,激活体积为6.75（±7.43）cm³·mol^-1;本次测量的电导率比低压下岩石的电导率要高,比矿物的电导率也要高.用本次的实验结果回归计算得到Fennoscandian地区的上地幔的一维电导率剖面,发现200 km以上本次实验计算的结果和大地电磁测深的电导率剖面吻合的比较好,在200 km以下本次实验得到的要比野外测量的电导率稍稍高一点,可能是因为实验过程中没有完全避免水的影响.本次的实验结果比用有效均匀介质方法计算得到的pyrolite矿物模型的电导率要高出两个数量级,这样的结果显示只用一种矿物的电导率或是几种矿物理论计算的结果有一定的不合理性.

     

高温高压下麻粒岩电导率研究

麻粒岩是下地壳的重要组成物质,模拟其在下地壳温压条件下的电导率对于认识下地壳的电导率分布具有十分重要的意义.本文用交流法在1 GPa下,373~1002 K范围内研究了麻粒岩样品的复阻抗,并且对测量结果进行了阻抗谱分析.研究结果表明,复阻抗对频率具有依赖性,且随温度的升高,复阻抗的实部、虚部均变小.在实验给出的温度范围内,电导率结果符合Arrhenius关系式.当温度在373~663 K范围内时,实验所获得的激化焓为0.31 eV,表明样品的电导率由低能带杂质离子所控制;当温度在673~1002 K范围内时,激化焓为0.67 eV,此时可能为小极化子导电.将所得电导率结果与西南峨边—马边地区以及华北应县―商河地区的大地电磁结果进行了对比,发现在所模拟的下地壳温压范围内,实验室测得的电导率位于野外MT数据范围内.     

Thermodynamically Constrained Averaging Theory Approach for Modeling Flow and Transport Phenomena in Porous Medium Systems: 8. Interface and Common Curve Dynamics

This work is the eighth in a series that develops the fundamental aspects of the thermodynamically constrained averaging theory (TCAT) that allows for a systematic increase in the scale at which multiphase transport phenomena is modeled in porous medium systems. In these systems, the explicit locations of interfaces between phases and common curves, where three or more interfaces meet, are not considered at scales above the microscale. Rather, the densities of these quantities arise as areas per volume or length per volume. Modeling of the dynamics of these measures is an important challenge for robust models of flow and transport phenomena in porous medium systems, as the extent of these regions can have important implications for mass, momentum, and energy transport between and among phases, and formulation of a capillary pressure relation with minimal hysteresis. These densities do not exist at the microscale, where the interfaces and common curves correspond to particular locations. Therefore, it is necessary for a well-developed macroscale theory to provide evolution equations that describe the dynamics of interface and common curve densities. Here we point out the challenges and pitfalls in producing such evolution equations, develop a set of such equations based on averaging theorems, and identify the terms that require particular attention in experimental and computational efforts to parameterize the equations. We use the evolution equations developed to specify a closed two-fluid-phase flow model.     

高温高压下地幔岩和苦橄质榴辉岩的电导率实验

为了探讨地幔岩模型和苦橄质榴辉岩模型在上地幔存在的合理性,建立上地幔的电性结构,本文利用YJ-3000t紧装式六面顶压机和Solartron IS-1260阻抗/增益-相位分析仪,在1.0~4.0 GPa、700~1150℃的条件下,采用交流阻抗谱法(频率范围10-1~106 Hz)分别测量了地幔岩和苦橄质榴辉岩的电导率.实验结果表明:随着温度的升高,地幔岩和苦橄质榴辉岩的电导率大幅增加;随着压力的增大,地幔岩的电导率略有增加,活化体积ΔV为-4.73 cm3·mol-1,而苦橄质榴辉岩的电导率几乎没有变化,活化体积ΔV为-0.11 cm3·mol-1;在电性方面,用苦橄质榴辉岩来表示深部的物质较为合理,地幔岩解释浅部可能更恰当,但浅部物质的分布不均匀,电导率随深度的变化主要受控于温度的影响,其次才是成分.