Calculated Phase Relations in the System Na2O-CaO-K2O-FeO-MgO-Al2O3-SiO2-H2O with Applications to UHP Eclogites and Whiteschists

Pressure–temperature grids in the system Na₂O–CaO–K₂O–FeO–MgO–Al₂O₃–SiO₂–H₂O and its subsystems have been calculatedin the range 15–45 kbar and 550–900°C, usingan internally consistent thermodynamic dataset and new thermodynamicmodels for amphibole, white mica, and clinopyroxene, with thesoftware THERMOCALC. Minerals considered for the grids includegarnet, omphacite, diopside, jadeite, hornblende, actinolite,glaucophane, zoisite, lawsonite, kyanite, coesite, quartz, talc,muscovite, paragonite, biotite, chlorite, and plagioclase. Compatibilitydiagrams are used to illustrate the phase relationships in thegrids. Coesite-bearing eclogites and a whiteschist from Chinaare used to demonstrate the ability of pseudosections to modelphase relationships in natural ultrahigh-pressure metamorphicrocks. Under water-saturated conditions, chlorite-bearing assemblagesin Mg- and Al-rich eclogites are stable at lower temperaturesthan in Fe-rich eclogites. The relative temperature stabilityof the three amphiboles is hornblende > actinolite > glaucophane(amphibole names used sensu lato). Talc-bearing assemblagesare stable only at low temperature and high pressure in Mg-and Al-rich eclogites. For most eclogite compositions, talccoexists with lawsonite, but not zoisite, in the stability fieldof coesite. Water content contouring of pressure–temperaturepseudosections, along with appropriate geotherms, provides newconstraints concerning dehydration of such rocks in subductingslabs. Chlorite and lawsonite are two important H₂O-carriersin subducting slabs. Depending on bulk composition and pressure–temperaturepath, amphibole may or may not be a major H₂O-carrier to depth.In most cases, dehydration to make ultrahigh-pressure eclogitestakes place gradually, with H₂O content controlled by divariantor higher variance assemblages. Therefore, fluid fluxes in subductionzones are likely to be continuous, with the rate of dehydrationchanging with changing pressure and temperature. Further, eclogitesof different bulk compositions dehydrate differently. Dehydrationof Fe-rich eclogite is nearly complete at relatively shallowdepth, whereas Mg- and Al-rich eclogites dehydrate continuouslydown to greater depth. KEY WORDS: dehydration; eclogites; phase relations; THERMOCALC; UHP metamorphism; whiteschists     

Calculated Phase Relations in the System NCKFMASH (Na2O-CaO-K2O-FeO-MgO-Al2O3-SiO2-H2O) for High-Pressure Metapelites

Petrogenetic grids in the system NCKFMASH (Na₂O–CaO–K₂O–FeO–MgO–Al₂O₃–SiO₂–H₂O)and the subsystems NCKMASH and NCKFASH calculated with the softwareTHERMOCALC 3.1 are presented for the P–T range 7–30kbar and 450–680°C, for assemblages involving garnet,chloritoid, biotite, carpholite, talc, chlorite, kyanite, staurolite,paragonite, glaucophane, jadeite, omphacite, diopsidic pyroxene,plagioclase, zoisite and lawsonite, with phengite, quartz/coesiteand H₂O in excess. These grids, together with calculated compatibilitydiagrams and P–T and T–X_Ca and P–X_Ca pseudosectionsfor different bulk-rock compositions, show that incorporationof Ca into the NKFMASH system leads to many of the NKFMASH invariantequilibria moving to lower pressure and/or lower temperature,which results, in most cases, in the stability of jadeite andgarnet being enlarged, but in the reduction of stability ofglaucophane, plagioclase and AFM phases. The effect of Ca onthe stability of paragonite is dependent on mineral assemblageat different P–T conditions. The calculated NCKFMASH diagramsare powerful in delineating the phase equilibria and P–Tconditions of natural pelitic assemblages. Moreover, contoursof the calculated phengite Si isopleths in P–T and P–X_Capseudosections confirm that phengite barometry in NCKFMASH isstrongly dependent on mineral assemblage. KEY WORDS: phase relations; metapelites; NCKFMASH; THERMOCALC; phengite geobarometry     

地震信号的相位恢复及误差分析

常见的几种地震子波有:零相位子波、最小相位子波、最大相位子波和混合相位子波。以不同类型的信号为激发子波,对不同地层模型合成地震道。仅用原始合成地震道记录的振幅谱来恢复原始记录,讨论不同方法下,由于激发子波的不同,对比分析原始记录与恢复记录的差异,及其恢复效果。     

微量相X射线衍射定量分析的一种简便方法

阐述了关于微量相X射线衍射定量分析的一种简便方法———近似线性关系法,它适用于二相体系(或基体吸收系数基本恒定的多相体系)材料中微量相的定量分析。该方法不需要在待测样品中加入参考物相,分析时仅需对一条衍射线进行测定,因而分析时间很短,其测定结果与实际含量的相对误差为-10.2%～ 21.4%,与外标法和内标法的测量结果的精度相当。     

三种类型烃-烃不混溶包裹体组合的特征和热力学条件的计算

在自然界广泛分布着烃-烃不混溶体系中捕获的流体包裹体,由于这些包裹体具有复杂的组成和相态,因此不混溶包裹体组合的判别和热力学参数的计算常常难以进行。根据烃-烃不混溶体系中两个端员组分流体包裹体室温下的相态特征和在温度-压力平面图上等容线交点显示的位置,划分成三种类型流体包裹体组合,本介绍了三种类型流体包裹体组合特征,叙述了不混溶烃-烃包裹体组合的测定和判别方法,并且阐述均一化包裹体相态方程和气-液平衡常数原理和方法与此同时列举了自然界简单的三种类型不混溶烃-烃包裹体组合的测定、判别和计算的几个实例,利用相态方程和气-液平衡常数,不但精确地计算出包裹体均一压力,并且精确地计算出流体密度和体积等热力学参数。最后,利用均一成气相和液相的两种包裹体在 p-T 平面图上等容线交点同样计算出流体包裹体组合的捕获温度和压力。     

水气两相流系统K-S-P模型参数反演的最优估计

应用渗流力学和参数反演辨识最优化理论,建立了水气两相流控制方程,给出了控制两相流系统K-S-P模型参数之间的本构关系,并对其中的相关参数进行数值反演辨识,得出待求参数的最优估计值。对不同质地土壤的实验测试值进行对比分析,其结果吻合较好,验证了模型的适定性和可靠性,解决了实验数据处理以及结果对初值的敏感性问题,并且该数值方法可通过多组初值参数进行数值实验,结合实验结果来满足实际工程的需求。这不仅对非饱和带内污染物传输、气藏工程试井计算、填埋气体的释放以及挥发性有机污染物的去污分析中两相流数值计算提供基础,而且对K-S-P本构关系之间参数的确定提供了可靠的理论依据。     

单频GPS整周模糊度动态快速求解的研究

提出了一种动态快速求解整周模糊度的方法,即先对系数阵进行QR分解,然后通过矩阵变换使模糊度参数和位置参数分离,从而降低矩阵的维数,满足实时动态求解的要求,最后应用LAMBDA方法搜索模糊度。为验证该算法,用单频GPS接收机进行了静态定位和动态定位两种试验。结果表明,静态定位误差小于1cm,动态定位误差小于3cm。由此可见,该方法能够为动态用户快速而精确地实施GPS单频动态定位。     

少量地面控制点的遥感图像精纠正与定位技术探讨

本文主要讨论了利用空间投影一空间斜墨卡托（SOM）投影进行卫星遥感图像精纠正与定位，解决缺少控制点地区卫星遥感图像精纠正与定位问题，并根据正反算公式讨论了卫星遥感图像纠正的空间投影数学基础。     

Numerical evaluation of the damping‐solvent extraction method in the frequency domain

The damping‐solvent extraction method for the analysis of unbounded visco‐elastic media is evaluated numerically in the frequency domain in order to investigate the influence of the computational parameters—domain size, amount of artificial damping, and mesh density—on the accuracy of results. An analytical estimate of this influence is presented, and specific questions regarding the influence of the parameters on the results are answered using the analytical estimate and numerical results for two classical problems: the rigid strip and rigid disc footings on a visco‐elastic half‐space with constant hysteretic material damping. As the domain size is increased, the results become more accurate only at lower frequencies, but are essentially unaffected at higher frequencies. Choosing the domain size to ensure that the static stiffness is computed accurately leads to an unnecessarily large domain for analysis at higher frequencies. The results improve by increasing artificial damping but at a slower rate as the total (material plus artificial) damping ratio ζ_t gets closer to 0.866. However, the results do not deteriorate significantly for the larger amounts of artificial damping, suggesting that ζ_t≈0.6 is appropriate; a larger value is not likely to influence the accuracy of results. Presented results do not support the earlier suggestion that similar accuracy can be achieved by a large bounded domain with small damping or by a small domain with larger damping. Copyright © 2002 John Wiley & Sons, Ltd.     

裂隙型单斜介质中弹性系数的计算及波的传播特性研究

根据Hudson等关于裂隙介质弹性系数计算的扰动理论及Bond变换矩阵原理，给出了各向同性介质中含多组垂直裂隙时等效弹性系数的计算方法。计算了含两组斜交的垂直裂隙形成的单斜各向异性介质中的等效弹性系数，并根据Christoffel方程推导、得出水平面内平面波传播的相速度和群速度随方位变化的特性。