The development of the fissures in soluble rock of karst areas directly affects the construction and operation safety of tunnel engineering. It is thus of theoretical and practical significance to study the characteristics of its corrosion and its influencing factors. Taking the Wulong tunnel as the research object, the numerical model of the study area was established to quantitatively analyze the corrosion range, corrosion ratio, and changes in the permeability and porosity of the fissures in soluble rock of karst areas of the tunnel over the past 100 years, and the simulation results were verified by field experiments. The results show that the main controlling factor of the fissure corrosion of the tunnel in the karst area is the flow rate. The corrosion range and corrosion ratio of the fissures of the tunnels in the karst area increased with temperature because the reaction rate constant increased with temperature, causing the reactions’ equilibrium to move towards the direction of the solution. The larger the initial permeability and the larger the porosity of the fissures, the faster the fissures corrode. In the same time period, the fissures with high permeability and large porosity will lead to the permeability and porosity being more enhanced, thus causing the corrosion of the fissures to exhibit secondary enhancement effects. The opening of the dead-end pores greatly enhanced the permeability and slightly increased the porosity, which caused the differential corrosion of fissures in the karst area. The protection of the tunnel should be strengthened, mainly in strong hydrodynamic conditions and in the fracture development zone.
China National Science and Technology2016ZX05016-005National Natural Science Foundation of China41372334Jilin University20171481. Introduction
Carbonate rocks account for 35% of the land area distribution of China [1, 2]. The karst development is very strong in Southwest China, which often led to geohazards and a variety of accidents of engineering construction [3–8]. Due to the distribution of abundant mountains in Southwest China, the ratio of bridges and tunnels relative to roads or railways is mostly over 70%, some up to 90% [2, 9, 10]. The karst caves and subterranean rivers that are broadly distributed in karst regions often lead to accidents during tunnel construction and increase the cost of the tunnels’ operation and maintenance [9–11]. The Wulong tunnel is one of the representative tunnel projects in the karst area of Southwest China. According to a preliminary investigation, karst doline, karst caves, and subterranean rivers were found in the tunnel crossing area [2, 9, 10].
Elucidating the corrosion mechanism of carbonate rocks helps us analyze the corrosion development of the fissures in soluble rock of karst areas. Scholars have carried out many studies on the influence factors of karst development through experiments, numerical simulations, and other methods [12–24]. These studies showed that the corrosion of carbonate rocks was influenced by internal factors, such as geological structures and rock compositions, as well as external factors, such as flow rates (rainfall), temperature, and pH.
In terms of the impact of fissures on karst corrosion, scholars have found that faults and fractures play a vital role in karst corrosion [25–27]. They found that in the unidirectional one-dimensional limestone fissures, the corrosion rate of the limestone fissures shows a nonlinear increase with time [16, 28–33]. At the same time, scholars have also studied the characteristics of fissures and pore evolution under two-dimensional [14] and three-dimensional conditions [17, 34–37]. In addition, by studying the changes in the microstructures of the carbonate rocks that were corroded, it was determined that corrosion is first determined by the lithology and its own structure [38–40].
As for the influence of external environmental factors, such as flow rates, on karst corrosion, which are based on the basic principles of mineral dissolution, scholars obtained the basic differential equation of calcite dissolution rates through experimental research and mathematical statistics [15, 41]; based on the corrosion law and surface corrosion theory of carbonate rocks, it was concluded that the flow rate, surface area, and concentration of major ions obviously affect the dissolution [22, 23, 42–44]. To assess the temperature factor, carbonate rock corrosion experiments between 0°C and 250°C were conducted, using rate equations for the corresponding temperature conditions to discuss the influence of temperature on the dissolution, precipitation, and surface corrosion of carbonate minerals [13, 19, 21, 23, 45–51]. In terms of pH, under natural conditions, the pH change in rain water is mostly caused by changes in the CO2 content of air [22, 23, 28, 52]. Scholars found that the variation of pCO2 in the range of 1~50 atm obviously affected the corrosion of carbonate minerals [22, 23], while at pCO2 values of lower than 0.1 atm, pCO2 had very little influence on the corrosion of carbonate minerals [28, 52, 53].
For tunnel engineering in a karst area, the geological and hydrogeological natural conditions must be taken into consideration [54]. The systematic study of the corrosion mechanism by numerical simulation and field experiments is of great importance for tunnel operation and maintenance.
In this paper, the Wulong tunnel, a typical tunnel project in Southwest China, was studied. The corrosion range, corrosion ratio, and changes in the porosity and permeability of karst cracks in the tunnel vault within a timeframe of 100 years were analyzed, and the differences in corrosion were discussed.
2. The Background of the Wulong Tunnel Engineering Area
The Wulong tunnel is located in Wulong County of China (Figure 1) with a full length of 9418 m; it stretches along with the Wujiang Canyon and goes through the Wuling Mountains, with a maximum burial depth of approximately 830 m. The middle part of the tunnel passes through 3800 m of soluble rock. Many karst caves and subterranean rivers have been found in the Wulong tunnel engineering area. Therefore, the engineering geological and hydrogeological conditions in the tunnel area are very complicated.
The location of the study area and tunnel project.
Wulong County is subtropical, with four distinct seasons, rain and heat over the same period, temperatures between −3.5°C and 41.7°C, and annual precipitation of approximately 1000–1200 mm.
Carbonate rocks and clastic rocks are predominant in the Wulong tunnel area. For the lithology of the tunnel from entrance to exit, see Table 1 and Figure 2.
Formation lithology of the Wulong tunnel area.
System
Formation
Lithology
Jurassic
Zhenzhuchong
J1z
Dark-brown shale and mudstone, light-gray sandstone
Triassic
Xujiahe
T3xj
Dark-brown shale, light-gray sandstone
Leikoupo
T2l
The upper part is dark-brown shale and marlstone. The middle part is dark-brown shale and light-gray sandstone. The lower part is gray limestone
Jialingjiang
T1j
Gray limestone intercalated with dolomitic limestone
Feixianguan
T1f
Gray limestone, dark-brown shale
Permian
Changxing
P2c
Dark-gray limestone
Wujiaping
P2w
Dark-gray limestone
Maokou
P1m
Gray and dark-gray limestone
Liangshan
P1q+l
Limestone intercalated with bioclastic limestone
Silurian
Luoreping
S1lr
Gray and dark-gray mudstone
Simplified hydrogeological cross section of the Wulong tunnel area.
The study site is located within the Maokou limestone formation (P1m), and it is close to the Wujiaping Formation (P2w), which is a set of coal-bearing strata near the contact zone of insoluble and soluble rock. A large “Tiankeng,” a karst doline (Yangaotuo doline), can be found on the ground above the tunnel (Figures 3 and 4(a)). The karst phenomenon is very developed near the area of the tunnel, and karst caves can be found in nearby strata (Figure 4(b)).
Geological map of the Wulong tunnel area.
Karst phenomena in the study area ((a) karst doline; (b) karst caves).
In the study area, groundwater mainly receives rainfall infiltration, and the recharge, runoff, and discharge are controlled by topography, lithology, and faults. Karst groundwater flows to the Wujiang River in the form of subterranean rivers and springs in the tunnel area. The 2# and 3# subterranean rivers pass through the middle section of the tunnel area, which is the most karst-developed section of the whole tunnel, having karst caves and karst funnels (Figures 2–4). The 2# river is located 4496 m from the tunnel entrance (Figure 2). The 3# river is located 4643 m from the tunnel entrance. The two subterranean rivers are both nearly orthogonal with the tunnel line. The 2# and 3# subterranean rivers inrushed into the tunnel many times during the tunnel construction, which caused great difficulties in the tunnel engineering construction and maintenance [55].
Karst fractures increase the contact area of water and carbonate rock, which is helpful to promote the corrosion of carbonate rock, and have a certain effect on the stability of tunnel bedrock. Moreover, the tunnel vault is one of the most vulnerable parts of a tunnel, and it is prone to cause tunnel accidents. Therefore, in this study, the most vulnerable area, the tunnel vault fracture zone, which is located under the Yangaotuo karst funnel and above the 2# and 3# subterranean rivers, was chosen as the key study area (Figures 2 and 3).
3. The Fissure Evolution Model and Simulation Scheme in the Karst Area of the Tunnel3.1. Conceptual Model of the Fissures in the Karst Area of the Tunnel
According to field investigations, the geological/hydrogeological conditions, length, hydrochemical characteristics, and flow rates of 2# and 3# subterranean rivers were similar. A conceptual model was established for the study area, the main fissure, and the damage zone, which is above the 2# subterranean river (Figure 5). Based on previous investigations about the impact zones of different sizes of faults and fissures [56–61], the widths of the impact zones are generally 5–10 meters. In this study, the length and width of the impact zone of the tunnel fracture zone are all set to 24 meters, and the height of the simulated area was set as 60 meters based on the saturated zone above the vault of the tunnel. In field conditions, groundwater infiltrates from the top fracture zone, flowing through fissures and finally converging into the subterranean rivers. In the model, the upper and lower boundaries of the model were set as constant flow boundaries, and the lateral boundaries were set as no flow boundaries (Figure 5). In the horizontal direction, the length of the long fissures is approximately 80 m and the length of the short fissure is approximately 5 m (Figure 5).
Conceptual model of the study area.
3.2. Numerical Model
The numerical model was constructed based on the conceptual model. The modeling area was divided into 10 layers in the vertical direction, with each layer being 6 meters high, and the horizontal direction was divided into 6 columns, with each column being 4 meters wide; overall, the numerical model area was divided into 360 grids (Figure 6). The central part of the model is the fracture core zone; moving outward from this zone are the damage zone and protolith zone successively. In addition, the permeability and porosity decrease from the central part of the fracture zone to both sides outward (Figure 6).
Numerical model.
3.3. Model Parameters and Simulation Scheme
According to previous research data [9, 10], the permeability of the karst fracture zone in the simulation area is between 1.8 × 10−11 m2 and 1.3 × 10−10 m2. The average permeability of 7.4 × 10−11 m2 was taken as the permeability of the fracture core zone. To study the impacts of porosity and permeability on corrosion, a high-permeability condition of 7.4 × 10−10 m2 and a low-permeability condition of 7.4 × 10−12 m2 were additionally assigned to the fracture core zone in our study based on its range of empirical values. The permeability of the protolith is approximately 1 × 10−17 m2~2.8 × 10−13 m2 according to the measured data of the porosity and permeability of the carbonatite (protolith) [62–64]. Therefore, in our model, the permeability of the protolith was set as 1.5 × 10−14 m2. Because the damage zone is distributed around the fracture core zone, the permeability in the model gradually decreases from the fracture core zone outward to the protolith zone.
According to the tunnel investigation department and at the same time the reference of the empirical value in Southwest China, the porosity of the protolith is approximately 7% and the porosity of the fracture core zone is approximately 20% to 35%. In the model, the porosity of the fracture core zone was set as 35%, 29%, and 20%, and the corresponding permeability values were set as 7.4 × 10−10 m2, 7.4 × 10−11 m2, and 7.4 × 10−12 m2, respectively. The porosity and permeability of the protolith were set as 7% and 1.5 × 10−14 m2, respectively. The mineral composition and density of the protolith were determined by the measured data of calcite (95%), illite (2%), and montmorillonite (3%). The chemical composition of groundwater was determined by water sample analysis, and its pH is 7.3. The water pressure range of the simulation area was assigned according to the field conditions (Table 1). The atmosphere temperature in the study area is between −3.5°C and ~41.7°C, and the temperature of the groundwater is between 10°C and 25°C. The temperatures of the injection water were set as 10°C, 15°C, 20°C, 25°C, 30°C, and 35°C in our model. The model parameters are shown in Table 2. The karst ratios of different tunnel sections (Table 3) were obtained from analyzing the number of the fissures in soluble rock of karst areas and the karstification level [9, 10].
Model parameters.
Category
Density (kg/m3)
Porosity
Water pressure (Pa)
Permeability (m2)
Protolith
2680
7%
1 × 105~6 × 105
1.5 × 10−14
Fracture core
2000~2280
20%~35%
7.9~790 × 10−12
Karst ratio of the Wulong tunnel.
The starting position of each section from the tunnel inlet (m)
The length of the section (m)
Karst ratio
0~3396
3396
0.2
3396~3892
496
0.4
3892~5276
1384
0.5
5276~7276
2000
0.3
7276~9418
2142
0.2
The simulation area is in the middle section of the tunnel (3892 m~5276 m), where the karst development is the strongest. Because the karst development depends on water-rock action, the flow rates directly influence the development of karst. According to the field investigation of the tunnel, based on the water balance principle [9, 10, 65], the annual rainfall infiltration amount of the whole tunnel area was calculated. The water quantity of the whole tunnel area was then allocated to the different sections of the tunnel using the karst ratios in Table 3; thus, the runoff distribution coefficient ki was obtained. ki can be calculated by
(1)ki=ri×li∑i=1nri×li,where ri is the karst ratio and li is the length of the tunnel section.
From the data in Table 3, the flow rates of different tunnel sections were calculated by (2). The flow rate of the simulated area in the model (Qz) was finally obtained by (3):
(2)Qi=ki×Q,(3)Qz=kz×Qi,where Q is the precipitation infiltration in the whole tunnel area and kz is the ratio of the simulation area size to the length of the most karst-developed parts of the tunnel.
According to the field conditions, the simulation time length was set as 100 years in consideration of the service period of the tunnel. The injection flow rates were based on the rainfall of the dry year, the normal flow year, and the wet year, that is, 800 mm, 1100 mm, and 1500 mm in Wulong, respectively. Using (1), (2), and (3), the injection flow rates of the weak hydrodynamic condition (a), the middle hydrodynamic condition (b), and the strong hydrodynamic condition (c) were calculated to be 1.9 kg/s, 2.9 kg/s, and 3.9 kg/s, respectively. There were 54 different conditions considered to evaluate the evolution of fissure corrosion (Table 4).
Simulation scenarios.
Scenario code
Permeability (m2)
Porosity
Temperature (°C)
Flow rates (kg/s)
a
b
c
Case 1
7.4 × 10−10
35%
10
1.9
2.9
3.9
Case 2
15
Case 3
20
Case 4
25
Case 5
30
Case 6
35
Case 7
7.4 × 10−11
29%
10
1.9
2.9
3.9
Case 8
15
Case 9
20
Case 10
25
Case 11
30
Case 12
35
Case 13
7.4 × 10−12
20%
10
1.9
2.9
3.9
Case 14
15
Case 15
20
Case 16
25
Case 17
30
Case 18
35
The simulation software used in this study was TOUGHREACT, which is widely used in the numerical simulation of water-rock interactions [66, 67].
4. Results and Discussion4.1. Corrosion Range under Different Conditions
The conditions of medium hydrodynamics and medium permeability were taken as examples (case 7(b)~case 12(b)) with which the corrosion process of fissures under different temperatures was analyzed (Figures 7–9). By comparing Figure 7 with Figure 9, during the same corrosion time, the corrosion range of 35°C was larger than that of 10°C. After 100 years of corrosion, the horizontal and vertical corrosion ranges of 35°C were 1.2 times and 1.3 times higher than those of 10°C, respectively.
The corrosion range of the fissure under medium hydrodynamic conditions-medium permeability—10°C.
The corrosion range of the fissure under medium hydrodynamic conditions-medium permeability—20°C.
The corrosion range of the fissure under medium hydrodynamic conditions-medium permeability—35°C.
The condition of medium hydrodynamics and the temperature of 20°C were taken as examples (case 3(b), case 9(b), and case 15(b)) with which the corrosion process of tunnel fractures under different porosity and permeability conditions was analyzed (Figures 10–12). By comparing Figure 10 with Figure 12, during the same corrosion time, the corrosion range of high-porosity and high-permeability conditions was larger than that of low-porosity and low-permeability conditions. After 100 years of corrosion, the horizontal and vertical corrosion ranges of high-porosity and high-permeability conditions were 1.4 times and 1.6 times higher than those of low-porosity and low-permeability conditions, respectively.
The corrosion range of the fissure under medium hydrodynamic conditions-low permeability—20°C.
The corrosion range of the fissure under medium hydrodynamic conditions-medium permeability—20°C.
The corrosion range of the fissure under medium hydrodynamic conditions-high permeability—20°C.
The conditions of 20°C and medium permeability were taken as examples (case 9(a), case 9(b), and case 9(c)) with which the corrosion process of the tunnel fissure under different flow rates was analyzed (Figures 13–15). By comparing Figure 13 with Figure 15, during the same corrosion time, the corrosion range of the strong hydrodynamic condition was larger than that of the weak hydrodynamic condition. After 100 years of corrosion, the corrosion ranges of the strong hydrodynamic condition were 1.8 times and 2.2 times higher than those of the weak hydrodynamic condition, respectively.
The corrosion range of the fissure under weak hydrodynamic conditions-medium permeability—20°C.
The corrosion range of the fissure under medium hydrodynamic conditions-medium permeability—20°C.
The corrosion range of the fissure under strong hydrodynamic conditions-medium permeability—20°C.
4.2. Corrosion Ratios under Different Conditions
Taking the high-porosity/permeability (case 1~case 6) and low-porosity/permeability conditions (case 13~case 18) as examples, we analyzed the corrosion ratios of the tunnel fissure under different flow rates and temperatures (Figures 16(a) and 17(b)). From Figures 16(a) and 17(b), we can see that the corrosion ratio increased with the time of corrosion. During the same corrosion time, the higher the temperature, the greater the corrosion ratio. After 100 years of corrosion, the corrosion ratio (strong hydrodynamic condition and high porosity and permeability) of 35°C (0.4032%) was 1.2 times higher than that of 10°C (0.3369%).
The corrosion ratio of the fissure under strong hydrodynamic conditions over time.
Low permeability
High permeability
The corrosion ratio of the fissure under weak hydrodynamic conditions over time.
Low permeability
High permeability
The comparisons of Figure 16(a) with Figure 16(b) and Figure 17(a) with Figure 17(b), during the same corrosion time, reveal that the corrosion ratio of the tunnel fissure under the high-permeability condition (7.4 × 10−10 m2) was greater than that under the low-permeability condition (7.4 × 10−12 m2). After 100 years of corrosion, the corrosion ratio (strong hydrodynamic condition, 35°C, and high porosity and permeability) (0.4032%) was 1.6 times higher than that of low porosity and permeability (0.2599%).
The comparison of Figure 16(a) with Figure 17(a) and Figure 16(b) with Figure 17(b), during the same corrosion time, reveals that under strong hydrodynamic conditions, the corrosion ratio was greater than that under the weak hydrodynamic condition. After 100 years of corrosion, the maximum corrosion ratio (0.4032%) under the strong hydrodynamic conditions (high porosity and permeability, 35°C) was 2.1 times than that under the weak hydrodynamic condition (0.1908%).
4.3. Evolution of Porosity and Permeability in the Fissure Zone
The structure of the fissures in soluble rock of karst areas is the foundation of the corrosion. The porosity and permeability conditions in fissure zones are the main factors affecting the development of karst. Through the analysis of the effects of different temperatures and flow rates on corrosion, it was found that the most significant corrosion occurred at the temperature of 35°C and strong hydrodynamic conditions. Therefore, we selected 35°C and strong hydrodynamic conditions as an example with which the variation of the porosity and permeability of the fissure zone within 100 years was discussed.
According to the simulation results of 10–100 years of corrosion, the porosity and permeability changes in the fracture core zone of the karst development area were obtained. The relation curves of the increase of permeability with time are shown in Figure 18. As shown in Figure 18, the rate of change of permeability also increased with time. Permeability increased much faster under the condition of high permeability than under the condition of low permeability. After 100 years of corrosion, the permeability increased from 7.40 × 10−12 m2 to 7.73 × 10−12 m2; the maximum increase of permeability was 0.33 × 10−12 m2 under the low-permeability condition. Meanwhile, under the high-permeability condition, the permeability increased from 740.00 × 10−12 m2 to 767.5 × 10−12 m2 and the maximum increase of permeability was 27.5 × 10−12 m2. It was 83.3 times higher than that under the low-permeability condition.
The rate of increase of permeability over time.
As shown in Figure 19, under the conditions of different porosities, the porosity increased with the corrosion time (linear relationship). Under the condition of high porosity, porosity increased much faster than that under the low-porosity condition. After 100 years of corrosion, the porosity of the low-porosity condition increased from 0.200 to 0.203, while under the high-porosity condition, the porosity increased from 0.350 to 0.354. The maximum increase of porosity under the high-porosity condition was 2 times higher than that under the low-porosity condition.
The increase of porosity over time.
The reasons for the changes in the corrosion range, porosity, and permeability of the fissure zone are discussed as follows:
The lithology of the study area is carbonate rocks comprising 95% calcite; the corrosion of carbonate rocks was mainly caused by calcite corrosion, and the chemical reaction equations are as follows:
(4)CO2g+H2O⇌CO2aq+H2OH2CO3⇌H++HCO3−HCO3−⇌H++CO32−H2O⇌H++OH−OH−+CO2aq⇌HCO3−CaCO3+H+⇌CaHCO3+CaHCO3+⇌Ca2++HCO3−CaCO3⇌Ca2++CO32−
Because calcite corrosion’s chemical balance will be affected by temperature, in the range of simulated temperatures (10°C~35°C), the reaction rate constant increased with temperature, causing the reaction’s equilibrium to move towards the direction of the solution. These results have the same corrosion trend as those obtained by Pokrovsky et al. [23] and Huang [68]. They used an experimental method to explore corrosion under neutral conditions and different temperatures (0.5°C~40°C).
According to (4), during the corrosion of calcite, Ca2+ was the main product of the chemical reaction. Under weak hydrodynamic conditions, the flow through the fissures was less and the travel time of water through the fissures was longer. Therefore, the concentration of dissolved Ca2+ was higher in the water, and when the water flow reached the lower part of the fissure, the concentration of Ca2+ in the water was close to saturation. Therefore, corrosion mainly occurred in the upper part of the fissure under the weak hydrodynamic condition, and the corrosion ratio was relatively small. In contrast, in the strong hydrodynamic conditions, the flow rate through the fissures was faster and the travel time of water through the fissures was shorter than that under the weak hydrodynamic condition. Therefore, the concentration of dissolved Ca2+ was lower in water under the strong hydrodynamic condition, and when the water seepage reached the lower part of the fracture, the concentration of Ca2+ was lower. The resulting corrosion of the fissures was more obvious than that of those in the weak hydrodynamic condition, and the corrosion range runs through the whole fissure with a higher corrosion ratio. In addition, due to the increased flow rate, the mechanical failure of the rock is enhanced, which makes the surface of the rock broken. The rock powder on the broken surface of the rock was carried away by the water, thus increasing the corrosion ratio and corrosion range.
Under the conditions of high porosity and permeability, the porosity and permeability increased faster than those under the low-porosity and low-permeability condition. After 100 years of corrosion, the permeability change under the high-permeability condition was 83.3 times higher than that under the low-permeability condition, while the change of porosity under the high-permeability condition was only 2 times higher than that under the low-permeability condition. It can be speculated that with continuing corrosion, more and more dead-end pores became effective pores. Due to the transformation of dead-end pores into effective pores, the connectivity between pores was significantly enhanced and the increase of fissure permeability became significantly larger, while the porosity increased only slightly. Therefore, the degree of change in permeability was much greater than that in porosity in the same period.
According to the simulation results, the larger the porosity and permeability were, the more significant the increases in porosity and permeability were, which led to the enhancement of corrosion in the high-porosity and high-permeability zones and enhancement the differential corrosion of karst.
5. Comparison of the Field Corrosion Experiment Results to the Numerical Simulation Results
To further explore the corrosion characteristics of rock samples at the tunnel’s key karst development parts and verify the reliability of the simulation results, a field corrosion experiment was carried out. The field experimental site was a gully in the Wulong tunnel construction area, which is in the key development parts of the tunnel karst area. The width of the gully is approximately 2.5 meters (Figure 20), and the lithology of the underlayer stratum is thick-layer limestone. The gully is covered with clay, gravel, and plants. The burial depth of groundwater was small.
The site of the field experiment.
The experimental rock samples were made of the limestone of the Changxing Formation, with a size of 30 mm × 10 mm × 10 mm (Figure 21(b)). The rock samples were placed into a plastic box with holes that allowed groundwater to flow through it. The plastic box was buried underground below the groundwater level (Figures 20 and 22).
Rock samples. (a) Original experimental samples, (b) rock samples after being cut, and (c) SEM photograph (×1500).
The process of the field experiment.
After one year of corrosion, the box was taken out and the degrees of corrosion of the rock samples (Figure 20) were measured. Because the corrosion of rock samples depends on the surface area of the water-rock interaction, we used (5) to calculate the surface corrosion ratio R (g/(m2·a)).
(5)R=ms∗t,where R=m/st is the corrosion ratio of the unit area (g/(m2·a)) ;s is the surface area of the water-rock interaction (m2) ;t is the corrosion time (a).
For numerical simulation, according to Figures 16(a) and 17(b), the corrosion ratio showed a linear variation with time. The corrosion ratio (after one year of corrosion) was obtained by linear interpolation; then, the surface corrosion ratio (g/(m2·a)) was calculated with (5). The surface corrosion ratios of rock samples per unit area in the field experiment and numerical simulation are listed in Table 5.
The surface corrosion ratios of the field experiment and numerical simulation.
From Table 5, we can see that the surface corrosion ratio of rock samples in the field experiments was less than that in the numerical simulation. The surface corrosion ratio of the weak hydrodynamic condition (numerical simulation) was close to that of the field experiments. Analyzing the hydrogeological conditions of the field experiment and the setting of hydrogeological conditions in numerical simulation, the water flow in the numerical simulation was mainly in the fracture core and the flow rates were faster than those in the field experiments. Thus, in the strong hydrodynamic condition and medium hydrodynamic condition, the surface corrosion ratios obtained in the simulation were higher than those obtained in the field experiment. However, overall, the differences in the corrosion ratio between the numerical simulation and field experiment were on the same order of magnitude. Therefore, the results of the numerical simulation model were reliable, and the model can be used to predict the corrosion evolution of the fissures in soluble rock of the tunnel in the future.
6. Conclusion
In this study, the numerical simulation model of the fissures in soluble rock of karst areas of the Wulong tunnel was established, and the corrosion characteristics of the fissures were evaluated. The conclusions were as follows:
The development of the fissures in soluble rock of karst areas of the tunnel was mostly affected by the flow rates. After 100 years of corrosion, under the influence of temperature and porosity and permeability, the maximum corrosion range and corrosion ratio were only 1.2 times~1.6 times higher than the minimum values, respectively. In contrast, the corrosion range of the strong hydrodynamic condition was approximately 1.8 times~2.2 times higher than that of the weak hydrodynamic condition, and the corrosion ratio of the strong hydrodynamic condition was 2.1 times higher than that of the weak hydrodynamic condition.
In the temperature range of 10°C~35°C, the corrosion range and the corrosion ratio were proportional to the temperature because the reaction rate constant increased with temperature, causing the reactions’ equilibria to move towards the direction of the solution. After 100 years of corrosion, the corrosion ranges of 35°C were 1.2 times~1.3 times higher than those of 10°C, and the corrosion ratio of 35°C was 1.2 times higher than that of 10°C.
The porosity and permeability of the fissure zone enhanced the corrosion of the fissure and enhanced the differential corrosion of karst. The higher the porosity and permeability were, the more significant the enhancement effects on the porosity and permeability were. Many of the dead-end pores became effective pores, which greatly increased the permeability, but the slightly increased porosity also accelerated the corrosion damage of the tunnel.
The simulation results were verified by the field experiment, and the simulation model can be used to quantitatively predict the corrosion evolution of a fissure zone within the acceptable range of error.
Data Availability
The data used to support the findings of this study are available from the corresponding author upon request.
Conflicts of Interest
The authors declare no conflict of interest.
Acknowledgments
This work was supported by the Graduate Innovation Fund of Jilin University (no. 2017148), the National Natural Science Foundation of China (no. 41372334), and the China National Science and Technology Major Project (no. 2016ZX05016-005).
LiD. T.LuoY.Measurement of carbonate rocks distribution area in China19832147150LiC. S.DingJ. F.LiaoY. K.2017Beijing, ChinaThe Science Publishing CompanyGutierrezF.Alcantara-AyalaI.GoudieA. S.Hazards associated with karst2010Cambridge, UKCambridge University Press16117610.1017/CBO9780511807527.0132-s2.0-84926113977GutierrezF.PariseM.De WaeleJ.JourdeH.A review on natural and human-induced geohazards and impacts in karst2014138618810.1016/j.earscirev.2014.08.0022-s2.0-84907284971PariseM.GunnJ.2007London, UKGeological Society Special Publication 279PariseM.ClossonD.GutierrezF.StevanovicZ.Anticipating and managing engineering problems in the complex karst environment201574127823783510.1007/s12665-015-4647-52-s2.0-84949098323QinX. Q.ZhuM. Q.JiangZ. C.A review on recent advances in rocky desertification in Southwest China karst region200625234238WangY.2003Kunming, ChinaScience and Technology PressLiC. S.GaoB.MeiZ. R.2013Chengdu, ChinaSouthwest Jiaotong University PressLiS. C.ZhouZ. Q.LiL. P.XuZ. H.ZhangQ. Q.ShiS. S.Risk assessment of water inrush in karst tunnels based on attribute synthetic evaluation system201338505810.1016/j.tust.2013.05.0012-s2.0-84878795805WalthamA. C.FookesP. G.Engineering classification of karst ground conditions200336210111810.1144/1470-9236/2002-332-s2.0-0041627683ArvidsonR. S.FischerC.LuttgeA.Calcite dissolution kinetics201521541542210.1007/s10498-015-9268-92-s2.0-84938286163BuhmannD.DreybrodtW.Calcite dissolution kinetics in the system H2O·CO2·CaCO3 with participation of foreign ions1987641-28910210.1016/0009-2541(87)90154-92-s2.0-0023523333DreybrodtW.GabrovšekF.RomanovD.2005Postojna, SloveniaKarst Research Institute, ZRC SAZUEisenlohrL.MetevaK.GabrovšekF.DreybrodtW.The inhibiting action of intrinsic impurities in natural calcium carbonate minerals to their dissolution kinetics in aqueous H2O–CO2 solutions1999637-8989100110.1016/S0016-7037(98)00301-92-s2.0-0032827544GabrovšekF.HuzjanM. Ž.2000Ljubljana, SloveniaZalozba ZRCHillerT.RomanovD.KaufmannG.EptingJ.HuggenbergerP.Karstification beneath the Birs weir in Basel/Switzerland: a 3D modeling approach2012448-44918119410.1016/j.jhydrol.2012.04.0402-s2.0-84861999129JeschkeA. A.DreybrodtW.Pitfalls in the determination of empirical dissolution rate equations of minerals from experimental data and a way out: an iterative procedure to find valid rate equations, applied to Ca-carbonates and -sulphates20021923-418319410.1016/S0009-2541(02)00135-32-s2.0-0037203024LasagaA. C.2014Princeton, NJ, USAPrinceton University PressMorseJ. W.BernerR. A.Dissolution kinetics of calcium carbonate in sea water; I, A kinetic origin for the lysocline1972272984085110.2475/ajs.272.9.840PlummerL. N.BusenbergE.The solubilities of calcite, aragonite and vaterite in CO2-H2O solutions between 0 and 90°C, and an evaluation of the aqueous model for the system CaCO3-CO2-H2O19824661011104010.1016/0016-7037(82)90056-42-s2.0-0019896215PokrovskyO. S.GolubevS. V.SchottJ.Dissolution kinetics of calcite, dolomite and magnesite at 25°C and 0 to 50 atm pCO220052173-423925510.1016/j.chemgeo.2004.12.0122-s2.0-18444407473PokrovskyO. S.GolubevS. V.SchottJ.CastilloA.Calcite, dolomite and magnesite dissolution kinetics in aqueous solutions at acid to circumneutral pH, 25 to 150°C and 1 to 55 atm pCO2: new constraints on CO2 sequestration in sedimentary basins20092651-2203210.1016/j.chemgeo.2009.01.0132-s2.0-60749095643ZhaoY. J.WangF. G.LiC. S.WuF. S.CaoY. Q.Carbonate karst erosion experiment research of Guizhou Huangguoshu tunnel2016357310312CaiH.ChenZ.ZhouH.Impact analysis of fracture on the karst development and groundwater flow of Wulong cave201522216KaufmannG.Modelling karst aquifer evolution in fractured, porous rocks2016543Part B79680710.1016/j.jhydrol.2016.10.0492-s2.0-85003443888SiemersJ.DreybrodtW.Early development of karst aquifers on percolation networks of fractures in limestone199834340941910.1029/97WR032182-s2.0-0031910280DreybrodtW.Principles of early development of karst conduits under natural and man-made conditions revealed by mathematical analysis of numerical models19963292923293510.1029/96WR013322-s2.0-0029663145DreybrodtW.1988Berlin, HeidelbergSpringerSpringer Series in Physical Environment10.1007/978-3-642-83352-6DreybrodtW.The role of dissolution kinetics in the development of karst aquifers in limestone: a model simulation of karst evolution199098563965510.1086/6294312-s2.0-0025584057DreybrodtW.GabrovšekF.KlimchoukA.FordD. C.PalmerA. N.DreybrodtW.Dynamics of the evolution of a single karst conduit2000Huntsville, AL, USANational Speleological Society184193HannaR. B.RajaramH.Influence of aperture variability on dissolutional growth of fissures in karst formations199834112843285310.1029/98WR015282-s2.0-0031730874PalmerA. N.Origin and morphology of limestone caves1991103112110.1130/0016-7606(1991)103<0001:OAMOLC>2.3.CO;2AnnableW. K.2003Waterloo, ON, CanadaUniversity of WaterlooHillerT.KaufmannG.RomanovD.Karstification beneath dam-sites: from conceptual models to realistic scenarios20113983-420221110.1016/j.jhydrol.2010.12.0142-s2.0-79551493475KaufmannG.Modelling karst geomorphology on different time scales20091061-2627710.1016/j.geomorph.2008.09.0162-s2.0-61849131188KaufmannG.RomanovD.HillerT.Modeling three-dimensional karst aquifer evolution using different matrix-flow contributions20103883-424125010.1016/j.jhydrol.2010.05.0012-s2.0-77953650526LiuQ.LuY. R.ZhangF. E.QiJ. X.ZhangS.Qualitative analysis of microcorrosion of limestone induced by temperature and hydrodynamic pressure201031S2149154YuH. D.YangD.WangD.HanM. Y.Top-down fabrication of calcite nanoshoot arrays by crystal dissolution201022293181318410.1002/adma.2010005032-s2.0-7995203410720535738ZhuW. H.QuX. Y.QiuL. W.ChenY.GongF. X.WuS. X.Characteristics and erosion mechanism of carbonate in acetic acid and hydrochloride solutions, an example from the Nanpu depression201534619625SvenssonU.DreybrodtW.Dissolution kinetics of natural calcite minerals in CO2-water systems approaching calcite equilibrium19921001-212914510.1016/0009-2541(92)90106-F2-s2.0-0027041124CubillasP.KöhlerS.PrietoM.CausserandC.OelkersE. H.How do mineral coatings affect dissolution rates? An experimental study of coupled CaCO3 dissolution—CdCO3 precipitation200569235459547610.1016/j.gca.2005.07.0162-s2.0-29544438004FinneranD. W.MorseJ. W.Calcite dissolution kinetics in saline waters20092681-213714610.1016/j.chemgeo.2009.08.0062-s2.0-70349400732FischerC.KurganskayaI.SchäferT.LüttgeA.Variability of crystal surface reactivity: what do we know?20144313215710.1016/j.apgeochem.2014.02.0022-s2.0-84896520215BusenbergE.Niel PlummerL.Kinetic and thermodynamic factors controlling the distribution of SO32− and Na+ in calcites and selected aragonites198549371372510.1016/0016-7037(85)90166-82-s2.0-0021890370GledhillD. K.MorseJ. W.Calcite dissolution kinetics in Na–Ca–Mg–Cl brines200670235802581310.1016/j.gca.2006.03.0242-s2.0-33750891822GongQ.DengJ.WangQ.YangL.SheM.Experimental determination of calcite dissolution rates and equilibrium concentrations in deionized water approaching calcite equilibrium201021440241110.1007/s12583-010-0103-32-s2.0-77955140187LasagaA. C.LüttgeA.Mineralogical approaches to fundamental crystal dissolution kinetics200489452754010.2138/am-2004-04072-s2.0-1942467061LüttgeA.Crystal dissolution kinetics and Gibbs free energy20061502-324825910.1016/j.elspec.2005.06.0072-s2.0-28444454825GongQ.DengJ.WangQ.YangL.SheM.Calcite dissolution in deionized water from 50°C to 250°C at 10 MPa: rate equation and reaction order2008825994100110.1111/j.1755-6724.2008.tb00655.xShirakiR.BrantleyS. L.Kinetics of near-equilibrium calcite precipitation at 100°C: an evaluation of elementary reaction-based and affinity-based rate laws19955981457147110.1016/0016-7037(95)00055-52-s2.0-0028794374AdamczykK.Premont-SchwarzM.PinesD.PinesE.NibberingE. T. J.Real-time observation of carbonic acid formation in aqueous solution200932659601690169410.1126/science.11800602-s2.0-7294909420819965381XuJ.FanC.TengH. H.Calcite dissolution kinetics in view of Gibbs free energy, dislocation density, and pCO22012322-323111810.1016/j.chemgeo.2012.04.0192-s2.0-84863494225ZhouW.BeckB. F.BeynenP. E.vanEngineering issues on karst2011Dordrecht, NetherlandsSpringer94510.1007/978-94-007-1207-2_2ZhangX. H.LiuQ. W.Analysis of the features of water outbursts from underground rivers and the countermeasures in Wulong tunnel2005425964Amann HildenbrandA.DietrichsJ. P.KroossB. M.Effective gas permeability of tight gas sandstones as a function of capillary pressure – a non-steady-state approach201616338310.1111/gfl.121552-s2.0-84952838389AngelierJ.BergeratF.DauteuilO.VilleminT.Effective tension-shear relationships in extensional fissure swarms, axial rift zone of Northeastern Iceland199719567368510.1016/S0191-8141(96)00106-X2-s2.0-0031147779FaulknerD. R.JacksonC. A. L.LunnR. J.SchlischeR. W.ShiptonZ. K.WibberleyC. A. J.WithjackM. O.A review of recent developments concerning the structure, mechanics and fluid flow properties of fault zones201032111557157510.1016/j.jsg.2010.06.0092-s2.0-78649764753GiwelliA.Delle PianeC.EstebanL.ClennellM. B.DautriatJ.RaimonJ.KagerS.KiewietL.Laboratory observations of fault transmissibility alteration in carbonate rock during direct shearing201616467210.1111/gfl.121832-s2.0-85006324385JohriM.DunhamE. M.ZobackM. D.FangZ.Predicting fault damage zones by modeling dynamic rupture propagation and comparison with field observations201411921251127210.1002/2013JB0103352-s2.0-84899734091MicarelliL.BenedictoA.WibberleyC. A. J.Structural evolution and permeability of normal fault zones in highly porous carbonate rocks20062871214122710.1016/j.jsg.2006.03.0362-s2.0-33745005539HeJ.KangY.LiuD.LeiM.ShuZ.The stress sensitivity research on porous and fractured porous carbonate reservoirs20052828486LiD. Q.KangY. L.YouL. J.Experimental study on permeability stress sensitivity of carbonate rocks201420231222WanY.ZanJ.TaoH.Study on pore structure of carbonate reservoir200827121314ZhuD.LiQ.Predication of outflow in tunnel based on experimental permeating coefficient of precipitation19951100102XuT.SonnenthalE.SpycherN.PruessK.TOUGHREACT—a simulation program for non-isothermal multiphase reactive geochemical transport in variably saturated geologic media: applications to geothermal injectivity and CO2 geological sequestration200632214516510.1016/j.cageo.2005.06.0142-s2.0-30944461755YangL.XuT.WeiM.FengG.WangF.WangK.Dissolution of arkose in dilute acetic acid solution under conditions relevant to burial diagenesis201554657310.1016/j.apgeochem.2015.01.0072-s2.0-84923931528HuangS. Y.The corrosion of carbonates and environment temperature19874287296