Method for determining strike-slip fracture relative activity intensity based on fracture piece length

Abstract

The invention discloses a method for determining strike-slip fracture relative activity intensity based on fracture piece length, and relates to the technical field of tectonic geology. Three-dimensional seismic data of a target strike-slip fracture zone are acquired; performing three-dimensional fine analysis on the strike-slip fault zone according to the three-dimensional seismic data, and identifying all fracture slices; determining the horizon of the strike-slip fracture zone according to the three-dimensional seismic data, and drawing a strike-slip fracture plane distribution diagram according to the horizon; according to the strike-slip fracture plane distribution diagram, numbering all fracture pieces; the lengths of all fracture pieces are measured according to the strike-slip fracture plane distribution diagram, and normalization processing is carried out; drawing a broken line graph according to the number of the fracture piece and the length of the fracture piece after normalization processing; and determining the relative activity intensity of the strike-slip fault zone according to the broken line graph. According to the method, the activity intensity of the strike-slip fracture system is determined by adopting the length of the fracture piece, so that the relative activity intensity of the strike-slip fracture can be accurately judged, the whole method is simple, and the estimation cost is reduced.

Description

Technical Field

[0001] This application relates to the technical field of structural geology, and specifically relates to a method for determining the relative activity intensity of strike-slip faults based on the length of fracture fragments. Background Art

[0002] Oil exploration is a crucial task in the field of energy development. It explores underground oil resources through scientific methods and technical means. With the continuous growth of global energy demand, oil exploration not only plays a key role in the economic development of each country, but also faces increasingly complex environmental and technical challenges.

[0003] Previous oil exploration was often carried out in basin thrust belts and uplift zones (e.g., the South Tianshan fold belt, the Tabei uplift, and the Tazhong uplift). However, with the increase in the scope and depth of oil and gas exploration, a large number of oil and gas accumulations have also been discovered within the basin, and these oil and gas accumulations are usually near strike-slip fault zones (e.g., the Tahe Oilfield and the Shunbei Oil and Gas Field). In recent years, strike-slip fault zones at the scale of cratonic basins have been identified in the three major cratonic basins in China. In the Tarim Basin, affected by the multi-stage activities of the strike-slip fault system, large-scale reservoirs have been formed in Ordovician carbonate rocks in multiple areas such as Tazhong, Shunbei, Tuoputai, and Halaha Tang.

[0004] With the continuous deepening of the understanding of the relationship between strike-slip fault systems and oil and gas, some scholars have pointed out that strike-slip fault systems play the roles of "controlling the source, transportation, storage, trapping, accumulation, and enrichment". In the study of the relationship between strike-slip faults and oil and gas, the activity intensity of strike-slip faults is an important part that cannot be ignored. Generally speaking, the stronger the activity of strike-slip faults, the stronger the transformation effect on the reservoir, and it is more likely to form a good reservoir. Therefore, determining the activity intensity of strike-slip faults has a very important impact on oil and gas exploration.

[0005] Currently, the methods for determining the relative activity intensity of strike-slip faults mainly include: (1) multi-layer vertical deformation amplitude analysis; (2) single-layer strike-slip fault zone width analysis; (3) deformation volume analysis of strike-slip fault overlapping segments; (4) quantitative analysis of horizontal slip distance using the method of changes in formation thickness on both sides of strike-slip faults. The above methods mainly rely on 3D seismic data and fracture analysis techniques to determine the activity intensity of strike-slip fault systems. However, the above several methods are relatively complex and it is difficult to accurately determine the relative activity intensity of strike-slip faults. Summary of the Invention

[0006] Therefore, this application provides a method for determining the relative activity intensity of strike-slip faults based on the length of fracture fragments to solve the problem that it is difficult to accurately determine the relative activity intensity of strike-slip faults in the prior art.

[0007] To achieve the above object, this application provides the following technical solutions:

[0008] In a first aspect, a method for determining the relative activity intensity of a strike-slip fault based on the length of fracture fragments includes:

[0009] Step 1: Obtain 3D seismic data of the target strike-slip fault zone;

[0010] Step 2: Conduct 3D fine analysis on the strike-slip fault zone according to the 3D seismic data, and identify all fracture fragments;

[0011] Step 3: Determine the horizons of the strike-slip fault zone according to the 3D seismic data, and draw a planar distribution map of the strike-slip fault based on the horizons;

[0012] Step 4: Number all fracture fragments according to the planar distribution map of the strike-slip fault;

[0013] Step 5: Measure the lengths of all fracture fragments according to the planar distribution map of the strike-slip fault, and perform normalization processing;

[0014] Step 6: Draw a line graph according to the fracture fragment numbers and the lengths of the fracture fragments after normalization processing;

[0015] Step 7: Determine the relative activity intensity of the strike-slip fault zone according to the line graph.

[0016] Preferably, Step 2 specifically includes:

[0017] Step 201: Conduct fracture interpretation on the strike-slip fault zone in the work area according to the 3D seismic data;

[0018] Step 202: Determine the strike of the strike-slip fault zone according to the fracture interpretation;

[0019] Step 203: Make a profile along the direction perpendicular to the strike of the strike-slip fault zone, and conduct fine interpretation on the fractures;

[0020] Step 204: Combine the fractures in 3D space according to the fine interpretation, and identify all fracture fragments.

[0021] Preferably, Step 204 specifically includes: Combining the fractures on each profile according to the fine interpretation through the change of the fracture penetration depth, the change of the fracture strike on the plane, and the different connection methods of the fractures in space, and dividing them into different fracture fragments.

[0022] Preferably, Step 3 specifically includes: Conduct horizon tracking according to the 3D seismic data, determine all horizons penetrated by the strike-slip fault zone according to the depth penetrated by the strike-slip fault zone, select appropriate horizons from all horizons, and draw a planar distribution map of the strike-slip fault.

[0023] Preferably, step 4 specifically includes: numbering from one end to the other end of the strike-slip fault along the strike of the strike-slip fault according to the strike-slip fault plane distribution map.

[0024] Preferably, in step 5, the minimum-maximum normalization method is used for the normalization process.

[0025] Preferably, step 6 specifically includes: taking the length of the fractured fragment after normalization as the Y-axis of the line graph. For the X-axis, mark the starting point of the strike-slip fault as 0, the ending point as the total length value of the strike-slip fault, and the distance from the midpoint of each fractured fragment to the starting point as D i , D i is the position of the corresponding fractured fragment on the X-axis. At this time, project the corresponding points onto the line graph and connect all the points with straight lines to obtain the line graph.

[0026] Preferably, step 7 is specifically: determining the relative activity intensity of the strike-slip fault zone according to the upward and downward change trend of the line graph.

[0027] In a second aspect, a device for determining the relative activity intensity of a strike-slip fault based on the length of fractured fragments includes:

[0028] A data acquisition module for acquiring three-dimensional seismic data of a target strike-slip fault zone;

[0029] A fractured fragment identification module for performing three-dimensional fine analysis on the strike-slip fault zone according to the three-dimensional seismic data and identifying all fractured fragments;

[0030] A strike-slip fault plane distribution map drawing module for determining the horizon of the strike-slip fault zone according to the three-dimensional seismic data and drawing a strike-slip fault plane distribution map according to the horizon;

[0031] A numbering module for numbering all fractured fragments according to the strike-slip fault plane distribution map;

[0032] A length determination module for measuring the lengths of all fractured fragments according to the strike-slip fault plane distribution map and performing normalization processing;

[0033] A line graph drawing module for drawing a line graph according to the fractured fragment numbers and the lengths of the fractured fragments after normalization processing;

[0034] A relative activity intensity determination module for determining the relative activity intensity of the strike-slip fault zone according to the line graph.

[0035] In a third aspect, a computer device includes a memory and a processor. The memory stores a computer program, and when the processor executes the computer program, the steps of a method for determining the relative activity intensity of a strike-slip fault based on the length of fractured fragments are implemented.

[0036] Compared with the prior art, the present application has at least the following beneficial effects:

[0037] The present application provides a method for determining the relative activity intensity of strike-slip faults based on the length of fracture fragments. The method includes obtaining three-dimensional seismic data of the target strike-slip fault zone; performing three-dimensional fine analysis on the strike-slip fault zone according to the three-dimensional seismic data and identifying all fracture fragments; determining the horizons of the strike-slip fault zone according to the three-dimensional seismic data and drawing a plane distribution map of the strike-slip fault based on the horizons; numbering all fracture fragments according to the plane distribution map of the strike-slip fault; measuring the lengths of all fracture fragments according to the plane distribution map of the strike-slip fault and performing normalization processing; drawing a line graph according to the fracture fragment numbers and the normalized fracture fragment lengths; and determining the relative activity intensity of the strike-slip fault zone according to the line graph. The present application uses the length of fracture fragments to determine the activity intensity of the strike-slip fault system, which can not only accurately judge the relative activity intensity of the strike-slip fault, but also has a simple and convenient method, reducing the calculation cost. BRIEF DESCRIPTION OF THE DRAWINGS

[0038] To more intuitively illustrate the prior art and the present application, exemplary drawings are given below. It should be understood that the specific shapes and structures shown in the drawings generally should not be regarded as limiting conditions when implementing the present application; for example, those skilled in the art are capable of making conventional adjustments or further optimizations to the addition / deletion / attribution division of certain units (components), specific shapes, positional relationships, connection methods, dimensional proportional relationships, etc. based on the technical concept disclosed in the present application and the exemplary drawings.

[0039] Figure 1 It is a flowchart of a method for determining the relative activity intensity of strike-slip faults based on the length of fracture fragments provided in Embodiment 1 of the present application;

[0040] Figure 2 It is a plane distribution map of the strike-slip fault profile interpretation provided in Embodiment 1 of the present application;

[0041] Figure 3 It is a strike-slip fault profile interpretation diagram provided in Embodiment 1 of the present application;

[0042] Figure 4 It is a plane distribution map of the strike-slip fault provided in Embodiment 1 of the present application;

[0043] Figure 5 It is a line graph of "fracture fragment length normalization - numbering" of the strike-slip fault provided in Embodiment 1 of the present application;

[0044] Figure 6 It is a line graph of "relative activity intensity of strike-slip fault - numbering" of the strike-slip fault provided in Embodiment 1 of the present application. DETAILED DESCRIPTION OF THE EMBODIMENTS

[0045] The following further details the present application through specific embodiments in conjunction with the accompanying drawings.

[0046] In the description of the present application: Unless otherwise specified, the meaning of "a plurality" is two or more. Terms such as "first", "second", "third", etc. in the present application are intended to distinguish the objects being referred to, without any special meaning in terms of technical connotations (for example, it should not be understood as emphasizing importance or order, etc.). Expressions such as "including", "comprising", "having", etc. also mean "not limited to" (certain units, components, materials, steps, etc.).

[0047] Terms such as "upper", "lower", "left", "right", "middle", etc. cited in the present application are usually indications of the general relative position relationship for the convenience of intuitively understanding with reference to the accompanying drawings, and are not absolute limitations on the position relationship in the actual product.

[0048] Embodiment 1

[0049] This embodiment provides a method for calibrating the relative activity intensity of strike-slip faults based on the length of fracture fragments, and uses a new geometric feature, that is, the length of fracture fragments to calibrate the activity intensity of the strike-slip fault system. Because the length L of the fracture fragment and the maximum displacement D of the fault max have a functional relationship: D max = cL n , where the value range of n is between 0.5 and 2. From the above formula, it can be obtained that the length L of the fracture fragment and the maximum displacement D of the fault max are positively correlated, that is, the longer the length of the fracture fragment, the larger the maximum displacement D of the fault max is, and the stronger the fault activity is. Therefore, the calibration of fault activity can be characterized by using the length of fracture fragments.

[0050] Please refer to Figure 1 , a method for calibrating the relative activity intensity of strike-slip faults based on the length of fracture fragments provided in this embodiment includes:

[0051] S1: Obtain the three-dimensional seismic data of the target strike-slip fault zone;

[0052] Specifically, through the three-dimensional seismic data, the morphology, scale, and distribution of the strike-slip fault zone can be accurately depicted. Compared with the traditional two-dimensional seismic data, the three-dimensional seismic data provides more detailed geological information and can better understand the geometric characteristics of the fault zone (such as strike, dip, depth, etc.).

[0053] S2: Conduct three-dimensional fine analysis on the strike-slip fault zone according to the three-dimensional seismic data, and identify all fracture fragments;

[0054] This step conducts three-dimensional fine analysis of the faults on the target strike-slip fault zone with three-dimensional seismic data, specifically including:

[0055] S201: Conduct fracture interpretation on the strike-slip fault zone in the work area based on 3D seismic data;

[0056] In this step, on the cross-section, conduct fracture interpretation on the strike-slip faults in the work area.

[0057] S202: Determine the trend of the strike-slip fault zone based on the fracture interpretation;

[0058] In this step, determine the general trend of the strike-slip fault based on the interpreted fractures in the work area.

[0059] S203: Make a cross-section along the trend perpendicular to the strike-slip fault zone and conduct fine interpretation on the fractures;

[0060] This step requires modifying it along the general trend perpendicular to the interpreted strike-slip fault to make the interpreted fractures more in line with the actual geological situation. In the embodiment, the interpreted cross-section should be consistent with the Figure 2 blue horizontal line shown in, and make more detailed modifications to the fractures accordingly.

[0061] S204: Combine in 3D space according to the fine interpretation and identify all fracture segments.

[0062] This step requires combining the interpreted fractures in 3D space. Combine the fractures on each cross-section through the change of the fracture penetration depth, the change of the fracture trend on the plane, and different connection methods (soft connection, hard connection) of the fractures in space, and divide them into different fracture segments. In the embodiment, some fractures are connected in space, but they should be hard connections, not one fracture, and are divided into two or more fracture segments.

[0063] S3: Determine the horizons of the strike-slip fault zone based on 3D seismic data and draw a plane distribution map of the strike-slip faults according to the horizons;

[0064] This step specifically includes:

[0065] S301: Conduct horizon tracking based on 3D seismic data and determine all horizons penetrated by the strike-slip fault zone according to the depth penetrated by the strike-slip fault zone;

[0066] This step requires conducting horizon tracking on the seismic data in the work area according to the geological background. Vertically, determine the horizons penetrated by it according to the interpreted strike-slip faults and the results after horizon tracking, laying a foundation for determining the appropriate horizons in the next step. In the embodiment, the horizons in the work area are divided into 12 layers, which are respectively Among them, the main trunk of the strike-slip fault penetrates 8 layers, which are respectively Such as Figure 3 shown.

[0067] S302: Select appropriate horizons from all horizons;

[0068] In this step, when selecting appropriate horizons, horizons that are penetrated by the largest number of faults in the same fault segment should be selected as much as possible. For each fault segment, the same horizon should be selected. Such a selection can provide a relatively accurate data basis when determining the relative activity intensity; for strike-slip faults, the shallowest horizon penetrated by the strike-slip fault, that is, the last penetrated horizon, should be selected.

[0069] S303: Draw a strike-slip fault plane distribution map based on the selected appropriate horizons.

[0070] In the embodiment, horizon T 7 4 is selected to draw a strike-slip fault plane distribution map, which is divided into 24 fault segments in total, as Figure 4 shown.

[0071] S4: Number all fault segments according to the strike-slip fault plane distribution map;

[0072] Specifically, when numbering the fault segments in this step, they should be numbered in sequence, that is, along the strike of the strike-slip fault, from one end of the strike-slip fault to the other end, and the number is denoted as i. In the embodiment, in the strike-slip fault plane distribution map, the fault segments are numbered from south to north as B1 to B24, as Figure 4 shown.

[0073] S5: Measure the lengths of all fault segments according to the strike-slip fault plane distribution map and perform normalization processing;

[0074] Specifically, in this step, first, according to the drawn strike-slip fault plane distribution map, measure the length of each fault segment, and the measured length is denoted as L i , with the unit of m. Among them, the maximum length is denoted as L max , and the minimum length is denoted as L min . Then record the measured results in a table and perform normalization processing on it. The normalization processing result is L i0-1 , dimensionless.

[0075] More specifically, when measuring the length of the fault segment in this step, the length change caused by the depth bending of the fault due to terrain changes does not need to be considered. Only the length of the fault segment on the plane needs to be measured, that is, on a horizon plane, only the length of the fault segment in the X and Y axis planes is considered, while the length change of the fault segment in the Z axis direction is ignored.

[0076] In the embodiment, by measuring the length of the fault segment in the X and Y axis planes, the length of each fault segment is obtained, which are respectively: L B1 : 1239.2m, L B2: 3181.29 m, L B3 : 7108.59 m, L B4 : 6999.23 m, L B5 : 9683.24 m, L B6 : 836.74 m, L B7 : 3565.93 m, L B8 : 4273.62, L B9 : 2597.62 m, L B10 : 4933.66, L B11 : 6922.28 m, L B12 : 5247.06 m, L B13 : 836.12 m, L B14 : 2683.69 m, L B15 : 5836.62 m, L B16 : 4227.87 m, L B17 : 8470.37 m, L B18 : 25890.81 m, L B19 : 8547.04 m, L B20 : 11037.12 m, L B21 : 14401.55 m, L B22 : 9355.41 m, L B23 : 5151.82 m, L B24 : 6396.18 m.

[0077] When normalizing the length of the fracture fragment, the min-max normalization method is adopted, and the formula is:

[0078]

[0079] Among them, L i is the length of the fracture fragment, with the unit of m; L max is the maximum length of the fracture fragments of this strike-slip fault, with the unit of m; L min is the minimum length of the fracture fragments of this strike-slip fault, with the unit of m; L i0-1 is the result after normalizing the corresponding fracture fragment length, dimensionless.

[0080] In the embodiment, through the above formula, min-max normalization processing is performed on each fracture fragment, and the processed results are respectively: L B1 : 0.0161, L B2 : 0.0936, L B3 : 0.2504, L B4 : 0.2460, L B5 : 0.3531, L B6 : 0, L B7: 0.1090, L B8 : 0.1372, L B9 : 0.0703, L B10 : 0.1635, L B11 : 0.2429, L B12 : 0.1761, L B13 : 0, L B14 : 0.0737, L B15 : 0.1996, L B16 : 0.1354, L B17 : 0.3047, L B18 : 1, L B19 : 0.3078, L B20 : 0.4071, L B21 : 0.5414, L B22 : 0.3400, L B23 : 0.1723, L B24 : 0.2219.

[0081] S6: Draw a line graph based on the fracture fragment number and the length of the fracture fragment after normalization;

[0082] Specifically, when drawing the line graph, use the result L of the normalization of the fracture fragment length i0-1 as the Y-axis of the line graph. For the X-axis, mark the starting point of the strike-slip fault as 0, the ending point as the total length value of the strike-slip fault, and the distance from the midpoint of each fracture fragment to the starting point as D i , and at this time D i is the position of the corresponding fracture fragment on the X-axis. At this time, project the corresponding points onto the line graph and connect all the points with a straight line. In the embodiment, the total length of the strike-slip fault is m, and the "fracture fragment length - distance" line graph is obtained by the above method, as Figure 5 shown.

[0083] S7: Determine the relative activity intensity of the strike-slip fault zone according to the line graph.

[0084] Specifically, for the "fault segment length - distance" broken line graph after drawing, it is considered that the up and down change trend of the broken line is the trend of the relative activity intensity change of the strike-slip fault, and the Y value of each corresponding point is the relative activity intensity of the strike-slip fault. The larger the value, the greater the activity intensity. After obtaining the broken line graph, the points on the X-axis representing the positions of each fault segment in the broken line graph are corresponding to the positions of the fault segments in the strike-slip fault plane distribution map. At this time, through the trend of the curve in the "fault segment length - distance" broken line graph, it is easy to obtain the relative activity intensity of each position on the strike-slip fault. In the embodiment, the obtained "fault segment length - distance" broken line graph is projected onto the strike-slip fault plane distribution map to obtain the relative activity intensity of the strike-slip fault. The rule is that the activity intensity is high in the middle and north, and low in the south, north, and branch faults, as Figure 6 shown.

[0085] The method for determining the relative activity intensity of strike-slip faults based on fault segment length provided in this embodiment can accurately determine the relative activity intensity of strike-slip faults. The method is simple and can be widely applied to the determination of the relative activity intensity of underground strike-slip faults. And only three-dimensional seismic data is required, and the estimation cost is low, with high operability.

[0086] It should be noted that this embodiment is mainly based on the influence of the strike-slip fault system activity on the fault segment length. In the strike-slip fault system, the position with greater activity intensity often has a greater maximum vertical displacement of the fault, and the fault segment length L and the maximum vertical displacement D of the fault max have a functional relationship, and the formula is D max = cL n . Therefore, there is a positive correlation between the fault segment length and the maximum vertical displacement of the fault. By analyzing the fault segment length in the strike-slip fault system, the relative activity intensity of the strike-slip fault can be obtained. Therefore, it is feasible and reasonable to quantitatively determine the relative activity intensity of the strike-slip fault using the fault segment length.

[0087] Embodiment Two

[0088] This embodiment provides a device for determining the relative activity intensity of strike-slip faults based on fault segment length, including:

[0089] A data acquisition module for acquiring three-dimensional seismic data of the target strike-slip fault zone;

[0090] A fault segment identification module for performing three-dimensional fine analysis on the strike-slip fault zone according to the three-dimensional seismic data and identifying all fault segments;

[0091] A strike-slip fault plane distribution map drawing module for determining the horizons of the strike-slip fault zone according to the three-dimensional seismic data and drawing a strike-slip fault plane distribution map according to the horizons;

[0092] A numbering module, configured to number all fracture segments according to the strike-slip fault plane distribution map;

[0093] A length determination module, configured to measure the lengths of all fracture segments according to the strike-slip fault plane distribution map and perform normalization processing;

[0094] A line graph drawing module, configured to draw a line graph according to the fracture segment numbers and the normalized lengths of the fracture segments;

[0095] A relative activity intensity determination module, configured to determine the relative activity intensity of the strike-slip fault zone according to the line graph.

[0096] For the specific implementation content of each module in a device for determining the relative activity intensity of a strike-slip fault based on the length of fracture segments, reference may be made to the limitations on the method for determining the relative activity intensity of a strike-slip fault based on the length of fracture segments in the foregoing text, and details are not described herein again.

[0097] Embodiment III

[0098] This embodiment provides a computer device, including a memory and a processor. The memory stores a computer program, and when the processor executes the computer program, the steps of a method for determining the relative activity intensity of a strike-slip fault based on the length of fracture segments are implemented.

[0099] The technical features of the above embodiments can be combined arbitrarily (as long as there is no contradiction in the combination of these technical features). For the sake of brevity of description, not all possible combinations of the technical features in the above embodiments are described; these embodiments that are not explicitly described should also be considered as within the scope described in this specification.

Claims

1. A method for determining the relative activity intensity of a strike-slip fault based on the length of a fault segment, characterized in that: include: Step 1: Acquire 3D seismic data of the target strike-slip fault zone; Step 2: Performing a three-dimensional fine analysis of the strike-slip fault zone according to the three-dimensional seismic data, and identifying all fault fragments; Step 3: determining the horizon of the strike-slip fault zone according to the three-dimensional seismic data, and drawing a strike-slip fault plane distribution map according to the horizon; Step 4: numbering all fault fragments according to the strike-slip fault plane distribution map; Step 5: measuring the lengths of all fault fragments according to the strike-slip fault plane distribution map and performing normalization processing; Step 6: Draw a line graph based on the fragment number and the normalized fragment length; Step 7: Determine the relative activity intensity of the strike-slip fault zone according to the line graph.

2. The method for determining the relative activity intensity of a strike-slip fault based on the length of a fault segment according to claim 1, characterized in that: Step 2 specifically includes: Step 201: performing fault interpretation on the strike-slip fault zone in the work area according to the three-dimensional seismic data; Step 202: determining the strike direction of the strike-slip fault zone according to the fault interpretation; Step 203: making a cross section along the direction perpendicular to the strike-slip fault zone, and performing a detailed interpretation of the fault; Step 204: Combine in three-dimensional space according to the detailed interpretation and identify all fractured pieces.

3. The method for determining the relative activity intensity of a strike-slip fault based on the length of a fault segment according to claim 2, characterized in that: Step 204 specifically includes: combining the faults on each section and dividing them into different fault slices according to the detailed interpretation through the change of the fault penetration depth, the change of the fault orientation on the plane and the different connection modes of the faults in space.

4. The method for determining the relative activity intensity of a strike-slip fault based on the length of a fault segment according to claim 1, characterized in that: Step 3 specifically includes: tracking the layers according to the three-dimensional seismic data, and determining all the layers of the strike-slip fault zone according to the depth of the strike-slip fault zone, selecting appropriate layers from all the layers, and drawing a plane distribution map of the strike-slip fault.

5. The method for determining the relative activity intensity of a strike-slip fault based on the length of a fault segment according to claim 1, characterized in that: Step 4 specifically includes: numbering the strike-slip faults from one end to the other end along the strike-slip fault direction according to the strike-slip fault plane distribution map.

6. The method for determining the relative activity intensity of a strike-slip fault based on the length of a fault segment according to claim 1, characterized in that: In step 5, the minimum and maximum value normalization method is used for the normalization process.

7. The method for determining the relative activity intensity of a strike-slip fault based on the length of a fault segment according to claim 1, characterized in that: Step 6 specifically includes: taking the normalized length of the fault piece as the Y axis of the line graph, for the X axis, the starting point of the strike-slip fault is recorded as 0, the end point is the total length of the strike-slip fault, and the distance from the midpoint of each fault piece to the starting point is recorded as D i , D i is the position of the corresponding fracture piece on the X-axis. At this time, the corresponding point is projected on the line graph, and all the points are connected with straight lines to obtain a line graph.

8. The method for determining the relative activity intensity of a strike-slip fault based on the length of a fault segment according to claim 1, characterized in that: Step 7 specifically includes: determining the relative activity intensity of the strike-slip fault zone according to the up and down changing trend of the line graph.

9. A device for determining the relative activity intensity of a strike-slip fault based on the length of a fault segment, characterized in that: include: A data acquisition module, used to acquire three-dimensional seismic data of a target strike-slip fault zone; A fault piece identification module, used for performing a three-dimensional fine analysis of the strike-slip fault zone according to the three-dimensional seismic data, and identifying all fault pieces; A strike-slip fault plane distribution map drawing module is used to determine the horizon of the strike-slip fault zone according to the three-dimensional seismic data, and draw a strike-slip fault plane distribution map according to the horizon; A numbering module, used for numbering all fault pieces according to the strike-slip fault plane distribution map; A length determination module, used to measure the lengths of all fractured pieces according to the strike-slip fracture plane distribution map and perform normalization processing; A line graph drawing module is used to draw a line graph according to the fracture piece number and the normalized fracture piece length; A relative activity intensity determination module is used to determine the relative activity intensity of the strike-slip fault zone according to the broken line graph.

10. A computer device comprising a memory and a processor, wherein the memory stores a computer program, characterized in that: When the processor executes the computer program, the steps of the method according to any one of claims 1 to 8 are implemented.

Images