Total and Effective Stresses in Backfilled Stopes during the Fill Placement on a Pervious Base for Barricade Design

Backfill is increasingly used in underground mines worldwide. Its successful application depends on the stability of the barricades built at the base of the stopes to hold the backfill in place, which in turn depends on the knowledge of the pore water pressure (PWP) and stresses during, or shortly after, the placement of the slurried backfill. Until now, self-weight consolidation is usually considered for the estimation of the PWP. There is no solution available to evaluate the total and effective stresses during, and shortly after, the filling operation. As excess PWP can simultaneously be generated (increased) and dissipated (decreased) during the backfilling operation, effective stresses can develop when the filling rate is low and/or hydraulic conductivity of the backfill is high. The arching effect has to be considered to evaluate the effective and total stresses in the backfilled stopes. In this paper, a pseudo-analytical solution is proposed to evaluate the effective and total stresses in backfilled stopes during the backfill deposition on a permeable base, by considering the self-weight consolidation and arching effect. The proposed solution is validated by numerical results obtained by Plaxis2D. A few sample applications of the proposed solution are shown.

When the mining backfill is transported by boreholes and pipelines, it must contain sufficient water to ensure workability and facilitate transportation.The construction of a confining structure, called barricade (bulkhead), in the draw-point at the base of the stope is necessary to hold the slurried backfill in place.The stability of the barricade can thus become a critical concern as the barricade failures, previously reported in the literature, have mostly resulted in serious consequences, such as damage of equipment, personal injury, and even loss of lives [21][22][23][24][25][26][27][28][29][30].
The stability analysis of barricades requires the proper estimation of the pore water pressure (PWP), and the total stresses in the backfilled stopes and on the barricades.It is important to know that the most critical moments for the stability of the barricades are during, and shortly after, the placement of the slurried backfill.After these critical moments, the total stresses in the backfilled stopes and on the barricades can dramatically diminish due to the development of arching effect and cementation of cemented backfill.
During the placement of the slurried backfill, excess PWP can be instantaneously generated (increased) with the addition of new backfill layers, and progressively dissipated (decreased).This process is known as self-weight consolidation [31][32][33][34].With the drainage and consolidation of slurried backfill, effective stresses and arching effect can develop in the backfilled stopes, especially when the filling rate is low and the hydraulic conductivity of the backfill is high.
For cemented backfill, the cementation can be an important process because it consumes pore water and reduces the PWP [35][36][37].The PWP and stresses can be further decreased by the generation of hydration products such as calcium silicate hydrate and calcium silicate, which, in turn, increase the strength and stiffness of the cemented backfill [15,18,19,[38][39][40].When the binder content is high (>7%) and the stope filling is slow, the effect of the cementation on the PWP and stresses can be very significant.The neglect of the cement hydration in pressure and stress estimation can lead to an overly conservative design of the barricades.
In most of the underground mines of Canada, the commonly used binding agents include, Portland cement, fly ash, and furnace slag, while the binder content is typically in the range of 3% to 7% by weight [18,38,41].As the cement content is low, the cementation during the stope filling can be insignificant.This has been confirmed by several in situ measurements [42][43][44][45][46].After a certain period of curing time (7 to 28 days), the effect of the cementation can become non-negligible with considerable gains in strength and stiffness [15,18,38,41].As the stope filling is usually completed within a few days, it is conservative to neglect the effect of cementation for the design of barricades.The knowledge of the total stresses is necessary for conducting the stability analysis for the barricades.
Until now, the self-weight consolidation is largely considered for the estimation of the PWP during the placement of the slurried backfill, without any consideration of the effective stress development [31,[46][47][48][49][50][51][52].The arching theory has been applied to evaluate the stresses, either for dry backfill [26,[53][54][55][56][57][58][59][60][61][62], or for submerged backfill with hydrostatic pressure [63][64][65].These conditions correspond to the final steady states of the slurried backfill when the drainage and self-weight consolidation are fully completed, with a pervious (for dry backfill), or an impervious (for submerged backfill) barricade.The conditions do not correspond to the most critical moment for the stability of the barricades during, and after, the filling operation.The question is how to evaluate the total and effective stresses during, and shortly after, the filling operation (the most critical moments for the stability of barricades).To reply to this question, the self-weight consolidation and arching effect should be jointly considered.
In this paper, a pseudo-analytical solution is proposed to calculate the effective and total stresses in backfilled stopes during the deposition of slurried backfill on a permeable base.The self-weight consolidation and arching effect are jointly taken into account.The PWP is first evaluated by using an analytical solution proposed by Zheng et al. [51].The effective and total stresses are then obtained by introducing the PWP expression in the arching model of Li and Aubertin [63].The proposed solution is validated by numerical modeling results obtained by Plaxis2D (v2016, Plaxis B.V. Delft, The Netherlands).

Proposed Solution
Figure 1 represents a typical vertical mine stope, which is continuously filled with a slurried backfill at a filling rate of m (m/h).The stope has a width of B (m) and filled to a final height of H (m). At a given time t, the backfill reaches a height h (=mt).In the figure, x (m) is the elevation of a studied point (0 ≤ x ≤ h).Other parameters shown in the figure will be presented on their application.As the base of the backfill is pervious with zero PWP (Figure 1), all the pore water in the backfill tends to move down and flow out of the stope.When the water reaches its final steady state, the water table falls down to the level of the base, resulting in zero PWP everywhere in the backfill (suction neglected here).Since the excess PWP is defined as the difference between the current total PWP and the final steady-state PWP [51,66,67], the excess PWP in a backfilled stope with a pervious base and barricade is always equal to the current total PWP.Therefore, a distinction between the excess PWP and current total PWP is not necessary.In the following sections, they are noted as (excess) PWP [51].

Solution to Estimate the (Excess) PWP
Gibson [31] proposed a solution to evaluate the (excess) PWP in slurried materials placed on a pervious base with the thickness increasing at a constant speed: where pw (kPa) is the (excess) PWP, t (h) is the filling time, γ (kN/m 3 ) is the saturated unit weight of the backfill, cv (m 2 /g) is consolidation coefficient of the backfill.
Using the boundary conditions pw = 0 at x = 0 and x = h, Gibson [31] gives the following solution: where m (m/h) is the filling rate; E is expressed as follows: where ξ (m) is an arbitrary integration parameter (0 < ξ < ∞).
The gradient of the (excess) PWP can be calculated as: As the base of the backfill is pervious with zero PWP (Figure 1), all the pore water in the backfill tends to move down and flow out of the stope.When the water reaches its final steady state, the water table falls down to the level of the base, resulting in zero PWP everywhere in the backfill (suction neglected here).Since the excess PWP is defined as the difference between the current total PWP and the final steady-state PWP [51,66,67], the excess PWP in a backfilled stope with a pervious base and barricade is always equal to the current total PWP.Therefore, a distinction between the excess PWP and current total PWP is not necessary.In the following sections, they are noted as (excess) PWP [51].

Solution to Estimate the (Excess) PWP
Gibson [31] proposed a solution to evaluate the (excess) PWP in slurried materials placed on a pervious base with the thickness increasing at a constant speed: where p w (kPa) is the (excess) PWP, t (h) is the filling time, γ (kN/m 3 ) is the saturated unit weight of the backfill, c v (m 2 /g) is consolidation coefficient of the backfill.Using the boundary conditions p w = 0 at x = 0 and x = h, Gibson [31] gives the following solution: where m (m/h) is the filling rate; E is expressed as follows: where ξ (m) is an arbitrary integration parameter (0 < ξ < ∞).
The gradient of the (excess) PWP can be calculated as: where F is given as Assuming: and substituting Equation (6) into Equation ( 3), E can be expressed as follows: where As g(y) is an even function, G(y) is also an even function and can be rewritten as follows: Goodwin [68] proposed the following transformation to approximate the calculation of the integral: where h 0 is step length of y; n is a series number in the range of −∞ to +∞.Goodwin [68] has shown that good approximations with high accuracy can be obtained with Equation (10) when the value of h 0 is in the range of 0 to 1 and n is large enough.
Equation (3) can thus be re-written as follows: Substituting Equation (11) into Equation (2) leads to: Equation ( 12) is the analytical solution to evaluate the (excess) PWP in the slurried material during the backfill deposition on a pervious base, which was proposed by Zheng et al. [51] based on the model of Gibson [31].Stable and reliable results of (excess) PWP can be obtained with Equation ( 12) when h 0 is taken as 0.5 and n is in the range of −55 and 55 for most cases [51].
Using the same derivation process as function E, function F can be expressed as follows: Substituting Equations (11) and (13) into Equation (4) yields: Minerals 2019, 9, 38 5 of 19 Equation ( 14) will be introduced into the below section to evaluate the effective and total stresses in the backfilled stopes by jointly taking into account the self-weight consolidation and arching effect.

Solution to Evaluate the Effective and Total Stresses
In Figure 1, the forces acting on the isolated layer element include its own weight W (kN), two compressive forces C (kN), and two lateral shearing forces S (kN) on two sides wall, a vertical total force V (kN) on the top and a vertical total force V + dV (kN) on the base.Considering the equilibrium state of the isolated layer element leads to: The layer element weight W can be expressed as: where l (=h − x) is the depth of the calculation point from the top of the slurried backfill, dl (m) is the thickness of the element layer.
Assuming the vertical total stress σ v is uniformly distributed across the stope width, the vertical force V can be expressed as follows: By applying the Mohr-Coulomb criterion and assuming the horizontal effective stress (σ' h , kPa) is proportional to the vertical effective stress (σ' v , kPa), the lateral shear force S can be calculated as follows: where δ' ( • ) is the effective friction angle along the backfill-wall contact interfaces, which is assumed here to be equal to the value of the effective friction angle of the backfill φ' ( • ); K is the lateral earth pressure coefficient for the placed backfill.Based on numerical and experimental results, some researchers suggested to use Rankine's active earth pressure coefficient K a [54,55,63,64,69]: while others proposed to take Jaky's at-rest earth pressure coefficient K 0 by arguing that the rock walls do not move [56,57,70,71]: Sobhi et al. [72] and Yang et al. [73] revealed that the wall immobilization is a necessary and sufficient condition for the backfill to be in an at-rest condition only when the natural soil (not backfill as shown in most textbook) before excavation (not backfilling) is initially in an at-rest state.For cases of the placement of a backfill in a pre-existing confining structure (rather than construction of a retaining wall later to hold a pre-existing soil), an active state is possible for the newly deposited backfill.Yang et al. [74] further showed why the state of the backfill along the vertical center line of the backfilled stope can be in an at-rest or an active state, depending on the relationship between the value of friction angle and Poisson's ratio.In the case of an at-rest state, the earth pressure coefficient will be expressed as K 0 = µ/(1 − µ), rather than the Jaky's expression.Further discussion on this aspect is beyond the scope of this paper.In the following sections, Rankine's active earth pressure coefficient of K a will be used.
The vertical (σ v , kPa) and horizontal (σ h , kPa) total stresses can be expressed as follows: Introducing Equation (21) to Equation ( 17) leads to: Substituting Equations ( 16), ( 18) and ( 23) into Equation ( 15) yields: Solving Equation ( 25) results in the following expression: where A is an arbitrary constant.By considering the boundary condition σ' v = 0 at l = 0, the vertical effective stress σ' v in the backfill can then be expressed as follows (see the detailed development in Appendix A): Considering that l = h − x (dp w /dl = −dp w /dx), Equation ( 27) can be re-written as: Knowing the vertical effective stress with Equation ( 28), the horizontal effective stress σ' h can then be calculated as: With the calculated PWP and effective stresses, the horizontal (σ h ) and vertical (σ v ) total stresses can be calculated as: Equations ( 28)-( 31) constitute the proposed pseudo-analytical solution to evaluate the vertical and horizontal effective and total stresses during the stope filling operation.For the barricade design, the key parameter is the horizontal total stress.A MATLAB program (R2014a, The MathWorks, Inc., Natick, MA, USA) to solve Equation (31) for the horizontal total stress has been presented in Appendix B by considering a case with It is important to use a large enough value of n max to ensure stable and accurate results of σ h .Sensitivity analysis should be done to study the influence of n max on the calculated σ h .
Figure 2 represents the variation of the calculated σ h at the base of the stope (x = 0) as a function of n max .It can be seen that a stable result of σ h can be obtained when the value of n max exceeds 101 for the given case (the results σ h = 167.37kPa for n max = 1001 and σ h = 167.27kPa for n max = 10,001 are not presented in the figure).This sensitivity analysis should be conducted for each new calculation to ensure stable and reliable results.
accurate results of σh.Sensitivity analysis should be done to study the influence of nmax on the calculated σh.
Figure 2 represents the variation of the calculated σh at the base of the stope (x = 0) as a function of nmax.It can be seen that a stable result of σh can be obtained when the value of nmax exceeds 101 for the given case (the results σh = 167.37kPa for nmax = 1001 and σh = 167.27kPa for nmax = 10,001 are not presented in the figure).This sensitivity analysis should be conducted for each new calculation to ensure stable and reliable results.

Validation of the Analytical Results by Numerical Results
In this section, the proposed pseudo-analytical solution will be validated by numerical results obtained by Plaxis2D [75].
The first validation is to compare the distribution and evolution of the vertical and horizontal total stresses in a backfilled stope calculated by the proposed solution (Equations ( 30) and ( 31)) with those reported by Fahey et al. Figure 3 shows the distribution of the vertical and horizontal total stresses along a vertical profile of 2 m from the vertical center line of the stope, calculated with the proposed solutions (Equations ( 30) and ( 31)) with K = Ka and obtained by the numerical modeling performed by Fahey et al. [70] with Plaxis.It can be seen that the analytical results compared fairly well with the numerical results.The slight difference between the two solutions can be attributed to the relatively large thickness of each backfill layer used in their numerical model.The 5 m for each layer in the numerical models is not thin enough to represent the continuous backfilling in the stope [51].

Validation of the Analytical Results by Numerical Results
In this section, the proposed pseudo-analytical solution will be validated by numerical results obtained by Plaxis2D [75].
The first validation is to compare the distribution and evolution of the vertical and horizontal total stresses in a backfilled stope calculated by the proposed solution (Equations ( 30) and ( 31)) with those reported by Fahey et al.  Figure 3 shows the distribution of the vertical and horizontal total stresses along a vertical profile of 2 m from the vertical center line of the stope, calculated with the proposed solutions (Equations ( 30) and ( 31)) with K = K a and obtained by the numerical modeling performed by Fahey et al. [70] with Plaxis.It can be seen that the analytical results compared fairly well with the numerical results.The slight difference between the two solutions can be attributed to the relatively large thickness of each backfill layer used in their numerical model.The 5 m for each layer in the numerical models is not thin enough to represent the continuous backfilling in the stope [51].30) and ( 31)) with K = K a and obtained by numerical simulation with Plaxis (numerical results taken from Fahey et al. [70]).
In order to further validate the proposed solutions, additional numerical simulations were conducted by using more representative backfill parameters.A series of very thin backfill layers was used in the numerical simulations to well-represent the continuous filling.
Figure 4a represents the physical model of an underground mine stope, which is 4 m wide and 12.4 m high.The stope is continuously filled with a saturated backfill at a filling rate of m = 0.5 m/h to a final height of 12 m with a gap of 0.4 m being left between the stope roof and top surface of the backfill.In practice, a gap that is 3 to 5 m high (the height of a drift) can be left to allow the passage of vehicles or workers.The backfill and rock walls are assumed to be elasto-plastic and obey the Mohr-Coulomb criterion.
The backfill is considered as cohesion-less and has a saturated unit weight of γ = 20 kN/m 3 , unit weight of water γ w = 10 kN/m 3 , drained Young's Modulus E' = 864 kPa, Poisson's ratio µ = 0.2, drained constrained modulus M' = 960 kPa, effective friction angle φ' = 10 • and effective dilation angle ψ' = 0 • .The consolidation coefficient (c v ) can be calculated with the permeability (k) and drained constrained modulus (M ' ) using the following equation: By considering k = 4.176 × 10 −2 m/h (i.e., 1.16 × 10 −5 m/s), one obtains c v = 4 m 2 /h.The surrounding rock mass is assumed to be impermeable.Its geotechnical properties are unit weight γ r = 27 kN/m 3 , Poisson's ratio µ r = 0.25, Young's modulus E r = 42 GPa, cohesion c r = 9.4 MPa, friction angle φ r = 38 • , and dilation angle ψ r = 0 • .Figure 4b presents the numerical model built by Plaxis2D, which shows only half the full numerical model by considering the vertical symmetry plane (VSP).The VSP and the right outer boundary are restrained in the horizontal direction, but can move freely in the vertical direction.The upper outer boundary of the model is free to move in the horizontal and vertical directions, while the lower outer boundary is restrained in all directions.The VSP and fill-wall interface are impermeable, while the top and bottom surfaces of the backfill, after adding a new layer, are permeable and imposed with zero PWP.
The numerical modeling was conducted in three steps.The first step was to build the numerical model with rock mass and attain an initial stress state in it without any excavations.The stope was then excavated in the rock mass.The last step was to fill the stope in layers.In practice, the backfill was continuously filled in the underground mine stope with a gradual increase in thickness.In the numerical model with Plaxis2D, this continuous and gradual increase of the backfill can only be approximately represented by instantaneously adding a series of very thin backfill layers.If the backfill layer thickness is dh, the waiting delay before adding a new layer will be dh/dm.Drainage and consolidation took place during this waiting time.In order to well represent the continuous backfilling in a mine stope, the thickness of each layer dh should be as thin as possible.Meanwhile, the calculation time increases exponentially as the layer thickness decrease and the number of layer increases.A compromise was made between the accuracy and time calculation.In this paper, each layer had a thickness of 0.4 m.The waiting delay before adding a new layer was 0.8 h for a filling rate of 0.5 m/h.The sensitivity analysis showed that the thickness was indeed small enough to ensure numerically stable results were achieved.
Figure 5 represents the distribution of the (excess) PWP (Figure 5a), total (Figure 5b), and effective (Figure 5c) stresses along the vertical center line (VCL) of the stope at the end of backfill deposition, obtained by numerical simulation with Plaxis2D and calculated with the proposed solution (Equations ( 12), ( 28)-( 31)) using K a .The (excess) PWP and total stresses based on the geostatic overburden solution (p w = σ v = σ h = γh) are also plotted on the figure.One sees that a good agreement is obtained between the analytical and numerical results.
Additional numerical simulations were also done.Good agreements were obtained again between the analytical and numerical results.These results are presented in the Ph.D. thesis of Zheng [76].The proposed solution is thus considered to be validated by numerical modeling results.It can thus be applied to calculate the stresses and PWP in backfilled stopes during the accreting deposition.

Sample Application
This section uses the proposed solution to investigate the influence of backfill properties, filling rate, and stope geometry on the total and effective stresses in a backfilled stope at the end of backfill deposition.The slurried backfill, having a saturated unit weight of γ = 20 kN/m 3 , was filled to a final height of H = 40 m.
Figure 6 represents the distribution of the total (Figure 6a) and effective (Figure 6b) stresses along the full height of the backfilled stope at the end of backfill deposition, predicted by the proposed

Sample Application
This section uses the proposed solution to investigate the influence of backfill properties, filling rate, and stope geometry on the total and effective stresses in a backfilled stope at the end of backfill deposition.The slurried backfill, having a saturated unit weight of γ = 20 kN/m 3 , was filled to a final height of H = 40 m.
Figure 6 represents the distribution of the total (Figure 6a) and effective (Figure 6b) stresses along the full height of the backfilled stope at the end of backfill deposition, predicted by the proposed solution when c v increases from 0.1 to 10 m 2 /h.The stope width B is taken as 6 m.The filling rate m is 0.2 m/h.The effective friction angle φ' is taken as 20 • .The Rankine's active earth pressure coefficient K a is 0.49.Without any surprise, the results indicate that the value of c v can significantly influence the effective and total stresses at the end of filling.With a very low value of c v (=0.1 m 2 /h), the drainage and consolidation can be very slow, resulting in total stresses almost equal to the geostatic overburden pressure at the upper part of the backfill.When the consolidation coefficient is high (c v = 10 m 2 /h), significant drainage and consolidation take place during the backfill accreting deposition, resulting in relatively low total stresses and relatively high effective stresses.
Similar results have been shown in previous publications, mostly through numerical modeling [34,70].This is now shown for the first time through pseudo-analytical solutions without numerical modeling.Similar results have been shown in previous publications, mostly through numerical modeling [34,70].This is now shown for the first time through pseudo-analytical solutions without numerical modeling.Figure 7 represents the distribution of the total (Figure 7a) and effective (Figure 7b) stresses along the full height of the backfilled stope when the deposition ceases as the filling rate m varies from 0.1 to 1 m/h, predicted by the proposed solution by considering a consolidation coefficient cv = 5 m 2 /h, an effective friction angle ϕ' = 20°, and an earth pressure coefficient of K = Ka = 0.49.One can see that high total stresses and low effective stresses result from high filling rate.Once again, these results are very similar to those shown in previous publications through numerical modeling [34].However, this is for the first time the results are shown through pseudo-analytical solution.
Near the bottom of the backfill, abnormally high effective stress at high filling rate is observed.The same phenomenon has been observed in previous cases with low consolidation coefficient (Figure 6b).These results are surprising, but easy to understand.When the filling rate is high or when the coefficient consolidation is low, the drainage and consolidation are slow and the arching effect is weak, resulting in high (excess) PWP, high total stresses, and low effective stresses along the full height of the backfill, except near the bottom.As the effective stresses are equal or very close to the total stresses near the pervious bottom along which the PWP is equal or close to zero, the high total stresses near the bottom due to the high total stresses in the upper portion result in high effective stresses near the bottom of the backfilled stope.This explains well the high total stresses along the full height of the stope and high total and effective stresses near the bottom of the stope, with low Figure 7 represents the distribution of the total (Figure 7a) and effective (Figure 7b) stresses along the full height of the backfilled stope when the deposition ceases as the filling rate m varies from 0.1 to 1 m/h, predicted by the proposed solution by considering a consolidation coefficient c v = 5 m 2 /h, an effective friction angle φ' = 20 • , and an earth pressure coefficient of K = K a = 0.49.One can see that high total stresses and low effective stresses result from high filling rate.Once again, these results are very similar to those shown in previous publications through numerical modeling [34].However, this is for the first time the results are shown through pseudo-analytical solution.
Near the bottom of the backfill, abnormally high effective stress at high filling rate is observed.The same phenomenon has been observed in previous cases with low consolidation coefficient (Figure 6b).These results are surprising, but easy to understand.When the filling rate is high or when the coefficient consolidation is low, the drainage and consolidation are slow and the arching effect is weak, resulting in high (excess) PWP, high total stresses, and low effective stresses along the full height of the backfill, except near the bottom.As the effective stresses are equal or very close to the total stresses near the pervious bottom along which the PWP is equal or close to zero, the high total stresses near the bottom due to the high total stresses in the upper portion result in high effective stresses near the bottom of the backfilled stope.This explains well the high total stresses along the full height of the stope and high total and effective stresses near the bottom of the stope, with low consolidation coefficient (Figure 6a) or high filling rate (Figure 7a).8a) and effective (Figure 8b) stresses along the full height of the backfilled stope at the end of filling, with a different effective friction angle of backfill ϕ', predicted by the proposed solution by using a consolidation coefficient of cv = 5 m 2 /h, a filling rate of m = 0.2 m/h, and an earth pressure coefficient of K = Ka.The iso-geostatic overburden pressure is also included in the figure.The results show that the total and effective stresses decrease as the value of the effective friction angle of the backfill increases.These results are straightforward because the high effective friction angle of backfill leads to the generation of a high degree of arching effect, and hence, reduces the effective and total stresses.Figure 9 represents the distribution of the total (Figure 9a) and effective (Figure 9b) stresses along the full height of the backfilled stope when the deposition ceases as stope width varies from 3 to 12 m, predicted by the proposed solution by considering a consolidation coefficient cv = 5 m 2 /h, a filling rate m = 0.2 m/h, an effective friction angle ϕ' = 20°, and an earth pressure coefficient K = Ka = 0.49.The results show that, high total and effective stresses are generated with increased stope width.This Figure 8 represents the distribution of the total (Figure 8a) and effective (Figure 8b) stresses along the full height of the backfilled stope at the end of filling, with a different effective friction angle of backfill φ', predicted by the proposed solution by using a consolidation coefficient of c v = 5 m 2 /h, a filling rate of m = 0.2 m/h, and an earth pressure coefficient of K = K a .The iso-geostatic overburden pressure is also included in the figure.The results show that the total and effective stresses decrease as the value of the effective friction angle of the backfill increases.These results are straightforward because the high effective friction angle of backfill leads to the generation of a high degree of arching effect, and hence, reduces the effective and total stresses.8a) and effective (Figure 8b) stresses along the full height of the backfilled stope at the end of filling, with a different effective friction angle of backfill ϕ', predicted by the proposed solution by using a consolidation coefficient of cv = 5 m 2 /h, a filling rate of m = 0.2 m/h, and an earth pressure coefficient of K = Ka.The iso-geostatic overburden pressure is also included in the figure.The results show that the total and effective stresses decrease as the value of the effective friction angle of the backfill increases.These results are straightforward because the high effective friction angle of backfill leads to the generation of a high degree of arching effect, and hence, reduces the effective and total stresses.Figure 9 represents the distribution of the total (Figure 9a) and effective (Figure 9b) stresses along the full height of the backfilled stope when the deposition ceases as stope width varies from 3 to 12 m, predicted by the proposed solution by considering a consolidation coefficient c v = 5 m 2 /h, a filling rate m = 0.2 m/h, an effective friction angle φ' = 20 • , and an earth pressure coefficient K = K a = 0.49.The results show that, high total and effective stresses are generated with increased stope width.This is because the degree of arching effect decreases as the stope width increases, resulting in high effective and total stresses.

Discussion
This paper proposed a new pseudo-analytical solution to calculate the (excess) PWP, as well as the total and effective stresses in a backfilled stope, by considering the occurrence of self-weight consolidation and arching effect during the placement of the backfill.The proposed solution has been validated by numerical simulation results obtained by Plaxis2D.It thus provides a useful tool to evaluate the pressure and stresses in backfilled stopes.However, one should keep in mind that the proposed solution, presented in this paper, contains some limitations due to the assumptions used in the development process.For instance, the validation of the proposed solution against numerical simulation was realized by considering Rankine's active earth pressure coefficient Ka.With the given problems and material parameters, one can conclude that the application of the proposed solution using Rankine's active earth pressure coefficient Ka is appropriate, while using Jaky's at-rest earth pressure coefficient K0, may lead to an overly conservative barricade design.However, this conclusion should not be generalized in all cases.In reality, the value of the lateral earth pressure coefficient depends on the relationship between the friction angle and Poisson's ratio.When the friction angle is higher than a critical value (determined by the value of Poisson's ratio), or when the Poisson's ratio is higher than a critical value (determined by the friction angle), the value of the lateral earth pressure coefficient is close to the at-rest earth pressure as defined by Poisson's ratio.Otherwise, the value of the lateral earth pressure coefficient should take Rankine's active earth pressure coefficient.More details can be seen in Yang et al. [73,74].
In this study, the (excess) PWP was evaluated by considering, one after the other, the self-weight consolidation and the arching effect.In reality, the arching effect tends to hold the backfill in the original place and has the effect of slowing down the consolidation.When the consolidation is slowed down, the arching development in turn slows down.The slow-down of the arching in turn favors the development of the consolidation.A coupled process is implicated between the consolidation and arching development.An equilibrium state can be found only after a number of iterations between the two processes.More work is required to consider the coupled processes in the future.
Another limitation is related to the application of the Gibson [31] model, which considers a onedimensional consolidation model, under small strain with constant physical, mechanical, and

Discussion
This paper proposed a new pseudo-analytical solution to calculate the (excess) PWP, as well as the total and effective stresses in a backfilled stope, by considering the occurrence of self-weight consolidation and arching effect during the placement of the backfill.The proposed solution has been validated by numerical simulation results obtained by Plaxis2D.It thus provides a useful tool to evaluate the pressure and stresses in backfilled stopes.However, one should keep in mind that the proposed solution, presented in this paper, contains some limitations due to the assumptions used in the development process.For instance, the validation of the proposed solution against numerical simulation was realized by considering Rankine's active earth pressure coefficient K a .With the given problems and material parameters, one can conclude that the application of the proposed solution using Rankine's active earth pressure coefficient K a is appropriate, while using Jaky's at-rest earth pressure coefficient K 0, may lead to an overly conservative barricade design.However, this conclusion should not be generalized in all cases.In reality, the value of the lateral earth pressure coefficient depends on the relationship between the friction angle and Poisson's ratio.When the friction angle is higher than a critical value (determined by the value of Poisson's ratio), or when the Poisson's ratio is higher than a critical value (determined by the friction angle), the value of the lateral earth pressure coefficient is close to the at-rest earth pressure as defined by Poisson's ratio.Otherwise, the value of the lateral earth pressure coefficient should take Rankine's active earth pressure coefficient.More details can be seen in Yang et al. [73,74].
In this study, the (excess) PWP was evaluated by considering, one after the other, the self-weight consolidation and the arching effect.In reality, the arching effect tends to hold the backfill in the original place and has the effect of slowing down the consolidation.When the consolidation is slowed down, the arching development in turn slows down.The slow-down of the arching in turn favors the development of the consolidation.A coupled process is implicated between the consolidation and arching development.An equilibrium state can be found only after a number of iterations between the two processes.More work is required to consider the coupled processes in the future.
Another limitation is related to the application of the Gibson [31] model, which considers a one-dimensional consolidation model, under small strain with constant physical, mechanical, and hydraulic properties.The one-dimensional model can be representative in tailing dams with large horizontal sections.However, the lateral drainage can take place in tailing dams when the waste rock inclusion is applied, or in mine stopes when the stope rock walls contain conductive fractures and joints.
It should be pointed out that the proposed solution (Equations ( 28)-( 31)) was developed by considering fully saturated backfill during the filling operation.This is only possible when the filling rate is large and/or when the consolidation coefficient is low.Otherwise, unsaturation can take place somewhere in the backfill due to the seepage through the pervious base.More work is necessary to take into account the base seepage and the subsequent unsaturation.
It is important to reiterate that the excess PWP is generally defined as the difference between the current total PWP and the PWP when the water reaches a steady state [66].Based on this general definition, two specific definitions of excess PWP result from two different situations [67].The first one applies when the PWP at the final (steady) state is equal to the PWP at the initial state; both are the hydrostatic pressure.The application of the general definition of the excess PWP leads to the difference between the current PWP and the hydrostatic pressure.This is a well-known definition of excess PWP, which is used in the traditional theory of consolidation in Terzaghi [77].The second specific definition applies when the PWP, at the final steady state, falls to zero.This is the case in this study with a pervious base.The PWP at the final steady state is zero when all water flows out of the stope through the pervious base (suction neglected).The application of the general definition of the excess PWP leads to the excess PWP always equaling the current PWP.No distinction is needed between them.This second specific definition of excess PWP is poorly known and not unanimously accepted.It would, however, be impossible to explain the subsidence or settlement of ground surface, which is associated with the drainage and consolidation due to the water pumping from a water well, if the initial hydrostatic pressure before water pumping in the well is not considered as excess PWP.This is because the dissipation of excess PWP is inseparable with the drainage and consolidation in soil mechanics [78].
Finally, further experimental work is required to validate or calibrate the proposed solutions.This work is ongoing and will be part of future publications.

Conclusions
In this paper, a new pseudo-analytical solution was presented to evaluate the (excess) PWP, as well as the effective and total stresses in backfilled stopes, by jointly considering the self-weight consolidation and arching effect during the placement of slurried backfill on a pervious base.The main conclusions are given as follows: 1.
One is used to define the excess PWP as the difference between the current PWP and the hydrostatic pressure.This definition is, however, only valid when the initial and final steady PWP are the same, all equaling the hydrostatic pressure.In this study, the final steady PWP is zero due to the pervious base (and barricade) with zero PWP.The PWP and the excess PWP in the backfilled stopes are always equal to each other.No distinction is needed between the PWP and excess PWP.

2.
The proposed solution has been validated against numerical results obtained by numerical modeling with Plaxis2D.The proposed solution constitutes a simple and useful tool to calculate the pressures and stresses in backfilled stopes for the design of barricades.It can be particularly useful for parametric sensitivity analysis in the preliminary stage of a project.

3.
The calculated effective and total stresses in the backfilled stopes, during the backfill deposition, are lower than the iso-geostatic overburden pressure, due to the occurrence of the consolidation and arching effects.4.
The sample applications of the proposed solution showed that high total stresses along the full height of backfilled stopes can result from low consolidation coefficient, high filling rate, low effective friction angle, or large stope width.Similar results have been reported in previous publications through numerical modeling.This is, however, shown for the first time through pseudo-analytical solutions without numerical modeling.

Figure 1 .
Figure 1.Schematic diagram of a vertical backfilled stope with a pervious base and continuously filled with slurried backfill.

Figure 1 .
Figure 1.Schematic diagram of a vertical backfilled stope with a pervious base and continuously filled with slurried backfill.
[70].The stope is 50 m high and 20 m wide.The backfill is saturated and characterized by the parameters: γ = 20 kN/m 3 (saturated unit weight), γw = 10 kN/m 3 (unit weight of water), k = 4.17 × 10 −3 m/h (permeability), cv = 4.625 m 2 /h (consolidation coefficient), c = 0 kPa (effective cohesion); E' = 10 MPa (drained Young's Modulus), μ = 0.2 (Poisson's ratio), ϕ' = 45° (effective friction angle) and ψ' = 0° (effective dilation angle).In the numerical model, a 2 m layer of backfill was first placed in the stope, followed by a 3 m layer of backfill in 0.1 days.The rest of the stope was filled with nine layers of backfill; each of them has 5 m and was filled in 0.1 days.The total filling time was 1.1 days, resulting an average filling rate of m = 1.89 m/h.
[70].The stope is 50 m high and 20 m wide.The backfill is saturated and characterized by the parameters: γ = 20 kN/m 3 (saturated unit weight), γ w = 10 kN/m 3 (unit weight of water), k = 4.17 × 10 −3 m/h (permeability), c v = 4.625 m 2 /h (consolidation coefficient), c = 0 kPa (effective cohesion); E' = 10 MPa (drained Young's Modulus), µ = 0.2 (Poisson's ratio), φ' = 45 • (effective friction angle) and ψ' = 0 • (effective dilation angle).In the numerical model, a 2 m layer of backfill was first placed in the stope, followed by a 3 m layer of backfill in 0.1 days.The rest of the stope was filled with nine layers of backfill; each of them has 5 m and was filled in 0.1 days.The total filling time was 1.1 days, resulting an average filling rate of m = 1.89 m/h.

Figure 3 .
Figure 3. Distribution of the vertical and horizontal total stresses along a vertical profile of 2 m from the vertical center line of the stope at the end of the filling operation, calculated by the proposed solutions (Equations (30) and (31)) with K = K a and obtained by numerical simulation with Plaxis (numerical results taken from Fahey et al.[70]).

Figure 4 .
Figure 4. (a) The physical model and (b) numerical model of the backfilled stope built by Plaxis2D after considering the vertical symmetry plane (VSP).

Figure 6 .
Figure 6.Distribution of the horizontal and vertical total (a) and effective (b) stresses with different consolidation coefficient c v along the full height of the stope at the end of backfill accreting deposition, predicted by the proposed solution (Equations (28)-(31)) by using: H = 40 m, B = 6 m, γ w = 10 kN/m 3 , γ = 20 kN/m 3 , m = 0.2 m/h, φ' = 20 • , and K = K a .

Figure 8
Figure 8  represents the distribution of the total (Figure8a) and effective (Figure8b) stresses along the full height of the backfilled stope at the end of filling, with a different effective friction angle of backfill ϕ', predicted by the proposed solution by using a consolidation coefficient of cv = 5 m 2 /h, a filling rate of m = 0.2 m/h, and an earth pressure coefficient of K = Ka.The iso-geostatic overburden pressure is also included in the figure.The results show that the total and effective stresses decrease as the value of the effective friction angle of the backfill increases.These results are straightforward because the high effective friction angle of backfill leads to the generation of a high degree of arching effect, and hence, reduces the effective and total stresses.

Figure 7 .
Figure 7. Distribution of the horizontal and vertical total (a) and effective (b) stresses with different filling rate m along the full height of the backfilled stope when the deposition ceases, predicted by the proposed solution (Equations (28)-(31)) by considering H = 40 m, B = 6 m, c v = 5 m 2 /h, γ = 20 kN/m 3 , γ w = 10 kN/m 3 , φ' = 20 • and K = K a .

Figure 8
Figure 8  represents the distribution of the total (Figure8a) and effective (Figure8b) stresses along the full height of the backfilled stope at the end of filling, with a different effective friction angle of backfill ϕ', predicted by the proposed solution by using a consolidation coefficient of cv = 5 m 2 /h, a filling rate of m = 0.2 m/h, and an earth pressure coefficient of K = Ka.The iso-geostatic overburden pressure is also included in the figure.The results show that the total and effective stresses decrease as the value of the effective friction angle of the backfill increases.These results are straightforward because the high effective friction angle of backfill leads to the generation of a high degree of arching effect, and hence, reduces the effective and total stresses.

Figure 9
Figure 9  represents the distribution of the total (Figure9a) and effective (Figure9b) stresses along the full height of the backfilled stope when the deposition ceases as stope width varies from 3 to 12 m, predicted by the proposed solution by considering a consolidation coefficient cv = 5 m 2 /h, a filling rate m = 0.2 m/h, an effective friction angle ϕ' = 20°, and an earth pressure coefficient K = Ka = 0.49.The results show that, high total and effective stresses are generated with increased stope width.This

Figure 8 .
Figure 8. Distribution of the horizontal and vertical total (a) and effective (b) stresses with different effective friction angle φ' along the full height of the backfilled stope when the deposition ceases, predicted by the proposed solution (Equations (28)-(31)) by using H = 40 m, B = 6 m, c v = 5 m 2 /h, m = 0.2 m/h, γ w = 10 kN/m 3 , γ = 20 kN/m 3 and K = K a .