首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 546 毫秒
1.
Computer vision provides several tools for analyzing and simulating textures. The principles of these techniques are similar to those in multiple-point geostatistics, namely, the reproduction of patterns and consistency in the results from a perceptual point of view, thus, ensuring the reproduction of long range connectivity. The only difference between these techniques and geostatistical simulation accounting for multiple-point statistics is that conditioning is not an issue in computer vision. We present a solution to the problem of conditioning simulated fields while simultaneously honoring multiple-point (pattern) statistics. The proposal is based on a texture synthesis algorithm where a fixed search (causal) pattern is used. Conditioning is achieved by adding a non-causal search neighborhood that modifies the conditional distribution from which the simulated category is drawn, depending on the conditioning information. Results show an excellent reproduction of the features from the training image, while respecting the conditioning information. Some issues related to the data structure and to the computer efficiency are discussed.  相似文献   

2.
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.  相似文献   

3.
Sequential Gaussian simulation is one of the most widespread algorithms for simulating regionalized variables in the earth sciences. Simplicity and flexibility of this algorithm are the most important reasons that make it popular, but its implementation is highly dependent on a screen effect approximation that allows users to use a moving neighborhood instead of a unique neighborhood. Because of this, the size of the moving neighborhood the number of conditioning data and the size of variogram range are important in the simulation process and should be chosen carefully. In this work, different synthetic and real case studies are presented to show the effect of the neighborhood size the number of conditioning data and the size of variogram range on the simulation result, with respect to the reproduction of the model first and second-order parameters. Results indicate that, in both conditional and non-conditional simulation cases, using a neighborhood with <50 conditioning data may lead to an inaccurate reproduction of the model statistics, and some cases require considering more than 200 conditioning data. It also can be understood from the result of example 3 that when the variogram range is beg compared to the simulation domain determination of inaccurate simulation program is harder.  相似文献   

4.
Gurdak JJ  McCray JE  Thyne G  Qi SL 《Ground water》2007,45(3):348-361
A methodology is proposed to quantify prediction uncertainty associated with ground water vulnerability models that were developed through an approach that coupled multivariate logistic regression with a geographic information system (GIS). This method uses Latin hypercube sampling (LHS) to illustrate the propagation of input error and estimate uncertainty associated with the logistic regression predictions of ground water vulnerability. Central to the proposed method is the assumption that prediction uncertainty in ground water vulnerability models is a function of input error propagation from uncertainty in the estimated logistic regression model coefficients (model error) and the values of explanatory variables represented in the GIS (data error). Input probability distributions that represent both model and data error sources of uncertainty were simultaneously sampled using a Latin hypercube approach with logistic regression calculations of probability of elevated nonpoint source contaminants in ground water. The resulting probability distribution represents the prediction intervals and associated uncertainty of the ground water vulnerability predictions. The method is illustrated through a ground water vulnerability assessment of the High Plains regional aquifer. Results of the LHS simulations reveal significant prediction uncertainties that vary spatially across the regional aquifer. Additionally, the proposed method enables a spatial deconstruction of the prediction uncertainty that can lead to improved prediction of ground water vulnerability.  相似文献   

5.
Application of the BME approach to soil texture mapping   总被引：3，自引：1，他引：3  
In order to derive accurate space/time maps of soil properties, soil scientists need tools that combine the usually scarce hard data sets with the more easily accessible soft data sets. In the field of modern geostatistics, the Bayesian maximum entropy (BME) approach provides new and powerful means for incorporating various forms of physical knowledge (including hard and soft data, soil classification charts, land cover data from satellite pictures, and digital elevation models) into the space/time mapping process. BME produces the complete probability distribution at each estimation point, thus allowing the calculation of elaborate statistics (even when the distribution is not Gaussian). It also offers a more rigorous and systematic method than kriging for integrating uncertain information into space/time mapping. In this work, BME is used to estimate the three textural fractions involved in a texture map. The first case study focuses on the estimation of the clay fraction, whereas the second one considers the three textural fractions (sand, silt and clay) simultaneously. The BME maps obtained are informative (important soil characteristics are identified, natural variations are well reproduced, etc.). Furthermore, in both case studies, the estimates obtained by BME were more accurate than the simple kriging (SK) estimates, thus offering a better picture of soil reality. In the multivariate case, classification error rate analysis in terms of BME performs considerably better than in terms of kriging. Analysis in terms of BME can offer valuable information to be used in sampling design, in optimizing the hard to soft data ratio, etc.  相似文献   

6.
Prestack depth imaging of seismic data in complex areas such as salt structures requires extensive velocity model updating. In many cases, salt boundaries can be difficult to identify due to lack of seismic reflectivity. Traditional amplitude based segmentation methods do not properly tackle this problem, resulting in extensive manual editing. This paper presents a selection of seismic attributes that can reveal texture differences between the salt diapirs and the surrounding geology as opposed to amplitude‐sensitive attributes that are used in case of well defined boundaries. The approach consists of first extracting selected texture attributes, then using these attributes to train a classifier to estimate the probability that each pixel in the data set belongs to one of the following classes: near‐horizontal layering, highly‐dipping areas and the inside of the salt that appears more like a low amplitude area with small variations in texture. To find the border between the inside of the salt and the highly‐dipping surroundings, the posterior probability of the class salt is input to a graph‐cut algorithm that produces a smooth, continuous border. An in‐line seismic section and a timeslice from a 3D North Sea data set were employed to test the proposed approach. Comparisons between the automatically segmented salt contours and the corresponding contours as provided by an experienced interpreter showed a high degree of similarity.  相似文献   

7.
Sand lenses at various spatial scales are recognized to add heterogeneity to glacial sediments. They have high hydraulic conductivities relative to the surrounding till matrix and may affect the advective transport of water and contaminants in clayey till settings. Sand lenses were investigated on till outcrops producing binary images of geological cross‐sections capturing the size, shape and distribution of individual features. Sand lenses occur as elongated, anisotropic geobodies that vary in size and extent. Besides, sand lenses show strong non‐stationary patterns on section images that hamper subsequent simulation. Transition probability (TP) and multiple‐point statistics (MPS) were employed to simulate sand lens heterogeneity. We used one cross‐section to parameterize the spatial correlation and a second, parallel section as a reference: it allowed testing the quality of the simulations as a function of the amount of conditioning data under realistic conditions. The performance of the simulations was evaluated on the faithful reproduction of the specific geological structure caused by sand lenses. Multiple‐point statistics offer a better reproduction of sand lens geometry. However, two‐dimensional training images acquired by outcrop mapping are of limited use to generate three‐dimensional realizations with MPS. One can use a technique that consists in splitting the 3D domain into a set of slices in various directions that are sequentially simulated and reassembled into a 3D block. The identification of flow paths through a network of elongated sand lenses and the impact on the equivalent permeability in tills are essential to perform solute transport modeling in the low‐permeability sediments.  相似文献   

8.
A number of statistical methods are typically used to effectively predict potential landslide distributions. In this study two multivariate statistical analysis methods were used （weights of evidence and logistic regression） to predict the potential distribution of shallow-seated landslides in the Kamikawachi area of Sabae City, Fukui Prefecture, Japan. First, the dependent variable （shallow-seated landslides） was divided into presence and absence, and the independent variables （environmental factors such as slope and altitude） were categorized according to their characteristics. Then, using the weights of evidence （WE） method, the weights of pairs comprising presence （w^＋（i）） or absence （w^-（i））, and the contrast values for each category of independent variable （evidence）, were calculated, Using the method that integrated the weights of evidence method and a logistic regression model, score values were calculated for each category of independent variable. Based on these contrast values, three models were selected to sum the score values of every gird in the study area. According to a receiver operating characteristic curve analysis （ROC）, model 2 yielded the best fit for predicting the potential distribution of shallow-seated landslide hazards, with 89% correctness and a 54.5% hit ratio when the occurrence probability （OP） of landslides was 70%. The model was tested using data from an area close to the study region, and showed 94% correctness and a hit ratio of 45.7% when the OP of landslides was 70%. Finally, the potential distribution of shallow-seated landslides, based on the OP, was mapped using a geographical information system.  相似文献   

9.
Using auxiliary information to improve the prediction accuracy of soil properties in a physically meaningful and technically efficient manner has been widely recognized in pedometrics. In this paper, we explored a novel technique to effectively integrate sampling data and auxiliary environmental information, including continuous and categorical variables, within the framework of the Bayesian maximum entropy (BME) theory. Soil samples and observed auxiliary variables were combined to generate probability distributions of the predicted soil variable at unsampled points. These probability distributions served as soft data of the BME theory at the unsampled locations, and, together with the hard data (sample points) were used in spatial BME prediction. To gain practical insight, the proposed approach was implemented in a real-world case study involving a dataset of soil total nitrogen (TN) contents in the Shayang County of the Hubei Province (China). Five terrain indices, soil types, and soil texture were used as auxiliary variables to generate soft data. Spatial distribution of soil total nitrogen was predicted by BME, regression kriging (RK) with auxiliary variables, and ordinary kriging (OK). The results of the prediction techniques were compared in terms of the Pearson correlation coefficient (r), mean error (ME), and root mean squared error (RMSE). These results showed that the BME predictions were less biased and more accurate than those of the kriging techniques. In sum, the present work extended the BME approach to implement certain kinds of auxiliary information in a rigorous and efficient manner. Our findings showed that the BME prediction technique involving the transformation of variables into soft data can improve prediction accuracy considerably, compared to other techniques currently in use, like RK and OK.  相似文献   

10.
本文采用三维应力-速度有限差分(SV-FD)方法,数值模拟了横向各向同性(TI)地层对称主轴与井轴斜交情况下正交偶极子声源激发的井孔声场.主要解决了与倾斜角有关的三维空间弹性模量矩阵的推导,柱坐标系下应力-速度有限差分方程组的建立,井轴上场点奇异性与内边界处理等几个关键问题,提高了计算精度.在横向各向同性地层对称主轴与井轴平行的情况下,与实轴积分法所得结果进行了对比,验证了本文方法的正确性.计算了不同倾角情况下xx和yy两分量的弯曲波,并用频域加权相似法提取了弯曲波频散曲线,结果显示了横向各向同性介质中不同方向偏振的弯曲波传播的分裂现象,其频散曲线在低频段分裂,随频率增大而逐渐重合.弯曲波低频截止频率处速度与理论公式得到的横波速度基本符合.  相似文献   

11.
Accurate and reliable methods for quantifying grain size are important for river science, management and in various other sedimentological settings. Remote sensing offers methods of quantifying grain size, typically providing; (a) coarse outputs (c. 1 m) at the catchment scale where individual grains are at subpixel level, or; (b) fine resolution outputs (c. 1 mm) at the patch scale. Recently, approaches using unmanned aerial vehicles (UAVs) have started to fill the gap between these scales, providing hyperspatial resolution data (< 10 cm) over reaches a few hundred metres in length, where individual grains are at suprapixel level. This ‘mesoscale’ is critical to habitat assessments. Most existing UAV‐based approaches use two‐dimensional (2D) textural variables to predict grain size. Validation of results is largely absent however, despite significant differences in platform stability and image quality obtained by manned aircraft versus UAVs. Here, we provide the first quantitative assessment of the accuracy and precision of grain size estimates produced from a 2D image texture approach. Furthermore, we present a new method which predicts subaerial gravel size using three‐dimensional (3D) topographic data derived from UAV imagery. Data is collected from a small gravel‐bed river in Cumbria, UK. Results indicate that our new topographic method gives more accurate measures of grain size (mean residual error ‐0.0001 m). Better results for the image texture method may be precluded by our choice of texture measure, the scale of analysis or the effects of image blur resulting from an inadequate camera gimbal. We suggest that at our scale of assessment, grain size is more strongly related to 3D variation in elevation than to the 2D textural patterns expressed within the imagery. With on‐going improvements, our novel method has potential as the first grain size quantification approach where a trade‐off between coverage and resolution is not necessary or inherent. Copyright © 2017 John Wiley & Sons, Ltd.  相似文献   

12.
Practical application of the power-law regression model with an unknown location parameter can be plagued by non-finite least squares parameter estimates. This presents a serious problem in hydrology, since stream flow data is mainly obtained using an estimated stage–discharge power-law rating curve. This study provides a set of sufficient requirements for the data to ensure the existence of finite least squares parameter estimates for a power-law regression with an unknown location parameter. It is shown that in practice, these requirements act as necessary for having a finite least squares solution, in most cases. Furthermore, it is proved that there is a finite probability for the model to produce data having non-finite least squares parameter estimates. The implications of this result are discussed in the context of asymptotic predictions, inference and experimental design. A Bayesian approach to the actual regression problem is recommended.  相似文献   

13.
Understanding the geological uncertainty of hydrostratigraphic models is important for risk assessment in hydrogeology. An important feature of sedimentary deposits is the directional ordering of hydrostratigraphic units (HSU). Geostatistical simulation methods propose efficient algorithm for assessing HSU uncertainty. Among different geostatistical methods to simulate categorical data, Bayesian maximum entropy method (BME) and its simplified version Markov-type categorical prediction (MCP) present interesting features. In particular, the zero-forcing property of BME and MCP can provide a valuable constrain on directional properties. We illustrate the ability of MCP to simulate vertically ordered units. A regional hydrostratigraphic system with 11 HSU and different abundances is used. The transitional deterministic model of this system presents lateral variations and vertical ordering. The set of 66 (11 × 12/2) bivariate probability functions is directly calculated on the deterministic model with fast Fourier transform. Despite the trends present in the deterministic model, MCP is unbiased for the HSU proportions in the non-conditional case. In the conditional cases, MCP proved robust to datasets over-representing some HSU. The inter-realizations variability is shown to closely follow the amount and quality of data provided. Our results with different conditioning datasets show that MCP replicates adequately the directional units arrangement. Thus, MCP appears to be a practical method for generating stochastic models in a 3D hydrostratigraphic context.  相似文献   

14.
Fluid flow behavior in a porous medium is a function of the geometry and topology of its pore space. The construction of a three dimensional pore space model of a porous medium is therefore an important first step in characterizing the medium and predicting its flow properties. A stochastic technique for reconstruction of the 3D pore structure of unstructured random porous media from a 2D thin section training image is presented. The proposed technique relies on successive 2D multiple point statistics simulations coupled to a multi-scale conditioning data extraction procedure. The Single Normal Equation Simulation Algorithm (SNESIM), originally developed as a tool for reproduction of long-range, curvilinear features of geological structures, serves as the simulation engine. Various validating criteria such as marginal distributions of pore and grain, directional variograms, multiple-point connectivity curves, single phase effective permeability and two phase relative permeability calculations are used to analyze the results. The method is tested on a sample of Berea sandstone for which a 3D micro-CT scanning image is available. The results confirm that the equi-probable 3D realizations obtained preserve the typical patterns of the pore space that exist in thin sections, reproduce the long-range connectivities, capture the characteristics of anisotropy in both horizontal and vertical directions and have single and two phase flow characteristics consistent with those of the measured 3D micro-CT image.  相似文献   

15.
阎维彰  郭跃宏 《地震研究》1991,14(3):215-222
本文利用山西地区的地震地质、地质构造、地震活动性、地球物理场等特征资料,对有关特征进行了图象识别,采用综合决策方案以提高识别的稳定性。综合判定结果分析表明,在研究区内的180个对象中,得到54个D类对象,占全部识别总面积的30％,其中向31个对象,历史上已发生过6级以上强震或在强震附近,另有23个对象为新判定的6级以上强震潜在震源区。得到山西强震潜在震源的预测图。  相似文献   

16.
A nonparametric resampling technique for generating daily weather variables at a site is presented. The method samples the original data with replacement while smoothing the empirical conditional distribution function. The technique can be thought of as a smoothed conditional Bootstrap and is equivalent to simulation from a kernel density estimate of the multivariate conditional probability density function. This improves on the classical Bootstrap technique by generating values that have not occurred exactly in the original sample and by alleviating the reproduction of fine spurious details in the data. Precipitation is generated from the nonparametric wet/dry spell model as described in Lall et al. [1995]. A vector of other variables (solar radiation, maximum temperature, minimum temperature, average dew point temperature, and average wind speed) is then simulated by conditioning on the vector of these variables on the preceding day and the precipitation amount on the day of interest. An application of the resampling scheme with 30 years of daily weather data at Salt Lake City, Utah, USA, is provided.  相似文献   

17.
The Earth??s gravity field is related to the topographic potential in medium and higher degrees, which is isostatically compensated. Hence, the topographic-isostatic (TI) data are indispensable for extending an available Earth Gravitational Model (EGM) to higher degrees. Here we use TI harmonic coefficients to construct a Synthetic Earth Gravitational Model (SEGM) to extend the EGMs to higher degrees. To achieve a high-quality SEGM, a global geopotential model (EGM96) is used to describe the low degrees, whereas the medium and high degrees are obtained from the TI or topographic potential. This study differes from others in that it uses a new gravimetric-isostatic model for determining the TI potential. We test different alternatives based on TI or only topographic data to determine the SEGM. Although the topography is isostatically compensated only to about degree 40?C60, our study shows that using a compensation model improves the SEGM in comparison with using only topographic data for higher degree harmonics. This is because the TI data better adjust the applied Butterworth filter, which bridges the known EGM and the new high-degree potential field than the topographic data alone.  相似文献   

18.
Probabilistic thresholds for triggering shallow landslides by rainfall are developed using two approaches: a logistic regression model and Iverson's physically based model. Both approaches are applied to a 180 km2 area in northern Italy. For the physically based model a Monte Carlo approach is used to obtain probabilities of slope failure associated with differing combinations of rainfall intensity and duration as well as differing topographic settings. For the logistic regression model hourly and daily rainfall data and split‐sample testing are used to explore the effect of antecedent rainfall on triggering thresholds. It is demonstrated that both the statistical and physically based models provide stochastic thresholds that express the probability of landslide triggering. The resulting thresholds are comparable, even though the two approaches are conceptually different. The physically based model also provides an estimate of the percentage of potentially unstable areas in which failure can be triggered with a certain probability. The return period of rainfall responsible for landslide triggering is studied by using a Gumbel scaling model of rainfall intensity–duration–frequency curves. It is demonstrated that antecedent rainfall must be taken into account in landslide forecasting, and a method is proposed to correct the rainfall return period by filtering the rainfall maxima with a fixed threshold of antecedent rainfall. This correction produces an increase of the return periods, especially for rainstorms of short duration. Copyright © 2009 John Wiley & Sons, Ltd.  相似文献   

19.
任意空间取向TI介质中体波速度特征分析   总被引：2，自引：0，他引：2       下载免费PDF全文
郝重涛  姚陈  王迅 《地震地质》2005,27(2):252-259
基于任意空间取向TI介质(简称ATI)中体波速度和偏振解析解,通过模型数值计算给出ATI介质中体波群速度和相速度的变化特征,说明TI空间取向与测线方位对速度的影响。研究表明,体波群速度图案和相速度图案相对TI对称轴固定,随TI对称轴倾角及其相对测线方位角的变化呈现出一定的对称性和重复性;可以针对ATI地区的地质情况,给出体波群速度和相速度变化图案,为进一步的理论研究提供便捷。此结果也可以直接用于VSP(垂直地震剖面)和井间地震资料的分析研究  相似文献   

20.
目的:筛选并分析影响肺结节良恶性的因素,建立预测模型、验证该模型并与梅奥模型、Brock模型对比。方法:回顾性分析2015年1月至2017年12月深圳市人民医院有病理结果的孤立性肺结节病例319例,其中229例作为建模组（A组）,90例作为验证组（B组）,分析A组病例性别、年龄、直径、吸烟史、毛刺、位于上叶、边界不清楚、分叶征、空泡征、血管集束征、胸膜凹陷征、含磨玻璃密度及钙化,通过单因素分析及Logistic回归分析,筛选出独立影响因子,并建立回归方程。用B组资料进行验证并将B组资料分别代入本研究模型、梅奥模型及Brock模型进行对比。结果:单因素分析示年龄、直径、毛刺、上叶、边界不清楚、分叶、空泡、血管集束征、胸膜凹陷征、是否含有磨玻璃密度在良恶性结节中的差异具有统计学意义（P<0.05）,Logistic回归分析示有毛刺、有分叶、边界不清楚和含有磨玻璃密度为恶性孤立性肺结节的独立影响因素,并据此建立的回归方程ROC曲线下面积为0.894,其灵敏度为91.3％,特异度为77.3％,阳性似然比为4.02,阴性似然比为0.11,阳性预测值为80.8％,阴性预测值为89.5％;本研究模型与梅奥模型的差异有统计学意义（P=0.0049）,与Brock模型差异没有统计学意义（P=0.79）。结论:有毛刺、有分叶、边界不清楚和含有磨玻璃密度为恶性孤立性肺结节的独立影响因素,据此建立的回归方程具有较高的诊断效能。本研究建立的模型诊断效能优于梅奥模型,与Brock模型诊断效能相当。   相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号