An efficient variance-based global sensitivity index based on C-DRM and Taylor expansion

Author(s):  
Wenxuan Wang ◽  
Xiaoyi Wang

Sensitivity analysis plays an important role in quantifying the impact of input uncertainty on model response uncertainty. Through sensitivity analysis, we can grasp the crucial parameters and gain a deeper understanding of the model behavior. In this paper, an effective analytical solution for solving the variance-based global sensitivity index is proposed. Firstly, the original performance function is approximated as the sum of a series of univariate functions using the conventional dimensional reduction method (C-DRM). Then, the Taylor series expansion is used to expand the univariate function as unary linear function and unary quadratic functions. Finally, the analytical solutions of the variance-based global sensitivity index based on unary linear function and unary quadratic functions are derived respectively. The computational cost of the proposed method is completely concentrated on the calculation of the partial derivative information of the performance function with respect to each variable. As long as the partial derivative information is obtained, the variance-based global sensitivity index can be obtained directly by the proposed method without any additional computational cost. For simple explicit performance functions, the derivative information can be directly derived analytically. For complex explicit or implicit performance functions, the derivative information can be estimated by some simple numerical difference methods. Five examples are studied to investigate the accuracy and efficiency of the proposed method.

SPE Journal ◽  
2013 ◽  
Vol 19 (04) ◽  
pp. 621-635 ◽  
Author(s):  
Cheng Dai ◽  
Heng Li ◽  
Dongxiao Zhang

Summary Reservoir simulations involve a large number of formation and fluid parameters, many of which are subject to uncertainties owing to the combination of spatial heterogeneity and insufficient measurements. Accurately quantifying the impact of varying parameters on simulation models can reveal the importance of the parameters, which helps in designing field-characterization strategies and determining parameterization for history matching. Compared with the commonly used local sensitivity analysis (SA), global SA considers the whole variation range of the parameters and can thus provide more-complete information. However, the traditional global sensitivity analysis that is derived from Monte Carlo simulation (MCS) is computationally too demanding for reservoir simulations. In this study, we propose an alternative approach that is both accurate and efficient. In the proposed approach, the model outputs such as pressure and reservoir production quantities are expressed by polynomial chaos expansions (PCEs). The probabilistic collocation method is used to determine the coefficients of the polynomial expansions by solving outputs at different sets of collocation points by means of the original partial-differential equations. Then, a proxy is constructed with such coefficients. Accurate statistical sensitivity indices of the uncertainty parameters can be obtained by running the proxy. We validate the approach with 2D examples by comparing with the MCS-based global SA. It is found that with only a small fraction of the computational cost required by the MCS approach, the new approach gives accurate global sensitivity for each parameter. The proposed approach is also demonstrated on a large-scale 3D black-oil model, for which the MCS-based global SA is found to be computationally infeasible. It is found that the developed approach possesses the following key advantages: It requires a much smaller number of reservoir simulations for accurate global SA; it is nonintrusive and can be implemented with existing codes or simulators; and it can accommodate arbitrary distributions of parameters encountered in realistic geological situations.


Author(s):  
Hyunkyoo Cho ◽  
Ujjwal Shrestha ◽  
Young-Do Choi ◽  
Jungwan Park

Abstract Global sensitivity analysis (GSA) estimates influence of design variables in the entire design domain on performance measures. Hence, using GSA, important design variables could be found for an engineering application with high dimension which require computationally expensive analyses. Then, similar engineering applications could use selected variables to carry out design process with smaller dimension and affordable computational cost. In this study, GSA has been carried out for the performance measures in design of stay vane and casing of reaction hydraulic turbines. Global sensitivity index method is used for GSA because it can fully capture the effect of interaction between the design variables. For efficiency, genetic aggregation surrogate models are constructed using the responses of computational fluid dynamic (CFD) analysis. Global sensitivity indices for the performance measures of stay vane and casing have been evaluated using the surrogate models. It is found that less than three design variables among 12 are effective in the design process of stay vane and casing in reaction hydraulic turbines.


2017 ◽  
Vol 21 (12) ◽  
pp. 6219-6234 ◽  
Author(s):  
Aronne Dell'Oca ◽  
Monica Riva ◽  
Alberto Guadagnini

Abstract. We propose new metrics to assist global sensitivity analysis, GSA, of hydrological and Earth systems. Our approach allows assessing the impact of uncertain parameters on main features of the probability density function, pdf, of a target model output, y. These include the expected value of y, the spread around the mean and the degree of symmetry and tailedness of the pdf of y. Since reliable assessment of higher-order statistical moments can be computationally demanding, we couple our GSA approach with a surrogate model, approximating the full model response at a reduced computational cost. Here, we consider the generalized polynomial chaos expansion (gPCE), other model reduction techniques being fully compatible with our theoretical framework. We demonstrate our approach through three test cases, including an analytical benchmark, a simplified scenario mimicking pumping in a coastal aquifer and a laboratory-scale conservative transport experiment. Our results allow ascertaining which parameters can impact some moments of the model output pdf while being uninfluential to others. We also investigate the error associated with the evaluation of our sensitivity metrics by replacing the original system model through a gPCE. Our results indicate that the construction of a surrogate model with increasing level of accuracy might be required depending on the statistical moment considered in the GSA. The approach is fully compatible with (and can assist the development of) analysis techniques employed in the context of reduction of model complexity, model calibration, design of experiment, uncertainty quantification and risk assessment.


2017 ◽  
Author(s):  
Aronne Dell'Oca ◽  
Monica Riva ◽  
Alberto Guadagnini

Abstract. We propose new metrics to assist global sensitivity analysis, GSA, of hydrological and Earth systems. Our approach allows assessing the impact of uncertain parameters on main features of the probability density function, pdf, of a target model output, y. These include the expected value of y, the spread around the mean and the degree of symmetry and tailedness of the pdf of y. Since reliable assessment of higher order statistical moments can be computationally demanding, we couple our GSA approach with a surrogate model, approximating the full model response at a reduced computational cost. Here, we consider the generalized Polynomial Chaos Expansion (gPCE), other model reduction techniques being fully compatible with our theoretical framework. We demonstrate our approach through three test cases, including an analytical benchmark, a simplified scenario mimicking pumping in a coastal aquifer, and a laboratory-scale conservative transport experiment. Our results allow ascertaining which parameters can impact some moments of the model output pdf while being uninfluential to others. We also investigate the error associated with the evaluation of our sensitivity metrics by replacing the original system model through a gPCE. Our results indicate that the construction of a surrogate model with increasing level of accuracy might be required depending on the statistical moment considered in the GSA. Our approach is fully compatible with (and can assist the development of) analysis techniques employed in the context of reduction of model complexity, model calibration, design of experiment, uncertainty quantification and risk assessment.


Processes ◽  
2019 ◽  
Vol 7 (9) ◽  
pp. 562 ◽  
Author(s):  
Satyajeet Bhonsale ◽  
Carlos André Muñoz López ◽  
Jan Van Impe

Spray drying is a key unit operation used to achieve particulate products of required properties. Despite its widespread use, the product and process design, as well as the process control remain highly empirical and depend on trial and error experiments. Studying the effect of operational parameters experimentally is tedious, time consuming, and expensive. In this paper, we carry out a model-based global sensitivity analysis (GSA) of the process. Such an exercise allows us to quantify the impact of different process parameters, many of which interact with each other, on the product properties and conditions that have an impact on the functionality of the final drug product. Moreover, classical sensitivity analysis using the Sobol-based sensitivity indices was supplemented by a polynomial chaos-based sensitivity analysis, which proved to be an efficient method to reduce the computational cost of the GSA. The results obtained demonstrate the different response dependencies of the studied variables, which helps to identify possible control strategies that can result in major robustness for the spray drying process.


Author(s):  
Souransu Nandi ◽  
Tarunraj Singh

The focus of this paper is on the global sensitivity analysis (GSA) of linear systems with time-invariant model parameter uncertainties and driven by stochastic inputs. The Sobol' indices of the evolving mean and variance estimates of states are used to assess the impact of the time-invariant uncertain model parameters and the statistics of the stochastic input on the uncertainty of the output. Numerical results on two benchmark problems help illustrate that it is conceivable that parameters, which are not so significant in contributing to the uncertainty of the mean, can be extremely significant in contributing to the uncertainty of the variances. The paper uses a polynomial chaos (PC) approach to synthesize a surrogate probabilistic model of the stochastic system after using Lagrange interpolation polynomials (LIPs) as PC bases. The Sobol' indices are then directly evaluated from the PC coefficients. Although this concept is not new, a novel interpretation of stochastic collocation-based PC and intrusive PC is presented where they are shown to represent identical probabilistic models when the system under consideration is linear. This result now permits treating linear models as black boxes to develop intrusive PC surrogates.


2014 ◽  
Vol 2014 ◽  
pp. 1-11 ◽  
Author(s):  
Lei Cheng ◽  
Zhenzhou Lu ◽  
Luyi Li

An extending Borgonovo’s global sensitivity analysis is proposed to measure the influence of fuzzy distribution parameters on fuzzy failure probability by averaging the shift between the membership functions (MFs) of unconditional and conditional failure probability. The presented global sensitivity indices can reasonably reflect the influence of fuzzy-valued distribution parameters on the character of the failure probability, whereas solving the MFs of unconditional and conditional failure probability is time-consuming due to the involved multiple-loop sampling and optimization operators. To overcome the large computational cost, a single-loop simulation (SLS) is introduced to estimate the global sensitivity indices. By establishing a sampling probability density, only a set of samples of input variables are essential to evaluate the MFs of unconditional and conditional failure probability in the presented SLS method. Significance of the global sensitivity indices can be verified and demonstrated through several numerical and engineering examples.


2021 ◽  
Author(s):  
Emilie Rouzies ◽  
Claire Lauvernet ◽  
Bruno Sudret ◽  
Arthur Vidard

Abstract. Pesticide transfers in agricultural catchments are responsible for diffuse but major risks to water quality. Spatialized pesticide transfer models are useful tools to assess the impact of the structure of the landscape on water quality. Before considering using these tools in operational contexts, quantifying their uncertainties is a preliminary necessary step. In this study, we explored how global sensitivity analysis can be applied to the recent PESHMELBA pesticide transfer model to quantify uncertainties on transfer simulations. We set up a virtual catchment based on a real one and we compared different approaches for sensitivity analysis that could handle the specificities of the model: high number of input parameters, limited size of sample due to computational cost and spatialized output. We compared Sobol' indices obtained from Polynomial Chaos Expansion, HSIC dependence measures and feature importance measures obtained from Random Forest surrogate model. Results showed the consistency of the different methods and they highlighted the relevance of Sobol' indices to capture interactions between parameters. Sensitivity indices were first computed for each landscape element (site sensitivity indices). Second, we proposed to aggregate them at the hillslope and the catchment scale in order to get a summary of the model sensitivity and a valuable insight into the model hydrodynamical behaviour. The methodology proposed in this paper may be extended to other modular and distributed hydrological models as there has been a growing interest in these methods in recent years.


Author(s):  
Rasool Koosha ◽  
Fatemeh Shahsavari

Abstract In the building energy performance simulation, the uncertainty analysis (UA) couples to the sensitivity analysis (SA) to handle ever-existing uncertainties; induced by the sources of uncertainty including random occupants behavior and degradation of building materials over time. As a building simulation tool reaches to a high level of complexity, it becomes more challenging for the sensitivity analysis to deliver reliable outputs; thus the accuracy of the SA results substantially depends upon the number of sample sets or the type of analysis performed. This paper describes a variance-based SA tool integrated into a building Resistance-Capacitance (RC) thermal model. Then, for a hypothetical residential building test case, three distinct first-order sensitivity index simulators and three total sensitivity index simulators are implemented and compared in terms of the dependency of results on the sample size, i.e., the demand for the computational cost.


2020 ◽  
Vol 2 (1) ◽  
Author(s):  
Elvis Felipe Elli ◽  
Neil Huth ◽  
Paulo Cesar Sentelhas ◽  
Rafaela Lorenzato Carneiro ◽  
Clayton Alcarde Alvares

Abstract Eucalyptus-breeding efforts have been made to identify clones of superior performance for growth and yield and how they will interact with global climate changes. This study performs a global sensitivity analysis for assessing the impact of genetic traits on Eucalyptus yield across contrasting environments in Brazil under present and future climate scenarios. The APSIM Next Generation Eucalyptus model was used to perform the simulations of stemwood biomass (t ha−1) for 7-year rotations across 23 locations in Brazil. Projections for the period from 2020 to 2049 using three global circulation models under intermediate (RCP4.5) and high (RCP8.5) greenhouse gas emission scenarios were performed. The Morris sensitivity method was used to perform a global sensitivity analysis to identify the influence of plant traits on stemwood biomass. Traits for radiation use efficiency, leaf partitioning, canopy light capture and fine root partitioning were the most important, impacting the Eucalyptus yield substantially in all environments under the present climate. Some of the traits targeted now by breeders for current climate will remain important under future climates. However, breeding should place a greater emphasis on photosynthetic temperature response for Eucalyptus in some regions. Global sensitivity analysis was found to be a powerful tool for identifying suitable Eucalyptus traits for adaptation to climate variability and change. This approach can improve breeding strategies by better understanding the gene × environment interactions for forest productivity.


Sign in / Sign up

Export Citation Format

Share Document