Uncertainty Analysis of Seepage-Induced Consolidation in a Fractured Porous Medium

Lingai Guo,Marwan Fahs,Hussein Hoteit,Rui Gao,⋆and Qian Shao,⋆

1School of Civil Engineering,Wuhan University,Wuhan,430072,China

2Institut Terre et Environnement de Strasbourg,Universit´e de Strasbourg,CNRS,ENGEES,UMR 7063,Strasbourg,France

3Physical Science and Engineering Division,King Abdullah University of Science and Technology,Thuwal,Saudi Arabia

ABSTRACT Numerical modeling of seepage-induced consolidation process usually encounters significant uncertainty in the properties of geotechnical materials.Assessing the effect of uncertain parameters on the performance variability of the seepage consolidation model is of critical importance to the simulation and tests of this process.To this end,the uncertainty and sensitivity analyses are performed on a seepage consolidation model in a fractured porous medium using the Bayesian sparse polynomial chaos expansion(SPCE)method.Five uncertain parameters including Young’s modulus,Poisson’s ratio,and the permeability of the porous matrix,the permeability within the fracture,and Biot’s constant are studied.Bayesian SPCE models for displacement,flow velocity magnitude,and fluid pressure at several reference points are constructed to represent the input-output relationship of the numerical model.Based on these SPCE models,the total and first-order Sobol’indices are computed to quantify the contribution of each uncertain input parameter to the uncertainty of model responses.The results show that at different locations of the porous domain,the uncertain parameters show different effects on the output quantities.At the beginning of the seepage consolidation process,the hydraulic parameters make major contributions to the uncertainty of the model responses.As the process progresses,the effect of hydraulic parameters decreases and is gradually surpassed by the mechanical parameters.This work demonstrates the feasibility to apply Bayesian SPCE approach to the uncertainty and sensitivity analyses of seepage-induced consolidation problems and provides guidelines to the numerical modelling and experimental testing of such problems.

KEYWORDS Fractured porous media;sensitivity analysis;polynomial chaos expansions;poroelasticity

1 Introduction

Seepage in porous media is a common phenomenon in civil engineering that provides both advantages and disadvantages.On one hand,land subsidence resulting from the seepage process is a critical problem that may induce foundation deformation and differential settlement,and can even damage the surface buildings[1–3].On the other hand,soil consolidation induced by seepage can improve the compressive strength and bearing capacity of the building foundation[4–6].During the compaction of a saturated porous medium,the seepage of fluid and the consolidation of solid happen simultaneously until the pore pressure dissipates and the effective stress of porous matrix balances with the external load.In this process,the fluid-solid interaction plays an important role and makes the problem difficult to be analyzed.Moreover,the existence of fractures and weak intercalations,which brings about the change of local flow distribution and deformation in the porous medium,makes the problem even more complex[7,8].

For decades,the seepage-induced consolidation of porous media has been extensively studied[9–13].The one-dimensional consolidation theory was firstly established by Terzaghi[14],which has been widely applied to the modeling and simulation of land subsidence.Then Biot[15]proposed the generalized three-dimensional consolidation theory based on Terzaghi’s work,and extended the theory to anisotropic and nonlinear porous materials[16,17].Biot’s constantαBwas introduced as a new physical constant to capture the influence of fluid pressure on the solid deformation.Other than saturated porous media,a one-dimensional consolidation theory for unsaturated soil was presented by Fredlund et al.[18].The expression ofαBwas given in a different manner by Zienkiewicz[19],and the basic formulation describing both static and dynamic behaviors of porous media was obtained.With the development of numerical techniques and computational methods,numerical modeling has been extensively used to study the seepage consolidation process and to predict the fluid flow and deformation in porous media during this process[20–22].In the simulation,the fluid is usually assumed to flow through the porous skeleton according to Darcy’s law.The solid deformation is modeled using the mechanical equilibrium equation with consideration of poroelasticity,which counts for the effect of fluid flow on solid deformation within a porous medium.Meanwhile,the influence of solid deformation on the fluid flow is taken into account by including the porosity change in governing equations.This fluid-solid interaction during seepage consolidation process can be simulated accurately by solving these governing equations using appropriate numerical methods in the premise that all the model parameters are predetermined by experiments,which is a difficult task in real geotechnical applications as the material properties are typically heterogeneous and anisotropic.Therefore,the numerical simulation of seepage-induced consolidation in porous media exhibits significant uncertainty resulting from the uncertain physical parameters[23,24].

Sensitivity analysis(SA)is a method that measures how the uncertainty of the model solution can be apportioned to the uncertainty of model parameters,and determines the most influential parameters[25].To understand the impact of uncertain parameters on the fluid flow and solid deformation in porous media,SA has been increasingly applied to geotechnical models.Morris’method was used by Jin et al.[24]and Chen et al.[26]to perform SA on the numerical model based on the Boit’s consolidation theory and the nonlinear rheological theory to study the influence of model parameters on land subsidence and ground fissures.Guo et al.[27]developed a three-dimensional hydraulic-thermal model of the Guide Basin,and assessed the uncertainty of the model output that originates from 9 uncertain physical parameters in this model.The Monte Carlo method was used by Li et al.[28]to perform uncertainty and sensitivity analyses so as to analyze the stability of the Bian Jia Gou tailings pond under flooding conditions.Soueid et al.[29]established a stochastic numerical framework for inverting self-potential data based on the Markov chains Monte Carlo method to localize seeps in dams and characterize their permeability and Darcy velocity.Recently,Liu et al.[30]performed a SA on the numerical model of layered periodic foundations based on the Gauss–Lobatto integration,and showed better convergence of this method than the Monte Carlo-based methods.Although some studies have been found in the literature to apply SA to geotechnical models,to date,there are only few SA studies regarding the seepage consolidation process in fractured porous media[23,31].

The most commonly used techniques to perform SA on mathematical models include the Morris’method,Monte Carlo-based methods,sparse polynomial chaos expansion-based methods.Polynomial chaos expansion has been demonstrated to be efficient and robust for the analyses of uncertainty and sensitivity on the computer models in engineering.In this work,we use a new algorithm proposed by Shao et al.[32]to perform SA and assess uncertainty on the seepage consolidation model in fractured porous media.This algorithm constructs sparse polynomial chaos expansion(SPCE),representing the input-output relation based on the Bayesian model averaging.This approach has been proven to be efficient and accurate in several engineering applications[33–36].Once constructing the SPCE of the original numerical model,the Variance-based Sobol’indices,which measure the impact of input uncertain parameters on the model output,are easy to be calculated analytically from SPCE.This approach not only quantitatively assesses the impact of each uncertain input parameter,but also captures its interaction with other parameters on the model output.The same physical model,as shown in literature[20],is adopted to simulate the seepage consolidation in fractured porous media,where the mechanical and hydraulic properties of the porous material are assumed isotropic.The physical parameters characterizing the porous medium such as the permeability,Young’s modulus,Poisson’s ratio,and Biot’s constant are considered as uncertain input parameters.The fluid pressure,velocity,and deformation of the porous matrix are the model responses of interest.Through this study,the critical parameters that have significant influence on the fluid flow and solid deformation are determined and their effects are quantified by Sobol’indices.The results are expected to identify the dominant parameters on the seepage consolidation process in fractured porous media,and to provide a guideline to the modeling and experiments in geotechnical engineering applications.

This paper is organized as follows.Section 2 describes the problem formulation of the seepage consolidation in fractured porous media.Section 3 briefly introduces the SPCE approach and Sobol’indices.Section 4 discusses the SA results for the seepage consolidation model in a fractured porous medium.Finally,conclusions are given in Section 5.

2 Problem Statement and Mathematical Model

The system under consideration is a two-dimensional domain shown in Fig.1,which is a typical example and widely used to study the seepage consolidation[20,37,38].In this domain,a quadrangle enclosure ofW= 10 m andH= 16 m is filled by a saturated heterogeneous porous medium.A fracture ofL= 8 m locates in the center of the domain with an inclination ofθ= 45◦and the fracture thickness is 0.05 m.The domain is assumed to be sealed on the right,left,and bottom boundaries,while the top boundary is allowed to drain freely.A uniformly distributed mechanical load of 10 kN/m2is applied to the top of the domain and kept constant throughout the seepage consolidation process.The bottom boundary of the domain is fixed,that is,no displacement is allowed on this boundary.On the left and right boundaries,the horizontal displacement is prescribed to be 0,and the fracture is assumed traction free.The plane strain condition is considered for the solid mechanics.The properties of the fluid and the porous medium are assumed to be constant during the consolidation process.

Darcy’s law and the continuity equation are used to describe the fluid flow in the fractured porous media:

Figure 1:Schematic diagram of the seepage consolidation model in a fractured porous medium under a mechanical load

3 Uncertainty and Sensitivity Analyses

4 Results and Discussion

4.1 Verification of the Seepage Consolidation Model

To evaluate the influence of uncertain parameters on the seepage-induced consolidation process in a fractured porous medium,a numerical model is developed using COMSOL Multiphysics and verified by comparing the solution against the results in the literature.The problem shown in Fig.1 is modelled using the poroelasticity multiphysics module,which couples the fluid flow and solid mechanics in porous media.The fracture is considered as an extremely narrow zone with weak mechanical properties(e.g.,very small Young’s modulus comparing to the porous matrix)and a large value of permeability.A mesh containing 4608 triangular elements and 2863 nodes is used to discretize the pressure and displacement fields simultaneously.The parameters adopted are listed in Tab.1,which are the same as in Lamb et al.[20]for the purpose of comparison and verification.A time-dependent study is performed,and the simulated displacement and pressure fields from 1 day to 100 days in the fractured porous domain are shown in Fig.2.To verify these results,the vertical displacement and pressure values at the reference points A and B(see Fig.1),respectively,are extracted at different time steps during the seepage consolidation process.The evolution of the vertical displacement at point A and the pressure at point B in 100 days are plotted and compared with the results in Lamb et al.[20],and good agreements are found in Fig.2.Notably,the vertical displacement at point A is negative in our study,while it is positive in Lamb et al.[20],as different coordinate systems are used.Hence,the absolute values of vertical displacement are compared in Fig.2b.

Table 1:Material properties for verification of the seepage consolidation model[20]

Figure 2:Verification of the seepage consolidation model.(a)Displacement field after 100 days(mm)(b)Surface displacement change during 100 days(c)Pressure field after 100 days(Pa)(d)Pressure change during 100 days

Subsequently,this verified numerical model is used to generate sample data for the uncertainty analysis of seepage-induced consolidation in a fractured porous medium.To this end,5 uncertain input parameters are analyzed,including Young’s modulus(E),Poisson’s ratio(v),the permeability of the porous matrix(k1),the permeability within the fracture(k2),and Biot’s constant(αB).These parameters are assumed to be statistically independent with each other and varying within the uncertain ranges as listed in the Tab.2.To evaluate the effect of these uncertain parameters,several quantities of interest are analyzed using the aforementioned numerical model.The flow velocity magnitude,the pressure,and the vertical displacement within the domain are considered to characterize the fluid flow and deformation of the porous media during the seepage consolidation.Noting that in some cases,the seepage consolidation progresses very slowly,and it can be time-consuming to obtain a steady-state solution using the numerical model.Thus,the seepage consolidation in a duration of 100 days is of special interest in this study.

Table 2:The range of values for uncertain input parameters

4.2 Training and Validation of the Bayesian SPCE Model

In order to perform uncertainty and sensitivity analyseson the seepage consolidation model,the Bayesian approach is used to build SPCE models that describe the relationship between the input parameters and the output responses of the numerical model.For this purpose,several output quantities are assessed to capture the characteristics of the model performance,including the vertical displacement(uA)and the flow velocity magnitude(QA)at point A and the pressure(PB)at point B(Fig.1).These quantities are typical for the characterization of the seepage consolidation in the considered porous domain.In this problem,the seepage and consolidation are induced by a mechanical load applied on the top surface,which is free to flow,and thus,the vertical displacement and flow velocity on the top surface(point A)is of interest.The pressure at point A is a constant equal to 0 during the consolidation process,but the pressure at the bottom boundary(point B)changes significantly,which is thereby selected as the model response of interest.

Figure 3:The flow chart of performing uncertainty and sensitivity analyses on the seepage consolidation model using Bayesian SPCE approach

Fig.4 shows the comparison between the Bayesian SPCE-predicted output quantities and the numerically calculated data from both the training and validation datasets.The calculated vertical displacement and the magnitude of flow velocity at point A varies from less than −30 to−5 mm and from 0 to 6.0 × 10−10m/s,respectively,due to the uncertainty of input parameters.Similarly,the calculated pressure at point B exhibits a large uncertainty and increases from 0 to approximately 1.0 × 104Pa with the variation of input parameters.For the displacement,velocity,and pressure,excellent agreements are found between the numerically calculated results and the Bayesian SPCE-predicted results on both training and validation datasets.Specifically,the coefficients of determination(R2),which measure the ability of a statistical model to fit or explain the observed data,on the training datasets of displacement,velocity,and pressure are larger than 0.99.The values of R2on the validation datasets(>0.98)are slightly smaller than those on the training datasets(Fig.4).For the three output metrics,the Bayesian SPCE model performs the best on the displacement of point A,exhibiting R2values of 0.9999 and 0.9998 for training and validation datasets,respectively.The comparative performances of the constructed Bayesian SPCE models on both the training and validation datasets show their good ability to well reproduce the numerical model responses and indicate their sufficient accuracy to calculate reliable sensitivity indices for assessing the effect of each uncertain input parameter.

Figure 4:Comparison of model outputs between the Bayesian SPCE-predicted values and modelcalculated data at 100 days,including the training dataset and the validation dataset:(a)the vertical displacement at point A,(b)the magnitude of flow velocity at point A,and(c)the pressure at point B

4.3 Uncertainty and Sensitivity Analyses of Model Responses

Based on the above constructed Bayesian SPCE models,the total and first-order Sobol’indices of each uncertain input parameter are calculated using Eqs.(11)–(14)and plotted in Fig.5.These Sobol’indices measure how much the uncertainty of each input parameter contributes to the uncertainty in the vertical displacement and flow velocity at point A and the pressure at point B after 100 days’consolidation.It is observed that for each output metric,the influences of the input parameters are different.For instance,the uncertainty of the vertical displacement at point A is mainly attributed to the uncertainty in the Young’s modulus of the porous matrix,while that of the pressure at point B is dominated by the uncertainty in the permeability of the porous matrix.More detailed discussions for each output metric are as follows.

Figure 5:The total and first-order Sobol’indices of the uncertain input parameters on the vertical displacement and flow velocity magnitude at point A and the pressure at point B after 100 days.(a)Vertical displacement(b)Flow velocity magnitude(c)Pressure

(1)Displacement of the porous domain

The variation of the displacement at point A arises mainly from the uncertainty in Young’s modulus of the porous matrix,where the total and first-order Sobol’indices ofEare larger than 0.85(Fig.5a).The Poisson’s ratio and the permeability of the porous matrix show slight influences on the vertical displacement at point A after 100 days,with Sobol’indices smaller than 0.1.The effects of the permeability within the fracture and Biot’s constant are negligible.Notably,the values of total and first-order Sobol’indices for each parameter have little difference,indicating that the interaction effect between different input parameters is not significant.Furthermore,the marginal effect describing the individual influence of each uncertain input parameter on the vertical displacement at point A is presented(Fig.6).Note that the Sobol’indices of the permeability within the fracture(k2)are 0,resulting in no marginal effect ofk2.The absolute values of the vertical displacement at point A decreases monotonously with the increase ofv,E,andαB,and increases with the increase ofk1(Fig.6).The increase of Poisson’s ratio means that the compressibility of the porous medium becomes smaller,leading to smaller displacements under a compressive load.The increase of Young’s modulus means that the porous medium is more resistant to deformation under a certain mechanical load,which yields a smaller absolute value of the vertical displacement.The increase of permeability in the porous matrix means that the fluid is easier to be discharged out of the domain,resulting in a decrease of pore pressure within the porous media,which in turn makes the domain easier to deform.Typically,the permeability only shows effect on the settlement rate,but has no impact on the final deformation of porous media when the seepage consolidation process reaches steady state.Therefore,the influence of permeability observed on displacement in Fig.6 is due to that we consider the vertical displacement after 100 days and the steady state may not be achieved yet.

Figure 6:Marginal effect of each uncertain input parameter on the vertical displacement of point A at 100 days.(a)Effect of Poisson’s ratio(b)Effect of Young’s modulus(c)Effect of matrix permeability(d)Effect of Biot’s constant

Furthermore,the uncertainty of the vertical displacement in the entire domain is analyzed using the Bayesian SPCE approach.To this end,the vertical displacement after 100 days on each node of the domain is considered as an output metric and exported from the COMSOL Multiphysics.For each output metric,a Bayesian SPCE model is trained,eventually resulting in 2863 Bayesian SPCE models,and based on which the Sobol’indices are calculated.Figs.7a–7c show respectively the spatial distribution of the mean,variance,and standard deviation of the vertical displacement resulting from the uncertain input parameters.The mean values of the vertical displacement range between −14 and 0 mm in the porous domain.Due to the fracture,a discontinuity is observed in the displacement-related distributions,and the absolute values of the vertical displacement above the fracture are significantly larger than that below the fracture.Also,the discrepancy of the vertical displacement from the mean value caused by the uncertain input parameters primarily happens in the upper part above the fracture.Figs.7d–7h show the spatial maps of sensitivity of the vertical displacement to each uncertain input parameter,i.e.,Poisson’s ratio,Young’s modulus,the permeability of the porous matrix,the permeability within the fracture,and Biot’s constant,where the total order Sobol’indices involving the individual effect of each uncertain input parameter and its interaction effects with the other parameters are shown.The influence of input parameters on the vertical displacement depends on the location lying in the domain.For instance,Poisson’s ratio shows negligible influence on the vertical displacement above the fracture(total order Sobol’indices of approximately 0),while has significant effect on that below the fracture(Fig.7d).Notably,Young’s modulus makes the major contribution to the uncertainty of the vertical displacement in the zone lying above the fracture,where the mean value and the standard deviation of the displacement are relatively large(Figs.7a–7c and 7e).In the zone lying below the fracture,although the complementary effect of Poisson’s ratio,Young’s modulus,the permeability of the porous matrix,and Biot’s constant on the uncertainty of vertical displacement seems noticeable,the mean,variance,and standard deviation values of the vertical displacement are almost zero.Therefore,the effect of all considered parameters on the vertical displacement below the fracture is negligible.The permeability within the fracture shows little influence on the vertical displacement distribution in the entire porous domain.

Figure 7:(Continued)

Figure 7:Spatial distribution of the(a)mean,(b)variance,and(c)standard deviation of the vertical displacement and spatial distribution of the total Sobol’indices of(d)the Poisson’s ratio v,(e)the Young’s modulus E,(f)the permeability in porous matrix k1,(g)the permeability within fracture k2,and(h)the Biot’s constant αB on the vertical displacement at 100 days

(2)Velocity of fluid flow

Unlike the vertical displacement,the uncertain input parameters show comparable Sobol’indices on the velocity magnitude of fluid flow at point A(Fig.5b).Young’s modulus and Poisson’s ratio are the most important parameters that contribute to the uncertainty in the flow velocity magnitude at point A,exhibiting Sobol’indices values of approximately 0.59 and 0.27,respectively.The permeability of the porous matrix and Biot’s constant also show considerable effects on the variation of the flow velocity.Notably,the total order Sobol’ indices of Poisson’s ratio,Young’s modulus,and permeability of porous matrix are significantly larger than the corresponding first-order Sobol’ indices,indicating that the coupling effects among them are important to the uncertainty of the flow velocity magnitude.Specifically,although the first-order sensitivity index of the permeabilityk1is small,the total Sobol’index shows a non-negligible effect ofk1on the flow velocity,which majorly originates from the coupling effect betweenk1and other parameters.The marginal effects shown in Fig.8 illustrate the univariate influence caused by each uncertain input parameter on the flow velocity magnitude at point A.Poisson’s ratio and Young’s modulus of the porous matrix show negative influences on the flow velocity magnitude(Figs.8a and 8b).As aforementioned,the increase of the Poisson’s ratio and Young’s modulus means smaller compressibility and stronger resistance to deformation of the porous medium,respectively,both yielding a smaller deformation of the domain,which then leads to a decrease of the fluid pressure and the flow velocity magnitude.In contrast tovandE,Biot’s constant shows an almost linearly positive effect on the flow velocity magnitude(Fig.8d).A turn point is observed in the marginal effect ofk1(Fig.8c).When the permeability of porous matrix is small,i.e.,k1<3 × 10−16m2,the velocity magnitude increases with the increase ofk1,while after the turn point(k1>3 × 10−16m2),the velocity magnitude declines gradually with the increase ofk1.

Figure 8:Marginal effect of each uncertain input parameter on the flow velocity magnitude of point A at 100 days.(a)Effect of Poisson’s ratio(b)Effect of Young’s modulus(c)Effect of matrix permeability(d)Effect of Biot’s constant

Furthermore,to analyze uncertainty in the flow velocity at the transient state during the seepage consolidation process,the total order Sobol’indices of the uncertain input parameters for each time step during the 100 days are calculated.The flow velocity magnitude at three reference points are selected as shown in Fig.1,including the central point of top boundary(point A)and the left and right tips of the fracture(points C and D).The evolution of the Sobol’indices for the flow velocity at the reference points A,C,and D is shown in Fig.9.At the beginning of the consolidation process,the permeability of the porous matrixk1makes the dominant contribution to the uncertainty of the flow velocity at the reference points.As the seepage consolidation progresses,the influence of permeabilityk1is gradually surpassed by that of Young’s modulus and Poisson’s ratio.At the reference point A,which is away from the fracture,the flow velocity is not affected by the permeability within the fracturek2.Contrarily,at reference points C and D,the permeability within the fracture shows significant effect on the flow velocity,especially at the beginning of the process.As shown in Figs.9b and 9c,during the first 20 days,the permeability of the porous matrix and the fracture jointly make the major contribution to the uncertainty in flow velocities near the fracture.In general,the hydraulic properties(k1andk2)have dominant effects on the fluid flow at the beginning,and the mechanical properties(Eandv)show more significant effects in the long run.Additionally,the properties within the fracture show local effect on the fluid flow near the fracture,but have little influence on the overall performance of the domain.

Figure 9:The evolution of total order Sobol’indices for flow velocity magnitude at the reference points A,C,and D in the transient regime.(a) for velocity at A(b) for velocity at C(c)for velocity at D

(3)Fluid pressure

The variation of the fluid pressure at point B arises mainly from the uncertainty in the permeability of the porous matrix,which exhibits the total and first-order Sobol’ indices of approximately 0.7(Fig.5c).Young’s modulus and Poisson’s ratio show moderate influence on the fluid pressure at point B,while the effect of the Biot’s constant is slight.The differences between the values of total and first-order Sobol’ indices forv,E,andk1indicate that the interaction effect exists between these input parameters.The marginal effect of Poisson’s ratio,Young’s modulus,and Biot’s constant on the fluid pressure at point B shows a similar tendency as that on fluid velocity at point A(Figs.8 and 10).As aforementioned,the increase of the Poisson’s ratio and Young’s modulus results in a smaller deformation of the porous domain under a certain mechanical load,which in turn leads to a lower fluid pressure at point B.The permeability of the porous matrix shows a positive marginal effect on the fluid pressure within a narrow range whenk1<5 × 10−17m2.As the permeability increases,the pressure at point B first encounters a sharp falling and then followed by a slow declination.

Figure 10:Marginal effect of each uncertain input parameter on the pressure of point B at 100 days.(a)Effect of Poisson’s ratio(b)Effect of Young’s modulus(c)Effect of matrix permeability(d)Effect of Biot’s constant

5 Conclusions

The seepage-induced consolidation process in a fractured porous medium is modeled numerically,and the uncertainty of the model is analyzed using the Bayesian SPCE approach.To perform the uncertainty and sensitivity analyses,Bayesian SPCE models are constructed to reproduce the relationship between the input parameters(i.e.,Young’s modulus,Poisson’s ratio,the permeability of the porous matrix,the permeability within the fracture,and Biot’s constant)and the output metrics(including the vertical displacement,the flow velocity magnitude,and the fluid pressure at reference points).The comparative performances of the constructed SPCE models on the training and validation datasets demonstrate their sufficient accuracy in calculating the Sobol’indices for uncertain parameters.The uncertainties in the model output metrics are quantitatively attributed to the five uncertain input parameters by the Sobol’indices.As a result,the following conclusions are drawn:

(1)For the settlement in the seepage consolidation process,Young’s modulus of the porous matrixEmakes the major contribution to the uncertainty in the settlement near the drainage outlet.While in the region below the fracture,the settlement is extremely small(almost zero),and thus,the effect of these parameters is negligible.

(2)For the flow velocity magnitude,the hydraulic properties(k1andk2)have dominant effects at the beginning of the seepage consolidation process,and the mechanical properties(Eandv)show more significant effects in the long run.The properties within the fracture show only local effect on the fluid flow near the fracture,but have little influence on the overall performance of the domain.

(3)For the fluid pressure,both hydraulic and mechanical properties have important contributions to the uncertainty in model responses,in which the permeability of the porous matrix shows a major effect on the pressure at the bottom boundary.

This work shows feasibility and effectiveness of the Bayesian SPCE approach for performing uncertainty and sensitivity analyses on the seepage consolidation problem.The influences of uncertain input parameters on the model performance are presented and discussed,which provides guidelines to the numerical modelling and experimental testing of such problems.In this work the fracture topology is assumed to be certain.More realistic analysis would include the effect of uncertain fracture topology on model outputs,which would require an efficient numerical model as the automatic mesh updating is needed to change the fracture topology.This would be an interesting topic for further studies.

Acknowledgement:This work was supported by the National Natural Science Foundation of China(11702199).

Funding Statement:National Natural Science Foundation of China(11702199),Qian Shao.National Natural Science Foundation of China(URL:nsfc.gov.cn).

Conflicts of Interest:The authors declare that they have no conflicts of interest to report regarding the present study.