Study on the Damage Model of Coal Rock Caused by Hydraulic Pressure and Electrical Impulse in Borehole

School of Civil Engineering, Inner Mongolia University of Science and Technology, Baotou Inner Mongolia 014010, China State Key Laboratory of Coal Resources and Safe Mining, China University of Mining and Technology, Xuzhou 221116, China School of Mining and Coal, Inner Mongolia University of Science and Technology, Baotou Inner Mongolia 014010, China School of Mining Technology, Taiyuan University of Technology, Taiyuan 030024, China


Introduction
The coalbed methane (CBM) is an associated mineral resource with methane as the main component stored in the coalbed. Its occurrence state is mainly adsorbed on the surface of coal matrix particles and partially freed from coal pores. It is unconventional natural gas. It is clean, highquality energy, and chemical raw material that has risen internationally in the past 1~2 decades. China is the third CBM reserve country in the world and has rich and evenly distributed coalbed gas [1,2]. Nevertheless, the CBM of China is generally characterized by low saturation, permeability, reservoir pressure and resource richness, and high metamorphism [3]. The current extraction rate is generally 30%~40%, and the mining efficiency is low. Some countries with a higher CBM extraction rate generally have an extraction rate of 30%~60% (the maximum extraction rate is 80%) [4]. At present, there are two main methods for improving the extraction rate of low-permeability coalbed methane. One of them is to change the occurrence state of CBM. The gas adsorbed on the inner surface of e pores of the coal matrix is changed from the adsorption state to the free state, thus improving the diffusion ability of CBM from the matrix and the micropores to the crack; You et al. [5] studied the effect of three surfactants in the fracturing fluid on the adsorption-desorption processes of coalbed methane (CBM). It was studied by conducting an isothermal adsorption-desorption experiment. The adsorption capacity of surfactants and the wettability of the coal surface were also systematically investigated to illustrate the impact of surfactants. The effect of fracturing fluid was proved to be of great significance for coalbed methane desorption and improved coalbed methane extraction rate. Zhu et al. [6] studied CBM in composite desorption mode and found that the desorption rate of CBM was negatively correlated with the water content, and the desorption amount showed an increasing trend with the decrease of the water content. The other is to expand the migration channel of gas; through the external load, the expansion of internal cracks of the coal seam is improved, and more complex network cracks are formed. It provides a channel for gas to flow inside the coal seam and allows the CBM to pass more cracks to the borehole (wellbore), which increases the gas permeability and improves the extraction efficiency. Li et al. [7] proposed to develop a method to increase the permeability of coal seam. It is realized by dissolving part of the coal matrix solvent and combining it with a physical destruction method to solve the problem of low permeability coal seam fracture connectivity. Lu et al. [8] proposed the high voltage electric pulse fracturing coal technology. He proved that this technology effectively and improved the coal seam fracture connectivity and increased its transport efficiency. It improved the coalbed methane production. Yan et al. [9] had proposed the synergistic technology of slitting and fracturing to increase the permeability, which has significantly improved the efficiency of gas extraction, in response to the problem of inadequate pressure release in the drill holes for extraction through the low permeability coal seam. Compared with the first method, the second method is more widely used and has a stronger guiding significance for engineering practice.
Among the various methods to expand the CBM migration channel and improve the extraction efficiency of CBM wells, the reservoir-cracking and permeability improvement technique has the most significant influence in increasing the production of CBM. Therefore, it has been highly valued and favored by CBM well extraction technicians in various countries [10,11]. More than 90% CBM wells in the United States use various cracking techniques to increase the production of CBM. The application and development of various cracking and antireflection technique in China are also very rapid. Most high-yield CBM wells have adopted this stimulation technique. This means that the artificial methods for the structural transformation of the coal-rock reservoir layer, the increase of crack extension range, the weakening of coal-rock strength, and the improvement of permeability are effective ways to improve the extraction efficiency. Among them, hydraulic fracturing is not only used for testing seismic sources [12,13] but also has been widely used in the development and exploitation of low-permeability oil and gas fields as a major measure to increase production in oil and gas wells. Zhong et al. [14] proposed that the stimulation of hydraulic fracturing on a combination of depressurization with thermal stimulation for natural gas hydrate (NGH) exploitation from challenging ocean hydrate reservoirs is investigated in this study. The influence of fracture permeability, well spacing on hydrate exploitation, and the enhancement effect of hydraulic fracturing on different well spacings is numerically investigated. The results show that hydraulic fracturing has a significant improvement in hydrate exploitation performance, including hydrate decomposition behavior, gas production behavior, gas recovery ratio, and energy ratio. Zhou et al. [15] conducted a simulation and experimental study on hydraulic fracturing of the low-permeability coal seam to solve the difficult problems of extraction, heavy workload, and low efficiency of gas extraction in a low-permeability coal seam. The results showed that the fracture damage range of coal seam before and after hydraulic fracturing was up to 14 m, and the permeability of the coal seam increased 67 times. Wu [16] studied the propagation law of hydraulic fracturing in the coal rock complexes, and the maximum opening of hydraulic fracturing in the parallel stratigraphic direction of the coal was 1.292 mm, the maximum opening in the vertical stratigraphic direction of the coal was 0.952 mm, and the maximum opening in the vertical stratigraphic direction of the sandy mudstone was 0.816 mm after the implementation of hydraulic fracturing.
Hydraulic fracturing has proven to be very effective in gas extraction and treatment, and it has also exposed some problems. Wang [17] conducted hydraulic fracturing test for coal seam wells with low permeability. It realized the purpose of pressure relief and permeability enhancement in lowpermeability coal seam, which greatly increased the extraction amount of gas from the low-permeability coal seam. Wang et al. [18] implemented the technique of directional hydraulic fracturing, where slits and hydraulic fracturing widen the range of gas pressure drop and reduce the time required to complete extraction. While hydraulic fracturing was carried out, some problems were also exposed. Gagnon et al. [19] suggested that environmental pollution problems such as methane pollution and other chemical pollution are the key problems to be solved based on hydraulic fracturing technology for shale gas wells. Bao [20] conducted a comparative experiment t based on high-pressure electric pulse hydraulic fracturing and traditional hydraulic fracturing. He found that traditional hydraulic fracturing technology exists limitations. They are insufficient fracturing capacity and long fracturing time in the application process. In order to further improve the traditional hydraulic fracturing effect, some new hydraulic fracturing techniques have emerged in recent years. Based on the electrohydraulic effect, Li et al. [21,22] proposed a high-voltage electric pulse hydraulic fracturing and permeability improvement technique of coal seam. The principle of this technique is to use the water shock vibration effect formed by the high-voltage electric pulse discharge in the crack tip of coal seam to break, develop, and penetrate the microcracks, thus achieving the purpose of reducing the resistance and improving the permeability and extraction efficiency of CBM. The specific principle was shown in Figure 1. The electrohydraulic effect produced by a high-voltage discharge pulse can also be referred to as the electrical explosion in the liquid (or electrohydraulic explosion). Li [23] and Yan et al. [24,25] studied the static breakdown, detonation, and heating effects of pulse discharge in hydrostatic pressure, the breakdown and detonation dynamic load, and the cracking characteristics of coal-rock mass through experiments and numerical calculations. In order to study the method in depth, Bian et al. [26,2 Geofluids 27] analyzed the pulse wave and expansion mechanism of pulse discharge with different hydrostatic pressures and discharge voltages under the action of coupling hydrostatic pressure and studied the fracturing effect of concrete specimen with different hydrostatic pressures and discharge parameters. Bao et al. [28,29] studied the coupling mechanism and the electro-solid-liquid law and used the high-voltage detection and CT scanning technique to preliminarily evaluate the damage characteristics and fracturing cracks effects of for high-voltage electric pulse hydraulic fracturing coal with different test parameters.
In recent years, in order to explore the micromechanical mechanism of macroscopic fracture of hydraulic fracturing coal-rock mass, domestic and foreign scholars have carried out much types of research on the microscopic damage laws of coal-rock mass and macroscopic fracture characteristics. Zhao et al. [30] studied the proposed equivalent Burgers model lays the foundations for revealing the rheological behavior of fractured rock mass subjected to the combined action of hydraulic pressure and stresses. Peng et al. [31] studied the macroscopic fracture characteristics of granite under different loading conditions based on different temperatures and depths. Wang et al. [32] studied the fracture characteristics of granite under impact loading and found that the fracture zone area and crater depth are variables depending on the impact angle at constant impact velocity. Based on the characteristics of probability and statistical models, Hu et al. [33] consider adsorbed-gas-induced swelling stress and erosion, pore/fracture-induced damage and failure to coal skeleton, the new effective stress equation, damage deformation control equation, and constitutive model for gas-bearing coal which are established. Based on the principle of statistics and fractal theory, we combined the microstructure characteristics of coal and established the control equation of the fracture field distribution. Based on the derivation of coal matrix strain, Liu et al. [34] established a coal matrix stress model under various stress conditions and a mechanical model for pulverized coal crushing. Tang et al. [35] studied the microdamage characteristics of sandstone in the loading process of different confining pressures and axial stress by using a triaxial compression test equipment. With the increase of axial stress, the degree of plastic deformation and the microdamage factor of sandstone increase continuously. The strength of sandstone gradually deteriorates, and the microdamage factor of the specimen basically obeys the Gauss function distribution. The electrohydraulic effect is produced by the discharge of high voltage and current transformer in water which converts the electromagnetic energy to the mechanical energy (strong shock wave). The shock wave amplitude pressure is more than 100 MPa, and the pulse width is several hundred microseconds.

Energy converter
Coal seam

Geofluids
The microdamage of specimen increases with the axial stress from the edge of the observation surface to the center. Voyiadjis [36] analyzed the damage variable in the coalrock damage process. Based on the evolution of the fatigue damage variable of coal-rock mass under the action of cyclic loading, Xiao [37] proposed a time-dependent damage variable and derived the damage constitutive model under the action of cyclic loading. Liu et al. [38] found that the gas permeability of coal seam was the main factor that determines the efficiency of gas extraction in mines. It was closely related to coal damage. The permeability of gasbearing coal increases exponentially with the development of damage. Li et al. [39] comprehensively discussed the effective stress of coal-rock mass under the action of external load and established a statistical constitutive model of elastoplastic damage. Based on the variable amplitude loading, Liu [40] derived the variation law of fatigue damage variable and deduced and established the constitutive model.
At present, the dynamic mechanical characteristics, dynamic damage mechanism, damage mode and crack initiation, and evolution mechanism of coal-rock mass under the action of high-voltage discharge pulse water shock have still been in an exploration stage. In order to study the fracturing effect of coal-rock mass with different electrohydraulic coupling parameters, the damage of microcracks in the damage process of coal-rock mass and the occurrence and expansion of macrocracks were quantitatively evaluated. Zhao et al. [41][42][43] investigated fracture toughness and subcritical crack growth of different rock and studied weakening the effects of different parameters. Based on Lemaitre damage strain equivalent principle [44][45][46], the damage variable of coalrock mass was defined, and the damage constitutive equation was established by combining the damaged area with effective stress of fine microfracture, which achieved the purpose of quantitatively evaluating the specific damage of coal-rock mass under the action of electrohydraulic coupling. The high-pressure electric pulse hydraulic fracturing test under the action of electrohydraulic coupling was designed and implemented. Based on the hydrostatic pressure and the electric pulse load electrohydraulic coupling, the numerical simulation was carried out by the ABAQUS/XFEM [47,48]. Based on the unit damage area and effective stress, the feasibility and correctness of damage variables and constitutive models were verified by the damage and deformation characteristics of coal-rock mass under the action of different loads.

Establishment of Electrohydraulic Coupling
Damage Model for Coal-Rock Mass

Mechanism of Electrohydraulic Coupling Damage.
Under the action of hydrostatic pressure and electric pulse load, the damage caused by coal-rock mass is mainly divided into two stages. The first stage is the damage of coal-rock mass under the action of hydrostatic pressure. There are a lot of microdefects inside the coal rock. When the hydrostatic pressure is applied to the coal-rock mass, the original microcracks existing in the coal-rock mass itself are initiated, expanded, and merged. This is the first stage of damage, that is, the solidliquid coupling damage stage. Under the action of hydro-static pressure, the electric pulse load is applied. As the pulse load repeatedly impacts, the crack expands and penetrates continuously, and the damage degree of coal rock increases continuously. This is the second stage, namely, the electrohydraulic coupling damage stage. The damage state effect can be represented by the damage variable as follows: where D 2 is the total damage variable of coal-rock mass under the action of hydrostatic pressure and electric pulse loading electrohydraulic coupling, D 1 is the damage variable under the action of hydrostatic pressure of coal-rock mass, D 1 ′ is the damage variables in coal-rock mass under electrical impulse loading, and D 1 D 1 ′ is the damage coupling term.

Damage Evolution Constitutive Relation.
In the process of high voltage electric pulse hydraulic fracturing of coalrock mass, the damage crack caused by the hydrostatic pressure and electric pulse pressure can be divided into macrocrack and micr-crack. The change of the macrofracture area is resulted from the generation and penetration of microcrack. Therefore, in the process of electric pulse hydraulic fracturing of coal-rock mass, the damage of the effective bearing area of coal-rock mass can be used to characterize the damage during the entire damage process. Based on the effective bearing area of coal-rock mass, the damage variable can be expressed as where A 1 is the original effective bearing area of coal-rock mass, m 2 , and A 2 is the effective bearing area of coal-rock under the action of electrohydraulic coupling, m 2 . As early as 1906, Professor Lemaitre proposed the strain equivalence criterion: where σ 1 is the original effective stress of coal-rock mass, MPa, σ 2 is the effective stress of coal-rock mass under the coupling of liquid and electricity, MPa, E 2 is the elastic modulus of coal-rock mass under the action of electrohydraulic coupling, MPa, and ε is the damage strain of coal-rock mass. According to the Lemaitre equivalent strain hypothesis, Zhang [49] proposed the strain equivalence principle: when the coal-rock material is subjected to the force F, two damage states of material are taken. The effective stress in the base damaged state acting on the strain caused by the second damage state is equivalent to the effective stress of material in the second damage state acting on the strain caused by the first damage state. Therefore, the following formula was obtained.
It can be seen from Formula (4) that if the coal-rock mass is defined as the base damage under the stress state, 4 Geofluids the original effective stress value of coal-rock mass can be obtained according to the derivation formula of Tang [50]: where σ x is the horizontal ground stress, MPa, σ y is the vertical ground stress, MPa, p 1 , the atmospheric pressure, is 0.1 MPa, and p 2 is the gas pressure in coal-rock mass, MPa.
The high-voltage electric pulse hydraulic fracturing coalrock mass is cyclically loaded by dynamic and static loads. According to Terzaghi's effective stress principle, the effective stress at a certain moment in the process of electric pulse hydraulic fracturing test can be obtained.
where σ is the total stress of hydrostatic pressure pulse stress wave. The peak stress of shock wave in the experiment can be taken.
Based on the stress change, the damage variables of electric pulse hydraulic fracturing coal-rock mass were obtained by combining Formulas (2), (4), (5), and (6): Combine Formula (2) with Formula (3), and the following formula was obtained: Therefore, the stress-strain constitutive relationship of electric pulse hydraulic fracturing was obtained: E 1 in the constitutive relation (9) is the initial elastic modulus of coal-rock mass, which is relatively easy to measure. Therefore, it is not necessarily to measure the elastic            Table 1.

Design of Test Scheme.
The high-voltage electric pulse hydraulic fracturing system independently developed by the research group was adopted. The test system was shown in Figure 2. The controlled trial of two groups with 8 coal samples was carried out. Coal samples of the first group were loaded into a rigid triaxial pressure chamber. After the filling was completed, the same ground stress was applied to the coal samples. The horizontal confining pressure is 8.66 MPa, and the vertical axial pressure is 7.28 MPa. After the stress was loaded, 3 MPa hydrostatic pressure was applied to the pipeline and coal samples firstly, and then the high-voltage discharges of 7 kV, 9 kV, 11 kV, and 13 kV    9 Geofluids were carried out, respectively. The coal sample was taken out after loading. The second group of coal samples was replaced. After the same stress loading was completed, the 11 kV dis-charge voltage was kept unchanged, and the system hydrostatic pressure was changed. The hydrostatic pressures of 1 MPa, 4 MPa, 6 MPa, and 8 MPa were applied, respectively.

Geofluids
During the test, the waveform and pressure of the stress wave in the water were collected by the high-frequency pressure sensor and the transient signal recorder. After the test, the CT scan detection equipment was used to detect the damage of the coal-rock mass. The specific scheme was shown in Table 2.

Test and Result Analysis of Stress Wave Peak Pressure.
The waveform and pressure of pulse stress wave (water shock wave) were collected and recorded during the experiment by using the high-frequency pressure sensor and the high-speed data-storage device. The results were shown in Figure 3. Then, the relationships of the stress wave peak pressure with the voltage and the water pressure were analyzed, respectively. The results were shown in Figure 4.
As shown in Figures 3 and 4, the maximum peak pressure of water shock after the discharge in water can be measured. After the high-pressure discharge in water, the rising edge of the water shock wave is steep. After reaching the maximum peak value, it decays in an approximate exponential form, and several peaks of the waveform appear. The oscillation decays into a vibration wave (seismic wave) finally. At the same time, it can be found that the maximum peak stress of water shock increases with the increase of voltage when the water pressure is constant. The higher the voltage, the greater the peak stress. When the discharge voltage is constant, the maximum peak stress of water shock also increases with the increase of water pressure, but the peak stress changes less than the amplitude of voltage change.

Damage Variable Calculation.
In order to quantitatively study the microdamage of coal samples under the action of coupling of different water pressures and discharge voltages, the damage variables of coal samples were calculated and analyzed. The ground stress values of coal samples and stress value of coal seam gas were substituted into formula (5), and the measured original effective stress of coal-rock mass was 8.1 MPa. The measured maximum peak stress of water shock wave was brought into formula (6), and the effective stress of coal rock under the action of liquid-electric coupling was obtained, as shown in Figure 5. Then, the damage variable value of coal samples with different liquid-electric parameters was obtained by substituting the effective stress into equation (7). The specific values were shown in Figure 6.
It can be seen from Figures 5 and 6 that under the condition of the constant hydrostatic pressure of 3 MPa, the discharge voltage value increases from 7 kV to 13 kV; the voltage growth rate is 185.7%, and the damage variable of coal sample increases from 0.575 to 0.666. The growth variable of the coal sample has a growth rate of 115.8%. The damage variable of the coal sample increases with the increase of discharge voltage. If the discharge voltage is 11 kV, the growth rate of hydrostatic pressure is 800%; the damage variable of the coal samples also increases with the increase of hydrostatic pressure, and its increase rate is only 106.1%. Relatively speaking, the voltage increases by about 0.86%, and the damage variable of the coal sample increases by 15.8%. The hydrostatic pressure increases by 700% while the damage variable of the coal sample only increases by 6.1%. It is indicated that the change of relative hydrostatic pressure and voltage has a more significant influence on the degree of damage of coal-rock mass, and the coal-rock mass has a better cracking effect.

Analysis of CT Test Results.
In order to further verify the microdamage inside the coal-rock mass, after the highvoltage electric pulse hydraulic fracturing of the coal-rock specimen, the tested coal-rock mass was sampled and analyzed, and then the sample was damaged by the CT test equipment. The CT scanning results were binarized by means of the Particle and Crack Analysis System (PCAS) developed by Nanjing University in which the black is the crack part and the white is the coal mass. The CT scanning results and the corresponding binary images obtained after PCAS treatment were shown in Figures 7 and 8 (the drilled specimens can be constructed with more than 1500 layers after scanning; only representative layers are selected, and the remaining layers are similar).
It can be seen from Figure 7 that after the CT scanning detection, binarization treatment, and 7 kV and 9 kV discharge under the condition of maintaining the hydrostatic pressure of 3 MPa, after, three irregular and intermittent cracks occur in the coal slice scanning section, and the length and the opening degree of cracks are small. Under the 11 kV discharge condition, five intermittent cracks occur in the coal seam scanning slice, and the number, length, and damage expansion degree of cracks increase. When the voltage increases to 13 kV, the penetration, opening degree, length, and roughness of cracks increase obviously, and the damage expansion degree is intensified. In general, under the condition of maintaining the hydrostatic pressure of 3 MPa, with the increase of discharge voltage, the number, total length, expansion degree, penetration, damage area, and crack rate of coal-rock specimens increase linearly. The trend indicates that the damage degree of coal rock increases with the increase of discharge voltage.
It can be seen from Figure 8 that when the discharge voltage is 11 kV and the hydrostatic pressure is 1 MPa, three intermittent cracks occur in the coal slice scanning section; when the hydrostatic pressure is 4 MPa, the crack penetrates to form the main crack, and its roughness and opening degree increase. Many short and irregular cracks occur around the main crack; when the hydrostatic pressure is 6 MPa, the main crack branches, and its length, opening x = 8.66 MPa y = 7.28 MPa Figure 9: Loading model. 12 Geofluids degree, and roughness further increase; when the pressure is increased to 8 MPa, three main cracks are generated. One has a long length, and another one has a large opening degree. The crack rate of coal-rock specimens further increases. Compared Figure 7 with Figure 8, it can be seen that the damage degree of coal-rock specimens is directly related to the change of hydrostatic pressure and discharge voltage, but the change of discharge pressure has a greater influence on the damage degree and crack rate of coal-rock specimens than that of hydrostatic pressure. It is consistent with the evaluation results of coal-rock damage caused by the macroscopic damage variables of coal-rock mass derived from the interaction of effective stress.

Microdamage Numerical Analysis
In order to further verify the damage model, it is more comprehensive to evaluate the damage degree of high-pressure electric pulse hydraulic fracturing technique on the coalrock mass quantitatively and to further deduct the damage variable of coal-rock mass. The damage variable was derived from the change of the effective bearing area of coal-rock mass by the numerical simulation, and the damage degree and expansion of microcrack were further evaluated. According to the definition of damage variable, Formula (10) was obtained: whereA 1 and A 2 have the same meaning as those in Formula (2); A D is the damage area of coal-rock under the electrohydraulic coupling.
The change of bearing area of the microjoint unit of the coal-rock mass was measured by the high-voltage electric pulse load. The change of damage variable was obtained by Formula (10), and the damage morphological change of coal-rock mass was further characterized. The length of the cell is 2 cm. Therefore, the original effective bearing area of the unit section is 4 × 10 −4 m 2 . 8.66 MPa and 7.28 MPa ground stress were applied in the horizontal direction and in the vertical direction, respectively. The model was shown in Figure 9. The model parameters were shown in Tables 3 and 4.

Damage Degree of Coal Rock at 3 MPa Hydrostatic
Pressure and Different Voltages. As shown in Figure 10, 7 kV high-voltage pulse discharge was performed in 3 MPa hydrostatic pressure. The maximum opening degree of cracks in coal-rock mass is 1:21 × 10 −2 m, and the maximum damage area of the crack unit is 1:21 × 10 −2 × 2 × 10 −2 = 2:42 × 10 −4 m 2 (where2 × 10 −2 is the original length of cell). The original effective bearing area of the unit section is 4 × 10 −4 m 2 . The damage variable is 0.605 by substituting the above data into Formula (10); the effective bearing area of the cell after loading is 2 × 10 −2 × 2 × 10 −2 − 2:42 × 10 −4 = 1:68 × 10 − 4 m 2 . At the same time, under the condition of constant hydrostatic pressure, the effective bearing area and damage variable at other discharge voltages can be obtained, as shown in Figure 11. It can be seen from Figure 11 that with the increase of electric pulse load, the effective bearing area of coal-rock mass decreases, and the damaged area and damage variable increase. It is indicated that with the increase of discharge voltage and stress peak generated in the water, the damage degree gradually increases.

Damage
Degree of Coal-Rock Mass at 11 kV Voltage and Different Hydrostatic Pressures. As shown in Figure 12, 11 kV high-voltage pulse discharge was performed in hydrostatic pressure of 1 MPa. The maximum opening degree of cracks in coal-rock mass is 1:39 × 10 −2 m, and the maximum damage area of the crack unit is 1:39 × 10 −2 × 2 × 10 −2 = 2:78 × 10 − 4 × m 2 (where 2 × 10 −2 is the original length of cell). The original effective bearing area of the unit section is 4 × 10 −4 m 2 . The above data can be substituted into the Formula (10) to obtain the damage variable of 0.695. At the same time, the effective bearing area of the loading cell is 2 × 10 −2 × 2 × 10 −2 − 2:78 × 10 −4 = 1:22 × 10 −4 m 2 . When the discharge voltage is constant, the effective bearing area and damage variable at other hydrostatic pressures can be calculated, as shown in Figure 13. It can be seen from Figure 13 when the discharge voltage remains unchanged, with the increase of hydrostatic pressure, the effective bearing area of coal-rock mass decreases, and the damaged area and damage variable increase. It is indicated that with the increase of discharge voltage and stress peak generated in the water, the damage degree gradually increases.

Discussion
When the high pressure is discharged in pressurized water, the stress wave generated by the electric explosion can make the defects crack, expand, and penetrate the coal-rock mass, such as voids and microcracks. Therefore, it is especially important to study the fracture and expansion behaviors of microcracks. However, how to test and evaluate the initiation, expansion behaviors, and effects of microcracks accurately and quantitatively has always been a problem that draws scholars' attention.    14 Geofluids From the above test and numerical simulation, combined with the observation results of CT scanning, the microscopic damage of coal-rock mass was accurately evaluated by calculating the damage variable of coal-rock mass. Then, the macroscopic damage behavior of the coal-rock mass structure was revealed. Compared with the damage variable at the same loading, the comparison diagram of damage variables was obtained, as shown in Figure 14.
It can be seen from Figure 14 that the stress-based damage variable obtained by the experiment and the damage    variable based on the damaged area obtained by numerical simulation has almost the same change trend at the same loading, and their values are not much different. This is mainly because the isotropic nature of material parameters during the simulation causes small changes of values in the test and simulation. These two research methods confirm each other, which indicates that the damage variable of coal-rock mass can reflect the internal microdamage evolution, the expansion, and the extension process of macrocrack. Therefore, the damage degree of coal-rock mass can be quantitatively characterized. Based on the relationships between the Newtonian fluid mechanics and the damage mechanics, and between the damage evolution characteristics and the macrocrack occurrence, the damage model was established by using the change of effective stress, damage area, and damage variable to better characterize the variation of effective stress, effective bearing area, and damage variable of coal-rock mass under the action electric pulse hydraulic fracturing electrohydraulic coupling effect. The damage evolution of coal-rock mass materials from the initial cracking to complete cracking under the action of electrohydraulic coupling was revealed, and the macrodamage degree of coal-rock mass was quantitatively characterized, which provides an important way and method for the application of engineering practice.
Despite the establishment of damage variables and damage models based on the effective stress and effective bearing area of coal-rock mass, the damage degree of coal-rock mass under the action of electrohydraulic coupling was quantitatively evaluated by the numerical simulation and experiment, but the current research is still in the preliminary exploration stage. Some difficult issues remain to be further studied.
(1) Due to the limitations of the detection technique, the damage width and area of microcracks in coal-rock mass at this stage cannot be quantified by test and detection. At this stage, only the numerical simulation can be used to derive the results (2) The quantitative test for the high-voltage electric pulse hydraulic fracturing is still in the research stage of single crack expansion, and the research on the composite crack change and mutual influence is slightly insufficient (3) The coal-rock mass expansion models adopted in most studies at home and abroad are based on the isotropic assumption of material, but in the engineering practice process, the coal-rock material and crack expansion are anisotropic. Therefore, the research on the anisotropy of materials and its influence on crack fracture and expansion should be strengthened in the future (4) Due to the nonuniformity and permeability of the coal-rock structure, the existing numerical simulation methods for high-voltage electric pulse hydraulic fracturing still have certain defects, such as difficulties in handling the multicrack expansion and acquiring the calculation parameters required by prefabrication cracks or numerical simulation. At the same time, the factors affecting the hydraulic crack expansion of high-voltage electric pulse in coal-rock mass need to be further studied systematically

Conclusion
(1) Based on the correlation between macroscopic defects and damage evolution of coal-rock mass, a dynamic damage constitutive model of coal-rock mass was constructed. The correctness of the model under the action of electric-liquid coupling was verified by experiments and numerical calculations (2) The effective bearing area, effective stress, and damage variable can be used to evaluate the damage degree of macroscopic cracks of coal-rock under the action of electric-liquid coupling quantitatively and effectively, thus revealing the damage law of coal-   16 Geofluids rock mass under the coupling effect of hydrostatic pressure and electric pulse load (3) Through tests and numerical simulations, it can be seen that compared with the change of water pressure, the change of voltage has more influence on the damage of coal-rock mass during the highvoltage pulse hydraulic fracturing process. The change of damage variable of coal-rock mass is more significant (4) With the increase of electric pulse load, the effective bearing area of the coal-rock bearing unit decreases, and its damage variable increases. The damage variable is positively correlated with the pulse load and negatively correlated with the effective bearing area of the unit. The damage variables obtained in the experiment are consistent with the calculation values of numerical simulation

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

Conflicts of Interest
The author(s) declare(s) that they have no conflicts of interest.