CN114139229A - Fault enhancement method, fault development interpretation method, storage medium, and electronic apparatus - Google Patents
Fault enhancement method, fault development interpretation method, storage medium, and electronic apparatus Download PDFInfo
- Publication number
- CN114139229A CN114139229A CN202010923056.XA CN202010923056A CN114139229A CN 114139229 A CN114139229 A CN 114139229A CN 202010923056 A CN202010923056 A CN 202010923056A CN 114139229 A CN114139229 A CN 114139229A
- Authority
- CN
- China
- Prior art keywords
- fault
- seismic data
- enhancement
- target interval
- dimensional
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/27—Design optimisation, verification or simulation using machine learning, e.g. artificial intelligence, neural networks, support vector machines [SVM] or training a model
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/004—Artificial life, i.e. computing arrangements simulating life
- G06N3/006—Artificial life, i.e. computing arrangements simulating life based on simulated virtual individual or collective life forms, e.g. social simulations or particle swarm optimisation [PSO]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T15/00—3D [Three Dimensional] image rendering
- G06T15/005—General purpose rendering architectures
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Artificial Intelligence (AREA)
- Computer Graphics (AREA)
- Computer Hardware Design (AREA)
- Life Sciences & Earth Sciences (AREA)
- Molecular Biology (AREA)
- Health & Medical Sciences (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Biomedical Technology (AREA)
- Biophysics (AREA)
- Computational Linguistics (AREA)
- Data Mining & Analysis (AREA)
- General Health & Medical Sciences (AREA)
- Medical Informatics (AREA)
- Computing Systems (AREA)
- Mathematical Physics (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention relates to the technical field of oil and gas exploration, in particular to a fault enhancement method, a fault development interpretation method, a storage medium and electronic equipment, which solve the problems that the method for carrying out fault enhancement on seismic data with poor quality and weak fault information in the prior art has large calculated amount, unobvious effect of enhancing tiny faults and poor fault continuity is difficult to overcome; the method comprises the following steps: acquiring three-dimensional seismic data of a target layer section, calculating an edge detection body and performing fault enhancement treatment to obtain a fault enhancement body; carrying out variance operation and ant tracking operation on the fault reinforcement in sequence, wherein angle control is carried out according to the fault development direction, and continuity reinforcement processing is carried out on the fracture body along the fault development direction to obtain a three-dimensional seismic attribute body of the target interval after the characteristics are reinforced; the purposes of improving the signal-to-noise ratio and the imaging quality of seismic data, further clearly depicting the fault and facilitating subsequent fault identification and extraction are achieved.
Description
Technical Field
The invention relates to the technical field of oil and gas exploration, in particular to a fault enhancement method, a fault development interpretation method, a storage medium and electronic equipment.
Background
Faults are tectonic forms that develop widely in tectonic movements, and because of the presence of tectonic stresses, the underground rock formations will inevitably fracture, thereby generating faults, which are one of the most important tectonic structures of the earth's crust. In the technical field of oil and gas exploration, the fault can be used as a channel for oil and gas migration and can also be used as the boundary of a block oil and gas field, so that the distribution of the oil and gas field is effectively controlled, and therefore, the fault is accurately identified, and the fault detection method has important significance for the exploration and development of the oil and gas field.
The seismic exploration generally comprises three steps of seismic data acquisition, seismic data processing and seismic data interpretation, wherein in the stage of acquiring and processing the seismic data, due to the influence of a plurality of factors, the signal-to-noise ratio of the seismic data is low and the imaging quality is poor, so that fault information in the seismic data is fuzzy, fault detection and identification are not facilitated, and the accuracy of seismic data interpretation is influenced. Therefore, the fault in the seismic data is subjected to signal enhancement while the noise is reduced through the fault enhancement technology, so that the fault depiction is clearer, and the enhanced fault is detected and identified, so that the accuracy of seismic data interpretation is improved.
However, most of the existing fault enhancement methods enhance the seismic coherence body through a filtering algorithm, so as to enhance fault information in the seismic coherence body. When the method is used for enhancing the fault in the seismic data with low signal-to-noise ratio, poor imaging quality and weak fault information, the problems of large calculated amount, unobvious effect on enhancing the tiny fault and poor fault continuity are difficult to solve, so that the accuracy of subsequent fault identification and detection is hindered.
Therefore, the present invention has been made in view of the above problems, and provides a fault enhancement method for performing optimization processing based on fault enhancement, and a fault development interpretation method, a storage medium, and an electronic apparatus based on the fault enhancement method.
Disclosure of Invention
The invention aims to: aiming at the problems, the invention provides an optimized fault enhancement processing method, an optimized fault enhancement processing device, a storage medium and electronic equipment, wherein the optimized processing step is added on the basis of a fault enhancement technology, the optimized processing is carried out on the basis of a fault enhancement body, and the variance operation and the ant tracking operation are sequentially carried out on the fault enhancement body obtained after the fault enhancement processing, so that the problems of larger calculated amount, unobvious enhancement effect of tiny faults and poor fault continuity which are difficult to overcome in the prior art of fault enhancement on seismic data with poor quality and weak fault information are solved, the signal-to-noise ratio and the imaging quality of the seismic data are improved, the faults are further clearly delineated, and the purposes of fault identification and extraction in the follow-up process are facilitated.
The technical scheme adopted by the invention is as follows:
to achieve the above object, in a first aspect, the present invention provides a fault enhancement method, including:
acquiring three-dimensional seismic data of a target interval;
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;
and sequentially carrying out variance operation and ant tracking operation on the fault reinforcement, wherein when the ant tracking operation is carried out, angle control is carried out according to the fault development direction so as to keep the fracture body along the fault development direction, and continuity enhancement treatment is carried out on the fracture body so as to obtain the three-dimensional seismic attribute body of the target interval after the characteristics are enhanced.
According to an embodiment of the present invention, optionally, in the fault enhancement method, the acquiring three-dimensional seismic data of the target interval includes:
acquiring seismic data of a target layer section and calculating a structural guide body of the seismic data;
and preprocessing the seismic data according to the structural guide body of the seismic data to acquire the three-dimensional seismic data of the target interval.
According to an embodiment of the present invention, optionally, in the above fault enhancement method, the calculating a structure director of the seismic data includes:
and analyzing the amplitude and waveform information of the seismic data through a scanning time window, and establishing a structural guide body representing the azimuth information of the seismic data as the structural guide body of the seismic data.
According to an embodiment of the present invention, optionally, in the above fault enhancement method, the preprocessing is performed on the seismic data according to the structural director of the seismic data to obtain three-dimensional seismic data of a target interval, and the method includes:
under the condition that the original geological form of the seismic data is kept unchanged according to the structural guide body of the seismic data, removing random noise interference of the seismic data to obtain denoised seismic data;
and performing dip angle guiding filtering processing on the denoised seismic data to highlight the fracture detail characteristics to obtain the filtered seismic data so as to obtain the three-dimensional seismic data of the target interval.
According to an embodiment of the present invention, optionally, in the fault enhancement method, the seismic data includes a seismic sub-volume;
prior to the step of removing random noise interference from the seismic data, the method further comprises:
and sequentially carrying out bad track removing treatment, shallow abnormal area removing treatment and cutting treatment on the acquired seismic data of the target interval to establish a seismic sub-body.
According to an embodiment of the present invention, optionally, in the above fault enhancement method, calculating an edge detection body according to the three-dimensional seismic data includes:
analyzing the amplitude and waveform information of the three-dimensional seismic data through a scanning time window, and establishing a construction guide body representing the azimuth information of the three-dimensional seismic data as the construction guide body of the three-dimensional seismic data;
and performing edge detection on the three-dimensional seismic data by adopting an edge detection algorithm along the stratum, and imaging discontinuity characteristics of the three-dimensional seismic data to obtain an edge detection body, wherein a scanning window in the edge detection is enabled to be along the spreading direction of the stratum according to a structural guide body of the three-dimensional seismic data.
According to an embodiment of the present invention, optionally, in the above fault enhancement method, the fault enhancement body is sequentially subjected to variance operation and ant tracking operation, wherein when the ant tracking operation is performed, angle control is performed according to a fault development direction so as to retain a fracture body along the fault development direction, and continuity enhancement processing is performed on the fracture body so as to obtain a three-dimensional seismic attribute body of a target interval after feature enhancement, the method includes:
calculating the similarity between adjacent seismic channels in the fault reinforcement body based on an error analysis method, and amplifying the region with the similarity lower than a given threshold value to reinforce the boundary characteristics of the region so as to obtain a variance body of a target interval;
determining the fault development direction according to the structural guide body and geological knowledge of the three-dimensional seismic data;
and carrying out ant tracking operation on the variance body for a plurality of times, wherein during the ant tracking operation, the angle control of an inclination angle and an azimuth angle is carried out according to the fault development direction so as to keep a fracture body along the fault development direction, and the continuity enhancement treatment is carried out on the fracture body by analyzing the curve closure of a time window so as to obtain the three-dimensional seismic attribute body of the target interval with enhanced characteristics.
In a second aspect, the present invention provides a fault development interpretation method, the method comprising:
acquiring three-dimensional seismic attribute bodies of all fault development directions in the target interval after characteristic enhancement by using the fault enhancement method;
and extracting the layered slices of all the three-dimensional seismic attribute bodies, fusing the layered slices into the layered slice of the whole target interval, and explaining the fault development condition of the target interval.
In a third aspect, the present invention provides a storage medium having stored thereon a computer program executable by one or more processors to perform the steps of the method as described above.
In a fourth aspect, the invention provides an electronic device comprising a memory and a processor, the memory having stored thereon a computer program which, when executed by the processor, performs the steps of the method as described above.
Compared with the prior art, one or more embodiments in the above scheme can have the following advantages or beneficial effects:
1. according to the fault enhancement method, the fault development interpretation method, the storage medium and the electronic device, the optimization processing step is added on the basis of the fault enhancement technology, optimization processing is carried out on the basis of the fault enhancement body, and compared with a coherent body, secondary fracture can be better depicted; the variance operation and the ant tracking operation are sequentially carried out on the fault reinforcement obtained after the fault reinforcement processing is carried out, so that the fault image is clearer and the resolution is higher; during ant tracking operation, angle control and continuity enhancement processing are carried out according to the fault development direction, a three-dimensional seismic attribute body of the target interval with enhanced characteristics is obtained, the signal-to-noise ratio and the imaging quality of seismic data are improved, fault continuity is enhanced, faults in the target interval are further clearly depicted, and accurate identification and evaluation of the faults are conveniently carried out subsequently; a set of complete fault enhancement method with good application effect is established, the processing flow of the fault development interpretation method is perfected, and the method has more advantages in the detection and identification of the tiny fault.
2. According to the method, after the seismic data of the target interval are obtained, the structure guide body of the seismic data is calculated, and before preprocessing, the step of calculating the structure guide body of the seismic data is added, so that the original seismic characteristics cannot be changed when the seismic data are preprocessed, and the condition that fault information is weakened is not easy to occur.
3. According to the invention, the seismic data are preprocessed according to the structural guide body of the seismic data to obtain the three-dimensional seismic data of the target interval, so that the signal-to-noise ratio of the seismic data is improved and the imaging quality of the seismic data is improved for the seismic data with weaker fault information.
4. According to the method, before the step of removing the random noise interference of the seismic data, the obtained seismic data of the target interval are sequentially subjected to bad track removing treatment, shallow abnormal area removing and cutting treatment to establish a seismic sub-body, so that the larger seismic data can be processed, the follow-up denoising and filtering are facilitated, and the processing efficiency and speed of the preprocessing are improved.
5. In the invention, the fault enhancement method is adopted to carry out fault development interpretation, after the three-dimensional seismic attribute bodies in all fault development directions are obtained, the layered slices of all the three-dimensional seismic attribute bodies are extracted and then fused to obtain the layered slice of the whole target interval, so that fault interpretation and extraction are convenient, and the accuracy of fault interpretation is improved.
Drawings
The invention will be described in more detail hereinafter on the basis of embodiments and with reference to the accompanying drawings.
Fig. 1 is a schematic flow chart of a fault enhancement method according to an embodiment of the present invention.
Fig. 2 is a time slice of seismic data in step 101.1 of a fault enhancement method according to an embodiment of the present invention.
Fig. 3 is a time slice diagram of a structure of a guide body in step 101.2 of a fault enhancement method according to an embodiment of the present invention.
Fig. 4 is a time slice of three-dimensional seismic data in step 101.3 of a fault enhancement method according to an embodiment of the present invention.
Fig. 5 is a time slice diagram of a structure of a guide body in step 102.1 of a fault enhancement method according to an embodiment of the present invention.
Fig. 6 is a time slice diagram of the edge detection object in step 102.2 of the fault enhancement method according to the first embodiment of the present invention.
Fig. 7 is a time slice diagram of a fault enhancement member broken in step 102.3 of a fault enhancement method according to an embodiment of the present invention.
FIG. 8 is a cross-sectional view of a fault enhancement member broken at step 103.1 of a fault enhancement method according to an embodiment of the present invention;
FIG. 9 is a slice along the horizon of the three-dimensional seismic attribute volume at step 103.3 of a fault enhancement method provided in accordance with an embodiment of the present invention;
fig. 10 is a flowchart illustrating a fault development interpretation method according to a second embodiment of the present invention.
In the drawings, like parts are designated with like reference numerals, and the drawings are not drawn to scale.
Detailed Description
The following detailed description of the embodiments of the present invention will be provided with reference to the accompanying drawings and examples, so that how to apply the technical means to solve the technical problems and achieve the corresponding technical effects can be fully understood and implemented. The embodiments of the present invention and the features of the embodiments can be combined with each other without conflict, and the formed technical solutions are within the scope of the present invention.
Example one
Referring to fig. 1 to 8, the present embodiment provides a fault enhancement method applicable to an electronic device, and when the method is applied to the electronic device, the method executes steps 101 to 103:
step 101: acquiring three-dimensional seismic data of a target interval;
step 101.1: acquiring seismic data of a target interval;
in this embodiment, taking actual seismic data of a certain area as an example, acquiring seismic data including a target interval, as shown in fig. 2, a time slice of the seismic data, where the seismic data is original seismic data, and it can be seen from fig. 2 that fault information in the original seismic data is weak, fault characteristics are not obvious, and contain much random noise, which brings great trouble to detection and identification of faults, where a fault refers to a structure in which rock strata or rock masses inside a crust generate planar destruction or planar rheological zones under stress, and rock masses on both sides of the fault generate significant displacement;
step 101.2: calculating a tectonic pilot of the seismic data;
specifically, information such as amplitude and waveform of the seismic data is analyzed through a scanning time window, and a structural guide body representing the azimuth information of the seismic data is established and used as the structural guide body of the seismic data;
in this embodiment, the construction guide body representing the seismic data azimuth information is established to keep the original geological form unchanged in the subsequent preprocessing denoising and edge detection body calculation processes, and as shown in fig. 3, a time slice of the construction guide body of the seismic data is shown; by adding the step of calculating the structure guide body of the original seismic data, the original seismic characteristics are not changed when the original seismic data are preprocessed, so that the condition of weakening fault information is not easy to occur in preprocessing; as can be seen from FIG. 3, the tectonic guide body contains the azimuth information of the original seismic data, which is expressed in the obvious fault characteristic place, but is still seriously influenced by noise;
step 101.3: preprocessing the seismic data according to the structural guide body of the seismic data to obtain three-dimensional seismic data of a target interval; preprocessing the seismic data to optimize the input seismic data, eliminate the influence of interference such as noise and the like as much as possible, improve the signal-to-noise ratio of the seismic data and highlight fracture detail characteristics;
in this embodiment, the seismic data containing the target interval obtained in step 101.1 is continuously preprocessed, so that the resolution is improved, and three-dimensional seismic data with higher quality containing the target interval is obtained;
step 101.3.1: sequentially carrying out bad track removing treatment, shallow abnormal area removing treatment and cutting treatment on the obtained seismic data of the target interval to establish a seismic sub-body;
in this embodiment, the seismic traces which are acquired in step 101.1 and contain the target interval and are missing and damaged are removed, then abnormal areas caused by acquisition or other reasons in the shallow layer are cut, then cutting is performed to cut off the top or bottom area, only the seismic data containing the target interval is reserved as a seismic daughter, the processing is performed on the larger original seismic data, and the processing efficiency and speed of the preprocessing can be improved;
it should be noted that, this step may be performed according to the actual situation of the size of the original seismic data, and when this step is not performed, other steps of preprocessing the seismic data acquired in step 101.1 are directly performed subsequently;
step 101.3.2: under the condition that the original geological form of the seismic subvolume is kept unchanged according to the structural guide body of the seismic data, removing random noise interference of the seismic subvolume to obtain the denoised seismic data;
in this embodiment, the structural director of the seismic data obtained in step 101.2 and the seismic subvolume obtained in step 101.3.1 are used as input data to remove the interference such as random noise therein, so as to obtain denoised seismic data;
step 101.3.3: performing dip angle guiding filtering processing on the denoised seismic data, and highlighting fracture detail characteristics to a certain extent while retaining original seismic information to obtain filtered seismic data so as to obtain three-dimensional seismic data of a target layer section;
in this embodiment, an appropriate calculation parameter is selected to perform filtering processing on the denoised seismic data, so as to obtain three-dimensional seismic data with noise removed and reflection event transverse continuity and actual discontinuity information enhanced, as shown in fig. 4, a time slice of the three-dimensional seismic data; as can be seen from fig. 4, compared with the time slice of the seismic data shown in fig. 2, the white speck random noise in fig. 4 is significantly reduced, and the signal-to-noise ratio of the slice is significantly improved, which indicates that the method is used for preprocessing the acquired original seismic data, so that the interferences such as random noise and the like are effectively removed, the signal-to-noise ratio of the seismic data is improved, and the effect of highlighting the fracture detail characteristics is also achieved.
Step 102: 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;
step 102.1: calculating a tectonic pilot of the three-dimensional seismic data;
specifically, information such as amplitude and waveform of the three-dimensional seismic data is analyzed through a scanning time window, and a construction guide body representing azimuth information of the three-dimensional seismic data is established and used as the construction guide body of the three-dimensional seismic data;
in this embodiment, the construction guide body of the three-dimensional seismic data azimuth information obtained in the characterization step 101.3.3 is established, so that the dominant azimuth of fault development can be conveniently and visually obtained from the construction guide body in the follow-up process, and the method is suitable for obvious trunk fracture; meanwhile, the method also provides for angle control in ant tracking operation in the subsequent optimization, and is shown as a time slice of the structure guide body of the three-dimensional seismic data in fig. 5; as can be seen in fig. 5, the imaging quality corresponding to the three-dimensional seismic data of fig. 4 is improved, as is the imaging quality of fig. 5;
step 102.2: performing edge detection on the three-dimensional seismic data by adopting an edge detection algorithm along a stratum, and imaging discontinuity characteristics of the three-dimensional seismic data to obtain an edge detection body, wherein a scanning window in the edge detection is enabled to be along the spreading direction of a stratum according to a structural guide body of the three-dimensional seismic data;
in this embodiment, an edge detection algorithm along a layer is adopted to perform edge detection on preprocessed three-dimensional seismic data, and a structural guide body of the three-dimensional seismic data is added in the calculation process, so as to ensure that a scanning window is along the spreading direction of a stratum, and an edge detection body is calculated, as shown in fig. 6, a time slice of the edge detection body is obtained;
step 102.3: carrying out fault enhancement processing on the edge detection body to obtain a fault enhancement body; the effective signal is enhanced, and the fault depicting capability is enhanced;
in this embodiment, the edge detection object of step 102.2 is used as input data to perform fault enhancement processing to obtain a fault enhancement body, wherein fig. 7 shows a time slice of the fault enhancement body, and fig. 8 shows a layered slice of the fault enhancement body;
as can be seen from fig. 7 and 8, the imaging quality of the slice along the slice of the slice enhancement is improved compared with the time slice, but the continuity of the slice is still poor, in order to further enhance the useful signal of the slice, the slice is depicted more clearly, the boundary feature of the slice is more obvious, the continuity is better, and then the optimization processing is performed on the basis of the slice enhancement;
step 103: carrying out variance operation and ant tracking operation on the fault reinforcement in sequence, wherein when the ant tracking operation is carried out, angle control is carried out according to the fault development direction so as to keep a fracture body along the fault development direction, and continuity enhancement processing is carried out on the fracture body so as to obtain a three-dimensional seismic attribute body of the target interval after the characteristics are enhanced;
in this embodiment, because the fault information of 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 a certain extent, and the imaging quality is improved, the obtained fault enhancement body still cannot well depict the fault, the fault information on the slice is weak, and the continuity is poor, so the optimization processing is continued through step 103;
step 103.1: calculating the similarity between adjacent seismic channels in the fault reinforcement body based on an error analysis method, and amplifying the region with the similarity lower than a given threshold value to reinforce the boundary characteristics of the region so as to obtain a variance body of a target interval;
in the embodiment, the fault reinforcement obtained in step 102.3 is used as input data to perform variance calculation, the similarity between adjacent seismic channels in the fault reinforcement is calculated, discontinuous areas, namely areas with the similarity lower than a given threshold value, are found, and the areas with poor similarity are amplified to enhance the boundary characteristics of the areas so as to obtain a variance of a target interval; the variance calculation is based on an error analysis theory, and is characterized in that the horizontal heterogeneity of strata, lithology and the like is described by utilizing the similarity between seismic signals of adjacent channels, when faults exist underground or strata in a certain local area are discontinuously changed, the reflection characteristics of some seismic channels are different from the reflection characteristics of the seismic channels nearby the seismic channels, so that the local discontinuity of the seismic channels is caused, the information of the faults or the discontinuous change can be detected by detecting the difference degree between the seismic channels, the variance calculation generates a new variance data body by quantizing the coherence attribute of a fault reinforcement body, the irrelevance of the seismic data is highlighted and emphasized, and cracks, faults, rivers, lithology boundaries and the like are reflected; the boundary feature in this embodiment is essentially a region with poor similarity in the fault reinforcement, i.e., a discontinuous place, which is a fracture region to be found, i.e., a fault in the target interval, and is amplified, i.e., the boundary feature is reinforced;
step 103.2: determining the fault development direction according to the structural guide body and geological knowledge of the three-dimensional seismic data;
in the embodiment, a fault development direction is determined according to the approximate fault development direction represented by the structure guide body of the three-dimensional seismic data and certain geological knowledge;
step 103.3: carrying out ant tracking operation on the variance cube for a plurality of times, so that the sectional image is clearer and the resolution ratio is higher; during ant tracking operation, angle control of an inclination angle and an azimuth angle is carried out according to the fault development direction so as to keep a fracture body along the fault development direction, and continuity enhancement processing is carried out on the fracture body by analyzing curve closure of a time window so as to obtain a three-dimensional seismic attribute body of a target interval with enhanced characteristics;
in this embodiment, the variance vector is subjected to ant tracking operation for three times, and when ant tracking operation is performed for the third time, angle control of an inclination angle and an azimuth angle is performed according to the analysis range of the fault development direction scanning time window determined in step 103.2, so that incoherence caused by changes of a stratum inclination angle and an azimuth angle is reduced, discontinuity caused by fault layers is highlighted, useful signals are enhanced, only a fracture vector in the fault development direction is retained, the purpose of running for three times is to highlight fault characteristics to the maximum extent and increase continuity of fault characteristics, and the number of ant tracking operation times can be selected according to actual conditions;
the continuity enhancement processing is carried out on the fracture body by analyzing the curve closure of the time window, the better continuity of the useful signal acquisition is enhanced to the maximum extent, the three-dimensional seismic attribute body of the target interval with enhanced characteristics is obtained, as shown in fig. 9, the slice along the layer of the three-dimensional seismic attribute body is shown, as can be seen from fig. 9, compared with the slice along the layer of the fault enhancement body of fig. 8, the useful signal of the slice along the layer of the three-dimensional seismic attribute body is enhanced to some extent, the continuity is better, and the imaging quality is also improved.
According to the fault enhancement method provided by the embodiment, the optimization processing step is added on the basis of the fault enhancement technology, and the optimization processing is carried out on the basis of the fault enhancement body, so that the secondary fracture can be better depicted compared with a coherent body; the variance operation and the ant tracking operation are sequentially carried out on the fault reinforcement obtained after the fault reinforcement processing is carried out, so that the fault image is clearer and the resolution is higher; during ant tracking operation, angle control and continuity enhancement processing are carried out according to the fault development direction, a three-dimensional seismic attribute body of the target interval with enhanced characteristics is obtained, the signal-to-noise ratio and the imaging quality of seismic data are improved, fault continuity is enhanced, faults in the target interval are further clearly depicted, and accurate identification and evaluation of the faults are conveniently carried out subsequently; a set of complete fault enhancement method with good application effect is established, the processing flow of the fault development interpretation method is perfected, and the fault enhancement method has more advantages in the detection and identification of the tiny fault; moreover, by adding the step of calculating the structure guide body of the original seismic data, the original seismic characteristics are not changed when the original seismic data are preprocessed, so that the condition of weakening fault information is not easy to occur in preprocessing; the method solves the problems that the method for enhancing the fault of the seismic data with poor quality and weak fault information in the prior art has large calculated amount, unobvious effect of enhancing the tiny fault and poor fault continuity and is difficult to overcome.
Example two
Referring to fig. 2 to 10, the present embodiment further provides a fault development interpretation method applicable to an electronic device based on the first embodiment, and when the method is applied to the electronic device, the steps 201 to 204 are executed:
step 201: acquiring three-dimensional seismic data of a target interval;
step 201.1: acquiring seismic data of a target interval;
in this embodiment, taking the actual seismic data of a certain area as an example, the seismic data including the target interval is obtained, and as shown in fig. 2, the seismic data is a time slice of the seismic data, and the seismic data is the original seismic data, and it can be seen from fig. 2 that the fault information in the original seismic data is weak, the fault characteristics are not obvious in representation, and contain more random noise, which brings great trouble to the detection and identification of the fault;
step 201.2: calculating a tectonic pilot of the seismic data;
specifically, amplitude and waveform information of the seismic data are analyzed through a scanning time window, and a structural guide body representing azimuth information of the seismic data is established and used as the structural guide body of the seismic data;
in this embodiment, the construction guide body representing the seismic data azimuth information is established to keep the original geological form unchanged in the subsequent preprocessing denoising and edge detection body calculation processes, and as shown in fig. 3, a time slice of the construction guide body of the seismic data is shown; by adding the step of calculating the structure guide body of the original seismic data, the original seismic characteristics are not changed when the original seismic data are preprocessed, so that the condition of weakening fault information is not easy to occur in preprocessing; as can be seen from FIG. 3, the tectonic guide body contains the azimuth information of the original seismic data, which is expressed in the obvious fault characteristic place, but is still seriously influenced by noise;
step 201.3: preprocessing the seismic data according to the structural guide body of the seismic data to obtain three-dimensional seismic data of a target interval;
in this embodiment, the seismic data including the target interval acquired in step 201.1 is continuously preprocessed, so that the resolution is improved, and three-dimensional seismic data including the target interval and having higher quality is obtained;
step 201.3.1: sequentially carrying out bad track removing treatment, shallow abnormal area removing treatment and cutting treatment on the obtained seismic data of the target interval to establish a seismic sub-body;
in this embodiment, the seismic traces which are acquired in step 201.1 and contain the target interval and are missing and damaged are removed, then abnormal areas caused by acquisition or other reasons in the shallow layer are cut, then cutting is performed to cut off the top or bottom area, only the seismic data containing the target interval is reserved as a seismic daughter, the processing is performed on the larger original seismic data, and the processing efficiency and speed of the preprocessing can be improved;
step 201.3.2: under the condition that the original geological form of the seismic subvolume is kept unchanged according to the structural guide body of the seismic data, removing random noise interference of the seismic subvolume to obtain the denoised seismic data;
in this embodiment, the structural director of the seismic data obtained in step 201.2 and the seismic subvolume obtained in step 201.3.1 are used as input data, and interferences such as random noise and the like in the input data are removed to obtain denoised seismic data;
step 201.3.3: performing dip angle guiding filtering processing on the denoised seismic data, and highlighting fracture detail characteristics to a certain extent while retaining original seismic information to obtain filtered seismic data so as to obtain three-dimensional seismic data of a target layer section;
in this embodiment, an appropriate calculation parameter is selected to perform filtering processing on the denoised seismic data, so as to obtain three-dimensional seismic data with noise removed and reflection event transverse continuity and actual discontinuity information enhanced, as shown in fig. 4, a time slice of the three-dimensional seismic data; as can be seen from fig. 4, compared to the time slice of the seismic data of fig. 2, the random noise of white spots in fig. 4 is significantly reduced, and the signal-to-noise ratio of the slice is significantly improved;
step 202: 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;
step 202.1: calculating a tectonic pilot of the three-dimensional seismic data;
specifically, amplitude and waveform information of the three-dimensional seismic data are analyzed through a scanning time window, and a construction guide body representing azimuth information of the three-dimensional seismic data is established and used as the construction guide body of the three-dimensional seismic data;
in this embodiment, the construction guide body of the three-dimensional seismic data azimuth information obtained in the characterization step 201.3.3 is established, so that the dominant azimuth of fault development can be conveniently and visually obtained from the construction guide body in the follow-up process, and the method is suitable for obvious trunk fracture; meanwhile, the method also provides for angle control in ant tracking operation in the subsequent optimization, and is shown as a time slice of the structure guide body of the three-dimensional seismic data in fig. 5; as can be seen in fig. 5, the imaging quality corresponding to the three-dimensional seismic data of fig. 4 is improved, as is the imaging quality of fig. 5;
step 202.2: performing edge detection on the three-dimensional seismic data by adopting an edge detection algorithm along a stratum, and imaging discontinuity characteristics of the three-dimensional seismic data to obtain an edge detection body, wherein a scanning window in the edge detection is enabled to be along the spreading direction of a stratum according to a structural guide body of the three-dimensional seismic data;
in this embodiment, an edge detection algorithm along a layer is adopted to perform edge detection on preprocessed three-dimensional seismic data, and a structural guide body of the three-dimensional seismic data is added in the calculation process, so as to ensure that a scanning window is along the spreading direction of a stratum, and an edge detection body is calculated, as shown in fig. 6, a time slice of the edge detection body is obtained;
step 202.3: carrying out fault enhancement processing on the edge detection body to obtain a fault enhancement body;
in this embodiment, the edge detection object of step 202.2 is used as input data to perform the slice enhancement processing to obtain a slice enhancement body, as shown in fig. 7, a time slice of the slice enhancement body, and as shown in fig. 8, a slice along the slice of the slice enhancement body;
step 203: carrying out variance operation and ant tracking operation on the fault reinforcement in sequence, wherein when the ant tracking operation is carried out, angle control is carried out according to the fault development direction so as to keep a fracture body along the fault development direction, and continuity enhancement processing is carried out on the fracture body so as to obtain a three-dimensional seismic attribute body of the target interval after the characteristics are enhanced;
step 203.1: calculating the similarity between adjacent seismic channels in the fault reinforcement body based on an error analysis method, and amplifying the region with the similarity lower than a given threshold value to reinforce the boundary characteristics of the region so as to obtain a variance body of a target interval;
in the embodiment, the fault reinforcement obtained in step 202.3 is used as input data to perform variance calculation, the similarity between adjacent seismic channels in the fault reinforcement is calculated, and the region with the similarity lower than a given threshold is amplified to enhance the boundary characteristics of the region so as to obtain a variance of a target interval;
step 203.2: determining the fault development direction according to the structural guide body and geological knowledge of the three-dimensional seismic data;
in the embodiment, a fault development direction is artificially determined according to the approximate fault development direction represented by the structure guide body of the three-dimensional seismic data and certain geological knowledge;
step 203.3: carrying out ant tracking operation on the variance body for a plurality of times, wherein during the ant tracking operation, angle control of an inclination angle and an azimuth angle is carried out according to the fault development direction so as to keep a fracture body along the fault development direction, and carrying out continuity enhancement treatment on the fracture body by analyzing curve closure of a time window so as to obtain a three-dimensional seismic attribute body of the target interval with enhanced characteristics;
in this embodiment, the variance body is subjected to ant tracking operation for three times, and when ant tracking operation is performed for the third time, angle control of an inclination angle and an azimuth angle is performed according to the analysis range of the fault development direction scanning time window determined in step 203.2, so that incoherence caused by changes of a stratum inclination angle and an azimuth angle is reduced, discontinuity caused by fault formation is highlighted, useful signals are enhanced, and only a fracture body in the fault development direction is reserved;
performing continuity enhancement treatment on the fracture body by analyzing curve closure of a time window, enhancing useful signals to the maximum extent to obtain better continuity, and obtaining a three-dimensional seismic attribute body of a target interval with enhanced characteristics, wherein as shown in fig. 9, a slab-along slice of the three-dimensional seismic attribute body is obtained, and as can be seen from fig. 9, compared with the slab-along slice of the fault enhancement body of fig. 8, the slab-along slice useful signals of the three-dimensional seismic attribute body are enhanced, the continuity is better, and the imaging quality is also improved;
step 204: acquiring three-dimensional seismic attribute bodies of all fault development directions in the target layer section after the characteristics are enhanced;
in this embodiment, returning to step 203.2, manually re-determining another fault development direction according to the approximate fault development direction represented by the structure guide body of the three-dimensional seismic data and combining with certain geological knowledge, and obtaining a three-dimensional seismic attribute body of the target interval after the enhancement feature of the other fault development direction;
repeating the steps 203.2 to 203.3 to obtain three-dimensional seismic attribute bodies in all fault development directions;
step 205: and extracting the layered slices of all the three-dimensional seismic attribute bodies, fusing the layered slices into the layered slice of the whole target interval, and explaining the fault development condition of the target interval so as to improve the accuracy of fault explanation.
In this embodiment, the fusion is a slab-wise slice of the entire target interval, that is, a planar diagram representing fault development of which the target interval has higher resolution, and the fault development condition of the target interval can be explained according to the diagram.
According to the fault development interpretation method provided by the embodiment, the fault enhancement method after the optimization processing step is added for fault development interpretation on the basis of the fault enhancement technology, firstly, fault enhancement processing and optimization processing are carried out on original seismic data, and then, after three-dimensional seismic attribute bodies in all fault development directions are obtained, the layered slices of all the three-dimensional seismic attribute bodies are extracted and fused to obtain the layered slice of the whole target interval, so that fault interpretation and extraction are facilitated, and the accuracy of fault interpretation is improved.
EXAMPLE III
This embodiment provides a computer readable storage medium, such as a flash memory, a hard disk, a multimedia card, a card type memory (e.g., 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 store, etc., on which a computer program is stored, which when executed by a processor can implement the following method steps:
step 301: acquiring three-dimensional seismic data of a target interval;
step 301.1: acquiring seismic data of a target interval;
step 301.2: calculating a tectonic pilot of the seismic data;
specifically, amplitude and waveform information of the seismic data are analyzed through a scanning time window, and a structural guide body representing azimuth information of the seismic data is established and used as the structural guide body of the seismic data;
step 301.3: preprocessing the seismic data according to the structural guide body of the seismic data to obtain three-dimensional seismic data of a target interval;
step 301.3.1: sequentially carrying out bad track removing treatment, shallow abnormal area removing treatment and cutting treatment on the obtained seismic data of the target interval to establish a seismic sub-body;
step 301.3.2: under the condition that the original geological form of the seismic subvolume is kept unchanged according to the structural guide body of the seismic data, removing random noise interference of the seismic subvolume to obtain the denoised seismic data;
step 301.3.3: performing dip angle guiding filtering processing on the denoised seismic data to highlight fracture detail characteristics to obtain filtered seismic data so as to obtain three-dimensional seismic data of a target interval;
step 302: 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;
step 302.1: calculating a tectonic pilot of the three-dimensional seismic data;
specifically, amplitude and waveform information of the three-dimensional seismic data are analyzed through a scanning time window, and a construction guide body representing azimuth information of the three-dimensional seismic data is established and used as the construction guide body of the three-dimensional seismic data;
step 302.2: performing edge detection on the three-dimensional seismic data by adopting an edge detection algorithm along a stratum, and imaging discontinuity characteristics of the three-dimensional seismic data to obtain an edge detection body, wherein a scanning window in the edge detection is enabled to be along the spreading direction of a stratum according to a structural guide body of the three-dimensional seismic data;
step 302.3: carrying out fault enhancement processing on the edge detection body to obtain a fault enhancement body;
step 303: carrying out variance operation and ant tracking operation on the fault reinforcement in sequence, wherein when the ant tracking operation is carried out, angle control is carried out according to the fault development direction so as to keep a fracture body along the fault development direction, and continuity enhancement processing is carried out on the fracture body so as to obtain a three-dimensional seismic attribute body of the target interval after the characteristics are enhanced;
step 303.1: calculating the similarity between adjacent seismic channels in the fault reinforcement body based on an error analysis method, and amplifying the region with the similarity lower than a given threshold value to reinforce the boundary characteristics of the region so as to obtain a variance body of a target interval;
step 303.2: determining the fault development direction according to the structural guide body and geological knowledge of the three-dimensional seismic data;
step 303.3: carrying out ant tracking operation on the variance body for a plurality of times, wherein during the ant tracking operation, angle control of an inclination angle and an azimuth angle is carried out according to the fault development direction so as to keep a fracture body along the fault development direction, and carrying out continuity enhancement treatment on the fracture body by analyzing curve closure of a time window so as to obtain a three-dimensional seismic attribute body of the target interval with enhanced characteristics;
step 304: returning to the step 303.2, re-determining another fault development direction, repeating the steps 303.2 to 303.3, and obtaining three-dimensional seismic attribute bodies of all fault development directions in the target layer section after the characteristics are enhanced;
step 305: extracting the layered slices of all the three-dimensional seismic attribute bodies, and fusing the layered slices into the layered slices of the whole target interval, namely a planar diagram representing fault development of the target interval with higher resolution, wherein the planar diagram is used for explaining the fault development condition of the target interval.
The specific embodiment process of the above method steps can be referred to in the first embodiment and the second embodiment, and the detailed description of the embodiment is not repeated herein.
Example four
The present embodiment provides an electronic device on the basis of the second embodiment, where the electronic device may be a mobile phone, a computer, a tablet computer, or the like, and includes a memory and a processor, where the memory stores a computer program, and when the computer program is executed by the processor, the method includes:
step 401: acquiring three-dimensional seismic data of a target interval;
step 401.1: acquiring seismic data of a target interval;
step 401.2: calculating a tectonic pilot of the seismic data;
specifically, amplitude and waveform information of the seismic data are analyzed through a scanning time window, and a structural guide body representing azimuth information of the seismic data is established and used as the structural guide body of the seismic data;
step 401.3: preprocessing the seismic data according to the structural guide body of the seismic data to obtain three-dimensional seismic data of a target interval;
step 401.3.1: sequentially carrying out bad track removing treatment, shallow abnormal area removing treatment and cutting treatment on the obtained seismic data of the target interval to establish a seismic sub-body;
step 401.3.2: under the condition that the original geological form of the seismic subvolume is kept unchanged according to the structural guide body of the seismic data, removing random noise interference of the seismic subvolume to obtain the denoised seismic data;
step 401.3.3: performing dip angle guiding filtering processing on the denoised seismic data to highlight fracture detail characteristics to obtain filtered seismic data so as to obtain three-dimensional seismic data of a target interval;
step 402: 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;
step 402.1: calculating a tectonic pilot of the three-dimensional seismic data;
specifically, amplitude and waveform information of the three-dimensional seismic data are analyzed through a scanning time window, and a construction guide body representing azimuth information of the three-dimensional seismic data is established and used as the construction guide body of the three-dimensional seismic data;
step 402.2: performing edge detection on the three-dimensional seismic data by adopting an edge detection algorithm along a stratum, and imaging discontinuity characteristics of the three-dimensional seismic data to obtain an edge detection body, wherein a scanning window in the edge detection is enabled to be along the spreading direction of a stratum according to a structural guide body of the three-dimensional seismic data;
step 402.3: carrying out fault enhancement processing on the edge detection body to obtain a fault enhancement body;
step 403: carrying out variance operation and ant tracking operation on the fault reinforcement in sequence, wherein when the ant tracking operation is carried out, angle control is carried out according to the fault development direction so as to keep a fracture body along the fault development direction, and continuity enhancement processing is carried out on the fracture body so as to obtain a three-dimensional seismic attribute body of the target interval after the characteristics are enhanced;
step 403.1: calculating the similarity between adjacent seismic channels in the fault reinforcement body based on an error analysis method, and amplifying the region with the similarity lower than a given threshold value to reinforce the boundary characteristics of the region so as to obtain a variance body of a target interval;
step 403.2: determining the fault development direction according to the structural guide body and geological knowledge of the three-dimensional seismic data;
step 403.3: carrying out ant tracking operation on the variance body for a plurality of times, wherein during the ant tracking operation, angle control of an inclination angle and an azimuth angle is carried out according to the fault development direction so as to keep a fracture body along the fault development direction, and carrying out continuity enhancement treatment on the fracture body by analyzing curve closure of a time window so as to obtain a three-dimensional seismic attribute body of the target interval with enhanced characteristics;
step 404: returning to the step 403.2, re-determining another fault development direction, repeating the steps 403.2 to 403.3, and acquiring three-dimensional seismic attribute bodies of all fault development directions in the target interval after the characteristics are enhanced;
step 405: extracting the layered slices of all the three-dimensional seismic attribute bodies, and fusing the layered slices into the layered slices of the whole target interval, namely a planar diagram representing fault development of the target interval with higher resolution, wherein the planar diagram is used for explaining the fault development condition of the target interval.
The specific embodiment process of the above method steps can be referred to in the first embodiment and the second embodiment, and the detailed description of the embodiment is not repeated herein.
It is understood that the electronic device may also include multimedia components, input/output (I/O) interfaces, and communication components.
The processor is configured to execute all or part of the steps of the application management method according to the first embodiment or the application management method according to the second embodiment.
The Processor may 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 execute the Application management method according to the first embodiment or the Application management method according to the second embodiment.
The memory is used to store various types of data, which may include, for example, instructions for any application or method in the electronic device, as well as application-related data.
The Memory may be implemented by any type of volatile or non-volatile Memory device or 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.
The multimedia components may include a screen, which may be a touch screen, and an audio component for outputting and/or inputting audio signals. For example, the audio component may include a microphone for receiving external audio signals. The received audio signal may further be stored in a memory or transmitted through a communication component. The audio assembly also includes at least one speaker for outputting audio signals.
The I/O interface provides an interface between the processor and other interface modules, such as a keyboard, a mouse, buttons, etc. These buttons may be virtual buttons or physical buttons.
The communication component is used for carrying out wired or wireless communication between the electronic equipment and other equipment. 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 that the corresponding Communication component may include: Wi-Fi module, bluetooth module, NFC module.
In summary, according to the fault enhancement method, the fault development interpretation method, the storage medium and the electronic device provided by the invention, the optimization processing step is added on the basis of the fault enhancement technology, and the optimization processing is carried out on the basis of the fault enhancement body, so that the secondary fracture can be better depicted compared with a coherent body; the variance operation and the ant tracking operation are sequentially carried out on the fault reinforcement obtained after the fault reinforcement processing is carried out, so that the fault image is clearer and the resolution is higher; during ant tracking operation, angle control and continuity enhancement processing are carried out according to the fault development direction, a three-dimensional seismic attribute body of the target interval with enhanced characteristics is obtained, the signal-to-noise ratio and the imaging quality of seismic data are improved, fault continuity is enhanced, faults in the target interval are further clearly depicted, and accurate identification and evaluation of the faults are conveniently carried out subsequently; a set of complete fault enhancement method with good application effect is established, the processing flow of the fault development interpretation method is perfected, and the fault enhancement method has more advantages in the detection and identification of the tiny fault; according to the method, after seismic data of a target interval is obtained, the structure guide body of the seismic data is calculated, and before preprocessing, the step of calculating the structure guide body of the seismic data is added, so that the original seismic characteristics cannot be changed when the seismic data are preprocessed, and the condition that fault information is weakened is not easy to occur; according to the invention, the seismic data are preprocessed according to the structural guide body of the seismic data to obtain the three-dimensional seismic data of the target interval, so that for the seismic data with weak fault information, the signal-to-noise ratio of the seismic data 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 method sequentially carries out bad channel removing treatment, shallow abnormal area removing and cutting treatment on the seismic data of the obtained target interval to establish a seismic sub-body, can treat the larger seismic data, is convenient for subsequent denoising and filtering, and improves the treatment efficiency and speed of pretreatment; according to the invention, the optimized fault enhancement method is adopted to carry out fault development interpretation, after the three-dimensional seismic attribute bodies in all fault development directions are obtained, the layered slices of all the three-dimensional seismic attribute bodies are extracted and then fused to obtain the layered slice of the whole target interval, so that fault interpretation and extraction are facilitated, and the accuracy of fault interpretation is improved.
It should be noted that, since the drawings in the specification should not be colored or modified, it is difficult to display the parts of the drawings in the present invention where the parts are clearly distinguished from each other, and if necessary, a color picture can be provided.
In the embodiments provided in the present invention, it should be understood that the disclosed system and method can be implemented in other ways. The system and method embodiments described above are exemplary only.
It should be noted that, 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 may include other elements not expressly listed or inherent to such process, method, article, or apparatus. Without further limitation, an element defined by the phrase "comprising an … …" does not exclude the presence of other identical elements in a process, method, article, or apparatus that comprises the element.
Although the embodiments of the present invention have been described above, the above descriptions are only for the convenience of understanding the present invention, and are not intended to limit the present invention. It will be understood by those skilled in the art that various changes in form and details may be made therein without departing from the spirit and scope of the invention as defined by the appended claims.
Claims (10)
1. A fault enhancement method, comprising the steps of:
acquiring three-dimensional seismic data of a target interval;
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;
and sequentially carrying out variance operation and ant tracking operation on the fault reinforcement, wherein when the ant tracking operation is carried out, angle control is carried out according to the fault development direction so as to keep the fracture body along the fault development direction, and continuity enhancement treatment is carried out on the fracture body so as to obtain the three-dimensional seismic attribute body of the target interval after the characteristics are enhanced.
2. The fault enhancement method of claim 1, wherein the acquiring three-dimensional seismic data for a target interval comprises:
acquiring seismic data of a target layer section and calculating a structural guide body of the seismic data;
and preprocessing the seismic data according to the structural guide body of the seismic data to acquire the three-dimensional seismic data of the target interval.
3. The fault enhancement method of claim 2, wherein computing a tectonic guide for the seismic data comprises:
and analyzing the amplitude and waveform information of the seismic data through a scanning time window, and establishing a structural guide body representing the azimuth information of the seismic data as the structural guide body of the seismic data.
4. The fault enhancement method of claim 2, wherein preprocessing the seismic data to obtain three-dimensional seismic data for a target interval based on a tectonic guide body of the seismic data comprises:
under the condition that the original geological form of the seismic data is kept unchanged according to the structural guide body of the seismic data, removing random noise interference of the seismic data to obtain denoised seismic data;
and performing dip angle guiding filtering processing on the denoised seismic data to highlight the fracture detail characteristics to obtain the filtered seismic data so as to obtain the three-dimensional seismic data of the target interval.
5. The fault enhancement method of claim 4, wherein the seismic data comprises seismic sub-volumes;
prior to the step of removing random noise interference from the seismic data, the method further comprises:
and sequentially carrying out bad track removing treatment, shallow abnormal area removing treatment and cutting treatment on the acquired seismic data of the target interval to establish a seismic sub-body.
6. The fault enhancement method of claim 1, wherein computing an edge detection volume from the three-dimensional seismic data comprises:
analyzing the amplitude and waveform information of the three-dimensional seismic data through a scanning time window, and establishing a construction guide body representing the azimuth information of the three-dimensional seismic data as the construction guide body of the three-dimensional seismic data;
and performing edge detection on the three-dimensional seismic data by adopting an edge detection algorithm along the stratum, and imaging discontinuity characteristics of the three-dimensional seismic data to obtain an edge detection body, wherein a scanning window in the edge detection is enabled to be along the spreading direction of the stratum according to a structural guide body of the three-dimensional seismic data.
7. The fault enhancement method according to claim 1, wherein the fault enhancement body is subjected to variance operation and ant tracking operation in sequence, wherein during ant tracking operation, angle control is performed according to a fault development direction so as to retain a fracture body along the fault development direction, and continuity enhancement processing is performed on the fracture body so as to obtain a three-dimensional seismic attribute body of a target interval after feature enhancement, and the method comprises the following steps:
calculating the similarity between adjacent seismic channels in the fault reinforcement body based on an error analysis method, and amplifying the region with the similarity lower than a given threshold value to reinforce the boundary characteristics of the region so as to obtain a variance body of a target interval;
determining the fault development direction according to the structural guide body and geological knowledge of the three-dimensional seismic data;
and carrying out ant tracking operation on the variance body for a plurality of times, wherein during the ant tracking operation, the angle control of an inclination angle and an azimuth angle is carried out according to the fault development direction so as to keep a fracture body along the fault development direction, and the continuity enhancement treatment is carried out on the fracture body by analyzing the curve closure of a time window so as to obtain the three-dimensional seismic attribute body of the target interval with enhanced characteristics.
8. A fault development interpretation method, comprising:
acquiring three-dimensional seismic attribute bodies of all fault development directions in the target interval after characteristic enhancement by using the fault enhancement method according to any one of claims 1 to 7;
and extracting the layered slices of all the three-dimensional seismic attribute bodies, fusing the layered slices into the layered slice of the whole target interval, and explaining the fault development condition of the target interval.
9. A storage medium having stored thereon a computer program executable by one or more processors to perform the steps of a method as claimed in any one of claims 1 to 8.
10. An electronic device, characterized in that the electronic device comprises a memory and a processor, the memory having stored thereon a computer program which, when executed by the processor, carries out the steps of the method according to any one of claims 1 to 8.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010923056.XA CN114139229A (en) | 2020-09-04 | 2020-09-04 | Fault enhancement method, fault development interpretation method, storage medium, and electronic apparatus |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010923056.XA CN114139229A (en) | 2020-09-04 | 2020-09-04 | Fault enhancement method, fault development interpretation method, storage medium, and electronic apparatus |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN114139229A true CN114139229A (en) | 2022-03-04 |
Family
ID=80438347
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010923056.XA Pending CN114139229A (en) | 2020-09-04 | 2020-09-04 | Fault enhancement method, fault development interpretation method, storage medium, and electronic apparatus |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114139229A (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116681840A (en) * | 2023-08-03 | 2023-09-01 | 中国科学院长春光学精密机械与物理研究所 | OSMO optimization improvement method for fault reconstruction and storage medium |
-
2020
- 2020-09-04 CN CN202010923056.XA patent/CN114139229A/en active Pending
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116681840A (en) * | 2023-08-03 | 2023-09-01 | 中国科学院长春光学精密机械与物理研究所 | OSMO optimization improvement method for fault reconstruction and storage medium |
| CN116681840B (en) * | 2023-08-03 | 2023-10-31 | 中国科学院长春光学精密机械与物理研究所 | OSMO optimization improvement method for fault reconstruction and storage medium |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US7844402B2 (en) | System and method for fault identification | |
| Cohen et al. | Detection and extraction of fault surfaces in 3D seismic data | |
| US10025971B2 (en) | Fracability measurement method and system | |
| Assous et al. | Automated detection of planar geologic features in borehole images | |
| CN102016642B (en) | Method and system for suppression of residual water bottom energy in surface seismic data | |
| US10705234B2 (en) | Line and edge detection and enhancement | |
| CN105431750A (en) | Fracability measurement method and system | |
| EP2926171B1 (en) | System and method for seismic velocity anomaly analysis | |
| EP3211594A1 (en) | Seismic modeling system providing seismic survey data inpainting based upon suspect region boundary comparisons and related methods | |
| CN114139229A (en) | Fault enhancement method, fault development interpretation method, storage medium, and electronic apparatus | |
| KR20150088816A (en) | System and method for velocity anomaly analysis | |
| Li et al. | Fracture extraction from FMI based on multiscale mathematical morphology | |
| CN112180447A (en) | Method and system for eliminating strong reflection shielding of reservoir | |
| CN113093274B (en) | Method, device, terminal and storage medium for identifying low-order faults | |
| CN113640876A (en) | Method for finely identifying trapping column by using chaotic body attribute | |
| Faiq Adenan et al. | Advancing subsurface analysis: Integrating computer vision and deep learning for the near real-time interpretation of borehole image logs in the Illinois Basin-Decatur Project | |
| Botter et al. | Seismic attribute analysis of a fault zone in the Thebe field, Northwest shelf, Australia | |
| CN113325474B (en) | Method for discriminating biological reef | |
| Sinha* et al. | Identification and Quantification of Parasequences Using Expectation Maximization Filter: Defining Well Log Attributes for Reservoir Characterization | |
| CN114415230B (en) | Linear fracture extraction method and device | |
| CN117148426A (en) | Seismic data feature boundary extraction method in image separation mode | |
| Montiel-Zafra et al. | Detection and classification of internal defects in limestone blocks based on a deconvolution technique with SI-PLCA applied to GPR signals | |
| CN113740914A (en) | Pre-stack depth migration method and device, storage medium and electronic equipment | |
| CN117828257A (en) | Fracture fine slice calculation method, fracture fine slice calculation device, computer equipment and storage medium | |
| Yong et al. | The most homogeneous dip-scanning method using edge-preserving smoothing for seismic noise attenuation |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination |