Prediction Method of Water Absorption of Soft Rock Considering the Influence of Composition, Porosity, and Solute Quantitatively

The study of water absorption of soft rock is of great significance to the prevention engineering disasters. However, the research on the prediction method of water absorption of soft rock considering the influence of composition, porosity, and solute is insufficient. Aiming to solve the above problem, water absorption tests are carried out by synthesizing soft rocks and water solutions. Then, the prediction model of water absorption of soft rock is established through quantitative analysis of water absorption data and compared with the water absorption characteristics of natural rock to verify the reliability of the model. The results show that the changes in water content and water absorption velocity of soft rock with time obey the second-order exponential attenuation function and the linear function (double logarithmic coordinates axis), respectively. The types of cations and anions and the type and content of clay minerals have the greatest influence on the process of rock water absorption. In addition, the water absorption prediction model could better predict the water absorption process of natural rock. The research results solve the problem of insufficient research on the soft rock water absorption prediction method considering the influence of composition, porosity and solute.


Introduction
Water has a great impact on the safety of the rock engineering. For example, seawater intrusion will threaten the construction and operation of subsea tunnels [1], damage the structure of underground sealed oil storage or reservoir [2,3] or cause collapse due to coastline erosion [4]. Moreover, water has a significant impact on the construction and operation of urban rock tunnels [5,6]. Combined with the influence of rock weathering and tunnel construction blasting, water is easier to enter the tunnel through some channels [7][8][9][10][11][12]. Water also affects the safety of slopes and chambers. For example, water will weaken the stability of red bed soft rock bank slope [13]. Water intrusion or erosion mainly occurs in areas with poor rock mass quality, which mainly depends on the mineral composition (especially clay minerals) and porosity of rock mass, solute composition and high salinity of water solutions, and the action of external environment [14,15]. The prediction of water absorption of soft rock is very important for the failure prediction of rock engineering. If the critical time point of rock damage due to the change in water absorption cannot be judged, the prevention and control of rock engineering disasters will be too conservative or too risky [16]. However, existing research studies mainly focus on the weakening effect of water on rock physical and mechanical properties [17,18], or natural rock samples and natural water solutions are used to study the water absorption of soft rock. Under the influence of anisotropic factors, the influence of composition, porosity and solute on the water absorption of soft rock cannot be quantitatively analyzed [19]. Therefore, this paper needs to study the test method and prediction model of water absorption of soft rock considering the effects of composition, porosity and solute quantitatively.
At present, there are many studies on rock water absorption. Some scholars have studied the effect of solute composition of water on rock water absorption. For example, Tinnacher et al. studied the adsorption of rocks rich in clay minerals in water with different mineralization [20]. Manaka et al. studied the effects of solute type and solution concentration in solution on chemical permeability, diffusion and hydraulic parameters of mudstone, and found that solute type has little effect on rock water absorption, while solution concentration has a great effect on rock water absorption [21]. Chai et al. studied the water absorption characteristics of rocks in solutions with different chemical properties and obtained the effects of different chemical paths on the water absorption process of clay rocks [22]. Some scholars mainly pay attention to the influence of rock porosity on its water absorption. For example, Binal and Wang et al. studied the relationship between physical properties such as porosity and water absorption of rocks and mechanical properties, and found that the greater the porosity, the greater the water absorption of rocks [23,24]. Tertre et al. studied the water (solution) absorption characteristics of natural clay-based materials with different porosities and found that osmotic swelling and its influence on pore structure are the key influencing factors of the rock water absorption process [25]. Simonyan studied the effect of alteration on rock water absorption and found that the formation of secondary facies changed the rock pore space (effective porosity) and played a leading role in controlling rock water absorption [26]. Moreover, the rock porosity can be obtained by ultrasonic velocity, image processing technology, neural network method, fractional flow theory, etc. [27][28][29][30][31][32][33]. Some scholars also mainly study the influence of clay minerals on rock water absorption. For example, Erguler et al. studied the change of water absorption over time of rocks containing different types and contents of clay minerals (siltstone, mudstone, marl, tuff) [34]. Rahromostaqim et al. studied the hydration expansion characteristics of pure clay and mixed clay and their effects on water absorption [35,36]. However, few scholars have quantitatively studied the effects of rock composition, porosity and solute on water absorption at the same time. Feng et al. studied the water absorption of shale with different clay mineral types and pore structures in water with different salinity, but did not quantitatively build a water absorption prediction model [37].
Therefore, according to the particularity of rock engineering, this paper studies the influence of soft rock composition, porosity and solute on water absorption, and establishes a water absorption prediction model. In order to eliminate the influence of anisotropy of rock and water solutions, firstly, rock samples with different clay mineral composition and porosity are synthesized and water solutions with different salinity are configured. Then, the water absorption tests of rock samples are carried out, and the effects of different factors on water absorption are quantitatively analyzed. Finally, the prediction model of water absorption of soft rock is established and its reliability is verified.

Test Materials
There are many kinds of ions in water solutions. The common cations are mainly Na + , Ca 2+ and K + , and the common anions are mainly Cl − , SO 2− 4 and CO 2− 3 [38]. The variation range of ion concentration and pH value in water solutions is large [39]. Rock engineering involves many sedimentary rocks (sandstone, mudstone and shale, etc.), which contain a large number of soft rocks. Soft rocks have special engineering properties, such as low strength, large porosity, poor cementation and high content of clay minerals (montmorillonite, illite and kaolinite, etc.), which have a great negative impact on engineering safety [40,41]. The interaction of the above conditions will have a complex impact on the water absorption of soft rock.
In order to eliminate the influence of anisotropy of soft rock and water solutions, soft rock analogues [42] and water solutions are synthesized to study. Synthetic soft rock analogues require clay minerals, quartz sand (particle size 0.075~0.1 mm) and pure Figure 1. The X-ray diffraction topographies of the clay minerals.

Test Methods
Soft rock analogues and water solutions are synthesized as shown in synthesizing soft rock analogues, the total mass of clay minerals and quar the mass of water is 5 g. Place the mixture of the above materials evenly i mold (diameter 50 mm) and hold the pressure for 10 min with an uniaxia testing machine according to the vertical stress shown in Table 1 (Figure sample out of the mold and put it into an oven heated to 105 °C until its constant value. The sample is then cooled to ambient temperature in a desi in order to prevent the sample turning muddy during water absorption, it porous stones (diameter 50 mm, thickness 5 mm) and heat shrinkable tub

Test Methods
Soft rock analogues and water solutions are synthesized as shown in Table 1. When synthesizing soft rock analogues, the total mass of clay minerals and quartz is 100 g and the mass of water is 5 g. Place the mixture of the above materials evenly in a cylindrical mold (diameter 50 mm) and hold the pressure for 10 min with an uniaxial compression testing machine according to the vertical stress shown in Table 1 (Figure 2a). Take the sample out of the mold and put it into an oven heated to 105 • C until its mass reaches a constant value. The sample is then cooled to ambient temperature in a desiccator. Finally, in order to prevent the sample turning muddy during water absorption, it is sealed with porous stones (diameter 50 mm, thickness 5 mm) and heat shrinkable tubes. The middle part of the bottom surface is made into a circle (diameter 38 mm) for the water absorption tests ( Figure 2b). In this study, 18 samples are synthesized for tests.
In order to determine the porosity of the samples by a nondestructive monitoring method, the P-wave velocity V p (km/s) of synthetic soft rock analogues is measured by C61 non-metallic ultrasonic detector. The porosity of the test sample could be obtained by the relationship between the P-wave velocity and the rock porosity n (%) (Equation (1), Figure 3) established by Yagiz [27].
Then, weigh a certain mass of salt according to the columns of CTs, ATs and CCs in Table 1. Put the salt in deionized water and add HCl and NaOH (considering the contents of Na + and Cl − ) to prepare solutions with certain pH values and mix evenly. Finally, 100 mL of solution supernatant is used as the test solution. In this study, 18 water solutions are synthesized for tests.
The water absorption test system is composed of a container, sample, water solution, high-precision electronic scale (Figure 2c), porous stone, etc. The sample is not in direct contact with the solution, but through contact with porous stone to reduce the saturation velocity of the sample (Figure 2d). Each sample is tested for regular intervals. A total of 32 tests are carried out for a total time of 4980 min (The samples are saturated or nearly saturated at 4980 min, and the extension of the test time has no effect on the main conclusions of this study). In each test, the total mass of the water absorption container including the sample should be weighed first and then the remaining mass of the water absorption container without the sample should be weighed. The change in water absorption in each test is calculated according to Equation (2), so as to eliminate the experimental error caused by solution evaporation. After the test, remove the heat shrinkable tube and porous stones on the surface of the sample and cut them from the middle. It can be seen that the water absorption of the sample is very uniform (Figure 2e).
where, ∆w i (%) is the change in water absorption when the test is conducted to i, i and i-1 are the test numbers, M (g) is the total mass of the water absorption system including the sample, m (g) is the residual mass of the water absorption system without the sample and M s (g) is the initial mass of the sample. In order to determine the porosity of the samples by a nondestructive monitoring method, the P-wave velocity (km/s) of synthetic soft rock analogues is measured by C61 non-metallic ultrasonic detector. The porosity of the test sample could be obtained by the relationship between the P-wave velocity and the rock porosity (%) (Equation (1), Figure 3) established by Yagiz [27]. In order to determine the porosity of the samples by a nondestructive monitoring method, the P-wave velocity (km/s) of synthetic soft rock analogues is measured by C61 non-metallic ultrasonic detector. The porosity of the test sample could be obtained by the relationship between the P-wave velocity and the rock porosity (%) (Equation (1), Figure 3) established by Yagiz [27].

Change in Water Absorption with Time
In order to analyze the water absorption process of soft rock, the results (Table S1) of 18 group tests in Table 1 are averaged to eliminate the influence of various factors. The water absorption velocity is obtained by time differentiation of the water content. Then, by fitting the values of water content and water absorption velocity, the best fitting curve is obtained and the rock water absorption process model is established (Figures 4 and 5).
where, ∆ (%) is the change in water absorption when the test is conducted to , and -1 are the test numbers, (g) is the total mass of the water absorption system including the sample, (g) is the residual mass of the water absorption system without the sample and (g) is the initial mass of the sample.

Change in Water Absorption with Time
In order to analyze the water absorption process of soft rock, the results (Table S1) of 18 group tests in Table 1 are averaged to eliminate the influence of various factors. The water absorption velocity is obtained by time differentiation of the water content. Then, by fitting the values of water content and water absorption velocity, the best fitting curve is obtained and the rock water absorption process model is established (Figures 4 and 5).  According to the average water content, a second-order exponential decay function is fitted (Equation (3), the corrected coefficient of determination = 0.992), and the rock water content model is obtained. The overall trend of the curve is that it first rises rapidly, then slows down and finally stabilizes near saturation. According to the parameters and in Equation (3), the whole curve can be divided into three sections. (1) Section I: extremely fast water absorption section; the water content exceeds 50% of the saturated water content in a very short time, and the curve is linear and almost vertical.
(2) Section II: deceleration water absorption section; within the time equivalent to several tens of times of Section I, the water content reaches 80% of the saturated water content and the curve is concave. (3) Section III: slow water absorption section; compared with the According to the average water content, a second-order exponential decay function is fitted (Equation (3), the corrected coefficient of determination R 2 adj = 0.992), and the rock water content model is obtained. The overall trend of the curve is that it first rises rapidly, then slows down and finally stabilizes near saturation. According to the parameters t 1 and t 2 in Equation (3), the whole curve can be divided into three sections. (1) Section I: extremely fast water absorption section; the water content exceeds 50% of the saturated water content in a very short time, and the curve is linear and almost vertical. (2) Section II: deceleration water absorption section; within the time equivalent to several tens of times of Section I, the water content reaches 80% of the saturated water content and the curve is concave. (3) Section III: slow water absorption section; compared with the first two sections, the water absorption time was significantly increased and the change rate of water absorption velocity is significantly reduced.
where, w is the water content. w 0 is the offset, which is close to the saturated water content, and the value here is 17.6. t is the water absorption time. According to the average water absorption velocity, a linear function is fitted (Equation (4), the corrected coefficient of determination R 2 adj = 0.983) in the double logarithmic coordinate axis (base 10), and the rock water absorption velocity model is obtained. Additionally, according to the parameters t 1 and t 2 in Equation (3), the whole curve is divided into three sections. (1) Section I: the water absorption velocity is greater than 0.4%/min. (2) Section II: the water absorption velocity is between 0.004-0.4%/min. (3) Section III: the water absorption velocity is less than 0.004 %/min. Equation (4) can also be written in the form of Equation (5).
where y represents lg(v). x represents lg(t). a is the intercept and is always positive, and the value here is 0.54. b is the slope and is always negative, and the value here is −1.14.
where v is the water absorption velocity. Under the premise of considering the scale problem, the above laboratory test results can be applied to engineering. Reasonable measures are taken to prevent and control engineering disasters by judging the section of rock water absorption. When the water absorption of the rock is in Section I, the water content is small and the water absorption velocity is large, and the water source needs to be cut off in time to prevent the water content from increasing rapidly. When the water absorption of the rock is in Section II, the water can be cut off while absorbing water. At this time, the operability time for taking engineering measures is sufficient. When the water absorption of the rock is in Section III, the influence of water can be ignored in a short period of time because the rock is basically saturated and the project has not been damaged. This kind of situation is suitable for short-term use of rock engineering. For example, on the premise of exploring the fault and other failure surfaces of Norwegian subsea tunnel, the rock water absorption model proposed in this study can be used to judge the water absorption stage of the surrounding rock during tunnel construction and certain measures can be taken to prevent and control tunnel water inflow [16].

Influence of Various Factors on the Water Absorption Process of Soft Rock
In rock engineering, water absorption by soft rock is the result of the combined action of internal and external factors. In different geological environments, the influence of different factors on the water absorption process needs to be considered. This section analyzes the experimental results for seven factors: clay mineral type, clay mineral content, vertical stress (porosity), cation type, anion type, concentration and pH value. Under the influence of various factors, the water content and water absorption velocity of rock are shown in Figures 6 and 7, respectively. The experimental design in Table 1 in this study was designed according to the seven factors and three levels in the orthogonal experiment, so the data sources in Figures 6 and 7 are the average of each factor and each level [43]. For example, the average value of water content for test numbers 1, 2, 3, 10, 11 and 12 in Table 1 is the effect of montmorillonite in the clay mineral type on the water absorption process. Meanwhile, the data in Figures 6 and 7 are fitted with Equations (3) and (4), respectively (most of the R 2 adj in the fitting results are greater than 0.95, indicating that the fitting results of water content and water absorption velocity over time are good). The extraction of the fitting results and comparison of the effects of each level on the water absorption characteristics are displayed in Table 2. The response data in Table 2 are normalized to obtain the degree of influence (DI L ) of each level on the rock water absorption process, as shown in Figure 8.          By analyzing Table 2 and Figure 8, it can be seen that seven factors have different degrees of influence on rock water absorption. Among the fitting parameters, A 1 and A 2 are less affected by seven factors, while other fitting parameters are greatly affected by seven factors. The specific surface area and swellability of montmorillonite are much higher than those of illite and kaolinite, which leads to slower water absorption velocity in rocks containing montmorillonite (Sections I and II last longer, t 1 and t 2 are larger, a is smaller) and the change rate of water absorption velocity (b) is smaller. The effect of clay mineral content on rock water absorption cannot be ignored. With the increase in clay mineral content, the duration of Section I (t 1 ) shortened, the duration of Section II (t 2 ) increased and the change rate of water absorption velocity decreased (b). The clay mineral type and clay mineral content have little effect on the saturated water content (w 0 ) of the rock. Generally speaking, the smaller the vertical stress is when the rock is synthesized, the higher the rock porosity, the better the rock permeability and the greater the water content and water absorption velocity [44,45]. However, when the vertical stress is too large, it will cause secondary fracture of the rock, and although its porosity is small, its permeability is very good [42]. Therefore, in this study, when the vertical stress is small and large (5 MPa and 15 MPa), the rock water absorption velocity is fast and the saturated water content is high. Regarding the influence of water solution on the water absorption process of rocks, the research mainly includes the interaction of ion type, concentration and pH value with clay minerals in rocks [22,26]. In this study, the saturated water content and the change rate of the water absorption velocity of the rock in the monovalent cation-containing solution are larger, and the solution containing Ca 2+ makes the rock water absorption velocity slower. Among the anions, the saturated water content of the rock in the SO 2− 4 -containing solution is the smallest, and the change rate of the water absorption velocity is the largest, while the water absorption velocity of the rock in the CO 2− 3 -containing solution is the slowest. The higher the solution concentration, the slower the rock water absorption velocity, the greater the saturated water content, and the smaller the change rate of the water absorption velocity. In the acidic solution, the rock has the fastest water absorption velocity, the largest saturated water content and the largest change rate of the water absorption velocity. The relevant indicators in the neutral and alkaline solutions are relatively close and small. Predecessors have also analyzed the above research contents and obtained fruitful results. However, there are still two problems in these studies: (1) the effects of clay minerals, porosity and solutes on the water absorption process are not comprehensively considered and (2) the water content model and water absorption velocity model are not used to quantitatively analyze the water absorption process [46][47][48].

Quantification of the Degree of Influence of Each Factor
The research in this section focuses on the degree of influence of seven factors (CMT, CMC, VSs, CTs, ATs, CC, and pH) on rock water content and water absorption velocity. According to the principle of orthogonal test analysis and the design of orthogonal table, the range analysis method is selected to analyze the degree of influence of various factors on rock water absorption [49]. The range values for all responses during water absorption process are calculated according to Equation (6), and the results are shown in Table 3.
where R is the range value. The greater the value, the greater the impact. K max and K min are the maximum and minimum values within each factor of the representative t 1 , t 2 , w 0 and b in the equation parameters in Table 2. Table 3. Range value of the seven factors and degree of influence on the water absorption process of rocks. The larger the range value of the factor, the higher the priority of a factor's degree of influence. PDI refers to the priority of degree of influence. From Table 3, the rock porosity (VSs) and the type of anions (ATs) in the solution have the greatest influence on the water absorption time (t 1 and t 2 ), and the solution concentration (CCs) has the least influence on the water absorption time. The concentration of the solution had the greatest influence on the saturated water content (w 0 ), and the clay mineral type (CMT) had the least influence on the saturated water content. The clay mineral type has the greatest influence on the change rate of water absorption velocity (b), and the rock porosity has the least influence on the change rate of water absorption velocity.

Factor
In order to eliminate the order of magnitude difference between the four parameters and facilitate the analysis, the range values are normalized and the degree of influence of the seven factors is drawn, and the overall degree of influence of each factor on the water absorption process is calculated (Figure 9, R n is the normalized degree of influence of each factor on each parameter, and R t is the sum of the degree of influence of each factor on the four parameters). It is found that seven factors have different effects on each parameter of the water absorption process. The total influence (R t ) of the seven factors on the water absorption process of soft rock is ATs, CMT, CMC, CTs, VSs, pH and CCs from the largest to the smallest. That is, the types of anions and cations in water solutions and clay mineral compositions in rocks have the greatest influence on rock water absorption process. Predecessors also studied the influence of these seven factors on rock water absorption, but did not use the water content model and water absorption velocity model to quantify its influence [37]. Appl. Sci. 2022, 12, x FOR PEER REVIEW 13 of 17 clay mineral compositions in rocks have the greatest influence on rock water absorption process. Predecessors also studied the influence of these seven factors on rock water absorption, but did not use the water content model and water absorption velocity model to quantify its influence [37].

Water Absorption Prediction Model of Soft Rock
On the basis of the above analysis, a prediction model of soft rock water absorption is established. The model can be used to predict the water absorption process of rocks, given the type and content of clay minerals, porosity and properties of water solutions. Firstly, build the factor and level matrices (Equations (7) and (8)). In the factor matrix, the sum of the values equals 7. If all factors are present, then each value is equal to 1. If only some factors exist, for example, only VSs and pH exist, then = = 3.5, and the remaining values are equal to 0. Similarly, in a level matrix, the values of each column sum to 3. For the first column, all values are 1 if all three clay minerals are present. If only one clay mineral, illite, exists, then = 3 and = = 0. Then, establish the degree of influence matrices of the factor and level matrices (Equations (9) and (10)). The values in the matrix are the results of the normalization of each factor in Figure  9, and the values in the matrix are the results of the normalization of the degrees of influence of each level in Figure 8. Finally, the comprehensive influence factor ( ) is obtained through matrix operations (Equation (11)).
where is the factor matrix, is the value of each factor and its value range is 0~7.
where is the level matrix, is the value of each level and the value range is 0~3.

Water Absorption Prediction Model of Soft Rock
On the basis of the above analysis, a prediction model of soft rock water absorption is established. The model can be used to predict the water absorption process of rocks, given the type and content of clay minerals, porosity and properties of water solutions. Firstly, build the factor and level matrices (Equations (7) and (8)). In the factor matrix, the sum of the f values equals 7. If all factors are present, then each f value is equal to 1. If only some factors exist, for example, only VSs and pH exist, then f VSs = f pH = 3.5, and the remaining f values are equal to 0. Similarly, in a level matrix, the l values of each column sum to 3. For the first column, all l values are 1 if all three clay minerals are present. If only one clay mineral, illite, exists, then l ILL = 3 and l MMT = l KAO = 0. Then, establish the degree of influence matrices of the factor and level matrices (Equations (9) and (10)). The values in the matrix DI F are the results of the normalization of each factor R t in Figure 9, and the values in the matrix DI L are the results of the normalization of the degrees of influence of each level in Figure 8. Finally, the comprehensive influence factor (p) is obtained through matrix operations (Equation (11)).
where F is the factor matrix, f is the value of each factor and its value range is 0~7.
where L. is the level matrix, l is the value of each level and the value range is 0~3.
where p is the comprehensive influence factor, and diag is the operator that takes the main diagonal value of the matrix. Since the interaction between factors and levels is not considered, only the diagonal of the matrix is meaningful. Substituting p into Equations (3) and (5) yields a multifactor water absorption prediction model (Equations (12) and (13)).
In order to verify the reliability of the soft rock water absorption prediction model, this model is compared with the related test results of the hard rock in Figure 1 of the literature [50] and the soft rock in Figure 5 of the literature [51], as shown in Figure 10. In general, the water absorption prediction model can better predict the water absorption process of rocks under the premise of known clay minerals, rock porosity and water solution properties. The curve calculated by the water absorption prediction model is basically consistent with the change trend of previous research results (except for the sample w-304). Moreover, the model can not only be used to predict the water absorption process of soft rock, but can also be used to predict the water absorption process of hard rock. Due to different experimental materials, methods and conditions, the data of some water absorption processes are different, but this does not affect the judgment of the critical time point when rocks are damaged due to changes in water content. Previous researchers also studied the water absorption prediction model, such as Zhu et al.'s entropy-based method to predict the water absorption process of rock, which has a very good prediction effect under certain conditions. However, the influence of clay minerals and solutes on the water absorption process was not considered in their study [52]. Appl where is the comprehensive influence factor, and is the operator that takes the main diagonal value of the matrix. Since the interaction between factors and levels is not considered, only the diagonal of the matrix is meaningful.
In order to verify the reliability of the soft rock water absorption prediction model, this model is compared with the related test results of the hard rock in Figure 1 of the literature [50] and the soft rock in Figure 5 of the literature [51], as shown in Figure 10. In general, the water absorption prediction model can better predict the water absorption process of rocks under the premise of known clay minerals, rock porosity and water solution properties. The curve calculated by the water absorption prediction model is basically consistent with the change trend of previous research results (except for the sample -304). Moreover, the model can not only be used to predict the water absorption process of soft rock, but can also be used to predict the water absorption process of hard rock. Due to different experimental materials, methods and conditions, the data of some water absorption processes are different, but this does not affect the judgment of the critical time point when rocks are damaged due to changes in water content. Previous researchers also studied the water absorption prediction model, such as Zhu et al.'s entropy-based method to predict the water absorption process of rock, which has a very good prediction effect under certain conditions. However, the influence of clay minerals and solutes on the water absorption process was not considered in their study [52].

Conclusions
In this study, the soft rock analogs and water solutions were synthesized, water absorption tests were carried out and a soft rock water absorption prediction model was established. The changes in water content with time in the process of rock water absorption obeys a second-order exponential decay function, and the change in water absorption velocity with time obeys a linear function in the double logarithmic coordinate axis. Among the seven factors, the type of anions and cations in solution and the type and content of clay minerals in rock have the greatest influence on rock water absorption. In addition, in comparison with the water absorption process of natural rocks, the reliability of the soft rock water absorption prediction model is verified and it is proved that the model can predict the water absorption characteristics of natural rocks. The related results solve the problem of insufficient research on the prediction method of soft rock water absorption that quantitatively considers the influence of composition, porosity and solute.
Supplementary Materials: The following supporting information can be downloaded at: www.mdpi.com/xxx/s1, Table S1: The water absorption experimental results of 18 samples in Table 1.

Conclusions
In this study, the soft rock analogs and water solutions were synthesized, water absorption tests were carried out and a soft rock water absorption prediction model was established. The changes in water content with time in the process of rock water absorption obeys a second-order exponential decay function, and the change in water absorption velocity with time obeys a linear function in the double logarithmic coordinate axis. Among the seven factors, the type of anions and cations in solution and the type and content of clay minerals in rock have the greatest influence on rock water absorption. In addition, in comparison with the water absorption process of natural rocks, the reliability of the soft rock water absorption prediction model is verified and it is proved that the model can predict the water absorption characteristics of natural rocks. The related results solve the problem of insufficient research on the prediction method of soft rock water absorption that quantitatively considers the influence of composition, porosity and solute.

Conflicts of Interest:
The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.