Fault seismic interpretation method, device and equipment and storage medium

Abstract

The invention discloses a fault seismic interpretation method, device and equipment and a storage medium, and relates to the technical field of seismic interpretation. The method comprises the following steps: firstly, acquiring a reference datum plane, and reading data from a fault attribute body according to the reference datum plane to form a near-geological isochronous slice; converting a fault line of the seismic time into relative time of a reference surface, and projecting the relative time to a near-geological isochronous slice; fault interpretation is carried out on the near-geological isochronous slice, and a near-geological isochronous fault line is obtained; recovering the near-geological isochronous fault line into a fault line of earthquake time; projecting a fault line of seismic time to a seismic section, referring to a projection point, performing conventional section fault interpretation, and storing a result; completing fault interpretation and closing, and obtaining an interpretation result of a final fault line. The method has the advantages of being high in operability and easy and rapid to implement, the workload of fault interpretation is greatly reduced, the recognition and interpretation precision is improved, and the method can be widely applied to fault interpretation work and software.

Description

Technical Field

[0001] The present invention relates to the technical field of structural interpretation of petroleum geophysical three-dimensional seismic data, and in particular to a fault seismic interpretation method, device, equipment and storage medium that can be used to improve the accuracy and efficiency of fault interpretation in geological conditions where the stratum structure is relatively large and the faults are complex. Background Art

[0002] Fault interpretation is the basis of seismic data interpretation. Good fault interpretation is directly related to structural accuracy. How to achieve accurate and rapid interpretation of faults through seismic data has always been a concern for interpreters. In recent years, with the development of computer technology, fault identification and interpretation technology has made great progress, and a variety of fault identification methods have emerged. Fault interpretation has gradually developed from single manual identification to auxiliary identification through multiple methods, which has greatly improved the accuracy of fault interpretation.

[0003] At present, fault identification methods can be generally divided into three categories: ① Conventional seismic profile fault identification and interpretation methods, including seismic profile identification and calibration of stratigraphic comparison breakpoints at well points, etc. For details, see the article "Three-dimensional Fault Combination" by Li Dong'an et al. in "Petroleum Geophysical Exploration" in 1992; ② Fault identification and interpretation methods based on seismic discontinuity attribute slicing. Commonly used discontinuity attributes that are more sensitive to faults include coherence attributes, variance attributes, curvature attributes, curvature, maximum likelihood volume, etc.; ③ Automatic fault tracking and interpretation methods based on seismic discontinuity attributes, which realize the carving of faults in three-dimensional space through ant tracking, neural network, image segmentation and other algorithms.

[0004] In actual production work, manual interpretation of fault lines through seismic and attribute profiles or slices is still the most basic and main working method. Other fault identification and interpretation methods are used as auxiliary methods in actual production.

[0005] In order to improve the accuracy of profile fault interpretation, Amoco proposed the coherent body fault interpretation method in 1995, which caused a sensation in the industry. For details, see S. Bahorich M's article "Using Three-Dimensional Coherence for Stratigraphic and Structural Interpretation" in the "Proceedings of the 65th Annual Meeting of the American Society of Exploration Geophysicists" published in 1995. Coherent slice interpretation includes: a new attribute data body obtained from the seismic data body through correlation calculation; horizontal slicing of coherent data, and interpretation of fault lines through the discontinuous lateral distribution features clearly shown on the slices; projecting the fault lines on the horizontal slices onto the seismic profile to assist and guide the profile fault interpretation; and realizing three-dimensional fault interpretation and closure through mutual projection and quality control of seismic profile and horizontal slice fault lines.

[0006] Coherent slice tomographic interpretation technology has been rapidly and widely used, and various software systems have integrated corresponding functions. However, the current coherent slice tomographic interpretation can only be implemented on horizontal slices, and this horizontal slice tomographic interpretation has the following shortcomings and problems.

[0007] 1. Horizontal slices do not take into account the characteristics and laws of geological strata, and are a pure seismic technology. They are isochronous in seismic bodies, but diachronous in geology, and one slice shows seismic response characteristics of different geological ages. Because they are not combined with geological structural characteristics, the regularity of fault characteristics on horizontal slices is poor, and the fault combination relationship is not good ( Figure 2 ), therefore, fault interpretation is difficult and the interpretation efficiency is low. At the same time, the fault plane combination has strong multi-solution characteristics and low interpretation accuracy, and the results may not be consistent with geological laws.

[0008] 2. When the structural undulations in the study area are large, the time range covered by the target layer also becomes larger. It is necessary to interpret multiple horizontal slices layer by layer from shallow to deep to achieve full coverage of the relevant faults in a certain layer, which requires a lot of interpretation work.

[0009] Therefore, in view of the shortcomings of coherent horizontal slice fault interpretation technology and the need for complex structural interpretation accuracy, the coherent fault interpretation method needs to be improved and optimized. Summary of the invention

[0010] In order to overcome the defects and shortcomings in the above-mentioned prior art, the present invention provides a fault seismic interpretation method. The purpose of the present invention is to provide a fault interpretation method along the fault attribute slices of the structural layer (trend surface) to solve the problems of poor consistency of horizontal time slice images, poor geological laws of the fault plane and low interpretation accuracy, unclear fault relationships and strong combination of multiple solutions, as well as large workload of high and steep structure interpretation in conventional structural interpretation.

[0011] In order to solve the above problems existing in the prior art, the present invention is implemented through the following technical solutions.

[0012] A first aspect of the present invention provides a fault seismic interpretation method, the method comprising the following steps:

[0013] S1. Select reference horizons, perform interpolation, smoothing and sampling, and obtain reference datum surfaces;

[0014] S2, reading data from the fault attribute body according to the reference datum surface obtained in step S1 to form a near geological isochronal slice;

[0015] S3, converting the fault lines of earthquake time into time relative to the reference datum and projecting them onto near geological isochronal slices;

[0016] S4, performing fault interpretation on the near geological isochronous slice obtained in step S3 to obtain near geological isochronous fault lines;

[0017] S5, restoring the near geological isochronous fault line obtained in step S4 to the fault line of the earthquake time;

[0018] S6, projecting the fault line of the earthquake time obtained in step S5 onto the seismic section, performing conventional section fault interpretation with reference to the projection point, and saving the result;

[0019] S7. Iterate steps S2 to S6 to complete fault interpretation and closure, and obtain the final fault line interpretation results.

[0020] Further preferably, in step S1, the selected reference layer is interpolated and smoothed according to the seismic grid to obtain a three-dimensional full-coverage continuous reference surface.

[0021] More preferably, in step S1, the method for smoothing the selected reference layer according to the seismic grid and the principles for determining the parameters include: the plane diagram of the structural mutation anomaly attributes has a certain signal-to-noise ratio, the anomaly characteristics are clear and easy to identify, the trend is consistent with the fault linear morphology, the anomaly and the fault interpreted on the seismic section have a good degree of consistency, and at the same time, the plane has a relatively consistent macroscopic law and conforms to the regional geomechanical background.

[0022] More preferably, in step S1, the following formula (1) is used to resample the continuous reference surface to obtain the reference datum surface, so that its time value is an integer multiple of the fault sensitive attribute sampling rate, so that data can be quickly extracted from the attribute body;

[0023] th′ x,y =round(th x,y / tr)×tr (1);

[0024] Where th′ represents the reference datum surface time, th represents the continuous reference surface time, tr represents the fault sensitive attribute volume sampling rate, x represents the line number, y represents the channel number, and round represents rounding.

[0025] Further preferably, in step S2, data is read from the fault sensitive attribute body according to the reference datum plane and the time shift; the data is displayed in the interpretation software in the form of horizontal slices to form near geological isochronal slices at a time Δt from the reference layer.

[0026] More preferably, the selected time shift amount is an integer multiple of the sampling rate.

[0027] Further preferably, in step S3, the fault line of the earthquake time is converted into the relative time of the reference datum plane, specifically, the fault line of the existing earthquake time is converted into the relative time of the reference datum plane at time Δt according to the following formula (2), and projected onto the near geological isochronous slice;

[0028] tf′ x,y =tf x,y +Δt-th′ x,y (2);

[0029] Where tf′ represents the time of the fault line relative to the reference layer at time Δt; tf represents the earthquake time of the fault line; Δt represents the time displacement; x represents the line number; y represents the track number; and th′ represents the reference datum plane time.

[0030] Further preferably, in step S5, the nearly geological isochronous fault line is restored to the fault line of the earthquake time, specifically, the nearly geological isochronous fault line is restored to the fault line of the earthquake time by calculating according to the following formula (3);

[0031] tf x,y =tf′ x,y -Δt+th′ x,y (3);

[0032] Where tf′ represents the time of the fault line relative to the reference layer at time Δt; tf represents the earthquake time of the fault line; Δt represents the time displacement; x represents the line number; y represents the track number; and th′ represents the reference datum plane time.

[0033] A second aspect of the present invention provides a fault seismic interpretation device, the device comprising:

[0034] The reference datum plane acquisition module is used to select a reference layer, perform interpolation, smoothing and sampling processing on the selected reference layer, and obtain the reference datum plane;

[0035] A near-geological isochronous slice acquisition module is used to read data from the fault attribute body according to the reference datum plane obtained by the reference datum plane acquisition module to form a near-geological isochronous slice;

[0036] A projection module is used to convert the fault lines of earthquake time into time relative to the reference datum and project them onto near geological isochron slices;

[0037] A geological isochronous fault line acquisition module is used to perform fault interpretation on the near geological isochronous slices obtained by the projection module to obtain near geological isochronous fault lines;

[0038] An earthquake time fault line acquisition module is used to restore the near geological isochronous fault lines acquired by the geological isochronous fault line acquisition module to fault lines of earthquake time;

[0039] A conventional profile fault interpretation module is used to project the fault line of the earthquake time obtained by the earthquake time fault line acquisition module onto the earthquake profile, refer to the projection point, perform conventional profile fault interpretation, and save the result;

[0040] The final interpretation module is used to complete fault interpretation and closure based on the iterative processing of the above-mentioned near-geological isochronous slice acquisition module, projection module, geological isochronous fault line acquisition module, seismic time fault line acquisition module and conventional profile fault interpretation module, and obtain the final fault line interpretation results.

[0041] A third aspect of the present invention provides a computer device, comprising a memory and a processor, wherein the processor is used to run a program stored in the processor, wherein the program, when running, executes the fault seismic interpretation method described in the first aspect of the present invention.

[0042] A fourth aspect of the present invention provides a readable storage medium, which is used to store a program, wherein when the program is running, the readable storage medium is controlled to execute the fault seismic interpretation method described in the first aspect of the present invention.

[0043] Compared with the prior art, the beneficial technical effects brought by the present invention are as follows:

[0044] 1. The present invention aims to solve the problems of poor consistency of horizontal time slice images, poor geological laws of fault planes, low interpretation accuracy, unclear fault relationships, strong combination of multiple solutions, and large workload of high-steep structure interpretation in conventional structural interpretation, and provide a fault interpretation method along the fault attribute slice of the structural layer (trend surface).

[0045] 2. The present invention adopts near-geological isochronal slices of sensitive fault attributes, which better reflects the position, trend and combination relationship of the overall fault, can quickly obtain the fault interpretation line near the target layer, and improves the slice fault recognition ability and interpretation accuracy.

[0046] 3. The present invention performs seismic time recovery on near-geological isochronal slice fault lines and projects them onto seismic sections in real time, thereby achieving dynamic interpretation of plane and section fault lines, reducing the workload of fault interaction adjustment and improving fault closure.

[0047] 4. The present invention only needs to carry out fault interpretation near the reference horizon, without considering the interpretation workload of other faults within other time windows such as non-target horizons, thereby reducing the number of slice interpretations and improving the efficiency of fault interpretation.

[0048] 5. The present invention makes full use of the basic functions of existing fault interpretation software and can realize the functions of the present invention without the development of a large amount of code. The implementation method is practical and effective.

[0049] 6. Based on the characteristics of fault line data, the present invention samples the reference surface according to the sampling rate of the fault attribute body, realizes direct extraction of data from the attribute body to form slices, avoids the calculation of layer attribute extraction, and meets the needs of real-time interactive interpretation;

[0050] 7. The present invention has the characteristics of strong operability, simple implementation and rapidity, which greatly reduces the workload of fault interpretation and improves the accuracy of identification and interpretation, and can be widely used in fault interpretation work and software. BRIEF DESCRIPTION OF THE DRAWINGS

[0051] Figure 1 A flow chart of the fault seismic interpretation method of the present invention;

[0052] Figure 2 It is a time coherence volume slice diagram of 1000ms (corresponding time of the main body of the buried mountain top) in the embodiment of the present invention;

[0053] Figure 3 It is a slice diagram of the geological isochronous coherence volume near the top of the buried mountain in the embodiment of the present invention;

[0054] Figure 4 This is a polygonal interpretation result diagram of the buried mountain top fault in an embodiment of the present invention. DETAILED DESCRIPTION

[0055] In order to enable those skilled in the art to better understand the scheme of the present invention, the technical scheme in the embodiments of the present invention will be clearly and completely described below in conjunction with the drawings in the embodiments of the present invention. Obviously, the described embodiments are only part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without creative work should fall within the scope of protection of the present invention.

[0056] Example 1

[0057] As a preferred embodiment of the present invention, refer to the attached specification Figure 1 As shown, this embodiment discloses a fault seismic interpretation method, which includes the following steps:

[0058] S1. Select reference horizons, perform interpolation, smoothing and sampling, and obtain reference datum surfaces;

[0059] S2, reading data from the fault attribute body according to the reference datum surface obtained in step S1 to form a near geological isochronal slice;

[0060] S3, converting the fault lines of earthquake time into time relative to the reference datum and projecting them onto near geological isochronal slices;

[0061] S4, performing fault interpretation on the near geological isochronous slice obtained in step S3 to obtain near geological isochronous fault lines;

[0062] S5, restoring the near geological isochronous fault line obtained in step S4 to the fault line of the earthquake time;

[0063] S6, projecting the fault line of the earthquake time obtained in step S5 onto the seismic section, performing conventional section fault interpretation with reference to the projection point, and saving the result;

[0064] S7. Iterate the above steps S2 to S6 to complete the fault interpretation and closure, and obtain the final fault line interpretation results.

[0065] Example 2

[0066] As another preferred embodiment of the present invention, this embodiment further supplements and elaborates on the technical solution of the present invention on the basis of the above-mentioned embodiment 1. In this embodiment, in step S1, the selected reference horizon is interpolated and smoothed according to the seismic grid to obtain a three-dimensional full-coverage continuous reference surface.

[0067] In step S1, the method and parameter determination principles for smoothing the selected reference horizon according to the seismic grid include: the plane map of the structural mutation anomaly attributes has a certain signal-to-noise ratio, the anomaly characteristics are clear and easy to identify, the strike is consistent with the fault linear morphology, the anomaly and the fault interpreted on the seismic section have a good degree of consistency, and at the same time, the plane has a relatively consistent macroscopic law and conforms to the regional geomechanical background.

[0068] In step S1, the continuous reference surface is resampled and calculated using the following formula (1) to obtain the reference datum surface, so that its time value is an integer multiple of the fault sensitive attribute sampling rate, so that data can be quickly extracted from the attribute body;

[0069] th′ x,y =round(th x,y / tr)×tr (1);

[0070] Where th′ represents the reference datum surface time, th represents the continuous reference surface time, tr represents the fault sensitive attribute volume sampling rate, x represents the line number, y represents the channel number, and round represents rounding.

[0071] As another implementation of this embodiment, in step S2, data is read from the fault sensitive attribute body according to the reference datum plane and the time shift; the data is displayed in the interpretation software in the form of horizontal slices to form a near geological isochronous slice at a time Δt from the reference layer. Preferably, the selected time shift is an integer multiple of the sampling rate.

[0072] As another implementation of this embodiment, in step S3, the fault line of the earthquake time is converted into the relative time of the reference reference plane, specifically, the fault line of the existing earthquake time is converted into the relative time of the reference reference plane at time Δt according to the following formula (2), and projected onto the near geological isochronal slice;

[0073] tf′ x,y =tf x,y +Δt-th′ x,y (2);

[0074] Where tf′ represents the time of the fault line relative to the reference layer at time Δt; tf represents the earthquake time of the fault line; Δt represents the time displacement; x represents the line number; y represents the track number; and th′ represents the reference datum plane time.

[0075] As another implementation of this embodiment, in step S5, the nearly geological isochronous fault line is restored to the fault line of the earthquake time, specifically, the nearly geological isochronous fault line is restored to the fault line of the earthquake time by calculating according to the following formula (3);

[0076] tf x,y =tf′ x,y -Δt+th′ x,y (3);

[0077] Where tf′ represents the time of the fault line relative to the reference layer at time Δt; tf represents the earthquake time of the fault line; Δt represents the time displacement; x represents the line number; y represents the track number; and th′ represents the reference datum plane time.

[0078] Example 3

[0079] This embodiment is an ancient buried hill reservoir in Block A in eastern China. The formation of the ancient buried hill is mainly affected by the tectonic background and fault activity. The buried hill type is mainly single-sided hill. The controlling fault is generally steep, with a cross-section dip angle of more than 50 degrees. The amplitude, scale and plane distribution characteristics of the buried hill trap are controlled by a series of faults dipping from northwest to north. Multi-phase tectonic movements have created a complex fault and fracture system in the area, and the dip angle of the buried hill top is steep, the tectonic amplitude is large, and the tectonic time range is between 500ms and 3000ms. The fault interpretation is difficult and the workload is huge.

[0080] Figure 2This is a horizontal time slice corresponding to the main body of the buried mountain top in Block A at a time of 1000ms. Obviously, the horizontal time slices are isochronous in the seismic body, but diachronous in geology. The same slice presents seismic response characteristics of different geological ages, making the fault characteristics on the horizontal slice unclear and the combination relationship poor. If the faults near the buried mountain top are to be traversed through horizontal slices, it is necessary to complete the fault interpretation of 625 horizontal slices (the 500-3000ms time window covered by the buried mountain top layer, calculated at a sampling interval of 4 milliseconds). In the case of low fault interpretation accuracy, there is still a huge workload of interpretation.

[0081] like Figure 1 As shown in Figure 2, methods to improve the accuracy of seismic inversion include:

[0082] 1) Select a reference layer, perform interpolation, smoothing and resampling to obtain a reference datum surface;

[0083] First, a reference horizon is selected, and the selected reference horizon is interpolated and smoothed according to the seismic grid to obtain a continuous reference surface with full three-dimensional coverage;

[0084] Specifically, the principles for determining the selected smoothing method and parameters include: the plane map of the tectonic mutation anomaly attributes has a certain signal-to-noise ratio, the anomaly features are clear and easy to identify, the strike is consistent with the fault linear morphology, the anomaly and the fault interpreted on the seismic profile have a good degree of consistency, and at the same time, the plane has a relatively consistent macroscopic law and conforms to the regional geomechanical background.

[0085] Secondly, for the continuous reference surface, the reference datum surface is obtained by resampling calculation using formula (1), so that its time value is an integer multiple of the sampling rate of the fault sensitive attribute body, so that data can be quickly extracted from the attribute body;

[0086] th′ x,y =round(th x,y / tr)×tr (1);

[0087] Where th′ represents the reference datum surface time, th represents the continuous reference surface time, tr represents the fault sensitive attribute volume sampling rate, x represents the line number, y represents the channel number, and round represents rounding.

[0088] In the embodiment of the present invention, the target layer is the top of the buried mountain, so the top of the buried mountain is selected as the reference layer, and the reference datum plane is obtained according to the above steps.

[0089] 2) Read data from the fault attribute volume according to the reference datum to form near-geological isochronal slices;

[0090] First, select the time shift of the slice, according to the reference plane and the time shift (t′ x,y+Δt) reads data from the fault sensitive attribute body; the selected time shift is an integer multiple of the sampling rate;

[0091] Secondly, the data are displayed in the interpretation software in the form of horizontal slices to form near-geological isochronal slices at the time of the reference layer;

[0092] In the embodiment of the present invention, the original sampling rate of seismic data is 4ms. Since the target layer is the top of the buried mountain, the time shift Δt is first set to 0, and data is read from the fault sensitive attribute body according to the reference datum plane of the top of the buried mountain to display the isochronous slice of the top of the buried mountain.

[0093] Figure 3 This is a near geological isochronous slice diagram along the top of the buried mountain in Block A of the embodiment of the present invention. Obviously, the fault features of this diagram are clearer, and the overall block structure of the block and the combination relationship between the faults are displayed. Therefore, the near geological isochronous slice diagram can greatly improve the accuracy of fault line interpretation.

[0094] 3) Convert the fault lines of earthquake time to time relative to the reference datum and project them onto near geological isochronal slices;

[0095] Among them, the fault line of the existing earthquake time is converted into the relative time of the reference datum surface at time Δt according to formula (2) and projected on the near geological isochronal slice:

[0096] tf′ x,y =tf x,y +Δt-th′ x,y (2);

[0097] Where tf′ represents the time of the fault line relative to the reference layer at time Δt; tf represents the earthquake time of the fault line; Δt represents the time displacement; x represents the line number; y represents the track number; and th′ represents the reference datum plane time.

[0098] 4) Perform fault interpretation on near geological isochronous slices to obtain near geological isochronous fault lines;

[0099] 5) Restore the nearly geological isochronous fault lines to the fault lines of the earthquake time and save the results;

[0100] The software system calculates according to formula (3) to restore the nearly geological isochronous fault line to the fault line of the earthquake time;

[0101] tf x,y =tf′ x,y -Δt+th′ x,y (3);

[0102] Where tf′ represents the time of the fault line relative to the reference layer at time Δt; tf represents the earthquake time of the fault line; Δt represents the time displacement; x represents the line number; y represents the track number; and th′ represents the reference datum plane time.

[0103] 6) Project the fault line of the earthquake time onto the seismic section, refer to the projection point, perform conventional section fault interpretation, and save the results;

[0104] 7) Iterate the above steps to complete the fault interpretation and closure and obtain the final fault line interpretation results.

[0105] In the embodiment of the present invention, the strata and fault systems above the buried hill top and the strata and faults inside the buried hill belong to two systems respectively, and are very different. Therefore, there is no need to pay too much attention to the faults above the buried hill top. On the premise of satisfying the fault interpretation work on the buried hill top, it is basically only necessary to interpret the faults in the time window of 20ms above the buried hill top and 80ms below the top surface. According to the 100ms time window and 4 millisecond sampling interval, the fault interpretation includes a total of 25 horizontal slices, which greatly reduces the workload.

[0106] Figure 4 This is the polygonal interpretation result diagram of the buried mountain top fault finally obtained in the embodiment of the present invention. Figure 3 From the slice diagram of the near geological isochronous coherence volume, it can be seen that the two have good similarity, which shows that the present invention has a good geological effect.

[0107] Example 4

[0108] As another preferred embodiment of the present invention, this embodiment discloses a fault seismic interpretation device, which includes:

[0109] The reference datum plane acquisition module is used to select a reference layer, perform interpolation, smoothing and sampling processing on the selected reference layer, and obtain the reference datum plane;

[0110] A near-geological isochronous slice acquisition module is used to read data from the fault attribute body according to the reference datum plane obtained by the reference datum plane acquisition module to form a near-geological isochronous slice;

[0111] A projection module is used to convert the fault lines of earthquake time into time relative to the reference datum and project them onto near geological isochron slices;

[0112] A geological isochronous fault line acquisition module is used to perform fault interpretation on the near geological isochronous slices obtained by the projection module to obtain near geological isochronous fault lines;

[0113] An earthquake time fault line acquisition module is used to restore the near geological isochronous fault lines acquired by the geological isochronous fault line acquisition module to fault lines of earthquake time;

[0114] A conventional profile fault interpretation module is used to project the fault line of the earthquake time obtained by the earthquake time fault line acquisition module onto the earthquake profile, refer to the projection point, perform conventional profile fault interpretation, and save the result;

[0115] The final interpretation module is used to complete the fault interpretation and closure according to the iterative processing of the above modules, and obtain the final fault line interpretation results.

[0116] Among them, the reference datum acquisition module interpolates and smoothes the selected reference horizon according to the seismic grid to obtain a continuous reference surface with full three-dimensional coverage. The method and parameter determination principles for smoothing the selected reference horizon according to the seismic grid include: the plane map of the structural mutation anomaly attribute has a certain signal-to-noise ratio, the anomaly characteristics are clear and easy to identify, the strike is consistent with the fault linear morphology, the anomaly and the fault interpreted on the seismic section have a good degree of consistency, and at the same time, the plane has a relatively consistent macroscopic law and conforms to the regional geomechanical background.

[0117] Furthermore, the reference datum plane acquisition module resamples the continuous reference plane using the following formula (1) to obtain the reference datum plane, so that its time value is an integer multiple of the fault sensitive attribute sampling rate, so as to be able to quickly extract data from the attribute body;

[0118] th′ x,y =round(th x,y / tr)×tr (1);

[0119] Where th′ represents the reference datum surface time, th represents the continuous reference surface time, tr represents the fault sensitive attribute volume sampling rate, x represents the line number, y represents the channel number, and round represents rounding.

[0120] The near-geological isochronous slice acquisition module reads data from the fault sensitive attribute body according to the reference datum and the time shift, and displays the data in the form of horizontal slices to form near-geological isochronous slices at the time Δt from the reference layer. The selected time shift is an integer multiple of the sampling rate.

[0121] The projection module converts the fault line of the earthquake time into the relative time of the reference datum plane. Specifically, the fault line of the existing earthquake time is converted into the relative time of the reference datum plane at time Δt according to the following formula (2), and projected onto the near geological isochronal slice;

[0122] tf′ x,y =tf x,y +Δt-th′ x,y (2);

[0123] Where tf′ represents the time of the fault line relative to the reference layer at time Δt; tf represents the earthquake time of the fault line; Δt represents the time displacement; x represents the line number; y represents the track number; and th′ represents the reference datum plane time.

[0124] The earthquake time fault line acquisition module calculates according to the following formula (3) to restore the nearly geological isochronous fault line to the fault line of earthquake time;

[0125] tf x,y =tf′ x,y -Δt+th′ x,y (3);

[0126] Where tf′ represents the time of the fault line relative to the reference layer at time Δt; tf represents the earthquake time of the fault line; Δt represents the time displacement; x represents the line number; y represents the track number; and th′ represents the reference datum plane time.

[0127] Example 5

[0128] In order to achieve the above-mentioned purpose, according to another aspect of the present application, a computer device is also provided, including a memory, a processor, and a computer program stored in the memory and executable on the processor, and when the processor executes the computer program, the steps of the above-mentioned hospital logistics equipment economic analysis method based on historical data are implemented.

[0129] In this embodiment, the processor may be a central processing unit (CPU). The processor may also be other general-purpose processors, digital signal processors (DSP), application-specific integrated circuits (ASIC), field-programmable gate arrays (FPGA) or other programmable logic devices, discrete gate or transistor logic devices, discrete hardware components, or a combination of the above chips.

[0130] The memory, as a non-transitory computer-readable storage medium, can be used to store non-transitory software programs, non-transitory computer executable programs and units, such as the corresponding program units in the above method embodiments of the present invention. The processor executes various functional applications of the processor and works data processing by running the non-transitory software programs, instructions and modules stored in the memory, that is, implementing the method in the above method embodiments.

[0131] The memory may include a program storage area and a data storage area, wherein the program storage area may store an operating system, an application required for at least one function; the data storage area may store data created by the processor, etc. In addition, the memory may include a high-speed random access memory, and may also include a non-volatile memory, such as at least one disk storage device, a flash memory device, or other non-volatile solid-state storage device. In some embodiments, the memory may optionally include a memory remotely arranged relative to the processor, and these remote memories may be connected to the processor via a network. Examples of the above-mentioned network include, but are not limited to, the Internet, an intranet, a local area network, a mobile communication network, and combinations thereof.

[0132] The one or more units are stored in the memory, and when executed by the processor, the steps in the above-mentioned embodiment 1, embodiment 2 or embodiment 3 are performed.

[0133] Example 6

[0134] As another preferred embodiment of the present invention, this embodiment discloses a computer-readable storage medium on which a computer program is stored. When the computer program is executed by a processor, the steps in the above-mentioned embodiment 1, embodiment 2 or embodiment 3 are implemented.

[0135] The above description is only a preferred specific implementation manner of the present invention, but the protection scope of the present invention is not limited thereto. Any technician familiar with the technical field can make equivalent replacements or changes according to the technical scheme and inventive concept of the present invention within the technical scope disclosed by the present invention, which should be covered by the protection scope of the present invention.

Claims

1. A fault seismic interpretation method, characterized in that: The method comprises the following steps, S1. Select reference horizons, perform interpolation, smoothing and sampling, and obtain reference datum surfaces; S2, reading data from the fault attribute body according to the reference datum surface obtained in step S1 to form a near geological isochronal slice; S3, converting the fault lines of earthquake time into time relative to the reference datum and projecting them onto near geological isochronal slices; S4, performing fault interpretation on the near geological isochronous slice obtained in step S3 to obtain near geological isochronous fault lines; S5, restoring the near geological isochronous fault line obtained in step S4 to the fault line of the earthquake time; S6, projecting the fault line of the earthquake time obtained in step S5 onto the seismic section, performing conventional section fault interpretation with reference to the projection point, and saving the result; S7. Iterate the above steps S2 to S6 to complete the fault interpretation and closure, and obtain the final fault line interpretation results.

2. A fault seismic interpretation method according to claim 1, characterized in that: In step S1, the selected reference layer is interpolated and smoothed according to the seismic grid to obtain a continuous reference surface with full three-dimensional coverage.

3. A fault seismic interpretation method according to claim 1 or 2, characterized in that: In step S1, the method and parameter determination principles for smoothing the selected reference horizon according to the seismic grid include: the plane map of the structural mutation anomaly attributes has a certain signal-to-noise ratio, the anomaly characteristics are clear and easy to identify, the strike is consistent with the fault linear morphology, the anomaly and the fault interpreted on the seismic section have a good degree of consistency, and at the same time, the plane has a relatively consistent macroscopic law and conforms to the regional geomechanical background.

4. A fault seismic interpretation method according to claim 1 or 2, characterized in that: In step S1, the continuous reference surface is resampled and calculated using the following formula (1) to obtain the reference datum surface, so that its time value is an integer multiple of the fault sensitive attribute sampling rate, so that data can be quickly extracted from the attribute body; th′ x,y =round(th x,y / tr)×tr (1); Where th′ represents the reference datum surface time, th represents the continuous reference surface time, tr represents the fault sensitive attribute volume sampling rate, x represents the line number, y represents the channel number, and round represents rounding.

5. A fault seismic interpretation method according to claim 1, characterized in that: In step S2, data are read from the fault-sensitive attribute body according to the reference datum plane and the time displacement; the data are displayed in the interpretation software in the form of horizontal slices to form near-geological isochronal slices at a time Δt from the reference layer.

6. A fault seismic interpretation method according to claim 5, characterized in that: , the selected time shift is an integer multiple of the sampling rate.

7. A fault seismic interpretation method according to claim 1, 2, 5 or 6, characterized in that: In step S3, the fault line of the earthquake time is converted into the relative time of the reference datum plane, specifically, the fault line of the existing earthquake time is converted into the relative time of the reference datum plane at time Δt according to the following formula (2), and projected onto the near geological isochronal slice; tf′ x,y =tf x,y +Δt-th′ x,y (2); Where tf′ represents the time of the fault line relative to the reference layer at time Δt; tf represents the earthquake time of the fault line; Δt represents the time displacement; x represents the line number; y represents the track number; and th′ represents the reference datum plane time.

8. A fault seismic interpretation method according to claim 1, 2, 5 or 6, characterized in that: In step S5, the nearly geological isochronous fault line is restored to the fault line of the earthquake time, specifically, the nearly geological isochronous fault line is restored to the fault line of the earthquake time by calculation according to the following formula (3); tf x,y =tf′ x,y -Δt+th′ x,y (3); Where tf′ represents the time of the fault line relative to the reference layer at time Δt; tf represents the earthquake time of the fault line; Δt represents the time displacement; x represents the line number; y represents the track number; and th′ represents the reference datum plane time.

9. A fault seismic interpretation device, characterized in that: The device includes The reference datum plane acquisition module is used to select a reference layer, perform interpolation, smoothing and sampling processing on the selected reference layer, and obtain the reference datum plane; A near-geological isochronous slice acquisition module is used to read data from the fault attribute body according to the reference datum plane obtained by the reference datum plane acquisition module to form a near-geological isochronous slice; A projection module is used to convert the fault lines of earthquake time into time relative to the reference datum and project them onto near geological isochron slices; A geological isochronous fault line acquisition module is used to perform fault interpretation on the near geological isochronous slices obtained by the projection module to obtain near geological isochronous fault lines; An earthquake time fault line acquisition module is used to restore the near geological isochronous fault lines acquired by the geological isochronous fault line acquisition module to fault lines of earthquake time; A conventional profile fault interpretation module is used to project the fault line of the earthquake time obtained by the earthquake time fault line acquisition module onto the earthquake profile, refer to the projection point, perform conventional profile fault interpretation, and save the result; The final interpretation module is used to complete fault interpretation and closure based on the iterative processing of the above-mentioned near-geological isochronous slice acquisition module, projection module, geological isochronous fault line acquisition module, seismic time fault line acquisition module and conventional profile fault interpretation module, and obtain the final fault line interpretation results.

10. A computer device, characterized in that: The method comprises a processor, an input device, an output device and a memory, wherein the processor, the input device, the output device and the memory are interconnected, wherein the memory is used to store a computer program, the computer program comprises program instructions, and the processor is configured to call the program instructions to execute the method according to any one of claims 1 to 8.

11. A computer-readable storage medium, characterized in that: The computer-readable storage medium stores a computer program, wherein the computer program includes program instructions, and when the program instructions are executed by a processor, the processor executes the method according to any one of claims 1 to 8.

Images