Analysis of Seismic Damage Zones: A Case Study of the Ordovician Formation in the Shunbei 5 Fault Zone, Tarim Basin, China

Fault damage zone has an important influence on subsurface fluid flow and petrophysical properties. Therefore, it is of great significance to study the characteristics of fault damage zone for oil and gas development of ultra-deep carbonate formation. This study uses seismic data and the derived variance attribute to identify two types of damage zones and analyze the spatial geometric characteristics of the damage zones. The results show that the type 1 damage zone is wider than the type 2 damage zone. The width of damage zones distributed on both sides of the Shunbei 5 fault core shows obvious asymmetry, and the damage zone width and throw conforms to the typical power-law distribution on the log-log plot. We discuss the factors affecting the width of the damage zone and its formation process. Finally, we discuss the influence of the damage zones on oil and gas exploration. It seems that the seismic variance attribute is a useful technique for characterizing the ultra-deep strike-slip fault damage zones.


Introduction
In recent years, a large number of field observations have shown that the fault zone is mainly composed of a narrow fault core and wide damage zone ( Figure 1) [1,2]. The fault core is the place where the displacement of the fault is concentrated, and the stratigraphy has experienced strong structural deformation. The fault core is usually composed of fault gouge and breccias. The damage zone is immediately adjacent to the fault core and usually develops multiple groups of fractures and secondary faults [3,4]. The damage zone is the product of the formation rheology and displacement on faults. Due to the complex fracture network within damage zones, the geometry size and architecture of damage zones have a significant effect on the migrations, and accumulation of subsurface fluids [5,6]. Therefore, the characterization of the structure and geometry of the damage zone is of great significance to evaluating the integrity of waste-disposal sites and CO 2 sequestration sites, and the exploitation of hydrocarbon in the areas that include the damage zone [1]. In the past, different methods were used to describe the damage zone at different positions. For the damage zone exposed to the surface, the functional relationship between the fracture density and the distance to the fault core was established to characterize the geometric surface, the functional relationship between the fracture density and the distance to the fault core was established to characterize the geometric characteristics of the damage zone [7]. Lin et al. made a field survey on active strike-slip fault damage zones with granite and rhyolitic tuff in southwest Japan and determined the width of the fault damage zone of different scales [8]. Nevertheless, the characterization of the underground damage zone mainly depends on geophysical logging, numerical simulation, seismic data, and analog models [9]. Lyu et al. used conventional logs and imaging logs to characterize the tight sandstone fracture zone in the upper Permian in southwest Ordos Basin, China, and analyzed the logging response characteristics of the tight sandstone damage zone [10]. With the development of seismic technology, a large number of researchers try to identify and describe the underground fault damage zone by using related seismic technology [11][12][13][14]. Botter et al. established discrete element models to simulate fracture zones of faults with different displacements and simulated the seismic response of fracture zones using ray-based prestack depth migration forward modeling technology [15]. Liao et al. used a variety of different seismic attributes (such as variance, coherence, curvature, etc.) and a clay model to describe the deep underground damage zone [16][17][18]. All the cases mentioned above prove the effectiveness of using seismic data and related attributes to characterize the subsurface fault damage zone.
Recently, the Northwest Oilfield Company (Sinopec) has found high-yield industrial oil and gas flow in the ultra-deep strike-slip fault damage zone in the Shunbei area [19][20][21][22]. However, the yield of oil wells deployed in different segments of the same strike-slip fault zone is quite different, and the characteristics of the fault damage zone in this area lack systematic study. In this paper, firstly, according to the original seismic profile, two structural styles of the Shunbei 5 fault zone are identified, and the geometric characteristics of the fracture zone corresponding to the two structural styles are analyzed by using the seismic variance attribute. Then, based on the fact that it is easier to identify the vertical throw of strike-slip fault from seismic data, we further analyze the relationship between the damage zone width and fault throw and propose a law that describes the evolution of the Shunbei 5 fracture zone in the study area. Finally, the possible effects of the two types of damage zones on the production of oil and gas wells are discussed.

Geological Setting
The Shunbei 5 fault zone is located in northern Shuntuo low uplift, which is adjacent to the Tabei Uplift in the north, Tazhong Uplift in the south, Manjiaer Depression in the East, and Awati Depression in the west (Figure 2a,b). The study area has undergone multistage tectonic movements. The main stages of structural evolution are as follows: In the Early Cambrian, due to the rapid tension and rifting around the Tarim plate, the sea level rose rapidly, and the Shuntuo uplift was located in the slope shelf sedimentary facies belt. At this time, the source rocks were widely developed [20]. In the Early Caledonian period, the uplift was high in the northwest and low in the southeast, due to the reverse boundary fault of the uplift. In the middle and late Caledonian, influenced by the closing of the North Kunlun ocean and the subduction of the South Tianshan ocean crust in the south margin of the Tarim Basin, the Tarim Basin is in a tectonic compressive stress environment [23]. In the Hercynian, influenced by the strong orogeny in the south of Tarim Basin, the Shunbei area continued to uplift, and a series of strike-slip faults developed in the uplift. In the Indosinian-Yanshanian period, the Shunbei uplift inherited the structural characteristics of the late Hercynian, and the local faults were further complicated by the influence of the structural movement. The Himalayan tectonic movement is the adjustment and finalization period of low uplift and fault in the study area. Since the Cenozoic, with the rapid uplift of surrounding mountains, the Tarim Basin has continued to subside, but the secondary structural units have remained stable [24].
As shown in Figure 3, the Ordovician strata in the study area are well-developed. According to the drilling data and published literature, the Ordovician strata can be divided into five formations from top to bottom. The upper Ordovician includes Sangtamu Fm., Lianglitage Fm., and Qiarbake Fm., and the lithology is mainly mudstone and marlstone; The middle and lower Ordovician includes Yijianfang Fm., Yingshan Fm., and Penglaiba Fm., and the lithology is mainly carbonate rock. The upper members of Yijianfang Fm. and Yingshan Fm. are limestone developed in open platform facies, and the lower members of Yingshan Fm. and Penglaiba Fm. are dolomite formed in restricted platform facies [25,26]. The middle-lower Ordovician and the upper Ordovician constitute a good reservoir cap assemblage. Oil and gas occur in the carbonate strata of Yijianfang Fm. and Yingshan Formation.

Data and Methods
In this study, high-resolution 3-d seismic data covering an area of about 90 Km 2 were used to characterize the Shunbei 5 fault zone. The time-migrated seismic data analyzed in this paper was collected by Northwest Oilfield Company (Sinopec) in 2014. The bin size of the observation system used to collect the data is 25 * 25 m, and the frequency range is 10~50 Hz. We used this data to carry out detailed horizon tracking and fault interpretation. According to the structural map of surface T7 4 (Figure 4a,b), the Shunbei 5 fault zone comprises several NW trending fault segments. On the seismic section, the typical positive flower structure and the near-vertical segment of the Shunbei 5 fault indicate its strike-slip property (Figure 5a,b). Because oil and gas are mainly distributed in the Ordovician strata, we mainly focus on the characteristics of the Ordovician fault damage zone. The Shunbei 5 fault is a strike-slip fault located in the deep underground; the study on the characteristics of its associated damage zones cannot be the same as that of outcrop fault damage zones. At present, the 3D geometry and structural characteristics of the subsurface fault damage zone are mainly revealed by drilling data and seismic data. In order to clarify the geometric and structural characteristics of the damage zone of the Shunbei 5 fault, the geological horizon on the seismic profile is calibrated by using a seismic well tie, and horizon tracking and fault interpretation are carried out on this basis. Three seismic reflection surfaces are traced, which are the top surface (T7 0 ), the Middle surface (T7 4 ), and the bottom surface (T8 0 ) of the Ordovician strata. Due to the limitation of the resolution of seismic data, it is a challenge to accurately describe the 3D structure characteristics of the deep underground fault damage zone. However, recent studies have shown that it is an effective and feasible method to analyze the characteristics of underground fault damage zone by using seismic amplitude data and its derived seismic attributes [17,18].
Then, we choose the seismic variance reflecting the discontinuities in the horizontal continuity of amplitude to characterize the geometric and structural characteristics of the damage zone. The calculation method of seismic variance volume is similar to that used

Data and Methods
In this study, high-resolution 3-d seismic data covering an area of about 90 Km 2 were used to characterize the Shunbei 5 fault zone. The time-migrated seismic data analyzed in this paper was collected by Northwest Oilfield Company (Sinopec) in 2014. The bin size of the observation system used to collect the data is 25 * 25 m, and the frequency range is 10~50 Hz. We used this data to carry out detailed horizon tracking and fault interpretation. According to the structural map of surface T 7 4 ( Figure 4a,b), the Shunbei 5 fault zone comprises several NW trending fault segments. On the seismic section, the typical positive flower structure and the near-vertical segment of the Shunbei 5 fault indicate its strike-slip property (Figure 5a,b). Because oil and gas are mainly distributed in the Ordovician strata, we mainly focus on the characteristics of the Ordovician fault damage zone. The Shunbei 5 fault is a strike-slip fault located in the deep underground; the study on the characteristics of its associated damage zones cannot be the same as that of outcrop fault damage zones. At present, the 3D geometry and structural characteristics of the subsurface fault damage zone are mainly revealed by drilling data and seismic data. In order to clarify the geometric and structural characteristics of the damage zone of the Shunbei 5 fault, the geological horizon on the seismic profile is calibrated by using a seismic well tie, and horizon tracking and fault interpretation are carried out on this basis. Three seismic reflection surfaces are traced, which are the top surface (T 7 0 ), the Middle surface (T 7 4 ), and the bottom surface (T 8 0 ) of the Ordovician strata. Due to the limitation of the resolution of seismic data, it is a challenge to accurately describe the 3D structure characteristics of the deep underground fault damage zone. However, recent studies have shown that it is an effective and feasible method to analyze the characteristics of underground fault damage zone by using seismic amplitude data and its derived seismic attributes [17,18].
Then, we choose the seismic variance reflecting the discontinuities in the horizontal continuity of amplitude to characterize the geometric and structural characteristics of the damage zone. The calculation method of seismic variance volume is similar to that used by Liao et al. [18]. In the calculation of variance volume, the length of the filter is 3 * 3 in the inline and crossline, and 15 samples in the vertical direction are used as smooth filtering parameters. Based on the characteristics of the seismic section and map variance property, we have identified two typical structural patterns in the study area, namely, type 1 and type 2 ( Figure 6). Previous studies have shown that the width of the damage zone is different, due to the different formation processes or overlapping mode of fault [1,27]. To analyze the characteristics of the damage zone corresponding to the two typical structural styles in this area, the variance extracted along the surface (T 7 4 ) is selected. In the direction perpendicular to the strike of the Shunbei 5 fault, a total of 17 sampling lines are arranged with an interval of 500 m ( Figure 4a). Then, the geometric characteristics of the two typical fault zones are analyzed by using the variance values extracted from these sampling lines ( Figure 7). Besides, to analyze the various characteristics of the damage zone width in the depth direction, L4 and L8 profiles are selected as the representatives of the two types of damage zones, and a set of variance values are obtained every 100 ms from 4300 ms to 5300 ms. Then, we interpreted these fracture zones at different depths according to the variance value and plotted the damage zone width versus depth ( Figure 8). Finally, because it is difficult to determine the horizontal displacement of deep strike-slip fault, to further explore the formation mechanism and evolution process of the Shunbei 5 fault, we made a statistical analysis on the throw and fracture zone width of the Shunbei 5 fault ( Figure 9). The throw of Shunbei 5 main fault is measured on seismic profiles, and the scaling relationship between the vertical throw and the width of damage zones is analyzed ( Figure 10). by Liao et al. [18]. In the calculation of variance volume, the length of the filter is 3 * 3 in the inline and crossline, and 15 samples in the vertical direction are used as smooth filtering parameters. Based on the characteristics of the seismic section and map variance property, we have identified two typical structural patterns in the study area, namely, type 1 and type 2 ( Figure 6). Previous studies have shown that the width of the damage zone is different, due to the different formation processes or overlapping mode of fault [1,27]. To analyze the characteristics of the damage zone corresponding to the two typical structural styles in this area, the variance extracted along the surface(T7 4 ) is selected. In the direction perpendicular to the strike of the Shunbei 5 fault, a total of 17 sampling lines are arranged with an interval of 500 m (Figure 4a). Then, the geometric characteristics of the two typical fault zones are analyzed by using the variance values extracted from these sampling lines (Figure 7). Besides, to analyze the various characteristics of the damage zone width in the depth direction, L4 and L8 profiles are selected as the representatives of the two types of damage zones, and a set of variance values are obtained every 100 ms from 4300 ms to 5300 ms. Then, we interpreted these fracture zones at different depths according to the variance value and plotted the damage zone width versus depth ( Figure  8). Finally, because it is difficult to determine the horizontal displacement of deep strikeslip fault, to further explore the formation mechanism and evolution process of the Shunbei 5 fault, we made a statistical analysis on the throw and fracture zone width of the Shunbei 5 fault (Figure 9). The throw of Shunbei 5 main fault is measured on seismic profiles, and the scaling relationship between the vertical throw and the width of damage zones is analyzed ( Figure 10).

Characteristics of Shunbei 5 Fault Damage Zone
We interpret the yellow dotted rectangle as the Shunbei 5 fault zone (Figure 6a,b), and Figure 6c,d show the seismic variance profiles corresponding to two typical damage zones, respectively. In the seismic amplitude profile (profile L8, Figure 6a), the type 1 damage zone is characterized by a disordered and discontinuous in-phase axis, while the variance section is characterized by a wide discontinuity area with a high variance value. Compared with type 1, the seismic amplitudes of type 2 (profile L4, Figure 6b) damage zone are more continuous, and the amplitudes only break at the fault core, showing a narrow area with a high variance value on the variance section. To analyze the characteristics of the two damage zones, we extract the variance values of the two fault

Characteristics of Shunbei 5 Fault Damage Zone
We interpret the yellow dotted rectangle as the Shunbei 5 fault zone (Figure 6a,b), and Figure 6c,d show the seismic variance profiles corresponding to two typical damage zones, respectively. In the seismic amplitude profile (profile L8, Figure 6a), the type 1 damage zone is characterized by a disordered and discontinuous in-phase axis, while the variance section is characterized by a wide discontinuity area with a high variance value. Compared with type 1, the seismic amplitudes of type 2 (profile L4, Figure 6b) damage zone are more continuous, and the amplitudes only break at the fault core, showing a narrow area with a high variance value on the variance section. To analyze the characteristics of the two damage zones, we extract the variance values of the two fault

Characteristics of Shunbei 5 Fault Damage Zone
We interpret the yellow dotted rectangle as the Shunbei 5 fault zone (Figure 6a,b), and Figure 6c,d show the seismic variance profiles corresponding to two typical damage zones, respectively. In the seismic amplitude profile (profile L8, Figure 6a), the type 1 damage zone is characterized by a disordered and discontinuous in-phase axis, while the variance section is characterized by a wide discontinuity area with a high variance value. Compared with type 1, the seismic amplitudes of type 2 (profile L4, Figure 6b) damage zone are more continuous, and the amplitudes only break at the fault core, showing a narrow area with a high variance value on the variance section. To analyze the characteristics of the two damage zones, we extract the variance values of the two fault segments on different sections. Because high variance coefficients indicate the discontinuities of subsurface geological objects [18]. Zones with higher variance values than a background value are interpreted as damage zones. In Figure 7, we interpret the grey area with the seismic variance greater than 0.1 as the damage zone, the yellow area with the variance greater than 0.6 as the fault core, and the blank area with variance lower than 0.1 as protolith. Based on the structural style of seismic sections and formation rheology, we separate the seismic damage zone into two types. The total width of type 1 (Figure 7a) [28]. To further investigate the relationship of the damage zones with depth, a series of variance profiles were drawn at intervals of 100 ms (Figure 8a,b). Type 1 has a broader damage zone width between 4500 ms and 4800 ms (Figure 8c), and the maximum width is 1400 m, while it is narrow above 4500 ms and below 4800 ms. Compared with type 1, type 2 was only slightly broader, around 4800 ms (Figure 8d). The maximum width is 970 m, the minimum width is 125 m, and the average width is 342 m.  [18]. Zones with higher variance valu than a background value are interpreted as damage zones. In Figure 7, we interpret t grey area with the seismic variance greater than 0.1 as the damage zone, the yellow ar with the variance greater than 0.6 as the fault core, and the blank area with variance low than 0.1 as protolith. Based on the structural style of seismic sections and formati rheology, we separate the seismic damage zone into two types. The total width of type (Figure 7a)

Relationship of Damage Zone Width with Throw
It is well known that there are various functional relationships among the various fault attributes, such as log-normal, exponential laws, power-law distribution, and so on [29][30][31]. The relationship between the fault displacement and the width of the damage zone is following the power-law distribution on the log-log diagram, which has been recognized by most researchers and justified in the frame of complexity theory for earth other planets, too [32][33][34][35]. However, the power-law distribution is still uncertain, due to the influence of sampling error or the complex fault evolution process. In this paper, because it is difficult to determine the horizontal displacement of deep strike-slip fault, we choose the vertical throw as the proxy of the displacement to analyze the relationship between the fault throw and the width of the damage zone. Based on 17 seismic profiles perpendicular to the strike of the Shunbei 5 fault and their seismic variance attributes, we measured the width of the damage zone and the throw at the surface (T7 4 ). To reduce the sampling error, only the damage zones and throw of the main fault are measured and

Relationship of Damage Zone Width with Throw
It is well known that there are various functional relationships among the various fault attributes, such as log-normal, exponential laws, power-law distribution, and so on [29][30][31]. The relationship between the fault displacement and the width of the damage zone is following the power-law distribution on the log-log diagram, which has been recognized by most researchers and justified in the frame of complexity theory for earth other planets, too [32][33][34][35]. However, the power-law distribution is still uncertain, due to the influence of sampling error or the complex fault evolution process. In this paper, because it is difficult to determine the horizontal displacement of deep strike-slip fault, we choose the vertical throw as the proxy of the displacement to analyze the relationship between the fault throw and the width of the damage zone. Based on 17 seismic profiles perpendicular to the strike of the Shunbei 5 fault and their seismic variance attributes, we measured the width of the damage zone and the throw at the surface (T7 4 ). To reduce the sampling error, only the damage zones and throw of the main fault are measured and

Relationship of Damage Zone Width with Throw
It is well known that there are various functional relationships among the various fault attributes, such as log-normal, exponential laws, power-law distribution, and so on [29][30][31]. The relationship between the fault displacement and the width of the damage zone is following the power-law distribution on the log-log diagram, which has been recognized by most researchers and justified in the frame of complexity theory for earth other planets, too [32][33][34][35]. However, the power-law distribution is still uncertain, due to the influence of sampling error or the complex fault evolution process. In this paper, because it is difficult to determine the horizontal displacement of deep strike-slip fault, we choose the vertical throw as the proxy of the displacement to analyze the relationship between the fault throw and the width of the damage zone. Based on 17 seismic profiles perpendicular to the strike of the Shunbei 5 fault and their seismic variance attributes, we measured the width of the damage zone and the throw at the surface (T 7 4 ). To reduce the sampling error, only the damage zones and throw of the main fault are measured and analyzed, while the secondary fault was excluded. At the same time, based on the average velocity of limestone strata in the study area, the time domain throw is converted into depth domain throw. It is found that the damage zone width and vertical throw of the Shunbei 5 fault vary along the strike of the fault. From Figure 9a, the width of the type 1 damage zone is positively correlated with throw, while that of the type 2 damage zone is not. The average width of the type 1 damage zone is about 831 m, and the average throw is about 270 m, while the average width of the type 2 damage zone is about 582 m, and the average throw is about 89 m (Figure 9b). In general, the type 1 damage zone is wider, and the throw is larger. Then, we compile the data of [36] in the Tazhong Uplift with our data (Figure 10). Although the Study area and the Tazhong Uplift are located in different structural units of Tarim Basin. It is found that the damage zone width and displacement (throw) in this area also obey the power-law distribution on the log-log plot, and the ratio of width to displacement is between 2-20. analyzed, while the secondary fault was excluded. At the same time, based on the average velocity of limestone strata in the study area, the time domain throw is converted into depth domain throw. It is found that the damage zone width and vertical throw of the Shunbei 5 fault vary along the strike of the fault. From Figure 9a, the width of the type 1 damage zone is positively correlated with throw, while that of the type 2 damage zone is not. The average width of the type 1 damage zone is about 831 m, and the average throw is about 270 m, while the average width of the type 2 damage zone is about 582 m, and the average throw is about 89 m (Figure 9b). In general, the type 1 damage zone is wider, and the throw is larger. Then, we compile the data of [36] in the Tazhong Uplift with our data (Figure 10). Although the Study area and the Tazhong Uplift are located in different structural units of Tarim Basin. It is found that the damage zone width and displacement (throw) in this area also obey the power-law distribution on the log-log plot, and the ratio of width to displacement is between 2-20.

Factors Affecting the Damage Zone Width of the Shunbei 5 Fault Zone
Previous studies have shown that the width of the fault damage zone is closely related to the buried depth of the fault, the stress field during the formation of the fault, the type of protolith, the interaction between fault segments, and the later diagenesis [37,38]. The Ordovician carbonate rocks in this study area are mainly composed of micritic  (Figure 9b). In general, the type 1 damage zone is wider, and the throw is larger. Then, we compile the data of [36] in the Tazhong Uplift with our data (Figure 10). Although the Study area and the Tazhong Uplift are located in different structural units of Tarim Basin. It is found that the damage zone width and displacement (throw) in this area also obey the power-law distribution on the log-log plot, and the ratio of width to displacement is between 2-20.

Factors Affecting the Damage Zone Width of the Shunbei 5 Fault Zone
Previous studies have shown that the width of the fault damage zone is closely related to the buried depth of the fault, the stress field during the formation of the fault, the type of protolith, the interaction between fault segments, and the later diagenesis [37,38]. The Ordovician carbonate rocks in this study area are mainly composed of micritic

Factors Affecting the Damage Zone Width of the Shunbei 5 Fault Zone
Previous studies have shown that the width of the fault damage zone is closely related to the buried depth of the fault, the stress field during the formation of the fault, the type of protolith, the interaction between fault segments, and the later diagenesis [37,38]. The Ordovician carbonate rocks in this study area are mainly composed of micritic limestone developed in open platform facies, and their buried depth ranges from 7300 to 7800 m. On the top of the Yijianfang Formation (surface T 7 4 ), the width of the damage zone along the strike of the Shunbei 5 fault shows segmentation. The width of the damage zone in the compressive uplifting segments (type 1) is larger than that in the strike-slip segments (type 2). The Shunbei 5 fault is a dextral strike-slip fault, which is subjected to strong compression and torsion stress in the compressive uplifting segments [39]. In the strike-slip segments, the mechanical interaction is mainly distributed at the end of the segment, and the overall compressive and torsional stress is relatively weak. On the profile, the width of the damage zone varies with the depth. According to the results of seismic variance analysis, the width of the Shunbei 5 fault damage zone is narrow at the root and upper part of the fault, and wide in the middle. This feature shows that the maximum displacement is mainly in the middle region of the Shunbei 5 fault [40]. Besides, the geometric characteristics of seismic profiles corresponding to the two damage zones show that the stress fields of different segments of the Shunbei 5 fault are different. Furthermore, the lithology between T 7 0 and T 7 4 is mainly mudstone, argillaceous limestone, and marl, while the lithology between T 7 4 and T 8 0 is mainly micritic limestone (Figure 3). The competency of micritic limestone is higher than that of mudstone and marl. During the formation and development of strike-slip fault, rocks with higher competency are more prone to stress concentration, which leads to rock deformation and forms a wider damage zone [41].

Implications for Growth of the Shunbei 5 Fault Damage Zone
A large number of studies have shown that there are various types of damage zones associated with strike-slip faults and their formation mechanisms. For the evolution and formation of damage zones, many researchers have proposed different models, such as the process zone model [42], fault wear model [43], and coseismic damage model. According to the relationship between the damage zones and their location along with the profile of a strike-slip fault, Kim et al. separated the damage zone into three types: Tip-damage zone, linking-damage zone, and wall-damage zone. Type 1 and type 2 in this paper correspond to the linking -and wall damage zones. Before forming a large strike-slip fault, it comprises a series of nearly parallel small faults. After multiple tectonic movements, these small faults gradually connected and overlapped, forming a large fault. The deformation characteristics of the seismic amplitude profile show that type 1 is a typical contraction zone. During forming strike-slip fault, the tips of two overstepping segments are pinned in the contractional zone, and strong mechanical interaction occurs, which results in the strong uplift of the strata in the contractional zone and wider damage zone. The type 2 damage zone mainly comes from the merger and coalescence between different fault segments. In addition, the heterogeneous stress field produced during the activity of strike-slip fault is geometrically shown as the ribbon effect of strike-slip fault, which may be the fundamental reason for the heterogeneous deformation of rocks around the fault and forming damage zones with different width and high seismic variance values along the strike of the Shunbei 5 fault [44]. Also, the study area has experienced multiple tectonic movements, and the multiple reactivations of the Shunbei 5 fault are the direct cause of forming the complex damage zones. Moreover, the roughness of the fault plane, the shear strength of the fault, and the spatial configuration of multiple segments also have a great influence on the dimension of the fault fracture zone. However, limited by the data, we will not discuss these issues.

Influence of Fault Damage Zone Characteristics to Hydrocarbon Exploration
Previous studies have shown that the damage zones associated with a strike-slip fault may be the main factor of reservoir quality and hydrocarbon distribution [22]. Therefore, the drilling targets in the Shunbei area are mainly fractured zones associated with large strike-slip faults. However, the drilled productivity of different types of damage zones associated with strike-slip faults shows a great difference, in general, the average unit pressure drop production of a single well deployed in the type 1 damage zone is 1189.7 t·MPa −1 , and that deployed in the type 2 damage zone is 524.5 t·MPa −1 [39,45]. It also directly indicates the difference in hydrocarbon distribution in different segments along a strike-slip fault. According to [45], pressure interference tests were carried out on wells deployed in the study area. In the process of testing, Shunbei 5 well is used as a stimulation well, Shunbei 5-B well, and Shunbei 5-C well are used as monitoring wells. The test results show that the unit propagation velocity of the type 1 damage zone is 9.8 m/h, and the type 2 damage zone is an independent reservoir unit controlled by a single well. Based on the above published dynamic data and the characteristic data of the damage zones, we can infer that the larger width of the type 1 damage zone leads to more reservoir space and better connectivity of fracture network, so it has greater production potential than the type 2 damage zone.

Conclusions
In this study, our results verify that it is an effective method to use seismic variance to characterize the deep fault damage zone. Based on seismic data and seismic variance attributes, the Ordovician strike-slip fault damage zones in the Shunbei area are analyzed, and the following conclusions are obtained. The results show that there is a positive correlation between the width of the type 1 damage zone and the vertical throw, while the correlation between the width of the type 2 damage zone and the vertical throw is not strong. Moreover, the ratio of the damage zone width to the vertical throw is mainly distributed in the range of 2-20, and the power-law distribution is shown in the log-log diagram. The damage zones developed in the middle and upper Ordovician strata are wider than those in the upper Ordovician strata, especially in the vicinity of the lithologic transition surface. Based on the published production data and the research results of the Shunbei 5 damage zone in this paper, we infer that the type 1 damage zone has greater oil and gas exploration potential than the type 2 damage zone.