成像光谱矿物填图精度兼与矿物丰度关系初探

成像光谱遥感技术在岩石矿物识别和矿物填图方面发挥着越来越大作用,而成像光谱矿物填图的精度问题一直是研究的焦点。文章利用高光谱数据HyMap对新疆东天山地区的矿物填图结果,通过引入矿物填绘分布强度的概念,探索性的分析了成像光谱矿物填图与填绘矿物丰度的关系。     

经验模式分解的尺度阈值滤波在变形分析中的应用

本文将信号处理领域的经验模式分解算法应用于变形信息的提取中,通过引入阈值函数,建立了基于经验模式分解的尺度阈值滤波模型,采用优化模型确定了经验模式分解的次数。分别通过模拟试验和实测数据与小波阈值法和多项式拟合法进行了比对,分析表明:在低噪声情况下,三种方法都有一定的滤波效果;在高噪声情况下,经验模式分解的尺度阈值滤波法具有与小波阈值去噪法相等的精度,而且瞬时强噪声识别能力更好,优于多项式拟合法。     

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.     

Zircon M257 ‐ a Homogeneous Natural Reference Material for the Ion Microprobe U‐Pb Analysis of Zircon

We introduce and propose zircon M257 as a future reference material for the determination of zircon U‐Pb ages by means of secondary ion mass spectrometry. This light brownish, flawless, cut gemstone specimen from Sri Lanka weighed 5.14 g (25.7 carats). Zircon M257 has TIMS‐determined, mean isotopic ratios (2s uncertainties) of 0.09100 ± 0.00003 for ²⁰⁶pb/²³⁸U and 0.7392 ± 0.0003 for ²⁰⁷pb/²³⁵U. Its ²⁰⁶pb/²³⁸U age is 561.3 ± 0.3 Ma (unweighted mean, uncertainty quoted at the 95% confidence level); the U‐Pb system is concordant within uncertainty of decay constants. Zircon M257 contains ～ 840 μg g^?1 U (Th/U ～ 0.27). The material exhibits remarkably low heterogeneity, with a virtual absence of any internal textures even in cathodoluminescence images. The uniform, moderate degree of radiation damage (estimated from the expansion of unit‐cell parameters, broadening of Raman spectral parameters and density) corresponds well, within the “Sri Lankan trends”, with actinide concentrations, U‐Pb age, and the calculated alpha fluence of 1.66 × 10¹⁸ g^?1. This, and a (U+Th)/He age of 419 ± 9 Ma (2s), enables us to exclude any unusual thermal history or heat treatment, which could potentially have affected the retention of radiogenic Pb. The oxygen isotope ratio of this zircon is 13.9%o VSMOW suggesting a metamorphic genesis in a marble or calc‐silicate skarn.     

Recommended Values of 239Pu, 240Pu and 239+240Pu Concentrations in Reference Material IAEA-315 (Marine Sediment) Estimated by Thermal Ionisation Mass Spectrometry, Inductively Coupled Plasma-Mass Spectrometry and Alpha Spectrometry

The recommended concentrations of ²³⁹Pu, ²⁴⁰Pu and ^239+240Pu in reference material IAEA‐315 (marine sediment) were estimated by three analytical methods: isotope dilution thermal ionisation mass spectrometry (TIMS), isotope dilution inductively coupled plasma‐mass spectrometry (ICP‐MS) and alpha spectrometry. The determination of ²³⁹Pu and ²⁴⁰Pu (^239+240Pu by alpha spectrometry) was carried out with samples from randomly selected bottles using each method. Plutonium‐238 was also measured by alpha spectrometry. A plutonium‐242 reference material was used as a spike for the quantitative analysis. The influence of ²⁴²Pu in the samples was therefore calculated; however, this contribution was less than the range of uncertainty and did not influence the final results. The obtained data were statistically analysed using variance component analysis and paired comparison. The combined standard uncertainties from “method/measurement”, “bottle” and “sub‐sample” were in the order of 3 to 6%. The main contributions to the uncertainty were from the material heterogeneity and from systematic differences between methods. Based on this study with twenty‐seven analyses using 10–14 g sample mass, concentrations of (38 ± 3) Bq kg^?1, (28 ± 3) Bq kg^?1 and (66 ± 4) Bq kg^?1 are proposed as recommended values for ²³⁹Pu, ²⁴⁰Pu and ^239+240Pu, respectively, and (9.5 ± 0.4) Bq kg^?1 for ²³⁸Pu as an information value in reference material IAEA‐315. In mass concentration units, these amount to (16.4 ± 1.2) ng kg^?1, (3.3 ± 0.4) ng kg^?1 and (0.015 ± 0.003) ng kg^?1 for ²³⁹Pu, ²⁴⁰Pu and ²³⁸Pu, respectively. The certified reference materials NIST 4350B and NIST 4354 were also analysed by TIMS for quality assurance of the method used in this study.     

A dynamic box model of bioactive elements in the southern Taiwan Strait

A dynamic box model was applied to study the characteristics of biogeochemical cycling of PO_4-P,NO_3-N,AOU,POC and PON in the southern Taiwan Strait region based on field data of the"Minnan Taiwan Bank Fishing Ground Upwelling Ecosystem Study" during the period of Dec.1987-Nov.1988.According to the unique hydrological and topographical features of the region,six boxesand three layers were considered in the model.The variation rates and fluxes of elements induced by hor-izontal current,upwelling,by diffusion,sinking of particles and biogeochemical processes were estimatedrespectively.Results further confirmed that upwellings had important effects in this region.Thenearshore upwelling areas had net input fluxes of nutrients brought by upwelling water,also had high de-pletion rates of nutrients and production rates of particulate organic matter and dissolved oxygen.Theabnormal net production of nutrients in the middle layer(10-30m) indicated the important role of bacte-ria in this high production region.Th