Abstract

Underground engineering, especially deep geotechnical engineering, is often affected alternately by groundwater infiltration and ventilation drying during construction and use. In addition, acid rain, mineral dissolution, acid deposition, and other factors make groundwater acidic. This caused the underground structure to erode and has also threatened its safety and durability. In this study, both the physical and mechanical properties under acid dry-wet (A-D-W) cycles were investigated. Deep sandstone was treated repeatedly under acid cycling, and its physical parameters were measured. Uniaxial compressive strength tests and microtests were carried out. Finally, using a combination of scanning electron microscopy (SEM) and backscatter electron images (BSE), the microstructural changes of sandstone under the combined action of an acidic environment and a dry-wet cycle were described. The test results show that the mass, P-wave velocity, peak stress, and elastic modulus of sandstone after the A-D-W cycles decreased by 0.43%, 7.87%, 70.20%, and 88.10%, respectively. With the increase in the number of cycles, the loss of these indicators increased. However, the peak strain increased with an average increase of 55.8%. In addition, 10 cycles are the critical point for physical and mechanical indicators of the A-D-W cycles. After 10 cycles, the changes of various indicators increase rapidly. Microscopic analysis shows that the reasons of this phenomenon were that the specimen was corroded by the sulfuric acid solution, which resulted in the development of pores and cracks and the decline of physical and mechanical properties. Conversely, the development of pores was hindered by sediment, which slowed down the decline rate of the physical and mechanical indexes of specimens, but this phenomenon disappeared with the increase of A-D-W cycles. Based on the analysis of the experimental phenomena, the constitutive model of uniaxial compression based on Weibull damage variable has been established in this article, and the model has the best effect in acid solution with less cycle times or pH value of solution greater than 6.

1. Introduction

Water-rock interaction is an important driving force for the evolution of the near-surface environment. It has a great influence on the geological structure of the strata and the stability of underground engineering under the coupling action of chemistry and mechanics [14]. As one of the important forms of water-rock interaction, dry-wet cycling studies the effects of water on the physical and mechanical properties of rocks by simulating periodic changes of water content. In recent years, sufficient research works have been done on the variation of static rock mechanical parameters under dry-wet cycling [58]. The results show that the dry-wet cycling has an obvious deterioration effect on the basic mechanical properties of rock, such as compressive strength, elastic modulus, cohesion, and the internal friction angle. Currently, the study of the variation of mechanical parameters of rock materials under a dynamic environment is becoming more and more abundant, Zhou et al. [9] found that the dynamic compressive strength and elastic modulus of sandstone decreased with the increase of the number of dry-wet cycles, and the empirical equation of the effect of the strain rate and dry-wet cycles on the dynamic compressive strength of the rock materials was established through experiments. Du et al. [10, 11] found that the dynamic compressive strength of red sandstone was positively correlated with the loading rate and negatively correlated with the number of dry-wet cycles.

The attenuation function model was established to evaluate the dynamic compressive strength of red sandstone under the influence of dry-wet cycles, whereas the constitutive model was developed for the coupling condition of impact load and dry-wet cycles. With the advancement of research, the testing range of mechanical parameters extends to tensile strength [12], shear properties [13, 14], and fracture characteristics [15]. With the development of microtechnology and the leap of detection technology, scanning electron microscopy (SEM) [16], nuclear magnetic resonance (NMR) [17], computed tomography (CT) [17], and acoustic emission (AE) [18] technologies are now often used to study the microchanges of rocks under dry-wet cycles.

In addition, due to acid rain, mineral dissolution, acid deposition, and other factors, the groundwater is acidic. The investigation of the rock damage law under the coupling action of the chemical solution and the dry-wet cycles is still in its initial stage. Because of the difference in chemical solutions, the experimental results will be slightly different, but the general law is consistent. Li et al. [19, 20] studied the variation of calcareous sandstone in an acid solution, and the chemical reaction rate was used to describe the degree of rock corrosion. The corresponding chemical damage model was established. Liu et al. [2126] took the sandstone in the ebb-flow zone of a reservoir in the Three Gorges as the research object. They summarized the changing law of mechanical properties and macro- and microdamage mechanisms of sandstone under dry-wet cycles [21, 22] as well as the change of the physical parameters and microdamage mechanisms of limestone under chemical solution immersion. They found that the basic reason for the decline of the mechanical properties of rocks was that the pH value of the chemical solution changed the porosity and microdamage evolution of the rocks [23]. In addition, they also analyzed the deterioration rule of sandstone mechanical parameters under acid-base conditions [24, 25] and developed a numerical simulation program based on an improved Duncan–Chang constitutive model [26]. Matsumoto et al. [27] found that the dissolution of some substances in rocks can promote the development of cracks and that physical and chemical weathering promote each other and jointly accelerate rock failure. Sun and Zhang [28] used a magnesium sulfate solution to treat sandstone specimens by dry-wet cycles.

The results show that 30 cycles is the threshold of the physical and mechanical parameters of sandstone. They pointed out that color lightness and thermal conductivity can effectively evaluate the tensile strength of sandstone. Li et al. [29] investigated the effects of acid solution immersion on the physical and mechanical properties of sandstone. It was found that the deformation characteristics of sandstone under immersion in an acid solution such as sulfuric acid were softening, and the softening degree was positively correlated with the immersion time. Based on the experimental results, the damage constitutive model of sandstone under an acid ring was established. In addition, Fang et al. [30] studied the failure mechanism of the rock mass under the action of dry-wet cycling from the micromorphological change of joint surface. His research found that the dry-wet cycling made the roughness coefficient of the rock surface, the compressive strength of the face wall, and the basic friction angle gradually decrease, which led to a decrease in the shear strength of the joint surface. Yang et al. [31] studied the influence of water-rock cycles on the macro- and microproperties of rocks.

Macroscopically, the uniaxial mechanical properties of the water-rock cycle decrease progressively and continuously. Microscopically, the grain size characteristics and pore distribution of the rock surface gradually increase. To simulate more extreme environments, Ni et al. [32] carried out dry-wet-freezing-weathering simulation tests on the distribution of two sandstones. It was found that water content and salt solution treatment have a great influence on rock properties. Although salt solution can alleviate rock weathering, the specific mechanism needs to be studied. Zhang et al. [33] treated expansive rocks by dry-wet cycling and studied the law of crack propagation. The results show that the expansibility of expansive rocks decreases with the dry-wet cycling, which shows that the degree of cracking, fractal dimension, and the strength of expansive rocks decrease. Hu et al. [34] studied the crack development mechanism of mudstone under dry-wet cycling. Expansion during drying resulted in the increase of crack width and the decrease of crack size due to the entry of water. However, the penetration area of water increased with the number of dry-wet cycles, which led to the rapid increase in crack number and accelerated the fracture process of the mudstone. All of the studies show that the repeated processes of dry expansion and wetting compaction of different rock masses are characterized by strength reduction and crack increase. However, the different water solutions and the change in wetting and drying cycles lead to a change of rock properties, and therefore, the wetting and drying cycles under specific conditions are still worthy of further study.

The construction of underground engineering in a coal mine is an example of typical deep underground engineering. The surrounding rock of the coal mine roadway is under high geostress for a long time. When disturbed by excavation, it usually produces expansion or deformation, which is in the state of dry-wet circulation under the combined action of roadway ventilation and acid groundwater infiltration. In this article, the sandstone with a buried depth of −965 m is processed into specimens. After soaking in a sulfuric acid solution with pH values of 4 to 7 for five days, the sandstone is dried and wetted 0, 5, 10, and 15 times, respectively, and uniaxial compressive strength tests and microtests are then carried out to study the change rules of the physical, mechanical, and the evolution process of the microstructure. Based on the test, the uniaxial damage model of the sandstones were established, and its applicability and accuracy were evaluated.

2. Materials and Methods

2.1. Specimen Preparation

Sandstone specimens were collected from a roadway with a depth of −965 m in the Huainan mining area, Anhui Province, China (see Figure 1(a)). The roadway has been affected by groundwater infiltration for a long time. The surrounding rock of the roadway roof has expanded, and the U-shaped steel support has deformed (see Figure 1(b)). The apparent characteristics of the rock specimens are grey-white, with a few dark brown spots (see Figure 1(c)).

According to the results of the X-ray diffraction (XRD) spectrum (see Figure 2(a)) and energy dispersive spectrometer (EDS) analysis, the main components of the specimens are quartz, albite, calcite, kaolinite, and iron compounds. The EDAX Octane plus spectrometer (see Figure 3(c)) was used to analyze the two detection points on the specimen surface (see Figure 2(b)).

The types of material elements on the surface of the specimens (see Figures 2(c) and 2(d)) and the results of the mass ratio analysis are shown in Table 1.

Combined with Figures 2(b)2(d), it can be seen that the surface of the specimens is relatively flat with a small number of white spots. Spots 1 and 2 were selected to analyze the types of elements and their mass ratios, and the relative number of elements was calculated according to the relative atomic mass (see Table 1). It can be seen that the brighter the color, the more the types of elements, and the number of elements and relative atoms are in accord with the XRD results. It is proven that the matrix is mainly quartz (SiO2). The speckles are metal element compounds and cementing materials (such as Fe2O3 and MgO). The metal element compounds are mainly Fe-containing compounds, and dark red spots can be observed on the surface of the macroscopic specimens.

In order to better understand the mineral composition, particle size, and structure of the rocks, specimens were taken and the thin sections were identified. Thin sections of rock with length × width × thickness of 30 mm × 20 mm × 0.03 mm were prepared by grinding. The thin sections were observed under a Leitz-ORTHOPLAN polarizing microscope in Germany (see Figure 4). The results show that the quartz content is 70%, the clay mineral is 15%, the siliceous rock is 10%, and the feldspar is 5%. The scattered distribution of quartz can be observed on the thin slices. The filler consists of muddy impurities and a small amount of biotite and kaolinite.

Specimens were selected based on similar P-wave velocities and made into cylindrical specimens with a diameter of 50 mm and height of 100 mm according to the ISRM (International Society for Rock Mechanics) standard [35]. The surface smoothness was controlled to within ±0.05 mm and the vertical deviation between upper and lower surfaces was controlled to within ±0.25 (see Figure 5).

The basic physical parameters of the prepared specimens were measured and the results are shown in Table 2.

2.2. Test Equipment

The mass of the sandstone specimens was measured by an electronic balance with an accuracy of 0.01 g and the P-wave velocity was measured by a nonmetallic ultrasonic testing analyzer of the NN-4B type. Drying and wetting cycles were dried in a 101-00 constant-temperature drying oven and immersed in a glass box. The uniaxial compression test used a WAW-2000 microcomputer-controlled electro-hydraulic servo universal testing machine with a maximum range of 2000 kN. The testing device is shown in Figure 3(a). In this experiment, the XRD (X-ray diffraction) instrument was a SmartLab high-resolution X-ray diffraction instrument (see Figure 3(b)), and the SEM (scanning electron microscope), BSE (back-scattered electron Imaging), and EDS (energy-dispersive spectrometer) were done using a Hitachi S-3400N scanning electron microscope (see Figure 3(c)).

2.3. Test Method
2.3.1. Solution Preparation

Solutions with pH of 4, 5, and 6 were prepared with concentrated H2SO4 and the solution with pH = 7 was prepared with tap water. Each solution measured 12 L.

2.3.2. The A-D-W Cycle Mechanism

The grouping method was as follows: the 48 specimens were divided into four groups, which were soaked in the solutions with a pH of 4, 5, 6, and 7, respectively. Each group was further divided into four additional groups, which were then treated in four stages: stage I (5-day solution soaking), stage II (0–5 cycles), stage III (5–10 cycles), and stage IV (10–15 cycles).

According to the groupings, the specimens were treated by soaking them in different pH solutions for 12 hours, and then drying them in the drying oven for 12 hours, which was regarded as a dry-wet cycle. The first drying temperature of the drying oven was set to 105°C, and then each temperature was maintained at 60°C. The schematic diagram of the A-D-W cycle process is shown in Figure 6.

2.3.3. Uniaxial Compressive Strength Test Method

In the uniaxial test, the displacement control method was used to load the specimens, and the loading speed was 2 mm/min. To reduce the influence of the surface smoothness of the specimens on the test and to ensure the stability of the test instrument, all specimens were preloaded at the beginning of the test, and the preload was 500 N.

2.3.4. Microscopic Test Method

In order to compare the microstructural changes of the specimens experiencing different A-D-W cycles, four groups of specimens in stages I-IV were selected for microtesting at pH = 4. The microevolution process of the specimens was observed by SEM and BSE images.

3. Results and Discussion

3.1. Variation in Mass

During the coupling process of the acid solution and the wet-dry cycle, the acid solution immersion will dissolve some substances in the sandstone or react with the acid solution. On the other hand, the specimen surface may be peeled off repeatedly due to wet-dry alternation, resulting in a decrease of specimen mass. Therefore, the variation of the sandstone mass is an important index of the influence of the A-D-W cycle on sandstone damage. Table 3 records the average mass change of each group of sandstone specimens after completion of different A-D-W cycles.

In Figure 7, the curve corresponds to the left coordinate axis and the column diagram corresponds to the right coordinate axis. The changes in total mass loss and stage mass loss of sandstone specimens in different pH solutions are shown It can be seen from the variation of the curve that when pH = 4, the mass loss gradually increases to 0.17%, 0.21%, 0.31%, and 0.43% with the progress of the cycle. By comparing different curves, it can be seen that the mass loss increases with the decrease of pH value. Taking the test group with pH = 4 as an example, the proportions of the mass loss in each stage to total loss are 40.4%, 22.9%, 8.7%, and 28.0%, respectively. It can be seen that the mass loss increases obviously in Stage IV. This is because the sulfuric acid solution and the calcium ions dissolved in the specimen will form gypsum (CaSO4), which has a certain inhibitory effect on the crack development of sandstone specimens with higher integrity. After 10 cycles, the crack development is complete, the inhibitory effect is weakened, so the internal material of the specimen continues to peel off and the decline in quality increases again.

The pH value of the acidic solution has a significant effect on the total mass loss of specimens. With the decrease in pH, the mass loss values are 1.05 g, 1.38 g, 1.74 g, and 2.18 g, respectively. It can be observed that the mass loss of the specimens increases with the increase in acid concentration at the same wet-dry cycle stage, and the difference becomes more obvious with the increase in cycle times.

3.2. Variation of P-Wave Velocity

P-wave velocity is an important index for measuring the compactness of a material. The higher the velocity, the higher the compactness of the sandstone. Ultrasound travels faster in solids than in air, and gas-solid coupling leads to low average velocity. This indirectly indicates the cracks and pores in the specimen. The rock integrity coefficient k is introduced to describe the degree of development of rock fractures, and the k value is expressed as shown in where Vm is the initial P-wave velocity of the specimen, m/s and Vr is the P-wave velocity of the specimen after different dry-wet cycles, m/s. The calculation results are shown in Table 4.

The smaller the k value in Table 4 indicates the poorer the rock integrity. It can be seen that with the increase of the number of dry-wet cycles, the k value of the specimen decreases gradually. At the same time, with the decrease of the pH value, the integrity of the sandstone specimen also decreases. The specific decline rule is shown in Figure 8.

In Figure 8, the curve corresponds to the left coordinate axis and the column diagram corresponds to the right coordinate axis. The variation of the total loss of P-wave velocity and stage loss of P-wave velocity of sandstone specimens in different pH solutions are shown.

It can be seen from the curve that the P-wave velocity of each group of specimens shows a downward trend with the progress of the A-D-W cycle. The most obvious decrease is at pH = 4, where the total decrease in P-wave velocity at each stage is 2.20%, 3.14%, 5.16%, and 7.87%, respectively. In addition, it can be seen from the column diagram that the P-wave velocity losses at stages I and II are equivalent, which is different from the law of continuous decline in mass loss analysis. This is because the ultrasonic wave is mainly emitted along the axis of the specimen during the P-wave velocity measurement, and the main measuring range is close to the center of the specimen. The spalling damage at the edge of the surface of the specimen has little effect on its interior. Also, it can be seen that the loss at stage IV is significantly higher than that of the first three stages, which agrees with the observation that the loss at stage IV increases sharply after 10 cycles found in Section 3.1.

3.3. Stress-Strain Curve

From the test results, the stress-strain curves of the treated and untreated specimens are drawn. The test curves of the treated specimens at stages I, II, III, and IV under different pH values are shown in Figure 9(a–d).

The typical σ − ε curves of the A-D-W cycle for each group of specimens show the following. (1) The peak stress of specimens experiencing the A-D-W cycle is significantly lower than that of untreated specimens, and with the increase of cycling stage, the decrease of stress increases gradually. (2) The slope of the experiencing curve, the elastic modulus, also decreases with the increase of the cyclic stage. (3) Figures 9(a)9(d) show that with the increase of the acid concentration, the peak strain increases gradually, and the more the cyclic stages are, the more obvious the phenomenon is.

According to the σ − ε curve, it is easy to get the σ (peak stress) and ε (peak strain) of sandstone under different states. The E (elastic modulus) of sandstone under different states can be calculated by fitting the slope of the curve in the straight line before the curve reaches the peak stress. Its data are recorded as shown in Table 5.

3.3.1. Peak Stress

The results shown in Figure 10(a) show that compared with the peak stress of 97.96 MPa in the control group, after 15 times of A-D-W cycles, the percent decrease in stress are 53.2%, 64.40%, 70.3%, and 70.2% with the increase in acid concentration. At the same pH value, pH = 5, the peak stress loss for stages I to IV at different cyclic stages can be found to be 29.68 MPa, 28.02 MPa, 2.89 MPa, and 8.30 MPa, respectively. The corresponding stress loss of specimens is 30.46 MPa, 19.20 MPa, 3.54 MPa, and 10.91 MPa at pH = 6. Two sets of data show that the loss of peak stress in each stage decreases gradually in the first three stages of the cycles and increases again in Stage IV.

3.3.2. Elastic Modulus

The results shown in Figure 10(b) show that as compared with the control group, the elastic modulus is 20.75 GPa. After 15 times of A-D-W cycles, the decreasing rates of the test group were 44.9%, 80.6%, 84.5%, and 88.1% for pH = 7, 6, 5, 4, respectively. The percent loss of the pH group is more than 80%, which is 1.80–1.96 times that of the neutral group. This shows that acidity has a significant effect on the elastic modulus of the specimens. When the pH value is the same, the elastic modulus of the specimens decreases continuously with the progression of cycles.

3.3.3. Peak Strain

The peak strain values of specimens with different pH groups at different cyclic stages are shown in Figure 11. In neutral solution, the peak strain of specimens decreases very gradually with the progression of cycles, which indicates that the brittleness of the specimens increases under normal dry-wet cycling conditions, and it is easier to produce sudden cracks leading to failure. However, with the addition of acidic solution, the peak stress of each group of specimens increases, which indicates that the acidic solution can soften the specimens and increase the ductility slightly, resulting in greater deformation and relatively mild failure mode. The average strain values of the different pH groups were also plotted. The results showed that the peak strain growth trend gradually slowed down from Stage I to Stage III, but by Stage IV, the total increase rate was 55.8%. The change rule also indicated that 10 cycles were the turning point for the physical and mechanical properties of the A-D-W cycle.

3.4. Microstructure Evolution

Specimens from stages I–IV of the pH = 4 groups were collected. Scanning electron microscopy (SEM) was used to observe and obtain the microscopic morphology and structure of the specimens. The images obtained are shown in Figure 12. The SEM and BSE images of the specimens in the natural state are shown in Figure 13 as a control group for comparison.

In Figure 13, it can be observed that the surfaces of the sandstone specimens in the natural state are relatively flat and compact, there are no obvious cracks, and there are a few scattered natural pores; at this time, the peak stress and elastic modulus of the specimens are the maximum values of each group, and their integrity and mechanical properties are good.

The main reason for the decrease of the mechanical parameters of the specimen during the dry-wet cycle is the development of pores and cracks, which will gradually reduce the strength and elastic modulus of the specimen. According to common sense, the corrosion of the acid solution will cause the development of cracks, and the expansion of the cracks in turn allows the acid solution to further enter the inside of the specimen. Under the repeated action of the dry-wet cycle, these two factors promote each other and accelerate the damage rate of the specimen. However, the test results show that the rate of decrease of the mechanical parameters of the specimen gradually decreases before 10 dry-wet cycles, and after 10 times, the specimen begins to deteriorate rapidly. It is found through the microstructure photography that the cause of this phenomenon is the development of pores or cracks affected. In the early stage of the dry-wet cycle, the surface of the specimen is dense, but the number of pores increased compared with Figure 13(a). The pores served as channels, which causes the acid solution to further corrode the inside of the specimen, causing the inside of the specimen to be damaged. The mechanical properties of the components are significantly reduced (see Figure 13(a)). As the cycle progresses (see Figure 12(b)), microcracks begin to develop, and the crystalline material gradually increases. When the specimen is under pressure, the cracks start to grow from the microcracks and promote the failure of the specimen. Mechanical parameters continue to decline. However, in Figure 12(c), we find that new precipitation occurs, which adheres to the surface of the specimen to a certain extent, hindering the entry of the acidic solution, forcing the acidic solution to react with the precipitated material, and then further accelerate the crack development. As a result, the rate of decline of the macro-mechanical parameters of the specimen is slowed down, but the amount of precipitation is limited after all. As the dry-wet cycle progresses, this phenomenon gradually disappears. After 10 dry-wet cycles, cracks are distributed in the network and connected. Crystal materials begin to gradually decrease (see Figure 12(d)), the macro-mechanical properties of the sample further deteriorated, and the rate of decline of the mechanical parameters of the specimen increased again.

In order to further observe the process of crack development, the surface of the specimen was ground and flattened, and a back-scattering electron (BSE) image of 200 μm was taken (Figure 14).

Figure 14(a) shows that the surface of the specimens is compact, the pores are small, most of them are in the cementing materials, and the matrix is relatively complete. With the enlargement of the pore area, new fine pores continue to appear on types of cement (see Figure 14(b)). Figure 14(c) shows that fissures first appear at the interface between the matrix and the cement, and the matrix in the original pore begins to curl in layers under the corrosion of the acidic substances. The cracks further develop and prolong, and microcracks begin to appear in the matrix. The matrix curling range in the pore is enlarged and the degree of curling is deepened, which results in the decrease of the strength of specimens at the macroscopic level (see Figure 14(d)).

The coal mine sandstone used in this test is mainly composed of quartz, albite, kaolinite, and calcite. According to the change of microstructure, the XRD diffraction patterns (see Figure 15) of the specimens before and after the test show that the main substances on the surface of the specimens are still quartz and feldspar, while kaolinite and calcite have disappeared, and new products, calcium silicate and gypsum, have appeared. It is inferred that the main ion reactions are as in Figure 15:

Quartz is almost insoluble in an acidic environment, so acidic solution has little influence on it. Albite (NaAlSi3O8) reacts with H+ in an acidic environment to form Na+, Al3+, and the precipitate orthosilicic acid (H2SiO4). The specific reaction equation is as follows:

The specific reaction equation of kaolinite (Al2Si2O5(OH)4) with acid is given by:and the specific reaction equation of calcite (CaCO3) with acid is as follows:

Reaction of metal compounds with acid solution is as follows:

Fe3+ causes the solution to be reddish brown after immersion.

In addition, the free Ca2+ ion binds with the ion to form gypsum (CaSO4), and free Ca2+ ion binds with the ion to form calcium silicate (CaSiO3), which are attached to the surface of the specimen with the original silicic acid (H2SiO4) and the silicic acid (H2SiO3). In the early stages of the A-D-W cycles, these substances can hinder the development of pores and slow down the damage rate of the physical and mechanical parameters. However, the structure and strength of the new substances are greatly different from those of sandstone, which ultimately leads to the failure of the specimen:

3.5. Damage Mechanism Analysis

The failure mechanism of sandstone in the acidic environment can be summarized by combining the microstructure evolution process of the specimens with the failure morphology of the macrospecimens as shown in Figure 16.

The effects of the acidic solutions on sandstone include physical, such as lubrication, softening, mudding, and strengthening of bound water; chemical, which includes ion exchange, dissolution, hydration, hydrolysis, redox, precipitation; and mechanical, including pore hydrostatic pressure and pore hydrodynamic pressure, effects.

3.5.1. Water Lubricates Sandstone

Lubrication occurs when water enters the structural planes such as cracks, joints, and fault planes of the sandstone. It reduces the friction resistance and enhances the shear stress effect of the structural plane. This results in the shear failure of the specimens after multiple dry-wet cycles, which is reflected in the reduction of the internal friction angle in the macrosense.

3.5.2. Chemical Damage

The influence of the chemical solution of H2SO4 as the main acidic substance on sandstone includes two aspects. The acidic solution reacts with mineral particles, which reduces the size of the mineral particles and changes the size and shape of the rock particles. This results in changes in the size and shape of the pores and cracks (see Figure 16(a)). But, as this happens, the sulfuric acid solution changes the mineral composition of the rocks. Gypsum (CaSO4) formed by the combination of and Ca2+ in solution, albite (NaAlSi3O8), and kaolinite (Al2Si2O5(OH)4) form orthosilicic acid and silicic acid in the acidic environment. These products can adhere to the pore of the specimen when the crack is small, and change the porosity, which prevents further reaction between the acid solution and the internal sandstone to a certain extent (see Figure 16(b)). As the number of cycles increases, the mixture dissolves and breaks down (see Figure 16(c)). Cracks develop further, and the acid solution promotes the development of cracks, which lead to the failure of the specimens (Figure 16(d)).

3.5.3. Effect of Pore Water Pressure

Because of the acid dry-wet cycling process, the pore passages in the specimens are opened and the water content of the rock mass rises. Under the high geostress environment, the pore water in the specimen cannot be eliminated. Because of the water wedge effect, it promotes the propagation of the crack tip, thus roughly destroying the rock.

4. Establishment of Uniaxial Damage Model

4.1. Constitutive Equation of Sandstone under Uniaxial Compression

Generally, sandstone is a brittle material, which has three stages in uniaxial compression, i.e., compaction stage, linear elastic stage, and failure stage. Because of the complexity of the change law of the compaction stage and the sudden change of the failure stage, it is difficult to summarize its change law. Therefore, the uniaxial constitutive equation of sandstone is mainly the linear elastic stage. Its constitutive equation is shown as follows:where σ is uniaxial compression stress value, MPa; E is modulus of elasticity, GPa; and ε is uniaxial strain value, 10−2.

4.2. Damage Constitutive Equation of Sandstone under Uniaxial Compression in the A-D-W Cycles
4.2.1. Determination of Weibull Damage Variable

Weibull distribution [36] is used to describe the strength change of sandstone under the A-D-W cycles, and its probability density function is as follows:

The cumulative probability of distribution is as follows:where m and F represent the inherent properties of rock materials [37].

From the microscopic point of view, the damage of rock material is caused by the damage of microunit structure in rock. It is assumed that for the specimens after n times of cycle, the number of microunits damaged under a certain load is c, the total number of microunits is N, and the statistical damage variable is defined as D.

During the [ε, ε + dε], the number of broken microunits is Np(x), and then,

Joint (10) and (11) calculations.

When D = 0, it means that the rock is undamaged; when D = 1, it means that all microunits are invalid; the size of D value indicates the degree of damage in the rock, D ∈ [0, 1].

The uniaxial damage constitutive model of sandstone with damage variable is shown as follows:

4.2.2. Calculation of Weibull Damage Variable Parameters

The damage constitutive model parameters F and m of rock under uniaxial compression can be calculated by the coordinate C(εc, σc) of the peak point on the uniaxial stress-strain curve. When the stress-strain curve reaches the peak value, the slope of the curve is 0, and equation (14) can be obtained.

The coordinate of point C satisfies equation (13) and can be substituted into equation (14).

Similarly for equations (10) and (11),

4.2.3. Determination of Damage Variables of A-D-W Cycle

The damage variable of acid dry wet cycle under uniaxial compression is defined as ,where En is the modulus of elasticity of n cycles under uniaxial compression, MPa; E0 is the elastic modulus of 0 cycle under uniaxial compression, MPa.

According to the mechanical test results (see Table 5), it is assumed that the change of elastic modulus E with the number of A-D-W cycles n is continuous. Introducing fitting function En = E0[1 − Kln(na + 1)], the change rule of elastic modulus under different pH values is fitted, and the results are shown in Figure 17.

4.2.4. Uniaxial Damage Constitutive Model of Sandstone under A-D-W Cycles

The uniaxial damage constitutive model of sandstone under different acid concentration is shown as follows:Where σ is the strength of sandstone, MPa; ε is the strain of sandstone; E0 is the elastic modulus of sandstone in its natural state, GPa; and K and a are the correlation coefficients of pH value. According to the fitting results in Table 6, specific figures of this value can be obtained under different pH values; m and F are inherent properties of rock materials.

4.3. Verification of Sandstone Uniaxial Damage Constitutive Model

The inherent properties of rock are changed due to the specimen has been immersed in different solutions for 5 days (stage I). Therefore, the values of m and F of the soaked specimen are calculated according to formulas (16) and (17), as shown in Table 6. According to the constitutive model (19), the uniaxial mechanical parameters of sandstone are calculated and compared with the test values. The comparison results are shown in Figure 17.

The test value of stage I is calculated according to m = 1.177, F = 1.268 × 10−3 (before immersion). Because immersion changes the inherent properties of the specimens, when in stages II, III, and IV, the test values and calculation values of each group of specimens are as shown in Figure 18, according to the m, F values corresponding to different pH values.

The predicted value of damage model and test value are shown in Table 7.

4.4. Evaluation of Sandstone Uniaxial Damage Constitutive Model

As shown in Figure 18, when pH = 7, the predicted value of peak stress is exact. However, when the pH value is lower than 7 and the number of dry-wet cycles is more than 10, the predicted value of the specimen is lower than the test value. Therefore, the correlation coefficient R is used to evaluate the calculated value curve and the theoretical value curve, and the calculation formula is shown in equation (20) [38].Where x and y, respectively, represent two groups of different variable values.

According to equation (20), calculate the correlation coefficient R between the peak stress test value and the calculated value of the specimen at different pH values, as shown in Table 8.

According to the R value of each experimental group in Table 8, the correlation coefficient of each group of specimens is greater than 0.8, showing a high correlation, indicating that the model prediction is accurate. Due to the cement produced on the surface of the specimen in the early stage of sulfuric acid solution (see Figure 14), further deterioration of the specimen is hindered to some extent, resulting in slightly higher test value than the predicted value.

According to equation (21), the relative error between the test value and the calculated value of each group of specimens is calculated.Where W is the relative error, σ is the test value and σp is the calculated value. The calculation results are shown in Table 9.

According to the calculation results in the table, the larger the W value, the larger the deviation of the predicted value. When the pH value is greater than 6, the W value can be controlled below 0.2, which is more applicable. At the same time, when the number of dry-wet cycles is less than 5, the W value can be controlled within 0.08, which is more applicable. The specific application scope is shown in Figure 19.

Therefore, the scope of application of the model can be summarized as the case that the pH value is more than 6 in the acid solution or the number of dry-wet cycles is less than 5.

5. Conclusions

(1)In the acidic environment, with the increase in the number of cycles, the quality of the sandstone and the P-wave velocity both decrease, and the amount decreases to 0.43% and 7.87%, respectively. At the same time, with the increase of the acidic concentration, the loss of both increases. In the first three stages of the cycle, the loss of the mass and the P-wave velocity in each stage decreases gradually and then increases rapidly after 10 cycles.(2)In an acidic environment, the peak stress and modulus of elasticity of sandstone decrease with the increase in the number of cycles, with decrease rates of 70.2% and 88.1%, respectively, whereas the peak strain increases to 55.8%, showing a softening phenomenon. The change rate of the three mechanical indicators in each stage has 10 cycles as the critical point and the change rate in Stage IV increases sharply.(3)Microscopic analysis shows that on the one hand, acid solution will corrode rock specimens, which increases the porosity of the specimen. Microcracks begin to develop at the interface of different substances, resulting in curling of the rock particles, which decreases the physical and mechanical properties of sandstone. On the other hand, the sulfuric acid reacts with some substances inside the specimen to form precipitation, which hinders the development of pores. The phenomenon disappears gradually with the increase of the number of cycles.(4)The damage model of uniaxial compression based on Weibull distribution function is established, and compared with the test results, and the model has the best effect in acid solution with less cycle times or pH value of solution greater than 6.

Data Availability

The datasets generated and analyzed during the current study are available from the corresponding author on reasonable request.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Authors’ Contributions

Xin Huang and Guangchen Liu performed the experiments and analyzed the test data; Xin Huang and Yu Chen wrote the manuscript; and Jianyong Pang conceived the experiment and provided the guidance and suggestion.

Acknowledgments

The authors acknowledge the School of Civil Engineering and Architecture, Anhui University of Science and Technology, for their support to perform this research. This research was funded by the Science and Technology Project Foundation of Key Technologies for Prevention and Cure of Major Accidents in Production Safety, General Administration of State Security Supervision (No. Anhui-0003-2016AQ).