Effects of different fracture parameters on microseisms induced by hydraulic fracturing

The microseism induced by hydraulic fracturing is of great significance to the development of geothermal reservoirs and the site selection of geothermal systems. In this study, taking geothermal field data from Qiabuqia as a geological reference, several models are developed for hydraulic fracturing simulations based on the models included in the commercial FracMan™ software suit. A series of numerical simulations are carried out to explore the effects of different fracture parameters, including the existence of faults, in-situ stress state, fracture occurrence, and fracture distribution near faults, on the induced microseisms. The results show that during hydraulic fracturing, the existence of faults does affect the propagation direction of the newly generated fractures, causing the fractures to extend toward the fault. The faults also increase the magnitude of microseismic events and make the distribution of microseismic events farther. Under three different in-situ stress states, the magnitudes of induced microseisms and the distance between microseismic events and injection wells are different, mainly due to the different forms of energy release. The total energy of induced microseisms under the reverse faulting stress state (RF) is the largest. Under the normal faulting stress state (NF), the number of microseismic events is the least, with only 28, and these microseismic events are concentrated near the well, with the farthest distribution distance of only 150 m. The fracture occurrence has a significant effect on induced microseisms, mainly affecting the number and distribution range of induced microseisms. With the same number of fractures, as the fracture concentration to the fault increases, both the maximum magnitude and the farthest distribution distance of the induced microseisms increase.


Introduction
Nowadays, hydraulic fracturing technology has been widely used in the development of geothermal energy. It not only improves the permeability of the rock in the reservoir, but also creates complex fracture networks in the reservoir as convective heat transfer channels, thereby increasing the productivity of the Enhanced Geothermal System (EGS) (Kruger 1976;Huenges and Ledru 2011;Zhu et al. 2015). However, many studies have pointed out that a large number of microseismic events occur during hydraulic fracturing (e.g., Warpinski et al. 2004;Vermylen and Zoback 2011;Anikiev et al. 2014). Especially when the fracture network of hydraulic fracturing is connected with the faults, it tends to produce some earthquakes with larger magnitude (Maxwell et al. 2009;Zhang et al. 2020). Therefore, it is necessary to study the effects of different fracture distributions, especially those near faults, on the microseismic induced by hydraulic fracturing, which can effectively avoid larger earthquakes and reduce damage to the environment.
Over the past few years, many studies have been conducted on the factors that influence hydraulic fracturing (e.g., Quosay et al. 2020;Ries et al. 2020;Liu et al. 2021;Qiao et al. 2022;Zheng et al. 2022). Generally, hydraulic fracture propagation proceeds through the following four steps: (1) the formation of new hydrofracture, (2) the reopening of pre-existing tensile fractures, (3) the opening of transverse natural fractures, and (4) the initiation of shear slip on the pre-existing fractures. Lamont and Jessen (1963) first investigated the effects of pre-existing fractures on the propagation of hydraulic fractures under tri-axial stresses in the laboratory, and they demonstrated that the width and orientation of existing fractures do not alter the extension or direction of the hydraulic fracture. Later, Cheng et al. (2015) also studied the effect of pre-existing fractures on hydraulic fracture propagation based on a series of hydraulic fracturing tests, and they proposed two empirical relationships to describe the interactions between pre-existing fractures and hydraulic fractures. De Pater and Beugelsdijk (2005) pointed out that when hydraulic fractures and natural fractures intersect, it may lead to the arrest of fracture propagation, the flow of fluid into discontinuities, the creation of multiple fractures, and the offset of fractures. Similarly, Dehghan et al. (2015) studied the effects of natural fracture dip and strike of pre-existing fractures on hydraulic fracture propagation when the pre-existing fractures intersect with hydraulic fractures, and they found that after interaction with the pre-fracture, hydraulic fracture geometry was significantly affected by differential horizontal principal stress and the pre-fracture dip and strike. Liu et al. (2014) developed a new experimental model to simulate the effect of a natural fracture network on the propagation geometry of hydraulic fractures in naturally fractured formations, and they gave the principle of hydraulic fracture propagation that it follows the least resistance, the most preferential propagation, and the shortest propagation path. The common conclusion from these studies shows that high horizontal stress difference and large approaching angle are conducive for the hydraulic fractures to cross the pre-existing fracture.
So far, most studies have only focused on the effect of natural fractures on hydraulic fracture propagation, while the effects of natural fractures on microseisms induced by hydraulic fracturing are less studied, especially the effects of the interaction between natural fractures and faults on induced microseisms (Dong and De Pater 2001;Zeng and Liu 2010;Behnia et al. 2015;Li et al. 2020). For example, Troiano et al. (2013) provided an interpretation of induced microseisms based on the computation of Coulomb stress changes resulting from fluid injection or withdrawal at depth, and they demonstrated that seismicity was more likely to occur where Coulomb stress varied more. Based on the microseismic monitoring system, Liu et al. (2016) found that the induced microseisms were mainly concentrated at the footwall of the fault, and the initiation and propagation of fractures were much faster under strong disturbance. Shan et al. (2021) studied the effects of prefabricated faults on microseisms induced by hydraulic fracturing through a series of hydraulic fracturing tests on granite samples and demonstrated that the prefabricated faults could lead to an increase in shear fractures and fracture energy. Since the characteristics of natural fractures are closely related to the site selection of geothermal systems, the optimization of reservoir reconstruction, and the measures of environmental protection, it is necessary to study the effects of natural fracture characteristics on the microseisms induced by hydraulic fracturing.
In this study, taking geothermal field data from Qiabuqia as a geological reference, a number of fracture models were developed by using the FracMan software to simulate hydraulic fracturing under different natural fracture distribution characteristics. The main objective of the present work is to explore the effects of different fracture parameters, including the existence of faults, in-situ stress state, fracture occurrence, and fracture distribution near faults, on the microseisms induced by hydraulic fracturing. The findings of this study not only fill the research gap of the effects of natural fractures on induced microseisms, but also provide a theoretical basis for the site selection of geothermal systems and the optimization of reservoir reconstruction.

Numerical simulation method
In this study, the fractures of rock mass are characterized using the models included in the commercial FracMan™ software suite (Golder Associates Inc.). These built-in models have been proven to be able to accurately simulate the regional in-situ stress field and volume pressure fracture network based on the given microseismic data (e.g., Crosta 1997;Johri and Zoback 2013;Nadimi et al. 2018;Zheng et al. 2019;Makedonska et al. 2020). According to the size, orientation, and strength characteristics of natural fractures in actual geothermal field data, the numerical models are established based on the discrete fracture network to simulate the hydraulic fracturing process under different fracture conditions.
The models used in this study have the following three main assumptions: (1) the flow in the fracture is an incompressible fluid and the seepage law conforms to Darcy's law; (2) the rock is an impermeable matrix and the water only flows in the fracture. The fluid flow in the rock matrix can be ignored; (3) the chemical reaction between fluid and rock matrix and the roughness of fracture are not considered.
The pore pressure of fractures is assumed to decrease linearly from the well-bore pressure in order to simulate the leak-off to the matrix or non-dilatable fractures and the pressure drop due to the resistance of fracture surface to flow. Based on Eq. (1) proposed by Secor Jr and Pollard (1975), the fracture pressure model is established, and the fracture pore pressure, P frac , can be calculated as follows: (1) where P pump is the pumping pressure; N is the fracture normal stress; d refers to the flow distance from the well intersection point to the fracture element, that is, the distance of water flowing along the hydraulic fracture generated in the injection well; d max represents the maximum flow distance, which can be automatically adjusted according to the length of hydraulic fracture at each time step; is the fluid density; and Δh is the vertical elevation difference from the well-intersection point to the fracture element. Equation (2), proposed by Secor Jr and Pollard (1975), is adopted to calculate the fracture aperture, e. This equation is similar to the width-pressure relation equation in elliptical elastic fracture, as follows: where v is the Poisson's ratio. The Mohr-Coulomb strength criterion is used to determine the fracture shear strength, f , as shown in Eq. (3) (Griffiths 1990): where S 0 is the cohesion of the fracture; N is the fracture normal stress, and f is the friction angle. ( The extension and propagation of hydraulic fractures are realized based on the critical stress analysis theory. Since hydraulic fracturing is achieved by maintaining a volume balance between the volume of injected water and natural fracture propagation and the volume of developed hydraulic fracture, the distribution of the active fracture network and the direction of fracturing fluid can be determined by judging the activation relationship between the main hydraulic fracture and the natural fracture. Generally, there are two activation relationships between artificial and natural fractures. One activation relationship occurs in the case that the tip pressure of fracturing fluid is less than the normal stress acting on the natural fracture surface, as shown in Fig. 1. As the hydraulic fracture gradually approaches the natural fracture, the fracturing fluid also enters the natural fracture from the hydraulic fracture (Fig. 1a). Since the tip pressure of fracturing fluid is less than the normal stress acting on the natural fracture surface, the fracturing fluid first forces the natural fracture to open wider (from Fig. 1b to c). When the tip pressure reaches a certain threshold, a new secondary fracture is generated (Fig. 1d). As the pressure of the fracturing fluid continues to increase, the fracturing fluid can cause both natural and secondary fractures to expand simultaneously until they can no longer expand (Fig. 1e), after which the fracturing fluid continues to flow along the hydraulic fracture (Fig. 1f). The other activation relationship occurs in the case that the tip pressure of fracturing fluid is greater than the normal stress acting on the natural fracture, as shown in Fig. 1 Activation relationship between artificial and natural fractures in the case that the tip pressure of fracturing fluid is less than the normal stress acting on the natural fracture surface, where the orange arrows represent the flow direction of the fracturing fluid, the green lines represent natural fractures, the orange lines represent newly created secondary fractures, the triangular shaped grey area surrounded by yellow lines represents the hydraulic fracture, and σ 3 is the minimum principal stress of in-situ stress Fig. 2. Different from the first activation relationship (Fig. 1), since the tip pressure of fracturing fluid is greater than the normal stress acting on the natural fracture surface, after filling with natural fracture (Fig. 2b), the fracturing fluid first creates the new secondary fracture instead of expanding the natural fracture (Fig. 2c). Until the fracturing fluid does not allow the secondary fracture to continue to expand, the fracturing fluid begins to expand the natural fracture (Fig. 2d), after which the fracturing fluid continues to flow along the hydraulic fracture (Fig. 2e).
In the numerical model, the microseismic points are generated by combining Eq. (4) (Kanamori 1977;Hanks and Kanamori 1979), Eq. (5) (Dershowitz et al. 1996), and Eq. (6) (Kanamori 1977;Dershowitz et al. 1996), as shown below: where G is the shear modulus; A is the shear slip surface area; Δu s is the shear displacement; K n is the fracture's normal stiffness; Δ is the change of shear stress; M 0 is the seismic moment, and M w is the seismic magnitude (instantaneous amplitude).

Hydraulic fracture simulation
As one of the three major geothermal fields that have been identified in the Gonghe Basin, the Qiabuqia geothermal field is located in the eastern part of the Gonghe basin on the Qinghai-Tibet Plateau. With a total area of about 246.9 square kilometers, the Qiabuqia geothermal field is now a very important area with exploration and development potential of hot dry rock geothermal resources in China. In recent years, many studies have been carried out on the geothermal characteristics and enhanced geothermal system of the Qiabuqia geothermal field (e.g., Zhang et al. 2018Zhang et al. , 2019Lei et al. 2019;Shan et al. 2020;Zhu et al. 2022), so there are abundant geothermal data related to the Qiaqia geothermal field available in the existing literature. In this study, taking geothermal field data from Qiabuqia as a geological reference, a series of numerical simulations were carried out to explore the effects of different fracture parameters on the microseisms induced by hydraulic fracturing.

Model setup
To study the effect of faults on hydraulic fracturing, a 3D geothermal reservoir model was first established using Frac-Man software. The parameters used for building the reservoir model are based on the Qiabuqia geothermal field Fig. 2 Activation relationship between artificial and natural fractures in the case that the tip pressure of fracturing fluid is greater than the normal stress acting on the natural fracture surface, where the orange arrows represent the flow direction of the fracturing fluid, the green lines represent natural fractures, the orange lines represent newly created secondary fractures, the triangular shaped grey area surrounded by yellow lines represents the hydraulic fracture, and σ 3 is the minimum principal stress of in-situ stress data in the existing literature (Lei et al. 2019) and the actual field survey data, as listed in Table 1. As shown in Fig. 3, this model is presented in a rectangular shape, with X, Y, and Z dimensions of 6000 m, 5000 m, and 2000 m, respectively, to avoid boundary effects. The total elevation of this reservoir model in the z-dimension is 2000 m, representing a depth from − 2000 to − 4000 m. The buried depth of the reservoir is from 2400 to 3600 m, which is located in the middle of the model. The influence of the overlying rock layers was simulated by applying stress to the upper surface  where the reservoir is actually buried, roughly located at the burial depth of 2400 m to 3600 m, where the color legend represents the distance from the fault and the dark green stripes represent the faults of the model. In the actual situation, there is no fault in the Qiabuqia geothermal reservoir. But to explore the effect of faults on the microseisms induced by hydraulic fracturing, two faults with different geometries were randomly generated and added to the model. In practice, a large number of numerical simulations were performed by adding various faults with randomly generated geometry, and it was found that the simulation results were not affected by the geometry of the fault. These two faults (dark green stripes in Fig. 3b) were distributed in the reservoir model, with a direction of North-North-East and an inclination of about 40°. It is noted that both the geothermal reservoir with faults and the fault-free geothermal reservoir were simulated in this study.
Although it is too ideal to run the simulation for the case without faults because faults usually exist in rock reservoirs, the non-fault reservoir scenario was still simulated for three main purposes: (1) to verify the accuracy and availability of the established model through this simplified situation; (2) indeed, there are no faults in the actual Qiabuqia geothermal reservoir; more importantly, (3) taking the case of the fault-free geothermal reservoir as a reference, the variation of fracture propagation in this case was compared to that in the case with faults to highlight the effect of the existence of faults on induced microseismic events. In this reservoir model (Fig. 3), the size of each grid unit to grid the domain was set to 50 × 50 × 50 m, and a total of 480,000 grid units were obtained for generation after the gridding. These grids representing the rock matrix contain the properties of rock strength and stress field. The discrete fracture network (DFN) was generated based on the fault characteristics. The input parameters used to establish the DFN are partly derived from Lei et al. (2019) and partly from field measurements of the Qiabuqia geothermal field, as listed in Table 2.

In-situ stress state
In fact, the rock matrix is often coupled with fractures. In this study, the selected fracture group was calculated based on the geometry-based Oda method (Wang et al. 2008).
Specifically, the compliance tensor for the fracture was first calculated and then combined with that for the rock matrix. The obtained tensor was inverted to give the corresponding stiffness tensor for the combination of the rock matrix and fracture sets. After that, the stress in the grid was calculated based on various input parameters, such as boundary conditions, material properties, and a strain field output from dynamic upscaling, by using FracStress that can provide a preliminary estimate of the likelihood of fracture slip under different stress conditions (Franzel et al. 2019). In this model, the stress boundary conditions were imposed by applying the stress field to the whole model. The constant gradient stress was also applied to the model, and the bottom of the model should be fixed. Given that the fluid density is 1 000 kg/m 3 and the gravitational acceleration is 9.81 m/s 2 , the main parameters of the strike-slip stress field used in the establishment of the model are shown in Table 3, including the gradients of minimum horizontal stress (S Hmin ), vertical stress (S v ), and maximum horizontal stress (S Hmax ) and fluid pore pressure, which are partly derived from Lei et al. (2019) and partly from field measurements of the Qiabuqia geothermal field.
After assigning these parameter values to the model, including rock matrix parameters (Table 1), fracture parameters (Table 2), and initial in-situ stress state (Table 3), a model specially designed for the Qiabuqia geothermal field was set up. On the basis of this model, the models used for different numerical simulations were established by changing the fracture parameters that may affect the induced microseisms. In practice, the appearance of these different models, derived from the model specially designed for the Qiabuqia geothermal field, did not show any change, similar to Fig. 3b. Therefore, the models used to explore the effects of the existence of faults (section "Effect of faults on induced microseisms"), in-situ stress state (section "Effect of in-situ stress state on induced microseisms"), fracture occurrence (section "Effect of fracture occurrence on induced microseisms"), and fracture distribution near faults (section " Effect of fracture distribution on induced microseisms") on the induced microseisms are no longer presented one by one in this paper.  It should be noted that the effect of the change of the groundwater level is negligible in the establishment of the models, as the groundwater level in the Qiabuqia geothermal field is relatively stable and its variation is considered to have little effect on hydraulic fracturing, although many studies have demonstrated that the groundwater level plays an important role in earthquakes (e.g., Steeples 2005; Grelle and Guadagno 2009).

Hydraulic fracturing simulation
The hydraulic fracturing simulations were carried out based on the critical stress analysis theory, which involves solving the constitutive relation of material mass in the conservative fracture network (Paris and Sih 1965). Hydraulic fracturing was simulated by maintaining a volume balance between the volume of injected water and natural fracture propagation and the volume of developed hydraulic fracture. In this study, the geothermal field data from Qiabuqia was taken as a geological reference for hydraulic fracturing simulations. According to the Qiabuqia geothermal field data in the existing literature (Lei et al. 2019) and the actual field survey data, the input parameters were chosen for hydraulic fracturing simulations, as listed in Table 4. It is noted that both Young's modulus and Poisson's ratio vary with depth. The established model can be divided into three layers: cap rock, reservoir, and substratum. From top to bottom, the values of Young's modulus are 35, 27, and 40 GPa, respectively, and the values of Poisson's ratio are 0.22, 0.24, and 0.21, respectively.
The detailed steps of the whole hydraulic fracturing simulation procedure are described as follows. Firstly, the appropriate model is applied to calculate the change in pore pressure due to water injection, which is usually affected by many factors such as pumping pressure, flow distance, and normal stress. The change of pore pressure is often accompanied by the variation of fracture characteristics, which causes the rock mass to break generating new fractures, and the expansion and propagation of the existing fractures. The fracture process of the rock mass follows the Mohr-Coulomb strength criterion (Eq. (3)). During the process of rock fracture and fracture propagation, the numerical values of the pore pressure and the normal stress of the fracture are compared by setting the model, which directly determines how the fracture extends. The magnitude of microseismic events generated during rock fracture and fracture propagation is calculated according to shear slip surface area, shear displacement, shear stress, etc. (Eqs. (4) ~ (6)). The final simulation results can show the expansion and propagation of existing fractures, the distribution of newly generated fractures, as well as the relevant microseismic events, as a result of hydraulic fracturing.
The simulation results of hydraulic fracturing simulation in the case of the fault-free geothermal reservoir are shown in Fig. 4. It can be found that the microseismic events generated by hydraulic fracturing are mainly concentrated near the well site. A total of 847 microseismic events occurred, among which the farthest microseismic event away from the fractured well occurred at 529 m, and the maximum magnitude of the microseismic event was 0.12. The simulation results are in good agreement with the field survey data, verifying the availability and correctness of the model.
In the simulation, the fractures generated by hydraulic fracturing can be generally divided into three types: main hydraulic fractures, activated fractures, and induced fractures. If hydraulic fracturing does not start from natural fractures, some new fractures are created, namely main hydraulic fractures, which are connected to natural fractures. Activated fractures are created when the fracturing fluid flows out of the main hydraulic fractures and communicates with natural fractures. Induced fractures are created due to the sliding and dislocation induced by activated fractures. In our case, the fractures generated by hydraulic fracturing are mainly activated fractures. In general, the propagation direction of these newly generated fractures is approximately perpendicular to the direction of the minimum principal stress of in-situ stress, σ 3 , (along the direction of the maximum principal stress of in-situ stress, σ 1 ), as shown in Fig. 5, which is consistent with the previous findings presented by Hubbert and Willis (1957).  The data of microseisms induced by hydraulic fracturing were statistically analyzed, and the results were presented in the form of a bar chart, as shown in Fig. 6. In Fig. 6a, the magnitude of microseismic events is generally small, mostly between − 2.0 and − 0.5. The maximum magnitude is 0.12, but the number of microseismic events with magnitudes greater than 0 accounts for less than 3% of the total, so hydraulic fracturing has little effect on microseisms. In Fig. 6b, the distance between microseismic events and injection wells is mostly between 200 and 450 m. The number of microseismic events with distances greater than 500 m only accounts for about 5% of the total, indicating that the influence scope of hydraulic fracturing is very small. All of these demonstrate that in the case of the geothermal reservoir without faults, hydraulic fracturing has little effect on microseismic events, not causing significant harm in actual engineering projects.

Model verification
Although the consistent performance of the hydraulic fracturing simulation results in the case of the fault-free reservoir and the actual field survey data proves the effectiveness of the models established in this study, the validity and accuracy of these models were further verified by comparing the numerical simulation results with the field test results of hydraulic fracturing in Qiabuqia geothermal field. The hydraulic fracturing test was conducted on site in the Qiabuqia geothermal field according to the testing procedures described in Haimson and Cornet (2003) and Courtier et al. (2017). In the field test, the operated fracturing pressure was 28.5 ~ 38.7 MPa, the fracturing duration was 88 min, the initial wellhead pressure was 8.31 MPa, the injection rate was 0.5 ~ 1.3 m 3 /min, and the pump-off pressure was 20.7 MPa. The field reservoir condition was the same as the reservoir model established in this study. The conditions described above for hydraulic fracturing field tests were assigned to the models established in this study for numerical simulation. The validity of the models was then verified by comparing the variation of the wellhead pressure and the induced microseismic events obtained by the field test and numerical simulation. The comparison of the variation of wellhead pressure obtained by the two methods is shown in Fig. 7. It can be seen that the simulation results are in good agreement with the field test results, with a maximum relative error of about 7.3%. During actual hydraulic fracturing, a total of about 160 microseismic events were detected, while a total of 147 microseismic events were statistically generated in the numerical simulation. The maximum magnitude of microseismic events and σ 3 are the maximum and minimum principal stresses of in-situ stress, respectively obtained by the two methods is both not more than 0.5. It can be found that the simulation results are coincident with the field hydraulic fracturing results, with a small relative error. This is consistent with the verification results of similar models established using the FracMan software in Crosta (1997), Johri and Zoback (2013), and Makedonska et al. (2020). Accordingly, the effects of different fracture parameters, such as fracture occurrence and fracture distribution, on the microseisms induced by hydraulic fracturing can be studied by using the models established based on

Effect of faults on induced microseisms
To explore the effect of faults on the microseisms induced by hydraulic fracturing, two faults were added to the reservoir model with other fracturing conditions unchanged. The distribution of these two faults can be seen in Fig. 3. The hydraulic fracturing simulation results in the case of the geothermal reservoir with faults (see Fig. 8) show that a total of 735 microseismic events occurred, among which the maximum magnitude is 0.22 and the maximum distance between the microseismic event and injection well is 682 m. Compared to the hydraulic fracturing simulation in the case of the fault-free geothermal reservoir described above, the microseismic points in this case not only gather near the injection well, but also tend to migrate to the fault, as shown in Fig. 8b. In addition, due to the existence of the fault, the propagation direction of the newly generated fractures is generally shifted toward that of the fault, as shown in Fig. 9, which is different from the case without faults (see Fig. 5). All of these indicate that the existence of faults has a significant influence on fracture propagation.
In Fig. 10a, the magnitude of microseismic events is mainly between 0 and 0.5, and the number of microseismic events with magnitudes greater than 0 accounts for about 85% of the total. The maximum magnitude of microseismic events is 0.22, which is higher than the maximum magnitude of 0.12 in the hydraulic fracturing simulation without fault.
For the case without faults, the magnitude of microseismic events is mostly concentrated at the middle magnitude and tends to the minimum magnitude (see Fig. 6a). However, affected by the fault, the magnitude of microseismic events, in this case, is mostly concentrated around the maximum magnitude. Therefore, it can be considered that the existence of faults increases the magnitude of microseismic events.
As shown in Fig. 10b, the maximum distance is 682 m, which is greater than that of 529 m in the case without faults. The distance between microseismic events and injection wells can be mainly divided into two parts: 50 ~ 200 m and 450 ~ 650 m. The microseismic points in the 50 ~ 200 m part are mainly concentrated near the well, while the microseismic points in the 450 ~ 650 m part mainly represent the offset microseismic events caused by faults. Compared to the microseismic points are mostly between 200 and 450 m in the case without faults (see Fig. 6b), combined with the observations in Fig. 10, it can be considered that the existence of faults can make the distribution of microseismic events more distant.
To sum up, the existence of faults does affect the propagation direction of the newly generated fractures, causing the fractures to extend in the direction of the fault rather than along the direction of the maximum principal stress of in-situ stress, σ 1 , as it would have been in the case without faults. The faults in the geothermal reservoir can increase the magnitude of microseismic events and make the distribution of microseismic events farther, which is probably because the existence of faults increases the shear fracture in the process of rock fracture, resulting in greater fracture energy and larger microseismic events (Shan et al. 2020).

Effect of in-situ stress state on induced microseisms
In general, the in-situ stress is represented by three orthogonal principal stresses, namely vertical stress (σ V ), maximum horizontal stress (σ H ), and minimum horizontal stress (σ h ). Based on Anderson's fault theory (Anderson 1951), the in-situ stress state can be divided into three types by comparing the relationship between these three principal stresses (Zhang et al. 2018): (1) normal faulting stress state (NF): the vertical stress is the greatest principal stress, i.e., σ V ≥ σ H ≥ σ h . In this stress state, the vertical stress drives normal faulting, and fault slip occurs when the minimum stress reaches a sufficiently low value; (2) strike-slip faulting stress state (SS): the vertical stress is the intermediate principal stress, i.e., σ H ≥ σ V ≥ σ h ; (3) reverse faulting stress state (RF): the vertical stress is the least principal stress, σ H ≥ σ h ≥ σ V .
In this study, all three different stress states were applied to the model to explore the effects of the in-situ stress state on the microseisms induced by hydraulic fracturing. The stress gradients under these three stress states used for simulation are given in Table 5. The hydraulic fracturing  Fig. 11, Fig. 12, and Fig. 13, respectively. It is found that there are 64 microseisms under NF, 735 microseisms under SS, and 1292 microseisms under RF. Obviously, the microseismic events occur the most under reverse faulting stress state (RF).
The magnitudes of microseismic events induced by hydraulic fracturing under NF, SS, and RF stress states are given in Figs. 14a, 15a, and 16a, respectively. It can be seen that the maximum magnitude of microseismic events under NF is 0.23, and there are 28 microseismic events with magnitudes greater than 0, accounting for 43.75% of the total. The maximum magnitude of microseismic events under SS is 0.21, and there are 623 microseismic events with magnitudes greater than 0, accounting for 84.76% of the total. The maximum magnitude of microseismic events under RF is 0.1, and there are 766 microseismic events with magnitudes greater than 0, accounting for 59.29% of the total. Figures 13b, 14b, and 15b, respectively, show the distance between microseismic events and injection wells under NF, SS, and RF. Under NF, all the microseismic events are relatively close to the injection wells, and the maximum distance is only 150 m. Under RF, the distance between microseismic events and injection wells is the farthest, which can reach 1400 m, followed by the maximum distance of 700 m under SS.
For natural tectonic earthquakes, the frequency and magnitude of earthquakes generated by the reverse fault are generally high (Ellsworth 2013;Coban and Sayil 2020). When the reverse fault is subjected to horizontal compression, there is an increase in the maximum static friction (Papazachos et al. 2004). When the fault plane is very close to vertical, the occurrence of relative sliding becomes very difficult (Papazachos et al. 2004). However, once the fault is strongly disturbed, relative sliding might occur at certain points of the fault, and the friction coefficient is reduced by sliding (Sayil 2013). The resulting fracture can gradually extend throughout the fault, accompanied by the release of energy, resulting in large earthquake magnitudes (Crowley and Bommer 2006). In other words, the earthquakes generated by the reverse fault require higher tectonic stress. When the accumulated tectonic stress is large enough, the fault system stores more strain energy at this time, which is more prone to large earthquakes (Crowley and Bommer 2006;Sayil 2013).
For the microseismic induced by hydraulic fracturing, under the premise of artificially applying the same energy required to generate formation fractures, the release forms are different under different stress states (Westwood et al. 2017). Under NF, there are a small number of microseismic events, but these microseismic events are concentrated near the well with relatively large magnitudes. Under RF, there are the most microseismic events. These microseismic events have the widest distribution range and the farthest influence distance, but the magnitudes are relatively small. The number of microseismic events under SS is in the middle of that under three stress states. The magnitude of microseismic events under SS is less than that under NF, Fig. 10 Statistical analysis of (a) microseismic magnitude and (b) the distance between microseismic events and injection wells, obtained from hydraulic fracturing simulation results in the case of the geothermal reservoir with faults and the distribution range and influence distance under SS are smaller than those under RF. In general, the total energy of microseisms under RF is higher than that under NF and SS, which is also consistent with the fact that earthquakes generated by reverse faults in natural tectonic earthquakes release more energy (Kelsey et al. 2008;Yang et al. 2021).
It should be noted that in the actual simulation, the hydraulic fracturing simulation was performed for only 24 h based on the established models. Due to the short simulation time, the accumulated energy was not enough to generate a large earthquake (microseismic). In addition, there are some differences between natural tectonic earthquakes and those induced by hydraulic fracturing. Tectonic earthquakes are caused entirely by tectonic stress (Sayil 2013). While for the earthquakes induced by hydraulic fracturing in this study, they are microseisms caused by rock fracture as a result of the change of internal rock pressure due to water injection. These two points result in a slight difference between the simulation results in this study and the actual situation of natural tectonic earthquakes described above. However, the simulation results of this study show that under RF, the total number of induced microseisms is much more than that under the other two stress conditions (NF and SS), and the total energy released (the sum of all the microseismic energy) is the largest, far greater than the other two stress conditions (NF and SS). Therefore, it can be considered that the simulation results are consistent with the actual occurrence of natural tectonic earthquakes.

Effect of fracture occurrence on induced microseisms
During hydraulic fracturing, fracture occurrence plays a very important role in the performance of hydraulic fracturing (Dehghan et al. 2015; Barthwal and van der Baan 2019). Through a series of laboratory rock tests, Dehghan et al. (2015) pointed out that the hydraulic fracture propagation behavior and geometry were affected by natural fracture inclination and strike. Based on this, the effects of fracture occurrence on the microseisms induced by hydraulic fracturing were studied by changing the dip angle of fractures and the positional relationship between fractures and faults.
The dip angles of the fracture used in this study include 0°, 30°, and 60°. The position relationship between fractures and faults can be roughly divided into three types: the parallel relationship between the fracture and the fault, the vertical relationship between the fracture and the fault, and the included angle between the fracture and the fault is 45°. Therefore, a total of nine different cases were used to study the effect of fracture occurrence on the induced microseisms. As shown in Fig. 17, the fracture joint strike in nine different cases is illustrated in the form of a combination of stereonet and rose diagram, where the deep red fan inside indicates the strike of natural fractures and the concentric circles indicate the dip angle of natural fractures, ranging from 0 to 90° from inside to outside. The legend on the right of each subfigure represents the density of natural fractures, and the numbers 1 ~ 9 next to the legend indicate that the density increases gradually. For each subfigure, from left to right, the dip angle of natural fractures are 0°, 30°, and 60°, respectively. Figure 17a, b, and c represents the three position relationships between fractures and faults described above: parallel, vertical, and the included angle of 45°. Figure 18 shows the variation of the number, the distribution range, and the maximum and average magnitude of the microseisms induced by hydraulic fracturing. It can be seen that the fracture occurrence has a significant influence on the microseisms induced by hydraulic fracturing. In Fig. 18a, when the fracture dip angle increases from 0° to 30° and then to 60°, the corresponding number of microseisms shows an upward trend. When the fracture is perpendicular to the fault, the number of microseisms is the least among the three positional relationships, regardless of the dip angle. As for the distribution range of microseisms, when the fracture dip angle increases from 0° to 30° and then to 60°, the corresponding farthest distance of the microseismic event decreases gradually (see Fig. 18b). When the included angle between the fracture and the fault is 45°, the microseismic farthest distribution distance is always the largest among the three positional relationships, with the same fracture dip angle. When the fracture is perpendicular to the fault, the microseismic farthest distribution distance is always the smallest among the three positional relationships, with the same fracture dip angle. Figure 18c shows that in all nine cases, the maximum magnitude of induced microseisms varies from − 0.5 to 0.5, suggesting that the fracture occurrence has little effect on the maximum magnitude of microseisms. However, the average magnitude varies greatly under different fracture dip angles and different positional relationships between fractures and faults, as shown in Fig. 18d. The maximum value of the average magnitude is about 0 when the fracture dip angle is 0° and the included angle between the fracture and the fault is 45°. When the fracture dip angle is 30° and fractures and faults are perpendicular to each other, the average magnitude reaches a minimum of about − 1.6. It is also found that when the fracture dip angle is 30°, the average magnitude is significantly smaller than that under the other two angles (0° and 60°). Even so, the maximum difference in the average magnitude in these nine cases is only about 1.5. Therefore, from the nine cases in this study, it can be considered that the effect of fracture occurrence on microseismic magnitude is not very obvious, especially on the maximum magnitude.

Effect of fracture distribution on induced microseisms
According to the possible aggregation of fractures near the fault in the actual situation, four different fracture distributions, A, B, C, and D, with the same number of fractures, were set in this study, from uniform to dense distribution, as shown in Fig. 19. Each subfigure in Fig. 19 represents a top view of the two faults and the distribution of fractures around them, where the green strip represents the fault, and the remaining scatters represent the fractures. A represents the uniform distribution of fractures, and D Fig. 15 Statistical analysis of (a) microseismic magnitude and (b) the distance between microseismic events and injection wells under strike-slip faulting stress state (SS) represents the fractures mainly distributed near the fault. From A to D, the fracture gradually concentrates on the fault. The effects of fracture distribution near the fault on induced microseisms were studied by comparing the characteristics of induced microseisms under these four fracture distributions.
As shown in Fig. 20a, the fracture distribution near faults has a great influence on the maximum magnitude and the farthest distribution distance of induced microseisms. With the same number of fractures, the maximum magnitude of induced microseisms increases as the fractures are more concentrated near the fault, although the maximum magnitude  Figure 20b shows that the variation trend of the average magnitude and average distance of microseismic events is similar to that of Under the fracture distribution D, the fractures are concentrated near the fault, and the induced microseismic events have the farthest distribution distance and the largest magnitude. This can be explained by the fact that during hydraulic fracturing, the fracturing fluid continues to expand outward, and after activating the natural fractures, some new fractures continue to be created, leading to more microseismic events. However, under the fracture distribution A, the fractures are almost uniformly distributed, and the fracture network enhances the hydraulic connectivity in the deep underground. During hydraulic fracturing, the fracturing fluid only needs to communicate with the natural fracture network. The effect of improving hydraulic connectivity can be achieved without or with few new fractures. Therefore, under the fracture distribution A, the microseismic events induced by hydraulic fracturing have the closest distribution distance and the smallest magnitude.

Conclusions
In this study, taking the Qiabuqia geothermal field data as a reference, the models were set up for different numerical simulations based on the built-in models in the commercial FracMan™ software suite. Through these established models, a series of numerical simulations were carried out to investigate the effects of different fracture parameters, including the existence of faults, in-situ stress state, fracture Under four different fracture distributions: a maximum magnitude and farthest distribution distance of microseismic events; b average magnitude and an average distance of microseismic events occurrence, and fracture distribution near faults, on the magnitude and distribution range of the microseisms induced by hydraulic fracturing. According to the numerical simulation results under different variables, the following conclusions are derived from this work: 1. Compared with the case of the fault-free geothermal reservoir, the existence of faults does affect the propagation direction of the newly generated fractures, causing the fractures to extend toward the fault. The faults in the geothermal reservoir increase the magnitude of microseismic events and make the distribution of microseismic events farther, which is probably because the existence of faults increases the shear fracture in the process of rock fracture, producing greater fracture energy. 2. Under three different in-situ stress states of NF, SS, and RF, the magnitude of induced microseisms and the distance between microseismic events and injection wells are different, mainly due to the different forms of energy release. The total energy of microseisms under RF is generally higher than that under NF and SS. Under RF, there are the most microseismic events, with a total of 1292, and these microseismic events have the widest distribution range and the farthest influence distance, but the magnitudes are relatively small, with a maximum magnitude of 0.1. Under NF, the number of microseismic events is the least, with only 28, and these microseismic events are concentrated near the well, with the farthest distribution distance of only 150 m, but there is the largest magnitude of microseismic events of 0.23. 3. Based on nine different cases in this study, the fracture occurrence, in terms of the fracture dip angle and the positional relationship between fractures and faults, can be considered to have a significant effect on the microseisms induced by hydraulic fracturing. When the fracture dip angle increases from 0 to 30° and then to 60°, the corresponding number of induced microseisms gradually increases, but the corresponding farthest distribution distance of microseismic events gradually decreases. With the same fracture dip angle, the microseismic farthest distribution distance is always the largest when the included angle between the fracture and the fault is 45°, and the microseismic farthest distribution distance is always the smallest when the fracture is perpendicular to the fault. The effect of fracture occurrence on microseismic magnitude is not very obvious, especially on the maximum magnitude. 4. With the same number of fractures, as the fracture concentration to the fault increases, both the maximum magnitude and the farthest distribution distance of the induced microseisms induced increase. 5. In the actual geothermal system site selection and reservoir reconstruction optimization, the fault, the area of RF stress state, and the excessive concentration of natural fractures should be avoided as far as possible. In this way, sufficient fracture networks generated by hydraulic fracturing can be satisfied to ensure the optimization of reservoir reconstruction, and meanwhile, the excessive magnitude of induced microseisms can be avoided by maintaining the magnitude of microseismic events within a controllable range.

Data Availability
The data sets generated during and/or analysed during the current study are available from the corresponding author on reasonable request.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http:// creat iveco mmons. org/ licen ses/ by/4. 0/.