A method for evaluating the activity of a carbonate rock depositional slope syn-sedimentary fault by three-dimensional seismic

Abstract

This invention discloses a three-dimensional seismic assessment method for the activity of syn-sedimentary faults on gentle slopes in carbonate rocks, specifically relating to the field of oil and gas geological exploration. Its key points are: interpreting the fault plane of the syn-sedimentary fault and interpreting the stratigraphic plane of the fault-related strata; selecting strata for assessing fault activity from the interpretation results and constructing a sequence stratigraphic framework using these strata; calculating the sedimentary time of the fault activity period using the time corresponding to the sequence stratigraphic boundaries of the sequence stratigraphic framework; selecting multiple seismic profiles from the fault-stratigraphic plane interpretation results and obtaining the fault plane trajectory and stratigraphic plane trajectory lines passing through each seismic profile; calculating the average fault activity rate of each seismic profile at different periods based on the sedimentary time of the fault activity period; and assessing the activity of the syn-sedimentary fault based on the average fault activity rate of each seismic profile at different periods to obtain the seismic assessment results.

Description

Technical Field

[0001] This invention relates to the field of petroleum and natural gas geological exploration, specifically to a three-dimensional seismic assessment method for the activity of co-sedimentary faults on gentle slopes of carbonate rocks. Background Technology

[0002] Co-sedimentary faults play a controlling role in the paleogeographic patterns and sedimentary facies distribution of tectonically active areas. The study of co-sedimentary fault activity rates is of great value in regional tectonic evolution, prediction of favorable sedimentary facies distribution, and characterization of lithological traps, and is a crucial step in the analysis of regional petroleum geological conditions.

[0003] Currently, the main methods for quantitatively studying the activity rate of syn-sedimentary faults include the growth index method, paleodisplacement method, fault activity rate method, fault displacement-length relationship analysis method, and paleoslip method (Lei Baohua. A review of the main methods for quantitatively studying the activity intensity of growth faults [J]. Advances in Earth Science, 2012, 27(9):947-956). These methods can reflect the activity intensity and changes of syn-sedimentary faults from different perspectives, but in general, they have problems such as a lack of effective constraints on fault activity time, complex paleotectonic reconstruction processes, and difficulty in applying them to studies related to gentle slopes of carbonate rocks. Compared with traditional two-dimensional seismic data, these problems are more prominent when using three-dimensional seismic data to study syn-sedimentary faults.

[0004] For syn-sedimentary carbonate rocks, due to the flat paleotopography and similar platform subsidence rates, the carbonate deposition rates in different regions are comparable when unaffected by syn-sedimentary fault activity. Therefore, the apparent fault distance between the two sides of a fault during tectonic activity and its corresponding time span can reflect the fault's activity rate during that period. To improve the accuracy of quantitative assessment of the activity of syn-sedimentary carbonate rocks, this invention optimizes the calculation method for the apparent fault distance and its corresponding time span during fault activity periods, and establishes a three-dimensional seismic assessment method for the activity of syn-sedimentary carbonate rocks. Summary of the Invention

[0005] The technical problem this invention aims to solve is that existing technical solutions lack effective methods for constraining fault activity time, have complex paleotectonic reconstruction processes, and are difficult to apply to the study of tectonic-sedimentary evolution of carbonate gentle slopes. The goal is to provide a three-dimensional seismic assessment method for the activity of syn-sedimentary faults on carbonate gentle slopes. By selecting strata without obvious sedimentary discontinuities and establishing a sequence stratigraphic framework to constrain the calculation of fault activity time, the accuracy of quantitative assessment of syn-sedimentary fault activity is improved. Furthermore, by applying the built-in seismic profile paleotectonic reconstruction module of the Petrel geological interpretation software to reconstruct paleotectonic patterns during fault activity, the accuracy of calculating apparent fault displacement during fault activity is improved, and the paleotectonic reconstruction process is simplified.

[0006] This invention is achieved through the following technical solution:

[0007] A three-dimensional seismic assessment method for the activity of syn-sedimentary faults on gentle slopes of carbonate rocks, the method comprising:

[0008] In the 3D seismic data volume, fault plane interpretation is performed on syn-depositional faults and stratigraphic plane interpretation is performed on fault-related strata to obtain fault plane-stratigraphic plane interpretation results.

[0009] Based on the pre-constructed evaluation selection rules, strata for evaluating fault activity are selected from the fault-stratigraphic interpretation results, and a sequence stratigraphic framework is constructed using the selected strata; and the sedimentary time of the strata during the fault activity period is calculated using the time corresponding to the sequence stratigraphic boundaries of the sequence stratigraphic framework.

[0010] Based on the pre-constructed seismic profile selection rules, multiple seismic profiles are selected from the fault-stratum layer interpretation results, and the fault trajectory line and stratum layer trajectory line passing through each seismic profile are obtained.

[0011] By using the fault plane trajectory lines and stratigraphic plane trajectory lines of each seismic profile, combined with the stratigraphic deposition time during the fault activity period, the average fault activity rate of each seismic profile at different periods is calculated.

[0012] The activity of syn-depositional faults is assessed based on the average fault activity rate at different times for each seismic profile, resulting in seismic assessment results.

[0013] Furthermore, fault plane interpretation is performed on syn-sedimentary faults and stratigraphic plane interpretation is performed on fault-related strata within the 3D seismic data volume, yielding fault plane-stratigraphic plane interpretation results, specifically:

[0014] In the three-dimensional seismic data volume, the pre-set seismic survey line grid is used to perform fault layer interpretation on the same sedimentary fault, and the fault layer interpretation results are obtained.

[0015] Based on the fault-related strata interpretation results, the fault-related strata are obtained, and the stratigraphic plane interpretation of the fault-related strata is performed using a preset seismic survey grid, resulting in fault-stratigraphic plane interpretation results.

[0016] Furthermore, based on pre-constructed evaluation selection rules, strata for assessing fault activity are selected from the fault-stratum interpretation results, and a sequence stratigraphic framework is constructed using the selected strata. The depositional time of the strata during fault activity is calculated using the time corresponding to the sequence stratigraphic boundaries of the sequence stratigraphic framework. Specifically:

[0017] In the fault-stratigraphic interpretation results, strata that have been continuously deposited during geological history are selected to assess fault activity.

[0018] Conduct comparative studies of seismic stratigraphy and well logging stratigraphy on the selected lithostratigraphic units to construct a sequence stratigraphic framework;

[0019] By using the time corresponding to the sequence stratigraphic boundaries of the sequence stratigraphic framework, the stratigraphic deposition time during the fault activity period can be calculated.

[0020] Furthermore, based on the pre-constructed seismic profile selection rules, multiple seismic profiles are selected from the fault-stratum plane interpretation results, and the fault plane trajectory line and stratum plane trajectory line passing through each seismic profile are obtained, specifically:

[0021] Multiple seismic profiles are selected from the fault-stratum interpretation results at preset horizontal intervals, and the fault trajectory lines and stratum trajectory lines passing through each seismic profile are obtained.

[0022] Furthermore, using the fault plane trajectory lines and stratigraphic plane trajectory lines of each seismic profile, combined with the stratigraphic deposition time during the fault activity period, the average fault activity rate of each seismic profile at different periods was calculated, specifically:

[0023] Using paleotectonic reconstruction software, fault plane trajectory lines and stratigraphic plane trajectory lines of each seismic profile are used to perform paleotectonic reconstruction interpretation of the strata at the end of the stratigraphic deposition period, and the horizontal and vertical apparent fault displacements after paleotectonic reconstruction are read from the reconstruction interpretation results.

[0024] By using the horizontal and vertical apparent fault displacements reconstructed from paleotectonic structures, as well as the stratigraphic deposition time during the fault activity period, the average fault activity rate of each seismic profile at different periods was calculated.

[0025] Furthermore, using the horizontal and vertical apparent fault displacements reconstructed from paleotectonic data, as well as the stratigraphic deposition time during fault activity, the average fault activity rate for each seismic profile at different periods was calculated, specifically: V = / T, where V represents the average fault activity rate. Indicates the horizontal apparent displacement of the fault. The vertical apparent fault displacement is represented by T, which represents the sedimentary time of the strata during the period of fault activity.

[0026] This invention also provides a three-dimensional seismic assessment system for the activity of syn-sedimentary faults in carbonate rocks with gentle slopes. This system is used in any of the three-dimensional seismic assessment methods for the activity of syn-sedimentary faults in carbonate rocks with gentle slopes described above. The system includes:

[0027] The layer interpretation module is used to interpret fault planes of syn-sedimentary faults and to interpret stratigraphic layers of fault-related strata in 3D seismic data volumes, so as to obtain fault plane-stratigraphic layer interpretation results.

[0028] The deposition time calculation module is used to select strata for assessing fault activity from the fault-stratum layer interpretation results based on pre-constructed assessment selection rules, and to construct a sequence stratigraphic framework using the selected strata; and to calculate the deposition time of each set of strata during the fault activity period using the time corresponding to the sequence stratigraphic boundaries of the sequence stratigraphic framework.

[0029] The trajectory acquisition module is used to select multiple seismic profiles from the fault-stratum layer interpretation results based on pre-built seismic profile selection rules, and to acquire the fault trajectory line and stratum layer trajectory line passing through each seismic profile.

[0030] The activity rate calculation module is used to calculate the average fault activity rate of each seismic profile at different periods by using the fault plane trajectory line and the stratigraphic plane trajectory line of each seismic profile, combined with the deposition time of each set of strata during the fault activity period.

[0031] The seismic assessment module is used to assess the activity of syn-depositional faults based on the average fault activity rate at different times for each seismic profile, and to obtain seismic assessment results.

[0032] The present invention also provides a computer device, including a system memory and a processor, wherein the system memory stores a computer program, and the processor executes the computer program to implement the steps of any of the methods described above.

[0033] The present invention also provides a computer-readable storage medium having a computer program stored thereon, wherein the computer program, when executed by a processor, implements the steps of any of the methods described above.

[0034] The present invention also provides a computer program product containing instructions that, when executed by a cluster of computer devices, cause the cluster of computer devices to perform the method described in any of the preceding claims.

[0035] Compared with the prior art, the present invention has the following advantages and beneficial effects:

[0036] In this invention, by selecting strata without obvious sedimentary discontinuities and establishing a sequence stratigraphic framework to constrain the calculation of fault activity time, the accuracy of quantitative assessment of syn-sedimentary fault activity is improved; by applying the seismic profile paleotectonic restoration module built into the Petrel geological interpretation software to restore the paleotectonic style during fault activity, the calculation accuracy of apparent fault displacement during fault activity is improved, and the paleotectonic restoration process is simplified. Attached Figure Description

[0037] To more clearly illustrate the technical solutions of the exemplary embodiments of the present invention, the accompanying drawings used in the embodiments will be briefly described below. It should be understood that the following drawings only show some embodiments of the present invention and should not be considered as a limitation of the scope. For those skilled in the art, other related drawings can be obtained based on these drawings without creative effort. In the drawings:

[0038] Figure 1 This is a schematic diagram of the process flow for a three-dimensional seismic assessment method for the activity of co-sedimentary faults on gentle slopes of carbonate rocks in this embodiment;

[0039] Figure 2 This is a three-dimensional seismic interpretation diagram of the top and bottom surfaces of the syn-depositional fault F and the I, J, and K strata provided in this embodiment;

[0040] Figure 3 This is a sequence stratigraphic framework of strata I, J, and K provided in this embodiment, and its correspondence with the regional sea-level change curve.

[0041] Figure 4 This is a seismic profile 1 of the I, J, and K strata at the end of their depositional period, as provided in this embodiment and after paleotectonic reconstruction interpretation by Petrel geological software. The reconstructed profile is marked with the two-way travel time A and the horizontal apparent displacement X corresponding to the vertical apparent displacement.

[0042] Figure 5 This is a line graph showing the average activity rate of the syn-depositional fault F during the depositional periods of strata I, J, and K.

[0043] Figure 6 This is a schematic diagram of the module connections of a three-dimensional seismic assessment system for the activity of carbonate rock slopes and sedimentary faults in this embodiment;

[0044] Figure 7 This is a schematic diagram of the structure of a computer device in this embodiment. Detailed Implementation

[0045] The exemplary embodiments of this disclosure are described below with reference to the accompanying drawings, including various details of the embodiments to aid understanding, and should be considered merely exemplary. Therefore, those skilled in the art will recognize that various changes and modifications can be made to the embodiments described herein without departing from the scope of this disclosure. Similarly, for clarity and brevity, descriptions of well-known functions and structures are omitted in the following description.

[0046] In this disclosure, unless otherwise stated, the use of terms such as "first," "second," etc., to describe various elements is not intended to limit the positional, temporal, or importance relationships of these elements; such terms are merely used to distinguish one element from another. In some examples, the first element and the second element may refer to the same instance of that element, while in other cases, based on the context, they may refer to different instances.

[0047] The terminology used in the description of the various examples in this disclosure is for the purpose of describing particular examples only and is not intended to be limiting. Unless the context explicitly indicates otherwise, an element may be one or more unless the number of elements is specifically limited. Furthermore, the term "and / or" as used in this disclosure covers any one of the listed items and all possible combinations thereof.

[0048] Example 1

[0049] See Figure 1 , Figure 1 A schematic flowchart of a three-dimensional seismic assessment method for the activity of syn-sedimentary faults on gentle slopes of carbonate rocks is shown, wherein the method includes:

[0050] S1: In the 3D seismic data volume, fault plane interpretation is performed on the syn-depositional faults and stratigraphic plane interpretation is performed on the fault-related strata to obtain the fault plane-stratigraphic plane interpretation results.

[0051] Specifically, in this embodiment, the fault-layer interpretation of the syn-depositional fault is performed using a preset seismic survey grid in the three-dimensional seismic data volume to obtain the fault-layer interpretation result; based on the fault-layer interpretation result, the fault-related strata are obtained, and the stratigraphic plane interpretation of the fault-related strata is performed using the preset seismic survey grid to obtain the fault plane-stratigraphic plane interpretation result.

[0052] It should be noted that in this embodiment, the fault plane of the same sedimentary fault F is interpreted in the 3D seismic data volume using a 10*10 seismic survey grid density (the spacing between adjacent survey lines is 25m, i.e., the sampling interval is 250m). The interpretation results show that fault F cuts through three strata: I, J, and K, meaning that I, J, and K are fault-related strata. Therefore, the same seismic survey grid is used to interpret the top and bottom surfaces of the three strata: I, J, and K. Since the three strata are continuously developed, the bottom surface of stratum I is the top surface of stratum J, and the bottom surface of stratum J is the top surface of stratum K. The final interpretation results of the fault plane-stratum plane are as follows: Figure 2 As shown.

[0053] S2: Based on the pre-constructed evaluation selection rules, select strata for evaluating fault activity from the fault-stratum interpretation results, and construct a sequence stratigraphic framework using the selected strata; and calculate the stratigraphic deposition time during the fault activity period using the time corresponding to the sequence stratigraphic boundaries of the sequence stratigraphic framework.

[0054] Specifically, in this embodiment, strata that have been continuously deposited during geological history are selected from the fault-stratum interpretation results to assess fault activity; seismic stratigraphy and well logging stratigraphy are compared between the selected strata and the lithostratigraphic units to construct a sequence stratigraphic framework; and the deposition time of the strata during the fault activity period is calculated using the time corresponding to the sequence stratigraphic boundaries of the sequence stratigraphic framework.

[0055] It should be noted that, in this embodiment, sedimentary analysis of the three strata (I, J, and K) shows that none of them have obvious sedimentary discontinuities and can be considered as continuously deposited during geological history. Therefore, they can all be used to assess the activity rate of the syn-sedimentary fault F. Simultaneously, using 3D seismic and drilling data, a comparative seismic stratigraphy and well logging study was conducted on the Upper Cretaceous lithostratigraphic units containing the three strata (I, J, and K). The resulting sequence stratigraphic framework is as follows: Figure 3 As shown, from bottom to top, strata K, J, and I correspond to the second and third transgressive systems tracts and the first regressive systems tract of the Upper Cretaceous, respectively; subsequently, referring to historical sea-level change data of the Cretaceous region, a correspondence between the sequence stratigraphic framework and the regional sea-level change curve was established. Figure 3 Based on the ages of each sequence boundary provided in historical data, the geological time T1 corresponding to the bottom of stratum K is 88 Ma, and the geological time T2 corresponding to the top of stratum K is 85 Ma; the geological time T1 corresponding to the bottom of stratum J is 85 Ma, and the geological time T2 corresponding to the top of stratum J is 81 Ma; the geological time T1 corresponding to the bottom of stratum I is 81 Ma, and the geological time T2 corresponding to the top of stratum I is 75 Ma. Using the formula for calculating the sedimentary time T: T = T2 - T1, we can obtain that the sedimentary time of stratum K is 3 Ma, the sedimentary time of stratum J is 4 Ma, and the sedimentary time of stratum I is 6 Ma.

[0056] S3: Based on the pre-constructed seismic profile selection rules, multiple seismic profiles are selected from the fault-stratum plane interpretation results, and the fault plane trajectory line and stratum plane trajectory line passing through each seismic profile are obtained.

[0057] Specifically, in this embodiment, multiple seismic profiles are selected from the fault-stratum interpretation results at preset horizontal intervals, and the fault trajectory line and stratum trajectory line passing through each seismic profile are obtained.

[0058] It should be noted that in this embodiment, the trace spacing of the three-dimensional seismic data volume of the study area is 25m, and the horizontal spacing of the selected seismic profiles is set to 1km, that is, starting from the profile cut through by the fault, one seismic profile is selected every 40 traces; and the trajectory lines of the same sedimentary fault F and the I, J, K strata that have picked up these seismic profiles.

[0059] S4: Using the fault plane trajectory line and stratigraphic plane trajectory line of each seismic profile, combined with the stratigraphic deposition time during the fault activity period, the average fault activity rate of each seismic profile at different periods is calculated.

[0060] Specifically, in this embodiment, paleotectonic reconstruction software is used to interpret the strata at the end of the stratigraphic deposition period by utilizing the fault plane trajectory lines and stratigraphic plane trajectory lines of each seismic profile. The horizontal and vertical apparent fault displacements after paleotectonic reconstruction are read from the reconstruction results. Using the horizontal and vertical apparent fault displacements after paleotectonic reconstruction and the stratigraphic deposition time during the fault activity period, the average fault activity rate of each seismic profile at different periods is calculated.

[0061] It should be noted that in this embodiment, the built-in paleotectonic reconstruction and interpretation module of the Petrel geological software is used to perform paleotectonic reconstruction and interpretation of a certain stratum at the end of its depositional period. The horizontal apparent fault displacement X and vertical apparent fault displacement Y are read from the reconstruction and interpretation results. Then, using the stratum depositional time T, the horizontal apparent fault displacement X, and the vertical apparent fault displacement Y, the average fault activity rate V across the seismic profile during the depositional period of this stratum is calculated, specifically: V = / T; then change to different seismic profiles and fault activity periods, repeat the above steps until the average fault activity rate of all selected strata and seismic profiles is calculated, so as to draw a line graph of the average activity rate of syn-sedimentary faults in different periods.

[0062] For example, taking seismic profile 1, which is cut through by the syn-depositional fault F, as an example, the trajectory lines of the syn-depositional fault F and the bedding planes of strata I, J, and K passing through the profile are input into the paleotectonic reconstruction and interpretation module built into the Petrel geological software. Paleotectonic reconstruction of the I, J, and K strata at the end of their depositional period is performed, using the top surfaces of strata I, J, and K as reference surfaces. The two-way travel time A and horizontal apparent fault distance X corresponding to the vertical apparent fault distance are measured on the reconstructed profile. Figure 4As shown in the figure, the two-way travel time A corresponding to the vertical apparent displacement at the end of the depositional period of stratum I is 9.0 ms, and the horizontal apparent displacement X is 14.0 m; the two-way travel time A corresponding to the vertical apparent displacement at the end of the depositional period of stratum J is 8.6 ms, and the horizontal apparent displacement X is 10.8 m; the two-way travel time A corresponding to the vertical apparent displacement at the end of the depositional period of stratum K is 10.0 ms, and the horizontal apparent displacement X is 11.8 m. Based on the analysis of the synthetic seismic records from 3D seismic data and wells, the average sound velocities S in strata I, J, and K are 3.9 m / ms, 4.3 m / ms, and 4.0 m / ms, respectively. Using the formula for calculating the vertical apparent displacement Y: Y = (A × S) / 2, where A is the two-way travel time corresponding to the vertical apparent displacement in the seismic profile, and S is the average sound velocity in the stratum, the calculated vertical apparent displacements Y for strata I, J, and K are 17.5 m, 18.5 m, and 20.0 m, respectively.

[0063] Based on the horizontal apparent displacement X, vertical apparent displacement Y, and sedimentary time T of strata I, J, and K, the formula V = / T, the average fault activity rates V of the I, J, and K strata during the depositional period through seismic profile 1 were obtained as 3.7 m / Ma, 5.4 m / Ma, and 7.7 m / Ma, respectively;

[0064] Then, horizontal intervals were selected for the seismic profiles at 1km intervals, resulting in a total of 9 seismic profiles passing through the same sedimentary fault F. After calculating the average activity rate V of the same sedimentary fault F during the depositional periods of strata I, J, and K in seismic profile 1, the seismic profiles were switched, and the same method as in step S103 was used to calculate the average activity rate of the same sedimentary fault F during the depositional periods of strata I, J, and K in seismic profiles 2-9. The statistics of the calculation results for all seismic profiles are shown in the table below.

[0065]

[0066] S5: The activity of syn-depositional faults is assessed based on the average fault activity rate at different times for each seismic profile, resulting in seismic assessment results.

[0067] For example, based on the calculation results of all seismic profiles in the table above, with the nine seismic profiles crossing faults as the horizontal axis and the activity rate of the syn-depositional fault F as the vertical axis, line graphs of the average activity rate of the syn-depositional fault during the depositional periods of strata I, J, and K are plotted, as follows. Figure 5As shown, during the deposition of strata K in geological history, fault F, which was deposited with the same sedimentary layer, exhibited the strongest relative activity, with a maximum average activity rate exceeding 20 m / Ma. Faults from the depositional periods of strata I and J showed relatively weaker activity, with a maximum average activity rate of only about 10 m / Ma. Furthermore, on the plane, the average activity rate at the intersection of seismic profile 5-8 and the fault is relatively higher, decreasing towards both sides of the fault. Therefore, it can be concluded that the main period of activity for the syn-depositional fault F was during the deposition of strata K, and the fault developed by extending outwards from the intersection of seismic profile 5-8 and the fault towards both flanks.

[0068] Specifically, in this embodiment, by selecting strata without obvious sedimentary discontinuities and establishing a sequence stratigraphic framework to constrain the calculation of fault activity time, the accuracy of quantitative assessment of syn-sedimentary fault activity is improved; by applying the seismic profile paleotectonic restoration module built into the Petrel geological interpretation software to restore the paleotectonic style during fault activity, the calculation accuracy of apparent fault displacement during fault activity is improved, and the paleotectonic restoration process is simplified.

[0069] Example 2

[0070] See Figure 6 As shown, the present invention also provides a three-dimensional seismic assessment system for the activity of carbonate rock syn-sedimentary faults on gentle slopes. This system is used in any of the three-dimensional seismic assessment methods for the activity of carbonate rock syn-sedimentary faults on gentle slopes described above. The system includes:

[0071] The layer interpretation module 100 is used to interpret fault planes of syn-depositional faults and to interpret stratigraphic layers of fault-related strata in the three-dimensional seismic data volume, so as to obtain fault plane-stratigraphic layer interpretation results.

[0072] The deposition time calculation module 200 is used to select strata for assessing fault activity from the fault-stratum layer interpretation results based on pre-constructed assessment selection rules, and to construct a sequence stratigraphic framework using the selected strata; and to calculate the deposition time of each set of strata during the fault activity period using the time corresponding to the sequence stratigraphic boundaries of the sequence stratigraphic framework.

[0073] The trajectory line acquisition module 300 is used to select multiple seismic profiles in the fault-stratum layer interpretation results based on pre-built seismic profile selection rules, and to acquire the fault trajectory line and stratum layer trajectory line passing through each seismic profile.

[0074] The activity rate calculation module 400 is used to calculate the average fault activity rate of each seismic profile at different periods by using the fault plane trajectory line and the stratigraphic plane trajectory line of each seismic profile, combined with the deposition time of each set of strata during the fault activity period.

[0075] The seismic assessment module 500 is used to assess the activity of syn-depositional faults based on the average fault activity rate at different times for each seismic profile, and to obtain seismic assessment results.

[0076] It should be noted that the modules in the system of Embodiment 2 correspond to the steps in the method of Embodiment 1. The steps in the method of Embodiment 1 have been described in detail in Embodiment 1, and the module content in the system will not be described in detail in this Embodiment 2.

[0077] Example 3

[0078] See Figure 7 As shown, this embodiment also provides a computer device, including a system memory 1005 and a processor 1001. The system memory 1005 stores a computer program, and the processor 1001 executes the computer program to implement the steps of any of the methods described above.

[0079] It should be noted that the processor 1001 is used to execute the steps in the above method embodiments according to the instructions in the program code. Alternatively, when the processor 1001 executes the computer program, it implements the functions of each module / unit in the above system / device embodiments.

[0080] Specifically, in this embodiment, the computer program can be divided into one or more modules / units. One or more modules / units are stored in the system memory 1005 and executed by the processor 1001 to complete this application. One or more modules / units can be a series of computer program instruction segments capable of performing specific functions, which describe the execution process of the computer program in the terminal device.

[0081] The terminal device can be a desktop computer, laptop, handheld computer, or cloud server, etc. The terminal device may include, but is not limited to, a processor 1001 and a system memory 1005. Those skilled in the art will understand that this does not constitute a limitation on the terminal device, which may include more or fewer components than shown, or a combination of certain components, or different components. For example, the terminal device may also include an input / output device 1003, a network access device 1002, a bus 1006, etc.

[0082] The processor 1001 can be a Central Processing Unit (CPU), or other general-purpose processors, digital signal processors (DSPs), application-specific integrated circuits (ASICs), field-programmable gate arrays (FPGAs), or other programmable logic devices, discrete gate or transistor logic devices, discrete hardware components, etc. A general-purpose processor can be a microprocessor or any conventional processor.

[0083] System memory 1005 can be an internal storage unit of the terminal device, such as a hard drive or RAM. System memory 1005 can also be a storage device 1004 of the terminal device, such as an external hard drive, SmartMedia Card (SMC), Secure Digital (SD) card, or FlashCard. Furthermore, system memory 1005 can include both internal storage units and storage device 1004. System memory 1005 is used to store computer programs and other programs and data required by the terminal device. System memory 1005 can also be used to temporarily store data that has been output or will be output.

[0084] Those skilled in the art will clearly understand that, for the sake of convenience and brevity, the specific working processes of the systems and units described above can be referred to the corresponding processes in the foregoing method embodiments, and will not be repeated here.

[0085] Example 4

[0086] This embodiment provides a computer-readable storage medium having a computer program stored thereon, which, when executed by a processor, implements the steps of any of the methods described above.

[0087] The computer-readable storage medium may be, for example, but not limited to, an electrical, magnetic, optical, electromagnetic, infrared, or semiconductor system, device, or any combination thereof. More specific examples of computer-readable storage media (a non-exhaustive list) include: an electrical connection having one or more wires, a portable computer disk, a hard disk, random access memory (RAM), read-only memory (ROM), erasable programmable read-only memory (EPROM), registers, hard disks, optical fibers, compact disc read-only memory (CD-ROM), optical storage devices, magnetic storage devices, or any suitable combination thereof, or any other form of computer-readable storage medium in the art.

[0088] An exemplary storage medium is coupled to a processor, enabling the processor to read information from and write information to the storage medium. Of course, the storage medium can also be a component of the processor. The processor and storage medium can reside within an application-specific integrated circuit (ASIC). In embodiments of the invention, the computer-readable storage medium can be any tangible medium containing or storing a program that can be used by or in conjunction with an instruction execution system, system, or device.

[0089] Example 5

[0090] This embodiment also provides a computer program product containing instructions that, when executed by a cluster of computer devices, cause the cluster of computer devices to perform the method described in Embodiment 1.

[0091] 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 description is only a specific embodiment of the present invention and is 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 three-dimensional seismic assessment method for the activity of co-sedimentary faults on gentle slopes of carbonate rocks, characterized in that, The methods include: In the 3D seismic data volume, fault plane interpretation is performed on syn-depositional faults and stratigraphic plane interpretation is performed on fault-related strata to obtain fault plane-stratigraphic plane interpretation results. Based on the pre-constructed evaluation selection rules, strata for evaluating fault activity are selected from the fault-stratum interpretation results, and the selected strata are used to construct a sequence stratigraphic framework. And by using the time corresponding to the sequence stratigraphic boundaries of the sequence stratigraphic framework, the stratigraphic deposition time during the fault activity period can be calculated. Based on the pre-constructed seismic profile selection rules, multiple seismic profiles are selected from the fault-stratum layer interpretation results, and the fault trajectory line and stratum layer trajectory line passing through each seismic profile are obtained. By using the fault plane trajectory lines and stratigraphic plane trajectory lines of each seismic profile, combined with the stratigraphic deposition time during the fault activity period, the average fault activity rate of each seismic profile at different periods is calculated. The activity of syn-depositional faults is assessed based on the average fault activity rate at different times for each seismic profile, resulting in seismic assessment results.

2. The three-dimensional seismic assessment method for the activity of syn-sedimentary faults on gentle slopes of carbonate rocks according to claim 1, characterized in that, In the 3D seismic data volume, fault plane interpretation is performed on syn-sedimentary faults, and stratigraphic plane interpretation is performed on fault-related strata, yielding fault plane-stratigraphic plane interpretation results, specifically: In the three-dimensional seismic data volume, the pre-set seismic survey line grid is used to perform fault layer interpretation on the same sedimentary fault, and the fault layer interpretation results are obtained. Based on the fault-related strata interpretation results, the fault-related strata are obtained, and the stratigraphic plane interpretation of the fault-related strata is performed using a preset seismic survey grid, resulting in fault-stratigraphic plane interpretation results.

3. The three-dimensional seismic assessment method for the activity of syn-sedimentary faults on gentle slopes of carbonate rocks according to claim 1, characterized in that, Based on pre-constructed evaluation selection rules, strata for assessing fault activity are selected from the fault-stratigraphic interpretation results, and a sequence stratigraphic framework is constructed using the selected strata. Furthermore, the depositional time of the strata during fault activity is calculated using the time corresponding to the sequence stratigraphic boundaries of the sequence stratigraphic framework. Specifically: In the fault-stratigraphic interpretation results, strata that have been continuously deposited during geological history are selected to assess fault activity. Conduct comparative studies of seismic stratigraphy and well logging stratigraphy on the selected lithostratigraphic units to construct a sequence stratigraphic framework; By using the time corresponding to the sequence stratigraphic boundaries of the sequence stratigraphic framework, the stratigraphic deposition time during the fault activity period can be calculated.

4. The three-dimensional seismic assessment method for the activity of syn-sedimentary faults on gentle slopes of carbonate rocks according to claim 1, characterized in that, Based on pre-constructed seismic profile selection rules, multiple seismic profiles are selected from the fault-stratum plane interpretation results, and the fault plane trajectory line and stratum plane trajectory line passing through each seismic profile are obtained, specifically: Multiple seismic profiles are selected from the fault-stratum interpretation results at preset horizontal intervals, and the fault trajectory lines and stratum trajectory lines passing through each seismic profile are obtained.

5. A three-dimensional seismic assessment method for the activity of syn-sedimentary faults on gentle slopes of carbonate rocks according to claim 1, characterized in that, Using the fault plane trajectory lines and stratigraphic plane trajectory lines of each seismic profile, combined with the stratigraphic deposition time during the fault activity period, the average fault activity rate of each seismic profile at different periods was calculated, specifically: Using paleotectonic reconstruction software, fault plane trajectory lines and stratigraphic plane trajectory lines of each seismic profile are used to perform paleotectonic reconstruction interpretation of the strata at the end of the stratigraphic deposition period, and the horizontal and vertical apparent fault displacements after paleotectonic reconstruction are read from the reconstruction interpretation results. By using the horizontal and vertical apparent fault displacements reconstructed from paleotectonic structures, as well as the stratigraphic deposition time during the fault activity period, the average fault activity rate of each seismic profile at different periods was calculated.

6. A three-dimensional seismic assessment method for the activity of syn-sedimentary faults on gentle slopes of carbonate rocks according to claim 5, characterized in that, Using the horizontal and vertical apparent fault displacements reconstructed from paleotectonic data, along with the stratigraphic deposition time during fault activity, the average fault activity rate for each seismic profile at different periods was calculated: V = / T, where V represents the average fault activity rate. Indicates the horizontal apparent displacement of the fault. The vertical apparent fault displacement is represented by T, which represents the sedimentary time of the strata during the period of fault activity.

7. A three-dimensional seismic assessment system for the activity of syn-sedimentary faults on gentle slopes of carbonate rocks, characterized in that, This system is used in a three-dimensional seismic assessment method for the activity of syn-sedimentary faults on gentle slopes of carbonate rocks as described in any one of claims 1-6, the system comprising: The layer interpretation module is used to interpret fault planes of syn-sedimentary faults and to interpret stratigraphic layers of fault-related strata in 3D seismic data volumes, so as to obtain fault plane-stratigraphic layer interpretation results. The deposition time calculation module is used to select strata for assessing fault activity from the fault-stratum layer interpretation results based on pre-constructed assessment selection rules, and to construct a sequence stratigraphic framework using the selected strata; and to calculate the deposition time of each set of strata during the fault activity period using the time corresponding to the sequence stratigraphic boundaries of the sequence stratigraphic framework. The trajectory acquisition module is used to select multiple seismic profiles from the fault-stratum layer interpretation results based on pre-built seismic profile selection rules, and to acquire the fault trajectory line and stratum layer trajectory line passing through each seismic profile. The activity rate calculation module is used to calculate the average fault activity rate of each seismic profile at different periods by using the fault plane trajectory line and the stratigraphic plane trajectory line of each seismic profile, combined with the deposition time of each set of strata during the fault activity period. The seismic assessment module is used to assess the activity of syn-depositional faults based on the average fault activity rate at different times for each seismic profile, and to obtain seismic assessment results.

8. A computer device comprising a system memory and a processor, wherein the system memory stores a computer program, characterized in that, When the processor executes the computer program, it implements the steps of the method according to any one of claims 1 to 6.

9. A computer-readable storage medium having a computer program stored thereon, characterized in that, When the computer program is executed by a processor, it implements the steps of the method described in any one of claims 1 to 6.

10. A computer program product containing instructions, characterized in that, When the instructions are executed by a cluster of computer devices, the cluster of computer devices causes the cluster of computer devices to perform the method as described in any one of claims 1 to 6.

Images