思林水电站通航建筑物高边坡稳定分析和评价

应用ANSYS结构分析软件,模拟了思林水电站通航建筑物高边坡的地形、地质条件以及开挖的过程,建立边坡开挖的有限元计算分析模型,分析边坡在自然状态、分级开挖以及锚杆支护加固等工况下的应力、变形及其整体稳定性.定义单元的局部安全系数,得到了边坡开挖加固后稳定安全系数的分布.根据计算成果,指出现设计加固措施可以保证边坡稳定,无需增加其他工程措施.     

拉拉铜矿黄铁矿微量元素地球化学特征及其成因意义

四川会理拉拉铜矿床是我国著名大型富铜矿床,针对该矿床中黄铁矿的微量元素、稀土元素地球化学分析表明:拉拉铜矿经历了早期火山喷发成岩成矿和晚期变质成岩成矿作用.条带状矿石中的黄铁矿Co/Ni比值集中于4.92～79.2之间,落入火山成因黄铁矿区,稀土元素分布具有Eu正异常和轻稀土富集的特征,反映矿床具有伴随河口群火山喷流沉积成岩过程的同生沉积成矿作用.脉状矿石中的黄铁矿Co/Ni比值集中于1.10～3.45,落在热液成因黄铁矿区,稀土元素较河口群岩石及其他典型块状硫化物矿床矿石稀土元素更加富集轻稀土元素,稀土含量变化范围更大,显著的负Eu异常,则又说明,矿床形成的主要成矿作用是伴随新元古代晋宁运动而发生的大规模的变质作用.     

Distribution Characteristics of Sulfur and the Main Harmful Trace Elements in China’s Coal

To promote the rational development and use of clean coal resources in China, data on the regional and age distribution of sulfur, arsenic and other harmful elements in Chinese coal was broadly collected, tested for content, and analyzed. Coal in northwestern China is characterized by low to extremely low levels of sulfur; the coal of the Taiyuan Formation in northern China mainly has high-sulfur content; that of the Shanxi Formation is mainly characterized by low sulfur coal; and the Late Permian coal in southern China has overall higher sulfur content; other regions have low sulfur coal. The average content of harmful trace elements in the bulk of China＇s coal is similar to the corresponding content in the coal of the North America and the rest of the world, whereas the content of various elements （Hg, Sb and Se） is different in magnitude to the corresponding percentage in the crust. The average content of the elements Cr, Se, Co, Be, U, Br in Late Permian coal in S China ranks first in the country whereas the average content of Hg and CI in the coals of Late Carboniferous to Early Permian age in N China are the highest. The average content of Mn in Early and Middle Jurassic coal is higher in NW China. The high content of harmful elements in some coal should cause particular concern both in the development and utilization of coal.     

U-Pb Age Determination for Seven Standard Zircons using Inductively Coupled Plasma–Mass Spectrometry Coupled with Frequency Quintupled Nd-YAG (λ = 213 nm) Laser Ablation System: Comparison with LA-ICP-MS Zircon Analyses with a NIST Glass Reference Material

This paper evaluates the analytical precision, accuracy and long‐term reliability of the U‐Pb age data obtained using inductively coupled plasma–mass spectrometry (ICP‐MS) with a frequency quintupled Nd‐YAG (λ = 213nm) laser ablation system. The U‐Pb age data for seven standard zircons of various ages, from 28 Ma to 2400 Ma (FCT, SL13, 91500, AS3, FC1, QGNG and PMA7) were obtained with an ablation pit size of 30 μm diameter. For ²⁰⁷Pb/²⁰⁶Pb ratio measurement, the mean isotopic ratio obtained on National Institute of Standards and Technology (NIST) SRM610 over 4 months was 0.9105 ± 0.0014 (n = 280, 95% confidence), which agrees well with the published value of 0.9096. The time‐profile of Pb/U ratios during single spot ablation showed no significant difference in shape from NIST SRM610 and 91500 zircon standards. These results encouraged the use of the glass standard as a calibration standard for the Pb/U ratio determination for zircons with shorter wavelength (λ = 213 nm) laser ablation. But ²⁰⁶Pb/²³⁸U and ²⁰⁷Pb/²³⁵U ages obtained by this method for seven zircon standards are systematically younger than the published U‐Pb ages obtained by both isotope dilution–thermal ionization mass spectrometry (ID‐TIMS) and sensitive high‐resolution ion‐microprobe (SHRIMP). Greater discrepancies (3–4% younger ages) were found for the ²⁰⁶Pb/²³⁸U ages for SL13, AS3 and 91500 zircons. The origin of the differences could be heterogeneity in Pb/U ratio on SRM610 between the different disks, but a matrix effect accuracy either in the ICP ion source or in the ablation‐transport processes of the sample aerosols cannot be neglected. When the ²⁰⁶Pb/²³⁸U (= 0.2302) newly defined in the present study is used, the measured ²⁰⁶Pb/²³⁸U and ²⁰⁷Pb/²³⁵U ages for the seven zircon standards are in good agreement with those from ID‐TIMS and SHRIMP within ±2%. This suggests that SRM610 glass standard is suitable for ICP‐MS with laser ablation sampling (LA‐ICP‐MS) zircon analysis, but it is necessary to determine the correction factor for ²⁰⁶Pb/²³⁸U by measuring several zircon standards in individual laboratories.     

Simulations of 3D DC borehole resistivity measurements with a goal-oriented <Emphasis Type="Italic">hp</Emphasis> finite-element method. Part II: through-casing resistivity instruments

We simulate direct current (DC) borehole resistivity measurements acquired in steel-cased deviated wells for the assessment of rock formation properties. The assumed data acquisition configuration considers one current (emitter) and three voltage (collector) electrodes that are utilized to measure the second difference of the electric potential along the well trajectory. We assume a homogeneous, 1.27-cm-thick steel casing with resistivity equal to 10^− 5 Ω· m. Simulations are performed with two different numerical methodologies. The first one is based on transferring two-dimensional (2D) axisymmetric optimal grids to a three-dimensional (3D) simulation software. The second one automatically produces optimal 3D grids yielded by a 3D self-adaptive goal-oriented algorithm. Both methodologies utilize high-order finite elements (FE) that are specially well-suited for problems with high-contrast coefficients and rapid spatial variations of the electric field, as it occurs in simulations that involve steel-cased wells. The method based on transferring 2D-optimal grids is efficient in terms of CPU time (few seconds per logging position). Unfortunately, it may produce inaccurate 3D simulations in deviated wells, even though the error remains below 1% for the axisymmetric (vertical) well. The method based on optimal 3D grids, although less efficient in terms of CPU time (few hours per logging position), produces more accurate results that are validated by a built-in a posteriori error estimator. This paper provides the first existing simulations of through-casing resistivity measurements in deviated wells. Simulated resistivity measurements indicate that, for a 30° deviated well, measurements in conductive layers 0.01 Ω· m) are similar to those obtained in vertical wells. However, in resistive layers (10,000 Ω· m), we observe 100% larger readings in the 30° deviated well. This difference becomes 3,000% for the case of a 60° deviated well. For this highly-deviated well, readings corresponding to the conductive formation layer are about 30% smaller in magnitude than those in a vertical well. Shoulder effects significantly vary in deviated wells.     

A coupling of mixed and discontinuous Galerkin finite-element methods for poroelasticity

In this paper, we formulate a finite-element procedure for approximating the coupled fluid and mechanics in Biot’s consolidation model of poroelasticity. We approximate the flow variables by a mixed finite-element space and the displacement by a family of discontinuous Galerkin methods. Theoretical convergence error estimates are derived and, in particular, are shown to be independent of the constrained specific storage coefficient, c _o . This suggests that our proposed algorithm is a potentially effective way to combat locking, or the nonphysical pressure oscillations, which sometimes arise in numerical algorithms for poroelasticity.     

Fractional crystallization models and B–Be–Li systematics at Mt Somma-Vesuvius volcano (Southern Italy)

A compilation of B–Be–Li data on rocks that cover the entire eruptive history of Somma-Vesuvius is presented and interpreted in the light of evolution models for the Somma-Vesuvius rocks. Using major and trace element data, fractional crystalllization models are presented for different geochemical units. These data were used to constrain the source mineralogy of the Somma-Vesuvius rocks (ol-opx-cpx-gar-amp of 0.4-0.3-0.1-0.1-0.1), the amount of sediment added (5–10%) and the melt fraction from batch partial melting computations (0.05–0.1). From the B–Li data it is inferred that the main process responsible for the B isotopic signature is sediment recycling. However, the B–Li data show a major variation in Li abundances respect to B which is explained with Li dehydration before the fluid enriched the mantle wedge that produced the arc magmas. The Somma-Vesuvius B isotope composition is intermediate between that of the Campi Flegrei and the broad field of the Eolian Island arc. A low Be isotopes in the recent volcanic rocks can be explained as: (a) the top 1–22 m of the incoming sediment is accreted, (b) large amounts of sediment erosion, (c) a slow rate of subduction which have provoked a long magmatic history for the Vesuvius magma, (d) the sediment component takes several Myr longer than the subducting plate to reach the magma source region beneath Italy.     

黑龙江省东部桦南隆起美作花岗岩的锆石U-Pb定年及其地质意义

为查明马家街群发生接触变质作用的时代及其构造背景,对出露于黑龙江省东部桦南隆起的美作花岗岩体进行了锆石LA-ICP-MS U-Pb定年研究.结果显示:锆石具有清晰的生长振荡环带,其Th/U比值为0.12～1.04,属于典型的岩浆成因锆石;获得的206Pb/238U加权平均年龄为(259.0±3.6) Ma(n=14,MSWD=5.9),代表了岩体结晶年龄.美作岩体与佳木斯地块南部的青山、楚山和柴河岩体形成时代一致,其轻稀土富集、重稀土亏损型式和Eu负异常以及明显的Nb、Ta、Sr、Ti亏损的特点显示为壳源成因的火山弧花岗岩;它们的形成可能与晚古生代古亚洲洋的消减作用有关,该期花岗岩的就位导致马家街群发生了接触变质作用.美作花岗岩体形成时代的厘定,不仅限定了马家街群的变质时代为晚二叠世,同时为进一步探讨佳木斯地块晚古生代构造演化提供了新证据.