A quantitative evaluation method for the intensity of strike-slip fault activity

Abstract

The present invention belongs to the technical field of oil and gas reservoir exploration and development, and discloses a method for quantitatively evaluating the activity intensity of a strike-slip fault. A fine structural interpretation is carried out based on high-precision three-dimensional seismic data. On this basis, the plane distribution range of the strike-slip fault crushing zone is obtained by fine calibration of the drilled wells, and the lateral deformation width and vertical deformation height difference data that can characterize the deformation degree of the strike-slip fault are further obtained, thereby obtaining an activity intensity coefficient that can characterize the activity intensity of the strike-slip fault. Finally, the activity intensity coefficient is used to quantitatively characterize the activity intensity of the strike-slip fault, which provides an operable quantitative evaluation index for the quantitative evaluation of the activity intensity of the strike-slip fault zone, and provides a strong technical support for the efficient exploration and development of fault-controlled fracture-cavity carbonate rocks.

Description

Technical Field

[0001] The invention belongs to the technical field of oil and gas reservoir exploration and development, and particularly relates to a method for quantitatively evaluating the intensity of strike-slip fault activity. Background Art

[0002] Exploration, development and production practice have confirmed that the intensity of strike-slip fault activity has an important influence on the distribution of Ordovician fracture-vuggy carbonate reservoirs, oil and gas reservoir properties and single well production. However, the Ordovician inner strike-slip faults are almost vertical. Existing technologies such as the growth index method, activity rate method, drop analysis method, and slip analysis method are difficult to achieve quantitative evaluation of the activity intensity of small slip-distance strike-slip faults within the inner plate of this type of craton due to their theoretical limitations. Summary of the invention

[0003] The purpose of the present invention is to provide a method for quantitatively evaluating the intensity of strike-slip fault activity, so as to achieve quantitative evaluation of the intensity of strike-slip fault activity with a small slip distance and nearly vertical inside the Ordovician system in a craton plate.

[0004] The present invention is achieved through the following technical solutions:

[0005] A method for quantitatively evaluating the intensity of strike-slip fault activity comprises the following steps:

[0006] Step 1: Based on the high-precision 3D seismic data in the depth domain, the microstructure map of the top surface of the Ordovician carbonate rock and the contour boundary of the strike-slip fault zone on the top surface of the Ordovician carbonate rock are obtained;

[0007] Step 2, using the contour boundary of the strike-slip fault zone on the top surface of the Ordovician carbonate rock obtained in step 1, to obtain a lateral deformation width measurement line of the strike-slip fault zone on the top surface of the Ordovician carbonate rock;

[0008] Step 3, using the lateral deformation width measurement line of the strike-slip fault zone on the top surface of the Ordovician carbonate rock obtained in step 2 to obtain the lateral deformation width data of the top surface of the Ordovician carbonate rock along the strike-slip fault direction;

[0009] Step 4, using the lateral deformation width data of the top surface of the Ordovician carbonate rock along the strike-slip fault obtained in step 3, to obtain the lateral deformation width coefficient of the top surface of the Ordovician carbonate rock along the strike-slip fault;

[0010] Step 5, using the microstructure map of the top surface of the Ordovician carbonate rock and the contour boundary of the top surface fracture zone of the Ordovician carbonate rock obtained in step 1, obtain the axial line of the strike-slip fault fracture zone on the top surface of the Ordovician carbonate rock;

[0011] Step 6, using the microstructure map of the top surface of the Ordovician carbonate rock obtained in step 1, the lateral deformation width measurement line of the strike-slip fault and crushing zone of the top surface of the Ordovician carbonate rock obtained in step 2, and the axial line of the strike-slip fault and crushing zone of the top surface of the Ordovician carbonate rock obtained in step 5, obtain the vertical deformation height difference data of the top surface of the Ordovician carbonate rock along the strike-slip fault direction;

[0012] Step 7, using the vertical deformation height difference data of the top surface of the Ordovician carbonate rock along the strike-slip fault obtained in step 6, to obtain the vertical deformation height difference coefficient of the top surface of the Ordovician carbonate rock along the strike-slip fault;

[0013] Step 8, using the lateral deformation width coefficient of the top surface of the Ordovician carbonate rock along the strike-slip fault obtained in step 4 and the vertical deformation height difference coefficient of the top surface of the Ordovician carbonate rock along the strike-slip fault obtained in step 7, obtain the activity intensity coefficient of the top surface of the Ordovician carbonate rock along the strike-slip fault;

[0014] The activity intensity of strike-slip faults is determined based on the activity intensity coefficient of the top surface of Ordovician carbonate rocks along the strike-slip faults.

[0015] Further, obtaining the microstructure map of the top surface of the Ordovician carbonate rock based on the high-precision three-dimensional seismic data in the depth domain includes the following steps:

[0016] 1.1. Based on the high-precision 3D seismic data in the depth domain, the top surface strata of the Ordovician carbonate rocks were carefully interpreted to obtain the top surface strata of the Ordovician carbonate rocks;

[0017] 1.2. Process the top surface of the Ordovician carbonate rock to obtain the microstructural layer of the top surface of the Ordovician carbonate rock;

[0018] 1.3. Using the microstructural horizons of the top surface of the Ordovician carbonate rocks, a microstructural map of the top surface of the Ordovician carbonate rocks is obtained;

[0019] The following steps are involved in obtaining the contour boundary of the strike-slip fault zone on the top of the Ordovician carbonate rock based on high-precision 3D seismic data in the depth domain:

[0020] 1.11. Extract coherent volume attributes based on high-precision 3D seismic data in depth domain;

[0021] 1.12, according to the coherent volume attributes extracted in step 1.11 and the top surface horizon of the Ordovician carbonate rock obtained in step 1.1, obtain the coherent plane attributes of the top surface of the Ordovician carbonate rock;

[0022] 1.13. Using the coherent plane attributes of the top surface of the Ordovician carbonate rock, the coherent plane attribute map of the top surface of the Ordovician carbonate rock is obtained;

[0023] 1.14. Use the coherent plane attribute map of the top surface of the Ordovician carbonate rock to obtain the contour boundary of the strike-slip fault zone on the top surface of the Ordovician carbonate rock.

[0024] Further, step 1.2 is specifically as follows: smoothing the top surface horizon of the Ordovician carbonate rock to obtain a smooth background horizon;

[0025] The microstructural horizons of the top surface of the Ordovician carbonate rock are obtained by subtracting the smooth background horizons of the top surface of the Ordovician carbonate rock from the top surface of the Ordovician carbonate rock.

[0026] Furthermore, in step 1.11, when extracting coherent plane attributes, the longitudinal window parameter setting is centered on the top surface of the Ordovician carbonate rock layer and includes a range of 20m above and below as the longitudinal extraction time window for coherent plane attributes.

[0027] Further, step 1.14 is specifically as follows: determining the coherent attribute threshold of the boundary of the fracture zone contour according to the coherent attribute value at the drilled leakage point on the fracture zone;

[0028] According to the coherent attribute threshold of the fracture zone contour boundary, in the coherent plane attribute map of the top surface of the Ordovician carbonate rock, the coherent attribute threshold contour line is drawn as the contour boundary of the strike-slip fracture zone on the top surface of the Ordovician carbonate rock.

[0029] Further, step 2 is specifically as follows: according to the contour boundary of the strike-slip fault zone on the top surface of the Ordovician carbonate rock obtained in step 1, a line connecting the midpoint of the contour boundary at the beginning of the strike-slip fault and the midpoint of the contour boundary at the end is used as the strike line of the strike-slip fault zone;

[0030] From the starting point of the strike-slip fault, along the strike line of the strike-slip fault zone, a straight line perpendicular to the strike line of the strike-slip fault zone and connecting the contour boundary of the strike-slip fault zone is drawn at intervals of 30m as the lateral deformation width measurement line of the strike-slip fault zone on the top surface of the Ordovician carbonate rock;

[0031] If there are multiple contour boundaries, the straight line connecting the outermost contour boundaries will be used as the measurement line of the lateral deformation width of the strike-slip fault zone on the top surface of the Ordovician carbonate rock.

[0032] Further, step 3 is specifically as follows: starting from the starting point of the strike-slip fault, along the strike-slip fault direction, the length of the lateral deformation width measurement line of the strike-slip fault zone on the top surface of the Ordovician carbonate rock obtained in step 2 is sequentially counted, so as to obtain the lateral deformation width data of the top surface of the Ordovician carbonate rock along the strike-slip fault direction.

[0033] Further, step 4 is specifically as follows: dividing the lateral deformation width data of the top surface of the Ordovician carbonate rock along the strike-slip fault obtained in step 3 by the maximum lateral deformation width data to obtain the lateral deformation width coefficient of the top surface of the Ordovician carbonate rock along the strike-slip fault;

[0034] Step 7 is specifically as follows: the absolute value of the vertical deformation height difference data of the top surface of the Ordovician carbonate rock along the strike-slip fault obtained in step 6 is divided by the maximum absolute vertical deformation height difference data, so as to obtain the vertical deformation height difference coefficient of the top surface of the Ordovician carbonate rock along the strike-slip fault.

[0035] Further, step 5 is specifically as follows: in the microstructure map of the top surface of the Ordovician carbonate rock obtained in step 1, and within the contour boundary of the strike-slip fault and crushing zone on the top surface of the Ordovician carbonate rock obtained in step 1, a maximum structural axial line is drawn along the strike-slip fault direction from the starting point of the strike-slip fault to obtain the axial line of the strike-slip fault and crushing zone on the top surface of the Ordovician carbonate rock.

[0036] Further, step 6 is specifically as follows: in the microstructure map of the top surface of the Ordovician carbonate rock, from the starting point of the strike-slip fault, along the strike-slip fault direction, sequentially read the structural values ​​at the intersection of the lateral deformation width measurement line of the strike-slip fault zone on the top surface of the Ordovician carbonate rock obtained in step 2 and the axial line of the fault and crushing zone on the top surface of the Ordovician carbonate rock obtained in step 5, and use the structural values ​​as the vertical deformation height difference data of the top surface of the Ordovician carbonate rock along the strike-slip fault direction.

[0037] Compared with the prior art, the present invention has the following beneficial technical effects:

[0038] The present invention discloses a method for quantitatively evaluating the intensity of strike-slip fault activity. Based on high-precision three-dimensional seismic data, a fine structural interpretation is carried out. On this basis, the plane distribution range of the strike-slip fault fracture zone is obtained by fine calibration of the drilled wells, and the lateral deformation width and vertical deformation height difference data that can characterize the degree of deformation of the strike-slip fault are further obtained. Then, the activity intensity coefficient that can characterize the intensity of strike-slip fault activity is obtained. Finally, the activity intensity coefficient is used to quantitatively characterize the intensity of strike-slip fault activity, which provides an operational quantitative evaluation index for the quantitative evaluation of the intensity of strike-slip fault zone activity, and provides strong technical support for the efficient exploration and development of fault-controlled fracture-cavity carbonate rocks. Good application effects have been achieved in the study of Ordovician fault-controlled fracture-cavity closures in the Tazhong No. 1 gas field in the Tarim Basin. The technology is used to quantitatively evaluate the intensity of strike-slip fault activity, providing a quantitative index for the division of Ordovician fault-controlled fracture-cavity closures. The design method and ideas of the present invention have good application effects and promotion value for the same industry. BRIEF DESCRIPTION OF THE DRAWINGS

[0039] Figure 1 It is a flow chart of a method for quantitatively evaluating the intensity of strike-slip fault activity according to the present invention;

[0040] Figure 2 To obtain the contour boundary of the strike-slip fault zone on the top surface of the Ordovician carbonate rock;

[0041] Figure 3The measurement line of the lateral deformation width of the strike-slip fault zone on the top surface of the Ordovician carbonate rock is obtained;

[0042] Figure 4 The axial line of the strike-slip fault zone on the top surface of the Ordovician carbonate rock is obtained;

[0043] Figure 5 To obtain the lateral deformation width and vertical deformation data of the top surface of the Ordovician carbonate rock along the strike-slip fault;

[0044] Figure 6 To obtain the lateral deformation width coefficient, vertical deformation height difference coefficient and activity intensity coefficient of the top surface of the Ordovician carbonate rock along the strike-slip fault. DETAILED DESCRIPTION

[0045] The present invention is further described in detail below in conjunction with specific embodiments, which are intended to explain the present invention rather than to limit it.

[0046] like Figure 1 As shown, the present invention discloses a method for quantitatively evaluating the intensity of strike-slip fault activity, comprising the following steps:

[0047] Step 1: Based on the high-precision 3D seismic data in the depth domain, the microstructure map of the top surface of the Ordovician carbonate rock and the contour boundary of the strike-slip fault zone on the top surface of the Ordovician carbonate rock are obtained;

[0048] The following steps are involved in obtaining the microstructure map of the top surface of the Ordovician carbonate rock based on high-precision 3D seismic data in the depth domain:

[0049] 1.1. Based on the high-precision 3D seismic data in the depth domain, the top surface strata of the Ordovician carbonate rocks were carefully interpreted to obtain the top surface strata of the Ordovician carbonate rocks;

[0050] 1.2. Process the top surface of the Ordovician carbonate rock to obtain the microstructural layer of the top surface of the Ordovician carbonate rock;

[0051] 1.3. Using the microstructural horizons of the top surface of the Ordovician carbonate rocks, a microstructural map of the top surface of the Ordovician carbonate rocks is obtained;

[0052] The following steps are involved in obtaining the contour boundary of the strike-slip fault zone on the top of the Ordovician carbonate rock based on high-precision 3D seismic data in the depth domain:

[0053] 1.11. Extract coherent volume attributes based on high-precision 3D seismic data in depth domain;

[0054] 1.12, according to the coherent volume attributes extracted in step 1.11 and the top surface horizon of the Ordovician carbonate rock obtained in step 1.1, obtain the coherent plane attributes of the top surface of the Ordovician carbonate rock;

[0055] 1.13. Using the coherent plane attributes of the top surface of the Ordovician carbonate rock, the coherent plane attribute map of the top surface of the Ordovician carbonate rock is obtained;

[0056] 1.14. Use the coherent plane attribute map of the top surface of the Ordovician carbonate rock to obtain the contour boundary of the strike-slip fault zone on the top surface of the Ordovician carbonate rock, such as Figure 2 As shown;

[0057] Step 2: Figure 3 As shown, the contour boundary of the strike-slip fault zone on the top surface of the Ordovician carbonate rock obtained in step 1 is used to obtain the lateral deformation width measurement line of the strike-slip fault zone on the top surface of the Ordovician carbonate rock;

[0058] Step 3: Figure 5 As shown, the lateral deformation width measurement line of the strike-slip fault zone on the top surface of the Ordovician carbonate rock obtained in step 2 is used to obtain the lateral deformation width data of the top surface of the Ordovician carbonate rock along the strike-slip fault direction;

[0059] Step 4, using the lateral deformation width data of the top surface of the Ordovician carbonate rock along the strike-slip fault obtained in step 3, to obtain the lateral deformation width coefficient of the top surface of the Ordovician carbonate rock along the strike-slip fault;

[0060] Step 5: Use the microstructure map of the top surface of the Ordovician carbonate rock and the contour boundary of the top surface of the Ordovician carbonate rock obtained in step 1, as shown in Figure 4 As shown, the axial line of the strike-slip fault zone on the top surface of the Ordovician carbonate rock is obtained;

[0061] Step 6: Using the microstructure map of the top surface of the Ordovician carbonate rock obtained in step 1, the lateral deformation width measurement line of the strike-slip fault zone of the top surface of the Ordovician carbonate rock obtained in step 2, and the axial line of the strike-slip fault zone of the top surface of the Ordovician carbonate rock obtained in step 5, as shown in FIG. Figure 5 As shown, the vertical deformation height difference data of the top surface of the Ordovician carbonate rock along the strike-slip fault is obtained;

[0062] Step 7, using the vertical deformation height difference data of the top surface of the Ordovician carbonate rock along the strike-slip fault obtained in step 6, to obtain the vertical deformation height difference coefficient of the top surface of the Ordovician carbonate rock along the strike-slip fault;

[0063] Step 8: Using the lateral deformation width coefficient of the top surface of the Ordovician carbonate rock along the strike-slip fault obtained in step 4 and the vertical deformation height difference coefficient of the top surface of the Ordovician carbonate rock along the strike-slip fault obtained in step 7, as follows: Figure 6 As shown, the activity intensity coefficient of the top surface of the Ordovician carbonate rock along the strike-slip fault is obtained;

[0064] The activity intensity of the strike-slip fault is determined based on the strike-slip fault activity intensity coefficient. The larger the value of the strike-slip fault activity intensity coefficient is, the greater the strike-slip fault activity intensity is.

[0065] Among them, in step 1.1, the accuracy of the detailed interpretation of the horizon is 1 line along the plane main survey line and the connecting line.

[0066] In step 1.11, the parameter settings in the coherent attribute window have a significant impact on the effect of fault identification. Through multiple tests and comparative analysis, the optimal parameters are to take 9 sample points along the main survey line and the connecting line in the plane, and 21 depth sample points in the vertical direction.

[0067] When extracting coherent plane attributes, the longitudinal window parameter setting is centered on the top surface of the Ordovician carbonate rock layer and includes a range of 20m above and below as the longitudinal extraction time window for coherent plane attributes.

[0068] Step 1.2 is specifically as follows: smoothing the top surface layer of the Ordovician carbonate rock to obtain a smooth background layer;

[0069] The microstructural horizons of the top surface of the Ordovician carbonate rock are obtained by subtracting the smooth background horizons of the top surface of the Ordovician carbonate rock from the top surface of the Ordovician carbonate rock.

[0070] The horizon smoothing parameter is 121 sample points along the plane main survey line and the connecting line.

[0071] Step 1.14 is specifically as follows: determine the coherent attribute threshold of the fracture zone contour boundary according to the coherent attribute value at the drilled leakage point on the fracture zone. Draw the coherent attribute threshold contour line in the coherent plane attribute map of the top surface of the Ordovician carbonate rock obtained in step 3 according to the coherent attribute threshold as the contour boundary of the strike-slip fracture zone on the top surface of the Ordovician carbonate rock (see Figure 2 ).

[0072] Step 2 is specifically as follows: according to the contour boundary of the strike-slip fault zone on the top surface of the Ordovician carbonate rock obtained in step 8, the line connecting the midpoint of the contour boundary at the beginning of the strike-slip fault and the midpoint of the contour boundary at the end is used as the strike line of the strike-slip fault zone. From the starting point of the strike-slip fault along the strike line of the strike-slip fault zone, a straight line perpendicular to the strike line of the strike-slip fault zone and connecting the contour boundary of the strike-slip fault zone is drawn at intervals of 30m as the lateral deformation width measurement line of the strike-slip fault zone on the top surface of the Ordovician carbonate rock. If there are multiple contour boundaries, the straight line connecting the outermost contour boundaries is used as the lateral deformation width measurement line of the strike-slip fault zone on the top surface of the Ordovician carbonate rock (see Figure 3 ).

[0073] Step 3 is specifically as follows: starting from the starting point of the strike-slip fault, the length of the measurement line of the lateral deformation width of the strike-slip fault zone on the top surface of the Ordovician carbonate rock obtained in step 9 is counted along the strike-slip fault direction, thus obtaining the lateral deformation width data of the top surface of the Ordovician carbonate rock along the strike-slip fault direction (see Figure 5 ).

[0074] Step 4 is specifically as follows: the lateral deformation width coefficient of the top surface of the Ordovician carbonate rock along the strike-slip fault direction is obtained by dividing the lateral deformation width data of the top surface of the Ordovician carbonate rock along the strike-slip fault direction obtained in step 3 by the maximum lateral deformation width data (see Figure 6 ).

[0075] Step 5 is specifically as follows: according to the microstructure map of the top surface of the Ordovician carbonate rock obtained in step 1, within the contour boundary of the strike-slip fault and crushing zone on the top surface of the Ordovician carbonate rock obtained in step 8, the maximum structural axial line is drawn along the strike-slip fault direction from the starting point of the strike-slip fault, which is the axial line of the strike-slip fault and crushing zone on the top surface of the Ordovician carbonate rock (see Figure 4 ).

[0076] Step 6 is specifically as follows: in the microstructure map of the top surface of the Ordovician carbonate rock, from the starting point of the strike-slip fault, along the strike-slip fault direction, sequentially read the structural values ​​at the intersection of the lateral deformation width measurement line of the strike-slip fault zone of the top surface of the Ordovician carbonate rock obtained in step 2 and the axial line of the fault zone of the top surface of the Ordovician carbonate rock obtained in step 5, and use the structural values ​​as the vertical deformation height difference data of the top surface of the Ordovician carbonate rock along the strike-slip fault direction. (See Figure 5 ).

[0077] Step 7 is specifically as follows: the vertical deformation height difference coefficient of the top surface of the Ordovician carbonate rock along the strike-slip fault direction is obtained by dividing the absolute value of the vertical deformation height difference data of the top surface of the Ordovician carbonate rock along the strike-slip fault direction obtained in step 6 by the maximum absolute vertical deformation height difference data (see Figure 6 ).

[0078] Step 8 is specifically as follows: the lateral deformation width coefficient of the top surface of the Ordovician carbonate rock along the strike-slip fault obtained in step 4 is added to the vertical deformation height difference coefficient of the top surface of the Ordovician carbonate rock along the strike-slip fault obtained in step 7 to obtain the activity intensity coefficient of the top surface of the Ordovician carbonate rock along the strike-slip fault (see Figure 6 ).

[0079] Application Example 1

[0080] In 2019, high-density 3D seismic data were re-collected in the south of the Tazhong No. 10 structural belt in the Tazhong No. 1 gas field in the Tarim Basin. High-quality 3D seismic data in the depth domain were obtained through pre-stack depth reverse time migration imaging, laying the foundation for the study of faults and reservoirs. Based on this data, fine structural interpretation and fault identification were carried out. The quantitative evaluation of the intensity of strike-slip fault activity provided a theoretical basis for the classification and segmentation of strike-slip faults, and further combined with reservoir prediction to provide a basis for the division and classification of traps. Finally, three Class I fault-controlled fracture-type pre-exploration traps with an area of ​​118km were submitted. 2 , resource oil 1120×10 4t, gas 86×10 8 m 3 .

[0081] Application Example 2

[0082] By reprocessing the existing 3D seismic data in the periphery of Block III of the Tazhong I gas field in the Tarim Basin, the accuracy of fault and reservoir identification has been greatly improved. Through fine structural interpretation, a strong quantitative evaluation of fault-traversing activities was carried out. Combined with reservoir prediction, a total of 2 Class I fault-controlled fracture-cavity pre-exploration traps with an area of ​​18.9 km were submitted. 2 , oil resources 513×10 4 t, gas 7×10 8 m 3 .

Claims

1. A quantitative evaluation method for the activity intensity of strike-slip faults, characterized in that, it includes the following steps: Step 1: Obtain the micro-structure map of the Ordovician carbonate rock top surface and the contour boundary of the strike-slip fault fracture zone on the Ordovician carbonate rock top surface based on high-precision 3D seismic data in the depth domain; Step 2: Use the contour boundary of the strike-slip fault fracture zone on the Ordovician carbonate rock top surface obtained in Step 1 to obtain the lateral deformation width measurement line of the strike-slip fault fracture zone on the Ordovician carbonate rock top surface; Step 3: Use the lateral deformation width measurement line of the strike-slip fault fracture zone on the Ordovician carbonate rock top surface obtained in Step 2 to obtain the lateral deformation width data of the Ordovician carbonate rock top surface along the strike of the strike-slip fault; Step 4: Use the lateral deformation width data of the Ordovician carbonate rock top surface along the strike of the strike-slip fault obtained in Step 3 to obtain the lateral deformation width coefficient of the Ordovician carbonate rock top surface along the strike of the strike-slip fault; Step 5: Use the micro-structure map of the Ordovician carbonate rock top surface and the contour boundary of the fracture zone on the Ordovician carbonate rock top surface obtained in Step 1 to obtain the axial line of the strike-slip fault fracture zone on the Ordovician carbonate rock top surface; Step 6: Use the micro-structure map of the Ordovician carbonate rock top surface obtained in Step 1, the lateral deformation width measurement line of the strike-slip fault fracture zone on the Ordovician carbonate rock top surface obtained in Step 2, and the axial line of the strike-slip fault fracture zone on the Ordovician carbonate rock top surface obtained in Step 5 to obtain the vertical deformation height difference data of the Ordovician carbonate rock top surface along the strike of the strike-slip fault; Step 7: Use the vertical deformation height difference data of the Ordovician carbonate rock top surface along the strike of the strike-slip fault obtained in Step 6 to obtain the vertical deformation height difference coefficient of the Ordovician carbonate rock top surface along the strike of the strike-slip fault; Step 8: Use the lateral deformation width coefficient of the Ordovician carbonate rock top surface along the strike of the strike-slip fault obtained in Step 4 and the vertical deformation height difference coefficient of the Ordovician carbonate rock top surface along the strike of the strike-slip fault obtained in Step 7 to obtain the activity intensity coefficient of the Ordovician carbonate rock top surface along the strike of the strike-slip fault; Judge the activity intensity of the strike-slip fault according to the activity intensity coefficient of the Ordovician carbonate rock top surface along the strike of the strike-slip fault.

2. The quantitative evaluation method for the activity intensity of strike-slip faults according to claim 1, characterized in that, obtaining the micro-structure map of the Ordovician carbonate rock top surface based on high-precision 3D seismic data in the depth domain includes the following steps: 1.1: Carry out fine interpretation of the Ordovician carbonate rock top surface layer based on high-precision 3D seismic data in the depth domain to obtain the Ordovician carbonate rock top surface layer; 1.2: Process the Ordovician carbonate rock top surface layer to obtain the micro-structure layer of the Ordovician carbonate rock top surface; 1.3: Use the micro-structure layer of the Ordovician carbonate rock top surface to obtain the micro-structure map of the Ordovician carbonate rock top surface; Obtaining the contour boundary of the strike-slip fault fracture zone on the Ordovician carbonate rock top surface based on high-precision 3D seismic data in the depth domain includes the following steps: 1.11: Extract the coherent body attributes according to the high-precision 3D seismic data in the depth domain; 1.12: Obtain the coherent plane attributes of the Ordovician carbonate rock top surface according to the coherent body attributes extracted in Step 1.11 and the Ordovician carbonate rock top surface layer obtained in Step 1.1; 1.

13. Using the coherent plane attributes of the top surface of the Ordovician carbonate rock, the coherent plane attribute map of the top surface of the Ordovician carbonate rock is obtained; 1.

14. Use the coherent plane attribute map of the top surface of the Ordovician carbonate rock to obtain the contour boundary of the strike-slip fault zone on the top surface of the Ordovician carbonate rock.

3. A method for quantitatively evaluating the intensity of strike-slip fault activity according to claim 2, It is characterized in that Step 1.2 is specifically as follows: smoothing the top surface layer of the Ordovician carbonate rock to obtain a smooth background layer; The microstructural horizons of the top surface of the Ordovician carbonate rock are obtained by subtracting the smooth background horizons of the top surface of the Ordovician carbonate rock from the top surface of the Ordovician carbonate rock.

4. A method for quantitatively evaluating the intensity of strike-slip fault activity according to claim 2, It is characterized in that In step 1.11, when extracting coherent plane attributes, the vertical window parameter setting is centered on the top surface of the Ordovician carbonate rock layer and includes a range of 20m above and below as the vertical extraction time window for coherent plane attributes.

5. A method for quantitatively evaluating the intensity of strike-slip fault activity according to claim 2, It is characterized in that Step 1.14 is specifically as follows: determining the coherent attribute threshold of the boundary of the fracture zone contour according to the coherent attribute value at the drilled leakage point on the fracture zone; According to the coherent attribute threshold of the fracture zone contour boundary, in the coherent plane attribute map of the top surface of the Ordovician carbonate rock, the coherent attribute threshold contour line is drawn as the contour boundary of the strike-slip fracture zone on the top surface of the Ordovician carbonate rock.

6. A method for quantitatively evaluating the intensity of strike-slip fault activity according to claim 1, It is characterized in that Step 2 is specifically as follows: according to the contour boundary of the strike-slip fault zone on the top surface of the Ordovician carbonate rock obtained in step 1, a line connecting the midpoint of the contour boundary at the beginning of the strike-slip fault and the midpoint of the contour boundary at the end is used as the strike line of the strike-slip fault zone; From the starting point of the strike-slip fault, along the strike line of the strike-slip fault zone, a straight line perpendicular to the strike line of the strike-slip fault zone and connecting the contour boundary of the strike-slip fault zone is drawn at intervals of 30m as the lateral deformation width measurement line of the strike-slip fault zone on the top surface of the Ordovician carbonate rock; If there are multiple contour boundaries, the straight line connecting the outermost contour boundaries will be used as the measurement line of the lateral deformation width of the strike-slip fault zone on the top surface of the Ordovician carbonate rock.

7. A method for quantitatively evaluating the intensity of strike-slip fault activity according to claim 1, It is characterized in that Step 3 is specifically as follows: starting from the starting point of the strike-slip fault, the length of the lateral deformation width measurement line of the strike-slip fault zone on the top surface of the Ordovician carbonate rock obtained in step 2 is sequentially counted along the strike-slip fault direction to obtain the lateral deformation width data of the top surface of the Ordovician carbonate rock along the strike-slip fault direction.

8. A method for quantitatively evaluating the intensity of strike-slip fault activity according to claim 1, It is characterized in that Step 4 is specifically as follows: the lateral deformation width data of the top surface of the Ordovician carbonate rock along the strike-slip fault obtained in step 3 is divided by the maximum lateral deformation width data to obtain the lateral deformation width coefficient of the top surface of the Ordovician carbonate rock along the strike-slip fault; Step 7 is specifically as follows: the absolute value of the vertical deformation height difference data of the top surface of the Ordovician carbonate rock along the strike-slip fault obtained in step 6 is divided by the maximum absolute vertical deformation height difference data, so as to obtain the vertical deformation height difference coefficient of the top surface of the Ordovician carbonate rock along the strike-slip fault.

9. A method for quantitatively evaluating the intensity of strike-slip fault activity according to claim 1, It is characterized in that Step 5 is specifically as follows: in the microstructure map of the top surface of the Ordovician carbonate rock obtained in step 1, and within the contour boundary of the strike-slip fault and crushing zone on the top surface of the Ordovician carbonate rock obtained in step 1, the maximum structural axial line is drawn along the strike-slip fault direction from the starting point of the strike-slip fault to obtain the axial line of the strike-slip fault and crushing zone on the top surface of the Ordovician carbonate rock.

10. A method for quantitatively evaluating the intensity of strike-slip fault activity according to claim 1, It is characterized in that Step 6 is specifically as follows: in the microstructure map of the top surface of the Ordovician carbonate rock, starting from the starting point of the strike-slip fault, along the strike-slip fault direction, sequentially read the structural values ​​at the intersection of the lateral deformation width measurement line of the strike-slip fault zone on the top surface of the Ordovician carbonate rock obtained in step 2 and the axial line of the fault and crushing zone on the top surface of the Ordovician carbonate rock obtained in step 5, and use the structural values ​​as the vertical deformation height difference data of the top surface of the Ordovician carbonate rock along the strike-slip fault direction.

Images