Improving Sequential Simulation with a Structured Path Guided by Information Content

Multiple-point simulation is a newly developed geostatistical method that aims at combining the strengths of two mainstream geostatistical methods: object-based and pixel-based methods. It maintains the flexibility of pixel-based algorithms in data conditioning, while enhancing its capability of reproducing realistic geological shapes, which is traditionally reserved to object-based algorithms. However, the current snesim program for multiple-point simulation has difficulty in reproducing large-scale structures, which have a significant impact on the flow response. To address this problem, we propose to simulate along a structured path based on an information content measure. This structured path accounts for not only the information from the data, but also some prior structural information provided by geological knowledge. Various case studies show a better reproduction of large-scale structures. This concept of simulating along a structured path guided by information content can be applied to any sequential simulation algorithms, including traditional variogram-based two-point geostatistical algorithms.     

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.     

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.     

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.     

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

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

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

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

合理设计岩屑规程保证正常钻进条件——金刚石钻进技术专题之一

论述了正常钻进的条件，提出了合理岩屑规程是保证正常钻进的关键和设计钻头的主要依据之一，讨论了岩屑规程与金刚石钻头唇面工作状态、孔底表面粗糙度和岩屑粒度的关系，提出了相应的建议。     

HJ-1复合型钻井液的研究及应用

分析和探讨了HJ-1复合型钻井液对复杂地层护壁的作用机理及护壁性能试验情况，介绍了其在寨上矿区复杂地层钻探中的成功应用。