共查询到20条相似文献，搜索用时 15 毫秒
1.
Numerical representations of multivariate natural phenomena, including characteristics of mineral deposits, petroleum reservoirs
and geo-environmental attributes, need to consider and reproduce the spatial relationships between correlated attributes of
interest. There are, however, only a few methods that can practically jointly simulate large size multivariate fields. This
paper presents a method for the conditional simulation of a non-Gaussian vector random field directly on block support. The
method is derived from the group sequential simulation paradigm and the direct block simulation algorithm which leads to the
efficient joint simulation of large multivariate datasets jointly and directly on the block support. This method is a multistage
process. First, a vector random function is orthogonalized with minimum/maximum autocorrelation factors (MAF). Blocks are
then simulated by performing LU simulation on their discretized points, which are later back-rotated and averaged to yield
the block value. The internal points are then discarded and only the block value is stored in memory to be used for further
conditioning through a joint LU, resulting in the reduction of memory requirements. The method is termed direct block simulation
with MAF or DBMAFSIM. A proof of the concept using an exhaustive data set demonstrates the intricacies and the performance
of the proposed approach. 相似文献
2.
Geostatistical Simulation of Regionalized Pore-Size Distributions Using Min/Max Autocorrelation Factors 总被引：1，自引：0，他引：1
In many fields of the Earth Sciences, one is interested in the distribution of particle or void sizes within samples. Like many other geological attributes, size distributions exhibit spatial variability, and it is convenient to view them within a geostatistical framework, as regionalized functions or curves. Since they rarely conform to simple parametric models, size distributions are best characterized using their raw spectrum as determined experimentally in the form of a series of abundance measures corresponding to a series of discrete size classes. However, the number of classes may be large and the class abundances may be highly cross-correlated. In order to model the spatial variations of discretized size distributions using current geostatistical simulation methods, it is necessary to reduce the number of variables considered and to render them uncorrelated among one another. This is achieved using a principal components-based approach known as Min/Max Autocorrelation Factors (MAF). For a two-structure linear model of coregionalization, the approach has the attractive feature of producing orthogonal factors ranked in order of increasing spatial correlation. Factors consisting largely of noise and exhibiting pure nugget–effect correlation structures are isolated in the lower rankings, and these need not be simulated. The factors to be simulated are those capturing most of the spatial correlation in the data, and they are isolated in the highest rankings. Following a review of MAF theory, the approach is applied to the modeling of pore-size distributions in partially welded tuff. Results of the case study confirm the usefulness of the MAF approach for the simulation of large numbers of coregionalized variables. 相似文献
3.
To speed up multivariate geostatistical simulation it is common to transform the set of attributes into spatially uncorrelated factors that can be simulated independently. Spatial decorrelation methods are usually based on the diagonalisation of the variance/covariance and semivariogram matrices of the set of attributes for a chosen family of lag spacings. These matrices are symmetric and there are several efficient methods for the approximate joint diagonalisation of a family of symmetric matrices. One of these is the uniformly weighted exhaustive diagonalisation with Gauss iterations (U-WEDGE) method. In contrast to the method of minimum/maximum autocorrelation factors (MAF), where a two structure linear model of coregionalisation is assumed, U-WEDGE can be applied directly to the set of experimental semivariogram matrices without having to place restrictions on the number of structures in the linear model of coregionalisation, thus removing one of the restrictions placed on the subsequent modelling of the spatial structure of the factors. We use an iron-ore data set to illustrate the method and present a comparison between the simulated attributes obtained from U-WEDGE and MAF with the full co-simulation of the attributes. 相似文献
4.
5.
We present a new model updating technique for the automatic calculation of the small strain shear modulus of geomaterials, based on bender element experiments. The technique searches for the shear modulus that maximizes the correlation between the output signal obtained from the experiment and the output signal acquired from its computational simulation. Instead of conventional extremum finding algorithms, a fixed point technique is used to update the shear modulus at each iteration. This option increases substantially the attraction basin of the absolute maximum correlation and improves the convergence of the algorithm. 相似文献
6.
7.
8.
9.
Pierre Petitgas Mathieu Woillez Mathieu Doray Jacques Rivoirard 《Mathematical Geosciences》2018,50(2):187-208
Marine research survey data on fish stocks often show a small proportion of very high-density values, as for many environmental data. This makes the estimation of second-order statistics, such as the variance and the variogram, non-robust. The high fish density values are generated by fish aggregative behaviour, which may vary greatly at small scale in time and space. The high values are thus imprecisely known, both in their spatial occurrence and order of magnitude. To map such data, three indicator-based geostatistical methods were considered, the top-cut model, min–max autocorrelation factors (MAF) of indicators, and multiple indicator kriging. In the top-cut and MAF approaches, the variable is decomposed into components and the most continuous ones (those corresponding to the low and medium values) are used to guide the mapping. The methods are proposed as alternatives to ordinary kriging when the variogram is difficult to estimate. The methods are detailed and applied on a spatial data set of anchovy densities derived from a typical fish stock acoustic survey performed in the Bay of Biscay, which show a few high-density values distributed in small spatial patches and also as solitary events. The model performances are analyzed by cross-validating the data and comparing the kriged maps. Results are compared to ordinary kriging as a base case. The top-cut model had the best cross-validation performance. The indicator-based models allowed mapping high-value areas with small spatial extent, in contrast to ordinary kriging. Practical guidelines for implementing the indicator-based methods are provided. 相似文献
10.
A compositional multivariate approach was used to analyse regional-scale soil geochemical data obtained as part of the Tellus Project generated by the Geological Survey of Northern Ireland. The multi-element total concentration data presented comprise X-ray fluorescence (XRF) analyses of 6862 rural soil samples collected at 20-cm depth on a non-aligned grid at one site per \(2\,\hbox {km}^{2}\). Censored data were imputed using published detection limits. Each soil sample site was assigned to the regional geology map, resulting in spatial data for one categorical variable and 35 continuous variables comprised of individual and amalgamated elements. This paper examines the extent to which soil geochemistry reflects the underlying geology or superficial deposits. Since the soil geochemistry is compositional, log-ratios were computed to adequately evaluate the data using multivariate statistical methods. Principal component analysis (PCA) and minimum/maximum autocorrelation factors (MAF) were used to carry out linear discriminant analysis (LDA) as a means to discover and validate processes related to the geologic assemblages coded as age bracket. Peat cover was introduced as an additional category to measure the ability to predict and monitor fragile ecosystems. Overall prediction accuracies for the age bracket categories were 68.4 % using PCA and 74.7 % using MAF. With inclusion of peat, the accuracy for LDA classification decreased to 65.0 and 69.9 %, respectively. The increase in misclassification due to the presence of peat may reflect degradation of peat-covered areas since the creation of superficial deposit classification. 相似文献
11.
12.
河北平原浅层地下水地球化学演化模拟 总被引：2，自引：0，他引：2
反应路径模拟技术是近十几年发展起来的定量研究地下水地球化学演化的重要方法。文中讨论了反应路径模拟与质量平衡模拟的差异，提出了对复杂水流系统采取分段模拟，以及对浅层成水形成的主要作用即蒸发浓缩作用和海侵的影响采用土壤盐类综合体的溶解来表征的处理方法，论述了研究区浅层地下水的地球化学演化过程。 相似文献
13.
14.
Obtaining accurate geological boundaries and assessing the uncertainty in these limits are critical for effective ore resource and reserve estimation. The uncertainty in the extent of an ore body can be the largest source of uncertainty in ore resource estimation when drilling is sparse. These limits are traditionally interpreted deterministically and it can be difficult to quantify uncertainty in the boundary and its impact on ore tonnage. The proposed methodology is to consider stochastic modeling of the ore boundary with a distance function recoding of the available data. This technique is modified to incorporate non-stationarities in the form of a locally varying anisotropy field used in kriging and sequential Gaussian simulation. Implementing locally varying anisotropy kriging retains the geologically realistic features of a deterministic model while allowing for a stochastic assessment of uncertainty. A case study of a gold deposit in Northern Canada is used to demonstrate the methodology. The proposed technique generates realistic, curvilinear geological boundary models and allows for an assessment of the uncertainty in the model. 相似文献
15.
ADINA软件是基于有限元技术的大型通用分析仿真平台,尤其针对结构非线性、流固耦合等复杂问题的求解具有强大优势,广泛应用于各个工业领域的工程仿真计算以及科学研究等各个领域。在机械领域的研究表明,有限元分析建模中的特征技术,对实现设计和有限元分析的紧密耦合,提高产品设计的质量、效率及自动化程度具有重要的意义。通过对几种种特征建模类型的分析和阐述,指出了ADINA建模技术的局限性,为其在各领域的应用提供参考。 相似文献
16.
In a probabilistic analysis of rock slope stability, the Monte Carlo simulation technique has been widely used to evaluate
the probability of slope failure. While the Monte Carlo simulation technique has many advantages, the technique requires complete
information of the random variables in stability analysis; however, in practice, it is difficult to obtain complete information
from a field investigation. The information on random variables is usually limited due to the restraints of sampling numbers.
This is why approximation methods have been proposed for reliability analyses. Approximation methods, such as the first-order
second-moment method and the point estimate method, require only the mean and standard deviation of the random variable; therefore,
it is easy to utilize when the information is limited. Usually, a single closed form of the formula for the evaluation of
the factor of safety is needed for an approximation method. However, the commonly used stability analysis method of wedge
failure is complicated and cumbersome and does not provide a simple equation for the evaluation of the factor of safety. Consequently,
the approximation method is not appropriate for wedge failure. In order to overcome this limitation, a simple equation, which
is obtained from the maximum likelihood estimation method for wedge failure, is utilized to calculate the probability of failure.
A simple equation for the direct estimation of the safety factors for wedge failure has been empirically derived from failed
and stable cases of slope, using the maximum likelihood estimation method. The developed technique has been applied to a practical
example, and the results from the developed technique were compared to the results from the Monte Carlo simulation technique. 相似文献
17.
考虑饱和-非饱和渗流、温度和应力耦合的三维有限元程序研制 总被引：13，自引：3，他引：10
在探索多场耦合理论基础上,研制出了一个考虑饱和-非饱和渗流场和温度场对应力场作用的三维弹塑性有限元程序,该程序采用了以下简化假设：不考虑力学过程和渗流过程对温度场产生的影响,而只考虑变温对力学过程产生的影响;不考虑力学过程和温度场对渗流场产生的影响,而只考虑渗流场对力学过程产生的影响。经与弹性力学解析解对比,验证了该程序计算结果是合理的。 相似文献
18.
Analysis of non climatic origins of floods in the downstream part of the Kura River,Azerbaijan 总被引：1，自引：1，他引：0
Over the past century, there is an increased contribution of non climatic factors to the flood formation processes in the
Kura River. Non climatic factors of floods refer to factors that are related to reductions in channel capacity and result
in floods. More recently, there are numerous non climatic factors occurring in and around the Kura River basin that have increased
the frequency of floods. Sediment accumulation in the riverbed over a long period of time has led to the reduction of channel
capacity and has raised the elevation of the riverbed above the surrounding territory. It is illustrated that construction
of dykes and levees do not actually prevent flooding, where hydrologic connections between groundwater and surface water are
high, since infiltrated waters from channel results in raising of ground waters, causing an effect of “underground flood.”
Since underground floods occur when water going from channels raises the level of ground waters, there is an urgent need to
carefully investigate the groundwater–surface water connections. With the purpose of predicting floods, the authors suggest
defining maximal acceptable flows (MAFs) rather than channel capacities. Results show that high rates of hydraulic conductivity
of soils will decrease MAF rates. MAF computations before high-water season allow for further regulation of outlets further
downstream in order to prevent flooding and enable flood forecasting. While the study focuses on a specific region, the overall
approach suggested is generic and may be applied elsewhere. 相似文献
19.
This paper presents a numerical technique for the analysis of the cone penetration test by means of the commercial finite element code ABAQUS. The von Mises yield criterion with its associated flow rule is assumed to model the plastic behaviour of elastoplastic undrained clays. An explicit finite element scheme is used to efficiently perform a large number of loading increments and to simplify the treatment of contact. An Arbitrary Langrangian–Eulerian (ALE) scheme is adopted to preserve the quality of mesh throughout the numerical simulation. A volumetric weighting algorithm adjusts the relative positions of nodes after each loading increment. This prevents mesh over distortion and allows the simulation to run continuously. The variation of the cone resistance is examined in relation to various parameters such as the in situ stress state, shaft and cone face roughness, and the material strength when steady state conditions have been reached. The trends of these variations are highlighted and compared with those found by other researchers. This technique can be extended to analyse the plastic behaviour of elastoplastic sands often modelled using either the Drucker–Prager yield criterion or a critical state model. 相似文献