A semianalytical approach for solving first-order perturbation equations of dissolution-timescale reactive infiltration instability problems in fluid-saturated rocks

This paper presents a semianalytical approach for solving first-order perturbation (FOP) equations, which are used to describe dissolution-timescale reactive infiltration instability (RII) problems in fluid-saturated rocks. The proposed approach contains two parts because the chemical dissolution reaction divides the whole problem domain into two subdomains. In the first part, the interface-condition substitution strategy is used to derive the analytical expressions of purely mathematical solutions for the FOP equations in the upstream subdomain, where the dissolution chemical reaction is ceased and the FOP equations are weakly coupled. In the second part, the finite element method (FEM) is used to derive the analytical expressions of numerical solutions for the FOP equations in the downstream subdomain, where the dissolution chemical reaction needs to be considered and the FOP equations are strongly coupled so that it is impossible to derive purely mathematical solutions for them. Particular attention is paid to the development of the element-by-element forward marching strategy, which is associated with the use of the FEM for solving this new kind of scientific problem. The related analytical results demonstrated that (1) both the dynamic characteristic of a reactive infiltration system and the dimensionless wavenumber can have pronounced influences on the distribution of the FOP dimensionless acid concentration within the entire domain of the dissolution-timescale RII problems in fluid-saturated rocks and (2) the FOP dimensionless acid concentration distribution exhibits two significantly different patterns in the upstream and downstream subdomains of the dissolution-timescale RII system.     

Development of a GPGPU-parallelized hybrid finite-discrete element method for modeling rock fracture

The hybrid finite-discrete element method (FDEM) is widely used for engineering applications, which, however, is computationally expensive and needs further development, especially when rock fracture process is modeled. This study aims to further develop a sequential hybrid FDEM code formerly proposed by the authors and parallelize it using compute unified device architecture (CUDA) C/C++ on the basis of a general-purpose graphics processing unit (GPGPU) for rock engineering applications. Because the contact detection algorithm in the sequential code is not suitable for GPGPU parallelization, a different contact detection algorithm is implemented in the GPGPU-parallelized hybrid FDEM. Moreover, a number of new features are implemented in the hybrid FDEM code, including the local damping technique for efficient geostatic stress analysis, contact damping, contact friction, and the absorbing boundary. Then, a number of simulations with both quasi-static and dynamic loading conditions are conducted using the GPGPU-parallelized hybrid FDEM, and the obtained results are compared both quantitatively and qualitatively with those from either theoretical analysis or the literature to calibrate the implementations. Finally, the speed-up performance of the hybrid FDEM is discussed in terms of its performance on various GPGPU accelerators and a comparison with the sequential code, which reveals that the GPGPU-parallelized hybrid FDEM can run more than 128 times faster than the sequential code if it is run on appropriate GPGPU accelerators, such as the Quadro GP100. It is concluded that the GPGPU-parallelized hybrid FDEM developed in this study is a valuable and powerful numerical tool for rock engineering applications.     

陕西蓝田灞源一带金矿地质特征及成因浅析

本文通过分析灞源地区的赋矿地层、岩浆岩及断裂构造与金矿形成的关系，查明了控制金矿形成与分布的地质条件，并重点论述了金矿床的产出、矿化类型、围岩蚀变和含金矿物等主要地质特征。最后根据矿床地球化学特征综合分析阐明本区金矿床的岩浆热液成因。     

华北陆块南缘富碱侵入岩岩石组合及时空分布

华北陆块南缘（河南段）出露的富碱侵入岩类集中分布于卢氏－确山一线。岩石类型主要分为三大类：①碱性岩类，即含有似长石或碱性暗色矿物的正长岩类；②碱性花岗岩类，包括钠铁闪石花岗岩以及孪生的钾长花岗岩类；③石英正长岩类，包括碱性长石为主的石英正长岩、英碱正长岩和花岗正长（斑）岩类。在空间上构成一条富碱岩浆岩带，据岩石组合自北而南划为三个亚带：北部碱性岩亚带、中部碱性花岗岩类亚带和南部石英正长（斑）岩亚带。在形成时间上富碱岩浆活动集中于晚元古代末期（７００Ｍａ±）和海西晚期（３００Ｍａ±）。     

大坂山公路隧道工程地质特征

大坂山公路隧道位于青海省西宁市—甘肃省张掖市公路(即227国道),自南向北穿越大坂山主峰。隧道轴线长约1.538km,轴线走向29°,系国内海拔最高的公路隧道。由于隧道址区所处的特殊地形地质背景条件,必须查明一些特殊的隧道工程地质问题。经工程地质勘探查明了:一、大坂山址区不存在多年冻土,仅发育季节冻土,对隧道工程危害不大;二、隧道围岩类型分为四类、11段,其中稳定性好的Ⅴ类围岩占隧道全长51％,Ⅳ类围岩占15％,Ⅲ类围岩占2％,稳定性差的Ⅱ类围岩占32％;三、隧道址区地壳稳定程度属次不稳定级,为6.5级地震潜在震源区,在五十年地震活动强度超越概率0.1的条件下,相应地震烈度值为7.8度;四、隧道址区位于分水岭地带,水文地质条件简单,大气降水和地表水不易补给地下水,预测隧道总涌水量仅668m~3/d。     

陕西潼关地区太华群变质岩原岩性质的恢复

潼关地区太华群地层是一套以角闪岩相为主(局部达麻粒岩相)的中深变质岩系。本区太华群变质岩的主要岩石类型:斜长片麻岩类的原岩为中酸性火山岩(二长安山质、英安质凝灰岩);浅(变)粒岩类的原岩为中酸性火山岩(英安岩,流纹岩)一中酸性火山碎屑岩(二长安山质、英安质凝灰岩),并以后者为主,还有少量粘土质沉积岩或长石质杂砂岩;斜长角闪岩的原岩为基性火山岩—拉斑玄武岩。     

Methodology for groundwater recharge assessment in carbonate aquifers: application to pilot sites in southern Spain

Aquifer recharge can be determined by conventional methods such as hydrodynamic or hydrologic balance calculations, or numerical, hydrochemical or isotopic models. Such methods are usually developed with respect to detrital aquifers and are then used on carbonate aquifers without taking into consideration their hydrogeological particularities. Moreover, such methods are not always easy to apply, sometimes requiring data that are not available. Neither do they enable determination of the spatial distribution of the recharge. For eight regions in southern Spain, the APLIS method has been used to estimate the mean annual recharge in carbonate aquifers, expressed as a percentage of precipitation, based on the variables altitude, slope, lithology, infiltration landform, and soil type. The aquifers are representative of a broad range of climatic and geologic conditions. Maps of the above variables have been drawn for each aquifer, using a geographic information system; thus they can be superimposed to obtain the mean value and spatial distribution of the recharge. The recharge values for the eight aquifers are similar to those previously calculated by conventional methods and confirmed by discharge values, which corroborates the validity of the method.     

Fate of individual sewage disposal system wastewater within regolith in mountainous terrain

In order to improve understanding of the fate of septic tank or individual sewage disposal system (ISDS) effluent in regolith overlying fractured-rock aquifers, effluent from an ISDS in such a setting was tracked via geophysical, hydrological, and geochemical methods. Under typical precipitation conditions, the effluent entered the fractured bedrock within 5 m of the boundary of the constructed infiltration area. During a period of unusually high spring recharge, the plume migrated between 50 and 100 m within the regolith before infiltrating the fractured bedrock. The chemical signature of the effluent is similar to that required to account for the decline in water quality, suggesting a causative relationship (as estimated from mass-balance models of the surface-water chemistry near the mouth of the basin). The elevated salt content of the effluent during periods of high natural recharge to the infiltration area correlates with elevated salt concentrations in surface and groundwater at the basin scale, suggesting that some of the effluent salt load may be stored in the unsaturated zone during dry periods and flushed during periods of elevated natural recharge.     

Fuzzy logic modeling of the resistivity parameter and topography features for aquifer assessment in hydrogeological investigation of a crystalline basement complex

A linguistic fuzzy logic system (LFLS)-based expert system model has been developed for the assessment of aquifers for the location of productive water boreholes in a crystalline basement complex. The model design employed a multiple input/single output (MISO) approach with geoelectrical parameters and topographic features as input variables and control crisp value as the output. The application of the method to the data acquired in Khondalitic terrain, a basement complex in Vizianagaram District, south India, shows that potential groundwater resource zones that have control output values in the range 0.3295–0.3484 have a yield greater than 6,000 liters per hour (LPH). The range 0.3174–0.3226 gives a yield less than 4,000 LPH. The validation of the control crisp value using data acquired from Oban Massif, a basement complex in southeastern Nigeria, indicates a yield less than 3,000 LPH for control output values in the range 0.2938–0.3065. This validation corroborates the ability of control output values to predict a yield, thereby vindicating the applicability of linguistic fuzzy logic system in siting productive water boreholes in a basement complex.     

北祁连黑矿型和塞浦路斯型硫化物矿床容矿火山岩的物质来源与形成环境

文章以白银厂石英角斑岩和石居里富钠玄武岩为代表,对产于北祁连造山带的"黑矿型"和"塞浦路斯型"铜多金属块状硫化物矿床的容矿火山岩做了系统的常量元素、微量元素、稀土元素及锶、钕、铅同位素分析。结果表明,这2种火山岩均表现出俯冲带岩浆作用的特征。所不同的是,白银厂"黑矿型"块状硫化物矿床容矿火山岩形成于活动大陆边缘的陆缘弧环境;石居里"塞浦路斯型"硫化物铜矿床形成于大洋板块内部的洋生弧环境。与各自的形成环境相对应,白银厂的石英角斑岩来源于原生地幔在俯冲板片流体参与下的部分熔融,岩浆在上升侵位过程中有部分大陆地壳物质混入。石居里的富钠玄武岩同样也是产生于原生地幔在板片流体作用下的部分熔融,但没有遭受陆壳物质混染。