Delineation of potential exploration targets based on 3D geological modeling: A case study from the Laoangou Pb-Zn-Ag polymetallic ore deposit,China

The Laoangou Pb-Zn-Ag polymetallic ore deposit (LAG) is a hydrothermal vein-type deposit in the East Qinling polymetallic metallogenic belt, which is controlled by NW-, NE- and nearly SN-trending faults and mainly hosted in the dolomite marble of Meiyaogou and Baishugou formations within the Luanchuan ore district (LOD). Although district-scale metallogenic prediction related to the Pb-Zn-Ag deposits in the LOD has been attempted, there have been no studies to formulate a detailed deposit-scale metallogenic prediction of the Pb-Zn-Ag deposit in 3D space. Here we selected the LAG to formulate a metallogenic prediction model for Pb-Zn-Ag deposit based on data from 1:10,000 scale geological map, and 66 boreholes and 39 sections of exploration lines on 1:1000 scale, to extract the salient spatial features of the deposit. We apply ordinary weights and weighted weights of evidence, followed by boosted weights of evidence, logistic regression and information entropy for integrating the features of the ore deposit for exploration targeting. The C–V fractal method is applied to classify the probabilities. Accordingly, three levels of exploration targets are delineated: the first level targets are mainly distributed in the periphery of known orebodies (M2, M3, II-3 and VI), which are also located in the ore-controlling strata and fault zones; the second level targets are located in the region of known orebodies, indirectly confirming that the prediction result is reliable; and the third level targets, which are mainly located at the contact between ore-controlling strata and metagabbro bodies, offer new targets for exploration in this ore district. Combining our results obtained in this study with those in previous studies, we provide some guidelines for exploration targets in the Pb-Zn-Ag deposits in the LOD, which might help in more effectively delineating the target zones.     

Pile driving formulas based on pile wave equation analyses

Pile driving formulas, which directly relate the pile set resulting from a hammer blow to the static load capacity of the pile, are often used to decide whether a pile will have the required design capacity. However, existing formulas do not consider soil or pile type, and do not, in general, reliably predict pile capacity. In this paper, an advanced model for dynamic pile driving analysis was used to develop accurate pile driving formulas. The proposed driving formulas are validated through well-documented case histories. Comparisons of predictions from the proposed formulas with results from static and dynamic load tests show that they produce reasonably accurate predictions of pile capacity based on pile set observations.     

“Transverse Action” – A model benchmark exercise for numerical analysis of the Callovo-Oxfordian claystone hydromechanical response to excavation operations

The Callovo-Oxfordian claystone (COx) is considered as a potential geological formation to host an industrial radioactive waste repository in France. A detailed understanding of the thermo-hydro-mechanical (THM) behavior of the COx is a key issue for design of different repository structures and the safety calculations of the project. More particularly, numerical modeling of induced fracture networks around drifts excavated at the main level of the Andra’s Meuse/Haut-Marne Underground Research Laboratory (M/HM URL) and short and long term behavior of the COx around these drifts are of great interest. Several constitutive models have been developed/used in the framework of the R&D and simulation programs of Andra. A model benchmark exercise has been launched since 2012 to provide an overall view of the developed models regarding the in situ observations. In this view, two series of test cases, respectively at material point scale and at drift excavation scale are defined. Different kinds of constitutive models based on the elasto-visco-plasticity concept, continuum damage mechanics, the rigid block spring method and two-scale computational homogenized model (CHM) are used within this exercise. The obtained results show that accounting for material anisotropy and strain localization treatment techniques can improve the obtained results when elasto-visco-plastic models are used. Damage mechanics based approaches and methods accounting for discontinuities through discrete elements provide also interesting insights especially when fracturing processes must be modeled. However, more efforts are necessary to improve the robustness of these kinds of approaches in the complex context of COx response to excavation works.     

Numerical 3D modeling of the effects of strain rate and confining pressure on the compressive behavior of Kuru granite

This paper deals with numerical modeling of the compressive behavior of granite rock under high strain rate dynamic loading and wide range of confining pressure. For this end, a constitutive model based on damage mechanics and viscoplasticity for rock is formulated and implemented in explicit dynamics FEM. The viscoplastic part is based on a simple power law type yield criterion that incorporates the rate-dependency with a linear viscosity term. Moreover, a Rankine type of tensile cut-off is employed. The damage part of the model is formulated with separate scalar damage variables in tension and compression. The model is calibrated for Kuru granite and validated with the experimental data from dynamic compression tests at the strain rate of 600 1/s up to 225 MPa of confining pressure. The numerical simulations demonstrate that, despite the underlying continuum modeling approach, the model captures the correct experimental failure modes, including the transition from single-to-multiple fragmentation, as well as the dynamic compressive strengths at different confining pressures.     

Cone penetration test in sand: A numerical-analytical approach

Separation of the effects of initial horizontal stress and relative density on cone tip resistance in sandy soils has been a complicated issue for many years. In order to overcome this problem, a numerical modeling of CPT which has been verified by calibration chamber tests, has been used in this paper to achieve a reliable analytical solution. The analytical solution has resulted in two relationships for sleeve friction and cone tip resistance in terms of the initial conditions of sandy soil. Based on the presented solution, the initial horizontal stress and relative density can be determined according to CPT measurements.     

塔里木盆地玉北地区埋藏史及热史分析

根据地层分层数据和古热流分析,对塔里木盆地玉北地区的5口代表性钻井进行了埋藏史、热史的模拟与分析,探讨了研究区烃源岩的热演化成熟过程,初步划分了烃源岩的成熟阶段和生烃期次。研究结果认为：玉北地区寒武系－泥盆系沉积时期,具有多次沉积埋藏和隆升剥蚀的交替;泥盆系沉积之后,石炭－二叠系呈现沉积埋藏的特征;中生代在玉北地区东部少量沉积三叠系,随后被剥蚀,而在中西部发生沉积间断;进入新生代,地层持续沉积,且沉积速率呈现增大趋势。研究区寒武系烃源岩自沉积以来成熟度演化相对连续,奥陶纪时进入生烃门限,成熟时期较早;志留－泥盆纪,受海西早期运动的影响,发生隆升剥蚀作用,烃源岩演化变缓;石炭－二叠纪受控于区域上持续稳定的地层沉积与广泛发育的火山活动,研究区寒武系烃源岩进入中高成熟阶段;由于喜马拉雅造山运动的远程效应,新生代研究区持续沉降,形成巨厚沉积,烃源岩成熟度达到高－过成熟,为天然气生成阶段     

基于谱元法的各向异性黏弹性介质地震波场模拟和分析

基于谱元法计算理论,通过坐标旋转算法实现任意各向异性黏弹性介质中的地震波场计算,研究复合介质中地震波场特征。根据计算结果,qSV波表现为剪切滞弹性形变,qP波表现为膨胀滞弹性形变;品质因子不仅影响地震波振幅,还影响其传播频率;黏弹性介质吸收高频快,在一定频率范围内,黏弹性介质没有频散现象。     

交通地理信息系统中的时空动态分段模型研究

动态分段技术是交通地理信息系统(GIS-T)中一项重要的线性要素动态显示与分析技术,该文提出了一种时空动态分段模型,使动态分段系统的组成成分融入时态因素,并将属性的时态信息作为属性信息表的字段存储,通过时态地理信息系统技术完成物理段的时态变化。利用统一建模语言(UML)的类图设计方式开发了物理模型,描述了时空动态分段模型中对象之间的关系,并定义了必要的属性和操作;同时基于线性参照系统(LRS)表达空间实体的方法,概括了模型中所涉及到的主要拓扑关系。实验表明,时空动态分段模型弥补其他模型将时空参考分为时间参考和空间参考所产生的不足,而且更明确地将多重属性和物理实体有机结合起来。     

基于交错网格Fourier伪谱微分矩阵算子的地震波场模拟GPU加速方案

作为高精度波形反演或逆时偏移的重要组成部分,地震波数值模拟对计算速度和效率提出了更高要求.GPU通用计算技术的产生及其内在数据并行性,为高效地震波数值模拟应用和研究得以有效开展奠定了基础.本文借助交错网格的Fourier伪谱微分矩阵算子和GPU上高效矩阵乘法,实现了复杂介质中地震波模拟的高效算法.数值试验表明,优化后的GPU计算相比CPU单核计算在大规模二维地震波场计算中获得至少100x以上的加速比.这对我们快速分析目标反射层在地震剖面中同相轴位置,制定优化采集方案具有重要意义.     

A mathematical model for steady‐state regolith production at constant erosion rate

It has been hypothesized that many soil profiles reach a steady‐state thickness. In this work, such profiles were simulated using a one‐dimensional model of reaction with advective and diffusive solute transport. A model ‘rock’ is considered, consisting of albite that weathers to kaolinite in the presence of chemically inert quartz. The model yields three different steady‐state regimes of weathering. At the lowest erosion rates, a local‐equilibrium regime is established where albite is completely depleted in the weathering zone. This regime is equivalent to the transport‐limited regime described in the literature. With an increase in erosion rate, transition and kinetic regimes are established. In the transition regime, both albite and kaolinite are present in the weathering zone, but albite does not persist to the soil–air interface. In the weathering‐limited regime, here called the kinetic regime, albite persists to the soil–air interface. The steady‐state thickness of regolith decreases with increasing erosion rate in the local equilibrium and transition regimes, but in the kinetic regime, this thickness is independent of erosion rate. Analytical expressions derived from the model are used to show that regolith production rates decrease exponentially with regolith thickness. The steady‐state regolith thickness increases with the Darcy velocity of the pore fluid, and in the local equilibrium regime may vary markedly with small variations in this velocity and erosion rate. In the weathering‐limited regime, the temperature dependences for chemical weathering rates are related to the activation energy for the rate constant for mineral reaction and to the ΔH of dissolution, while for local equilibrium regimes they are related to the ΔH only. The model illustrates how geochemical and geomorphological observations are related for a simple compositional system. The insights provided will be useful in interpreting natural regolith profiles. Copyright © 2010 John Wiley & Sons, Ltd.