改进时间分数阶模型模拟非Fick溶质运移

通过将经典时间分数阶对流-弥散方程的等待时间分布函数的尾部修改为指数型,推导出了改进时间分数阶对流-弥散方程,并提出有效的时空算子分裂数值求解方法。对两个理想算例和一个实际算例进行计算,结果表明,改进的时间分数阶对流-弥散方程继承了时间分数阶对流-弥散方程能模拟穿透曲线幂率型拖尾分布的优点,还可模拟穿透曲线尾部由幂率型转换到指数型的过程;特征时间λ、分数阶指数γ和两相容量比例系数β共同决定了运移行为。改进的新模型可以区分非均质介质中流动相和非流动相中的溶质浓度, 更细微地模拟非Fick溶质运移行为。     

Permian High Ba-Sr Granitoids: Geochemistry, Age and Tectonic Implications of Erlangshan Pluton, Urad Zhongqi, Inner Mongolia

Abstract: Erlangshan Pluton from Urad Zhongqi, central Inner Mongolia, is located in the middle segment of the northern margin of the North China Plate. The rocks consist mainly of diorites with gneissic structure. Petrochemical characteristics reveal that the diorites belong to metaluminous, high-potassium calc-alkaline series, with chemical signatures of I-type granites. They are characterized by low SiO2 contents (56.63%–58.53%) and A/CNK (0.90–0.96), high Al2O3 contents (17.30%–17.96%) and Na2O/K2O ratios (1.20–1.70), enrichment in large ion lithophile elements (LILE, e.g., Ba=556–915 ppm, Sr=463–595 ppm), and relative depletion in high field strength elements (HFSE, e.g., Nb, Ta, Ti) in primitive mantle-normalized spidergram, and right-declined rare earth element patterns with slightly negative Eu anomalies (d?Eu=0.72–0.90). They have Sr/Y ratios (20–25) evidently less than Kebu Pluton (49–75) to its east. Sensitive high resolution ion micro-probe U-Pb zircon dating of the diorites has yielded an intrusive age of 270±8?Ma. This leads us to conclude that Erlangshan diorites were formed by mixing between the middle or lower crustal-derived magma and minor mantle-derived mafic magma, followed by fractional crystallization, which was trigged by crustal extension and fault activity in post-collisional setting.     

A garnet population in Yellowknife schist, Canada

Abstract Data are presented on a garnet population in a specimen of garnet-biotite-plagioclase-quartz schist from the cordierite zone of an Archaean thermal dome in the Southern Slave Province of the Canadian Shield. Garnet crystals are bounded by planar dodecahedral faces and by trapezohedral faces which on the 10-μm scale are corrugated. Crystal distribution, as revealed by dissection of a small cubic volume of rock, is random. The size distribution is normal, with a mean diameter of 0.81 mm and a standard deviation of 0.32 mm. In the largest crystal of the population (mean radius 0.83 mm), [Mn] = 100 Mn/(Fe + Mg + Mn + Ca) decreases from 14.5 at the centre to 7.5 and then increases in the outer margin to 8.5; [Fe] increases continuously from 67 at the centre to 77 at the surface; [Mg] increases from 12.5 to 13.5 and then falls sharply to 11; [Ca] remains unchanged at 4.0 and then drops to 3.3. Progressively smaller crystals have progressively lower [Mn] and higher [Fe] concentrations at their centres, while all crystals have the same margin composition. Growth vectors extending from given concentration contours to crystal surfaces are of equal length regardless of the size of the crystal in which the vector is located. A garnet-forming model is presented in which reaction was initiated by a rise in temperature. Nucleation sites were randomly selected. The nucleation rate increased with time and then declined. Crystal faces advanced at a constant linear rate, which implies an increase in volume proportional to surface area. Initially, the composition of garnet deposited on crystal surfaces was determined by van Laar equations of equilibrium, which demanded the withdrawal of Mn and Fe from within chlorite crystals. This transfer reaction was then accompanied by an ion exchange reaction which moved Mn and Fe to garnet surfaces from biotite, in exchange for Mg. The exchange reaction provides an explanation for the high overall concentration of Mn and Fe in garnet and for the observed Mn and Mg reversals in the margins of crystals. The increase of garnet volume in the garnet population is found to be parabolic, i.e. Vαα⁵.     

二维分形技术在地震勘探构造解释中的应用

作者在文中探讨了二维分形布朗运动在地震勘探构造解释中的应用,并针对地震资料,提出了具体实用的算法,并以此算用Ｃ＾＋＋语言编制了相应处理软件,对实测地震资料进行了处理,取得了明显的效果。     

An iterative pressure‐stabilized fractional step algorithm in saturated soil dynamics

This paper presents an iterative, incremental pressure‐stabilized fractional step algorithm for coupled hydro‐mechanical problems with mixed formulations of the displacement–pressure ( u –p) model in saturated soil dynamics that allows the use of finite elements with equal low order of interpolation approximation of u and p. In comparison with the original fractional step algorithm, the distinct features of the proposed algorithm lie in its enhanced stability owing to the introduction of both an iteration procedure and a finite increment calculus (FIC) process into the algorithm. The introduction of the iterative procedure makes the velocity term satisfy the momentum conservation equation in an implicit sense and allows much larger time step sizes to be used than those limited in existing explicit and semi‐implicit versions of the algorithm. The introduction of the FIC process removes the dependence of the stability of the proposed algorithm on the time step size, as a result it allows to using the incremental version of the algorithm and evades the minimum time step size requirement presented in the existing versions of the fractional step algorithm that restricts the application of the algorithm to saturated soil dynamics problems with high frequencies. Numerical experiments demonstrate the effectiveness and improved performance of the proposed iterative pressure‐stabilized fractional step algorithm. Copyright © 2009 John Wiley & Sons, Ltd.     

平衡剖面研究中的新方法

随着现代化计算手段的运用的新的地学思维方法的影响,在平衡剖面研究中引入计算机模拟技术,分数维方法,三维模式等新方法,使平衡剖面技术具有了更加广阔的应用前景。另外,本文还对平衡剖面技术的局限性和发展趋势进行了探讨。     

Phase Equilibria Constraints on the Chemical and Physical Evolution of the Campanian Ignimbrite

The Campanian Ignimbrite is a > 200 km³ trachyte–phonolitepyroclastic deposit that erupted at 39·3 ± 0·1ka within the Campi Flegrei west of Naples, Italy. Here we testthe hypothesis that Campanian Ignimbrite magma was derived byisobaric crystal fractionation of a parental basaltic trachyandesiticmelt that reacted and came into local equilibrium with smallamounts (5–10 wt%) of crustal rock (skarns and foid-syenites)during crystallization. Comparison of observed crystal and magmacompositions with results of phase equilibria assimilation–fractionationsimulations (MELTS) is generally very good. Oxygen fugacitywas approximately buffered along QFM + 1 (where QFM is the quartz–fayalite–magnetitebuffer) during isobaric fractionation at 0·15 GPa ( 6km depth). The parental melt, reconstructed from melt inclusionand host clinopyroxene compositions, is found to be basaltictrachyandesite liquid (51·1 wt% SiO₂, 9·3 wt%MgO, 3 wt% H₂O). A significant feature of phase equilibria simulationsis the existence of a pseudo-invariant temperature, 883 °C,at which the fraction of melt remaining in the system decreasesabruptly from 0·5 to < 0·1. Crystallizationat the pseudo-invariant point leads to abrupt changes in thecomposition, properties (density, dissolved water content),and physical state (viscosity, volume fraction fluid) of meltand magma. A dramatic decrease in melt viscosity (from 1700Pa s to 200 Pa s), coupled with a change in the volume fractionof water in magma (from 0·1 to 0·8) and a dramaticdecrease in melt and magma density acted as a destabilizingeruption trigger. Thermal models suggest a timescale of 200kyr from the beginning of fractionation until eruption, leadingto an apparent rate of evolved magma generation of about 10^–3km³/year. In situ crystallization and crystal settling in density-stratifiedregions, as well as in convectively mixed, less evolved subjacentmagma, operate rapidly enough to match this apparent volumetricrate of evolved magma production. KEY WORDS: assimilation; Campanian Ignimbrite; fractional crystallization; magma dynamics; phase equilibria     

西藏驱龙斑岩铜矿床岩浆演化过程中的Cu同位素分馏

随着同位素测试技术的发展,高温下重同位素分馏逐渐运用到地球化学领域.为了探索岩浆演化过程中是否会发生Cu同位素分馏,文章对西藏驱龙斑岩铜矿床同来源、不同演化阶段的两套岩浆岩(闪长岩包体和花岗闪长岩)进行了Cu同位素测试.测试结果显示,闪长岩包体的δ65Cu值集中在+0.08‰-+0.35‰之间,平均值为+0.25‰,花岗闪长岩的δ65Cu值为+0.42‰～+0.87‰,平均值为+0.60‰.△65Cu花岗闪长岩-闪长岩包体≈+0.4‰,并且δ65Cu值与样品w(SiO)变化存在一定的正相关性.结合地质学及同位素方面考虑,花岗闪长岩和闪长岩包体的Cu同位素组成差异可能是岩浆演化过程中发生了Cu同位素分馏所致.在驱龙矿区中新世岩浆演化过程中,随着岩浆去气作用,63Cu随HS-、Cl-等挥发分优先进入气相,导致残留岩浆熔体相对富集65Cu.此外,两者的Cu同位素组成差异也可能与岩浆演化过程中斜长石、角闪石、黑云母等矿物的不断结晶分离有关,矿物结晶分离时,基性矿物富集63Cu,而残余熔体则富集65Cu.     

岩浆分离结晶作用过程微量元素反演

介绍岩浆分离结晶作用中微量元素系统反演的基本原理和方法,反演计算中极小化问题的求解,反演所需要的优选模式参数估计方法,进行反演的计算机程序,利用模拟数据和实测数据检验了方法和程序的有效性.