共查询到20条相似文献，搜索用时 437 毫秒
1.
Conditional Simulation with Patterns 总被引：17，自引：0，他引：17
An entirely new approach to stochastic simulation is proposed through the direct simulation of patterns. Unlike pixel-based
(single grid cells) or object-based stochastic simulation, pattern-based simulation simulates by pasting patterns directly
onto the simulation grid. A pattern is a multi-pixel configuration identifying a meaningful entity (a puzzle piece) of the
underlying spatial continuity. The methodology relies on the use of a training image from which the pattern set (database)
is extracted. The use of training images is not new. The concept of a training image is extensively used in simulating Markov
random fields or for sequentially simulating structures using multiple-point statistics. Both these approaches rely on extracting
statistics from the training image, then reproducing these statistics in multiple stochastic realizations, at the same time
conditioning to any available data. The proposed approach does not rely, explicitly, on either a statistical or probabilistic
methodology. Instead, a sequential simulation method is proposed that borrows heavily from the pattern recognition literature
and simulates by pasting at each visited location along a random path a pattern that is compatible with the available local
data and any previously simulated patterns. This paper discusses the various implementation details to accomplish this idea.
Several 2D illustrative as well as realistic and complex 3D examples are presented to showcase the versatility of the proposed
algorithm. 相似文献
2.
Fast FILTERSIM Simulation with Score-based Distance 总被引：5，自引：3，他引：2
FILTERSIM is a pattern-based multiple-point geostatistical algorithm for modeling both continuous and categorical variables.
It first groups all the patterns from a training image into a set of pattern classes using their filter scores. At each simulation
location, FILTERSIM identifies the training pattern class closest to the local conditioning data event, then samples a training
pattern from that prototype class and pastes it onto the simulation grid. In the original FILTERSIM algorithm, the selection
of the closest pattern class is based on the pixel-wise distance between the prototype of each training pattern class and
the local conditioning data event. Hence, FILTERSIM is computationally intensive for 3D simulations, especially with a large
and pattern-rich training image. In this paper, a novel approach is proposed to accelerate the simulation process by replacing
that pixel-wise distance calculation with a filter score comparison, which is the difference between the filter score of local
conditioning data event and that of each pattern prototype. This score-based distance calculation significantly reduces the
CPU consumption due to the tremendous data dimension reduction. The results show that this new score based-distance calculation
can speed up FILTERSIM simulation by a factor up to 10 in 3D applications. 相似文献
3.
Snehamoy Chatterjee Roussos Dimitrakopoulos Hussein Mustapha 《Mathematical Geosciences》2012,44(3):343-374
A pattern-based simulation technique using wavelet analysis is proposed for the simulation (wavesim) of categorical and continuous variables. Patterns are extracted by scanning a training image with a template and then storing
them in a pattern database. The dimension reduction of patterns in the pattern database is performed by wavelet decomposition
at certain scale and the approximate sub-band is used for pattern database classification. The pattern database classification
is performed by the k-means clustering algorithm and classes are represented by a class prototype. For the simulation of categorical variables,
the conditional cumulative density function (ccdf) for each class is generated based on the frequency of the individual categories at the central node of the template. During
the simulation process, the similarity of the conditioning data event with the class prototypes is measured using the L
2-norm. When simulating categorical variables, the ccdf of the best matched class is used to draw a pattern from a class. When continuous variables are simulated, a random pattern
is drawn from the best matched class. Several examples of conditional and unconditional simulation with two- and three- dimensional
data sets show that the spatial continuity of geometric features and shapes is well reproduced. A comparative study with the
filtersim algorithm shows that the wavesim performs better than filtersim in all examples. A full-field case study at the Olympic Dam base metals deposit, South Australia, simulates the lithological
rock-type units as categorical variables. Results show that the proportions of various rock-type units in the hard data are
well reproduced when similar to those in the training image; when rock-type proportions between the training image and hard
data differ, the results show a compromise between the two. 相似文献
4.
Fast direct sampling for multiple-point stochastic simulation 总被引：1，自引：0，他引：1
Multiple-point statistics simulation has recently attracted significant attention for the simulation of complex geological structures. In this paper, a fast direct sampling (FDS) algorithm is presented based on a fast gradient descent pattern matching strategy. The match is directly extracted from the training image (TI) and so the method does not require intensive preprocessing and database storage. The initial node of the search path is selected randomly but the following nodes are selected in a principled manner so that the path is conducted to the right match. Each node is selected based on the matching accuracy and the behavior of the TI in the previous node. A simple initialization strategy is presented in this paper which significantly accelerates the matching process at the expense of a very naïve preprocessing stage. The proposed simulation algorithm has several outstanding advantages: it needs no (or very limited) preprocessing, does not need any database storage, searches for the match directly in the TI, is not limited to fixed size patterns (the pattern size can be easily changed during simulation), is capable of handling both continuous and categorical data, is capable of handling multivariate data, and finally and more importantly, is a fast method while maintaining high standards for the matching quality. Experiments on different TIs reveal that the simulation results of FDS and DS are comparable in terms of pattern reproduction and connectivity while FDS is far faster than DS. 相似文献
5.
Weights of evidence modeling for combining indicator patterns in mineral resource evaluation is based on an application of Bayes' rule. Two weights are defined for each indicator pattern and Bayes' rule is applied repeatedly to combine indicator patterns. If all patterns are conditionally independent with respect to deposits, the logit of the posterior probability can be calculated as the sum of the logit of the prior probability plus the weights of the overlay patterns. The information to be integrated for gold exploration in Xiong-er Mountain Region comes from a geological map, an interpreted map of a Thematic Mapper (TM) image, and the locations of known gold deposits. Favorable stratigraphic units, structural control factors, and alteration factors are considered. The work was conducted on an S600 I2S image-processing system. FORTRAN programs were developed for creating indicator patterns, statistical calculations, and pattern integration. Six indicator patterns were selected to predict mineral potential. They are conditionally independent according to pairwiseG
2 tests, and an overall chi-square test. The potential area predicted using the 32 known deposits generally coincides with the prospect areas determined by geological fieldwork. 相似文献
6.
以随机函数理论为基础,采用相控-多参数协同的随机建模方法,建立塔河油气田AT1区块凝析气藏三维地质模型,实现气藏精细三维表征。首先,以钻井和岩芯资料为基础构建储层构造模型;然后,以小层界面为控制条件建立储层结构模型;接着,在沉积相、地质条件的约束下,采用序贯指示模拟法来建立砂体骨架模型;随后,在砂体骨架模型内进行优势相计算,形成最终有效砂体骨架模型;最后,以有效砂体骨架模型为约束,采用序贯高斯模拟法建立储层物性参数模型。结果表明：将物性参数变量与微相分布结合的序贯高斯模拟法建立孔隙度等物性参数的分布模型,以及采用地质分析类比、地质统计分析等方法优选最佳模型是有效的地质建模方法;所建地质模型精确细致地表征了塔河油气田AT1区块凝析气藏构造格架及储层、流体三维分布,反映了辫状水道复合连片,东北向展布,储层物性受相控较明显,非均质性较强。 相似文献
7.
Luis Manuel de Vries Jesus Carrera Oriol Falivene Oscar Gratacós Luit Jan Slooten 《Mathematical Geosciences》2009,41(1):29-42
Simulation of flow and solute transport through aquifers or oil reservoirs requires a precise representation of subsurface
heterogeneity that can be achieved by stochastic simulation approaches. Traditional geostatistical methods based on variograms,
such as truncated Gaussian simulation or sequential indicator simulation, may fail to generate the complex, curvilinear, continuous
and interconnected facies distributions that are often encountered in real geological media, due to their reliance on two-point
statistics. Multiple Point Geostatistics (MPG) overcomes this constraint by using more complex point configurations whose
statistics are retrieved from training images. Obtaining representative statistics requires stationary training images, but
geological understanding often suggests a priori facies variability patterns. This research aims at extending MPG to non-stationary
facies distributions. The proposed method subdivides the training images into different areas. The statistics for each area
are stored in separate frequency search trees. Several training images are used to ensure that the obtained statistics are
representative. The facies probability distribution for each cell during simulation is calculated by weighting the probabilities
from the frequency trees. The method is tested on two different object-based training image sets. Results show that non-stationary
training images can be used to generate suitable non-stationary facies distributions. 相似文献
8.
The time and temperature evolution of twinning in cordierite is simulated using three computer models. The orientation of
walls between twin domains in natural cordierite follows mainly the ferroelastic pattern which minimises the strain energy
of the walls between twin-related domains. Such ferroelastic twinning is simulated in an elastic three-states Potts model
in which each structural six-membered ring is represented by a three state pseudo-spin. The resulting twin pattern in a sample
with 3169 structural rings shows sector trilling and fine scale ferroelastic wall patterns which coarsen with increasing annealing
time. The poorly defined wall directions observed in cordierite were found to be related to twin walls which do not minimise
the strain energy. Instead, these walls are located along the corners of pseudo-hexagonal rings and appear as the consequence
of local rather than global interatomic interactions. Simulations using two-dimensional (38028 atoms) and three-dimensional
(408 228 atoms) structural models show a predominance of these topological walls over the strain walls at early stages in
the ordering process. The domain structure in the simulation is patchy rather than corresponding to repeated stripe structures
found in other ferroelastic and co-elastic materials. In all models, a strong tendency for sector trilling is observed. In
kinetic tweed patterns a novel 60o tweed is found at atomic length scales while the usual strain-mediated 90o tweed appears
at mesoscopic length scales. An unusual surface tension effect in domain formation and ’writhing’, fluid-like motion was found
in the three-dimensional structural model. This motion, along with the existence of non strain-mediated walls may contribute
to cordierite’s poorly defined domain wall directions at the early stages of domain coarsening.
Received: 11 March 1998 / Revised, accepted: 28 August 1998 相似文献
9.
Simulation of Geological Contacts from Interpreted Geological Model Using Multiple-Point Statistics 总被引：2，自引：1，他引：1
Alexandre Boucher Joao Felipe Costa Luis Gustavo Rasera Eduardo Motta 《Mathematical Geosciences》2014,46(5):561-572
Applications of multiple-point statistics (mps) algorithms to large non-repetitive geological objects such as those found in mining deposits are difficult because most mps algorithms rely on pattern repetition for simulation. In many cases, an interpreted geological model built from a computer-aided design system is readily available but suffers as a training image due to the lack of patterns repetitiveness. Porphyry copper deposits and iron ore formations are good examples of such mining deposits with non-repetitive patterns. This paper presents an algorithm called contactsim that focuses on reproducing the patterns of the contacts between geological types. The algorithm learns the shapes of the lithotype contacts as interpreted by the geologist, and simulates their patterns at a later stage. Defining a zone of uncertainty around the lithological contact is a critical step in contactsim, because it defines both the zones where the simulation is performed and where the algorithm should focus to learn the transitional patterns between lithotypes. A larger zone of uncertainty results in greater variation between realizations. The definition of the uncertainty zone must take into consideration the geological understanding of the deposit, and the reliability of the contact zones. The contactsim algorithm is demonstrated on an iron ore formation. 相似文献
10.
11.
Solving Speed and Memory Issues in Multiple-Point Statistics Simulation Program SNESIM 总被引：2，自引：2，他引：0
In the last 10 years, Multiple-Point Statistics (MPS) modeling has emerged in Geostatistics as a valuable alternative to traditional variogram-based and object-based modeling. In contrast to variogram-based simulation, which is limited to two-point correlation reproduction, MPS simulation extracts and reproduces multiple-point statistics moments from training images; this allows modeling geologically realistic features, such as channels that control reservoir connectivity and flow behavior. In addition, MPS simulation works on individual pixels or small groups of pixels (patterns), thus does not suffer from the same data conditioning limitations as object-based simulation. The Single Normal Equation Simulation program SNESIM was the first implementation of MPS simulation to propose, through the introduction of search trees, an efficient solution to the extraction and storage of multiple-point statistics moments from training images. SNESIM is able to simulate three-dimensional models; however, memory and speed issues can occur when applying it to multimillion cell grids. Several other implementations of MPS simulation were proposed after SNESIM, but most of them manage to reduce memory demand or simulation time only at the expense of data conditioning exactitude and/or training pattern reproduction quality. In this paper, the original SNESIM program is revisited, and solutions are presented to eliminate both memory demand and simulation time limitations. First, we demonstrate that the time needed to simulate a grid node is a direct function of the number of uninformed locations in the conditioning data search neighborhood. Thus, two improvements are proposed to maximize the ratio of informed to uniformed locations in search neighborhoods: a new multiple-grid approach introducing additional intermediary subgrids; and a new search neighborhood designing process to preferentially include previously simulated node locations. Finally, because SNESIM memory demand and simulation time increase with the size of the data template used to extract multiple-point statistics moments from the training image and build the search tree, a simple method is described to minimize data template sizes while preserving training pattern reproduction quality. 相似文献
12.
V. I. Znak 《Mathematical Geology》2005,37(2):207-221
Median filters are nonlinear and a theoretical analysis of their behavior is very difficult and so are rarely used for the processing of seismic data. However, they are able to preserve steps, sharp discontinuities and edges that are lost using most other standard filters. As seismic data are mostly harmonic or frequency modulated signals in the frequency range 1 to 10 Hz, the median filter must be adapted so as to suppress noise but not unduly distort the signal. This can be accomplished by using a median filter whose length is n times that of signal period where n is even. By use of weighted-order statistics, samples of the signals closest in phase can be obtained. To illustrate the method, a signal from a frequency band of 6.3–7.5 Hz was processed and the quality of the signal enhanced two fold over the quality of the signal without processing. 相似文献
13.
G. Agrosì G. C. Capitani E. Scandale G. Tempesta 《Physics and Chemistry of Minerals》2011,38(2):101-109
Transmission Electron Microscopy (TEM) and High Resolution TEM (HRTEM) investigations were used to study a complex micro-structure
of twin-related lamellae of 15R SiC polytype, in syntactic coalescence with 6H-SiC. The analysis of the local stacking microstructure by means of HRTEM revealed that the perfect structure of 15R polytype, (23)3 was locally interrupted by numerous adjoining stacking faults parallel to (0001) with stacking of the coupled (22) and (33)
bilayers superimposed on the twin boundaries. HR images taken exactly across both the twin boundaries showed a zig-zag pattern
(23) that switched to (32) by a twin coherent interfaces or, alternatively, a zig-zag pattern (32) which passed to (23) through
an isolated (33), 6H like sequence. The selected area electron diffraction (SAED) patterns taken exactly above both the twin interfaces indicate
classifying of the twin found in this study as a “Friedelian” reticular merohedric twinning. However, two indistinguishable
twin operations matched the observed features: a reflection through rational plane (0001), and 180°-rotation around [0001].
Since individual Si and C atoms and even the SiC bilayer polarity could not be established from these HR images, the real
twin law was deduced by taking into account that the coherent structural match at the interface could be guaranteed only by
a 180°-rotation around [0001]. In order to explain the origin of the planar defects found in this sample, the growth mechanism
and the influence of the low-energy stacking faults were considered. 相似文献
14.
The evolution of a subduction channel and orogenic wedge is simulated in 2D for an active continental margin, with P-T paths being displayed for selected markers. In our simulation, subduction erosion affects the active margin and a structural pattern develops within a few tens of millions of years, with four zones from the trench into the forearc: (i) an accretionary complex of low grade metamorphic sedimentary material, (ii) a wedge of nappes with alternating upper and lower crustal provenance, and minor interleaving of oceanic or hydrated mantle material, (iii) a megascale melange composed of high pressure (HP) and ultra-high pressure (UHP) metamorphic rocks extruded from the subduction channel, and (iv) the upward tilted frontal part of the remaining lid. The P–T paths and time scales correspond to those typically recorded in orogenic belts. The simulation shows that HP/UHP metamorphism of continental crust does not necessarily indicate collision, but that the material can be derived from the active margin by subduction erosion and extruded from the subduction channel beneath the forearc during ongoing subduction. 相似文献
15.
区域尺度上含水层非均质具有复杂的结构性和随机性,难以准确刻画,造成非均质对区域地下水流系统的影响机制研究不够深入。本文以鄂尔多斯盆地白垩系地下水流系统为研究实例,选择典型剖面,采用剖面二维随机数值模拟方法,通过对比不同非均质刻画方法下地下水流场的变化,探讨含水层层状非均质对地下水流系统的影响机制。结果显示,均质条件下模型各向异性（含水层水平和垂向渗透系数比值Kh/Kv）取值为1000时,地下水流场与实际条件较为接近;非均质条件下,渗透系数方差取值0.91,水平相关长度取值5000 m,Kh/Kv取值150时,接近实际条件。研究表明,在大尺度地下水流模拟研究中,采用水平相关长度、渗透系数方差和各向异性值三个变量生成的随机场能很好地刻画含水层的层状非均质特征及其对水流系统的影响控制作用。由于含水层不同尺度层状非均质的叠加效应,采用均质各向异性介质等效概化含水层层状非均质性会造成等效各向异性值偏大失真的效应。 相似文献
16.
Iwao Kawabe 《Mathematical Geology》1977,9(1):39-53
A stochastic model of the magmatic differentiation by fractional crystallization is given. The probability distribution of the quantity of crystallized solid removed from silicate melt was assumed to be binomial at any stage of the differentiation. According to this model, when concentrations of an element are transformed into their powers by using the reciprocal of the difference between bulk-partition coefficient for the element and unity as an exponent, the resultant frequency-distribution pattern becomes absolutely normal. The patterns of concentrations, however, cannot be expressed by a particular type of function but are variable according to the values of bulk-partition coefficients. Frequency distributions of element concentrations were computed under selected conditions on the basis of the model. The result shows that the skewed pattern often observed in frequency distributions of minor-element concentrations is explained by this fractional crystallization process. The frequency distributions of Ni and Cr concentrations in the geosynclinal basalt of southwestern Japan were examined in terms of this model. It was concluded that the model can be applied to the formation of the basalt concerning at least these two elements. 相似文献
17.
Memory-Efficient Categorical Multi-point Statistics Algorithms Based on Compact Search Trees 总被引：1，自引：1，他引：0
Tuanfeng Zhang Stein Inge Pedersen Christen Knudby David McCormick 《Mathematical Geosciences》2012,44(7):863-879
Multi-point statistics (MPS) has emerged as an advanced geomodeling approach. A practical MPS algorithm named snesim (simple normal equations simulation), which uses categorical-variable training images, was proposed in 2001. The snesim algorithm generates a search tree to store the occurrence statistics of all patterns in the training image within a given set of search templates before the simulation proceeds. The snesim search tree concept makes MPS simulation central processing unit efficient but consumes large amounts of memory, particularly when three-dimensional training images contain complex patterns and when a large search template is required to ensure optimal reproduction of the image patterns. To crack the memory-restriction bottleneck, we have developed a compact search tree that contains the same information but reduces memory cost by one order of magnitude. Furthermore, the compact structure also accelerates MPS simulation significantly. Such remarkable improvement makes MPS a more practical tool to use in building the large and complex three-dimensional facies models required in the oil and gas industry. 相似文献
18.
Because modeling the flow of gas and oil in fractured media is increasingly important, multipoint geostatistics approaches, such as single normal equation simulation and FILTERSIM, have become a focus of research in the simulation of sedimentary faces. However, there are problems when these methods are applied to the simulation of discrete fracture networks (DFNs). The biggest problem is that the generated fractures are discontinuous, behaving as separated and isolated points. This paper proposes an improved approach we call CCL-FILTERSIM, which combines FILTERSIM with a run-based, connected-component labeling algorithm to generate more continuous fractures. This paper focuses on the continuity of the simulated fracture network and adopts the scan of the training image, classification of the patterns, and the sequential simulation proposed by FILTERSIM. The principle innovation of the CCL-FILTERSIM approach is the choice of data pattern. Only patterns that form better continuous object patterns are chosen. This is realized by a connected-component labeling algorithm that calculates the connected runs of selected patterns and chooses the pattern with the least runs. Some cases are presented to compare realizations from CCL-FILTERSIM and FILTERSIM, which show that CCL-FILTERSIM yields more continuous objects than FILTERSIM. The DFNs generated by this method provide better conditions for subsequent fluid flow research. 相似文献
19.
One of the common sources of error in triaxial tests is the penetration of the latex membrane into the peripheral voids of the specimen as confining pressure is applied. Existing analytical solutions developed to predict the magnitude of this membrane penetration have all assumed that the peripheral sand particles in contact with the surrounding membrane are arranged so that the centroids of the sand–membrane contact areas form a square pattern of side length D50 on a plane parallel to the undeflected membrane. Measurements made on digital images of sand–membrane contact patterns for triaxial specimens of varying characteristics have shown this assumption to be inadequate. The analysis presented herein shows how actual contact pattern can be reasonably replicated by simulated patterns. These simulated patterns are generated by adding randomly selected distances to the centroid co‐ordinates of contact areas initially placed in a systematic grid pattern. The results indicate that it is possible to systematically model sand–membrane contact patterns as a function of triaxial specimen characteristics such as relative density and preparation method. The ability to simulate the real patterns eliminates the reliance of existing analytical solutions on a single assumed square contact pattern and hence can serve as the basis for improving the accuracy of future solutions to account for the effects of membrane penetration. Copyright © 2000 John Wiley & Sons, Ltd. 相似文献
20.
离散裂隙随机网络模型中结构面的分布特征仅代表其统计规律,并非地质成因、结构模式、构造形迹的客观反映,以致后续力学计算及渗流模拟结果可信度偏低。针对国家高放废物处置库甘肃北山预选区出露良好的花岗岩体,基于实测结构面统计规律、内在成生关系及水力联系,在前期纯随机模型关键部位修正特征参数,并融入人工辨识的确定结构模式而构建“随机-确定”耦合模型。图形和渗流模拟两种定量检验结果表明：耦合模型结构面数量更接近客观实际,模型准确度提升约48.8%;耦合模型渗流路径与流量较之前更显客观真实,模拟结果与试验数据的接近程度比随机模型大约1/3。另外,不同孔位渗流结果显示：相比于随机结构面,确定性结构面对区域渗流控制作用更加明显,在渗流模拟中扮演更为重要角色。此耦合模型有益于拓宽结构面网络模拟的发展方向。 相似文献