Simulations of a coupled hydro-chemo-mechanical system in rocks

A coupled hydro-chemo-mechanical numerical model is developed for these coupled phenomena in many engineering fields. The model has been applied to predicting the response of a stressed rockmass column to an injected reactive fluid (reagent) flow. The response includes evolutions of porosity, permeability, reagent and mineral concentrations during dissolution. In the model, the progress of dissolution is defined by the change in porosity ratio and the porosity increases with dissolution assuming there is no precipitation. The numerical evolutions of porosity, permeability, reagent and mineral concentrations during dissolution are validated against steady state solutions. The model results show that these evolutions are regulated to a certain extent by the applied external loadings: an applied extensional stress enhances the progress of the dissolution process while an applied compression stress slows the progress of the dissolution process. This revised version was published online in July 2006 with corrections to the Cover Date.     

Preliminary Characterisation of New Glass Reference Materials (GSA-1G, GSC-1G, GSD-1G and GSE-1G) by Laser Ablation-Inductively Coupled Plasma-Mass Spectrometry Using 193 nm, 213 nm and 266 nm Wavelengths

New glass reference materials GSA-1G, GSC-1G, GSD-1G and GSE-1G have been characterised using a prototype solid state laser ablation system capable of producing wavelengths of 193 nm, 213 nm and 266 nm. This system allowed comparison of the effects of different laser wavelengths under nearly identical ablation and ICP operating conditions. The wavelengths 213 nm and 266 nm were also used at higher energy densities to evaluate the influence of energy density on quantitative analysis. In addition, the glass reference materials were analysed using commercially available 266 nm Nd:YAG and 193 nm ArF excimer lasers. Laser ablation analysis was carried out using both single spot and scanning mode ablation. Using laser ablation ICP-MS, concentrations of fifty-eight elements were determined with external calibration to the NIST SRM 610 glass reference material. Instead of applying the more common internal standardisation procedure, the total concentration of all element oxide concentrations was normalised to 100%. Major element concentrations were compared with those determined by electron microprobe. In addition to NIST SRM 610 for external calibration, USGS BCR-2G was used as a more closely matrix-matched reference material in order to compare the effect of matrix-matched and non matrix-matched calibration on quantitative analysis. The results show that the various laser wavelengths and energy densities applied produced similar results, with the exception of scanning mode ablation at 266 nm without matrix-matched calibration where deviations up to 60% from the average were found. However, results acquired using a scanning mode with a matrix-matched calibration agreed with results obtained by spot analysis. The increased abundance of large particles produced when using a scanning ablation mode with NIST SRM 610, is responsible for elemental fractionation effects caused by incomplete vaporisation of large particles in the ICP.     

Three dimensional numerical simulation of residential building on shrink–swell soils in response to climatic conditions

Shrink–swell soils can cause distresses in buildings, and every year, the economic loss associated with this problem is huge. This paper presents a comprehensive system for simulating the soil–foundation–building system and its response to daily weather conditions. Weather data include rainfall, solar radiation, air temperature, relative humidity, and wind speed, all of which are readily available from a local weather station or the Internet. These data are used to determine simulation flux boundary conditions. Different methods are proposed to simulate different boundary conditions: bare soil, trees, and vegetation. A coupled hydro‐mechanical stress analysis is used to simulate the volume change of shrink–swell soils due to both mechanical stress and water content variations. Coupled hydro‐mechanical stress‐jointed elements are used to simulate the interaction between the soil and the slab, and general shell elements are used to simulate structural behavior. All the models are combined into one finite element program to predict the entire system's behavior. This paper first described the theory for the simulations. A site in Arlington, Texas, is then selected to demonstrate the application of the proposed system. Simulation results are shown, and a comparison between measured and predicted movements for four footings in Arlington, Texas, over a 2‐year period is presented. Finally, a three‐dimensional simulation is made for a virtual residential building on shrink–swell soils to identify the influence of various factors. Copyright © 2015 John Wiley & Sons, Ltd.     

Elastic lateral dynamic impedance functions for a rigid cylindrical shell type foundation

Elastic lateral dynamic impedance functions are defined as the ratio of the lateral dynamic force/moment to the corresponding lateral displacement/rotation at the top ending of a foundation at very small strains. Elastic lateral dynamic impedance functions have a defining influence on the natural frequencies of offshore wind turbines supported on cylindrical shell type foundations, such as suction caissons, bucket foundations, and monopiles. This paper considers the coupled horizontal and rocking vibration of a cylindrical shell type foundation embedded in a fully saturated poroelastic seabed in contact with a seawater half‐space. The formulation of the coupled seawater–shell–seabed vibration problem is simplified by treating the shell as a rigid one. The rigid shell vibration problem is approached by the integral equation method using ring‐load Green's functions for a layered seawater‐seabed half‐space. By considering the boundary conditions at the shell–soil interface, the shell vibration problem is reduced to Fredholm integral equations. Through an analysis of the corresponding Cauchy singular equations, the intrinsic singular characteristics of the problem are rendered explicit. With the singularities incorporated into the solution representation, an effective numerical method involving Gauss–Chebyshev method is developed for the governing Fredholm equations. Selected numerical results for the dynamic contact load distributions, displacements of the shell, and lateral dynamic impedance functions are examined for different shell length–radius ratio, poroelastic materials, and frequencies of excitation. Copyright © 2016 John Wiley & Sons, Ltd.     

Modelling landslides in unsaturated slopes subjected to rainfall infiltration using material point method

This paper presents a dynamic fully coupled formulation for saturated and unsaturated soils that undergo large deformations based on material point method. Governing equations are applied to porous material while considering it as a continuum in which the pores of the solid skeleton are filled with water and air. The accuracy of the developed method is tested with available experimental and numerical results. The developed method has been applied to investigate the failure and post‐failure behaviour of rapid landslides in unsaturated slopes subjected to rainfall infiltration using two different bedrock geometries that lie below the top soil. The models show different failure and post‐failure mechanisms depending on the bedrock geometry and highlight the negative effects of continuous rain infiltrations. Copyright © 2016 John Wiley & Sons, Ltd.     

利用LA-ICP-MS在16 μm和10 μm激光束斑条件下测定独居石U-Th-Pb年龄

我国锂钾主要矿产资源大多分布在西部偏远地区,其勘查找矿或综合利用迫切需要现场快速分析技术的支持。本文介绍了自行开发的便携式Li-K分析仪的主要性能及其在锂辉石中锂的分析应用。Li-K分析仪是一种基于大气压的液体阴极辉光放电光谱仪,以待测液体为放电阴极,实现样品中Li、K等元素的原子化和激发。该仪器以光纤传导CCD光谱仪作为检测器,其波长范围为345～1015 nm,分辨率3 nm;以Li 670.78 nm和K 769.90 nm分析谱线,在最佳仪器工作条件下测定Li、K的精密度（RSD,n=14）均低于2%,检出限为0.03 μg/mL,检测范围0.1～10 μg/mL。研究表明,不同酸度和酸的类型对谱线强度影响较大,且存在显著的样品基体效应,标准曲线法分析锂辉石中Li的结果偏差高达267%;而标准加入法可克服基体效应的影响,获得与ICP-OES一致性较好的分析结果,为现场开展固体样品中Li的测定奠定了基础。     

Sources of anthropogenic fire ignitions on the peat-swamp landscape in Kalimantan,Indonesia

Fire disturbance in many tropical forests, including peat swamps, has become more frequent and extensive in recent decades. These fires compromise a variety of ecosystem services, among which mitigating global climate change through carbon storage is particularly important for peat swamps. Indonesia holds the largest amount of tropical peat carbon globally, and mean annual CO₂ emissions from decomposition of deforested and drained peatlands and associated fires in Southeast Asia have been estimated at ∼2000 Mt y-1. A key component to understanding and therefore managing fire in the region is identifying the land use/land cover classes associated with fire ignitions. We assess the oft-asserted claim that escaped fires from oil palm concessions and smallholder farms near settlements are the primary sources of fire in a peat-swamp forest area in Central Kalimantan, Indonesia, equivalent to around a third of Kalimantan's total peat area. We use the MODIS Active Fire product from 2000 to 2010 to evaluate the fire origin and spread on the land use/land cover classes of legal, industrial oil palm concessions (the only type of legal concession in the study area), non-forest, and forest, as well as in relation to settlement proximity. We find that most fires (68–71%) originate in non-forest, compared to oil palm concessions (17%–19%), and relatively few (6–9%) are within 5 km of settlements. Moreover, most fires started within oil palm concessions and in close proximity to settlements stay within those boundaries (90% and 88%, respectively), and fires that do escape constitute only a small proportion of all fires on the landscape (2% and 1%, respectively). Similarly, a small proportion of fire detections in forest originate from oil palm concessions (2%) and within close proximity to settlements (2%). However, fire ignition density in oil palm (0.055 ignitions km⁻²) is comparable to that in non-forest (0.060 km-2 ignitions km-²), which is approximately ten times that in forest (0.006 ignitions km⁻²). Ignition density within 5 km of settlements is the highest at 0.125 ignitions km⁻². Furthermore, increased anthropogenic activity in close proximity to oil palm concessions and settlements produces a detectable pattern of fire activity. The number of ignitions decreases exponentially with distance from concessions; the number of ignitions initially increases with distance from settlements, and, around from 7.2 km, then decreases with distance from settlements. These results refute the claim that most fires originate in oil palm concessions, and that fires escaping from oil palm concessions and settlements constitute a major proportion of fires in this study region. However, there is a potential for these land use types to contribute substantially to the fire landscape if their area expands. Effective fire management in this area should therefore target not just oil palm concessions, but also non-forested, degraded areas where ignitions and fires escaping into forest are most likely to occur.     

电感耦合等离子体质谱法(ICP-MS)测定水系沉积物中的铌钽锆铪

建立电感耦合等离子体质谱法(ICP-MS)测定水系沉积物中铌、钽、锆、铪四种元素的分析方法。将样品与氢氧化钠、过氧化钠混合物放入高温炉中,熔融分解完全,用热水提取,过滤后弃去滤液,将滤纸及沉淀用酒石酸-盐酸溶液溶解,稀释至刻度测定。方法检出限(3s)为:LD(Nb)=0.08μg/g,LD(Ta)=0.04μg/g,LD(Zr)=0.5μg/g,LD(Hf)=0.04μg/g,精密度(RSD%,n=6)为:0.84%~4.21%。该测定方法具有灵敏度高、精密度好、分析速度快、线性范围宽、操作性强等优势。采用该方法对国家一级标准物质进行测定表明,其结果与标准值吻合。此方法已在实际地质调查样品分析中得到应用。     

均生函数与BP神经网络耦合模拟预测模型(MGF-BP-I)的建立与应用

水文气候因子模拟预测对气候变化研究、农业墒情预报、生态环境改善、水资源合理开发利用等具有一定参考意义。均生函数、BP神经网络及其结合改进方式在模拟预测中各有优点,被广泛应用,但仍有进一步改进空间。针对MGF、MGF-OSR、MGF-OSR-BP等方法粗选因子集、粗选集组合筛选、收敛适应性、精度控制等可改进空间,进一步发挥均生函数和BP神经网络优势,建立了MGF-BP-I模拟预测模型。利用MGF-OSR、MGF-OSR-BP、MGF-BP-I对科尔沁沙地区域平均年降水进行了模拟预测。结果表明,建模期MGF-OSR-BP、MGF-BP-I拟合效果均较好,MGF-BP-I建模阶段最优模式精度优于MGF-OSR-BP,MGF-BP-I整体同时最优模式结果也非常好。检验期,MGF-BP-I检验阶段最优及整体同时最优两种模式拟合效果最好,相比其他模式精度有所提高。MGF-BP-I考虑更加全面,充分发挥了均生函数和BP神经网络优势,精度远高于MGF-OSR和MGF-OSR-BP,MGF-BP-I整体同时最优模式更符合实际应用,效果理想,可用于水文气候因子模拟预测。     

An Investigation of Digestion Methods for Trace Elements in Bauxite and Their Determination in Ten Bauxite Reference Materials Using Inductively Coupled Plasma‐Mass Spectrometry

Trace elements from samples of bauxite deposits can provide useful information relevant to the exploration of the ore‐forming process. Sample digestion is a fundamental and critical stage in the process of geochemical analysis, which enables the acquisition of accurate trace element data by ICP‐MS. However, the conventional bomb digestion method with HF/HNO₃ results in a significant loss of rare earth elements (REEs) due to the formation of insoluble AlF₃ precipitates during the digestion of bauxite samples. In this study, the digestion capability of the following methods was investigated: (a) ‘Mg‐addition’ bomb digestion, (b) NH₄HF₂ open vessel digestion and (c) NH₄F open vessel digestion. ‘Mg‐addition’ bomb digestion can effectively suppress the formation of AlF₃ and simultaneously ensure the complete decomposition of resistant minerals in bauxite samples. The addition of MgO to the bauxite samples resulted in (Mg + Ca)/Al ratios ≥ 1. However, adding a large amount of MgO leads to significant blank contamination for some transition elements (V, Cr, Ni and Zn). The NH₄HF₂ or NH₄F open vessel digestion methods can also completely digest resistant minerals in bauxite samples in a short period of time (5 hr). Unlike conventional bomb digestion with HF/HNO₃, the white precipitates and the semi‐transparent gels present in the NH₄HF₂ and NH₄F digestion methods could be efficiently dissolved by evaporation with HClO₄. Based on these three optimised digestion methods, thirty‐seven trace elements including REEs in ten bauxite reference materials (RMs) were determined by ICP‐MS. The data obtained showed excellent inter‐method reproducibility (agreement within 5% for REEs). The relative standard deviation (% RSD) for most elements was < 6%. The concentrations of trace elements in the ten bauxite RMs showed agreement with the limited certified (Li, V, Cr, Cu, Zn, Ga, Sr, Zr and Pb) and information values (Co, Ba, Ce and Hf) available. New trace element data for the ten RMs are provided, some of which for the first time.