共查询到20条相似文献，搜索用时 46 毫秒
1.
Corentin Faucher Antoine Saucier Denis Marcotte 《Stochastic Environmental Research and Risk Assessment (SERRA)》2013,27(1):253-273
We present a new stochastic simulation method that builds two-dimensional images by assembling together square image pieces called blocks. The blocks are taken from a reference image. Our method, called patchwork simulation method (PSM), enforces pattern continuity in the image. Moreover, PSM allows to control the image local-mean histogram. This histogram bin-frequencies can be set to user-defined target values that may differ from the reference image local-mean histogram. This flexibility enhances the PSM generality by enlarging the set of all possible simulations. The local-mean histogram control is achieved by adjusting suitably the transition probabilities that associate a new block to an existing neighborhood in the partly simulated image. For several types of synthetic images and one polymer blend image, we show that PSM reproduces faithfully the reference image visual appearance (i.e. patterns are correctly shaped) and that simulated images are statistically compatible with the target local-mean histogram. Moreover, we show that our method has the ability to produce simulations that respect conditional hard data as well as a target local-mean histogram. 相似文献
2.
Subsurface flow and solute transport simulations are performed using different scenarios of permeability fields generated from the sequential Gaussian simulation method (SGS), the multiple-point FILTERSIM algorithm and a new multiple-point wavelet-based simulation method (SWS). The SWS method is a multiple-point pattern-based simulation method which uses discrete wavelet transformation for the representation of geologic heterogeneity. For pattern-based simulation, patterns are generated by scanning a training image with a spatial template. The pattern classifications were performed after reducing the dimension of patterns by wavelet decomposition at the suitable scale and by taking only scaling components of wavelet decomposed patterns. The simulation is performed in a sequential manner by finding the best-matched class corresponding to the conditioning data and by randomly sampling a pattern from the best-matched class. The developed method is compared with two other multi-point simulation algorithms, FLTERSIM and SIMPAT. The comparative results revealed that the proposed method is computationally faster than the other two methods while the simulation maps are comparable. Numerical simulations of two flow problems are performed using SGS, SWS and FILTERSIM realizations. The numerical results show a superiority of the SWS method over SGS and FILTERSIM in terms of reproduction of the reference images main features, and agreement with flow and transport results obtained on reference images. 相似文献
3.
S. De Iaco M. Palma 《Stochastic Environmental Research and Risk Assessment (SERRA)》2002,16(5):333-341
In geostatistics, stochastic simulation is often used either as an improved interpolation algorithm or as a measure of the
spatial uncertainty. Hence, it is crucial to assess how fast realization-based statistics converge towards model-based statistics
(i.e. histogram, variogram) since in theory such a match is guaranteed only on average over a number of realizations. This
can be strongly affected by the random number generator being used. Moreover, the usual assumption of independence among simulated
realizations of a random process may be affected by the random number generator used. Simulation results, obtained by using
three different random number generators implemented in Geostatistical Software Library (GSLib), are compared. Some practical
aspects are pointed out and some suggestions are given to users of the unconditional LU simulation method. 相似文献
4.
Multiple-point geostatistical simulation is used to simulate the spatial structures of geological phenomena. In contrast to conventional two-point variogram based geostatistical methods, the multiple-point approach is capable of simulating complex spatial patterns, shapes, and structures normally observed in geological media. A commonly used pattern based multiple-point geostatistical simulation algorithms is called FILTERSIM. In the conventional FILTERSIM algorithm, the patterns identified in training images are transformed into filter score space using fixed filters that are neither dependent on the training images nor on the characteristics of the patterns extracted from them. In this paper, we introduce two new methods, one for geostatistical simulation and another for conditioning the results. At first, new filters are designed using principal component analysis in such a way to include most structural information specific to the governing training images resulting in the selection of closer patterns in the filter score space. We then propose to combine adaptive filters with an overlap strategy along a raster path and an efficient conditioning method to develop an algorithm for reservoir simulation with high accuracy and continuity. We also combine image quilting with this algorithm to improve connectivity a lot. The proposed method, which we call random partitioning with adaptive filters simulation method, can be used both for continuous and discrete variables. The results of the proposed method show a significant improvement in recovering the expected shapes and structural continuity in the final simulated realizations as compared to those of conventional FILTERSIM algorithm and the algorithm is more than ten times faster than FILTERSIM because of using raster path and using small overlap specially when we use image quilting. 相似文献
5.
Characterization of shale reservoirs, which are typically of low permeability, is very difficult because of the presence of multiscale structures. While three-dimensional (3D) imaging can be an ultimate solution for revealing important complexities of such reservoirs, acquiring such images is costly and time consuming. On the other hand, high-quality 2D images, which are widely available, also reveal useful information about shales’ pore connectivity and size. Most of the current modeling methods that are based on 2D images use limited and insufficient extracted information. One remedy to the shortcoming is direct use of qualitative images, a concept that we introduce in this paper. We demonstrate that higher-order statistics (as opposed to the traditional two-point statistics, such as variograms) are necessary for developing an accurate model of shales, and describe an efficient method for using 2D images that is capable of utilizing qualitative and physical information within an image and generating stochastic realizations of shales. We then further refine the model by describing and utilizing several techniques, including an iterative framework, for removing some possible artifacts and better pattern reproduction. Next, we introduce a new histogram-matching algorithm that accounts for concealed nanostructures in shale samples. We also present two new multiresolution and multiscale approaches for dealing with distinct pore structures that are common in shale reservoirs. In the multiresolution method, the original high-quality image is upscaled in a pyramid-like manner in order to achieve more accurate global and long-range structures. The multiscale approach integrates two images, each containing diverse pore networks – the nano- and microscale pores – using a high-resolution image representing small-scale pores and, at the same time, reconstructing large pores using a low-quality image. Eventually, the results are integrated to generate a 3D model. The methods are tested on two shale samples for which full 3D samples are available. The quantitative accuracy of the models is demonstrated by computing their morphological and flow properties and comparing them with those of the actual 3D images. The success of the method hinges upon the use of very different low- and high-resolution images. 相似文献
6.
Shu-Wei Gan Shou-Dong Wang Yang-Kang Chen Jiang-Long Chen Wei Zhong Cheng-Lin Zhang 《应用地球物理》2016,13(1):127-134
Conventional f?x empirical mode decomposition (EMD) is an effective random noise attenuation method for use with seismic profiles mainly containing horizontal events. However, when a seismic event is not horizontal, the use of f?x EMD is harmful to most useful signals. Based on the framework of f?x EMD, this study proposes an improved denoising approach that retrieves lost useful signals by detecting effective signal points in a noise section using local similarity and then designing a weighting operator for retrieving signals. Compared with conventional f?x EMD, f?x predictive filtering, and f?x empirical mode decomposition predictive filtering, the new approach can preserve more useful signals and obtain a relatively cleaner denoised image. Synthetic and field data examples are shown as test performances of the proposed approach, thereby verifying the effectiveness of this method. 相似文献
7.
根据材料内部真实结构建立数值计算模型是岩土力学研究的一个热点。本文针对不同材料内部结构特征,介绍了等效图形法与边缘矢量化两种建模方法,并以自制孔隙岩土材料和混凝土的CT图像为例,说明了两种建模步骤。分析表明:当图像的灰度直方图具有明显双峰或多峰结构时,特征区域易分割;而当灰度直方图为单峰结构时,可以通过数学形态学的top-hat方法进行图像增强,以方便特征区域的提取。对于孔隙结构,可以利用等效圆或等效椭圆表征,建立相应的数值模型;而当内部结构复杂时,可利用特征边缘矢量化方法进行建模。本研究为如何建立反映岩土材料内部真实结构数值模型提供了方法。 相似文献
8.
Potential of time‐lapse photography for identifying saturation area dynamics on agricultural hillslopes 
下载免费PDF全文
下载免费PDF全文 Rasmiaditya Silasari Juraj Parajka Camillo Ressl Peter Strauss Günter Blöschl 《水文研究》2017,31(21):3610-3627
Mapping saturation areas during rainfall events is important for understanding the dynamics of overland flow. In this study, we evaluate the potential of high temporal resolution time‐lapse photography for mapping the dynamics of saturation areas (i.e., areas where water is visually ponding on the surface) on the hillslope scale during natural rainfall. We take 1 image per minute over a 100 × 15 m2 depression area on an agricultural field in the Hydrological Open Air Laboratory, Austria. The images are georectified and classified by an automated procedure, using grey intensity as a threshold to identify saturation area. The optimum threshold T is obtained by comparing saturation areas from the automated analysis with the manual analysis of 149 images. T is found to be highly correlated with an image brightness characteristic defined as the greyscale image histogram mode M (Pearson correlation r = 0.91). We estimate T as T = M + C where C is a calibration parameter assumed to be constant during each event. The automated procedure estimates the total saturation area close to the manual analysis with mean normalized root mean square error of 9% and 21% if C is calibrated for each event and taken constant for all events, respectively. The spatial patterns of saturation are estimated with a geometric mean accuracy index of 94% as compared to the manual analysis of the same photos. The patterns are tested against field observations for one date as a preliminary demonstration, which yields a root mean square error of the shortest distance between the measured boundary points and the automatically classified boundary as 23 cm. The usefulness of the patterns is illustrated by exploring run‐off generation processes of an example event. Overall, the proposed classification method based on grey intensity is found to process images with highly varying brightnesses well. It is more efficient than the manual tracing for a large number of images, which allows the exploration of surface flow processes at high temporal resolution. 相似文献
9.
I. V. Nizkous I. A. Sanina E. Kissling L. I. Gontovaya 《Izvestiya Physics of the Solid Earth》2006,42(4):286-296
Arrival times of P and S waves from local earthquakes in the Kamchatka area of the Kurile-Kamchatka Island Arc are used for calculating a spatial model of the elastic wave velocity distribution to a depth of 200 km. The lithosphere is shown to be strongly stratified in its velocity properties and laterally heterogeneous within the mantle wedge and seismic focal zone. A lower velocity layer (an asthenospheric wedge) is identified at depths of 70–130 km beneath the Eastern Kamchatka volcanic belt. The morphology of the Moho interface and the velocity properties of the crust are studied. The main tectonic structures of the region are shown to be closely interrelated with deep velocity heterogeneities. Regular patterns in the statistics of the earthquakes are analyzed in relation to variations in the elastic wave velocities in the focal layer. A mechanism of lithospheric block displacements along weakened zones of the lower crust and upper mantle is proposed. 相似文献
10.
微米X射线计算机断层扫描作为一种数字岩心探测手段,被广泛应用于研究含天然气水合物沉积物赋存形态,但由于水合物与水对X射线的衰减系数相近,二者在CT图像中灰度区间存在交集,导致在从CT图像上对水合物与水进行分割时存在强非唯一性。为提高对CT图像中水合物与水阈值分割的准确性,本文通过分析天然气水合物生长过程中不同时刻CT图像直方图特征,提出归一化CT图像及其直方图的方法。首先选定甲烷气与石英砂的峰值灰度基准;然后用高斯函数分别对当前CT图像直方图中的甲烷气与石英砂曲线进行拟合,得到当前CT图像直方图中的甲烷气与石英砂峰值灰度;再将当前CT图像直方图中的甲烷气峰值灰度与石英砂峰值灰度归一化到选定的峰值灰度基准;进而用归一化的直方图对CT图像进行归一化;最后根据归一化灰度直方图的变化趋势,定量统计得到CT图像中水合物增加和气-水减少的灰度区间,完成图像中不同组分的阈值划分。实验结果表明,提出的阈值分割方法能够为天然气水合物CT图像中水合物与水边界的确定和水合物饱和度计算提供定量的依据,具有实际的工程应用价值。 相似文献
11.
Mathias Raschke 《Journal of Seismology》2013,17(4):1157-1182
We introduce a new approach for ground motion relations (GMR) in the probabilistic seismic hazard analysis (PSHA), being influenced by the extreme value theory of mathematical statistics. Therein, we understand a GMR as a random function. We derive mathematically the principle of area equivalence, wherein two alternative GMRs have an equivalent influence on the hazard if these GMRs have equivalent area functions. This includes local biases. An interpretation of the difference between these GMRs (an actual and a modeled one) as a random component leads to a general overestimation of residual variance and hazard. Beside this, we discuss important aspects of classical approaches and discover discrepancies with the state of the art of stochastics and statistics (model selection and significance, test of distribution assumptions, extreme value statistics). We criticize especially the assumption of logarithmic normally distributed residuals of maxima like the peak ground acceleration (PGA). The natural distribution of its individual random component (equivalent to exp(ε 0) of Joyner and Boore, Bull Seism Soc Am 83(2):469–487, 1993) is the generalized extreme value. We show by numerical researches that the actual distribution can be hidden and a wrong distribution assumption can influence the PSHA negatively as the negligence of area equivalence does. Finally, we suggest an estimation concept for GMRs of PSHA with a regression-free variance estimation of the individual random component. We demonstrate the advantages of event-specific GMRs by analyzing data sets from the PEER strong motion database and estimate event-specific GMRs. Therein, the majority of the best models base on an anisotropic point source approach. The residual variance of logarithmized PGA is significantly smaller than in previous models. We validate the estimations for the event with the largest sample by empirical area functions, which indicate the appropriate modeling of the GMR by an anisotropic point source model. The constructed distances like the Joyner–Boore distance do not work well for event-specific GMRs. We discover also a strong relation between magnitude and the squared expectation of the PGAs being integrated in the geo-space for the event-specific GMRs. One of our secondary contributions is the simple modeling of anisotropy for a point source model. 相似文献
12.
Reconstruction of porous media using multiple-point statistics with data conditioning 总被引：1，自引：0，他引：1
Ting Zhang Yi Du Tao Huang Xue Li 《Stochastic Environmental Research and Risk Assessment (SERRA)》2015,29(3):727-738
The reconstruction of porous media is of great importance in predicting fluid transport properties, which are widely used in various fields such as catalysis, oil recovery, medicine and aging of building materials. The real three-dimensional structural data of porous media are helpful to describe the irregular topologic structures of porous media. By using multiple-point statistics (MPS) to extract the characteristics of real porous media acquired from micro computed tomography (micro-CT) scanning, the probabilities of structural characteristics of pore spaces are obtained first, and then reproduced in the reconstructed regions. One solution to overcome the anisotropy of training images is to use real 3D volume data as a training image (TI). The CPU cost and memory burden brought up by 3D simulations can be reduced greatly by selecting the optimal multiple-grid template size that is determined by the entropy of a TI. Moreover, both soft data and hard data are integrated in MPS simulation to improve the accuracy of reconstructed images. The variograms and permeabilities, computed by lattice Boltzmann method, of the reconstructed images and the target image obtained from real volume data are compared, showing that the structural characteristics of reconstructed porous media using our method are similar to those of real volume data. 相似文献
13.
R.C. Acharya M.I.J. van Dijke K.S. Sorbie S.E.A.T.M. Van der Zee A. Leijnse 《Advances in water resources》2007
We present a 3D network model with particle tracking to upscale 3D Brownian motion of non-reactive tracer particles subjected to a velocity field in the network bonds, representing both local diffusion and convection. At the intersections of the bonds (nodes) various jump conditions are implemented. Within the bonds, two different velocity profiles are used. At the network scale the longitudinal dispersion of the particles is quantified through the coefficient DL, for which we evaluate a number of methods already known in the literature. Additionally, we introduce a new method for derivation of DL based on the first-arrival times distribution (FTD). To validate our particle tracking method, we simulate Taylor’s classical experiments in a single tube. Subsequently, we carry out network simulations for a wide range of the characteristic Péclet number Pe? to assess the various methods for obtaining DL. Using the new method, additional simulations have been carried out to evaluate the choice of nodal jump conditions and velocity profile, in combination with varying network heterogeneity. In general, we conclude that the presented network model with particle tracking is a robust tool to obtain the macroscopic longitudinal dispersion coefficient. The new method to determine DL from the FTD statistics works for the full range of Pe?, provided that for large Pe? a sufficiently large number of particles is used. Nodal jump conditions should include molecular diffusion and allow jumps in the upstream direction, and a parabolic velocity profile in the tubes must be implemented. Then, good agreement with experimental evidence is found for the full range of Pe?, including increased DL for increased porous medium heterogeneity. 相似文献
14.
Properties and limitations of sequential indicator simulation 总被引：2，自引：0，他引：2
The sequential indicator algorithm is a widespread geostatistical simulation technique that relies on indicator (co)kriging and is applicable to a wide range of datasets. However, such algorithm comes up against several limitations that are often misunderstood. This work aims at highlighting these limitations, by examining what are the conditions for the realizations to reproduce the input parameters (indicator means and correlograms) and what happens with the other parameters (other two-point or multiple-point statistics). Several types of random functions are contemplated, namely: the mosaic model, random sets, models defined by multiple indicators and isofactorial models. In each case, the conditions for the sequential algorithm to honor the model parameters are sought after. Concurrently, the properties of the multivariate distributions are identified and some conceptual impediments are emphasized. In particular, the prior multiple-point statistics are shown to depend on external factors such as the total number of simulated nodes and the number and locations of the samples. As a consequence, common applications such as a flow simulation or a change of support on the realizations may lead to hazardous interpretations. 相似文献
15.
Many types of structural systems that undergo cycles of inelastic deformation under severe natural hazard loadings exhibit ‘pinching’ of hysteresis loops. In this paper, a generally pinching hysteretic restoring force model—an extension of the Bouc–Wen differential hysteresis model—is used in stochastic equivalent linearization of single-degree-of-freedom structural systems. The severity and rate of pinching are controlled by the hysteretic energy dissipation and the pinching level can be specified to match experimental data. Under white noise excitations, estimates of reponse statistics from linearization are shown to compare favourably with those from Monte Carlo simulation. Numerical studies on the sensitivity of the accuracy of response statistics obtained by linearization to changes in the hysteresis parameters showed that, for a range of practical cases, the linearization method can be used in lieu of simulation and that, in low-frequency systems, some hysteresis parameters may be set to a constant value a priori to reduce the number of model parameters that needs to be estimated or identified, and to simplify further random vibration analysis and/or performance evaluation studies. 相似文献
16.
17.
特征值相干属性是一种高效、成熟的裂缝预测技术。但即使在连续、稳定的目的层段,部分裂缝隐没于相近高灰度值区域而难以识别;在破碎带发育区域,由于成像能量难以完美聚焦,中、小尺度的裂缝区域呈现相近低灰度值的云雾状模糊而无法区分。本文提出一种以直方图均衡化处理为核心技术的新型裂缝增强方法。该方法能够强化相干图像中不连续性信息与背景的差异性,突出裂缝的线性结构:采用相干图像阈值逐级调节的方式,实现对不同尺度裂缝的预测。该方法同时提升了显性、隐性裂缝识别能力和精度。 相似文献
18.
The statistics of the relative variations (δfoF2) in the critical frequency of the ionospheric F 2 region during catastrophic earthquakes with magnitude M ≥ 6 at 70 vertical sounding (VS) automatic ionospheric stations (AISs) of the global network is considered. Five-day periods including 70 such situations have been selected for 1981 and 1983–1985. Five sets of 24 foF2 measurements from 3 days before to a day after the earthquake instant are formed for each AIS. The number of AISs, for which the fourth statistical invariant is above the significance level (E ≥ 1), has been analyzed. It has been indicated that the histogram deformation at the set of AISs is substantial a day before earthquakes. 相似文献
19.
Testing the correctness of the sequential algorithm for simulating Gaussian random fields 总被引：2，自引：2，他引：0
The sequential algorithm is widely used to simulate Gaussian random fields. However, a rigorous application of this algorithm is impractical and some simplifications are required, in particular a moving neighborhood has to be defined. To examine the effect of such restriction on the quality of the realizations, a reference case is presented and several parameters are reviewed, mainly the histogram, variogram, indicator variograms, as well as the ergodic fluctuations in the first and second-order statistics. The study concludes that, even in a favorable case where the simulated domain is large with respect to the range of the model, the realizations may poorly reproduce the second-order statistics and be inconsistent with the stationarity and ergodicity assumptions. Practical tips such as the multiple-grid strategy do not overcome these impediments. Finally, extending the original algorithm by using an ordinary kriging should be avoided, unless an intrinsic random function model is sought after. 相似文献