Abstract

The creep characteristics of soft clays have been studied for decades. However, the lateral deformation of soils is not allowed during the commonly used one-dimensional consolidation tests, which cannot describe the real deformation features of soils in practice. On the other hand, the influence of drainage distance on the mechanical properties of soil is still controversial, classified as hypothesis A and hypothesis B. For a better understanding of deformation characteristics of soft clay, especially which in long-terms, a series of conventional oedometer tests as well as novel geometric confined consolidation tests was conducted on soft dredger fill. The results show that the secondary consolidation coefficient of the soil sample Cα would increase firstly, followed by a small decrease with the increase of consolidation pressure generally. Cα would decrease with the consolidation time and also be reduced by preloading. The strain at the completion of primary consolidation would increase with the drainage distance, but the Cα would be affected little. Both compression index Cc and Cα of soft clay would reduce after preconsolidation, in which two parameters show an approximate linear relationship. The creep coefficient of soft clay under the geometric confinement is larger than that under the oedometer test Cαε. However, the trends of the relationship between the creep coefficient and loading are consistent regardless of the confinement conditions.

1. Introduction

Soft clay could be widely found in nature, which has high water content, low permeability, high compressibility, and high viscosity. Such features could result in an excessive deformation, especially for a long period for constructions on soft clay. The deformation of soft clay has a time effect involving both consolidation and creep, which interact with each other and contribute to long-term deformation. Oedometer tests have been widely used to characterize the deformation of soft clay in long-term, which is closely related to mineral composition, magnitude of load, loading rate, temperature, water content, stress history, etc. [110].

Due to the different conditions such as sedimentary environment, mineral composition, and stress history, the physical characteristics and mechanical behaviour of soft clay in different regions would vary a lot. For example, secondary consolidation coefficient Cα for soft clay in different regions might increase [11, 12], remain the same [13, 14], or increase first then decrease [5, 6, 1517] with increasing of loading. When the soil state changes from overconsolidated to normally consolidated, Cα would increase firstly followed by a decreasing with the increasing of load. The upper load corresponding to the maximum Cα may appear near the preconsolidation pressure [6, 18] or equal to 1.5 to 4 times the preconsolidation pressure [1517]. When the soil is in the normal consolidation state, Cα would not change with the consolidation pressure anymore [14]. Cα is also closely related to the loading ratio. When the loading ratio decreases, Cα would decrease due to the quasi-overconsolidation phenomenon caused by creep [16, 19, 20]. Mesri [21] revealed a positive correlation between Cα and compression index Cc for several typical soft clays. Subsequently, Santagate et al. [22] and Jesmani et al. [4] also reached the consistent conclusion in their studies.

Mainstream theories to explain whether the creep effect occurs after the excess pore water pressure has dissipated or it already begins at the primary consolidation stage were classified into Hypothesis A and Hypothesis B, according to Ladd et al. [23]. The scholars supporting Hypothesis A [2429] believe that the creep deformation occurs after the completion of primary consolidation. For the same soil mass, the strain value at the completion of primary consolidation is unique and independent of thickness of soil sample or drainage boundary conditions. Mesri and Choi [26] conducted a series of oedometer tests under different drainage conditions, and they agreed that creep existed in the primary consolidation stage. However, they proposed that the relationship between the end-of-primary (E-O-P) void ratio and effective vertical stress of soil samples with different thicknesses was still unique due to the interaction between consolidation effects and creep effects in the primary consolidation stage, which was a proof of the Hypothesis A. In Hypothesis B, creep deformation has occurred in the primary consolidation stage. For the same kind of soil, the strain value at the completion of primary consolidation is related to the thickness of soil sample as well as the drainage boundary conditions, e.g., more creep deformation occurs in the primary consolidation stage if the distance of drainage path is larger. The majority of laboratory as well as in situ test results support the Hypothesis B [3036]. For the contradiction of the two hypotheses, Kabbaj et al. [35] pointed out that the soil sample used in tests for Mesri and Choi [26] was obtained by using a soil collector with a diameter of 50 mm. In this case, the soil sample was disturbed significantly, which could lead to an underestimate of the practical settlement prediction by laboratory test results. However, an agreement is not being reached on this issue yet.

For a better understanding of deformation characteristics for soft clay, especially for that in long term, this study conducted a series of one-dimensional confined consolidation tests on soft clay in the eastern Shantou City of China, including the conventional oedometer tests and novel geometric confined consolidation tests, which could simulate the real deformation characteristics of soft clay in practice. During the improved geometric confined consolidation tests, vertical and limited lateral deformations of soil would occur at the same time. Thus, the deformation parameters considering the influence of lateral deformation could be obtained. Based on the test results, the effect of loading magnitude, preloading, loading time, thickness of soil sample, and lateral constraint conditions on the long-term deformation of soft clay was systematically studied and discussed.

2. Materials, Apparatus, and Testing Procedure

The soil tested in this study is soft dredger fill, which is taken from the eastern economic zone in Shantou City, Guangdong Province, PRC (Figure 1). A thin-walled soil sampler with a height of 50 cm and a diameter of 10 cm was used to obtain the undisturbed soil samples at a distance of 7.0 to 9.0 m below the ground, which was below the groundwater level (Figure 2). The sampling and transportation procedures were carefully done to avoid any disturbance to the intact soil samples [37, 38]. The soil samples were muddy soft clay with grey-black colour. It could be found that the soft clay was mixed with a small amount of fine sand or chips of shells locally, which should be removed from soil samples in order to ensure the accuracy and uniformity of test results [39, 40]. The results of particle size analysis showed that the content of clay (d < 5 μm) accounted for about 37.4%, the content of cement (d < 2 μm) accounted for about 18.7%, the content of sand particles (d > 75 μm) was about 14%, and the content of organic matter was approximately 4%. The physical properties of the soil samples are shown in Table 1. The soil samples were saturated with a high initial void ratio and water content. According to the Casagrande method, the preconsolidation pressures of the samples are 48 to 53 kPa, while the overburden pressure is about 160 kPa, so that the soil is in an unconsolidated state [41].

A series of oedometer tests was conducted on soil samples with a same area of 30 cm2 and thickness of 2 cm, 4 cm, 6 cm, and 10 cm, respectively. The samples with double-drainage boundaries were tested in a constant temperature of 24 ± 1°C to avoid the scatters of test results caused by temperature change. Details of all the oedometer tests under step loading are summarized in Table 2. It is noted that, for test no. 2, the soil sample was firstly preloaded for 3 days under 100 kPa to investigate the effect of preloading on the mechanical behaviour of the sample. Besides, preloading of 25 kPa was also applied on test no.7–no.10 for two days to avoid the scatter caused by the sampling procedure of remoulded samples. After the preloading was removed and rebounding completed, the step loading listed in Table 2 was performed on the soil samples.

Another kind of test conducted in this study was the geometric confined consolidation test using an improved GJY K0 consolidation apparatus made by Nanjing Soil Instrument, PRC (Figure 3). For the improved apparatus, the soil sample was separated from the pressure chamber by a rubber membrane with a thickness of 3 mm. Meanwhile, the pressure chamber was filled with deaired water and connected with the confining pressure controlling system to obtain a constant lateral confining pressure. Mechanical behaviour of soil sample, which has a thickness of 4 cm and an area of 30 cm3 with double-drainage boundaries, was studied during step loading. In this case, a parameter named “geometric confined coefficient” was proposed, which was defined as  = , where σ1 is the vertical stress and σ3 is the confining pressure. remained a constant throughout the test, e.g., 0.6, smaller than the earth pressure coefficient at rest K0, i.e., 0.74∼0.78, to allow a certain lateral deformation of soil sample under loading. The temperature during the tests was consistent with that for oedometer tests as 24 ± 1°C. The axial displacement and lateral pressure are recorded by a dial gauge and load cells, respectively. The details of the two geometric confined consolidation tests are listed in Table 2.

3. Results and Discussion

The Casagrande method is commonly used to distinguish primary consolidation and secondary consolidation for one-dimensional consolidation tests (Figure 4). In this method, the slope of the e − logt curve for secondary consolidation is considered as the value of Cα. According to Leroueil et al. [42] and Mesri and Godlewski [43], the values of the Cα would decrease with time. However, Cα could remain constant at least in the first two logarithmic cycles. Therefore, in order to avoid the influence of different logarithmic cycles on the determination of Cα, the first logarithmic cycle after the completion of primary consolidation was selected to determine Cα in this study. The equation to obtain Cα is as follows:where t1 is the time when the primary consolidation completed, t2 is the time when the first logarithmic cycle completed, and e1 and e2 are the void ratios of the soil sample at t1 and t2, respectively.

In oedometer tests, lateral deformation of soil sample is not allowed and the area of the soil sample remains constant. Therefore, the void ratio e could be simply derived from vertical cumulative strain. If equation (1) is described in terms of vertical strain ε, a parameter named creep coefficient for one-dimensional consolidation Cαε is proposed:where t1 is the time when the primary consolidation completed for one-dimensional consolidation, t2 is the time when the first logarithmic cycle completed for one-dimensional consolidation, and ε1, ε2 are the vertical cumulative strains at t1 and t2, respectively.

In geometric confined consolidation tests (e.g., test nos. 13 and 14), the lateral constraint of soil sample is flexible, which allows a certain lateral deformation of soil samples. In this case, the area of soil sample has been increasing and void ratio is difficult to measure or calculated from vertical strain. Therefore, the secondary consolidation coefficient Cα is difficult to be obtained and equation (1) is invalid. The creep coefficient for geometric confined consolidation is introduced, which is derived from the vertical strain as follows:where is the time when the primary consolidation completed for geometric confined consolidation, is the time when the first logarithmic cycle completed for geometric confined consolidation and and are the vertical cumulative strains at and , respectively. The introduction of creep coefficients Cαε and makes the creep properties of soil samples under different lateral constraint more comparable.

Figure 5 shows the relationship between the void ratio e and lgt for oedometer test no. 6. It could be observed that when the consolidation pressure is small, e.g., less than 50 kPa, the soil sample is in the overconsolidation state so that the inflection point on the e − lgt curve is not obvious and the value of Cα is small. As the consolidation pressure gradually increases, the soil sample turns to normal consolidation state and the vertical deformation of the sample would increase significantly. In this case, the e − lgt curve shows a reverse “S” and an obvious inflection point occurs. Moreover, the slope of e − lgt curve in the secondary consolidation stage has little change with the load in the normal consolidation state.

Chen’s method [44, 45] was used to analyze the results of soil samples under the step-loading, in order to obtain the deformation-time curves under a single-loading and then the isochronous curves of stress-strain relationship (Figure 6). It could be observed that the stress-strain relationships of soil samples are a series of curves with a same trend and the deformation of the soil exhibits a nonlinear characteristic. However, when the upper load is small (less than 63 kPa) indicating the soil is in the overconsolidated state, the stress-strain relationship of the soil sample is approximately linear. With the increase of the load, the soil turns to the normal consolidation stage and the stress-strain relationship curves of the soil samples are deflected to the stress axis and present significant nonlinear characteristics.

3.1. Relationship between Secondary Consolidation Coefficient and Loading for Oedometer Tests

In this study, for test nos. 1, 11, 13, and 14, the value of Cα increases first and then slightly decreases with the increasing of loading. It could reach a maximum value when the consolidation pressure is about η· (η ≈ 3.2–3.5; is the preconsolidation pressure), as shown in Figure 7(a), taking test no. 1 as an example. For the other oedometer tests, e.g., test no. 4, Cα would increase continually with the consolidation pressure increasing and the slope of the Cα −  curve is related to the preconsolidation pressure (Figure 7(b)). It could be seen that when the loading is less than η· (η ≈ 2.3), Cα increased dramatically with the load increasing. Whereas when the loading exceeds η· (η ≈ 2.3), the rise of Cα would decline.

Figure 7 shows that the inflection point of Cα −  curve does not occur near the preconsolidation pressure but the value of η· does. This result is consistent with the previous research such as [15] and [46], in which the value of η was about 1.5 to 3.5. However, Chen et al. [3] proposed that the inflection point of Ca −  curve would occur near based on their test results. The reasons for this contradiction need to be further studied, but it might be attributed to the quasi-overconsolidation due to the loading time being longer than the primary consolidation completion time or the load increment ratio.

According to the mechanism of the variation of Ca with the consolidation pressure, the Ca- curves in Figure 7 could be divided into three stages. In stage I ( < ), the soil is in the overconsolidated state. The soil structure formed during natural deposition would restrict the development of secondary consolidation. In this case, the rearrangement of soil particles and the creep deformation of the bound water would be restricted so that the value of Cα is small. Based on Bjerrum’s theory of equitime e-lg curve, Yin et al. [19] developed a method to calculate the secondary consolidation deformation of overconsolidated soils as follows:where is the time increment, and the equivalent creep time can be calculated bywhere is the time when the primary consolidation completed, is the preconsolidation pressure, is the pressure after loading, which is less than , and is the compression index.

It can be found that the secondary consolidation deformation in the overconsolidation stage is related to the consolidation pressure. When the consolidation pressure increases to , the Cα −  curve would turn into stage II, in which the cementations between the soil particles are gradually destroyed with the loading increasing. Then, soil particles would be rearranged along the force direction to their new equilibrium positions. Meanwhile, the intermolecular force between bound water and particles would decrease, leading to an increase of the secondary consolidation effect. However, when the load continues to increase to exceed η·, the structures between soil particles and between bound water and particles could be destroyed completely, resulting in a decrease of the volume of voids inside the soil. The soil sample approaches to the remoulded state, and the increment of secondary consolidation deformation of the sample declines with the load and tends to be stable gradually, which is defined as stage III. For test nos. 1 and 11, due to their relatively high water content and macrovoid ratio, the intermolecular forces between soil particles and between the bound water and particles would decease dramatically. Therefore, the rearrangement of soil particles in stage II is more remarkable and the secondary consolidation effect is more obvious, leading to a longer stage II and a larger Cα even compared with that for stage III. In test nos. 13 and 14, a certain lateral deformation was allowed due to the flexible lateral confinement during the test, resulting in more rearrangements of soil particles along with a longer stage II.

3.2. Effect of Preloading on Secondary Consolidation Coefficient

For test no. 2, the soil sample was preloaded for 3 days under 100 kPa and then unloaded. The step loading in a sequence of 25 kPa–50 kPa–100 kPa–200 kPa–400 kPa–800 kPa would not be performed on the soil sample until the rebounding displacement did not change any more. Figure 8 illustrates that the relationship between Cα and for test no. 2 compared with that for the soil sample without preloading (test no. 3).

It could be concluded that Cα of soil sample after preloading tends to be smaller than that for soil sample without preloading under the same consolidation pressure, whereas the reduction differs under the different loading levels. For example, under the load of 100 kPa, Cα for the sample after preloading reduces by nearly 50.9% compared with that for the sample without preloading, while this reduction is only 3.84% under 800 kPa. The effect of preloading on Cα increases first and then decreases with the increasing of the consolidation pressure. It could be observed from Figure 8 that when consolidation pressure is 100 kPa, the preloading has the greatest impact on the secondary consolidation effect. The above results could be explained by theory of equitime e − lg curve. If , the soil without preloading was in the overconsolidation state. The value of Cα was small, and the preloading did not change the consolidation state of the soil sample. Therefore, the change of the secondary consolidation effect caused by the preloading is not remarkable, e.g., from point A to point B in Figure 9 (point B is on the rebound recompression curve). If the soil sample was preloaded for 3 days under 100 kPa, the deformation of the soil would creep from point C to point D, and it would rebound along a straight line l when unloaded. When it was reloaded again, the yield point of sample would then develop to point E. It is noted that the preconsolidation pressure at point E was greater than preloading of 100 kPa, resulting in a quasi-overconsolidation of the soil sample. For a soil sample without preloading, it was in a normal consolidation state under a load of 100 kPa. While after preloading, the consolidation state of the soil changed to overconsolidation state due to quasi-overconsolidation phenomenon caused by the creep effect under the load of 100 kPa. The difference of consolidation state of the soil sample caused by preloading leads to a notable variation of . When the consolidation pressure was greater than , the soil samples were in a normal consolidation state no matter with or without preloading. As the consolidation pressure increased, the relationship between e and would develop along the normal consolidation line (line tc) as shown in Figure 9. In this situation, it could be considered that the decrease of the initial void ratio of the soil sample after the preloading is the reason of the decrease of Cα.

The results also show that the reduction of the compression index Cc of soil sample caused by preloading at the step loading of 100 kPa was about 79.2%, which was about 1.5 times of the reduction of Cα under the same condition. In practice, preloading is commonly used for soft clay in foundation treatment. This study could explain the reasonability and effectiveness of this method from the mechanism and confirmed that preloading could prevent more deformations of primary consolidation under the same loading as the prepressure.

3.3. Variation of Secondary Consolidation Coefficient with Consolidation Time

For test nos. 4, 5, and 6, soil samples were consolidated for 1 day, 2 days, and 3 days, respectively. The values of Cα under different consolidation pressures as well as different consolidation time periods are listed in Table 3. It could be seen from Table 3 that the consolidation time cannot affect the trend of the variation of Cα with the consolidation pressure. However, Cα would decrease slightly with the increasing of consolidation time under the same consolidation pressure and the effect of consolidation time on Cα would increase with the increasing of the consolidation pressure. The longer the consolidation time, the larger the creep deformation and the smaller the void ratio of soil sample, resulting in a decreasing of the value of Cα. As mentioned above, the value of Cα for soil sample in overconsolidated state is smaller than that for soil sample in the normal consolidated state. The step loading was conducted in this study so that the impact of the consolidation time on Cα would be accumulated. Therefore, the reduction of Cα for soil sample consolidated for 3 days compared with that for 1 day under 800 kPa is about 2.8 times of the reduction under 50 kPa.

3.4. Effect of Drainage Distance on Secondary Consolidation Coefficient

In order to study the influence of drainage distance on Cα, a series of oedometer tests were conducted on soil samples with different thicknesses of 2 cm, 4 cm, 6 cm, and 10 cm (no. 7 to no. 10 soil samples, respectively). The samples were remoulded soils, and the test results were analyzed using Chen’s method. Figure 10 shows the ε − lgt curves of soil samples with different thicknesses under a load of 100 kPa. It can be found that, with the decrease of the thickness of the sample, ε − lgt curve tends to move downward, despite of coincidence of curves during the secondary consolidation state. Meanwhile, the strain at the completion of primary consolidation would also decrease with the decrease of the thickness of the sample, i.e., 7.15%, 6.6%, 6.25%, and 5.75% for the thickness of 10 cm, 6 cm, 4 cm, and 2 cm, respectively. The deformation characteristics are consistent with Hypothesis B [23], in which the creep already occurs in the primary consolidation state. The strain at the completion of primary consolidation would vary with the thickness of the soil sample, i.e., the strain would decrease with the decreasing of the thickness of the sample. This is because the thinner the soil sample, the shorter completion time of the primary consolidation and the less creep deformation during the primary consolidation. This study has further confirmed that the E-O-P curve of the soil is related to the thickness of the soil sample. Therefore, if the deformation parameters obtained from the laboratory tests would be used in the assessment of the settlement in practice, the effect of thickness of soil layers should be taken into consideration.

Figure 11 shows the relationship between Cα and the thickness of the soil sample under different consolidation pressure conditions. It is noted that Cα for soil samples with different thicknesses changed little under the same consolidation pressure. This is because, although the strain of the thick soil sample is larger at the completion of the primary consolidation, the primary consolidation time is also longer. According to equation (3), the effect of thickness of soil samples on the value of Cα is slightly under the same loading step.

3.5. Effect of Confinement on the Creep Coefficient

As previously described, the creep coefficients of soil samples in one-dimensional consolidation as well as in geometric confined consolidation could be obtained from equations (2) and (3), respectively. For oedometer tests, the relationship between the creep coefficients Cαε and the secondary consolidation coefficient Cα is shown as follows:where e0 is the initial void ratio of the soil sample.

Figure 12 illustrates the relationship between the creep coefficient and consolidation loading for both undisturbed soil samples and remoulded soil samples under different confinements (specimens used in test nos. 11 and 13 were undisturbed soil while those in test nos. 12 and 14 were remoulded soil). It is noted that both the creep coefficient for oedometer tests Cαε and that for geometric confined consolidation test would increase firstly with the load increasing, followed by a slight decrease and then tend to be stable for undisturbed soil (Figure 12(a)). The same trend could be found for remoulded soil as shown in Figure 12(b). This is because, under geometric confinement, the lateral deformation of soil sample results in more rearrangements of soil particles and distortions of bound water inside the specimen. Besides, the creep deformation of the sample in this case concludes two parts: volume creep and shear creep. By contrast, for oedometer tests, as the lateral deformation is restricted rigidly, the creep deformation is only associated with volume creep, i.e., the secondary consolidation deformation.

In geometric confined consolidation tests, the ratio of confining pressure to the vertical load remained a constant as 0.6. Although the confinement was not rigid, there was still a constraint of lateral deformation of the specimen to a certain extent. Therefore, the volume creep still dominated the total creep while the shear creep played a secondary role and no shear creep failure occurred, so that the trend of creep coefficient with load increasing for geometric confined consolidation tests changed little compared with that in oedometer tests (Figure 10).

3.6. Relationship between Secondary Consolidation Coefficient and Compression Index

Walker [47] conducted a series of oedometer tests on Leda soft clay. The results indicated that the ratio of the secondary consolidation coefficient Cα to the compression index Cc was approximately a constant about 0.025. Subsequently, Mesri and Choi [26] investigated the relationship between Cα and Cc for 22 kinds of soft clays. They proposed that the value of Cα/Cc remained almost a constant for the same kind of clay. This ratio could range from 0.025 to 0.1 for different kinds of soft clays. In this study, the results of 15 oedometer tests (thickness of 2 cm for soil samples) were analyzed, and then the ratios of Cα/Cc for soft clay in the eastern Shantou City were obtained. The value of Cα/Cc ranges from 0.021 to 0.045, which roughly agrees with the conclusion of Mesri and Choi [26]. Little scatters of data might contribute to experimental errors or errors during secondary consolidation analysis. The relationships between Cα and Cc for the 15 tests were plotted in Figure 13. The results present a linear relationship between Cα and Cc approximately, and the average ratio of Cα/Cc is about 0.032.

3.7. Characteristics of Microstructure of Soil Samples after Consolidation

Compaction characteristics in terms of microstructure of soil samples after consolidation for 15 days under a load of 100 kPa could be observed by the scanning electron microscope (SEM). Figure 14 presents the images of soil samples before and after oedometer tests (at a magnification of 1000 times and 5000 times, respectively). It could be seen from Figures 14(a) and 14(c) that there are more pores and higher void ratio inside the undisturbed soft clay. The microstructure is mainly loose flocculation structure and soil particles in a kind of slice state distributed loosely. The soil particles are aligned along the direction of gravity slightly but not markedly. In this case, the interconnections between soil particles are mostly edge-to-edge or edge-to-surface contacts, for which the cementations are weak. By contrast, the images of soil samples with different magnifications after consolidation (Figures 14(b) and 14(d)) show a compact microstructure of samples, resulting in an obvious reduction of void ratio in terms of macroscopic scales. The soil particles rearrange in the direction perpendicular to the consolidation pressure to be more close and even overlap with each other. The connection between particles turns to surface-to-surface contact and the microstructure of the soil sample tends to be a stable clot-like structure.

4. Conclusions

In this study, a series of oedometer tests and geometric confined consolidation tests were conducted on soft clay samples in eastern Shantou City, China. The deformation characteristics of consolidation, especially in long-term consolidation, were systematically investigated for soft clay. The results could be concluded as follows:(1)Generally, the secondary consolidation coefficient Cα of the soft clay would continuously increase with the load increasing, although the increment could reduce gradually. However, Cα would increase firstly to the maximum value, followed by a decrease with the increase of the load for soil samples with a high initial void ratio. A turning point occurs at the point where the loading approaches to 1.5–3.5 times the preconsolidation pressure.(2)Both compression index and secondary consolidation coefficient of soft clay would reduce after preconsolidation while the reductions could be different under different loadings. The effect of preconsolidation is more dramatic on the compression index than that on the secondary consolidation coefficient.(3)The secondary consolidation coefficient would decrease with the increase of the consolidation time, being hardly affected by the drainage distance.(4)The creep coefficient of soft clay under the geometric confinement is larger than that under the oedometer test. However, the trends of the relationship between the creep coefficient and loading are consistent regardless of the confinement conditions. The main reason is that, although shear creep occurs during geometric confined consolidation, the majority of creep deformation is still caused by the volumetric creep, which contributes to the entire creep deformation of soil samples in one-dimensional consolidation.(5)An approximate linear correlation can be observed between the secondary consolidation coefficient and compression index. The value of Cα/Cc ranging from 0.021 to 0.045 falls roughly within the variation range proposed by Mesri and Choi [26].

Data Availability

The data used to support the findings of this study are included within the article.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

This work was supported by the National Natural Science Foundation of China (Grant nos. 51809050, 41902288, and 51978177) and Natural Science Foundation of Guangdong Province (Grant no. 2020A1515010811).