温度作用下岩石的热黏弹塑性模型研究

在西原体模型的基础上,通过引入弹簧的热膨胀系数,将温度效应引入到模型中,推导了单轴应力状态下微分形式的热粘塑性本构方程。通过蠕变曲线和卸载曲线的分析,考虑了温度效应的西原体模型当应力水平较低时为稳定蠕变,存在瞬时弹性、弹性后效和由温度引起的粘性流动;当应力水平较高时为不稳定蠕变,存在瞬时弹性、弹性后效和由温度和应力共同引起的粘性流动,该模型较全面地反映了岩石在温度作用下的蠕变性质。     

一种基于修正西原模型的冻土蠕变本构关系

冻土在高应力水平下通常会表现出加速蠕变特性。经典西原（Nishihara）模型可以反映衰减蠕变和稳态蠕变,但不适用于描述加速蠕变。通过对Nishihara模型中黏性元件进行修正,将黏弹性元件中的定常的黏滞系数修正为时间的函数,并在黏塑性元件中采用非定常、非牛顿黏性元件,得到的蠕变柔量可以考虑应力水平和时间的影响,使其能够反映其加速蠕变的过程。用提出的模型分别对一维和三维应力、不同应力水平下的冻土蠕变试验数据进行了拟合,拟合结果能够反映冻土从衰减蠕变和稳定蠕变状态进入加速蠕变的过程,说明了该模型对于描述冻土在不同应力下不同蠕变状态的适用性。     

VES-18清洁压裂液性能试验与应用

用阳离子季胺盐和有机盐等复配制备VES-18清洁压裂液,对其主要性能进行了研究,并进行了现场应用试验。结果表明:VES-18具有很好的耐低温性和时间稳定性;在常温、170s-1剪切速率下,表观粘度为11mPa·s;在200~1000s-1区间内,随着剪切速率的增大,粘度缓慢下降,且下行粘度高于上行粘度;在应力扫描范围内,粘弹曲线中弹性模量一直大于粘性模量,属于弹性流体;模拟动态沉降时间大于粘度高的胍胶基液和丙三醇,体现了很好的携砂性能;少量的原油和自制破胶剂PJJ都可使VES-18在4h内完全破胶,破胶液粘度小于3mPa·s,破胶液表面张力较低,无残渣,破胶效果好,与地层配伍性良好。      

动扭剪试验中砂土液化后流动特性分析

土体在初始液化后仍然可能承受动荷载作用发生大变形。引入流体力学中的剪应变率和表观黏度的概念,对振动扭剪试验中饱和砂土液化后的流动特性进行了分析。分析中将砂土的状态分为0有效应力状态和非0有效应力状态。结果表明,砂土在0有效应力状态下表现出与静扭剪试验类似的“剪切稀化”非牛顿流体的特征,表观黏度随着剪应变率的增大而减小。加载周数对“剪切稀化”状态下的剪应变率幅值有影响,随着加载周数的增加,剪应变率幅值逐渐增大,而对流动曲线的形状没有影响。在非0有效应力状态下,砂土的表观黏度随应变的增大而增大,随孔压比的减小而增大,且所有的试验得到的表观黏度与孔压比具有一致关系。     

An improved hardening-damage creep model of lean coal: a theoretical and experimental study

Creep is an important rheological behavior of coal, impacting pillar stability, gateway support performance, and construction costs in underground mining. It is of great significance to study the creep mechanism of coal to prevent creep failure. In this paper, uniaxial compressive creep tests of raw lean coal under multiple stages are conducted, the increasing instantaneous elastic modulus and decreasing viscosity coefficient are calculated to analyze the hardening-damage creep mechanism of coal, and some new conclusions are as follows. When the stress is low, the hardening effect exists, and coal strain shows an instantaneous response. At high stress, both the damage effect and hardening effect work and coal strain shows time-dependent transient creep and steady creep. During multistage creep, coal is first hardened, then weakened, and finally fails due to well-developed cracks. Based on the strain evolution laws of lean coal at different stress levels, an elastic viscous-plastic model is established, and the hardening function and damage function are introduced into the model to obtain an improved hardening-damage creep model. This model is used to fit the creep data, and the model curves match the testing data very well, showing much better accuracy than other models. This model is simple with clear physical meanings of the elements and can reflect the nonlinear hardening-damage creep of raw lean coal under uniaxial compression very well.     

Computation of the filtration laws through porous media for a non-Newtonian fluid obeying the power law

We consider a stationary flow of an incompressible non-Newtonian flow through a porous medium, induced by an injection velocity when inertial effects are negligible. At the pore scale, the governing equations are based on a nonlinear relation between the stress and the rate of deformation. In such a situation, the limit problem obtained when the pore size tends to zero, is called the homogenized problem that leads to the filtration law. This filtration law is given by a non-linear system coupling a local problem on a typical cell of the porous medium to a global problem at the scale of the whole porous medium. We propose, in this work, a numerical method to solve this homogenized problem and apply this method when the velocity dependent viscosity is given by the power law. Finally, we propose some numerical experiments to illustrate our approach.     

A new approach to the equation of state of silicate melts: An application of the theory of hard sphere mixtures

A comparison of compressional properties of silicate solids, glasses, and liquids reveals the following fundamental differences: (1) Liquids have much smaller bulk moduli than solids and glasses and the bulk moduli of various silicate melts have a narrow range of values; (2) Liquids do not follow the Birch’s law of corresponding state as opposed to solids and glasses; (3) The Grüneisen parameter increases with increasing pressure for liquids but decreases for solids; (4) The radial distribution functions of liquids show that the interatomic distances in liquids do not change upon compression as much as solids do. The last observation indicates that the compression of silicate melts occurs mostly through the geometrical arrangement of various units whose sizes do not change much with compression, i.e., the entropic mechanism of compression plays a dominant role over the internal energy contribution. All of the other three observations listed above can be explained by this point of view. In order to account for the role of the entropic contribution, we propose a new equation of state for multi-component silicate melts based on the hard sphere mixture model of a liquid. We assign a hard sphere for each cation species that moves in the liquid freely except for the volume occupied by other spheres. The geometrical arrangement of these spheres gives the entropic contribution to compression, while the Columbic attraction between all ions provides the internal energy contribution to compression. We calibrate the equation of state using the experimental data on room-pressure density and room-pressure bulk modulus of liquids. The effective size of a hard sphere for each component in silicate melts is determined. The temperature and volume dependencies of sphere diameters are also included in the model in order to explain the experimental data especially the melt density data at high pressures. All compressional properties of a silicate melt can be calculated using the calibrated sphere diameters. This equation of state provides a unified explanation for most of compressional behaviors of silicate melts and the experimental observations cited above including the uniformly small bulk moduli of silicate melts as well as the pressure dependence of Grüneisen parameters. With additional data to better constrain the key parameters, this equation of state will serve as a first step toward the unified equation of state for silicate melts.     

硅酸盐熔体包裹体均匀化过程的动力学性质

硅酸盐熔体包裹体的研究,和以模拟实验为主的实验岩石学,实验矿物学一样,已经成为现代岩浆岩石学,尤其是理论岩浆岩石学研究的一个必不可少的组成部分。均匀化测温研究,是硅酸盐熔体包裹体研究的一个重要组成部分。其方法是对于呈封闭系统的硅酸盐熔体包裹体进行加热,在反方向上,再造这些包裹体从结晶到熔融的所有阶段,直至其完全均匀化,也就是恢复到它们被生长着的矿物包含、捕获时的状态,借以测量矿物近液相线或固相线的结晶温度。     

Triaxial tests on the fluidic behavior of post-liquefaction sand

Liquefaction-induced ground deformation is a major cause of structural damage during earthquakes. However, a better understanding of seismic liquefaction is needed to improve earthquake hazard analyses and mitigate structural damage. In this paper, a dynamic triaxial test apparatus was employed to investigate the fluidic characteristics of post-liquefaction sand. The specimens were vibrated to the point of liquefaction by dynamic loading, and then the liquefied sand was further sheared by triaxial compression in an undrained manner. It was found that a non-Newtonian fluid model can accurately describe the shear stress and the shear strain rate of post-liquefaction sand during undrained triaxial compression. The apparent viscosity, a major parameter in a constitutive model of a non-Newtonian fluid, decreases with an increase in the shear strain rate.     

Manganese partitioning during hydrous melting of peridotite

Manganese contents and the iron/manganese ratio of igneous rocks have been used as a method of probing the heterogeneity in the Earth’s mantle during melting of peridotite and pyroxenite lithologies. Most previous work has assumed that changes in these parameters require differences in either source lithology or composition based on experiments indicating that manganese is slightly incompatible during melting and that the iron/manganese ratio is fixed by the presence of olivine. However, the presence of volatiles in the mantle drives melting at lower temperatures and with different compositions than in volatile-free systems, and thus the partitioning of Fe and Mn may in fact vary. We have produced silicate liquids in equilibrium with a peridotite assemblage under hydrous conditions at 3 GPa that show that Mn can also be unexpectedly compatible in garnet at 1375 °C and that Mn partitioning between solids and liquids can be strongly affected by temperature and liquid composition. The compatibility of Mn in garnet provides a mechanism for large variations of Mn contents and the Fe/Mn ratio in silicate melts that solely involves melting of mantle peridotite with only small compositional changes. Correlations between Mn variations and other indices indicative of melting in the presence of garnet may provide a means of more completely understanding the role of garnet at high pressures in peridotite melting.     

高温硅酸盐熔体粘度与网络分数维值的相关性研究

对高温硅酸盐熔体粘度的估算一直是国际地学界热点问题之一,本文在研究了高温硅酸盐熔体网络分数维值的基础上,建立了估算熔体粘度的新模式 (简称 FD模式 ) ,阐明了熔体粘度值除了与温度成反比外,还与分子网络介观尺度的自相似比 ri 和分数维值 D相关,即与分子网络 (纳米尺度以上 )中的桥氧数 Ni成正比 (Ni∝ ) , 与单位硅氧四面体中的非桥氧数成反比.经对 4个硅酸盐系列高温熔体的粘度测定证实, FD模式的理论计算值与实测值吻合,且优于现今国际通用的 VTF模式.为探索岩浆迁移演化规律和完善新型低维材料的性能提供了新的理论依据.     

岩体的加、卸载状态与能量的分配关系

岩体是流变体,在高地应力的长时间作用下其内部非均匀应力随着时间自我松弛,逐渐趋向于静水压力状态。在这样的状态下,虽然其体积仍维持弹性,但丧失了进一步抵抗剪切变形的能力,处于流动的状态,其能量储存具有体积的特征。采用平面压缩模型,研究岩石在不同应力水平长期作用下的变形性状及能量分配关系,给出在快速卸荷条件下发生自持续动力断裂破坏条件,以期揭示深部岩体初始地应力水平（能量在体变及形变上的分配关系）对卸荷条件下岩体运动过程中表现的惯性、黏性及弹塑性变形的影响规律,为建立深部岩体加卸荷本构关系打下基础。     

A computer simulation study of natural silicate melts. Part II: High pressure properties

The thermodynamic, structural and transport properties of natural silicate melts under pressure are investigated by molecular dynamics simulation with the help of a force field recently introduced by us [Guillot B. and Sator N. (2007) A computer simulation study of natural silicate melts. Part I: low pressure properties. Geochim. Cosmochim. Acta71, 1249-1265]. It is shown that the simulation reproduces accurately the bulk moduli of a large variety of silicate liquids as evaluated from ultrasonic studies. The equations of state (EOS) of the simulated melts are in good agreement with the density data on mid-ocean ridge basalt, komatiite, peridotite and fayalite as obtained either by sink/float experiments or by shock-wave compression. From the structural point of view it is shown that the population of ^[5]Al and ^[6]Al species increases rapidly upon initial compression (0-50 kbar) whereas for Si these highly coordinated species are found in a significant abundance (>5%) only above ∼50 kbar for ^[5]Si and ∼100-150 kbar for ^[6]Si. This increase of the coordination of network formers is not the only response of the melt structure to the densification: there is also a large redistribution of the T-O-T (T = Si, Al) bond angles with the pressure and noticeably upon initial compression in rhyolitic and basaltic liquids. Furthermore, a detailed analysis of the population of bridging oxygens (BO) and nonbridging oxygens (NBO) points out that the polymerization of the melt generally increases when the pressure increases, the magnitude of this polymerization enhancement being all the more important that the melt is depolymerized at low pressure. The role of triclusters (threefold coordinated oxygens to network former cations) is particularly emphasized in acidic and basaltic liquids. The pressure-induced redistribution of the oxygen atoms through the melt structure is also stressed. Finally, the simulation predicts a nonmonotonic behavior of the diffusivity of network former ions when the pressure increases, a feature with depends on the melt composition. This could have a counterpart in the electrical conductivity at sufficiently high temperature when the viscosity of the liquid is low.     

北京西山膨缩石香肠流变性质构造研究初探

本文探索了利用膨缩石香肠构造定量估测古岩石流变参数的方法,获得北京西山云母大理岩中的方解石石英脉的应力指数在3.06-7.69之间,与云母大理岩的粘度比在8.74-20.27之间;板岩中的石英脉的应力指数在2.57-4.75之间,与板岩的粘度比在94.8-175.8之间;白云岩中硅质条带的应力指数在2.25-3.93之间,与白云岩的粘度比在2.38-4.26之间.      

硅酸盐熔体结构的研究进展和问题

概括地介绍了硅酸盐熔体结构的三种研究途径和方法,从硅酸盐熔体的分子聚合结构单元测试、阳离子和挥发份的结构作用、物理化学性质测试及量子化学研究等五个方面,阐明了硅酸直熔体结构的主要成就,指出了该研究领域中存在的主要问题。     

Iron-bearing silicate melts: Relations between pressure and redox equilibria

Mossbauer spectroscopy has been used to determine the redox equilibria of iron and structure of quenched melts on the composition join Na₂Si₂O₅-Fe₂O₃ to 40 kbar pressure at 1400° C. The Fe³⁺/ΣFe decreases with increasing pressure. The ferric iron appears to undergo a gradual coordination transformation from a network-former at 1 bar to a network-modifier at higher (≧10 kbar) pressure. Ferrous iron is a network-modifier in all quenched melts. Reduction of Fe³⁺ to Fe²⁺ and coordination transformation of remaining Fe³⁺ result in depolymerization of the silicate melts (the ratio of nonbridging oxygens per tetrahedral cations, NBO/T, increases). It is suggested that this pressure-induced depolymerization of iron-bearing silicate liquids results in increasing NBO/T of the liquidus minerals. Furthermore, this depolymerization results in a more rapid pressure-induced decrease in viscosity and activation energy of viscous flow of iron-bearing silicate melts than would be expected for iron-free silicate melts with similar NBO/T.     

Silicate Liquid Immiscibility in Tholeiitic Basalts

Immiscible silicate liquids, preserved as fresh glasses, constituteat least 32 per cent of a tholeiitic basalt in the Triassic-Jurassicbasin of Southbury, Connecticut. This apparently uncommon abundanceof immiscible liquids cannot be attributed to an unusual magmacomposition, for the rock is similar to many quartz-normativetholeiites associated with other basins of similar age in easternNorth America. Nor can peculiar conditions of fractional crystallizationbe involved, for experiments carried out at one atmosphere totalpressure reveal that two silicate liquids coexist over a widetemperature interval during crystallization, and that this intervalis increased by raising the oxygen fugacity. Experiments onthe standard diabase, W-1, which differs only slightly in compositionfrom the Connecticut basalt, indicates the presence of immiscibleliquids during the crystallization of this magma compositionas well. It is concluded that liquid immiscibility is to beexpected during the crystallization of many tholeiitic magmas,and at a sufficiently early stage of solidification to playa role in igneous differentiation.     

A back-of-the-envelope approach to boudinage mechanics

Smith has shown from linear hydrodynamic stability that pinch-and-swell boudinage is possible only in a non-Newtonian material. By assuming uniaxial stress and strain and large viscosity contrast between the competent layer and the incompetent matrix, Smith's result can be derived in a very simple way. The one-dimensional force balance also predicts growth rates, although not dominant wavelengths. However, we question the ability of even linear stability analysis to accurately predict boudin widths.     

The coordination number of ferrous ions in silicate glasses

Mössbauer and optical absorption measurements on a number of silicate glasses show that most of the ferrous ions in these glasses are distributed over two types of octahedral sites. Less than 0.5% of all ferrous ions occupy tetrahedral sites. It is shown that the Bernal liquid model of random close packing does not apply to silicate glasses. Attention is drawn to the fact that measurements on silicate glasses cannot provide reliable information about coordination numbers in silicate liquids.     

Decoupling of paired elements,crossover REE patterns,and mirrored spider diagrams: Fingerprinting liquid immiscibility in the Tapira alkaline–carbonatite complex,SE Brazil

Tapira is an alkaline silicate–carbonatite complex belonging to the kamafugite-carbonatite association in the Late-Cretaceous Alto Paranaíba Igneous Province (APIP). It is dominated by coarse-grained plutonic rocks (bebedourite – a phlogopite-, apatite-, and perovskite-rich clinopyroxenite – with subordinated dunites, wehrlites, carbonatites and phoscorites). The plutonic rocks are crosscut by fine-grained ultramafic alkaline rocks (phlogopite picrites, bebedouritic dikes) and fine-grained carbonatites. Both types of dike-rocks show petrographic evidence of the coexistence of immiscible silicate and carbonatite liquids, such as carbonate ocelli present in the silicate rocks and, more rarely, silicate ocelli within carbonatites. A detailed geochemical study of the rock types in the complex, with emphasis on the fine-grained varieties, showed that whilst some rocks may be related to each other through crystal fractionation (e.g. phlogopite picrites and bebedouritic dikes), others display anomalous trace-element behaviour that cannot be readily explained by the fractionation of a particular phase or combination of phases. We interpret such anomalous geochemical signatures as produced by silicate–carbonate liquid immiscibility, on the basis of available experimental data on partition coefficients between coexisting immiscible liquids. The immiscibility signatures comprise: (a) decoupling of geochemical pairs, such as Nb–Ta and Zr–Hf; (b) rotation of REE patterns, which cross over the patterns of the primitive liquids; and (c) matching and opposite enrichment-depletion trace elements relationships in spider diagrams of conjugate immiscible liquids. We suggest that, once established, such geochemical signatures are very difficult to erase during the subsequent petrogenetic evolution processes, which may result in superimposed conflicting signatures.