Indicator Simulation Accounting for Multiple-Point Statistics

Geostatistical simulation aims at reproducing the variability of the real underlying phenomena. When nonlinear features or large-range connectivity is present, the traditional variogram-based simulation approaches do not provide good reproduction of those features. Connectivity of high and low values is often critical for grades in a mineral deposit. Multiple-point statistics can help to characterize these features. The use of multiple-point statistics in geostatistical simulation was proposed more than 10 years ago, on the basis of the use of training images to extract the statistics. This paper proposes the use of multiple-point statistics extracted from actual data. A method is developed to simulate continuous variables. The indicator kriging probabilities used in sequential indicator simulation are modified by probabilities extracted from multiple-point configurations. The correction is done under the assumption of conditional independence. The practical implementation of the method is illustrated with data from a porphyry copper mine.     

Weathering: Toward a Fractal Quantifying

Weathering occurs over a wide range of scales. To link features through these scales is a major challenge for interdisciplinary weathering studies. Fractal approach seems to be specially useful for this purpose. We introduce a multistep fractal weathering assessment scheme devoted to extract fractal weathering classifiers from texture analysis of the mineral's image. Our scheme enables to quantitatively estimate the global and local information about the geometry of the weathering pattern. This information is basic to develop geometrical indices of weathering, which can significantly enrich the common qualitative and semiquantitative weathering assessment schemes. To justify the fractal approach, a strong statistical self-similarity has been documented for both the weathering and fresh features of two common silica minerals: quartz and biogenic A-opal (phytolith) over four orders of length scales. The procedure is fast, drastically reduces thresholding bias, promises to be universal, it is valid for genetically different minerals and rock types, scale independent, and specially useful for monitoring the changes in the mineral's roughness during the alteration. Two of the proposed classifiers seem to be potentially useful for direct application in the field and be used by nonspecialist.     

Applicability of a Backprojection Algorithm to Reconstruct Images of Subsurface Horizontal Planes for Laboratory Experiments in Electrical Resistance Tomography

This paper describes the use of a backprojection algorithm to reconstruct subsurface images of the electrical resistivity in horizontal planes parallel to the surface. The algorithm can be applied to detect buried objects such as tanks or pipes and possible leakages from them. Two imaging strategies are compared: juxtaposition of vertical planes, and 3D reconstruction from the sensitivity matrix corresponding to the entire volume whose surface is explored. The electrode arrays used for voltage measurement are the dipole–dipole array and a modified Schlumberger array. A personal computer controls current injection, electrode switching, and voltage detection. The system injects 1 kHz, 20 V peak-to-peak square waveforms, thus avoiding electrode polarization effects. Experimental laboratory measurements show that the algorithm detects localized objects such as an insulating sphere and a conductive cylinder immersed in water. Furthermore, covering half of the cylinder by a rubber sleeve to simulate a nonconductive leak, yields a distinct image for the leak. The backprojection algorithm does not need any regularization parameter and it is very fast in inverting the sensitivity matrix because it approximates the inverse matrix by its transposed. The dipole–dipole array usually yields a lower overall pixel error than the modified Schlumberger array but both allow the detection of simulated underground leaks.     

Genesis of Kanggur Gold Deposit in Eastern Tianshan Orogenic Belt, NW China: Fluid Inclusion and Oxygen and Hydrogen Isotope Constraints

Abstract. Primary fluid inclusions in quartz and carbonates from the Kanggur gold deposit are dominated by aqueous inclusions, with subsidiary CO₂-H₂O inclusions that have a constant range in CO₂ content (10–20 vol %). Microthermometric results indicate that total homogenization temperatures have a wide but similar range for both aqueous inclusions (120 to 310C) and CO₂-H₂O inclusions (140 to 340C). Estimates of fluid salinity for CO₂-H₂O inclusions are quite restricted (5.9∼10.3 equiv. wt% NaCl), whereas aqueous inclusions show much wider salinity ranging from 2.2 to 15.6 equivalent wt %NaCl.
The 6D values of fluid inclusions in carbonates vary from -45 to -61 %, in well accord with the published δD values of fluid inclusions in quartz (-46 to -66 %). Most of the δ¹⁸O and δD values of the ore-forming fluids can be achieved by exchanged meteoric water after isotopic equilibration with wall rock by fluid/rock interaction at a low water/rock ratio. However, the exchanged meteoric water alone cannot explain the full range of δ¹⁸O and δD values, magmatic and/or meta-morphic water should also be involved. The wide salinity in aqueous inclusions may also result from mixing of meteoric water and magmatic and/or metamorphic water.     

Methane-rich fluid inclusions in skarn near the giant REE–Nb–Fe deposit at Bayan Obo, Northern China

We report fluid inclusion data for skarn, formed at the contact between Hercynian granitoids and dolomite of the Proterozoic Bayan Obo Group, in the vicinity of Bayan Obo REE–Nb–Fe deposit, Inner Mongolia, China. Three types of fluid inclusions are identified: two-phase CH₄-rich, three-phase liquid–vapour–solid and two-phase aqueous inclusions. Using microthermometry and laser Raman microprobe analysis to calculate isochores for CH₄-bearing inclusions, we estimate fluid trapping conditions at T=280 to 344 °C and P<1 to 2.3 kbar. Such conditions are compatible with formation of CH₄ inclusions as a result of reaction between graphite in the country rocks (black slate sequence) and fluids derived from magma. The lack of carbonaceous material in the inclusions supports the hypothesis that CH₄ was generated during fluid migration rather than by in situ reaction. In contrast to the skarn, and despite the fact that similar graphite-bearing slates are found in the host rocks, no CH₄-bearing inclusions have been so far reported from Bayan Obo REE ores. We therefore conclude that the skarn-forming fluids in the contact aureole of the Hercynian granitoids were not involved at any stage in the formation of the Bayan Obo deposit.     

Transformation of Residuals to Avoid Artifacts in Geostatistical Modelling with a Trend

Trend modelling is an important part of natural resource characterization. A common approach to account for a variable with a trend is to decompose it into a relatively smoothly varying trend and a more variable residual component. Then, the residuals are stochastically modelled independent of the trend. This decomposition can result in values outside the plausible range of variability, such as grades below zero or ratios that exceed 1.0. We transform the residuals conditional to the trend component to explicitly remove these complex features prior to geostatistical modelling. Back transformation of the modelled residual values allows the complex relations to be reproduced. A petroleum-related application shows the robustness of the proposed transformation. Furthermore, a mining application shows that when this conditional transformation is applied to the original variable, instead of the residual, simulated values are assured to be nonnegative.     

北京延庆北部地区后城组火山岩岩石化学特征与成因探讨

根据TAS图解 ,北京延庆北部地区后城组火山属于玄武粗安岩 流纹岩 (和粗面岩 )的双峰式组合。该套火山岩不属于典型的钾玄岩系列 ,在基性岩和中性岩中总体上以碱性系列和向碱性系列过渡的高钾钙碱性系列、钾玄岩系列为主 ,流纹岩类属于富钾质的拉斑系列。该套火山岩轻稀土元素明显富集。玄武粗安岩类岩石没有负铕异常或负铕异常不明显 ,而流纹岩存在明显的负铕异常。延庆北部地区后城组火山岩形成于燕山运动早期或中期阶段的造山变形事件趋于结束时的后碰撞构造环境中。     

An attenuation relationship based on Turkish strong motion data and iso-acceleration map of Turkey

This paper presents an attenuation relationship of peak ground acceleration (PGA) derived from Turkish strong motion data for rock, soil and soft soil sites and an iso-acceleration map of Turkey based on this relationship. For the purpose, among all the three-component accessible records, 221 records from 122 earthquakes that occurred in Turkey between 1976 and November 2003 were selected. The database was compiled for earthquakes with moment magnitudes (M_w) and PGA values ranging between 4.1 and 7.5, and 20 and 806 gal, and distances to epicenter considered in the database were between 5 and 100 km. From the regression analysis of the data, an attenuation equation of PGA considering rock, soil and soft soil conditions was developed. The PGA values predicted from the equation suggested in this study and those both from a few domestic equations and some imported equations were compared. In addition, an iso-acceleration map of Turkey was constructed using the suggested attenuation equation and considering both known active faults and epicenter locations of the earthquakes that have occurred in Turkey.     

复杂环境条件下深基坑综合支护措施的应用

北京某交易中心一期基坑工程属大型深基坑工程，各种地下管线和地上建(构)筑物复杂，介绍了针对其复杂环境条件的综合支护措施的设计思路与实际应用效果。     

管井降水技术在多种含水介质基坑工程中的应用

结合蚌埠城市中心广场会展中心工程，介绍了管井降水技术在多种含水介质的复杂地质环境下的应用。