Fault activity period determination method and device

Abstract

The application discloses a fault activity period determination method and device, and the method comprises the following steps: determining the position of a target fault and the stratum range of the target fault based on obtained seismic data; determining the standard interface of the stratum range of the target fault and a plurality of preset key interfaces which do not penetrate the stratum range of the target fault according to the position of the target fault and the stratum range of the target fault; performing curve sampling on the standard interface and the plurality of preset key interfaces, calculating the similarity coefficient of each preset key interface and the standard interface based on the sampling point set of the standard interface and the sampling point set of the plurality of preset key interfaces; comparing the similarity coefficients between adjacent preset key interfaces, and determining whether the preset key interfaces are affected by fault activity according to the comparison result; converting the geological time corresponding to the preset key interfaces affected by the fault activity into time sequence data, and determining the fault activity period according to the time sequence data. The application can improve the accuracy of fault activity period determination.

Description

Technical Field

[0001] This invention relates to the field of oil and gas exploration technology, and in particular to a method and apparatus for determining fault activity phases. Background Technology

[0002] This section is intended to provide background or context for the embodiments of the invention set forth in the claims. The description herein is not an admission that it is prior art simply because it is included in this section.

[0003] Large oil and gas basins contain various types of faults, which are significant for hydrocarbon accumulation, influencing reservoir development, trap preservation conditions, and hydrocarbon migration pathways. Normal faults are numerous and their formation stages are complex. Determining the fault's activity stage is a crucial task in structural analysis and geological exploration, and a prerequisite for confirming the fault's role in controlling hydrocarbon growth. Current techniques only consider the influence of faults on penetrating strata, neglecting the potential influence on non-penetrating strata, leading to inaccurate determinations of fault activity stages. Therefore, a method for determining fault activity stages is urgently needed. Summary of the Invention

[0004] This invention provides a method for determining the phases of fault activity, addressing the problem in existing technologies that do not consider the influence of fault activity on non-penetrating strata, thereby improving the accuracy of fault activity phase determination. The method includes:

[0005] Based on the acquired seismic data, the location of the target fault and the extent of the strata it penetrates are determined.

[0006] Based on the location and stratigraphic extent of the target fault, the standard interface of the strata that the fault penetrates and several preset key interfaces of the strata that it does not penetrate are determined. The standard interface is the strata within the penetrating stratigraphic extent that are affected by fault activity; the multiple preset key interfaces are the strata outside the penetrating stratigraphic extent.

[0007] Curve sampling is performed on the standard interface and multiple preset standard interfaces to obtain the sampling point set of the standard interface and the sampling point set of multiple preset key interfaces respectively. Based on the sampling point set of the standard interface and the sampling point set of multiple preset key interfaces, the similarity coefficient between each preset key interface and the standard interface is calculated.

[0008] The similarity coefficients between adjacent preset standard interfaces are compared, and the preset key interfaces are determined to be affected by fault activity based on the comparison results.

[0009] The geological time corresponding to the preset key interface affected by fault activity is converted into time series data, and the fault activity period is determined based on the time series data.

[0010] This invention also provides a device for determining the phase of fault activity, which addresses the problem in the prior art that does not consider the influence of fault activity on non-penetrating strata, thereby improving the accuracy of fault activity phase determination. The device includes:

[0011] The target fault location and penetration range determination module is used to determine the location of the target fault and the penetration range of the target fault based on the acquired seismic data.

[0012] The standard interface and preset key interface determination module is used to determine the standard interface of the strata that the fault penetrates and multiple preset key interfaces of the strata that it does not penetrate, based on the location and stratigraphic extent of the target fault. The standard interface is the strata within the penetrated stratigraphic extent that are affected by fault activity; the multiple preset key interfaces are the strata outside the penetrated stratigraphic extent.

[0013] The similarity coefficient calculation module is used to perform curve sampling on the standard interface and multiple preset standard interfaces to obtain the sampling point set of the standard interface and the sampling point set of multiple preset key interfaces respectively. Based on the sampling point set of the standard interface and the sampling point set of multiple preset key interfaces, the similarity coefficient between each preset key interface and the standard interface is calculated.

[0014] The fault activity impact determination module is used to compare the similarity coefficients between adjacent preset standard interfaces and determine whether the preset key interfaces are affected by fault activity based on the comparison results.

[0015] The fault activity period determination module is used to convert the geological time corresponding to the preset key interface affected by fault activity into time series data, and determine the fault activity period based on the time series data.

[0016] This invention also provides a computer device, including a memory, a processor, and a computer program stored in the memory and executable on the processor. When the processor executes the computer program, it implements the above-described method for determining the fault activity period.

[0017] This invention also provides a computer-readable storage medium storing a computer program that, when executed by a processor, implements the above-described method for determining fault activity periods.

[0018] This invention also provides a computer program product, which includes a computer program that, when executed by a processor, implements the above-described method for determining fault activity periods.

[0019] In this embodiment of the invention, the location of the target fault and the extent of the strata it penetrates are determined using acquired seismic data. Based on the location and extent of the target fault, a standard interface of the strata it penetrates and multiple preset key interfaces of the strata it does not penetrate are determined. The standard interface is the stratum within the penetrated stratum that is affected by fault activity; the multiple preset key interfaces are the strata outside the penetrated stratum. Curve sampling is performed on the standard interface and the multiple preset standard interfaces to obtain the sampling point sets of the standard interface and the multiple preset key interfaces, respectively. Based on the sampling point sets of the standard interface and the multiple preset key interfaces, a similarity coefficient between each preset key interface and the standard interface is calculated. The similarity coefficients between adjacent preset standard interfaces are compared, and the preset key interfaces are determined to be affected by fault activity based on the comparison results. The geological time corresponding to the preset key interfaces affected by fault activity is converted into time series data, and the fault activity period is determined based on the time series data. In the above process, the embodiments of the present invention consider the impact of fault activity on unpenetrated strata, convert the preset key interfaces and standard interfaces into curves for curve sampling, calculate the similarity coefficient between each preset key interface and the standard interface, and determine whether each unpenetrated preset key interface has been affected by fault activity. This solves the problem in the prior art that the impact of fault activity on unpenetrated strata is not considered. The fault activity period is determined based on the unpenetrated preset key interfaces affected by fault activity, thereby improving the accuracy of fault activity period determination. Attached Figure Description

[0020] To more clearly illustrate the technical solutions in the embodiments of the present invention or the prior art, the drawings used in the description of the embodiments or the prior art will be briefly introduced below. Obviously, the drawings described below are only some embodiments of the present invention. For those skilled in the art, other drawings can be obtained based on these drawings without creative effort. In the drawings:

[0021] Figure 1 This is a flowchart of the method for determining fault activity phases in an embodiment of the present invention;

[0022] Figure 2 This is a schematic diagram of the preset key interface and standard interface in an embodiment of the present invention;

[0023] Figure 3 This is a schematic diagram showing the standard interface and preset key interface curves after sampling in an embodiment of the present invention.

[0024] Figure 4 This is a flowchart of the optimal distance calculation in an embodiment of the present invention;

[0025] Figure 5 This is a flowchart illustrating the process of determining whether a formation interface is affected by fault activity in an embodiment of the present invention.

[0026] Figure 6 This is a schematic diagram of the fault activity period determination device in an embodiment of the present invention. Detailed Implementation

[0027] To make the objectives, technical solutions, and advantages of the embodiments of the present invention clearer, the embodiments of the present invention will be further described in detail below with reference to the accompanying drawings. Here, the illustrative embodiments of the present invention and their descriptions are used to explain the present invention, but are not intended to limit the present invention.

[0028] The acquisition, storage, use, and processing of data in this application all comply with the relevant provisions of national laws and regulations.

[0029] Figure 1 This is a flowchart of a method for determining fault activity phases in an embodiment of the present invention. The method includes:

[0030] Step 101: Based on the acquired seismic data, determine the location of the target fault and the extent of the strata it penetrates;

[0031] Step 102: Based on the location and stratigraphic extent of the target fault, determine the standard interface of the strata that the fault penetrates and multiple preset key interfaces of the strata that it does not penetrate. The standard interface is the strata within the penetrating stratigraphic extent that are affected by fault activity; the multiple preset key interfaces are the strata outside the penetrating stratigraphic extent.

[0032] Step 103: Perform curve sampling on the standard interface and multiple preset standard interfaces to obtain the sampling point set of the standard interface and the sampling point set of multiple preset key interfaces respectively. Based on the sampling point set of the standard interface and the sampling point set of multiple preset key interfaces, calculate the similarity coefficient between each preset key interface and the standard interface.

[0033] Step 104: Compare the similarity coefficients between adjacent preset standard interfaces, and determine the preset key interfaces affected by fault activity based on the comparison results.

[0034] Step 105: Convert the geological time corresponding to the preset key interface affected by fault activity into time series data, and determine the fault activity period based on the time series data.

[0035] Each step is explained in detail below.

[0036] In step 101, based on the acquired seismic data, the location of the target fault and the extent of the strata it penetrates are determined. In step 102, based on the location of the target fault and the extent of the strata it penetrates, the standard interface of the strata it penetrates and multiple preset key interfaces of the strata it does not penetrate are determined. The standard interface is the strata within the extent of the strata it penetrates that are affected by fault activity. The multiple preset key interfaces are the strata outside the extent of the strata it penetrates.

[0037] Figure 2 This is a schematic diagram of the preset key interface and standard interface in an embodiment of the present invention. In a specific embodiment, such as... Figure 2 As shown, a high-quality 2D seismic dataset from a certain region was selected to obtain the location, attitude, and extent of the target normal fault. Within the hanging wall of the fault and the strata penetrated by the normal fault, a clearly imaged stratigraphic interface with continuous phase axes was selected as the standard interface SL. Within the hanging wall of the fault and the strata not penetrated by the normal fault, several pre-defined key interfaces for the unpenetrated strata were determined, namely, interfaces A, B, C, and D overlying the hanging wall of the normal fault. By studying whether stratigraphic interfaces A, B, C, and D are affected by fault activity, the phase of fault activity can be determined. Figure 2 The red vertical lines in the middle indicate faults whose formation period needs to be determined.

[0038] In step 103, curve sampling is performed on the standard interface and multiple preset standard interfaces to obtain the sampling point set of the standard interface and the sampling point set of multiple preset key interfaces, respectively. Based on the sampling point set of the standard interface and the sampling point set of multiple preset key interfaces, the similarity coefficient between each preset key interface and the standard interface is calculated.

[0039] Figure 3 This is a schematic diagram of the standard interface and preset key interface curves after sampling in an embodiment of the present invention. In a specific embodiment, before sampling the standard interface and multiple preset standard interfaces, a coordinate transformation is first performed, and the standard interface curve l is selected. SL Middle point l SL0 Draw a vertical line from its location in the cross-section, and the vertical line intersects the stratigraphic interface curve l. A l B The intersection of ... is denoted as l. A0 l B0 ...; using point l in the standard interface curve SL0 A two-dimensional coordinate system is established with the origin as the origin, and the horizontal and vertical directions are used as the x and y axes, respectively. Axe and y coordinates are established based on fixed horizontal and vertical distances. In this two-dimensional coordinate system, samples are taken at appropriate intervals along the horizontal direction to determine the coordinates of the curve corresponding to the standard interface on the corresponding x-axis, forming a point set representing the curve. For curves corresponding to other preset key interfaces, a set of points is formed using l... A0 l B0Using , ... as the origin, establish two-dimensional coordinate systems identical to the standard interface curves. Sample at the same intervals as the standard interface curves and determine the coordinates of each point in the curve point set corresponding to the preset key interfaces to form a series of point sets representing the curves; thus obtaining the sampling point set of the standard interface and the sampling point sets of multiple preset key interfaces.

[0040] In one embodiment, based on a set of sampling points from a standard interface and a set of sampling points from multiple preset key interfaces, a similarity coefficient between each preset key interface and the standard interface is calculated, including:

[0041] The Pearson correlation algorithm is used to calculate the similarity coefficient between each preset key interface and the standard interface based on the sampling point set of the standard interface and the sampling point set of multiple preset key interfaces.

[0042] In one embodiment, the Pearson correlation algorithm is used to calculate the similarity coefficient between each preset key interface and the standard interface based on the sampling point set of the standard interface and the sampling point set of multiple preset key interfaces, including:

[0043] The similarity coefficient is calculated using the following formula:

[0044]

[0045] Where r(A,SL) is the Pearson correlation coefficient, A is the preset key interface, SL is the standard interface, and d DTW (A,SL) represents the optimal distance between the preset key interface and the standard interface calculated using the dynamic programming algorithm, and α is the weighting coefficient. By combining the Pearson correlation coefficient (global similarity) with the optimal distance (local morphological difference) of the dynamic programming algorithm, the deformation characteristics of the stratigraphic interface are comprehensively captured, avoiding the one-sidedness of a single index, which could lead to inaccurate similarity coefficient calculation.

[0046] Figure 4 This is a flowchart illustrating the optimal distance calculation in an embodiment of the present invention. In one embodiment, a dynamic programming algorithm is used to calculate the optimal distance between a preset key interface and a standard interface, including:

[0047] Step 401: Construct the distance matrix between the sampling point set of the preset key interface and the sampling point set of the standard interface;

[0048] Step 402: Using a dynamic programming algorithm, the optimal distance between the preset key interface and the standard interface is calculated based on the distance matrix between points.

[0049] In a specific embodiment, a dynamic programming algorithm is used to calculate the curves of each preset key interface. A l B ...and standard interface curve l SL The optimal distance d DTWGiven (A,SL), analyze the differences in curve shape, construct a distance matrix, calculate the distance between points on the two curves, and find the optimal path d. DTW (A,SL) is used to calculate the similarity coefficient.

[0050] In step 104, the similarity coefficients between adjacent preset standard interfaces are compared, and the preset key interfaces are determined to be affected by fault activity based on the comparison results.

[0051] Figure 5 This is a flowchart illustrating the process of determining whether a formation interface is affected by fault activity in an embodiment of the present invention. In one embodiment, the similarity coefficients between adjacent preset standard interfaces are compared, and the preset key interface is determined to be affected by fault activity based on the comparison results, including:

[0052] Step 501: Extract the similarity coefficients corresponding to adjacent preset standard interfaces, and calculate the relative rate of change between the similarity coefficients of adjacent preset standard interfaces.

[0053] Step 502: Compare the relative rate of change with a preset threshold. When the relative rate of change exceeds the preset threshold, determine that the preset key interface corresponding to the relative rate of change is affected by fault activity.

[0054] In a specific embodiment, similarity coefficients R1, R2, R3, and R4 corresponding to at least four adjacent preset standard interfaces are extracted from bottom to top, where R1 is the similarity coefficient of the bottommost interface and R4 is the similarity coefficient of the topmost interface; the relative rate of change between the similarity coefficients of adjacent preset standard interfaces is calculated.

[0055]

[0056] Where, ΔR 1-2 Let ΔR be the relative rate of change between the first preset standard interface similarity coefficient R1 and the second preset standard interface similarity coefficient R2. 2-3 ΔR represents the relative rate of change between the second preset standard interface similarity coefficient R2 and the third preset standard interface similarity coefficient R3. 3-4 Let R be the relative rate of change between the third preset standard interface similarity coefficient R3 and the fourth preset standard interface similarity coefficient R4. If ΔR 1-2 and ΔR 3-4 All are less than the preset threshold, and ΔR 2-3If the deformation exceeds a preset threshold, a deformation difference is determined to exist between the second and third preset standard interfaces. The third preset standard interface corresponding to the deformation difference represents strata unaffected by fault activity, while the second preset standard interface represents strata affected by fault activity. The geological time corresponding to the second preset standard interface represents the fault activity period. By analyzing the deformation differences at the interfaces of overlying strata not penetrated by the fault, the technical problem of traditional methods being unable to effectively determine the impact of fault activity on strata not penetrated by the fault is solved.

[0057] In step 105, the geological time corresponding to the preset key interface affected by fault activity is converted into time series data, and the fault activity period is determined based on the time series data.

[0058] In a specific embodiment, structural features corresponding to the identified fault activity periods are identified in the seismic profile, including angular unconformities, stratigraphic thickening, or local folding deformation, in order to confirm the reliability of the fault activity period determination results.

[0059] This invention also provides a device for determining the phase of fault activity, as described in the following embodiments. Since the principle behind this device is similar to that of the method for determining the phase of fault activity, the implementation of this device can be referred to the implementation of the method for determining the phase of fault activity, and repeated details will not be elaborated further.

[0060] Figure 6 This is a schematic diagram of a fault activity phase determination device according to an embodiment of the present invention. The device includes:

[0061] The target fault location and penetration range determination module 601 is used to determine the location of the target fault and the penetration range of the target fault based on the acquired seismic data.

[0062] The standard interface and preset key interface determination module 602 is used to determine the standard interface of the strata that the fault penetrates and multiple preset key interfaces of the strata that it does not penetrate, based on the location of the target fault and the stratigraphic range of the target fault. The standard interface is the strata that are affected by fault activity within the strata that penetrate the fault. The multiple preset key interfaces are the strata that are outside the strata that penetrate the fault.

[0063] The similarity coefficient calculation module 603 is used to perform curve sampling on the standard interface and multiple preset standard interfaces to obtain the sampling point set of the standard interface and the sampling point set of multiple preset key interfaces respectively. Based on the sampling point set of the standard interface and the sampling point set of multiple preset key interfaces, the similarity coefficient between each preset key interface and the standard interface is calculated.

[0064] The fault activity impact determination module 604 is used to compare the similarity coefficients between adjacent preset standard interfaces and determine the preset key interfaces affected by fault activity based on the comparison results.

[0065] The fault activity period determination module 605 is used to convert the geological time corresponding to the preset key interface affected by fault activity into time series data, and determine the fault activity period based on the time series data.

[0066] In one embodiment, the similarity coefficient calculation module 603 is specifically used for:

[0067] The Pearson correlation algorithm is used to calculate the similarity coefficient between each preset key interface and the standard interface based on the sampling point set of the standard interface and the sampling point set of multiple preset key interfaces.

[0068] In one embodiment, the similarity coefficient calculation module 603 is specifically used for:

[0069] The similarity coefficient is calculated using the following formula:

[0070]

[0071] Where r(A,SL) is the Pearson correlation coefficient, A is the preset key interface, SL is the standard interface, and d DTW (A,SL) represents the optimal distance between the preset key interface and the standard interface calculated using the dynamic programming algorithm, where α is the weight coefficient.

[0072] In one embodiment, an optimal distance calculation module is further included, specifically for:

[0073] Construct a distance matrix between the sampling point set of the preset key interface and the sampling point set of the standard interface;

[0074] The optimal distance between the preset key interface and the standard interface is calculated using a dynamic programming algorithm based on the distance matrix between points.

[0075] In one embodiment, the fault activity influence determination module 604 is specifically used for:

[0076] Extract the similarity coefficients of adjacent preset standard interfaces and calculate the relative rate of change between the similarity coefficients of adjacent preset standard interfaces;

[0077] The relative rate of change is compared with a preset threshold. When the relative rate of change exceeds the preset threshold, it is determined that the preset key interface corresponding to the relative rate of change is affected by fault activity.

[0078] This invention also provides a computer device, including a memory, a processor, and a computer program stored in the memory and executable on the processor. When the processor executes the computer program, it implements the above-described method for determining the fault activity period.

[0079] This invention also provides a computer-readable storage medium storing a computer program that, when executed by a processor, implements the above-described method for determining fault activity periods.

[0080] This invention also provides a computer program product, which includes a computer program that, when executed by a processor, implements the above-described method for determining fault activity periods.

[0081] In this embodiment of the invention, the location of the target fault and the extent of the strata it penetrates are determined using acquired seismic data. Based on the location and extent of the target fault, a standard interface of the strata it penetrates and multiple preset key interfaces of the strata it does not penetrate are determined. The standard interface is the stratum within the penetrated stratum that is affected by fault activity; the multiple preset key interfaces are the strata outside the penetrated stratum. Curve sampling is performed on the standard interface and the multiple preset standard interfaces to obtain the sampling point sets of the standard interface and the multiple preset key interfaces, respectively. Based on the sampling point sets of the standard interface and the multiple preset key interfaces, a similarity coefficient between each preset key interface and the standard interface is calculated. The similarity coefficients between adjacent preset standard interfaces are compared, and the preset key interfaces are determined to be affected by fault activity based on the comparison results. The geological time corresponding to the preset key interfaces affected by fault activity is converted into time series data, and the fault activity period is determined based on the time series data. In the above process, the embodiments of the present invention consider the impact of fault activity on unpenetrated strata, convert the preset key interfaces and standard interfaces into curves for curve sampling, calculate the similarity coefficient between each preset key interface and the standard interface, and determine whether each unpenetrated preset key interface has been affected by fault activity. This solves the problem in the prior art that the impact of fault activity on unpenetrated strata is not considered. The fault activity period is determined based on the unpenetrated preset key interfaces affected by fault activity, thereby improving the accuracy of fault activity period determination.

[0082] Those skilled in the art will understand that embodiments of the present invention can be provided as methods, systems, or computer program products. Therefore, the present invention can take the form of a completely hardware embodiment, a completely software embodiment, or an embodiment combining software and hardware aspects. Furthermore, the present invention can take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, etc.) containing computer-usable program code.

[0083] This invention is described with reference to flowchart illustrations and / or block diagrams of methods, apparatus (systems), and computer program products according to embodiments of the invention. It will be understood that each block of the flowchart illustrations and / or block diagrams, and combinations of blocks in the flowchart illustrations and / or block diagrams, can be implemented by computer program instructions. These computer program instructions can be provided to a processor of a general-purpose computer, special-purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, generate instructions for implementing the flowchart illustrations and / or block diagrams. Figure 1 One or more processes and / or boxes Figure 1 A device that provides the functions specified in one or more boxes.

[0084] These computer program instructions may also be stored in a computer-readable storage medium that can direct a computer or other programmable data processing device to function in a particular manner, such that the instructions stored in the computer-readable storage medium produce an article of manufacture including instruction means, which are implemented in a process Figure 1 One or more processes and / or boxes Figure 1 The function specified in one or more boxes.

[0085] These computer program instructions may also be loaded onto a computer or other programmable data processing equipment to cause a series of operational steps to be performed on the computer or other programmable equipment to produce a computer-implemented process, thereby providing instructions that execute on the computer or other programmable equipment for implementing the process. Figure 1 One or more processes and / or boxes Figure 1 The steps of the function specified in one or more boxes.

[0086] The specific embodiments described above further illustrate the purpose, technical solution, and beneficial effects of the present invention. It should be understood that the above descriptions are merely specific embodiments of the present invention and are not intended to limit the scope of protection of the present invention. Any modifications, equivalent substitutions, improvements, etc., made within the spirit and principles of the present invention should be included within the scope of protection of the present invention.

Claims

1. A method for determining the phases of fault activity, characterized in that, include: Based on the acquired seismic data, the location of the target fault and the extent of the strata it penetrates are determined. Based on the location and stratigraphic extent of the target fault, determine the standard interface of the strata that the fault penetrates and several preset key interfaces of the strata that it does not penetrate; The standard interface is the strata affected by fault activity within the range of the penetrating strata; Several preset key interfaces are for strata outside the scope of the penetrating strata; Curve sampling is performed on the standard interface and multiple preset standard interfaces to obtain the sampling point set of the standard interface and the sampling point set of multiple preset key interfaces respectively. Based on the sampling point set of the standard interface and the sampling point set of multiple preset key interfaces, the similarity coefficient between each preset key interface and the standard interface is calculated. The similarity coefficients between adjacent preset standard interfaces are compared, and the preset key interfaces are determined to be affected by fault activity based on the comparison results. The geological time corresponding to the preset key interface affected by fault activity is converted into time series data, and the fault activity period is determined based on the time series data.

2. The method as described in claim 1, characterized in that, Based on the sampling point set of the standard interface and the sampling point sets of multiple preset key interfaces, the similarity coefficient between each preset key interface and the standard interface is calculated, including: The Pearson correlation algorithm is used to calculate the similarity coefficient between each preset key interface and the standard interface based on the sampling point set of the standard interface and the sampling point set of multiple preset key interfaces.

3. The method as described in claim 2, characterized in that, The Pearson correlation algorithm is used to calculate the similarity coefficient between each preset key interface and the standard interface, based on the sampling point set of the standard interface and the sampling point set of multiple preset key interfaces. This includes: The similarity coefficient is calculated using the following formula: Where r(A,SL) is the Pearson correlation coefficient, A is the preset key interface, SL is the standard interface, and d DTW (A,SL) represents the optimal distance between the preset key interface and the standard interface calculated using the dynamic programming algorithm, where α is the weight coefficient.

4. The method as described in claim 3, characterized in that, The optimal distance between the preset key interface and the standard interface is calculated using a dynamic programming algorithm, including: Construct a distance matrix between the sampling point set of the preset key interface and the sampling point set of the standard interface; The optimal distance between the preset key interface and the standard interface is calculated using a dynamic programming algorithm based on the distance matrix between points.

5. The method as described in claim 1, characterized in that, The similarity coefficients between adjacent preset standard interfaces are compared, and the preset key interfaces are determined to be affected by fault activity based on the comparison results, including: Extract the similarity coefficients of adjacent preset standard interfaces and calculate the relative rate of change between the similarity coefficients of adjacent preset standard interfaces; The relative rate of change is compared with a preset threshold. When the relative rate of change exceeds the preset threshold, it is determined that the preset key interface corresponding to the relative rate of change is affected by fault activity.

6. A device for determining the phases of fault activity, characterized in that, include: The target fault location and penetration range determination module is used to determine the location of the target fault and the penetration range of the target fault based on the acquired seismic data. The standard interface and preset key interface determination module is used to determine the standard interface of the strata that the fault penetrates and multiple preset key interfaces of the strata that it does not penetrate, based on the location of the target fault and the stratigraphic extent of the target fault. The standard interface is the strata affected by fault activity within the stratigraphic range. Several preset key interfaces are strata outside the scope of the penetrating strata; The similarity coefficient calculation module is used to perform curve sampling on the standard interface and multiple preset standard interfaces to obtain the sampling point set of the standard interface and the sampling point set of multiple preset key interfaces respectively. Based on the sampling point set of the standard interface and the sampling point set of multiple preset key interfaces, the similarity coefficient between each preset key interface and the standard interface is calculated. The fault activity impact determination module is used to compare the similarity coefficients between adjacent preset standard interfaces and determine whether the preset key interfaces are affected by fault activity based on the comparison results. The fault activity period determination module is used to convert the geological time corresponding to the preset key interface affected by fault activity into time series data, and determine the fault activity period based on the time series data.

7. The apparatus as claimed in claim 6, characterized in that, The similarity coefficient calculation module is specifically used for: The Pearson correlation algorithm is used to calculate the similarity coefficient between each preset key interface and the standard interface based on the sampling point set of the standard interface and the sampling point set of multiple preset key interfaces.

8. The apparatus as claimed in claim 7, characterized in that, The similarity coefficient calculation module is specifically used for: The similarity coefficient is calculated using the following formula: Where r(A,SL) is the Pearson correlation coefficient, A is the preset key interface, SL is the standard interface, and d DTW (A,SL) represents the optimal distance between the preset key interface and the standard interface calculated using the dynamic programming algorithm, where α is the weight coefficient.

9. The apparatus as claimed in claim 8, characterized in that, It also includes an optimal distance calculation module, specifically used for: Construct a distance matrix between the sampling point set of the preset key interface and the sampling point set of the standard interface; The optimal distance between the preset key interface and the standard interface is calculated using a dynamic programming algorithm based on the distance matrix between points.

10. The apparatus as claimed in claim 6, characterized in that, The fault activity impact determination module is specifically used for: Extract the similarity coefficients of adjacent preset standard interfaces and calculate the relative rate of change between the similarity coefficients of adjacent preset standard interfaces; The relative rate of change is compared with a preset threshold. When the relative rate of change exceeds the preset threshold, it is determined that the preset key interface corresponding to the relative rate of change is affected by fault activity.

11. A computer device, comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, characterized in that, When the processor executes the computer program, it implements the method of any one of claims 1 to 5.

12. A computer-readable storage medium, characterized in that, The computer-readable storage medium stores a computer program that, when executed by a processor, implements the method of any one of claims 1 to 5.

13. A computer program product, characterized in that, The computer program product includes a computer program that, when executed by a processor, implements the method of any one of claims 1 to 5.

Images