Fault enhancement method and fault development interpretation method, storage medium, and electronic device

Abstract

The present application relates to the technical field of oil and gas exploration, in particular to a fault enhancement method and a fault development interpretation method, a storage medium and an electronic device, which solve the problems of large amount of calculation, non-obvious enhancement effect of small faults and poor fault continuity in the prior art method for fault enhancement of seismic data with poor quality and weak fault information; the method comprises the following steps: obtaining three-dimensional seismic data of a target layer section, calculating an edge detection volume and performing fault enhancement processing to obtain a fault enhancement volume; sequentially performing variance operation and ant tracking operation on the fault enhancement volume, wherein the angle is controlled according to the fault development direction, the continuity of the fault body along the fault development direction is enhanced, and the three-dimensional seismic attribute volume of the target layer section after enhancement is obtained; the purpose of improving the signal-to-noise ratio and imaging quality of seismic data, further clearly depicting faults and facilitating subsequent fault identification and extraction is achieved.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of oil and gas exploration, in particular to a fault enhancement method and a fault development interpretation method, a storage medium and an electronic device. BACKGROUND

[0002] Faults are widely developed structural forms in tectonic movement. Because of the existence of tectonic stress factors, underground rock layers will inevitably be broken, thereby generating faults. Faults are one of the most important structures of the crust. In the field of oil and gas exploration technology, faults can not only serve as a channel for oil and gas migration but also as a boundary of fault block oil and gas fields, thereby effectively controlling the distribution of oil and gas fields. Therefore, accurate identification of faults is of great significance to the exploration and development of oil and gas fields.

[0003] Seismic exploration generally includes three steps of seismic data acquisition, seismic data processing and seismic data interpretation. During the acquisition and processing of seismic data, due to the influence of many factors, the signal-to-noise ratio of the seismic data may be low and the imaging quality may be poor, thereby leading to the fault information in the seismic data being relatively blurred, which is not conducive to the detection and identification of faults, thereby affecting the accuracy of seismic data interpretation. Therefore, through the fault enhancement technology, the fault in the seismic data is enhanced in signal while reducing noise, so that the fault is depicted more clearly, and then the enhanced fault is detected and identified, thereby improving the accuracy of seismic data interpretation.

[0004] However, the current fault enhancement method mostly enhances the seismic coherence volume through a filtering algorithm, and then enhances the fault information in the seismic coherence volume. When this method is used to enhance the fault in the seismic data with low signal-to-noise ratio, poor imaging quality and weak fault information, there are problems of large amount of calculation, not obvious enhancement effect on small faults and difficulty in solving the problem of poor fault continuity, thereby hindering the accuracy of subsequent fault identification and detection.

[0005] Therefore, based on the above problems, the present application provides a fault enhancement method for optimizing processing on the basis of fault enhancement, and a fault development interpretation method, a storage medium and an electronic device based on the fault enhancement method. SUMMARY

[0006] The application aims at the above problems, and provides an optimized fault enhancement processing method, device, storage medium and electronic equipment, which is based on fault enhancement technology, increases an optimization processing step, performs optimization processing based on a fault enhancement body, and sequentially performs variance operation and ant tracking operation on the fault enhancement body obtained after fault enhancement processing, so as to solve the problems of large calculation amount, non-obvious enhancement effect of small faults and poor fault continuity in the prior art, improve the signal-to-noise ratio and imaging quality of seismic data, further clearly depict faults, and facilitate subsequent fault identification and extraction.

[0007] The technical scheme adopted by the application is as follows.

[0008] To achieve the above object, in a first aspect, the application provides a fault enhancement method, which comprises the following steps.

[0009] acquiring three-dimensional seismic data of a target layer section;

[0010] calculating an edge detection body according to the three-dimensional seismic data, and performing fault enhancement processing on the edge detection body to obtain a fault enhancement body;

[0011] sequentially performing variance operation and ant tracking operation on the fault enhancement body, wherein, when performing ant tracking operation, angle control is performed according to a fault development direction to retain a fault body along the fault development direction, and continuity enhancement processing is performed on the fault body to obtain a three-dimensional seismic attribute body of the target layer section after enhancement.

[0012] According to the embodiments of the application, in the fault enhancement method, the acquisition of the three-dimensional seismic data of the target layer section comprises the following steps.

[0013] acquiring seismic data of a target layer section and calculating a structure guide body of the seismic data;

[0014] performing preprocessing on the seismic data according to the structure guide body of the seismic data to obtain three-dimensional seismic data of the target layer section.

[0015] According to the embodiments of the application, in the fault enhancement method, the calculation of the structure guide body of the seismic data comprises the following steps.

[0016] amplitude and waveform information of the seismic data are analyzed by scanning a time window, a structure guide body representing azimuth information of the seismic data is established, and the structure guide body is taken as the structure guide body of the seismic data.

[0017] According to the embodiment of the present application, optionally, in the fault enhancement method, the three-dimensional seismic data of the target layer section is obtained by preprocessing the seismic data according to the structure guide volume of the seismic data, comprising:

[0018] In the case that the original geological shape of the seismic data is kept unchanged according to the structure guide volume of the seismic data, the random noise interference of the seismic data is removed to obtain the denoised seismic data;

[0019] The denoised seismic data is subjected to the dip angle guide filtering processing to highlight the fracture detail features, and the filtered seismic data is obtained to obtain the three-dimensional seismic data of the target layer section.

[0020] According to the embodiment of the present application, optionally, in the fault enhancement method, the seismic data comprises a seismic subvolume;

[0021] Before the step of removing the random noise interference of the seismic data, the method further comprises:

[0022] The obtained seismic data of the target layer section is sequentially subjected to the bad channel processing, the cutting of the abnormal shallow layer region and the cutting processing to establish the seismic subvolume.

[0023] According to the embodiment of the present application, optionally, in the fault enhancement method, the edge detection volume is calculated according to the three-dimensional seismic data, comprising:

[0024] The amplitude and waveform information of the three-dimensional seismic data is analyzed by scanning the time window, the structure guide volume representing the azimuth information of the three-dimensional seismic data is established as the structure guide volume of the three-dimensional seismic data;

[0025] The edge detection algorithm along the layer edge is adopted to perform the edge detection on the three-dimensional seismic data, and the discontinuity features of the three-dimensional seismic data are imaged to obtain the edge detection volume, wherein the scanning window during the edge detection is along the distribution direction of the stratum according to the structure guide volume of the three-dimensional seismic data.

[0026] According to the embodiment of the present application, optionally, in the fault enhancement method, the fault enhancement volume is sequentially subjected to the variance operation and the ant tracking operation, wherein during the ant tracking operation, the angle control is performed according to the fault development direction to retain the fracture volume along the fault development direction, and the continuity enhancement processing is performed on the fracture volume to obtain the three-dimensional seismic attribute volume of the target layer section after the enhancement features, comprising:

[0027] The similarity between the adjacent seismic traces inside the fault enhancement volume is calculated based on the error analysis method, the region with the similarity lower than the given threshold value is enlarged to enhance the boundary features, and the variance volume of the target layer section is obtained;

[0028] Determine the fault development direction according to the structure-oriented volume and geological understanding of the three-dimensional seismic data;

[0029] Perform several ant tracking operations on the variance volume, wherein, when performing the ant tracking operation, angle control of the dip angle and the azimuth angle is performed according to the fault development direction to retain a fracture volume along the fault development direction, and the fracture volume is subjected to continuity enhancement processing by analyzing the curve closure of the time window to obtain the three-dimensional seismic attribute volume of the target layer section after enhancement of the features.

[0030] In a second aspect, the present application provides a fault development interpretation method, which comprises:

[0031] Using the fault enhancement method as described above, obtain the three-dimensional seismic attribute volume of all fault development directions in the target layer section after enhancement of the features;

[0032] Extract the along-layer slices of all three-dimensional seismic attribute volumes, and fuse them into the along-layer slices of the entire target layer section, which are used for interpreting the fault development of the target layer section.

[0033] In a third aspect, the present application provides a storage medium having a computer program stored thereon, which can be executed by one or more processors to implement the steps of the method as described above.

[0034] In a fourth aspect, the present application provides an electronic device comprising a memory and a processor, wherein the memory has a computer program stored thereon, and the computer program is executed by the processor to perform the steps of the method as described above.

[0035] Compared with the prior art, one or more embodiments in the above solution can have the following advantages or beneficial effects:

[0036] 1. The fault enhancement method and the fault development interpretation method, the storage medium and the electronic device provided by the present application are based on the fault enhancement technology, increase the optimization processing step, perform the optimization processing based on the fault enhancement volume, and can better depict the secondary fractures compared with the coherence volume; the variance operation and the ant tracking operation are sequentially performed on the fault enhancement volume obtained after the fault enhancement processing, so that the fault image is clearer and has higher resolution; and when the ant tracking operation is performed, the angle control and the continuity enhancement processing are performed according to the fault development direction to obtain the three-dimensional seismic attribute volume of the target layer section after enhancement of the features, improve the signal-to-noise ratio and the imaging quality of the seismic data, enhance the fault continuity, further clearly depict the faults in the target layer section, and facilitate the subsequent accurate identification and evaluation of the faults; a complete and well-applied fault enhancement method is established, the processing flow of the fault development interpretation method is improved, and the method has more advantages in the detection and identification of the small faults.

[0037] 2. The application, after obtaining the seismic data of the target layer section, calculates the structure guide body of the seismic data, and adds the step of calculating the structure guide body of the seismic data before preprocessing, so that the original seismic characteristics will not be changed when the seismic data is preprocessed, thereby avoiding the weakening of fault information.

[0038] 3. The application, according to the structure guide body of the seismic data, pre-processes the seismic data to obtain three-dimensional seismic data of the target layer section, improves the signal-to-noise ratio of the seismic data, and improves the imaging quality of the seismic data.

[0039] 4. The application, before removing the random noise interference of the seismic data, sequentially performs bad channel processing, cuts off the abnormal area of the shallow layer, and performs cutting processing to establish the seismic sub-body, so that larger seismic data can be processed, the subsequent denoising and filtering are facilitated, and the processing efficiency and speed of preprocessing are improved.

[0040] 5. The application, using the fault enhancement method of the application for fault development interpretation, after obtaining the three-dimensional seismic attribute body of all fault development directions, extracting the along-layer slice of all three-dimensional seismic attribute bodies, and fusing the along-layer slice of the entire target layer section, the fault interpretation and extraction are facilitated, and the accuracy of fault interpretation is improved. BRIEF DESCRIPTION OF DRAWINGS

[0041] In the following, the application will be described in more detail based on the embodiments and with reference to the accompanying drawings.

[0042] Figure 1 A flowchart of a fault enhancement method provided for the first embodiment of the application.

[0043] Figure 2 A time slice diagram of seismic data in step 101.1 of the fault enhancement method provided for the first embodiment of the application.

[0044] Figure 3 A time slice diagram of the structure guide body in step 101.2 of the fault enhancement method provided for the first embodiment of the application.

[0045] Figure 4 A time slice diagram of three-dimensional seismic data in step 101.3 of the fault enhancement method provided for the first embodiment of the application.

[0046] Figure 5 A time slice diagram of the structure guide body in step 102.1 of the fault enhancement method provided for the first embodiment of the application.

[0047] Figure 6The time slice graph of the edge detection body in step 102.2 of the fault enhancement method provided by the first embodiment of the present application.

[0048] Figure 7 The time slice graph of the fault enhancement body in step 102.3 of the fault enhancement method provided by the first embodiment of the present application.

[0049] Figure 8 The along-layer slice graph of the fault enhancement body in step 103.1 of the fault enhancement method provided by the first embodiment of the present application.

[0050] Figure 9 The along-layer slice graph of the three-dimensional seismic attribute body in step 103.3 of the fault enhancement method provided by the first embodiment of the present application.

[0051] Figure 10 The flowchart of the fault development interpretation method provided by the second embodiment of the present application.

[0052] In the drawings, the same components are denoted by the same reference numerals, and the drawings are not drawn according to the actual scale. DETAILED DESCRIPTION

[0053] The embodiments of the present application will be described in detail below with reference to the accompanying drawings and embodiments, so that how the present application applies technical means to solve technical problems and achieves corresponding technical effects can be fully understood and implemented. The embodiments of the present application and various features in the embodiments can be combined with each other without conflict, and the formed technical solutions are all within the protection scope of the present application.

[0054] Example One

[0055] Please refer to Figures 1 to 8 The present embodiment provides a fault enhancement method applicable to an electronic device. When the method is applied to the electronic device, steps 101 to 103 are performed.

[0056] Step 101: Obtain three-dimensional seismic data of a target layer section.

[0057] Step 101.1: Obtain seismic data of the target layer section.

[0058] In the present embodiment, actual seismic data of a certain area is taken as an example to obtain seismic data containing a target layer section, such as Figure 2 The time slice graph of the seismic data is shown in the figure, and the seismic data is original seismic data, which is obtained from Figure 2It can be seen from the original seismic data that the fault information is weak, the fault characteristics are not obvious, and there are a lot of random noises, which brings great trouble to the detection and identification of the fault. The fault refers to the surface damage or surface rheological zone of the rock layer or rock mass in the crust under the action of stress, and the rock blocks on both sides of the fault have obvious displacement;

[0059] Step 101.2: Calculate the structure guide body of the seismic data;

[0060] Specifically, the structure guide body representing the azimuth information of the seismic data is established by scanning the amplitude and waveform information of the seismic data, and is used as the structure guide body of the seismic data.

[0061] In this embodiment, the structure guide body representing the azimuth information of the seismic data is established to maintain the original geological form unchanged in the subsequent preprocessing denoising and edge detection body calculation process, such as Figure 3 The time slice diagram of the structure guide body of the seismic data is shown in the figure; by adding the step of calculating the structure guide body of the original seismic data, the original seismic characteristics will not be changed when the original seismic data is preprocessed, so that the situation that the fault information is weakened is not easy to occur when the original seismic data is preprocessed; from Figure 3 It can be seen from the structure guide body that the structure guide body contains the azimuth information of the original seismic data, and the fault characteristics are obvious, but are still seriously affected by the noise;

[0062] Step 101.3: Preprocessing the seismic data according to the structure guide body of the seismic data to obtain three-dimensional seismic data of the target layer; wherein, the seismic data is preprocessed to optimize the input seismic data, eliminate the influence of noise and other interference as much as possible, improve the signal-to-noise ratio of the seismic data, and highlight the fracture detail characteristics;

[0063] In this embodiment, the seismic data containing the target layer obtained in step 101.1 is further preprocessed to improve the resolution, and three-dimensional seismic data containing the target layer with high quality is obtained;

[0064] Step 101.3.1: sequentially performing bad channel removal, cutting off the abnormal area in the shallow layer, and cutting processing on the obtained seismic data of the target layer to establish a seismic sub-body;

[0065] In this embodiment, the seismic channel containing the target layer obtained in step 101.1 is first removed, then the abnormal area in the shallow layer caused by acquisition or other reasons is cut off, and then cutting is performed to cut off the top or bottom area, and only the seismic data containing the target layer is reserved as the seismic sub-body. The above processing can be performed on the larger original seismic data to improve the processing efficiency and speed of the preprocessing.

[0066] It should be noted that this step can be selected according to the actual size of the original seismic data, and when this step is not performed, the subsequent other steps of pre-processing the seismic data obtained in step 101.1 are directly performed;

[0067] Step 101.3.2: removing random noise interference of the seismic sub-volume while keeping the original geological form of the seismic sub-volume unchanged according to the structure guide volume of the seismic data, to obtain denoised seismic data;

[0068] In this embodiment, the structure guide volume of the seismic data obtained in step 101.2 and the seismic sub-volume obtained in step 101.3.1 are taken as input data together, and random noise and other interferences are removed therefrom to obtain denoised seismic data;

[0069] Step 101.3.3: performing dip guide filtering processing on the denoised seismic data, retaining the original seismic information while highlighting the fracture detail features to a certain extent, to obtain filtered seismic data, so as to obtain three-dimensional seismic data of the target layer section;

[0070] In this embodiment, appropriate calculation parameters are selected to perform filtering processing on the denoised seismic data, so as to obtain three-dimensional seismic data in which noise is removed and reflection events have enhanced lateral continuity and actual discontinuity information, as shown in FIG. 6. Figure 4 As can be seen from FIG. 6, compared with the time slice graph of the seismic data shown in FIG. 5, Figure 4 the white spot-like random noise in FIG. 6 is obviously reduced, and the signal-to-noise ratio of the slice is obviously improved, which indicates that the method for pre-processing the obtained original seismic data effectively removes random noise and other interferences, improves the signal-to-noise ratio of the seismic data, and has the effect of highlighting fracture detail features. Figure 2 Figure 4

[0071] Step 102: calculating an edge detection volume according to the three-dimensional seismic data, and performing fault enhancement processing on the edge detection volume to obtain a fault enhancement volume;

[0072] Step 102.1: calculating a structure guide volume of the three-dimensional seismic data;

[0073] Specifically, the amplitude and waveform information of the three-dimensional seismic data are analyzed by scanning a time window, and a structure guide volume representing azimuth information of the three-dimensional seismic data is established as the structure guide volume of the three-dimensional seismic data;

[0074] ​​In this embodiment, the structural guide volume representing the azimuth information of the three-dimensional seismic data obtained in step 101.3.3 is established, which facilitates the subsequent intuitive acquisition of the dominant azimuth of the fault development from the structural guide volume, and is suitable for the main fault which is relatively obvious; meanwhile, it also prepares for the angle control in the subsequent ant tracking operation in the optimization, as shown in Figure 5 FIG. 6 is a time slice diagram of the structural guide volume of the three-dimensional seismic data. Figure 5 As can be seen from Figure 4 the imaging quality of the three-dimensional seismic data corresponding to Figure 5 is improved,

[0075] Step 102.2: edge detection is performed on the three-dimensional seismic data by using a layer-following edge detection algorithm, and the discontinuity features of the three-dimensional seismic data are imaged to obtain an edge detection volume, wherein the scanning window in the edge detection is made to follow the distribution direction of the strata according to the structural guide volume of the three-dimensional seismic data.

[0076] In this embodiment, the layer-following edge detection algorithm is used to perform edge detection on the preprocessed three-dimensional seismic data, and the structural guide volume of the three-dimensional seismic data is added in the calculation process to ensure that the scanning window follows the distribution direction of the strata, and the edge detection volume is calculated, as shown in Figure 6 FIG. 7 is a time slice diagram of the edge detection volume.

[0077] Step 102.3: fault enhancement processing is performed on the edge detection volume to obtain a fault enhancement volume; the effective signal is enhanced, and the ability to depict the fault is increased.

[0078] In this embodiment, the fault enhancement processing is performed on the edge detection volume of step 102.2 to obtain the fault enhancement volume, as shown in Figure 7 FIG. 8 is a time slice diagram of the fault enhancement volume. Figure 8 FIG. 9 is a layer slice diagram of the fault enhancement volume.

[0079] As can be seen from Figure 9 and Figure 9 , the imaging quality of the layer slice of the fault enhancement volume is improved compared with the time slice, but the continuity of the fault is still poor. In order to further enhance the useful signal of the fault, make the fault depiction clearer, and make the boundary features of the fault more obvious and the continuity better, the optimization processing is performed on the basis of the fault enhancement volume.

[0080] Step 103: variance operation and ant tracking operation are sequentially performed on the fault enhancement volume, wherein in the ant tracking operation, the angle control is performed according to the fault development direction to retain the fault body along the fault development direction, and the continuity enhancement processing is performed on the fault body to obtain the three-dimensional seismic attribute volume of the target layer section after the enhancement features.

[0081] In this embodiment, since the fault information in the area is weak, although the preprocessing is performed before the fault enhancement processing, the signal-to-noise ratio of the seismic data is improved to some extent, and the imaging quality is improved, but the obtained fault enhancement body still cannot well depict the fault, the fault information is weak on the slice, and the continuity is poor, therefore, the optimization processing is continuously performed through step 103;

[0082] Step 103.1: Similarity between adjacent seismic traces in the fault enhancement body is calculated based on the error analysis method, and the area with similarity lower than a given threshold is enlarged to enhance the boundary feature, so as to obtain the variance body of the target layer section;

[0083] In this embodiment, the fault enhancement body obtained in step 102.3 is taken as the input data, the variance operation is performed, the similarity between adjacent seismic traces in the fault enhancement body is calculated, the discontinuous area, i.e., the area with similarity lower than a given threshold, is found out, and the area with poor similarity is enlarged to enhance the boundary feature, so as to obtain the variance body of the target layer section; wherein, the theoretical basis of the variance operation is the error analysis theory, and the similarity between adjacent seismic traces is used to describe the lateral heterogeneity of strata and lithology. When there is a fault or a local discontinuous change in the strata, the reflection characteristics of some seismic traces will be different from those of the nearby seismic traces, so as to cause the local discontinuity of the seismic traces. Thus, the fault or the discontinuous change information can be detected by detecting the difference degree between the seismic traces. The variance operation quantitatively processes the coherence attribute of the fault enhancement body, generates a new variance data body, highlights and emphasizes the irrelevance of the seismic data, and reflects the cracks, faults, river channels, lithological boundaries, etc. The boundary feature in this embodiment is essentially the area with poor similarity in the fault enhancement body, i.e., the discontinuous place. Enlarging the discontinuous place is to enhance the boundary feature of the discontinuous place, and the discontinuous place is the fracture area to be found, i.e., the fault in the target layer section;

[0084] Step 103.2: The fault development direction is determined according to the structure guide body of the three-dimensional seismic data and the geological understanding;

[0085] In this embodiment, the approximate fault development direction is determined according to the structure guide body of the three-dimensional seismic data and the geological understanding;

[0086] Step 103.3: Perform several ant-tracing operations on the variance volume to make the fault image clearer and have higher resolution. During the ant-tracing operation, the dip and azimuth angles are controlled according to the fault development direction to retain the fault body along the fault development direction. The fault body is then subjected to continuity enhancement processing by closing the curve in the analysis time window to obtain a 3D seismic attribute volume of the target layer segment after the enhanced features are obtained.

[0087] In this embodiment, three ant-tracing operations are performed on the variance volume. During the third ant-tracing operation, the dip and azimuth angles are controlled based on the analysis range of the scanning time window of the fault development direction determined in step 103.2. This reduces the incoherence caused by changes in the formation dip and azimuth, highlights the discontinuity caused by the fault, enhances the useful signal, and retains only the fracture volume in the fault development direction. The purpose of running the operation three times is to maximize the highlighting of the fault characteristics and increase their continuity. The number of ant-tracing operations can be selected according to actual conditions.

[0088] By analyzing the curve closure of the time window, the fracture body is subjected to continuity enhancement processing, and the useful signal is enhanced to the greatest extent to obtain better continuity, thus obtaining the 3D seismic attribute body of the target layer segment after the enhanced features, such as Figure 8 The following is a slice of the 3D seismic attribute volume along the layer. Example Two It can be seen that compared with Figures 2 to 10 The useful signals of the layer-by-layer slices of the fault enhancement body and the three-dimensional seismic attribute body are enhanced, the continuity is better, and the imaging quality is also improved. This shows that the application of this method in this implementation has achieved the purpose of improving the signal-to-noise ratio and imaging quality, further clearly depicting the fault, and improving the accuracy of subsequent fault identification and evaluation.

[0089] The fault enhancement method provided by the embodiment is based on the fault enhancement technology, an optimization processing step is added, and the optimization processing is performed based on the fault enhancement volume, so that the secondary faults can be better described than the coherent volume; the variance operation and the ant tracking operation are sequentially performed on the fault enhancement volume obtained after the fault enhancement processing, so that the fault image is clearer and has higher resolution; and when the ant tracking operation is performed, the angle control and the continuity enhancement processing are performed according to the fault development direction, the three-dimensional seismic attribute volume of the target layer section after the enhanced features is obtained, the signal-to-noise ratio and the imaging quality of the seismic data are improved, the fault continuity is enhanced, the faults in the target layer section are further clearly described, and the subsequent accurate identification and evaluation of the faults are facilitated; a complete fault enhancement method with good application effect is established, the processing flow of the fault development interpretation method is improved, and the method has more advantages in the detection and identification of the small faults; and by adding the step of calculating the structure-oriented volume of the original seismic data, the original seismic characteristics are not changed when the original seismic data is preprocessed, and the situation that the fault information is weakened is not easily caused when the preprocessing is performed; the embodiment solves the problems that the fault enhancement method for the seismic data with poor quality and weak fault information in the prior art has large calculation amount, the enhancement effect of the small faults is not obvious, and the poor fault continuity is difficult to overcome.

[0090] Figure 2

[0091] Referring to Figure 2 , the embodiment further provides a fault development interpretation method applicable to an electronic device, when the method is applied to the electronic device, steps 201 to 204 are performed.

[0092] Step 201: Obtain three-dimensional seismic data of a target layer section.

[0093] Step 201.1: Obtain seismic data of the target layer section.

[0094] In the embodiment, actual seismic data of a certain area is still taken as an example, and seismic data containing the target layer section is obtained, for example, Figure 3 The time slice diagram of the seismic data is shown in the figure, the seismic data is original seismic data, and from Figure 3 it can be seen that the fault information in the original seismic data is weak, the fault characteristics are not obvious, and a large amount of random noise is contained, which brings great trouble to the detection and identification of the faults.

[0095] Step 201.2: Calculate a structure-oriented volume of the seismic data.

[0096] Specifically, the amplitude and waveform information of the seismic data is analyzed by scanning a time window to establish a structure guide volume representing the azimuth information of the seismic data as a structure guide volume of the seismic data.

[0097] In this embodiment, the structure guide volume representing the azimuth information of the seismic data is established to maintain the original geological form unchanged in the subsequent preprocessing denoising and edge detection volume calculation process, for example, Figure 4 The time slice diagram of the structure guide volume of the seismic data is shown in FIG. 6. By adding the step of calculating the structure guide volume of the original seismic data, the original seismic characteristics are not changed when the original seismic data is preprocessed, so that the fault information is not easily weakened during the preprocessing. As shown in Figure 4 It can be seen from FIG. 6 that the structure guide volume contains the azimuth information of the original seismic data, and the fault features are obvious, but are still seriously affected by noise.

[0098] Step 201.3: Preprocessing the seismic data according to the structure guide volume of the seismic data to obtain three-dimensional seismic data of the target layer section;

[0099] In this embodiment, the seismic data containing the target layer section obtained in step 201.1 is further preprocessed to improve the resolution, and three-dimensional seismic data containing the target layer section with high quality is obtained.

[0100] Step 201.3.1: sequentially performing bad channel removal, cutting off abnormal areas in shallow layers and cutting processing on the obtained seismic data of the target layer section to establish a seismic sub-volume;

[0101] In this embodiment, the seismic channels containing the target layer section obtained in step 201.1 are first removed, then the abnormal areas in the shallow layers caused by acquisition or other reasons are cut off, and then cutting is performed to cut off the top or bottom area and only keep the seismic data containing the target layer section as a seismic sub-volume. The above processing can be performed on the larger original seismic data to improve the processing efficiency and speed of the preprocessing;

[0102] Step 201.3.2: removing random noise interference of the seismic sub-volume while maintaining the original geological form of the seismic sub-volume unchanged according to the structure guide volume of the seismic data to obtain denoised seismic data;

[0103] In this embodiment, the structure guide volume of the seismic data obtained in step 201.2 and the seismic sub-volume obtained in step 201.3.1 are taken as input data to remove random noise and other interference in the seismic data to obtain denoised seismic data.

[0104] Step 201.3.3: Perform dip-guided filtering on the denoised seismic data to retain the original seismic information while highlighting the detailed features of the fault to a certain extent, thereby obtaining filtered seismic data to obtain 3D seismic data of the target interval;

[0105] In this embodiment, appropriate calculation parameters are selected to filter the de-noised seismic data, thereby obtaining three-dimensional seismic data with noise removed and enhanced lateral continuity and actual discontinuity information of reflection events, such as Figure 2 The three-dimensional seismic data is shown as a time slice diagram; Figure 4 It can be seen that compared with Figure 5 A time slice diagram of the seismic data, Figure 5 The random noise in the form of white spots in the middle was significantly reduced, and the signal-to-noise ratio of the slices was significantly improved;

[0106] Step 202: calculating an edge detection volume according to the 3D seismic data, and performing fault enhancement processing on the edge detection volume to obtain a fault enhancement volume;

[0107] Step 202.1: Calculate the structural guidance volume of the 3D seismic data;

[0108] Specifically, by analyzing the amplitude and waveform information of the three-dimensional seismic data through a scanning time window, a structural guide body representing the azimuth information of the three-dimensional seismic data is established as the structural guide body of the three-dimensional seismic data;

[0109] In this embodiment, a structural guide body is established to represent the azimuth information of the three-dimensional seismic data obtained in step 201.3.3, so as to facilitate the subsequent intuitive acquisition of the dominant orientation of fault development from the structural guide body, which is suitable for more obvious main faults; at the same time, it also prepares for the angle control in the ant tracking operation during the subsequent optimization, such as Figure 4 Shown is a time slice diagram of the structural guide volume of the three-dimensional seismic data; Figure 5 It can be seen that the corresponding Figure 6 Improved imaging quality of 3D seismic data, Figure 7 The imaging quality has also been significantly improved;

[0110] Step 202.2: Perform edge detection on the 3D seismic data using a layer-wise edge detection algorithm, image discontinuity features in the 3D seismic data, and obtain an edge detection volume. The scanning window during edge detection is aligned with the direction of stratum distribution based on the structural guide volume of the 3D seismic data.

[0111] In this embodiment, the edge detection algorithm along layer edge is used to detect the edge of the preprocessed 3D seismic data, and the structure guide volume of the 3D seismic data is added in the calculation process to ensure that the scanning window is along the distribution direction of the stratum, and the edge detection volume is calculated, as shown in the following figure: Figure 8 The time slice graph of the edge detection volume is shown in the following figure:

[0112] Step 202.3: The fault enhancement processing is performed on the edge detection volume to obtain a fault enhancement volume;

[0113] In this embodiment, the fault enhancement processing is performed on the edge detection volume of step 202.2 to obtain a fault enhancement volume, as shown in the following figure: Figure 9 The time slice graph of the fault enhancement volume is shown in the following figure: Figure 9 The layer slice graph of the fault enhancement volume is shown in the following figure:

[0114] Step 203: The variance operation and the ant tracking operation are sequentially performed on the fault enhancement volume, wherein, when the ant tracking operation is performed, the angle control is performed according to the fault development direction to retain the fracture volume along the fault development direction, and the continuity enhancement processing is performed on the fracture volume to obtain the 3D seismic attribute volume of the target layer section after the enhancement feature;

[0115] Step 203.1: The similarity between adjacent seismic traces in the fault enhancement volume is calculated based on the error analysis method, the region with the similarity lower than the given threshold is enlarged to enhance the boundary feature, and the variance volume of the target layer section is obtained;

[0116] In this embodiment, the fault enhancement volume obtained in step 202.3 is taken as the input data to perform the variance operation, the similarity between adjacent seismic traces in the fault enhancement volume is calculated, the region with the similarity lower than the given threshold is enlarged to enhance the boundary feature, and the variance volume of the target layer section is obtained;

[0117] Step 203.2: The fault development direction is determined according to the structure guide volume of the 3D seismic data and the geological understanding;

[0118] In this embodiment, the approximate fault development direction represented by the structure guide volume of the 3D seismic data is combined with certain geological understanding to artificially determine a fault development direction;

[0119] Step 203.3: The ant tracking operation is performed on the variance volume for several times, wherein, when the ant tracking operation is performed, the angle control of the dip angle and the azimuth angle is performed according to the fault development direction to retain the fracture volume along the fault development direction, and the continuity enhancement processing is performed on the fracture volume by analyzing the curve closure of the time window to obtain the 3D seismic attribute volume of the target layer section after the enhancement feature;

[0120] In this embodiment, three ant tracking operations are performed on the variance volume. During the third ant tracking operation, the dip and azimuth angles are controlled based on the analysis range of the scanning time window in the fault development direction determined in step 203.2. This reduces the incoherence caused by changes in the formation dip and azimuth, highlights the discontinuity caused by the fault, enhances the useful signal, and retains only the fault volume in the fault development direction.

[0121] By analyzing the curve closure of the time window, the fracture body is subjected to continuity enhancement processing, and the useful signal is enhanced to the greatest extent to obtain better continuity, thus obtaining the 3D seismic attribute body of the target layer segment after the enhanced features, such as Figure 8 The following is a slice of the 3D seismic attribute volume along the layer. Example Three It can be seen that compared with Example Four The useful signals of the slices along the fault enhancement body and the slices along the 3D seismic attribute body are enhanced, the continuity is better, and the imaging quality is also improved;

[0122] Step 204: obtaining a 3D seismic attribute volume of all fault development directions in the target layer segment after feature enhancement;

[0123] In this embodiment, returning to step 203.2, another fault development direction is manually determined based on the approximate fault development direction represented by the structural guide body of the 3D seismic data and in combination with certain geological knowledge, to obtain a 3D seismic attribute body of the target layer segment with enhanced features of the other fault development direction;

[0124] Repeat steps 203.2 to 203.3 to obtain the three-dimensional seismic attribute volume of all fault development directions;

[0125] Step 205: extracting layer slices of all 3D seismic attribute volumes and merging them into layer slices of the entire target layer segment for interpreting the fault development of the target layer segment, thereby improving the accuracy of fault interpretation.

[0126] In this embodiment, the fusion is a layer-wise slice of the entire target layer segment, that is, a planar map of the target layer segment with a relatively high resolution for characterizing the development of faults. The fault development of the target layer segment can be interpreted based on the map.

[0127] This embodiment provides a fault development interpretation method, which interprets fault development by using a fault enhancement method based on fault enhancement technology and adding an optimization processing step. First, the original seismic data is subjected to fault enhancement and optimization processing. Then, after obtaining three-dimensional seismic attribute volumes of all fault development directions, layer-by-layer slices of all three-dimensional seismic attribute volumes are extracted and fused to obtain layer-by-layer slices of the entire target layer segment, which facilitates fault interpretation and extraction and improves the accuracy of fault interpretation.

[0128] ​

[0129] The embodiment provides a computer readable storage medium, such as a flash memory, a hard disk, a multimedia card, a card memory (for example, an SD or DX memory, etc.), a random access memory (RAM), a static random access memory (SRAM), a read only memory (ROM), an electrically erasable programmable read only memory (EEPROM), a programmable read only memory (PROM), a magnetic memory, a magnetic disk, an optical disk, a server, an App application market, etc., on which a computer program is stored, and the computer program can implement the following method steps when executed by a processor.

[0130] Step 301: obtaining three-dimensional seismic data of a target layer section;

[0131] Step 301.1: obtaining seismic data of a target layer section;

[0132] Step 301.2: calculating a structure guide volume of the seismic data;

[0133] Specifically, the amplitude and waveform information of the seismic data are analyzed by scanning a time window, a structure guide volume representing azimuth information of the seismic data is established as a structure guide volume of the seismic data;

[0134] Step 301.3: pre-processing the seismic data according to the structure guide volume of the seismic data to obtain three-dimensional seismic data of a target layer section;

[0135] Step 301.3.1: sequentially performing bad channel processing, cutting off an abnormal area in a shallow layer and cutting processing on the obtained seismic data of a target layer section to establish a seismic sub-volume;

[0136] Step 301.3.2: removing random noise interference of the seismic sub-volume without changing the original geological form of the seismic sub-volume according to the structure guide volume of the seismic data, to obtain denoised seismic data;

[0137] Step 301.3.3: performing dip angle guide filtering processing on the denoised seismic data to highlight fracture detail features, to obtain filtered seismic data, to obtain three-dimensional seismic data of a target layer section;

[0138] Step 302: calculating an edge detection volume according to the three-dimensional seismic data, and performing fault enhancement processing on the edge detection volume to obtain a fault enhancement volume;

[0139] Step 302.1: calculating a structure guide volume of the three-dimensional seismic data;

[0140] Specifically, by scanning the amplitude and waveform information of the three-dimensional seismic data, a structure guide volume representing the azimuth information of the three-dimensional seismic data is established as a structure guide volume of the three-dimensional seismic data;

[0141] Step 302.2: edge detection is performed on the three-dimensional seismic data by using a layer edge detection algorithm to image the discontinuity features of the three-dimensional seismic data, and an edge detection volume is obtained, wherein the scanning window during edge detection is along the distribution direction of the stratum according to the structure guide volume of the three-dimensional seismic data;

[0142] Step 302.3: a fault enhancement processing is performed on the edge detection volume to obtain a fault enhancement volume;

[0143] Step 303: a variance operation and an ant tracking operation are sequentially performed on the fault enhancement volume, wherein during the ant tracking operation, angle control is performed according to the fault development direction to retain a fracture volume along the fault development direction, and a continuity enhancement processing is performed on the fracture volume to obtain a three-dimensional seismic attribute volume of the target layer section after enhancement features;

[0144] Step 303.1: similarity between adjacent seismic traces inside the fault enhancement volume is calculated based on an error analysis method, and a region with a similarity lower than a given threshold is enlarged to enhance its boundary features, and a variance volume of the target layer section is obtained;

[0145] Step 303.2: a fault development direction is determined according to the structure guide volume of the three-dimensional seismic data and geological understanding;

[0146] Step 303.3: a plurality of ant tracking operations are performed on the variance volume, wherein during the ant tracking operation, angle control of dip angle and azimuth angle is performed according to the fault development direction to retain a fracture volume along the fault development direction, and a continuity enhancement processing is performed on the fracture volume by analyzing curve closure of a scanning window to obtain a three-dimensional seismic attribute volume of the target layer section after enhancement features;

[0147] Step 304: return to step 303.2, and another fault development direction is determined again, and steps 303.2 to 303.3 are repeated to obtain three-dimensional seismic attribute volumes of all fault development directions in the target layer section after enhancement features;

[0148] Step 305: an along-layer slice of all three-dimensional seismic attribute volumes is extracted and fused into an along-layer slice of the entire target layer section, i.e., a plan view representing fault development with high resolution of the target layer section, which is used for interpreting the fault development of the target layer section.

[0149] The specific embodiment process of the above method steps can be seen in Embodiment One and Embodiment Two, which will not be repeated here.

[0150] ​

[0151] The embodiment provides an electronic device based on Embodiment Two, which can be a mobile phone, a computer, or a tablet computer, etc., comprising a memory and a processor, wherein the memory stores a computer program, and the computer program is executed by the processor to implement the following method steps:

[0152] Step 401: obtaining three-dimensional seismic data of a target layer section;

[0153] Step 401.1: obtaining seismic data of a target layer section;

[0154] Step 401.2: calculating a structure guide volume of the seismic data;

[0155] Specifically, the amplitude and waveform information of the seismic data are analyzed by scanning a time window, a structure guide volume representing azimuth information of the seismic data is established as a structure guide volume of the seismic data;

[0156] Step 401.3: pre-processing the seismic data according to the structure guide volume of the seismic data to obtain three-dimensional seismic data of a target layer section;

[0157] Step 401.3.1: sequentially performing bad channel processing, cutting off an abnormal area in a shallow layer, and cutting processing on the obtained seismic data of a target layer section to establish a seismic sub-volume;

[0158] Step 401.3.2: removing random noise interference of the seismic sub-volume under the condition that the original geological form of the seismic sub-volume is not changed according to the structure guide volume of the seismic data, to obtain denoised seismic data;

[0159] Step 401.3.3: performing dip angle guide filtering processing on the denoised seismic data to highlight fracture detail features, to obtain filtered seismic data, to obtain three-dimensional seismic data of a target layer section;

[0160] Step 402: calculating an edge detection volume according to the three-dimensional seismic data, and performing fault enhancement processing on the edge detection volume to obtain a fault enhancement volume;

[0161] Step 402.1: calculating a structure guide volume of the three-dimensional seismic data;

[0162] Specifically, by scanning the amplitude and waveform information of the three-dimensional seismic data, a structure guide volume representing the azimuth information of the three-dimensional seismic data is established as a structure guide volume of the three-dimensional seismic data;

[0163] Step 402.2: edge detection is performed on the three-dimensional seismic data by using a layer edge detection algorithm to image the discontinuity features of the three-dimensional seismic data, and an edge detection volume is obtained, wherein the scanning window during edge detection is along the distribution direction of the stratum according to the structure guide volume of the three-dimensional seismic data;

[0164] Step 402.3: a fault enhancement processing is performed on the edge detection volume to obtain a fault enhancement volume;

[0165] Step 403: a variance operation and an ant tracking operation are sequentially performed on the fault enhancement volume, wherein during the ant tracking operation, angle control is performed according to the fault development direction to retain a fracture volume along the fault development direction, and a continuity enhancement processing is performed on the fracture volume to obtain a three-dimensional seismic attribute volume of the target layer section after enhancement features;

[0166] Step 403.1: similarity between adjacent seismic traces inside the fault enhancement volume is calculated based on an error analysis method, and a region with a similarity lower than a given threshold is enlarged to enhance its boundary features, and a variance volume of the target layer section is obtained;

[0167] Step 403.2: a fault development direction is determined according to the structure guide volume of the three-dimensional seismic data and geological understanding;

[0168] Step 403.3: a plurality of ant tracking operations are performed on the variance volume, wherein during the ant tracking operation, angle control of dip angle and azimuth angle is performed according to the fault development direction to retain a fracture volume along the fault development direction, and a continuity enhancement processing is performed on the fracture volume by analyzing curve closure of a scanning window to obtain a three-dimensional seismic attribute volume of the target layer section after enhancement features;

[0169] Step 404: return to step 403.2, and another fault development direction is determined again, and steps 403.2 to 403.3 are repeated to obtain three-dimensional seismic attribute volumes of all fault development directions in the target layer section after enhancement features;

[0170] Step 405: an along-layer slice of all three-dimensional seismic attribute volumes is extracted and fused into an along-layer slice of the entire target layer section, i.e., a plan view representing fault development with high resolution of the target layer section, which is used for interpreting the fault development of the target layer section.

[0171] The specific embodiment process of the above method steps can be seen in Embodiment One and Embodiment Two, which will not be repeated here.

[0172] It can be understood that the electronic device can further include a multimedia component, an input / output (I / O) interface, and a communication component.

[0173] The processor is configured to perform all or part of the steps of the application management method described in Embodiment One or the application management method described in Embodiment Two.

[0174] The processor can be an Application Specific Integrated Circuit (ASIC), a Digital Signal Processor (DSP), a Digital Signal Processing Device (DSPD), a Programmable Logic Device (PLD), a Field Programmable Gate Array (FPGA), a controller, a microcontroller, a microprocessor, or other electronic components, and is configured to perform the application management method described in Embodiment One or the application management method described in Embodiment Two.

[0175] The memory is configured to store various types of data, which can include, for example, instructions of any application program or method in the electronic device, and application program related data.

[0176] The memory can be implemented by any type of volatile or non-volatile storage device or a combination thereof, such as Static Random Access Memory (SRAM), Electrically Erasable Programmable Read-Only Memory (EEPROM), Erasable Programmable Read-Only Memory (EPROM), Programmable Read-Only Memory (PROM), Read-Only Memory (ROM), magnetic memory, flash memory, magnetic disk, or optical disk.

[0177] The multimedia component can include a screen and an audio component, the screen can be a touch screen, and the audio component is used to output and / or input audio signals. For example, the audio component can include a microphone for receiving external audio signals. The received audio signals can be further stored in the memory or transmitted through the communication component. The audio component also includes at least one speaker for outputting audio signals.

[0178] The I / O interface provides an interface between the processor and other interface modules, which can be a keyboard, a mouse, a button, etc. These buttons can be virtual buttons or physical buttons.

[0179] The communication component is used for wired or wireless communication between the electronic device and other devices. Wireless communication, such as Wi-Fi, Bluetooth, Near Field Communication (NFC), 2G, 3G or 4G, or a combination of one or more of them, so the corresponding communication component can include: Wi-Fi module, Bluetooth module, NFC module.

[0180] In summary, the fault enhancement method and the fault development interpretation method, the storage medium and the electronic device provided by the application are based on the fault enhancement technology, increase the optimization processing step, and perform optimization processing based on the fault enhancement volume, so that the secondary fractures can be better described compared with the coherent volume; the variance operation and the ant tracking operation are sequentially performed on the fault enhancement volume obtained after the fault enhancement processing, so that the fault image is clearer and has higher resolution; and during the ant tracking operation, the angle control and the continuity enhancement processing are performed according to the fault development direction, the three-dimensional seismic attribute volume of the target layer section after the enhanced features is obtained, the signal-to-noise ratio and the imaging quality of the seismic data are improved, the fault continuity is enhanced, the faults in the target layer section are further clearly described, and the subsequent accurate identification and evaluation of the faults are facilitated; a complete fault enhancement method with good application effect is established, the processing flow of the fault development interpretation method is improved, and the detection and identification of the small faults are more advantageous; after the seismic data of the target layer section is obtained, the structural guide volume of the seismic data is calculated, the step of calculating the structural guide volume of the seismic data is added before the preprocessing, so that the original seismic features are not changed when the seismic data is preprocessed, and the weakening of the fault information is not easily caused; according to the structural guide volume of the seismic data, the seismic data is preprocessed to obtain the three-dimensional seismic data of the target layer section, the signal-to-noise ratio of the seismic data with weak fault information is improved, and the imaging quality of the seismic data is improved; before the step of removing the random noise interference of the seismic data, the seismic data of the target layer section obtained is sequentially subjected to the bad channel processing, the cutting of the abnormal shallow layer region and the cutting processing to establish the seismic sub-volume, the larger seismic data can be processed, the subsequent denoising and filtering are facilitated, and the processing efficiency and speed of the preprocessing are improved; the optimized fault enhancement method is used for the fault development interpretation, after the three-dimensional seismic attribute volume of all the fault development directions is obtained, the along-layer slice of the entire target layer section is obtained by fusing the along-layer slices of all the three-dimensional seismic attribute volumes, the fault interpretation and extraction are facilitated, and the accuracy of the fault interpretation is improved.

[0181] It should be noted that, due to the fact that the drawings of the specification cannot be colored and repainted, some parts that are obviously different in the drawings of the application are relatively difficult to display, and if necessary, color pictures can be provided.

[0182] In the several embodiments of the application provided in the embodiments, it should be understood that the disclosed system and method can also be implemented by other manners. The system and method embodiments described above are only exemplary.

[0183] It should be noted that, as used in this document, the terms "comprises", "comprising", or any other variation thereof, are intended to cover a non-exclusive inclusion, such that a process, method, article, or apparatus that comprises a list of elements does not include only those elements but can include other elements not expressly listed or inherent to such process, method, article, or apparatus. An element proceeded by "comprises... a" does not, without more constraints, exclude the presence of additional identical elements in the process, method, article, or apparatus that comprises the element.

[0184] While the present application has been disclosed with reference to the embodiments described above, it should be understood that these have been disclosed by way of example only and that numerous modifications, additions and deletions could be made without departing from the spirit and scope of the application.

Claims

1. A fault enhancement method characterized by, The method comprises the following steps: acquiring three-dimensional seismic data of a target interval; calculating an edge detection volume according to the three-dimensional seismic data, and performing fault enhancement processing on the edge detection volume to obtain a fault enhancement volume; performing variance operation and ant tracking operation on the fault enhancement volume in sequence, wherein, when performing ant tracking operation, angle control is performed according to a fault development direction to retain a fracture volume along the fault development direction, and continuity enhancement processing is performed on the fracture volume to obtain a three-dimensional seismic attribute volume of the target interval after enhancement of features; the acquiring of the three-dimensional seismic data of the target interval comprises: acquiring seismic data of a target interval and calculating a structure guide volume of the seismic data; performing preprocessing on the seismic data according to the structure guide volume of the seismic data to acquire three-dimensional seismic data of a target interval; calculating the structure guide volume of the seismic data comprises: analyzing amplitude and waveform information of the seismic data through scanning time window, and establishing a structure guide volume representing azimuth information of the seismic data as the structure guide volume of the seismic data; calculating an edge detection volume according to the three-dimensional seismic data comprises: analyzing amplitude and waveform information of the three-dimensional seismic data through scanning time window, and establishing a structure guide volume representing azimuth information of the three-dimensional seismic data as the structure guide volume of the three-dimensional seismic data; performing edge detection on the three-dimensional seismic data by using a layer edge detection algorithm, and imaging discontinuity features of the three-dimensional seismic data to obtain an edge detection volume, wherein the scanning window during edge detection is made to follow the distribution direction of the stratum according to the structure guide volume of the three-dimensional seismic data.

2. The fault enhancement method of claim 1, wherein, performing preprocessing on the seismic data according to the structure guide volume of the seismic data to acquire three-dimensional seismic data of a target interval comprises: removing random noise interference of the seismic data to obtain denoised seismic data while keeping the original geological form of the seismic data unchanged according to the structure guide volume of the seismic data; performing dip guide filtering processing on the denoised seismic data to highlight fracture detail features, obtaining filtered seismic data, and acquiring three-dimensional seismic data of a target interval.

3. The fault enhancement method of claim 2, wherein, The seismic data comprises seismic subvolumes; before the step of removing random noise interference of the seismic data, the method further comprises: performing bad channel removal processing, cutting off abnormal shallow layer regions and cutting processing in sequence on the acquired seismic data of a target interval to establish seismic subvolumes.

4. The fault enhancement method of claim 1, wherein, performing variance operation and ant tracking operation on the fault enhancement volume in sequence, wherein, when performing ant tracking operation, angle control is performed according to a fault development direction to retain a fracture volume along the fault development direction, and continuity enhancement processing is performed on the fracture volume to obtain a three-dimensional seismic attribute volume of the target interval after enhancement of features, comprising: calculating similarity between adjacent seismic traces inside the fault enhancement volume based on an error analysis method, and enlarging a region with a similarity lower than a given threshold to enhance boundary features of the region, to obtain a variance volume of a target interval; determining a fault development direction according to the structure guide volume of the three-dimensional seismic data and geological understanding; The variance volume is subjected to several ant tracking operations, wherein, during the ant tracking operation, the angle control of the dip angle and the azimuth angle is performed according to the fault development direction to retain a fracture volume along the fault development direction, and the continuity of the fracture volume is enhanced by analyzing the curve closure of the time window to obtain a three-dimensional seismic attribute volume of the target layer section after the enhancement.

5. A fault development interpretation method characterized by, The method comprises: a three-dimensional seismic attribute volume of all fault development directions in the target layer section after the enhancement is obtained by using the fault enhancement method according to any one of claims 1 to 4; a layer slice of all three-dimensional seismic attribute volumes is extracted and fused into a layer slice of the entire target layer section, which is used for interpreting the fault development of the target layer section.

6. A storage medium having stored thereon a computer program, characterized in that The computer program can be executed by one or more processors to implement the steps of the method according to any one of claims 1 to 5.

7. An electronic device, comprising: The electronic device comprises a memory and a processor, and the memory stores a computer program which, when executed by the processor, implements the steps of the method according to any one of claims 1 to 5. The electronic device comprises a memory and a processor, and the memory stores a computer program which, when executed by the processor, implements the steps of the method according to any one of claims 1 to 5.

Images