Analysis of Factors Affecting Vacuum Formation and Drainage in the Siphon-Vacuum Drainage Method for Marine Reclamation

: Traditional drainage methods for marine reclamation typically consume large amounts of energy and have a negative environmental impact. The siphon-vacuum drainage method (SVD) automatically forms a vacuum and drains using less energy. It has significant potential for research and application. In this study, a theoretical model is used to calculate the vacuum formation process and drainage rate. Qualitative analysis and global sensitivity analysis were conducted to investigate the effect of various factors in the SVD on vacuum formation and drainage. The qualitative analysis suggests that modifying the length and diameter of the siphon pipe and the thickness of the sealing soil column to increase the siphon rate can improve the vacuum degree and drainage efficiency. Sobol global sensitivity analysis reveals that the sealing soil column thickness is the main factor affecting the vacuum, with a first-order sensitivity index accounting for up to 79.48%. The impact of cylinder diameter and the local resistance coefficient (0.43%) can be almost neglected. A fitting equation for estimating the maximum achievable vacuum is provided. Calculations show that the vacuum formed by the SVD can reach over 80 kPa. This work can help optimize SVD design and advance environmentally friendly marine reclamation projects.


Introduction
Many developed regions are situated in coastal areas, benefiting from access to maritime resources and trade routes.The demand for land is escalating due to population growth and economic expansion in these areas.The marine reclamation method can be used to increase the land area in coastal regions [1].Marine dredged soil, commonly utilized in reclamation, has characteristics including high compressibility, high water content, low bearing capacity, and low permeability [2][3][4].The drainage consolidation method is economical and efficient, making it a prevalent choice for large-scale marine reclamation [5][6][7][8].Techniques such as surcharge preloading, vacuum preloading, and electroosmosis are often employed for drainage [9][10][11].The surcharge preloading applies an external load to the soil surface to promote water expulsion.The vacuum preloading enhances drainage by creating a vacuum under a membrane, reducing pore water pressure.Electroosmosis uses electric currents to drive water and ion movement, facilitating soil dewatering.Among them, vacuum preloading is widely used.It has been implemented in various projects, including the Tianjin Lingang Economic Zone in China [12], the expansion of Fisherman Island in Australia [13], and land reclamation in Singapore [14].
However, vacuum preloading requires a considerable amount of electricity.Traditional vacuum pumps require 75 kW for every 10,000 m 2 of soft soil, and a new water and gas separation vacuum pump demands 25 kW [15].In most cases, vacuum preloading lasts at least 90 days.If the electricity originates from coal, this process will result in significant carbon dioxide emissions.Furthermore, the transmission and maintenance of the vacuum will consume a considerable amount of materials [16].Taking the maintenance of the vacuum as an example, a significant amount of plastic film is required for sealing.The production and disposal of plastic film can lead to greenhouse gas emissions.Additionally, the abandonment of plastic film can lead to soil and water pollution.Today, climate change and environmental pollution have become global concerns [17,18].Consequently, exploring low-carbon and environmentally friendly methods is an effective measure to alleviate these issues.
Siphon is a gravity-driven drainage method that operates without the need for external power [19].However, the siphon drainage at high lift will be interrupted due to air accumulation.Cai et al. [20] found that small-diameter siphon pipes could solve this problem, allowing the siphon to work effectively for a long time.These small-diameter siphon pipes have been successfully applied in slope drainage engineering.Sun et al. [21] applied the siphon drainage method to soft soil drainage first.The pore water within the soft soil seeps into the vertical drainage wells under the action of a siphon.The model tests have shown that siphon drainage could reduce the water level in soft soil, enhance drainage efficiency, lower the pore water pressure, and thereby increase the effective stress of the soft soil.Shen et al. [22] conducted a full-scale field test on improving soft soil using the siphon drainage method.The results indicated a 2.78 m drop in groundwater levels and a settlement of 17.82 cm after 16.5 days of employing siphon.This method eliminates the necessity for vacuum pumps, sealing membranes, and sand cushions.Furthermore, it is easy to construct and maintain in engineering projects.In summary, the siphon drainage method emerges as an economical, low-carbon, and environmentally friendly method for improving marine dredged soil.However, groundwater seepage occurs due to its own weight in siphon drainage.The drainage effect of the siphon drainage method on surface soil is not as effective as vacuum preloading, which relies on vacuum pressure.
In vacuum preloading, the vacuum pressure is used to discharge pore water, thereby inducing soil consolidation [23].In addition, the drainage rate has an impact on both project duration and costs.The increase in vacuum pressure can facilitate drainage.Therefore, the effectiveness and economy of the siphon drainage method for improving soft soils can be enhanced by the addition of vacuum pressure.The siphon-vacuum drainage method (SVD) is shown in Figure 1.This method uses a siphon to automatically create vacuum pressure innovatively in soft soil.Shu et al. [24] have studied the mechanism of vacuum formation and developed a theoretical model to calculate the vacuum pressure and drainage rate.The vacuum pressure and drainage are influenced by factors in the SVD.In order to achieve greater vacuum pressure and better drainage effectiveness, it is necessary to analyze the factors in the SVD.
This research aims to improve understanding of how SVD factors affect vacuum formation and drainage for better marine reclamation.Firstly, the model for calculating the vacuum pressure and drainage rate is summarized.Then, the qualitative effects of each factor on vacuum formation and drainage are analyzed by the calculation results.Moreover, a global sensitivity analysis was conducted to obtain the proportion that each factor affects vacuum formation.Finally, the maximum vacuum degree that can be formed theoretically by the SVD was calculated.The results of this study can serve as a basis for the optimized design of the SVD.This research aims to improve understanding of how SVD factors affect vacuum formation and drainage for better marine reclamation.Firstly, the model for calculating the vacuum pressure and drainage rate is summarized.Then, the qualitative effects of each factor on vacuum formation and drainage are analyzed by the calculation results.Moreover, a global sensitivity analysis was conducted to obtain the proportion that each factor affects vacuum formation.Finally, the maximum vacuum degree that can be formed theoretically by the SVD was calculated.The results of this study can serve as a basis for the optimized design of the SVD.

Cylindrical Model
When dealing with large-scale marine dredged soil, a large number of vertical drainage wells are typically inserted into the dredged soil.The strength of marine dredged soils is improved, as pore water flows into vertical drainage wells.The vertical drainage bodies are arranged in a certain form.As shown in Figure 2, the common forms of arrangement are square and equilateral triangles.Square and hexagonal units are obtained by uniformly dividing the soil layers.The units are simplified to a circle for ease of calculation.The conversion equations between the diameter of the circle (de) and the spacing of the vertical drainage wells (l) are, respectively:  When dealing with large-scale marine dredged soil, a large number of vertical drainage wells are typically inserted into the dredged soil.The strength of marine dredged soils is improved, as pore water flows into vertical drainage wells.The vertical drainage bodies are arranged in a certain form.As shown in Figure 2, the common forms of arrangement are square and equilateral triangles.Square and hexagonal units are obtained by uniformly dividing the soil layers.The units are simplified to a circle for ease of calculation.The conversion equations between the diameter of the circle (d e ) and the spacing of the vertical drainage wells (l) are, respectively: When improving marine dredged soil with the SVD, the siphon pipes and permeable chambers will be inserted into the soil layer.The simplified cylindrical model of a single permeable chamber is shown in Figure 3.The bottom and sides of the cylinder are impermeable.The initial groundwater level is located at the surface of the soil layer.The bottom of the permeable chamber is situated at the bottom of the cylinder, and the permeable chamber is sealed by saturated dredged soil.The thickness of the sealing soil column above the chamber is represented by hs.The center point at the bottom of the permeable chamber serves as the origin (0, 0) for a two-dimensional axisymmetric coordinate system.When improving marine dredged soil with the SVD, the siphon pipes and permeable chambers will be inserted into the soil layer.The simplified cylindrical model of a single permeable chamber is shown in Figure 3.The bottom and sides of the cylinder are imper-meable.The initial groundwater level is located at the surface of the soil layer.The bottom of the permeable chamber is situated at the bottom of the cylinder, and the permeable chamber is sealed by saturated dredged soil.The thickness of the sealing soil column above the chamber is represented by h s .The center point at the bottom of the permeable chamber serves as the origin (0, 0) for a two-dimensional axisymmetric coordinate system.When improving marine dredged soil with the SVD, the siphon pipes and permeable chambers will be inserted into the soil layer.The simplified cylindrical model of a single permeable chamber is shown in Figure 3.The bottom and sides of the cylinder are impermeable.The initial groundwater level is located at the surface of the soil layer.The bottom of the permeable chamber is situated at the bottom of the cylinder, and the permeable chamber is sealed by saturated dredged soil.The thickness of the sealing soil column above the chamber is represented by hs.The center point at the bottom of the permeable chamber serves as the origin (0, 0) for a two-dimensional axisymmetric coordinate system.

Calculation Equations
After the SVD starts, both the air pressure and the water level in the permeable chamber will decrease.The air pressure at point A (pA, t) will become negative.The pore water in the dredged soil will seep into the permeable chamber and then be discharged through the siphon.A theoretical model of vacuum formation has been developed by Shu et al. [24].The model involves detailed equations derivations for the three processes of siphon drainage, seepage drainage, and changes in air pressure.The calculations are based on the principle of energy conservation for siphon drainage The Dupuit assumption is used to describe the distribution of seepage velocities.The non-steady seepage process is broken

Calculation Equations
After the SVD starts, both the air pressure and the water level in the permeable chamber will decrease.The air pressure at point A (p A, t ) will become negative.The pore water in the dredged soil will seep into the permeable chamber and then be discharged through the siphon.A theoretical model of vacuum formation has been developed by Shu et al. [24].The model involves detailed equations derivations for the three processes of siphon drainage, seepage drainage, and changes in air pressure.The calculations are based on the principle of energy conservation for siphon drainage The Dupuit assumption is used to describe the distribution of seepage velocities.The non-steady seepage process is broken down into a sequence of steady flows, which are solved iteratively.It is assumed that the gas within the permeable chamber dissolves and releases in water almost immediately.
In this study, the equations in the model are generalized and simplified.The SVD factors involved in the three processes are summarized.The flow rate of siphon drainage (U out, t ) and the amount of gas dissolved in it (V out, t ) can be expressed by Equations ( 3)-( 5): where a out, t is the acceleration of siphon at time t; h w, t−1 is the water level in the permeable chamber at time t−1; p A, t−1 is the air pressure of point A at time t−1; u out, t−1 is the siphon drainage velocity at time t−1; H 0 is the depth of the permeable chamber; p s is the standard atmospheric pressure (101.3 kPa); ζ i and λ are the local and frictional resistance coefficients of the siphon pipe, respectively; γ w is the weight of water; d and r are the diameter and radius of the siphon pipe, respectively; d pc is the diameter of the permeable chamber; ∆t is the time interval; h c is the distance between the inlet of the siphon pipe and the bottom of the permeable chamber; and V m is the solubility of gas dissolved in water under standard atmospheric pressure.According to Equations (3)-( 5), the factors L, d, d pc , ζ i , and λ in the SVD are related to the siphon drainage.The water level (h ′ w, t ) and the total volume of gas (V ′ all, t ) in the permeable chamber after siphon drainage are as follows: where r pc is the radius of the permeable chamber, and V all, t−1 is the total volume of gas in the chamber at time t−1.
The gas in the chamber will reach a balance between dissolution and release.The relationship between the total volume of gas in the chamber (V ′ all, t ) and the air pressure at point A (p ′ A, t ) can be expressed as follows: where h s is the thickness of the sealing soil column.The first term on the left side of the equation represents the volume of the gaseous gas in the cavity, while the second term represents the volume of the dissolved gas in the water within the chamber.According to Equations ( 6)-( 8), the factors h s and d pc in the SVD are related to the changes in air pressure.Due to the short time of vacuum formation, the seepage at this stage is dominated by the horizontal direction.Therefore, only seepage in the horizontal direction is considered in this model.The flow rate of the seepage (U in, t ) and the amount of gas dissolved in it (V in, t ) can be expressed by Equations ( 9) and (10).
V in, t = 2πK h ∆t where H(r pc , z, t) represents the boundary conditions at the outlet of seepage, and K h is the horizontal hydraulic conductivity.According to Equations ( 9) and ( 10), the factors K h , h s , d e , and d pc in the SVD are related to the seepage drainage.

Results and Discussion of Vacuum Pressure and Drainage
The above equations were compiled using MATLAB.Time was discretized into multiple steps, and the iterative calculation method is based on the results of the previous step to compute the results of the next step.The calculations within each step follow the sequence of siphon drainage, resulting in changes in air pressure and water level, subsequent seepage, and the corresponding variations of air pressure and water level.When the air pressure reaches its minimum value, the process of vacuum formation is considered complete.

Vacuum Pressure
The theoretical model has been validated by full-scale field test results [24], with one such test result depicted in Figure 4, including both measured and calculated results.Table 1 provides the parameters of the field test.In the test, permeable chambers were arranged in a square pattern with a spacing of 0.9 m.The diameter (d e ) of the cylinder was calculated to be 1.015 m according to Equation (1).As shown in Figure 4, the model calculation results are in good agreement with the measured results.At 351 s, the measured results show that the air pressure reaches its minimum value at −42.7 kPa, the model calculation result is −41.5 kPa.The air pressure rapidly decreases in the first 50 s.After 100 s, the rate of air pressure decline significantly slows down.Before 400 s, the air pressure stabilizes essentially.
such test result depicted in Figure 4, including both measured and calculated results.Table 1 provides the parameters of the field test.In the test, permeable chambers were arranged in a square pattern with a spacing of 0.9 m.The diameter (de) of the cylinder was calculated to be 1.015 m according to Equation (1).As shown in Figure 4, the model calculation results are in good agreement with the measured results.At 351 s, the measured results show that the air pressure reaches its minimum value at −42.7 kPa, the model calculation result is −41.5 kPa.The air pressure rapidly decreases in the first 50 s.After 100 s, the rate of air pressure decline significantly slows down.Before 400 s, the air pressure stabilizes essentially.Table 1.The parameters and their values required for the model calculation.

Parameters Values
Thickness of sealing soil column, h s (m) Although the theoretical model disregards the time for gas dissolution and release in water, it still effectively simulates the vacuum formation process.Thus, the model can be used to analyze the impact of various factors on vacuum formation.If more accurate equations could be proposed to represent seepage and siphon drainage, the calculated results of the vacuum formation process would align more closely with the measured results.

Drainage
The water storage within the chamber is influenced by both siphon drainage and seepage drainage.The siphon drainage reduces the water storage within the chamber, while seepage drainage increases it.The flow rates of siphon drainage and seepage drainage are calculated based on the theoretical model described above.The model parameters are as given in Table 1.The model calculation results for the flow rate of siphon drainage and seepage drainage are given in Figure 5.The area of the shaded part is the result of the change in water storage within the chamber.The process of vacuum formation is divided into an early stage and a later stage.In the early stage, the flow rate of the siphon experiences a brief, rapid increase, and it is significantly greater than that of seepage.Hence, siphon drainage plays a dominant role in the variability of water storage within the chamber during this early stage.Subsequently, the flow rate of the siphon decreased while the flow rate of seepage gradually increased.The impact of seepage drainage in the later stage is strengthened.At 400 s, the flow rate of seepage drainage is very similar to that of siphon drainage.The flow change trends of siphon and seepage by calculations are consistent with theoretical analysis.Considering the agreement between the calculated and measured results of air pressure, it is believed that the model can be utilized to calculate the seepage flow rate.Thus, the calculated results can also be used to analyze the effects of various factors on seepage.However, for computational efficiency, the model uses stable well flow theory instead of actual unstable seepage, which makes it unsuitable for calculating seepage over extended periods.
The effectiveness of the SVD in improving marine dredged soil is related to the amount of pore water discharge.The pore water is discharged through seepage.Therefore, the flow rate of seepage drainage was selected to evaluate the effectiveness of the SVD in discharging pore water from marine dredged soil.

Results and Discussion on the Impact of Factors
The discharge of pore water from marine dredged soil using the SVD involves three main processes: siphoning, air pressure changes in the chamber, and seepage.Therefore, the siphon, chamber size, and seepage flow rate are the primary aspects influencing vacuum formation and drainage.Factors affecting the siphon flow rate include the siphon pipe length (L), the minor head loss coefficient of the siphon pipe (ζ), and the diameter of the siphon pipe (d).Factors influencing chamber size include the thickness of the sealing soil column (hs) and the diameter of the permeable chamber (dpc).Additionally, factors affecting the flow rate of seepage include the horizontal permeability coefficient (Kh), the diameter of the permeable chamber (dpc), and the diameter of the cylindrical model (de).
Based on the engineering conditions and the dimensions of existing materials, different values for these factors are selected to calculate vacuum formation and drainage.Furthermore, the effects of these factors on vacuum formation and drainage were analyzed.The factor values for the reference scenario (0-0) are shown in Table 1.Fourteen scenarios have been designed to investigate the influence of seven factors.Table 2 provides the values for these fourteen scenarios.The flow change trends of siphon and seepage by calculations are consistent with theoretical analysis.Considering the agreement between the calculated and measured results of air pressure, it is believed that the model can be utilized to calculate the seepage flow rate.Thus, the calculated results can also be used to analyze the effects of various factors on seepage.However, for computational efficiency, the model uses stable well flow theory instead of actual unstable seepage, which makes it unsuitable for calculating seepage over extended periods.
The effectiveness of the SVD in improving marine dredged soil is related to the amount of pore water discharge.The pore water is discharged through seepage.Therefore, the flow rate of seepage drainage was selected to evaluate the effectiveness of the SVD in discharging pore water from marine dredged soil.

Results and Discussion on the Impact of Factors
The discharge of pore water from marine dredged soil using the SVD involves three main processes: siphoning, air pressure changes in the chamber, and seepage.Therefore, the siphon, chamber size, and seepage flow rate are the primary aspects influencing vacuum formation and drainage.Factors affecting the siphon flow rate include the siphon pipe length (L), the minor head loss coefficient of the siphon pipe (ζ), and the diameter of the siphon pipe (d).Factors influencing chamber size include the thickness of the sealing soil column (h s ) and the diameter of the permeable chamber (d pc ).Additionally, factors affecting the flow rate of seepage include the horizontal permeability coefficient (K h ), the diameter of the permeable chamber (d pc ), and the diameter of the cylindrical model (d e ).
Based on the engineering conditions and the dimensions of existing materials, different values for these factors are selected to calculate vacuum formation and drainage.Furthermore, the effects of these factors on vacuum formation and drainage were analyzed.The factor values for the reference scenario (0-0) are shown in Table 1.Fourteen scenarios have been designed to investigate the influence of seven factors.Table 2 provides the values for these fourteen scenarios.Figure 6 presents the variation curves of vacuum formation and drainage rate (total seepage rate) for different diameters of siphon pipes (d).As shown in Figure 6a, there are significant differences observed in the pressure variation among the three scenarios during the early stage.In the later stage, the pressure change for all three scenarios becomes relatively gradual.The larger the diameter of the siphon pipe (d), the smaller the final air pressure formed (greater vacuum degree).The impact of the diameter (d) on the final vacuum degree is relatively small.A comparison between Figures 6b and 6a reveals that the impact of the diameter of the siphon pipe (d) on drainage is similar to its effect on vacuum formation.The effect of the siphon pipe diameter (d) on drainage and vacuum formation is similar when comparing Figure 6b with Figure 6a.

Length of Siphon Pipe
As the siphon pipe length (L) increases, the resulting vacuum degree and drainage rate decrease, as shown in Figure 7.The siphon pipe length (L) affects the frictional resistance along the siphon pipe.As a result, it also only influences the early stages of vacuum formation and drainage through the siphon rate.The frictional resistance increases due to the increase in siphon pipe length (L).This resistance diminishes the siphon rate, thereby resulting in a lower vacuum degree.It is observed that the siphon pipe length (L) has a relatively minor effect on the final vacuum degree and drainage rate, according to Figure 7.The diameter of the siphon pipe (d) primarily affects the siphon rate and does not directly influence seepage.A larger siphon rate results in a faster decrease in water storage and air pressure within the chamber (cavity enlargement).The seepage rate is increased by the smaller water pressure (volume of water storage) and air pressure in the chamber.Additionally, the siphon primarily impacts the early stage of drainage.Therefore, the diameter of the siphon pipe (d) mainly affects the vacuum formation and drainage in the early stages.Although the vacuum degree and drainage rate can be increased by using a larger diameter of siphon pipe (d), the increase is small.

Length of Siphon Pipe
As the siphon pipe length (L) increases, the resulting vacuum degree and drainage rate decrease, as shown in Figure 7.The siphon pipe length (L) affects the frictional resistance along the siphon pipe.As a result, it also only influences the early stages of vacuum formation and drainage through the siphon rate.The frictional resistance increases due to the increase in siphon pipe length (L).This resistance diminishes the siphon rate, thereby resulting in a lower vacuum degree.It is observed that the siphon pipe length (L) has a relatively minor effect on the final vacuum degree and drainage rate, according to Figure 7.

Length of Siphon Pipe
As the siphon pipe length (L) increases, the resulting vacuum degree and drainage rate decrease, as shown in Figure 7.The siphon pipe length (L) affects the frictional resistance along the siphon pipe.As a result, it also only influences the early stages of vacuum formation and drainage through the siphon rate.The frictional resistance increases due to the increase in siphon pipe length (L).This resistance diminishes the siphon rate, thereby resulting in a lower vacuum degree.It is observed that the siphon pipe length (L) has a relatively minor effect on the final vacuum degree and drainage rate, according to Figure 7.

Local Resistance Coefficient of Siphon Pipe
The local resistance coefficient (ζ) affects vacuum formation and drainage in the same way as the siphon pipe length (L).The siphon rate can be reduced by increasing the local resistance coefficient (ζ). Figure 8 shows that the process of vacuum formation and the increase in drainage rate slows down at a local resistance coefficient of 700.However, the final vacuum degree and drainage rate are almost unaffected.A local resistance coefficient (ζ) of 700 is an unusually large value.This unusual local resistance may occur when the siphon pipe is accidentally bent.However, the local resistance hardly discounts the final vacuum degree and drainage rate formed by the SVD.It shows that the siphon pipe has some resistance to risk.The siphon pipe can effectively perform its intended function when it is neither leaking nor clogging.

Local Resistance Coefficient of Siphon Pipe
The local resistance coefficient (ζ) affects vacuum formation and drainage in the same way as the siphon pipe length (L).The siphon rate can be reduced by increasing the local resistance coefficient (ζ). Figure 8 shows that the process of vacuum formation and the increase in drainage rate slows down at a local resistance coefficient of 700.However, the final vacuum degree and drainage rate are almost unaffected.A local resistance coefficient (ζ) of 700 is an unusually large value.This unusual local resistance may occur when the siphon pipe is accidentally bent.However, the local resistance hardly discounts the final vacuum degree and drainage rate formed by the SVD.It shows that the siphon pipe has some resistance to risk.The siphon pipe can effectively perform its intended function when it is neither leaking nor clogging.

Thickness of Sealing Soil Column
The thickness of the sealing soil column (hs) determines the height of the permeable chamber.The cavity formation and air pressure changes are influenced by the volume of the permeable chamber.The surface area of the chamber determines the cross-sectional area of the seepage.As shown in Figure 9a, the thickness (hs) has a significant impact on

Thickness of Sealing Soil Column
The thickness of the sealing soil column (h s ) determines the height of the permeable chamber.The cavity formation and air pressure changes are influenced by the volume of the permeable chamber.The surface area of the chamber determines the cross-sectional area of the seepage.As shown in Figure 9a, the thickness (h s ) has a significant impact on the early stage of vacuum formation and the final vacuum degree.A lower thickness (h s ) results in more water storage in the chamber.More time is required to discharge more water storage in order to achieve a balance between siphon and seepage rates.Hence, the process of vacuum formation will be more gradual and prolonged in scenarios where the thickness (h s ) is smaller.
the permeable chamber.The surface area of the chamber determines the cross-sectional area of the seepage.As shown in Figure 9a, the thickness (hs) has a significant impact on the early stage of vacuum formation and the final vacuum degree.A lower thickness (hs) results in more water storage in the chamber.More time is required to discharge more water storage in order to achieve a balance between siphon and seepage rates.Hence, the process of vacuum formation will be more gradual and prolonged in scenarios where the thickness (hs) is smaller.
The effect of the thickness of the sealing soil column (hs) on the drainage rate is similar to its effect on the vacuum formation (Figure 9a,b).The thickness (hs) affects both the vacuum formation and the cross-sectional area.In scenario 4-1, a greater vacuum degree results in a greater seepage velocity.Additionally, the cross-sectional area in scenario 4-1 is larger, resulting in a greater seepage rate.Consequently, the difference in seepage rates between the three scenarios gradually becomes significant in the later stages.
The effect of the sealing soil column thickness (hs) on both vacuum formation and drainage rate is significant.Decreasing the thickness of the sealing soil column (hs) is beneficial for increasing both the final vacuum degree and the drainage rate.Therefore, it is an important factor that requires particular attention in engineering applications.The effect of the thickness of the sealing soil column (h s ) on the drainage rate is similar to its effect on the vacuum formation (Figure 9a,b).The thickness (h s ) affects both the vacuum formation and the cross-sectional area.In scenario 4-1, a greater vacuum degree results in a greater seepage velocity.Additionally, the cross-sectional area in scenario 4-1 is larger, resulting in a greater seepage rate.Consequently, the difference in seepage rates between the three scenarios gradually becomes significant in the later stages.
The effect of the sealing soil column thickness (h s ) on both vacuum formation and drainage rate is significant.Decreasing the thickness of the sealing soil column (h s ) is beneficial for increasing both the final vacuum degree and the drainage rate.Therefore, it is an important factor that requires particular attention in engineering applications.

Diameter of Permeable Chamber
The vacuum formation curves at different diameters of the permeable chambers (d pc ) are presented in Figure 10a.The diameter of the chambers (d pc ) affects both the permeable chamber volume and the cross-sectional area of the seepage, which is similar to the thickness of the sealing soil column (h s ).However, the effects of the two on vacuum formation are not the same.The diameter of the chamber (d pc ) changes the time required to achieve a stable vacuum degree but has very little effect on the final vacuum degree.This is because the increase in the diameter of the permeable chamber has a greater promoting effect on seepage, leading to the seepage flow rate approaching the siphon flow rate more quickly.
(Figure 10b).The increase in the diameter (dpc) slows down the vacuum formation, resulting in a slow growth of the seepage rate in the early stages.However, the vacuum degrees of the three scenarios are essentially similar in the later stage.A larger diameter of the chamber (dpc) results in a greater cross-sectional area, enabling a higher seepage rate in the later stage.Therefore, increasing the inner diameter of the permeable chamber (dpc) appropriately can expedite drainage and reduce construction time.

Horizontal Hydraulic Conductivity
The seepage rate is influenced directly by the horizontal hydraulic conductivity (Kh).Therefore, the horizontal hydraulic conductivity (Kh) primarily influences the later stages of vacuum formation and drainage.As shown in Figure 11, the influence of the horizontal hydraulic conductivity (Kh) on the seepage rate is significant, while its effect on the vacuum degree is minor.It shows that the vacuum formation in low-permeability soil is mainly affected by siphon drainage and permeable chamber size.The hydraulic conductivity is an inherent property of the soil and can be challenging to adjust.However, a higher hydraulic conductivity contributes to drainage.There is an intersection in the drainage rate curves between scenarios 5-1 and 5-2 (Figure 10b).The increase in the diameter (d pc ) slows down the vacuum formation, resulting in a slow growth of the seepage rate in the early stages.However, the vacuum degrees of the three scenarios are essentially similar in the later stage.A larger diameter of the chamber (d pc ) results in a greater cross-sectional area, enabling a higher seepage rate in the later stage.Therefore, increasing the inner diameter of the permeable chamber (d pc ) appropriately can expedite drainage and reduce construction time.

Horizontal Hydraulic Conductivity
The seepage rate is influenced directly by the horizontal hydraulic conductivity (K h ).Therefore, the horizontal hydraulic conductivity (K h ) primarily influences the later stages of vacuum formation and drainage.As shown in Figure 11, the influence of the horizontal hydraulic conductivity (K h ) on the seepage rate is significant, while its effect on the vacuum degree is minor.It shows that the vacuum formation in low-permeability soil is mainly affected by siphon drainage and permeable chamber size.The hydraulic conductivity is an inherent property of the soil and can be challenging to adjust.However, a higher hydraulic conductivity contributes to drainage.

Cylinder Diameter
The cylinder diameter affects drainage by influencing the seepage rate.Its impact pattern is similar to that of the horizontal hydraulic conductivity (Kh).However, the effect of cylinder diameter on seepage is not as significant as the effect of the horizontal hydraulic conductivity (Kh).As shown in Figure 12, the cylinder diameter has almost no impact

Cylinder Diameter
The cylinder diameter affects drainage by influencing the seepage rate.Its impact pattern is similar to that of the horizontal hydraulic conductivity (K h ).However, the effect of cylinder diameter on seepage is not as significant as the effect of the horizontal hydraulic conductivity (K h ).As shown in Figure 12, the cylinder diameter has almost no impact on vacuum formation and only a minor influence on drainage.The cylinder diameter can be determined based on the cost and duration of the project.

Cylinder Diameter
The cylinder diameter affects drainage by influencing the seepage rate.Its impact pattern is similar to that of the horizontal hydraulic conductivity (Kh).However, the effect of cylinder diameter on seepage is not as significant as the effect of the horizontal hydraulic conductivity (Kh).As shown in Figure 12, the cylinder diameter has almost no impact on vacuum formation and only a minor influence on drainage.The cylinder diameter can be determined based on the cost and duration of the project.

Global Sensitivity Analysis Method
The global sensitivity analysis method has the ability to analyze the sensitivity of multiple factors simultaneously.It can extend the range of factor change to the entire domain while allowing each factor to have different ranges of variation and achieve simultaneous changes.Sobol is a widely applied global sensitivity analysis method.Several scholars have already employed this method in various fields such as geotechnical engineering, agriculture, materials science, and environmental studies [25][26][27][28][29].The impact proportion of each factor on the final vacuum degree can be calculated using Sobol.This enables the optimization of design factors in the SVD.

Global Sensitivity Analysis Method
The global sensitivity analysis method has the ability to analyze the sensitivity of multiple factors simultaneously.It can extend the range of factor change to the entire domain while allowing each factor to have different ranges of variation and achieve simultaneous changes.Sobol is a widely applied global sensitivity analysis method.Several scholars have already employed this method in various fields such as geotechnical engineering, agriculture, materials science, and environmental studies [25][26][27][28][29].The impact proportion of each factor on the final vacuum degree can be calculated using Sobol.This enables the optimization of design factors in the SVD.
The core principle of the Sobol is to decompose the function model into individual factors and their combinations.The corresponding sensitivity indices are obtained by calculating the impact of the variance of individual or multiple input factors on the total output variance.The calculation equations for Sobol are as follows [29,30]: The function model f (X) can be decomposed into a sum of 2 p increasing terms to obtain the Equation (11).(11) where f (0) is a constant; X = [x 1 , x 2 ,. .., x i ,. .., x p ], x i is a performance indicator or influence factor, and i = 1, 2,. .., p.
The total variance (V) and bias-variance (V i ) of Sobol are denoted as: The expression for the decomposition of the total variance into the bias-variance is: The Sobol sensitivity indices can be obtained by dividing both sides of Equation ( 14) by V: where S i is the first-order sensitivity index, indicating the influence of the variation of factor i on the model output; S i,j is the second-order sensitivity index, indicating the joint influence of simultaneous variations in factors i and j on the model output; and S i T is the total-order sensitivity index, indicating the sum of sensitivity indices for factor i across all orders.
The sensitivity indices can serve as indicators to evaluate the importance of various factors.Before calculating the sensitivity indices for each factor, a reasonable range of values (minimum to maximum) must be defined for the seven factors under study, as shown in Table 3.Then, the factors are sampled using the sobolset function in MATLAB.The random numbers generated by this sampling method form a low-discrepancy sequence, which is uniform and efficient.Figure 13 shows the 200 randomly collected data points for the thickness of the sealing soil column (h s ), with a uniform data distribution.The sampling dimension for this study is seven, meaning that each set of samples contains seven data points (one for each factor).To eliminate the influence of sample size on the results, a series of sample libraries were obtained by incrementing the sample size at intervals of 50 (50, 100, 150, etc.).First-order and total-order sensitivity indices were calculated for each sample library.

Results and Discussion
As shown in Figure 14, the sensitivity indices calculated in this study converge well.The first-order and total-order sensitivity indices obtained are essentially stable when the sample size exceeds 2000 sets.Table 4 provides the first-order and total-order sensitivity indices for the seven parameters when the sample size is 5000.The first-order sensitivity indices are ranked from highest to lowest as follows: 0.7772 (hs), 0.0731 (d), 0.0618 (L), 0.0419 (Kh), 0.0196 (dpc), 0.0024 (de), and 0.0018 (ζ).The sum of the first-order sensitivity

Results and Discussion
As shown in Figure 14, the sensitivity indices calculated in this study converge well.The first-order and total-order sensitivity indices obtained are essentially stable when the sample size exceeds 2000 sets.Table 4 provides the first-order and total-order sensitivity indices for the seven parameters when the sample size is 5000.The first-order sensitivity indices are ranked from highest to lowest as follows: 0.7772 (h s ), 0.0731 (d), 0.0618 (L), 0.0419 (K h ), 0.0196 (d pc ), 0.0024 (d e ), and 0.0018 (ζ).The sum of the first-order sensitivity indices has reached 0.9779.The total-order sensitivity indices are ranked from highest to lowest as follows: 0.7825 (h s ), 0.0861 (d), 0.0749 (L), 0.0483 (K h ), 0.0273 (d pc ), 0.0027 (d e ), and 0.0024 (ζ).The thickness of the sealing soil column (h s ) has an extremely high impact on the vacuum degree.Its first-order sensitivity index accounts for 79.48%, and the total-order sensitivity index accounts for 76.40%.This significantly exceeds the impact of other factors.Therefore, the thickness of the sealing soil column (h s ) is the factor that needs the most attention in the SVD design.The results of the global sensitivity analysis are consistent with the analysis results in Section 2.3.The diameter of the permeable chamber (d pc ) and the local resistance coefficient of the siphon pipe (ζ) have a very small impact on the final vacuum degree.

Calculation of Maximum Vacuum Degree
The thickness of the sealing soil column (hs) has an impact of approximately 80% on the vacuum degree.Therefore, exploring the quantitative relationship between vacuum degree and sealing soil column thickness (hs) is valuable for assessing the maximum vacuum degree that the SVD can form.

Calculation of Maximum Vacuum Degree
The thickness of the sealing soil column (h s ) has an impact of approximately 80% on the vacuum degree.Therefore, exploring the quantitative relationship between vacuum degree and sealing soil column thickness (h s ) is valuable for assessing the maximum vacuum degree that the SVD can form.The range of values for the thickness of the sealing soil column (h s ) is [0.1, 0.9] m.Secondly, based on the theoretical model, calculate the maximum vacuum degree (Max|p A, t |) under different thicknesses of the sealing soil column (h s ), as shown in Figure 15.Finally, a fitting equation for the calculation results of the maximum vacuum degree is obtained:

Conclusions
Siphon-vacuum drainage is a promising method for marine reclamation, with the potential to reduce the environmental impact and operational costs.In this study, the factors influencing vacuum formation and drainage in the SVD have been comprehensively analyzed.The findings provide a basis for optimizing the design of the SVD in marine reclamation projects, leading to more sustainable and efficient land development along coastal regions.The following conclusions can be drawn:

•
After the SVD is started, a vacuum gradually forms inside the permeable chamber.The duration of the vacuum formation process is generally less than 400 s.During this process, the seepage drainage rate gradually increases.The siphon rate (Uout, t) initially rises rapidly and then decreases slowly, gradually approaching the seepage rate (Uin, t).The siphon drainage plays a dominant role in the water storage change within the permeable chamber during the early stages.

•
The length (L), local resistance coefficient (ζ), and diameter of the siphon pipe (d) influence the early stages of vacuum formation and drainage by affecting the siphon.
Increasing the siphon rate (Uout, t) can lead to a greater vacuum degree and drainage rate.Reducing the thickness of the sealing soil column (hs) can create a larger vacuum degree, resulting in a larger cross-sectional area and drainage rate.The diameter of the permeable chamber (dpc) primarily affects the drainage rate through the crosssectional area.

•
The first-order sensitivity indices are as follows: 0.7772 (the sealing soil column thickness, hs), 0.0731 (the siphon pipe diameter, d), 0.0618 (the siphon pipe length, L), 0.0419 (the horizontal hydraulic conductivity, Kh), 0.0196 (the permeable chamber diameter, dpc), 0.0024 (the cylinder diameter, de), and 0.0018 (the local resistance coefficient, ζ).The sealing soil column thickness (hs) has an extremely high impact on the vacuum degree, accounting for 79.48%.The combined impact of the cylinder diameter (de) and local resistance coefficient (ζ) on the vacuum degree accounts for only The determination coefficient (R 2 ) of the fitting curve is 1, indicating a good fit.When h s = 0.1 m, the vacuum degree can reach 80.49 kPa.The smaller the sealing soil column thickness (h s ), the greater the vacuum degree that can be automatically formed in the SVD.Additionally, it can shorten the path for pore water from the upper soil to seep into the vertical drainage well (permeable chamber).Thus, using smaller thicknesses of the sealing soil column (h s ) can enhance the reinforcement effect on soft soil.Naturally, this should be done while ensuring proper sealing.

Conclusions
Siphon-vacuum drainage is a promising method for marine reclamation, with the potential to reduce the environmental impact and operational costs.In this study, the factors influencing vacuum formation and drainage in the SVD have been comprehensively analyzed.The findings provide a basis for optimizing the design of the SVD in marine reclamation projects, leading to more sustainable and efficient land development along coastal regions.The following conclusions can be drawn: • After the SVD is started, a vacuum gradually forms inside the permeable chamber.
The duration of the vacuum formation process is generally less than 400 s.During this process, the seepage drainage rate gradually increases.The siphon rate (U out, t ) initially rises rapidly and then decreases slowly, gradually approaching the seepage rate (U in, t ).The siphon drainage plays a dominant role in the water storage change within the permeable chamber during the early stages.• The length (L), local resistance coefficient (ζ), and diameter of the siphon pipe (d) influence the early stages of vacuum formation and drainage by affecting the siphon.
Increasing the siphon rate (U out, t ) can lead to a greater vacuum degree and drainage rate.Reducing the thickness of the sealing soil column (h s ) can create a larger vacuum degree, resulting in a larger cross-sectional area and drainage rate.The diameter of the permeable chamber (d pc ) primarily affects the drainage rate through the cross-sectional area.• The first-order sensitivity indices are as follows: 0.7772 (the sealing soil column thickness, h s ), 0.0731 (the siphon pipe diameter, d), 0.0618 (the siphon pipe length, L), 0.0419 (the horizontal hydraulic conductivity, K h ), 0.0196 (the permeable chamber diameter, d pc ), 0.0024 (the cylinder diameter, d e ), and 0.0018 (the local resistance coefficient, ζ).The sealing soil column thickness (h s ) has an extremely high impact on the vacuum degree, accounting for 79.48%.The combined impact of the cylinder diameter (d e ) and local resistance coefficient (ζ) on the vacuum degree accounts for only 0.43%.• A fitting curve was obtained that describes the relationship between the maximum vacuum degree (Max|p A, t |) and different sealing soil column thicknesses (h s ).Calculations indicate that the maximum vacuum degree (Max|p A, t |) can exceed 80 kPa.
Reducing the sealing soil column thickness (h s ) can also shorten the path for pore water in the upper soil to seep into the vertical drainage well.Therefore, using the smallest possible sealing soil column thickness (h s ) while ensuring sealing can enhance the improvement effect.

Figure 3 .
Figure 3. Schematic diagram of a cylindrical model of a single permeable chamber.

Figure 3 .
Figure 3. Schematic diagram of a cylindrical model of a single permeable chamber.

Figure 4 .
Figure 4. Test results and calculations of air pressure changes.

Figure 4 .
Figure 4. Test results and calculations of air pressure changes.

J 19 Figure 5 .
Figure 5. Drainage flow rate of siphon and seepage.

Figure 5 .
Figure 5. Drainage flow rate of siphon and seepage.

JFigure 6 .
Figure 6.(a) Vacuum formation and (b) drainage rate under the influence of the siphon pipe diameter.

Figure 6 .
Figure 6.(a) Vacuum formation and (b) drainage rate under the influence of the siphon pipe diameter.

Figure 6 .
Figure 6.(a) Vacuum formation and (b) drainage rate under the influence of the siphon pipe diameter.

Figure 7 .
Figure 7. (a) Vacuum formation and (b) drainage rate under the influence of the siphon pipe length.

Figure 7 .
Figure 7. (a) Vacuum formation and (b) drainage rate under the influence of the siphon pipe length.

JFigure 8 .
Figure 8.(a) Vacuum formation and (b) drainage rate under the influence of the minor head loss coefficient of siphon pipe.

Figure 8 .
Figure 8.(a) Vacuum formation and (b) drainage rate under the influence of the minor head loss coefficient of siphon pipe.

Figure 9 .
Figure 9. (a) Vacuum formation and (b) drainage rate under the influence of the thickness of sealing soil column.

Figure 10 .
Figure 10.(a) Vacuum formation and (b) drainage rate under the influence of the permeable chamber diameter.

Figure 10 .
Figure 10.(a) Vacuum formation and (b) drainage rate under the influence of the permeable chamber diameter.

JFigure 11 .
Figure 11.(a) Vacuum formation and (b) drainage rate under the influence of the horizontal permeability coefficients.

Figure 11 .
Figure 11.(a) Vacuum formation and (b) drainage rate under the influence of the horizontal permeability coefficients.

Figure 11 .
Figure 11.(a) Vacuum formation and (b) drainage rate under the influence of the horizontal permeability coefficients.

Figure 12 .
Figure 12.(a) Vacuum formation and (b) drainage rate under the influence of the cylinder diameter.

Figure 12 .
Figure 12.(a) Vacuum formation and (b) drainage rate under the influence of the cylinder diameter.

Figure 13 .
Figure 13.Distribution of randomly collected data points for sealing soil column thickness.

Figure 14 .
Figure 14.Sensitivity indices of various factors corresponding to different sample sizes: (a) first-order; (b) total-order.
First, the values of the remaining six factors are determined to be L = 50 m, ζ = 0.131, d = 5 mm, d pc = 5 cm, K h = 3.75 × 10 −9 m/s, and d e = 1.015 m.

Figure 15 .
Figure 15.Maximum vacuum degree achievable with different thicknesses of sealing soil column.

Figure 15 .
Figure 15.Maximum vacuum degree achievable with different thicknesses of sealing soil column.

Table 2 .
Factor values of each scenario.

Table 3 .
Range of factor values.

Table 3 .
Range of factor values.

Diameter of Siphon Pipe, d (mm) Length of Siphon Pipe, L (m) Local Re- sistance Co- efficient of Siphon Pipe, ζ Thickness of Sealing Soil Col- umn, hs (m) Diameter of Permeable Chamber, dpc (cm) Horizontal Hydrau- lic Conductivity, Kh (×10 −9 m/s) Cylinder Diameter, de (m)
Figure 13.Distribution of randomly collected data points for sealing soil column thickness.

Table 4 .
The first-order and total-order sensitivity indices of various factors when the sample size is 5000.

Table 4 .
The first-order and total-order sensitivity indices of various factors when the sample size is 5000.