scholarly journals Identification of high-permeability subsurface structures with multiple point geostatistics and normal score ensemble Kalman filter

2017 ◽  
Vol 548 ◽  
pp. 208-224 ◽  
Author(s):  
Francesco Zovi ◽  
Matteo Camporese ◽  
Harrie-Jan Hendricks Franssen ◽  
Johan Alexander Huisman ◽  
Paolo Salandin
2011 ◽  
Vol 44 (2) ◽  
pp. 169-185 ◽  
Author(s):  
Haiyan Zhou ◽  
Liangping Li ◽  
Harrie-Jan Hendricks Franssen ◽  
J. Jaime Gómez-Hernández

2011 ◽  
Vol 8 (4) ◽  
pp. 6749-6788 ◽  
Author(s):  
L. Li ◽  
H. Zhou ◽  
H. J. Hendricks Franssen ◽  
J. J. Gómez-Hernández

Abstract. The normal-score ensemble Kalman filter (NS-EnKF) is tested on a synthetic aquifer characterized by the presence of channels with a bimodal distribution of its hydraulic conductivities. Fourteen scenarios are analyzed which differ among them in one or various of the following aspects: the prior random function model, the boundary conditions of the flow problem, the number of piezometers used in the assimilation process, or the use of covariance localization in the implementation of the Kalman filter. The performance of the NS-EnKF is evaluated through the ensemble mean and variance maps, the connectivity patterns of the individual conductivity realizations and the degree of reproduction of the piezometric heads. The results show that (i) the localized NS-EnKF can identify correctly the channels when a large number of conditioning piezometers are used even when an erroneous prior random function model is used, (ii) localization plays an important role to prevent filter inbreeding and results in a better logconductivity characterization, and (iii) the NS-EnKF works equally well under very different flow configurations.


2012 ◽  
Vol 2012 ◽  
pp. 1-18 ◽  
Author(s):  
Haiyan Zhou ◽  
Liangping Li ◽  
J. Jaime Gómez-Hernández

The localized normal-score ensemble Kalman filter is shown to work for the characterization of non-multi-Gaussian distributed hydraulic conductivities by assimilating state observation data. The influence of type of flow regime, number of observation piezometers, and the prior model structure are evaluated in a synthetic aquifer. Steady-state observation data are not sufficient to identify the conductivity channels. Transient-state data are necessary for a good characterization of the hydraulic conductivity curvilinear patterns. Such characterization is very good with a dense network of observation data, and it deteriorates as the number of observation piezometers decreases. It is also remarkable that, even when the prior model structure is wrong, the localized normal-score ensemble Kalman filter can produce acceptable results for a sufficiently dense observation network.


2012 ◽  
Vol 16 (2) ◽  
pp. 573-590 ◽  
Author(s):  
L. Li ◽  
H. Zhou ◽  
H. J. Hendricks Franssen ◽  
J. J. Gómez-Hernández

Abstract. The normal-score ensemble Kalman filter (NS-EnKF) is tested on a synthetic aquifer characterized by the presence of channels with a bimodal distribution of its hydraulic conductivities. This is a clear example of an aquifer that cannot be characterized by a multiGaussian distribution. Fourteen scenarios are analyzed which differ among them in one or various of the following aspects: the prior random function model, the boundary conditions of the flow problem, the number of piezometers used in the assimilation process, or the use of covariance localization in the implementation of the Kalman filter. The performance of the NS-EnKF is evaluated through the ensemble mean and variance maps, the connectivity patterns of the individual conductivity realizations and the degree of reproduction of the piezometric heads. The results show that (i) the localized NS-EnKF can characterize the non-multiGaussian underlying hydraulic distribution even when an erroneous prior random function model is used, (ii) localization plays an important role to prevent filter inbreeding and results in a better logconductivity characterization, and (iii) the NS-EnKF works equally well under very different flow configurations.


2012 ◽  
Vol 132 (10) ◽  
pp. 1617-1625
Author(s):  
Sirichai Pornsarayouth ◽  
Masaki Yamakita

Author(s):  
Nicolas Papadakis ◽  
Etienne Mémin ◽  
Anne Cuzol ◽  
Nicolas Gengembre

2021 ◽  
Vol 14 (6) ◽  
Author(s):  
Jinming Yang ◽  
Chengzhi Li

AbstractSnow depth mirrors regional climate change and is a vital parameter for medium- and long-term numerical climate prediction, numerical simulation of land-surface hydrological process, and water resource assessment. However, the quality of the available snow depth products retrieved from remote sensing is inevitably affected by cloud and mountain shadow, and the spatiotemporal resolution of the snow depth data cannot meet the need of hydrological research and decision-making assistance. Therefore, a method to enhance the accuracy of snow depth data is urgently required. In the present study, three kinds of snow depth data which included the D-InSAR data retrieved from the remote sensing images of Sentinel-1 synthetic aperture radar, the automatically measured data using ultrasonic snow depth detectors, and the manually measured data were assimilated based on ensemble Kalman filter. The assimilated snow depth data were spatiotemporally consecutive and integrated. Under the constraint of the measured data, the accuracy of the assimilated snow depth data was higher and met the need of subsequent research. The development of ultrasonic snow depth detector and the application of D-InSAR technology in snow depth inversion had greatly alleviated the insufficiency of snow depth data in types and quantity. At the same time, the assimilation of multi-source snow depth data by ensemble Kalman filter also provides high-precision data to support remote sensing hydrological research, water resource assessment, and snow disaster prevention and control program.


2021 ◽  
Vol 11 (7) ◽  
pp. 2898
Author(s):  
Humberto C. Godinez ◽  
Esteban Rougier

Simulation of fracture initiation, propagation, and arrest is a problem of interest for many applications in the scientific community. There are a number of numerical methods used for this purpose, and among the most widely accepted is the combined finite-discrete element method (FDEM). To model fracture with FDEM, material behavior is described by specifying a combination of elastic properties, strengths (in the normal and tangential directions), and energy dissipated in failure modes I and II, which are modeled by incorporating a parameterized softening curve defining a post-peak stress-displacement relationship unique to each material. In this work, we implement a data assimilation method to estimate key model parameter values with the objective of improving the calibration processes for FDEM fracture simulations. Specifically, we implement the ensemble Kalman filter assimilation method to the Hybrid Optimization Software Suite (HOSS), a FDEM-based code which was developed for the simulation of fracture and fragmentation behavior. We present a set of assimilation experiments to match the numerical results obtained for a Split Hopkinson Pressure Bar (SHPB) model with experimental observations for granite. We achieved this by calibrating a subset of model parameters. The results show a steady convergence of the assimilated parameter values towards observed time/stress curves from the SHPB observations. In particular, both tensile and shear strengths seem to be converging faster than the other parameters considered.


Sign in / Sign up

Export Citation Format

Share Document