基于数值模拟的再生水对地下水质影响实例

为了探讨在再生水地表回灌条件下,再生水对地下水中氨素和主要盐分浓度变化的影响,选取现实存在的清河再生水湿地公园为例,以氨氮、Cl-离子为主要研究对象,利用Hydrus-1D溶质运移模型和Modflow、MT3DMSD耦合地下水模型,模拟了再生水长期利用后的地下水氨氮和Cl-指标变化情况。研究结果显示,再生水到达包气带底部时,水中氮氨、Cl-离子含量高于研究区潜水地下水中的相应指标含量,氨氮超出地下水质量Ⅲ类标准。预测5~20年后,潜水中氨氮超标面积不断扩展,20年后将达14km2,氨氮超标垂向延伸至第一层弱透水层;Cl-离子在潜水中不超标,但在湿地公园周边浓度有所增加。湿地公园再生水的长期地表回灌将给当地地下水质量带来负面影响。     

塔里木陆块西北缘萨热克砂岩型铜矿床构造-流体演化对成矿的制约

塔里木陆块西北缘萨热克砂岩型铜矿床构造演化、流体演化与成矿之间具有密切关系,处于一个统一系统中。矿床成岩期方解石中包裹体水的δD值为-65.3‰~-99.2‰,改造成矿期石英包裹体水的δD值为-77.7‰~-96.3‰,成岩成矿期成矿流体δ~(18)OH_2O变化范围为-3.22‰~1.84‰,改造成矿期成矿流体δ~(18)OH_2O变化范围为-4.26‰~5.14‰,指示萨热克铜矿成岩期、改造期成矿流体主要为中生代大气降水及其经水岩作用而成的盆地卤水。矿石中辉铜矿δ~(34)S值为-24.7‰~-15.4‰,指示硫主要源自硫酸盐细菌与有机质还原,部分源于有机硫。构造与成矿流体演化对砂岩铜矿成矿起关键制约作用。盆地发展早期强烈的抬升运动使盆地周缘基底与古生界剥蚀,为富铜矿源层的形成提供了丰富物源,至晚侏罗世盆地发展晚期,长期演化积聚的巨量含矿流体在库孜贡苏组砾岩胶结物及裂隙中富集,在萨热克巴依盆地内形成具有经济意义的砂岩型铜矿床。     

Development and verification of the comprehensive model for physical properties of hydrate sediment

Natural gas hydrate is widely distributed all over the world and may be a potential resource in the near future, whereas hydrate dissociation during the development affects wellbore stability and drilling safety. However, the present modeling of hydrate reservoir parameters ignored the influence of effective stress and only considered the hydrate saturation. In this paper, a series of stress sensitivity experiments for the unconsolidated sandstone were carried out, and the influence of mean effective stress on physical parameters was obtained; a comprehensive model for the physical parameters of hydrate reservoir was developed subsequently. With the help of ABAQUS finite element software, the established comprehensive model was verified by the use of the wellbore stability numerical model of hydrate reservoir. The verification results show that ignoring the effect of mean effective stress on the parameters of hydrate formation aggravates the invasion of drilling fluid into the hydrate formation. Besides, ignoring the stress sensitivity of reservoir physical parameters will underestimate the wellbore instability during hydrate drilling, which will be a threat to the safety of gas hydrate drilling. At the end of the drilling operation, the maximum plastic strain of the model for considering and not considering stress sensitivity was 0.0145 and 0.0138, respectively. Therefore, the established comprehensive model will provide a theoretical support for accurately predicting the engineering geological disasters in hydrate development process.     

Evolution of the shape parameter in the Weibull distribution for brittle rocks under uniaxial compression

Applying the statistical damage theory based on the Weibull distribution to describe rock deformation and failure processes is an important development in rock mechanics. The shape parameter of the Weibull distribution, m, determines the basic shape of the distribution curve; additionally, it also represents a physical characteristic which can be applied when constructing rock constitutive models. To study the evolution of m during rock failure when applying the Weibull distribution to rock mechanics, uniaxial cyclic loading tests of shale specimens were conducted and previous rock mechanics experiments under different temperatures and loading rates were reviewed. The results indicate that m varied with the accumulation of damage but was almost constant between the volume expansion point and the peak strength point of each specimen. Combined with previous studies about the accelerated failure behavior of rocks, we conclude that between the volume expansion point and the peak strength point, the mechanical behavior of the rock fracture process did not change significantly. Based on the characteristics of m at different damage stages during the rock failure process, ranges of m values at different damage stages are proposed. The conclusions reached in this study may be used as an important reference for theoretical research on rock mechanics.     

Porosity and wave velocity evolution of granite after high-temperature treatment: a review

The evolution of porosity and changes in wave velocity in granite after high-temperature treatment has been experimentally investigated in different studies. Statistical analysis of the test results shows that there is a temperature threshold value that leads to variations in porosity and wave velocity. At a temperature that is less than 200 °C, the porosity of granite slowly increases with increases in temperature, while the wave velocity decreases. When the temperature is greater than 200 °C (especially between 400 and 600 °C), the porosity quickly increases, while the wave velocity substantially decreases. The temperature ranges of room temperature to 200 and 200–400 °C correspond to the undamaged state and the micro-damage state, respectively. The results confirm that there is an important link between the variations of physical and mechanical properties in response to thermal treatment. By studying the relationships among rock porosity, wave velocity and temperature, this provides the basis for solving multi-variable coupling problems under high temperatures for the thermal exploitation of petroleum and safe disposal of nuclear waste.     

A GIS-based method of risk assessment on no. 11 coal-floor water inrush from Ordovician limestone in Hancheng mining area,China

The no. 11 coal seam in the deep area of Hancheng mining area is mining in recent years, which is threatened by the water inrush from the Ordovician limestone aquifer. Coal-floor water inrush is governed by the water abundance of coal-floor aquifer, the water-resisting performance of coal-floor aquitard, and the pathway connecting the water source and the working face. To make an accuracy risk assessment of water inrush from the no. 11 coal seam floor, a GIS-based vulnerability index method (VIM) is adopted for its superior comprehensive consideration of more controlling factors, powerful spatial analysis, and intuitively display functions. This study firstly established an index system including the water pressure of the coal-floor aquifer, the unit water inflow, the thickness, the core recovery percentage, the thickness ratio of brittle rocks to ductile rocks, the thickness of effective aquitard, and the accumulated length of faults and folds, of which the former six indexes governed the water abundance of the coal-floor aquifer which was combined with the last two factors to determine the risk of coal-floor water inrush. Secondly, the thematic map of each controlling factor is established by GIS using the geological prospecting data, and the weight of each factor is determined by the analytic hierarchy process (AHP) after consulting the expert review panel. At last, a vulnerability index is obtained and used to assess the risk of coal-floor water inrush of the no. 11 coal seam. The risk of water inrush of the no. 11 coal seam of the study area was ranked to three zones: the southeastern shallow area in red color is the dangerous zone, the wide northwestern area in green color is the safe zone, and the transition area in yellow color is the moderate-risk zone. Compared with the actual water-inrush incidents, the risk assessment result was verified to achieve an accuracy of 82.35%, which is proved to be a dependable reference for the prevention and controlling of coal-floor water inrush of the no. 11 coal seam in Hancheng mining area.     

Stability study on the northern batter of MBC Open Pit using Plaxis 3D

Cracks appeared on the northern batter at Maddingley Brown Coal Open Pit Mine, Victoria, Australia, on 8 November 2013 and a 2-day rainfall event happened 5 days later. This study models the stability of the northern batter considering the effect of the rainfall event and an emergency buttress using finite element method (FEM) encoded in Plaxis 3D. It is found that the batter tended to lead to block sliding after overburden removal. The observed vertical crack would be a combined action of the overburden removal and groundwater flow. The simulated location of cracks agrees well with the actual location, and the simulated heave of the coal seam is in good agreement with the experience in Victoria brown coal open pit mining. The rainfall accelerated the development of the cracks. With the construction of the emergency buttress, the batter became stable that is in good agreement with the monitored data.     

湘西-黔东下寒武统清虚洞组藻灰岩沉积特征及其找矿意义

铅锌矿是湘西-黔东地区的优势矿种,资源丰富,开发历史悠久,找矿潜力巨大。综合分析前人有关湘西-黔东生物礁的资料,介绍了藻灰岩的沉积特征,并结合大量实际地质资料,对藻灰岩控矿进行了论述。结果表明: 成矿物质主要来源于藻类及碳酸盐泥对Pb²⁺、Zn²⁺离子的吸取,矿床成因属于沉积成岩矿床类型,兼有成岩期后矿床性质; 藻礁灰岩与不同岩性的接合部位及其附近等微地球化学障区,往往就是铅锌富矿体产出部位; 藻礁灰岩与礁间通道的薄层泥质白云质灰岩呈指状交叉接触处,也常有铅锌富矿体产出; 铅锌富矿体基本上产于清虚洞组灰岩段中,且明显受该段藻礁灰岩控制。为满足国内外对铅锌矿日益增长的消费需要,深入研究区内铅锌矿的地层岩相岩石控矿因素及其富集规律,指导该地区藻灰岩中铅锌矿找矿工作,具有重要的现实意义。     

全国地热资源调查评价与勘查示范工程进展

地热作为一种清洁能源,在我国经济发展与能源转型过程中正发挥着越来越重要的作用。“全国地热资源调查评价与勘查示范工程”自实施以来,主要开展了京津冀地区地热资源评价以及东南沿海和青藏高原东部干热岩的勘查等工作。这些工作的开展在获得一系列科研突破成果的同时,也解决了大量地热开发的工程问题,积累了大量成果和经验。该文对工程取得的进展及成果进行了总结,对工作开展中出现的问题进行了探讨,可用于指导我国地热资源勘查开发工作。     

羌塘盆地晚侏罗世夏里组磁学参数指示的沉积环境及其成盐意义

羌塘盆地位于青藏高原中部,是一个具有沉积钠盐、钾盐远景的蒸发岩盆地。最近在羌塘盆地上侏罗统夏里组中发现多处具有钾、钠异常的盐泉这一重要的找钾盐线索,然而目前对上侏罗统夏里组的古环境演变仍然不甚了解,限制了对晚侏罗世羌塘盆地的成盐(成钾)条件和可能的层位等重大科学问题的理解。因此本文以热磁分析为主,对羌塘盆地雁石坪地区上侏罗统夏里组沉积物进行高分辨率岩石磁学研究,并讨论其古环境意义。最后,根据成盐(成钾)理论,综合分析羌塘盆地上侏罗统夏里组成盐(成钾)的构造、气候和物源条件。κ-T曲线、磁滞回线和热退磁结果共同表明磁铁矿和赤铁矿是夏里组的主要载磁矿物;磁学参数分析表明夏里组上段赤铁矿含量相对增加,揭示夏里组上段形成于干旱气候环境;综合分析构造、气候和物源条件表明,羌塘盆地上侏罗统夏里组上段最为具备良好的成盐(成钾)条件。