Hydraulic Effects of Shales in Fluvial-Deltaic Deposits: Ground-Penetrating Radar, Outcrop Observations, Geostatistics, and Three-Dimensional Flow Modeling for the Ferron Sandstone, Utah   总被引：1，自引：0，他引：1  

Djuro Novakovic  Christopher D. White  Rucsandra M. Corbeanu  William S. Hammon III  Janok P. Bhattacharya  George A. McMechan 《Mathematical Geology》2002,34(7):857-893

Ground-penetrating radar (GPR) surveys, outcrop measurements, and cores provide a high-resolution 3D geologic model to investigate the hydraulic effects of shales in marine-influenced lower delta-plain distributary channel deposits within the Cretaceous-age Ferron Sandstone at Corbula Gulch in central Utah, USA. Shale statistics are computed from outcrop observations. Although slight anisotropy was observed in mean length and variogram ranges parallel and perpendicular to pale of low , the anisotropy is not statistically significant and the estimated mean length is 5.4 m. Truncated Gaussian simulation was used to create maps of shales that are placed on variably dipping stratigraphic surfaces interpreted from high-resolution 3D GPR surveys, outcrop interpretations, and boreholes. Sandstone permeability is estimated from radar responses calibrated to permeability measurements from core samples. Experimentally designed flow simulations examine the effects of variogram range, shale coverage fraction, and trends in shale coverage on predicted upscaled permeability, breakthrough time, and sweep efficiency. Approximately 1500 flow simulations examine three different geologic models, flow in the 3 coordinate directions, 16 geostatistical parameter combinations, and 10 realizations for each model. ANOVA and response models computed from the flow simulations demonstrate that shales decrease sweep, recovery, and permeability, especially in the vertical direction. The effect on horizontal flow is smaller. Flow predictions for ideal tracer displacements at Corbula Gulch are sensitive to shale-coverage fraction, but are relatively insensitive to twofold variations in variogram range or to vertical trends in shale coverage. Although the hydraulic effects of shale are statistically significant, the changes in flow responses rarely exceed 20%. As a result, it may be reasonable to use simple models when incorporating analogous shales into models of reservoirs or aquifers.  相似文献   

Numerical Method for Conditional Simulation of Levy Random Fields   总被引：2，自引：0，他引：2  

Scott Painter 《Mathematical Geology》1998,30(2):163-179

Stochastic simulations of subsurface heterogeneity require accurate statistical models for spatial fluctuations. Incremental values in subsurface properties were shown previously to be approximated accurately by Levy distributions in the center and in the start of the tails of the distribution. New simulation methods utilizing these observations have been developed. Multivariate Levy distributions are used to model the multipoint joint probability density. Explicit bounds on the simulated variables prevent nonphysical extreme values and introduce a cutoff in the tails of the distribution of increments. Long-range spatial dependence is introduced through off-diagonal terms in the Levy association matrix, which is decomposed to yield a maximum likelihood type estimate at unobserved locations. This procedure reduces to a known interpolation formula developed for Gaussian fractal fields in the situation of two control points. The conditional density is not univariate Levy and is not available in closed form, but can be constructed numerically. Sequential simulation algorithms utilizing the numerically constructed conditional density successfully reproduce the desired statistical properties in simulations.  相似文献   

Bayesian area-to-point kriging using expert knowledge as informative priors     

《International Journal of Applied Earth Observation and Geoinformation》2014

Area-to-point (ATP) kriging is a common geostatistical framework to address the problem of spatial disaggregation or downscaling from block support observations (BSO) to point support (PoS) predictions for continuous variables. This approach requires that the PoS variogram is known. Without PoS observations, the parameters of the PoS variogram cannot be deterministically estimated from BSO, and as a result, the PoS variogram parameters are uncertain. In this research, we used Bayesian ATP conditional simulation to estimate the PoS variogram parameters from expert knowledge and BSO, and quantify uncertainty of the PoS variogram parameters and disaggregation outcomes. We first clarified that the nugget parameter of the PoS variogram cannot be estimated from only BSO. Next, we used statistical expert elicitation techniques to elicit the PoS variogram parameters from expert knowledge. These were used as informative priors in a Bayesian inference of the PoS variogram from BSO and implemented using a Markov chain Monte Carlo algorithm. ATP conditional simulation was done to obtain stochastic simulations at point support. MODIS (Moderate Resolution Imaging Spectroradiometer) atmospheric temperature profile data were used in an illustrative example. The outcomes from the Bayesian ATP inference for the Matérn variogram model parameters confirmed that the posterior distribution of the nugget parameter was effectively the same as its prior distribution; for the other parameters, the uncertainty was substantially decreased when BSO were introduced to the Bayesian ATP estimator. This confirmed that expert knowledge brought new information to infer the nugget effect at PoS while BSO only brought new information to infer the other parameters. Bayesian ATP conditional simulations provided a satisfactory way to quantify parameters and model uncertainty propagation through spatial disaggregation.  相似文献   

2011年长江中下游旱涝急转及汛期暴雨的对流条件研究     

黄明策  沈新勇  刘会鹏  李小凡 《热带气象学报》2020,36(5):590-602

利用ERA-Interim及雨量和土壤水分观测资料,对比诊断了2011年5—6月长江中下游梅汛前极旱期急转为梅汛期洪涝的极端天气事件的对流条件（水汽、不稳定、抬升作用）差异及特征,并研究条件性湿位涡垂直通量（CMF）指数与暴雨之间的定量关系。结果表明:在极旱期,干冷的东北气流控制,西太平洋副热带高压偏东,低层水汽通量弱且以偏北风输送为主,中低层为下沉气流,无低空急流,等θ_se线稀疏,边界层抬升机制缺乏,是干旱加剧的主要因子;在梅汛期,西南气流增强,西太平洋副热带高压西伸,低层气流在长江地区辐合,低层水汽通量增加且转为西南和东南风输送为主,伴随高低空急流耦合和深厚的上升运动,等θ_se线密集形成梅雨锋,增强不稳定暖湿空气强迫抬升和垂直输送,造成暴雨频发,引起区域性洪涝。暴雨中心600 hPa以下为负湿位涡的不稳定层,对流不稳定与条件性对称不稳定共同作用是强降水发生的不稳定机制。CMF指数与旱涝变化、暴雨过程演变非常一致,在极旱（梅汛）期,CMF指数低（高）,变化平缓（剧烈）,CMF指数在暴雨开始时逐步剧增,结束时迅速减小。   相似文献   

A bivariate gamma distribution for use in multivariate flood frequency analysis     

Sheng Yue 《水文研究》2001,15(6):1033-1045

A gamma distribution is one of the most frequently selected distribution types for hydrological frequency analysis. The bivariate gamma distribution with gamma marginals may be useful for analysing multivariate hydrological events. This study investigates the applicability of a bivariate gamma model with five parameters for describing the joint probability behavior of multivariate flood events. The parameters are proposed to be estimated from the marginal distributions by the method of moments. The joint distribution, the conditional distribution, and the associated return periods are derived from marginals. The usefulness of the model is demonstrated by representing the joint probabilistic behaviour between correlated flood peak and flood volume and between correlated flood volume and flood duration in the Madawask River basin in the province of Quebec, Canada. Copyright © 2001 John Wiley & Sons, Ltd.  相似文献   

Geostatistical analysis of field hydraulic conductivity in compacted clay     

A. S. Rogowski  D. E. Simmons 《Mathematical Geology》1988,20(4):423-446

Hydraulic conductivity(K) of fractured or porous materials is associated intimately with water flow and chemical transport. Basic concepts imply uniform flux through a homogeneous cross-sectional area. If flow were to occur only through part of the area, actual rates could be considerably different. Because laboratory values ofK in compacted clays seldom agree with field estimates, questions arise as to what the true values ofK are and how they should be estimated. Hydraulic conductivity values were measured on a 10×25 m elevated bridge-like platform. A constant water level was maintained for 1 yr over a 0.3-m thick layer of compacted clay, and inflow and outflow rates were monitored using 10×25 grids of 0.3-m diameter infiltration rings and outflow drains subtending approximately 1×1 m blocks of compacted clay. Variography of inflow and outflow data established relationships between cores and blocks of clay, respectively. Because distributions of outflow rates were much less and bore little resemblance to the distributions of break-through rates based on tracer studies, presence of macropores and preferential flow through the macropores was suspected. Subsequently, probability kriging was applied to reevaluate distribution of flux rates and possible location of macropores. Sites exceeding a threshold outflow of 100×10^–9 m/s were classified as outliers and were assumed to probably contain a significant population of macropores. Different sampling schemes were examined. Variogram analysis of outflows with and without outliers suggested adequacy of sampling the site at 50 randomly chosen locations. Because of the potential contribution of macropores to pollutant transport and the practical necessity of extrapolating small plot values to larger areas, conditional simulations with and without outliers were carried out. Simulated scenarios based on all available data compared well with conditional simulations based on randomly chosen locations.This paper was presented at MGUS 87 Conference, Redwood City, California, 14 April 1987.  相似文献   

Excel在GPS基线解算精度统计中的应用     

张宏 《测绘与空间地理信息》2012,(6):140-143

GPS基线解算就是利用GPS观测值,通过数据处理,得到测站间的基线向量值。基线解算是GPS网平差的基础和关键,基线解算精度直接关系着GPS网的精度,所以,基线解算精度统计是一项很有必要的工作。Excel中大量的公式函数可以应用选择,主要用来进行有繁重计算任务的预算、财务、数据处理等工作,借助于Excel的数据筛选和条件格式功能,可以方便、高效、准确无误地进行GPS精度统计,这在GPS广泛大量使用的今天,有一定的现实意义。  相似文献   

The Roles of Spatial Locations and Patterns of Initial Errors in the Uncertainties of Tropical Cyclone Forecasts   总被引：1，自引：0，他引：1  

CHEN Boyu  MU Mu 《大气科学进展》2012,29(1):63-78

In this study,a series of sensitivity experiments were performed for two tropical cyclones (TCs),TC Longwang (2005) and TC Sinlaku (2008),to explore the roles of locations and patterns of initial errors in uncertainties of TC forecasts.Specifically,three types of initial errors were generated and three types of sensitive areas were determined using conditional nonlinear optimal perturbation (CNOP),first singular vector (FSV),and composite singular vector (CSV) methods.Additionally,random initial errors in randomly selected areas were considered.Based on these four types of initial errors and areas,we designed and performed 16 experiments to investigate the impacts of locations and patterns of initial errors on the nonlinear developments of the errors,and to determine which type of initial errors and areas has the greatest impact on TC forecasts.Overall,results from the experiments indicate the following:(1) The impact of random errors introduced into the sensitive areas was greater than that of errors themselves fixed in the randomly selected areas.From the perspective of statistical analysis,and by comparison,the impact of random errors introduced into the CNOP target area was greatest.(2) The initial errors with CNOP,CSV,or FSV patterns were likely to grow faster than random errors.(3) The initial errors with CNOP patterns in the CNOP target areas had the greatest impacts on the final verification forecasts.  相似文献   

Algorithm Studies on How to Obtain a Conditional Nonlinear Optimal Perturbation (CNOP)     

SUN Guodong  MU Mu  ZHANG Yale 《大气科学进展》2010,27(6):1311-1321

The conditional nonlinear optimal perturbation (CNOP), which is a nonlinear generalization of the linear singular vector (LSV), is applied in important problems of atmospheric and oceanic sciences, including ENSO predictability, targeted observations, and ensemble forecast. In this study, we investigate the computational cost of obtaining the CNOP by several methods. Differences and similarities, in terms of the computational error and cost in obtaining the CNOP, are compared among the sequential quadratic programming (SQP) algorithm, the limited memory Broyden-Fletcher-Goldfarb-Shanno (L-BFGS) algorithm, and the spectral projected gradients (SPG2) algorithm. A theoretical grassland ecosystem model and the classical Lorenz model are used as examples. Numerical results demonstrate that the computational error is acceptable with all three algorithms. The computational cost to obtain the CNOP is reduced by using the SQP algorithm. The experimental results also reveal that the L-BFGS algorithm is the most effective algorithm among the three optimization algorithms for obtaining the CNOP. The numerical results suggest a new approach and algorithm for obtaining the CNOP for a large-scale optimization problem.  相似文献   

利用星载测雨雷达探测结果对夏季中国南方对流和层云降水气候特征的分析   总被引：10，自引：5，他引：5  

刘鹏  傅云飞 《大气科学》2010,34(4):802-814

本文利用热带测雨卫星(TRMM)上搭载的测雨雷达(PR)十年的探测结果, 对夏季中国南方对流降水和层云降水的气候特征进行了分析。研究结果表明:夏季中国南方层云降水频次较对流降水频次高出两倍以上, 而对流降水强度至少是层云降水强度的4倍; 就整个中国南方而言, 这两种类型的降水对总降水量贡献相当。日变化分析表明夏季中国南方大部分地区的对流降水主要出现在午后, 层云降水出现时间并不集中, 但这两类降水的频次日变化均显示了明显的地域性特征; 对降水廓线日变化的分析结果表明, 对流降水和层云降水廓线的日变化主要表现在“雨顶”高度的日变化, 即对流降水云的厚度有明显的日变化变化特征, 不同地区的降水廓线存在明显的差异。降水率剖面分析结果显示了对流降水的“雨顶” 高度日变化较层云降水剧烈, 降水率的日变化则相反, 且层云降水率的地域性特征更强。  相似文献