4.1. Vertical Settlement Characteristics
The loading stages for various cross-sections were divided into four to eight levels to reliably consider the difference in the ground surface elevation and the soft soil settlement. For example, in section N5 + 850, the riprap loading process was divided into seven stages, as shown in
Figure 6. The vertical settlements at the central axis of section N5 + 850 were obtained by CJP-3 and F-2 (
Figure 4) monitoring data and plotted versus the results obtained numerically, as illustrated in
Figure 6.
Figure 6 illustrates that the numerical simulation results agreed well with the monitoring data. Consequently, it can be concluded that the deployed model can accurately simulate the vertical settlement of soft soil ground treated by PVD-assisted staged riprap loading. The settlement of the surface and shallow soil layers had a similar stepped pattern following the applied staged loading. However, the gradual settlement pattern faded deeper in the soil profiles. Furthermore, according to the layered settlement monitoring data from three in situ observation sections, by accumulating the layered settlements in the PVD-treated zone (elevation from −1.8 m to −21.8 m) and the non-PVD zone (elevation from −21.8 m to −30.8 m), the vertical settlement in the PVD zone (
) and the non-PVD zone (
) were determined. It was found that
reached about 85% of the total settlement (
), while
was limited to about 15% of
, as summarized in
Table 5. The comparison suggested that the consolidation settlement occurred mainly in the PVD-treated zone.
4.2. Horizontal Displacement Characteristics
The horizontal displacements at different elevations with time in CX-2 (
Figure 4) of section N5 + 850, obtained from the monitoring data using an inclinometer and numerically, are plotted in
Figure 7. It can be observed that the simulation results were in good agreement with the monitored data. Furthermore, the horizontal displacement mainly occurred in the PVD-treated zone, with a maximum displacement of 246 mm located at around −10 m elevation after 715 days of loading. This −10 m depth was equalt to about half of the used PVD length, while the horizontal displacement decreased sharply for the deeper soil layers and converged to zero at around −40 m elevation.
The filling height of the crest of the Lingni Seawall was higher than the berms, as shown in
Figure 2. Before 130 d, the staged riprap loading was carried out simultaneously at the center and both sides (berms) of the embankment. When the berms reached the intended elevation, the riprap filling was stopped for the berms. After that, only the center part of the embankment continued to be riprap-loaded. Therefore, with the passage of construction time, the horizontal displacement curves in
Figure 7 showed different variation trends due to different loading stages and time periods.
For −10 m (corresponding to the maximum displacement depth), the horizontal displacement increased with time and converged to a constant value after around 600 days. While for −26 m (representing the non-PVD zone), the horizontal displacement converged to a constant value only after 350 days, which mainly corresponded to the initial preloading stage, as shown in
Figure 8.
4.3. Excess Pore Water Pressure Dissipation Regulations
The excess pore water pressure (Δ
u) at different elevations with time at K-2 of section N5 + 850, monitored using pore water pressure piezometers and obtained numerically, is plotted in
Figure 9. It was found that the pore water pressure simulation results were consistent with the monitoring data. For the shallow soil ground in the PVD-treated zone, the pore water pressure increased quickly during the initial preloading stage and converged to zero once the loading became constant. Meanwhile, for the non-PVD zone, the excess pore water pressure increased slowly until reaching a constant value when the loading became constant for the deep soil layers and did not dissipate during the observation period. This phenomenon demonstrated the efficiency of utilizing PVDs to facilitate drainage for deep soft soils characterized by low permeability.
The deployed regulations of the excess pore water pressure dissipation utilizing the PVDs in the Lingni Seawall were investigated. Initially, the original consolidation state of the soft soil ground was analyzed under the hydrostatic water level and compared to the initial pore water pressure observations, as shown in
Figure 10a. It can be observed that the initial pore water pressure almost coincided with the hydrostatic water pressure, which verified that the excess pore water pressure of the natural soil ground was almost zero. Then, by accumulating the increment in excess pore water pressure (Δ
u) at each staged loading (Δ
p), the relation curves between the cumulative excess pore water pressure (
u) and the load (
P =
p) could be drawn, as shown in
Figure 10b. It can be seen that
u in the muck layer increased almost linearly with the load increase. The regulations for the excess pore water pressure dissipation induced by loading with PVD were analogous to previously reported field observations and model testing reports [
33,
34]. Further analysis showed a remarkable variation between the pore pressure coefficient (Δ
u/Δ
p) and the cumulative load (
p), as shown in
Figure 10c.
In
Figure 10b,
u/P shows a three-period variation that was found to be similar to the Δ
u/Δ
p variations illustrated in
Figure 10c. In period 1 and period 3, the excess pore water pressure increment was less than the external load increment. During period 2, the excess pore water pressure increment was greater than or equal to the external load increment, as presented in
Table 6. This was associated with the loading increment and rate during the staged loading. Among all the staged loadings, the fourth stage had the largest increment, resulting in a sharp increase in the excess pore water pressure within the muck layer, as demonstrated in
Figure 9. Finally, it can be concluded that the staged loading rate and the riprap height can be controlled to ensure the safety and stability of the construction.
4.4. Influence of the PVD on the Degree of Consolidation
In this study, the average radial degree of consolidation (
) was calculated using the measured excess pore water pressure (
) [
35]. However,
corresponds to a particular point and at a specific time (
t) and might not represent the average excess pore water pressure (
). Therefore, to calculate
, the coefficient
was introduced for conversion, as shown in Equation (
5) below:
where
is the initiative excess pore water pressure;
is the average radial excess pore water pressure at time
t;
is the measured excess pore water pressure at a particular point in time
t; and
=
/
.
Then, coefficient
was obtained based on Barron’s analytical solution (1948) [
36] under a constant strain condition following Equations (6) and (7):
where
is the radial excess pore water pressure at time
t at a radial distance
r;
t is the consolidation time;
r is the radial distance from the central axis of the vertical drain well;
is the radius of the effective zone of drainage, as shown in
Figure 2;
is the radius of the vertical drain well; and
n is the drain spacing ratio, expressed as
n =
/
.
Furthermore, by transforming Equation (
6), we can obtain Equation (
8). In the Lingni Seawall monitoring plan, the pore water pressure cells were embedded in the center of the effective zone of the PVDs; so, when
r =
,
=
, coefficient
can be calculated using Equation (
9) as follows:
where the PVD equivalent radius
= 0.033 m adopted from Hansbo’s solution (1979) [
37]. For PVDs in a quincunx equilateral triangle pattern with interval spacing (
l) of 1.5 m, the radius of the effective zone of the PVDs [
38]
= 0.788 m; thus,
n = 24. Finally, the coefficient
was found to equal 0.908.
Using Equation (
5),
at the different elevation (
H) was calculated for section N5 + 850 and is plotted in
Figure 11. It can be observed that
obviously decreased with increase in the muck depth. The consolidation effect of the muck layer in the PVD zone was much better than that for the non-PVD zone. In particular, above −11.8 m elevation, corresponding to the half-length of the PVD,
reached up to 75–100%. The results showed that the PVD was efficient for accelerating the consolidation drainage for the soft soil layers; however, the effective depth of the PVD may be different due to the complexity of the prevailing engineering conditions and of the construction process. For example, in section N5 + 850, at the range extending from −11.8 m to −21.8 m elevation,
was limited to about 20–50%. However, in section N3 + 850 at −18 m elevation,
reached up to 80%. The result for the consolidation effectiveness agreed, in general, with the findings of previous research [
39].
4.5. Consolidation Coefficient under PVD Performance
The PVD-assisted staged riprap preloading method involved draining the excess pore water from the surrounding soil layer during the consolidation associated with the hydraulic gradients induced by the embankment preloading. Thus, the water filling the pores flowed easily toward the PVD in the horizontal direction, and then travelled freely along the PVD vertically toward the permeable drainage layer on the ground surface. The consolidation coefficient of the soil layer will change during the loading and draining process, and the actual value of the consolidation coefficient will be different from that measured in the laboratory before construction. Therefore, back-analysis was used to optimize the design and the model parameters based on the monitoring data recorded during the embankment construction [
40]. According to Carrillo’s analytical solution (1942) [
41] of the vertical drain model, the consolidation coefficient of the soil layers can be inversely calculated by the three-point method expressed in Equations (10) and (11) [
42]. The three-point method is illustrated in
Figure 12.
where
is the coefficient of the drain conditions;
,
and
are the settlements at time
,
and
, respectively, and
−
=
−
;
and
are annotated in Equations (6) and (7);
and
are the vertical and horizontal consolidation coefficients, respectively; and
is the vertical drain distance.
Furthermore, by analyzing the experimental data from the oedometer tests using the square-root time method, the relationship between the vertical consolidation coefficient
and the horizontal consolidation coefficient
was determined with
= 1.18
from our previous research [
43]. Thus, the
and
values can be obtained using Equation (
11) for each section, as summarized in
Table 7. Comparing the back-calculation (bc) values of
(bc) and
(bc) to the values determined in the laboratory
(test) and
(test) [
43], it was found that the inversely calculated value was larger than the test value. A good linear relationship was obtained by plotting the relation between the measured values versus the inversely calculated values, as shown in
Figure 13.
Based on the relation illustrated in
Figure 13, it was found that the drainage effect of the PVD and the riprap loading in the soft soil layers was accelerated, with the consolidation coefficient increasing by 1.4 times. The consolidation coefficient is a key parameter for consolidation calculation. It reflects the consolidation rate of soils. After PVD treatment, the consolidation rate of the soft soil ground was also obviously increased with increase in the consolidation coefficient in comparison to the non-PVD zone. Similar conclusions were also drawn in previous studies using back-analyses of the consolidation coefficient of PVD-improved soft clay from settlement data [
44].
4.6. Strength Improvement under PVD Performance
A series of in situ vane shear tests were conducted on the natural muck soil before loading and at six different holes after staged loading to investigate the strength development of the reinforced ground. Holes VST-1, VST-2, and VST-3 correspond to after the berm loading stage, while holes VST-4, VST-5, and VST-6 (
Figure 4) correspond to the full riprap loading conditions. The in situ vane shear strength results from section N5 + 850 are depicted in
Figure 14. The strength development in the PVD zone was found to be greater than that in the non-PVD zone. Moreover, the strength improvement at the central axis was greater than that at the berms. Especially, at shallower than −11.8 m elevations, corresponding to the half-length of the PVD, the soil strength increased significantly up to 200–700%. The strength increased by 60–150% in the range of −11.8 m to −21.8 m elevation in the PVD zone. This pattern was analogous to the verified consolidation and settlement results, reflecting the strength improvement and ground stability. It must be noted that similar strength variation trends in vertical drain-improved clay deposits subjected to vacuum or surcharge loading have been reported in the literature [
45,
46].
4.7. Error and Evaluation Metrics Analysis of Numerical Simulation
To further evaluate the accuracy of the numerical simulation model, the mean absolute error (MAE) and the coefficient of determination or R-squared (R
) were calculated, as presented in
Table 8. All the simulation data had smaller MAE (1.870–7.461) and higher R
(0.962–0.998), which indicated that the model used and the assigned boundaries accurately estimated the vertical settlement, horizontal displacement, and excess pore water pressure. Compared with the totally solid element model, with the R
from 0.914 to 0.977 [
30], the model accuracy and calculating efficiency were increased, and the numbers of elements and nodes were reduced by using the PVD linear element.