Deformation Intensity Distribution and Zoning of Seismic Deformation Fields in Near-Fault Areas

: To describe the distribution law of the degree of stratigraphic deformation in the cross-fault region under the inﬂuence of fault dislocation, the concept of deformation intensity is introduced, and a computational model of the coseismic deformation ﬁeld in the cross-fault region is established based on the dislocation model theory and the rich fault inversion parameters in the Sichuan–Tibet region. According to the deformation law of a typical earthquake in Yushu, key sections and parameters are selected for analysis. The results show that the maximum dislocation of the deformation ﬁeld in the cross-fault region decays from the fault rupture to the surface, and the deformation inﬂuence extends from the upper edge of the fault to the surface, with an overall funnel-shaped symmetric distribution. The deformation intensity is elliptically distributed in the vertical proﬁle. The maximum deformation of the strata across the fault within the burial depth of the structure is ﬁtted to ﬁll the lack of the existing statistical equation of coseismic deformation without considering the burial depth factor. Finally, the strata are zoned by deformation intensity according to the degree of deformation to provide a reference for the seismic protection of structures in the near-fault area. Author Contributions: Conceptualization, formal analysis, investigation, curation, preparation, visualization, supervision, B.Z.; funding acquisition, B.Z.


Introduction
The coseismic deformation field is caused by fault dislocation, distributed in the area near the fault, and is the dominant factor causing the seismic response of underground engineering structures in the near-fault region. Wells and Coppersmith [1] gave the surface displacement-earthquake magnitude fitting relationship through the statistics and analysis of surface rupture zone observation data worldwide, which was used as a basis for assessing the maximum deformation of the coseismic deformation field at different earthquake magnitudes. J. Cheng [2] described the relationship between rupture and earthquake magnitude in China and compared it with the worldwide pattern. Further, Ran H.L. [3] presented the surface deformation data for western China.
Underground projects such as tunnels play an important role as lifeline projects in daily production life and post-earthquake disaster relief. The existing studies mainly focus on the surface area of the co-earthquake deformation field. The ground deformation caused by the occultation rupture [4] extends from the fault surface, and the difference between surface and subsurface deformation becomes non-negligible with increasing depth. Callisto [5] found that the same significant seismic damage exists in deeply buried tunnels in a statistical study of seismic damage. A similar typical situation was also found in the Longxi tunnel [6] in the Wenchuan earthquake, which was a typical deeply buried tunnel that also suffered severe damage during the earthquake. In a study of the Kumamoto earthquake, Zhang X. [7] found a correlation between seismic damage and surface coseismic deformation in the Tawarayama tunnel. Huang, Lu Y. [8] studied the significance of coseismic displacement for tunnel design. K. Ma [9] found that the Symmetry 2022, 14, 239 2 of 18 stability of the surrounding rock around the underground structure is closely related to the distance between the structure and the fault. Shrestha Rajyaswori [10] found that the factors influencing the seismic damage of tunnels in the near-fault region are fault distance, earthquake magnitude, and burial depth of the structure. With the increase in the burial depth of engineering structures in near-fault regions, the existing knowledge of the depth of the coseismic deformation field lags behind the current engineering practice, and it is necessary to extend the study of the coseismic deformation field from the surface to the deeply buried strata in the burial depth range of underground engineering.
For the coseismic deformation field inside the strata, scholars have researched from the perspectives of experiments, numerical simulations, and theories, etc. Anastasopoulos [11] studied the surface deformation and rupture law caused by fault dislocation based on model tests; Ashtiani [12] studied the effect of fault dislocation on building foundation by centrifugal tests. Model tests are limited by the research tools and cannot consider both structural scale and fault scale. Loukidis [13] studied the process of fault rupture propagation in sandy soil by numerical simulation. Baziar [14] evaluated the effect of fault dislocation on tunnel structure by a numerical model. Li H. [15] studied the propagation law of fault dislocations and the influencing factors of surface deformation based on the finite element method. Compared with model tests, numerical methods can study a larger scale area, but the large computational volume prevents them from being able to have a sufficiently fine simulation of the strata at the structural scale while considering faults. At present, the research on the interior of the strata is mainly focused on the fault scale, and the research on the deformation field of the strata within the burial depth of the structure at the engineering scale is rarely reported.
By studying the distribution of the deformation degree of the deformation field inside the strata, it is of great significance to quickly and accurately evaluate the ground shaking damage force and guide the post-earthquake emergency rescue. Most of the current seismic intensities are classified by the PGA of ground shaking, while there is no relevant index for the evaluation of the degree of stratigraphic dislocation. Zhang C. J. [16] found that the surface ground shaking intensity has a strong correlation with the earthquake coseismic deformation field. The existing code [17] has only a rough suggested value for the avoidance of the fault-affected zone for above-ground structures, and there is no clear method and standard for the consideration of anti-staggering for underground works in the near-fault area. Concerning the concept of ground vibration intensity distribution at the surface, it is necessary to introduce a corresponding concept to quantitatively describe the distribution of the ground deformation degree inside the strata, defining this indicator as to the deformation intensity.
Because of the large spatial location difference between the stratum in which the deeply buried tunnel is located and the surface and the need to recognize the influence of the deformation field inside the stratum on the underground engineering, this paper investigates the intensity of stratum deformation inside the stratum. Firstly, based on the fault parameters obtained from the inversion of the typical Yushu earthquake in the Sichuan-Tibet region of China, the coseismic deformation field is calculated orthogonally in combination with the stratigraphic model of the Sichuan-Tibet region. The calculated results are compared with the surface observed deformation fields to verify the accuracy of the calculated results and the applicability of the adopted method. Based on the validated model, the deformation distribution law of the internal deformation field of the strata is analyzed according to the deformation propagation law with the surface deformation field and fault parameters as constraints. Typical monitoring locations are selected to statistically regress the equation of the relationship between the maximum deformation of the strata and the seismic intensity within a certain depth underground. By introducing the concept of deformation intensity and referring to the relevant structural deformation index, the intensity division of the subsurface deformation field is performed to form the intra-stratigraphic deformation intensity partition. The results obtained in this paper are intended to provide a basis for future construction of underground structures, such as tunnels near faults, and to provide a useful reference for rapid damage assessment and rescue planning in the relevant areas after earthquakes.

Calculation Method Based on Dislocation Model
The coseismic displacement field caused by fault dislocation can be described by the dislocation theory. J.A. Steketee et al. [18] were the first to introduce the elastic dislocation theory into the calculation of the earthquake coseismic deformation field, and then scholars worldwide made further developments and refinements on this basis. Y. Okada et al. [19,20] generalized the equation for the stratigraphic deformation caused by faults in the elastic half-space and formed one of the most widely used equations, which is shown in Equation (1) where u i (x 1 , x 2 , x 3 ) denotes the dislocation deformation field generated by the dislocation ∆u j (ξ 1 , ξ 2 , ξ 3 ) over the entire surface Σ in an isotropic medium. Hashima A. et al. [21] studied stratified half-space moderate effect strain fields and demonstrated by comparison that the results based on stratified stratigraphic models are more accurate than those obtained by uniform half-space stratigraphic models. Subsequently, Wang R.J. et al. [22,23] improved the transfer stability of Green's function in stratigraphic strata to control the propagation of computational errors in stratigraphic strata calculations. Zhou J. et al. [24] analyzed the calculation of the internal deformation field of a spherical formation and found that the influence of the earth curvature effect on the formation deformation is small in the near field within 100 km of the fault. Dong J. et al. [25] also found that the effect of geocurvature can be neglected in the local area of the near field. Taking into account the above research results, this paper proposes to use the Green's function based on the horizontal stratigraphy in the near-fault zone at the engineering scale to calculate the co-seismic deformation field inside and at the surface of the strata by combining the fault model and stratigraphic model parameters, and the overall workflow is shown in Figure 1. performed to form the intra-stratigraphic deformation intensity partition. The results obtained in this paper are intended to provide a basis for future construction of underground structures, such as tunnels near faults, and to provide a useful reference for rapid damage assessment and rescue planning in the relevant areas after earthquakes.

Calculation Method Based on Dislocation Model
The coseismic displacement field caused by fault dislocation can be described by the dislocation theory. J.A. Steketee et al. [18] were the first to introduce the elastic dislocation theory into the calculation of the earthquake coseismic deformation field, and then scholars worldwide made further developments and refinements on this basis. Y. Okada et al. [19,20] generalized the equation for the stratigraphic deformation caused by faults in the elastic half-space and formed one of the most widely used equations, which is shown in Equation (1) where ( , , ) denotes the dislocation deformation field generated by the dislocation Δ ( , , ) over the entire surface Σ in an isotropic medium.
Hashima A. et al. [21] studied stratified half-space moderate effect strain fields and demonstrated by comparison that the results based on stratified stratigraphic models are more accurate than those obtained by uniform half-space stratigraphic models. Subsequently, Wang R.J. et al. [22,23] improved the transfer stability of Green's function in stratigraphic strata to control the propagation of computational errors in stratigraphic strata calculations. Zhou J. et al. [24] analyzed the calculation of the internal deformation field of a spherical formation and found that the influence of the earth curvature effect on the formation deformation is small in the near field within 100 km of the fault. Dong J. et al. [25] also found that the effect of geocurvature can be neglected in the local area of the near field. Taking into account the above research results, this paper proposes to use the Green's function based on the horizontal stratigraphy in the near-fault zone at the engineering scale to calculate the co-seismic deformation field inside and at the surface of the strata by combining the fault model and stratigraphic model parameters, and the overall workflow is shown in Figure 1.

Typical Earthquake Example Verification
At 7:49 a.m. on 14 April 2010, an earthquake of magnitude 7.1 on the Richter scale occurred in Yushu County, Yushu Tibetan Autonomous Prefecture, Qinghai Province, with a typical left-rotation almost purely strike-slip fault. The intensity of the epicenter reached nine degrees, causing a large number of casualties and property damage and was another earthquake that triggered serious disasters in China after the Wenchuan earthquake.

Typical Earthquake Example Verification
At 7:49 a.m. on 14 April 2010, an earthquake of magnitude 7.1 on the Richter scale occurred in Yushu County, Yushu Tibetan Autonomous Prefecture, Qinghai Province, with a typical left-rotation almost purely strike-slip fault. The intensity of the epicenter reached nine degrees, causing a large number of casualties and property damage and was another earthquake that triggered serious disasters in China after the Wenchuan earthquake.
With the progress of measurement technology, geodetic data, such as GPS and INSAR, have a high surface coverage, which provides good conditions for studying the surface coseismic deformation field. Wang Q.C. [26] found that the Yushu M S 7.3 earthquakes occurred on the Ganzi-Yushu fault zone, and the source fault was a north-west-trending high-angle fault. Shan B. [27] simplified the finite fault model of the Yushu earthquake by ignoring the bending of the fault rupture surface in the northwest section, but the simplification of the model had less effect on the results because the coseismic dislocations of the Yushu earthquake were mainly located in the southeast section of the rupture surface. The above studies and the corresponding fault inversion data [28] were synthesized, and the average slip of the Yushu earthquake-prone fault was determined to be 1.1 m. The specific parameters are shown in Table 1. To obtain a more accurate inversion, detailed stratigraphic parameters were collected for the Yushu area. Wu J. P. [29] gave a model of the crustal velocity structure in the Sichuan and Yunnan regions in western China, which is shown in Table 2. Based on the Yushu fault and stratigraphic parameters, the amount of stratigraphic deformation and its distribution at different depths on the surface and within the strata are calculated, and the results are shown in Figure 2. With the progress of measurement technology, geodetic data, such as GPS and INSAR, have a high surface coverage, which provides good conditions for studying the surface coseismic deformation field. Wang Q.C. [26] found that the Yushu MS7. 3 earthquakes occurred on the Ganzi-Yushu fault zone, and the source fault was a northwest-trending high-angle fault. Shan B. [27] simplified the finite fault model of the Yushu earthquake by ignoring the bending of the fault rupture surface in the northwest section, but the simplification of the model had less effect on the results because the coseismic dislocations of the Yushu earthquake were mainly located in the southeast section of the rupture surface. The above studies and the corresponding fault inversion data [28] were synthesized, and the average slip of the Yushu earthquake-prone fault was determined to be 1.1 m. The specific parameters are shown in Table 1. To obtain a more accurate inversion, detailed stratigraphic parameters were collected for the Yushu area. Wu J. P. [29] gave a model of the crustal velocity structure in the Sichuan and Yunnan regions in western China, which is shown in Table 2. Based on the Yushu fault and stratigraphic parameters, the amount of stratigraphic deformation and its distribution at different depths on the surface and within the strata are calculated, and the results are shown in Figure 2.   To verify the accuracy of the calculation results, previous observations of the Yushu earthquake in the post-earthquake period were collected. Ma Y.S. [30] found that the surface rupture in the south-central section was around 30 cm in the observation of the surface deformation of the Yushu earthquake. Xun Wang [31] observed the deformation of the surface isoseismic displacement map generated by InSAR to be around 0.3m, with the maximum deformation occurring in the upper disk, which is close to the simulation results. Zhou H.F. [32] also concluded that the surface deformation was around 30 cm. After comparison, it is found that the surface deformation calculated by the dislocation model is close to the results obtained from the surface observation, which verifies the feasibility of calculating the surface coseismic deformation by this method.
It can be seen from Figure 2 that the deformation on both sides of the stratum is smaller near the fault center because the coordinates are established at the origin. The stratigraphy on both sides of the fault is asymmetrically distributed about the upper edge of the fault, with a slight upper plate effect. From the comparison of the deformation field at different depths, it can be seen that the deformation volume increases with the increase of burial depth, while the deformation influence range decreases. The morphology of the overall deformation field is funnel-shaped.

Analysis of Coseismic Deformation Field of Sichuan-Tibet Strata
The amount of coseismic deformation of surface strata is an intuitive indicator for observing seismic deformation, while in the internal regions of strata, structures are affected by stratigraphic dislocation deformation. To accurately describe the affected degree of different regions within the stratum, the indicator of relative stratigraphic dislocation deformation is introduced. The deformation degree of the stratum is defined as the relative deformation between different points of the stratum in the threedimensional deformation field, taking two adjacent points i and i + 1 inside the stratum, and the deformation difference is defined in Figure 3. To verify the accuracy of the calculation results, previous observations of the Yushu earthquake in the post-earthquake period were collected. Ma Y.S. [30] found that the surface rupture in the south-central section was around 30 cm in the observation of the surface deformation of the Yushu earthquake. Xun Wang [31] observed the deformation of the surface isoseismic displacement map generated by InSAR to be around 0.3 m, with the maximum deformation occurring in the upper disk, which is close to the simulation results. Zhou H.F. [32] also concluded that the surface deformation was around 30 cm. After comparison, it is found that the surface deformation calculated by the dislocation model is close to the results obtained from the surface observation, which verifies the feasibility of calculating the surface coseismic deformation by this method.
It can be seen from Figure 2 that the deformation on both sides of the stratum is smaller near the fault center because the coordinates are established at the origin. The stratigraphy on both sides of the fault is asymmetrically distributed about the upper edge of the fault, with a slight upper plate effect. From the comparison of the deformation field at different depths, it can be seen that the deformation volume increases with the increase of burial depth, while the deformation influence range decreases. The morphology of the overall deformation field is funnel-shaped.

Analysis of Coseismic Deformation Field of Sichuan-Tibet Strata
The amount of coseismic deformation of surface strata is an intuitive indicator for observing seismic deformation, while in the internal regions of strata, structures are affected by stratigraphic dislocation deformation. To accurately describe the affected degree of different regions within the stratum, the indicator of relative stratigraphic dislocation deformation is introduced. The deformation degree of the stratum is defined as the relative deformation between different points of the stratum in the three-dimensional deformation field, taking two adjacent points i and i + 1 inside the stratum, and the deformation difference is defined in Figure 3.
Using the deformation indexes defined in Figure 3, the deformation difference cloud maps for different depth strata are calculated and the general shape of the subsurface deformation intensity field distribution is obtained. To better show the misalignment in the middle region of the results, logarithmic coordinates are used to display the deformation cloud map, see Figure 4.  Using the deformation indexes defined in Figure 3, the deformation difference cloud maps for different depth strata are calculated and the general shape of the subsurface deformation intensity field distribution is obtained. To better show the misalignment in the middle region of the results, logarithmic coordinates are used to display the deformation cloud map, see Figure 4.   Using the deformation indexes defined in Figure 3, the deformation difference cloud maps for different depth strata are calculated and the general shape of the subsurface deformation intensity field distribution is obtained. To better show the misalignment in the middle region of the results, logarithmic coordinates are used to display the deformation cloud map, see  Compared with surface deformation, the distribution of subsurface deformation is not easy to observe and has an important impact on the seismic protection of structures in this area. The ground deformation caused by the fault is mainly concentrated in the middle of the fault, and the displacement at the two ends of the fault is almost zero. From the most unfavorable situation, the monitoring profile is set to the vertical tangent of the fault. It can be observed from Figure 4 that the width range of deformation influence decreases gradually from shallow to deep, and the change rate is non-uniform. In addition, the variation of deformation intensity along the fault strike direction is small and almost uniformly distributed. Therefore, for simplification, we are mainly concerned with the variation of the deformation degree of the strata along the vertical and fault strike direction and depth direction, ignoring the variation of the deformation field along the fault strike direction, and simplifying the three-dimensional deformation field to two dimensional. The horizontal range of the study is limited to the near-field engineering scale within 60 km on both sides of the fault, and the depth range is limited to the stratigraphic region between the upper edge of the fault and the ground surface. The deformation characteristics of the stratigraphic deformation field are clarified through the analysis of the typical Yushu earthquake in the Sichuan-Tibet region. To further analyze the law of stratigraphic deformation in depth, monitoring profiles are established through the stratigraphic deformation characteristics clarified above, as shown in Figure 5. Compared with surface deformation, the distribution of subsurface deformation is not easy to observe and has an important impact on the seismic protection of structures in this area. The ground deformation caused by the fault is mainly concentrated in the middle of the fault, and the displacement at the two ends of the fault is almost zero. From the most unfavorable situation, the monitoring profile is set to the vertical tangent of the fault. It can be observed from Figure 4 that the width range of deformation influence decreases gradually from shallow to deep, and the change rate is non-uniform. In addition, the variation of deformation intensity along the fault strike direction is small and almost uniformly distributed. Therefore, for simplification, we are mainly concerned with the variation of the deformation degree of the strata along the vertical and fault strike direction and depth direction, ignoring the variation of the deformation field along the fault strike direction, and simplifying the three-dimensional deformation field to two dimensional. The horizontal range of the study is limited to the near-field engineering scale within 60 km on both sides of the fault, and the depth range is limited to the stratigraphic region between the upper edge of the fault and the ground surface. The deformation characteristics of the stratigraphic deformation field are clarified through the analysis of the typical Yushu earthquake in the Sichuan-Tibet region. To further analyze the law of stratigraphic deformation in depth, monitoring profiles are established through the stratigraphic deformation characteristics clarified above, as shown in Figure 5. Based on the analytical laws obtained in the Yushu earthquake, parametric modeling is performed by combining the range of parameters collected in the Sichuan-Tibet region.
Using the basic stratigraphic model, fault parameters, and dislocation modeling methods, the coseismic deformation field of the strata at the time of the earthquake is calculated, and the index of the deformation degree of the coseismic deformation field is used to estimate the intra-stratigraphic deformation intensity. A sensitivity analysis of the parameters is also carried out to obtain general patterns in the Sichuan-Tibet region. The current and more accurate global stratigraphic model Crust1.0 [33] is selected, which has a high degree of confidence and applicability. This results in a more general stratigraphic model of formation, as shown in Table 3. Based on the analytical laws obtained in the Yushu earthquake, parametric modeling is performed by combining the range of parameters collected in the Sichuan-Tibet region.
Using the basic stratigraphic model, fault parameters, and dislocation modeling methods, the coseismic deformation field of the strata at the time of the earthquake is calculated, and the index of the deformation degree of the coseismic deformation field is used to estimate the intra-stratigraphic deformation intensity. A sensitivity analysis of the parameters is also carried out to obtain general patterns in the Sichuan-Tibet region. The current and more accurate global stratigraphic model Crust1.0 [33] is selected, which has a high degree of confidence and applicability. This results in a more general stratigraphic model of formation, as shown in Table 3. Based on the Yushu earthquake, the magnitude was kept the same as that of the Yushu earthquake, and the study was carried out according to the magnitude seven earthquakes. To further determine the characteristics of stratigraphic deformation at different stratigraphic depths, a typical depth of 2000 m is selected for parametric sensitivity analysis of fault parameters. A single variable was controlled to analyze the effects of the sliding angle, dip angle, and stratigraphic depth on stratigraphic deformation perpendicular to the fault strike direction. In a survey of faults in the Sichuan-Tibet region, it was found that strike-slip faults account for 40% of the seismogenic faults in the Sichuan-Tibet region [34]. Therefore, taking the strike-slip fault as an example for illustration, the deformation rates of strata at different depths under a magnitude 7 earthquake were analyzed, and the calculated results are shown in Figure 6. Based on the Yushu earthquake, the magnitude was kept the same as that of the Yushu earthquake, and the study was carried out according to the magnitude seven earthquakes. To further determine the characteristics of stratigraphic deformation at different stratigraphic depths, a typical depth of 2000 m is selected for parametric sensitivity analysis of fault parameters. A single variable was controlled to analyze the effects of the sliding angle, dip angle, and stratigraphic depth on stratigraphic deformation perpendicular to the fault strike direction. In a survey of faults in the Sichuan-Tibet region, it was found that strike-slip faults account for 40% of the seismogenic faults in the Sichuan-Tibet region [34]. Therefore, taking the strike-slip fault as an example for illustration, the deformation rates of strata at different depths under a magnitude 7 earthquake were analyzed, and the calculated results are shown in Figure 6. It can be seen from Figure 6 that the effect of stratigraphic depth on the fault deformation is more significant, and therefore the stratigraphic depth is listed as the main consideration. Based on previous analytical experience, it is known that the dip angle of the fault produces a hanging wall effect, which affects the degree of symmetry of the distribution of deformation. Further analysis of the influence on the degree of stratigraphic deformation when the dip angle of the strike-slip fault varies from 15° to 90° is shown in Figure 7. It can be seen from Figure 6 that the effect of stratigraphic depth on the fault deformation is more significant, and therefore the stratigraphic depth is listed as the main consideration. Based on previous analytical experience, it is known that the dip angle of the fault produces a hanging wall effect, which affects the degree of symmetry of the distribution of deformation. Further analysis of the influence on the degree of stratigraphic deformation when the dip angle of the strike-slip fault varies from 15 • to 90 • is shown in Figure 7.
It can be observed from Figure 7 that there is a slight hanging wall effect of stratigraphic deformation, as manifested by the instability of the values near the hanging wall because the dip angle of the Yushu earthquake fault is 80 • . From the deformation rate of the deeply buried strata, the influence range and deformation values of the strata deformation under different dip angles do not differ much. Therefore, for the sake of simplification, the dip angle of the fault is not lowered as the main consideration for the distribution of the intensity of stratigraphic deformation in the further analysis of stratigraphic deformation. It can be observed from Figure 7 that there is a slight hanging wall effect of stratigraphic deformation, as manifested by the instability of the values near the hanging wall because the dip angle of the Yushu earthquake fault is 80°. From the deformation rate of the deeply buried strata, the influence range and deformation values of the strata deformation under different dip angles do not differ much. Therefore, for the sake of simplification, the dip angle of the fault is not lowered as the main consideration for the distribution of the intensity of stratigraphic deformation in the further analysis of stratigraphic deformation.

Study of the Distribution Law of Deformation Field
Through the previous analysis, how to predict the possible dislocation deformation of faults before the occurrence of earthquakes and delineate the "intensity" of stratigraphic deformation within the strata is an important guideline for the seismic protection of underground structures in near-fault areas. Among the fault parameters, the amount of sliding is, therefore, the main influencing factor of stratigraphic dislocation deformation. To carry out the ground deformation analysis, it is necessary to clarify the relationship between the earthquake level and the fault sliding amount first. According to Aki K. and Richards P.G. [35], the relationship between the average sliding volume u on the fault surface and the seismic moment M is obtained statistically as in Equation (2) M = μAu (2) where u is the shear modulus, and A is the rupture area of the fault. The relationship between seismic moment M and moment magnitude M , obtained statistically by Hanks T.C. and Kanamori H. [36], is given in Equation (3) The unit of seismic moment M in Equation (3) is N-m, and the relationship between the average slip on the fault surface and the matrix level is obtained as in Equation (4

Study of the Distribution Law of Deformation Field
Through the previous analysis, how to predict the possible dislocation deformation of faults before the occurrence of earthquakes and delineate the "intensity" of stratigraphic deformation within the strata is an important guideline for the seismic protection of underground structures in near-fault areas. Among the fault parameters, the amount of sliding is, therefore, the main influencing factor of stratigraphic dislocation deformation. To carry out the ground deformation analysis, it is necessary to clarify the relationship between the earthquake level and the fault sliding amount first. According to Aki K. and Richards P.G. [35], the relationship between the average sliding volume u on the fault surface and the seismic moment M 0 is obtained statistically as in Equation (2) M 0 = µAu (2) where u is the shear modulus, and A is the rupture area of the fault. The relationship between seismic moment M 0 and moment magnitude M w , obtained statistically by Hanks T.C. and Kanamori H. [36], is given in Equation (3) The unit of seismic moment M 0 in Equation (3) is N-m, and the relationship between the average slip on the fault surface and the matrix level is obtained as in Equation (4) The average sliding volume of the fault surface can be obtained from the rupture area by the determined matrix level. Wells and Coppersmith [1] gave an empirical relationship among the rupture area, rupture length, and matrix level as in Equations (5) Through the above analysis and derivation, the relationship between the amount of fault sliding and the seismic level is given in Table 4. In this section, the amount of dislocation deformation is taken as the main influencing factor of fault dislocation under different earthquake levels, and the influence of fault dislocation on stratigraphic deformation is analyzed, and the variation law of stratigraphic deformation along with depth and the variation law of the influence range of different deformation intensities under different earthquake levels are further quantified.

Equation of Maximum Displacement of the Same Seismic Deformation Field
The maximum deformation of the coseismic deformation field is an important reference indicator for analyzing the influence of stratigraphic deformation on the cross-fault structure. For the attenuation of stratigraphic deformation when crossing the fault directly, only the strata on both sides of the fault as well as the surface rupture values are considered, in fact, there is a large number of faults that are hidden within the strata and do not even form ruptures at the surface [37]. In this case, the attenuation of stratigraphic deformation occurs not only in the horizontal direction perpendicular to the fault strike but also in the depth direction of the strata. At the same time, the existing statistical formulas for stratigraphic deformation mainly consider the earthquake magnitude and count only the maximum displacement deformation at the surface. From the analysis in the previous section, it is clear that the burial depth of the fault is an important influencing factor among the influencing factors regarding the ground deformation. Guo Yao [38] proposed the concept of the burial depth of the earthquake source relative to the thickness of the earth's crust in the analysis of ground stress. Since both the burial depth of the fault and the depth of the stratum in which the project is located influence the degree of deformation of the strata here, in this study, the two are considered together, and the concept of relative fault burial depth is proposed, namely.
In Equation (7), H indicates the burial depth of the upper edge of the fault, and h indicates the burial depth of the structure. When the relative burial depth is 0, it is the surface; when the relative burial depth is 1, it is the depth of burial of the fault. The study area is defined as the stratum between the upper edge of the fault and the surface, In the analysis of the previous section, the parameters of the fault itself are also important influencing factors of the ground deformation. For simplicity, the fault misalignment is considered as the main factor by sensitivity analysis. The projection of the fault midpoint on the surface is used as the origin to establish the coordinate axis. The X-axis is perpendicular to the fault strike direction, and the Y-axis is along the depth direction. The maximum stratigraphic deformation at different depths of the epicenter is collected for each relative burial position under different fault misalignment amounts, and the misalignment deformation of the monitored strata is regressed with the relative burial depth and the change of misalignment amount to form an empirical formula. The ground deformation within 4 km from the fault is taken care of to establish a zoning for the degree of ground deformation. The monitoring location and the shape of the deformation field are shown schematically in Figure 8. direction. The maximum stratigraphic deformation at different depths of the epicenter is collected for each relative burial position under different fault misalignment amounts, and the misalignment deformation of the monitored strata is regressed with the relative burial depth and the change of misalignment amount to form an empirical formula. The ground deformation within 4 km from the fault is taken care of to establish a zoning for the degree of ground deformation. The monitoring location and the shape of the deformation field are shown schematically in Figure 8. The maximum dislocations and rupture zone widths caused by the slip fault at different depths for different earthquake levels = 7.0~8.5 are shown in Table 5. From the data in the table, it can be seen that the maximum dislocation of the ground surface increases with increasing depth at the same seismic level. A multiple regression analysis of the calculated results in Table 5 is fitted to obtain the equation of the relationship between the maximum dislocation of the surface under the action of the strike-slip fault and the earthquake magnitude and the burial depth of the fault, which is shown in Equation (8) = 0.607 + 0.041 − 4.517 R = 0.984 (8) where D is the maximum dislocation of the formation. is the moment magnitude, and is the relative burial depth of the formation. The correlation coefficient is 0.984, which indicates that the formula is a good fit for this law.
From the above analysis, the fault sliding volume and the relative burial depth of strata are selected as the main influencing factors, and the study is set to the cross-fault area of the occult fault. The maximum deformation of strata with different relative burial depths under different earthquake levels was monitored, and the data within this region were statistically obtained from the regression equation. After calculation, it is found that The maximum dislocations and rupture zone widths caused by the slip fault at different depths for different earthquake levels M w = 7.0~8.5 are shown in Table 5. From the data in the table, it can be seen that the maximum dislocation of the ground surface increases with increasing depth at the same seismic level. A multiple regression analysis of the calculated results in Table 5 is fitted to obtain the equation of the relationship between the maximum dislocation of the surface under the action of the strike-slip fault and the earthquake magnitude and the burial depth of the fault, which is shown in Equation (8) lgD = 0.607M W + 0.041δ h − 4.517 R 2 = 0.984 (8) where D is the maximum dislocation of the formation. M w is the moment magnitude, and δ h is the relative burial depth of the formation. The correlation coefficient is 0.984, which indicates that the formula is a good fit for this law. From the above analysis, the fault sliding volume and the relative burial depth of strata are selected as the main influencing factors, and the study is set to the cross-fault area of the occult fault. The maximum deformation of strata with different relative burial depths under different earthquake levels was monitored, and the data within this region were statistically obtained from the regression equation. After calculation, it is found that the difference of maximum deformation between surface and deeply buried strata is 10ˆ0.04 ≈ 110% at the fault burial depth of 2.5 km; i.e., the maximum deformation increases by 10% with increasing depth, which is on the dangerous side. Comparing the maximum deformation at different depths proves the necessity of considering the depth.

Deformation Field Intensity Zoning
Underground structures located within the near-fault region need to consider the attenuation and distribution of the deformation intensity in addition to the influence of the maximum deformation of the strata. In a study by Zhang Y. [39], it was found that the distribution characteristics of the coseismic deformation field of the Wenchuan earthquake were highly consistent with the characteristics of the intensity distribution of that earthquake, and the coseismic deformation could respond to the seismic intensity to some extent. In this section, based on the concept of intensity, a theoretical calculation method is adopted to study the distribution of the degree of stratigraphic deformation. By defining the stratigraphic intensity thresholds and establishing the relationship with the deformation intensity according to the characteristics of the coseismic deformation field, different deformation degrees in different regions within the strata are classified.
For the study of deformation within the strata, there are three directional components, namely the depth direction, the strike of the fault, and the direction perpendicular to the strike of the fault. Here, the three-dimensional deformation field is simplified to a two-dimensional plane deformation field for consideration based on the principle of two dimensions with relatively large change gradients, and the soil displacement field under different earthquake levels, see Figure 9.
by 10% with increasing depth, which is on the dangerous side. Comparing the maximum deformation at different depths proves the necessity of considering the depth.

Deformation Field Intensity Zoning
Underground structures located within the near-fault region need to consider the attenuation and distribution of the deformation intensity in addition to the influence of the maximum deformation of the strata. In a study by Zhang Y. [39], it was found that the distribution characteristics of the coseismic deformation field of the Wenchuan earthquake were highly consistent with the characteristics of the intensity distribution of that earthquake, and the coseismic deformation could respond to the seismic intensity to some extent. In this section, based on the concept of intensity, a theoretical calculation method is adopted to study the distribution of the degree of stratigraphic deformation. By defining the stratigraphic intensity thresholds and establishing the relationship with the deformation intensity according to the characteristics of the coseismic deformation field, different deformation degrees in different regions within the strata are classified.
For the study of deformation within the strata, there are three directional components, namely the depth direction, the strike of the fault, and the direction perpendicular to the strike of the fault. Here, the three-dimensional deformation field is simplified to a two-dimensional plane deformation field for consideration based on the principle of two dimensions with relatively large change gradients, and the soil displacement field under different earthquake levels, see Figure 9.  As can be seen from Figure 9, the main deformation occurred within 1.5 km of the upper edge of the fault and its projection on the surface. Therefore, after the earthquake, the main affected area of the earthquake can be quickly identified. The deformation distribution map obtained from the Yushu earthquake parameter calculation shows that the displacement deformation field is close to a funnel-shaped Symmetric distribution. By observing the deformation cloud map, it can be found that the degree of stratigraphic deformation intensifies as the horizontal center and burial depth of the fault are approached, so the two dimensions of depth and epicenter distance are taken as the main considerations for stratigraphic zoning, and the influence law of both on stratigraphic deformation is shown in Figure 10.
upper edge of the fault and its projection on the surface. Therefore, after the earthquake, the main affected area of the earthquake can be quickly identified. The deformation distribution map obtained from the Yushu earthquake parameter calculation shows that the displacement deformation field is close to a funnel-shaped Symmetric distribution. By observing the deformation cloud map, it can be found that the degree of stratigraphic deformation intensifies as the horizontal center and burial depth of the fault are approached, so the two dimensions of depth and epicenter distance are taken as the main considerations for stratigraphic zoning, and the influence law of both on stratigraphic deformation is shown in Figure 10. The deformation cloud diagram shows that the intensity of deformation gradually increases within the stratum in the process of approaching the fault. The stratum deformation will have an important influence on the structural deformation force as the displacement boundary of structural deformation, and the reference index of stratum deformation intensity is established concerning the limit value of structural deformation in the Code for Seismic Design of Buildings [40]. Combined with the range of displacement difference obtained from the calculation, the deformation difference is graded to form the intensity grading standard, and the ground deformation field is defined into three grades according to the deformation difference: no influence, weak influence, and strong influence. The intensity grades of stratigraphic deformation and the corresponding thresholds of stratigraphic deformation indicators are shown in Table 6.

Stratigraphic Deformation Rate
Intensity 0-0.002 Little Impact 0.002-0.005 Weak Impact >0.005 Strong Impact The maximum value of stratigraphic deformation doubles with the doubling of the fault dislocation. The distribution trends are the same. To further clarify the deformation difference within the soil layer, the overall deformation difference cloud diagram and local enlargement diagram of the soil body under different earthquake levels are made, which are shown in Figure 11. The deformation cloud diagram shows that the intensity of deformation gradually increases within the stratum in the process of approaching the fault. The stratum deformation will have an important influence on the structural deformation force as the displacement boundary of structural deformation, and the reference index of stratum deformation intensity is established concerning the limit value of structural deformation in the Code for Seismic Design of Buildings [40]. Combined with the range of displacement difference obtained from the calculation, the deformation difference is graded to form the intensity grading standard, and the ground deformation field is defined into three grades according to the deformation difference: no influence, weak influence, and strong influence. The intensity grades of stratigraphic deformation and the corresponding thresholds of stratigraphic deformation indicators are shown in Table 6. The maximum value of stratigraphic deformation doubles with the doubling of the fault dislocation. The distribution trends are the same. To further clarify the deformation difference within the soil layer, the overall deformation difference cloud diagram and local enlargement diagram of the soil body under different earthquake levels are made, which are shown in Figure 11.
As can be seen in Figure 11, the overall distribution pattern on the deformation difference distribution pattern is close to that of a funnel-shaped distribution from the fault to the surface. In addition, with the increase in earthquake magnitude, the influence range expands. At the same earthquake magnitude, the intensity of stratigraphic deformation increases with the proximity to the fault. Specifically, at the same horizontal position, the intensity of deformation increases with the increase in stratigraphic depth. At the same depth position, the intensity of deformation increases with the proximity to the epicenter. . Distribution of the degree of stratigraphic deformation. Figure 11. Distribution of the degree of stratigraphic deformation.
In order to further quantify the description of the impact area of different deformation degrees, this paper introduces quantitative index parameters to describe the distribution of the deformation intensity field. The shape of the deformation field is approximated by an ellipse, where the range of intensity field in the vertical direction is defined as b, and the range of intensity field in the horizontal direction is defined as a. The specific location relationship is shown in Figure 12. In this fix, the intensity distribution parameters of the deformation field under different earthquake levels are statistically described in Table 7.
difference distribution pattern is close to that of a funnel-shaped distribution from the fault to the surface. In addition, with the increase in earthquake magnitude, the influence range expands. At the same earthquake magnitude, the intensity of stratigraphic deformation increases with the proximity to the fault. Specifically, at the same horizontal position, the intensity of deformation increases with the increase in stratigraphic depth. At the same depth position, the intensity of deformation increases with the proximity to the epicenter.
In order to further quantify the description of the impact area of different deformation degrees, this paper introduces quantitative index parameters to describe the distribution of the deformation intensity field. The shape of the deformation field is approximated by an ellipse, where the range of intensity field in the vertical direction is defined as b, and the range of intensity field in the horizontal direction is defined as a. The specific location relationship is shown in Figure 12. In this fix, the intensity distribution parameters of the deformation field under different earthquake levels are statistically described in Table 7.

Discussion
Aiming at the objective of qualitative and quantitative research on the degree of deformation and deformation distribution of the coseismic deformation field in the nearfault region, this paper establishes a model of stratigraphic deformation in the Sichuan-Tibet region based on the stratigraphic and fault parameters and uses dislocation model theory to calculate the degree of deformation of the stratigraphic deformation field around the fault after an earthquake. By introducing depth-related parameters, a more refined consideration of the coseismic displacement field is achieved. The calculated deformation law is fitted by multiple regression to form a formula, which extends the analysis field to within the burial depth scale of underground engineering compared with the existing formula. The applicability of statistical formulas in deeply buried strata is improved.

Discussion
Aiming at the objective of qualitative and quantitative research on the degree of deformation and deformation distribution of the coseismic deformation field in the nearfault region, this paper establishes a model of stratigraphic deformation in the Sichuan-Tibet region based on the stratigraphic and fault parameters and uses dislocation model theory to calculate the degree of deformation of the stratigraphic deformation field around the fault after an earthquake. By introducing depth-related parameters, a more refined consideration of the coseismic displacement field is achieved. The calculated deformation law is fitted by multiple regression to form a formula, which extends the analysis field to within the burial depth scale of underground engineering compared with the existing formula. The applicability of statistical formulas in deeply buried strata is improved.
In the depth dimension within the engineering scale, the seismic intensity gradually decays from the location of the source depth to the surface, and the magnitude of the decay is first fast and then slow, showing a nonlinear change characteristic. The concept of "relative depth of burial" is proposed by considering the comprehensive influence of faults and depth of burial on deformation intensity. Under the same earthquake level and different burial depths, the closer the stratigraphic position is to the fault, the higher the intensity; the further away from the fault, the lower the intensity. The influence range of deformation is a funnel-shaped symmetric distribution. The influence range of the stratigraphic deformation field caused by fault dislocation gradually expands with the depth away from the fault, and the displacement maximum is far from the fault location, and the stratigraphic deformation is mainly bending deformation. As it is closer to the fault, the influence range is concentrated, and the ground deformation is mainly shear deformation. Both the range of stratigraphic deformation and the degree of deformation are closely related to the depth.
The model is based on the elasticity theory for the study and has good applicability at the location of deeply buried rock. As the depth of the formation becomes shallower, the mechanical properties of the formation show plasticity. Therefore, the correspondence of the calculated results of formation deformation in the surface area is poorer than that of the deeply buried strata. In addition, since the multiple influencing factors of formation deformation are not independent of each other, there is a coupling effect. In this paper, we ignore the influence of the fault sliding direction, fault dip angle, sliding surface sliding distribution, and other factors on the distribution of the ground deformation field and only consider the fault misalignment under the upright slip fault and the depth of the strata relative to the fault as two main factors and consider the fault as a uniform sliding on a rectangular plane so that the ground deformation profile is symmetrically distributed. The amount of misalignment corresponds to the energy released by the actual seismic magnitude, while in practice there is inhomogeneity in the fault surface sliding. There are some methods to consider the inhomogeneity and analyze the effect of this inhomogeneity, and in further research, more abundant factors are considered to improve the existing model.

Conclusions
In this paper, according to the crustal stratification model of the Sichuan-Tibet region, the distribution pattern of the stratigraphic deformation field in the vicinity of the fault is obtained by combining theoretical analysis and surface survey data, using multiple data joint inversion to obtain fault information and combining with stratigraphic parameter orthorectification calculation. Compared with previous studies, it makes up for the lack of description of the deformation degree within the strata. The effects of various factors are considered, and their parameter influence laws are analyzed separately and targeted for simplification in the subsequent quantitative analysis. For the first time, the concept of intensity is introduced to describe the degree of stratigraphic dislocation deformation, and the research scope is extended from the surface to the scale where the burial depth of underground works is located, and the results of the study can make a preliminary estimation of the degree of stratigraphic deformation after the occurrence of earthquakes. Finally, formulas are fitted to the maximum displacement based on the main influencing factors, and a practical partitioning of the deformation influence range is carried out.
Synthesizing the above work, the main conclusions obtained are as follows.
(1) The stratigraphic deformation field is significantly different in terms of deformation volume and deformation range within a certain depth at the surface and underground. When analyzing the seismic intensity, it is necessary to consider the distribution of the intensity of the ground deformation field at different locations within the stratum under two factors: the burial depth of the fault and the burial depth of the structure. (2) The application range of the existing surface rupture formula is extended to a certain depth below the ground surface, and after considering the relative burial depth factor, the formula is fitted with the ground deformation law to facilitate the preliminary estimation of the maximum displacement of the strata under different earthquake levels and depths. (3) Based on the concept of seismic intensity, the deformation intensity index is introduced to describe the deformation degree and distribution in different areas of the subsurface deformation field, and the practical zoning is carried out according to the ground deformation law, which provides a useful reference for structural deformation prevention and post-earthquake damage prediction at different depths and ranges in the near-fault area under the corresponding seismic magnitude.