Combined mine tremors source location and error evaluation in the Lubin Copper Mine (Poland)

A modified method of mine tremors location used in Lubin Copper Mine is presented in the paper. In mines where an intensive exploration is carried out a high accuracy source location technique is usually required. The effect of the flatness of the geophones array, complex geological structure of the rock mass and intense exploitation make the location results ambiguous in such mines. In the present paper an effective method of source location and location's error evaluations are presented, combining data from two different arrays of geophones. The first consists of uniaxial geophones spaced in the whole mine area. The second is installed in one of the mining panels and consists of triaxial geophones. The usage of the data obtained from triaxial geophones allows to increase the hypocenter vertical coordinate precision. The presented two-step location procedure combines standard location methods: P-waves directions and P-waves arrival times. Using computer simulations the efficiency of the created algorithm was tested. The designed algorithm is fully non-linear and was tested on the multilayered rock mass model of the Lubin Copper Mine, showing a computational better efficiency than the traditional P-wave arrival times location algorithm. In this paper we present the complete procedure that effectively solves the non-linear location problems, i.e. the mine tremor location and measurement of the error propagation.     

Crustal structure at the easternmost termination of the Variscan belt based on CELEBRATION 2000 and ALP 2002 data

The eastern margin of the Variscan belt in Europe comprises plate boundaries between continental blocks and terranes formed during different tectonic events. The crustal structure of that complicated area was studied using the data of the international refraction experiments CELEBRATION 2000 and ALP 2002. The seismic data were acquired along SW–NE oriented refraction and wide-angle reflection profiles CEL10 and ALP04 starting in the Eastern Alps, passing through the Moravo-Silesian zone of the Bohemian Massif and the Fore-Sudetic Monocline, and terminating in the TESZ in Poland. The data were interpreted by seismic tomographic inversion and by 2-D trial-and-error forward modelling of the P waves. Velocity models determine different types of the crust–mantle transition, reflecting variable crustal thickness and delimiting contacts of tectonic units in depth. In the Alpine area, few km thick LVZ with the Vp of 5.1 km s^− 1 dipping to the SW and outcropping at the surface represents the Molasse and Helvetic Flysch sediments overthrust by the Northern Calcareous Alps with higher velocities. In the Bohemian Massif, lower velocities in the range of 5.0–5.6 km s^− 1 down to a depth of 5 km might represent the SE termination of the Elbe Fault Zone. The Fore-Sudetic Monocline and the TESZ are covered by sediments with the velocities in the range of 3.6–5.5 km s^− 1 to the maximum depth of 15 km beneath the Mid-Polish Trough. The Moho in the Eastern Alps is dipping to the SW reaching the depth of 43–45 km. The lower crust at the eastern margin of the Bohemian Massif is characterized by elevated velocities and high Vp gradient, which seems to be a characteristic feature of the Moravo-Silesian. Slightly different properties in the Moravian and Silesian units might be attributed to varying distances of the profile from the Moldanubian Thrust front as well as a different type of contact of the Brunia with the Moldanubian and its northern root sector. The Moho beneath the Fore-Sudetic Monocline is the most pronounced and is interpreted as the first-order discontinuity at a depth of 30 km.     

激光拉曼光谱法分析多种显微组分荧光变化及其应用

利用激发488 nm激光的氩离子激光器作为RENISHAW inVia型激光拉曼光谱仪的光源,建立了多种显微组分荧光变化(FAMM)分析方法,并对镜质体反射率明显抑制的东营凹陷有效烃源岩进行了测定。结果表明,东营凹陷有效烃源岩的有机质类型越好,镜质体反射率抑制程度越高,有效烃源岩的真实成熟度应主要处于0.64%~1.30%,而不是实测镜质体反射率所反映的0.37%~1.10%。     

西藏驱龙超大型斑岩铜矿床:地质、蚀变与成矿

驱龙超大型矿床是一个产于后碰撞伸展环境下、与大洋俯冲无关的新型斑岩铜矿。文章通过对驱龙铜矿床地质、蚀变与矿化的详细研究,建立了驱龙中新世岩浆演化序列,初步查明了岩浆浅成侵位的构造控制要素,厘定了主要的围岩蚀变类型及空间展布规律,查明了引起各期蚀变事件的地质记录及矿化的空间分布规律,并探讨了成矿物质沉淀的机制,初步建立了该矿床的成矿模型。研究表明,驱龙铜矿中新世斑岩是闪长质深部岩浆房不断演化的产物,花岗闪长岩中新发现的、结晶时间为22.2Ma左右的闪长质包体可近似代表深部岩浆房组分,依次产出的花岗闪长岩、呈岩株或岩枝产出的P斑岩、X斑岩及最晚期的闪长玢岩(15.7±0.2)Ma,均为深部岩浆房连续演化的产物,岩浆持续6Ma左右。岩浆演化过程中角闪石、斜长石不断的结晶分异,导致了岩石常量元素、稀土元素及微量元素组成的规律性变化,斑岩埃达克质的特征也因岩浆演化过程中角闪石等矿物的不断结晶分异而引起。X斑岩中锆石的Hf同位素特征表明,岩石可能形成于新生下地壳的部分熔融。大面积产出的花岗闪长岩为驱龙铜矿最主要的含矿围岩,容纳了驱龙矿床70%以上的矿体,主要由斜长石、钾长石和石英组成,具花岗结构-似斑状结构,近EW向产出,其浅成就位可能受背斜控制,其后的各期斑岩均沿该侵位中心上侵,而冈底斯地壳中新世的快速抬升与剥蚀是导致含矿斑岩浅成侵位的根本原因;矿区内的SN向裂隙带既不控岩,也不控矿。浅成侵位的斑岩及深部岩浆房均发生了流体出溶。发生了大量流体出溶的深部岩浆房,是矿区早期蚀变流体的主要来源,显微晶洞构造及单向固结结构(UST)是流体出溶的地质记录。蚀变主要有3种类型,分别为早期的钾硅酸盐化、青磐岩化以及晚期的长石分解。钾硅酸盐化可分为2个阶段,即蚀变矿物以次生钾长石为主的早期钾硅酸盐化和以次生黑云母为主的晚期钾硅酸盐化。青磐岩化因产出的岩石类型不同,蚀变矿物组合具有明显差异性:产于叶巴组地层中的青磐岩化相对较强,蚀变矿物以绿帘石为主;产于花岗闪长岩中的青磐岩化相对较弱,蚀变矿物以绿泥石为主。晚期长石分解蚀变以破坏长石类矿物为特征,蚀变矿物主要为绢云母-绿泥石-粘土等。石英和硬石膏贯穿于上述各种蚀变中。空间上,钾硅酸盐化位于斑岩体及其周围地区,青磐岩化位于钾硅酸岩化外侧。后期形成的长石分解蚀变强烈叠加了早期钾硅酸盐化,介于钾硅酸盐化带与青磐岩化带之间。与早期钾长石化有关的脉体主要为不规则石英-钾长石脉,与晚期黑云母化有关的脉体主要为不规则至板状的石英-硬石膏脉、黑云母脉,与青磐岩化有关的脉体主要为板状的绿帘石-石英脉,与晚期长石分解蚀变有关的脉体主要为板状黄铜矿-黄铁矿脉及黄铁矿脉;在早期钾硅酸盐蚀变与晚期长石分解蚀变转换阶段,发育一组板状的石英-硫化物脉。早期不规则的脉体形成于斑岩结晶早期、矿区裂隙小规模发育阶段;晚期的板状脉体形成于斑岩弱固结或固结之后、矿区大规模连通裂隙发育阶段。驱龙矿区的铜矿化分布较为均一,主体产于花岗闪长岩中,其中,铜矿化主体形成于黑云母化蚀变阶段,转变阶段及长石分解阶段也有大量铜的形成;钼主要形成于转换阶段,长石分解蚀变阶段也有产出。黑云母化阶段,铜的沉淀与角闪石黑云母化、斜长石钾长石化过程中Ca2 的大量释放有关;转换阶段,铜钼矿化可能与压力和(或)温度骤降有关;晚期铜矿化与长石矿化蚀变阶段,斜长石绿泥石化、黑云母绿帘石化过程中Ca2 及Fe2 的释放有关。     

豫西前河构造蚀变岩型金矿成矿过程中的流体-岩石反应

前河金矿区位于华北地台南缘,赋存在中元古界熊耳群安山岩和英安岩中,矿体受断裂破碎带控制。含矿热液在迁移过程中与围岩发生了广泛的流体-岩石反应而引起热液蚀变。本区石英中有4种类型的流体包裹体,均一温度范围为145～331℃,其中含CO2包裹体的完全均一温度主要分布在中-高温区。成矿流体的密度和压力变化范围分别是0.68～0.94g/cm3和(367.01～896.55)×105Pa。金大量沉淀成矿时的流体特征为:温度213～260℃、密度0.80～0.89g/cm3和压力(502.86～710.57)×105Pa。流体相为SO42->Na >Cl->K 型,CO2/H2O比值及N2、H2S、Ar、C2H6等挥发分的含量明显增高,f(CO2)、f(H2S)、f(CH4)和Eh值增大;f(O2)、f(H2O)和pH值减小。在青磐岩化安山岩的基础上发生的流体-岩石反应是造成本矿床金沉淀成矿的最主要原因。     

Three-dimensional discontinuous deformation analysis (3-D DDA) using a new contact resolution algorithm

This paper presents a new contact calculating algorithm for contacts between two polyhedra with planar boundaries in the three-dimensional discontinuous deformation analysis (3-D DDA). In this algorithm, all six type contacts in 3-D (vertex-to-face, vertex-to-edge, vertex-to-vertex, face-to-face, edge-to-edge, and edge-to-face) are simply transformed into the form of point-to-face contacts. The presented algorithm is a simple and efficient method and it can be easily coded into a computer program. In this paper, formulations of normal contact, shear contact and frictional force submatrices based on the new method are derived and the algorithm has been programmed in VC⁺⁺. Examples are provided to demonstrate the new contact rule between two blocks.     

A comparative study of four vulnerability mapping methods in a detritic aquifer under mediterranean climatic conditions

This paper presents the results of a comparative study relating to the application of four vulnerability mapping methods, GOD, AVI, DRASTIC and SINTACS, in a pilot detritic aquifer situated in NW Morocco, known as the Martil–Alila aquifer. The principal objective of this work is to determine the most suitable such methods for this aquifer type within a Mediterranean context, and to show the effect of the rainfall variations that are characteristic of the Mediterranean climate on the degree of vulnerability. The methods applied distinguish five classes of vulnerability, these being irregularly divided up in space, with the division varying according to the method in question. The vulnerability maps obtained by the different methods strongly suggest that the eastern half of the aquifer is more vulnerable to contamination than the western half, for all hydrological situations. The effect of climatic conditions on the degree of vulnerability is well represented by the DRASTIC, according to which the aquifer is moderately to strongly vulnerable during humid hydrological years and weakly to moderately vulnerable during dry ones. For the other methods, this climatic effect is limited to the area occupied by the two predominant classes (“High” and “Low” for GOD and “High” and “Moderate” for SINTACS) while it is null for AVI. In conclusion, DRASTIC appears the most suitable for mapping the vulnerability to contamination of Mediterranean coastal detritic aquifers such as the Martil–Alila aquifer.     

Mineralogy and geochemistry of the Texeo Cu–Co mine site (NW Spain): screening tools for environmental assessment

The Cu–Co–Ni Texeo mine has been the most important source of Cu in NW Spain since Roman times and now, approximately 40,000 m³ of wastes from mine and metallurgical operations, containing average concentrations of 9,263 mg kg⁻¹ Cu, 1,100 mg kg⁻¹ As, 549 mg kg⁻¹ Co, and 840 mg kg⁻¹ Ni, remain on-site. Since the cessation of the activity, the abandoned works, facilities and waste piles have been posing a threat to the environment, derived from the release of toxic elements. In order to assess the potential environmental pollution caused by the mining operations, a sequential sampling strategy was undertaken in wastes, soil, surface and groundwater, and sediments. First, screening field tools were used to identify hotspots, before defining formal sampling strategies; so, in the areas where anomalies were detected in a first sampling stage, a second detailed sampling campaign was undertaken. Metal concentrations in the soils are highly above the local background, reaching up to 9,921 mg kg⁻¹ Cu, 1,373 mg kg⁻¹ As, 685 mg kg⁻¹ Co, and 1,040 mg kg⁻¹ Ni, among others. Copper concentrations downstream of the mine works reach values up to 1,869 μg l⁻¹ and 240 mg kg⁻¹ in surface water and stream sediments, respectively. Computer-based risk assessment for the site gives a carcinogenic risk associated with the presence of As in surface waters and soils, and a health risk for long exposures; so, trigger levels of these elements are high enough to warrant further investigation.     

Determination of historical emission of heavy metals into the atmosphere: Celje case study

The research shows that in the Celje area (Slovenia), the historical anthropogenical emissions are 1,712 tons of Zn and 9.1 tons of Cd. For Zn, this value represents approximately 0.3% of the total Zn production in that area. Close to the former zinc smelting plant, the “Zn precipitation” has been estimated to be up to 0.036 mm. The 100-year Zn production left behind a heavily contaminated area with maximum concentrations of Zn of up to 5.6% in attic dust and 0.85% in the soil, and 456 mg/kg of Cd in attic dust and 59.1 mg/kg in the soil. The calculation of historical emissions is based on the data of heavy metals concentration in the attic dust at 98 sampling points and on the data from 19 measurement sites of the weight of total monthly air deposit. The main idea behind determining past emissions is that when the weight of the deposited dust on a small area is multiplied by the concentration of the element in that area, the mass of the polluter which has been transported to the place of interest by air can be calculated. If we sum up all the weight over the whole geochemical anomaly, we get the quantity of historical emissions.