CN114185091B - Ant body crack tracking method and device based on frequency spectrum decomposition and electronic equipment - Google Patents
Ant body crack tracking method and device based on frequency spectrum decomposition and electronic equipment Download PDFInfo
- Publication number
- CN114185091B CN114185091B CN202010968413.4A CN202010968413A CN114185091B CN 114185091 B CN114185091 B CN 114185091B CN 202010968413 A CN202010968413 A CN 202010968413A CN 114185091 B CN114185091 B CN 114185091B
- Authority
- CN
- China
- Prior art keywords
- frequency
- band
- seismic data
- ant
- crack
- 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.)
- Active
Links
- 238000000354 decomposition reaction Methods 0.000 title claims abstract description 91
- 238000001228 spectrum Methods 0.000 title claims abstract description 67
- 238000000034 method Methods 0.000 title claims abstract description 56
- 238000004364 calculation method Methods 0.000 claims abstract description 16
- 230000003595 spectral effect Effects 0.000 claims description 24
- 238000004590 computer program Methods 0.000 claims description 5
- 238000005516 engineering process Methods 0.000 description 14
- 238000011160 research Methods 0.000 description 10
- 238000011161 development Methods 0.000 description 7
- 238000003384 imaging method Methods 0.000 description 7
- 238000004458 analytical method Methods 0.000 description 6
- 230000000694 effects Effects 0.000 description 6
- BVKZGUZCCUSVTD-UHFFFAOYSA-L Carbonate Chemical compound [O-]C([O-])=O BVKZGUZCCUSVTD-UHFFFAOYSA-L 0.000 description 4
- 241000257303 Hymenoptera Species 0.000 description 4
- 230000008901 benefit Effects 0.000 description 4
- 238000004422 calculation algorithm Methods 0.000 description 4
- 230000001427 coherent effect Effects 0.000 description 4
- 238000001514 detection method Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 4
- 238000012545 processing Methods 0.000 description 4
- 239000003016 pheromone Substances 0.000 description 3
- 230000001629 suppression Effects 0.000 description 3
- 230000009466 transformation Effects 0.000 description 3
- 241000220317 Rosa Species 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 230000000007 visual effect Effects 0.000 description 2
- 230000002776 aggregation Effects 0.000 description 1
- 238000004220 aggregation Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 230000037237 body shape Effects 0.000 description 1
- 238000009933 burial Methods 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 230000002349 favourable effect Effects 0.000 description 1
- 230000002431 foraging effect Effects 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 235000001968 nicotinic acid Nutrition 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000005325 percolation Methods 0.000 description 1
- 230000000704 physical effect Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/30—Analysis
- G01V1/301—Analysis for determining seismic cross-sections or geostructures
- G01V1/302—Analysis for determining seismic cross-sections or geostructures in 3D data cubes
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/30—Analysis
- G01V1/307—Analysis for determining seismic attributes, e.g. amplitude, instantaneous phase or frequency, reflection strength or polarity
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/63—Seismic attributes, e.g. amplitude, polarity, instant phase
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/64—Geostructures, e.g. in 3D data cubes
- G01V2210/642—Faults
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention discloses an ant body crack tracking method and device based on spectrum decomposition and electronic equipment, wherein the method comprises the following steps: performing spectrum decomposition on the full-band post-stack seismic data to obtain a seismic tuning body; determining a main frequency band of the effective signal of the seismic data based on the seismic tuning body; performing spectrum decomposition on the seismic data in the main frequency band to obtain a discrete single-frequency energy body; and carrying out ant body tracking calculation on the discrete single-frequency energy bodies of the low frequency band and the discrete single-frequency energy bodies of the high frequency band respectively to obtain a low-frequency-band ant body crack and a high-frequency-band ant body crack. The ant body crack tracking method based on frequency spectrum decomposition eliminates mutual interference of different frequency components in a time domain, so that an interpretation result higher than the traditional resolution is obtained, fine cracks and crack boundaries can be clearly identified, meanwhile, the influence of man-made subjectivity in broken layer interpretation is reduced to a certain extent, and the efficiency and the accuracy of broken layer interpretation are effectively improved.
Description
Technical Field
The invention belongs to geophysical exploration, and particularly relates to an ant body crack tracking method, device, electronic equipment and medium based on spectrum decomposition.
Background
The crack is not only a main percolation channel of the carbonate reservoir, but also an important reservoir space, the crack has important control effect on the development of the carbonate reservoir, and the crack development zone is often the most developed zone of the carbonate reservoir and is also the most favorable zone for oil gas aggregation, so the crack development zone prediction is the gravity center of the carbonate reservoir research work. Large faults can be identified by 2D or 3D seismic data, and small scale faults or fractures can be identified by log data. For faults with intermediate dimensions of a few decimeters to 30 meters in break distance, neither seismic nor logging data are generally well identifiable. But the faults have very important influence on the reservoir, so that the oil and gas recovery ratio can be improved, and the economic benefit can be obtained to the maximum extent.
The common crack tracking and identifying method mainly comprises a coherent data volume technology, a curvature attribute analysis technology and an ant tracking algorithm. The coherent data volume technology does not need horizon restraint, and the resolution in the longitudinal and transverse directions is improved. The method is widely popularized and applied, and geological phenomena such as faults, river channels and the like can be well identified. But this method is not very capable of identifying small faults. The fine-phase correlation technique utilizes the multi-attribute of coherence, dip angle, azimuth angle and the like to research and identify microstructure, cracks, special lithology and the like. After the seismic data is subjected to relevant processing, an inclination angle body, an azimuth angle body and a coherent body can be obtained and displayed in an overlapping manner, and the subtle changes of the geologic body shape can be more clearly described, so that powerful technical means are provided for geologist research on deformation, fold, crack, lithology change and the like of the structure.
The curvature attribute solves the trouble of geologist, and compared with a coherent algorithm, the curvature attribute has better identification capability on small characteristics such as disturbance, fold, bulge and the like on seismic data. However, curvature analysis is applied to predict the fracture, and only the final structural morphology of the stratum is considered in curvature analysis, and the structural event or the specific structural evolution process of the folds experienced by the reservoir are not considered. For example, for a box-folded wing, there is little change in pitch, the curvature is substantially zero, and the magnitude of the curvature value indicates that the crack density is also low in these places. However, when we consider the evolution of the folds, it is clear that the deformation of the wing is very strong, which will become a region of very high crack density. Thus, in analyzing the curvature properties, consideration must be given to whether the formation undergoes structural deformation or not, otherwise, the interpretation of the curvature properties may deviate greatly.
The ant tracking algorithm is a very effective fracture and crack identification technology which is currently recognized, a boundary map is identified more clearly, the boundary map comprises all small cracks and earthquake discontinuous interfaces, the algorithm utilizes the bionics concept, namely, ants can always climb along the shortest path from nest to foraging, and the ants can communicate with each other by utilizing a chemical substance of a pheromone. Ants along the shortest path will arrive at the destination earlier and then the other ants will arrive successively according to the pheromone on the road. There will be more pheromones on the shortest path. The technology can automatically identify and analyze the fracture system, has the characteristics of high speed, less human intervention and the like, and has higher requirements on the quality of seismic data. However, as the depth of burial of the exploration target area increases, geological conditions tend to be complex, and exploration difficulty increases, and the main expression is as follows: the quality of the seismic data is reduced, the signal-to-noise ratio is low, and the resolution can not meet the requirement of directly using ant tracking to identify the tiny fracture zone.
Therefore, how to suppress noise while keeping effective signals to the greatest extent and guaranteeing the quality of deep seismic data becomes a key for identifying deep cracks by using an ant tracking technology, and in this regard, an ant tracking method capable of identifying small crack zones is particularly needed.
Disclosure of Invention
The invention aims to provide an ant body crack tracking method based on spectrum decomposition, which solves the problem that a tiny crack zone cannot be identified by directly applying the ant tracking method to the existing seismic data.
In view of the above, the invention provides an ant body crack tracking method, an ant body crack tracking device, an electronic device and a medium based on spectrum decomposition, which at least solve the problem that the existing seismic data cannot identify a tiny crack zone by directly applying the ant tracking method.
In a first aspect, the present invention provides a method for tracking ant body cracks based on spectral decomposition, including: performing spectrum decomposition on the full-band post-stack seismic data to obtain a seismic tuning body; determining a main frequency band of a seismic data effective signal based on the seismic tuning body; performing spectrum decomposition on the seismic data in the main frequency band to obtain a discrete single-frequency energy body; and carrying out ant body tracking calculation on the discrete single-frequency energy bodies of the low frequency band and the discrete single-frequency energy bodies of the high frequency band respectively to obtain a low-frequency-band ant body crack and a high-frequency-band ant body crack.
Optionally, performing spectrum decomposition on the full-band post-stack seismic data, and acquiring the seismic tuning body includes: explaining the full-band post-stack seismic data in a time domain to obtain an explained three-dimensional seismic data volume; acquiring a data subvolume of a target interval based on the interpreted three-dimensional seismic data volume; and converting the data subvolumes of the target interval from a time domain to a frequency domain, and acquiring the seismic tuning body.
Optionally, performing discrete fourier transform on the seismic data in the main frequency band to obtain the discrete single-frequency energy body.
Optionally, the frequency interval of the discrete single frequency energy bodies is 1Hz.
In a second aspect, the present invention also provides an electronic device, including: a memory storing executable instructions; and the processor runs the executable instructions in the memory to realize the ant body crack tracking method based on frequency spectrum decomposition.
In a third aspect, the present invention also provides a computer readable storage medium storing a computer program which, when executed by a processor, implements the above-described ant body crack tracking method based on spectral decomposition.
In a fourth aspect, the present invention further provides an ant body crack tracking device based on spectral decomposition, including: the frequency spectrum decomposition module is used for performing frequency spectrum decomposition on the full-band post-stack seismic data to obtain a seismic tuning body; a main frequency band determining module for determining a main frequency band of the effective signal of the seismic data based on the seismic tuning body; the energy body acquisition module is used for carrying out frequency spectrum decomposition on the seismic data in the main frequency band to acquire discrete single-frequency energy bodies; and the ant body tracking module is used for respectively carrying out ant body tracking calculation on the discrete single-frequency energy bodies of the low frequency band and the discrete single-frequency energy bodies of the high frequency band to obtain a low-frequency-band ant body crack and a high-frequency-band ant body crack.
Optionally, performing spectrum decomposition on the full-band post-stack seismic data, and acquiring the seismic tuning body includes: explaining the full-band post-stack seismic data in a time domain to obtain an explained three-dimensional seismic data volume; acquiring a data subvolume of a target interval based on the interpreted three-dimensional seismic data volume; and converting the data subvolumes of the target interval from a time domain to a frequency domain, and acquiring the seismic tuning body.
Optionally, performing discrete fourier transform on the seismic data in the main frequency band to obtain the discrete single-frequency energy body.
Optionally, the frequency interval of the discrete single frequency energy bodies is 1Hz.
The invention has the beneficial effects that: when the method for tracking the ant body cracks based on frequency spectrum decomposition is used for carrying out geological discontinuity imaging and interpretation on three-dimensional seismic data, the frequency spectrum decomposition technology is used for converting full-band post-stack seismic data from a time domain to a frequency domain, a seismic tuning body and a discrete single-frequency energy body are calculated, then the low-frequency and high-frequency discrete single-frequency energy bodies are optimized for carrying out ant body tracking calculation, low-frequency-band ant body cracks and high-frequency-band ant body cracks are obtained, mutual interference of different frequency components in a time domain is eliminated, so that interpretation results higher than the traditional resolution are obtained, fine cracks and crack boundaries can be identified more clearly, meanwhile, the influence of man-made subjectivity in broken layer interpretation is reduced to a certain extent, and the efficiency and the accuracy of broken layer interpretation are effectively improved.
The invention has other features and advantages which will be apparent from or are set forth in detail in the accompanying drawings and the following detailed description, which are incorporated herein, and which together serve to explain certain principles of the invention.
Drawings
The foregoing and other objects, features and advantages of the invention will be apparent from the following more particular descriptions of exemplary embodiments of the invention as illustrated in the accompanying drawings wherein like reference numbers generally represent like parts throughout the exemplary embodiments of the invention.
Fig. 1 shows a flowchart of a method for tracking ant body cracks based on spectral decomposition according to an embodiment of the present invention.
Fig. 2 shows a flowchart of a method for spectrum decomposition-based ant body fracture tracking that applies full-band post-stack seismic data to generate a seismic tuning body, in accordance with one embodiment of the invention.
Fig. 3 shows a seismic tuning volume diagram of an ant body fracture tracking method based on spectral decomposition according to an embodiment of the invention.
Fig. 4 shows a flowchart of a method for spectrum decomposition-based ant body fracture tracking using full-band post-stack seismic data to generate discrete single-frequency energy bodies, in accordance with one embodiment of the present invention.
Fig. 5 shows a10 Hz single frequency energy volume plan view of an ant body crack tracking method based on spectral decomposition according to an embodiment of the present invention.
Fig. 6 shows a 20Hz single frequency energy volume plan view of an ant body crack tracking method based on spectral decomposition according to an embodiment of the present invention.
Fig. 7 shows a 30Hz single frequency energy volume plan view of an ant body crack tracking method based on spectral decomposition according to an embodiment of the present invention.
Fig. 8 shows a 40Hz single frequency energy volume plan view of an ant body crack tracking method based on spectral decomposition according to an embodiment of the present invention.
Fig. 9 illustrates a10 Hz ant body crack recognition plan view of a spectrum decomposition-based ant body crack tracking method according to an embodiment of the present invention.
Fig. 10 illustrates a 30Hz ant body crack recognition plan view of a spectrum decomposition-based ant body crack tracking method according to an embodiment of the present invention.
Fig. 11 shows an ant body crack identification plan view calculated from full-band seismic data of a top surface of a target layer based on a spectrum decomposition ant body crack tracking method according to an embodiment of the present invention.
FIG. 12 shows a fracture strike rose plot for well imaging data statistics in a work area.
Fig. 13 shows a block diagram of an ant body crack tracking device based on spectral decomposition according to an embodiment of the present invention.
102. A spectrum decomposition module; 104. a primary band determination module; 106. an energy body acquisition module; 108. And the ant body tracking module.
Detailed Description
Preferred embodiments of the present invention will be described in more detail below. While the preferred embodiments of the present invention are described below, it should be understood that the present invention may be embodied in various forms and should not be limited to the embodiments set forth herein.
The invention provides an ant body crack tracking method based on spectrum decomposition, which comprises the following steps: performing spectrum decomposition on the full-band post-stack seismic data to obtain a seismic tuning body; determining a main frequency band of the effective signal of the seismic data based on the seismic tuning body; performing spectrum decomposition on the seismic data in the main frequency band to obtain a discrete single-frequency energy body; and carrying out ant body tracking calculation on the discrete single-frequency energy bodies of the low frequency band and the discrete single-frequency energy bodies of the high frequency band respectively to obtain a low-frequency-band ant body crack and a high-frequency-band ant body crack.
Specifically, spectrum decomposition technology is used for carrying out spectrum decomposition on the full-band post-stack seismic data to generate a seismic tuning body, a main frequency band of a seismic data effective signal is determined in the seismic tuning body, spectrum decomposition technology is used for carrying out spectrum decomposition on the seismic data in a main frequency band range, discrete single-frequency energy bodies are generated in the main frequency band range, and a plurality of single-frequency energy bodies are preferentially selected in a low frequency band and a high frequency band to carry out ant body tracking calculation to obtain a low-frequency-band ant body crack and a high-frequency-band ant body crack.
Through the research of the invention, small cracks and crack development zones which are difficult to identify by naked eyes in an original seismic section are clearly reflected on an ant tracking section, the morphology and the spread of the crack development zones are clear and visual, the application effect of noise suppression is good, the definition and the accuracy of crack detection are improved compared with those of the traditional ant tracking, and the research result of the invention provides an effective technical means for the seismic crack prediction of an oil field.
According to the embodiment, when geological discontinuity imaging and interpretation are carried out on three-dimensional seismic data, the method for tracking the ant body cracks based on frequency spectrum decomposition utilizes the frequency spectrum decomposition technology to convert full-band post-stack seismic data from a time domain to a frequency domain, calculates a seismic tuning body and a discrete single-frequency energy body, and then carries out ant body tracking calculation by optimizing the discrete single-frequency energy bodies with low frequency and high frequency to obtain low-frequency-band ant body cracks and high-frequency-band ant body cracks, so that mutual interference of different frequency components in a time domain is eliminated, interpretation results higher than the traditional resolution are obtained, fine cracks and crack boundaries can be identified more clearly, meanwhile, artificial subjective influence in broken layer interpretation is reduced to a certain extent, and the efficiency and accuracy of broken layer interpretation are effectively improved.
As an alternative, performing spectral decomposition on the full-band post-stack seismic data, the obtaining the seismic tuning body includes: explaining the full-band post-stack seismic data in a time domain to obtain an explained three-dimensional seismic data volume; acquiring a data subvolume of the target interval based on the interpreted three-dimensional seismic data volume; and converting the data subvolumes of the target interval from a time domain to a frequency domain, and acquiring the seismic tuning volume.
The tuning body is a data body composed of a plurality of frequency components in a single time window, the vertical upper scale represents continuously variable frequency, and the plane upper scale represents the tuning resonance amplitude corresponding to each single frequency. The heterogeneity of the target layer can be reflected by the trend of the amplitude.
And interpreting the full-band post-stack seismic data in a time domain, acquiring an interpreted three-dimensional seismic data volume, acquiring data subvolumes of a target interval in the interpreted three-dimensional seismic data volume, and converting the data subvolumes of the target interval from the time domain to a frequency domain through Fourier transformation to acquire a seismic tuning volume.
Alternatively, discrete Fourier transform is performed on the seismic data in the main frequency band to obtain discrete single-frequency energy volumes.
Alternatively, the frequency interval of the discrete single frequency energy bodies is 1Hz.
Specifically, discrete Fourier computation is performed on the seismic data volume within the main frequency band, and a discrete single-frequency energy volume with a 1Hz interval is generated.
The present invention also provides an electronic device including: a memory storing executable instructions; and the processor runs executable instructions in the memory to realize the ant body crack tracking method based on frequency spectrum decomposition.
The invention also provides a computer readable storage medium storing a computer program which when executed by a processor implements the ant body crack tracking method based on spectrum decomposition.
The invention also provides an ant body crack tracking device based on frequency spectrum decomposition, which comprises: the frequency spectrum decomposition module is used for performing frequency spectrum decomposition on the full-band post-stack seismic data to obtain a seismic tuning body; the main frequency band determining module is used for determining a main frequency band of the effective signal of the seismic data based on the seismic tuning body; the energy body acquisition module is used for carrying out frequency spectrum decomposition on the seismic data in the main frequency band to acquire discrete single-frequency energy bodies; and the ant body tracking module is used for respectively carrying out ant body tracking calculation on the discrete single-frequency energy bodies of the low frequency band and the discrete single-frequency energy bodies of the high frequency band to obtain a low-frequency-band ant body crack and a high-frequency-band ant body crack.
Specifically, spectrum decomposition technology is used for carrying out spectrum decomposition on the full-band post-stack seismic data to generate a seismic tuning body, a main frequency band of a seismic data effective signal is determined in the seismic tuning body, spectrum decomposition technology is used for carrying out spectrum decomposition on the seismic data in a main frequency band range, discrete single-frequency energy bodies are generated in the main frequency band range, and a plurality of single-frequency energy bodies are preferentially selected in a low frequency band and a high frequency band to carry out ant body tracking calculation to obtain a low-frequency-band ant body crack and a high-frequency-band ant body crack.
Through the research of the invention, small cracks and crack development zones which are difficult to identify by naked eyes in an original seismic section are clearly reflected on an ant tracking section, the morphology and the spread of the crack development zones are clear and visual, the application effect of noise suppression is good, the definition and the accuracy of crack detection are improved compared with those of the traditional ant tracking, and the research result of the invention provides an effective technical means for the seismic crack prediction of an oil field.
According to the embodiment, when the geological discontinuity imaging and interpretation are carried out on the three-dimensional seismic data, the full-band stacked seismic data are transformed from a time domain to a frequency domain by utilizing a frequency spectrum decomposition technology, the seismic tuning body and the discrete single-frequency energy body are calculated, then the low-frequency and high-frequency discrete single-frequency energy bodies are optimized to carry out ant tracking calculation, the low-frequency ant body cracks and the high-frequency ant body cracks are obtained, the mutual interference of different frequency components in a time domain is eliminated, so that interpretation results higher than the traditional resolution are obtained, the tiny cracks and crack boundaries can be identified more clearly, meanwhile, the artificial subjective influence in the broken layer interpretation is reduced to a certain extent, and the efficiency and the accuracy of the broken layer interpretation are effectively improved.
As an alternative, performing spectral decomposition on the full-band post-stack seismic data, the obtaining the seismic tuning body includes: explaining the full-band post-stack seismic data in a time domain to obtain an explained three-dimensional seismic data volume; acquiring a data subvolume of the target interval based on the interpreted three-dimensional seismic data volume; and converting the data subvolumes of the target interval from a time domain to a frequency domain, and acquiring the seismic tuning volume.
The tuning body is a data body composed of a plurality of frequency components in a single time window, the vertical upper scale represents continuously variable frequency, and the plane upper scale represents the tuning resonance amplitude corresponding to each single frequency. The heterogeneity of the target layer can be reflected by the trend of the amplitude.
And interpreting the full-band post-stack seismic data in a time domain, acquiring an interpreted three-dimensional seismic data volume, acquiring data subvolumes of a target interval in the interpreted three-dimensional seismic data volume, and converting the data subvolumes of the target interval from the time domain to a frequency domain through Fourier transformation to acquire a seismic tuning volume.
Alternatively, discrete Fourier transform is performed on the seismic data in the main frequency band to obtain discrete single-frequency energy volumes.
Alternatively, the frequency interval of the discrete single frequency energy bodies is 1Hz.
Specifically, discrete Fourier computation is performed on the seismic data volume within the main frequency band, and a discrete single-frequency energy volume with a 1Hz interval is generated.
Example 1
Fig. 1 shows a flowchart of a method for tracking ant body cracks based on spectral decomposition according to an embodiment of the present invention. Fig. 2 shows a flowchart of a method for spectrum decomposition-based ant body fracture tracking that applies full-band post-stack seismic data to generate a seismic tuning body, in accordance with one embodiment of the invention. Fig. 3 shows a seismic tuning volume diagram of an ant body fracture tracking method based on spectral decomposition according to an embodiment of the invention. Fig. 4 shows a flowchart of a method for spectrum decomposition-based ant body fracture tracking using full-band post-stack seismic data to generate discrete single-frequency energy bodies, in accordance with one embodiment of the present invention. Fig. 5 shows a 10Hz single frequency energy volume plan view of an ant body crack tracking method based on spectral decomposition according to an embodiment of the present invention. Fig. 6 shows a 20Hz single frequency energy volume plan view of an ant body crack tracking method based on spectral decomposition according to an embodiment of the present invention. Fig. 7 shows a 30Hz single frequency energy volume plan view of an ant body crack tracking method based on spectral decomposition according to an embodiment of the present invention. Fig. 8 shows a 40Hz single frequency energy volume plan view of an ant body crack tracking method based on spectral decomposition according to an embodiment of the present invention. Fig. 9 illustrates a 10Hz ant body crack recognition plan view of a spectrum decomposition-based ant body crack tracking method according to an embodiment of the present invention. Fig. 10 illustrates a 30Hz ant body crack recognition plan view of a spectrum decomposition-based ant body crack tracking method according to an embodiment of the present invention. Fig. 11 shows an ant body crack identification plan view calculated from full-band seismic data of a top surface of a target layer based on a spectrum decomposition ant body crack tracking method according to an embodiment of the present invention. FIG. 12 shows a fracture strike rose plot for well imaging data statistics in a work area.
As shown in fig. 1, 2,3, 4,5, 6, 7, 8, 9, 10, 11 and 12, the ant body crack tracking method based on spectrum decomposition includes:
Step 1: performing spectrum decomposition on the full-band post-stack seismic data to obtain a seismic tuning body;
The method for obtaining the seismic tuning body comprises the following steps of: explaining the full-band post-stack seismic data in a time domain to obtain an explained three-dimensional seismic data volume; acquiring a data subvolume of the target interval based on the interpreted three-dimensional seismic data volume; and converting the data subvolumes of the target interval from a time domain to a frequency domain, and acquiring the seismic tuning volume.
The tuning body is a data body composed of a plurality of frequency components in a single time window, the vertical upper scale represents continuously variable frequency, and the plane upper scale represents the tuning resonance amplitude corresponding to each single frequency. The heterogeneity of the target layer can be reflected by the trend of the amplitude.
And interpreting the full-band post-stack seismic data in a time domain, acquiring an interpreted three-dimensional seismic data volume, acquiring data subvolumes of a target interval in the interpreted three-dimensional seismic data volume, converting the data subvolumes of the target interval from the time domain to a frequency domain through Fourier transformation, and acquiring a seismic tuning volume, as shown in figure 2.
Step 2: determining a main frequency band of the effective signal of the seismic data based on the seismic tuning body;
Taking a room group stratum of a Tarim basin YB work area as an example, taking a top interface T74+40ms as a center, setting a time window to 140ms (basically including a room group stratum), and performing discrete Fourier transform on the post-stack seismic data to obtain a tuning body of the post-stack seismic data in a frequency domain. Fig. 3 is a cross-sectional view of a tuning body of any measuring line in the north-west direction of an over-work area, and it can be seen that data above 60Hz are very noise-like, the effective signal energy of a stratum of a room group is mainly concentrated at 10-20 Hz, and a part of the area is about 40Hz, so that the main frequency band range of the seismic data of a target interval of the work area is determined to be 10-40 Hz.
Step 3: performing spectrum decomposition on the seismic data in the main frequency band to obtain a discrete single-frequency energy body;
and performing discrete Fourier transform on the seismic data in the main frequency band to obtain a discrete single-frequency energy body.
Wherein the frequency interval of the discrete single frequency energy bodies is 1Hz.
Taking a room group stratum of a Tarim basin YB work area as an example, selecting a seismic data body with a spectral decomposition frequency band range of 10-40HZ for discrete Fourier calculation to generate a discrete single-frequency energy body with an interval of 1 Hz. Fig. 4 is a workflow of forming discrete single frequency energy volumes. The frequency analysis method has the advantages that the influence of the structural form on interpretation can be eliminated by utilizing a sliding time window method along a layer, and the control and influence of the horizon can be avoided by utilizing an isochronous window analysis method. Fig. 5-8 are single frequency energy volume plan views of 10Hz, 20Hz, 30Hz, 40Hz, respectively. The characteristics of the seismic profile are obviously different, and the small differences of the seismic amplitudes of different frequency bands caused by lithology, physical properties and other factors among the strata of one room group can be amplified by a spectrum decomposition technology, and the small differences are displayed and studied.
Step 4: and carrying out ant body tracking calculation on the discrete single-frequency energy bodies of the low frequency band and the discrete single-frequency energy bodies of the high frequency band respectively to obtain a low-frequency-band ant body crack and a high-frequency-band ant body crack.
And preliminarily checking crack characteristics reflected by each frequency through the discrete single-frequency energy bodies, and performing ant body calculation on the single-frequency energy bodies in the low frequency band and the high frequency band respectively. Through experimental analysis, single-frequency ant bodies of 10Hz and 30Hz are finally selected. In fig. 9, for the 10Hz ant body crack detection plan of the target layer of the research area, gray is a small-scale fault and cracks which can be identified on the core, and it can be seen that the 10Hz ant body can finely depict the middle-large scale cracks in the large fault, and fig. 10 is the 30Hz ant body crack detection plan of the target layer of the research area, the darker the color in the drawing, the higher the reliability of crack identification is represented, and the 30Hz ant body depicts more details of the small fault and micro cracks near the fault. 11 is the ant body of the full frequency band of the seismic data, can see that the mutual interference of different frequency components has better identification effect on large fracture and is obvious in small-scale crack suppression. The invention provides more fault and crack information than the ant body calculated by full-band seismic data.
And verifying the application effect of the invention by combining the fracture trend characteristics revealed by imaging logging, ancient stress field and other data in the research area. The developed fracture group lines of YJ1-9 and YC1 wells revealed by the imaging well logging of FIG. 12 mainly run in the northeast direction, and are consistent with the fracture identification result of the invention.
Example two
Fig. 13 shows a block diagram of an ant body crack tracking device based on spectral decomposition according to an embodiment of the present invention.
As shown in fig. 13, the ant body crack tracking device based on spectrum decomposition includes:
the spectrum decomposition module 102 performs spectrum decomposition on the full-band post-stack seismic data to obtain a seismic tuning body;
The method for obtaining the seismic tuning body comprises the following steps of: explaining the full-band post-stack seismic data in a time domain to obtain an explained three-dimensional seismic data volume; acquiring a data subvolume of the target interval based on the interpreted three-dimensional seismic data volume; and converting the data subvolumes of the target interval from a time domain to a frequency domain, and acquiring the seismic tuning volume.
A main band determination module 104 that determines a main band of the seismic data effective signal based on the seismic tuning volume;
the energy body acquisition module 106 is used for carrying out frequency spectrum decomposition on the seismic data in the main frequency band to acquire discrete single-frequency energy bodies;
and performing discrete Fourier transform on the seismic data in the main frequency band to obtain a discrete single-frequency energy body.
Wherein the frequency interval of the discrete single frequency energy bodies is 1Hz.
The ant tracking module 108 performs ant tracking calculation on the low-frequency discrete single-frequency energy and the high-frequency discrete single-frequency energy to obtain a low-frequency ant crack and a high-frequency ant crack.
The present disclosure provides an electronic device including: a memory storing executable instructions; and the processor runs executable instructions in the memory to realize the ant body crack tracking method based on frequency spectrum decomposition.
An electronic device according to an embodiment of the present disclosure includes a memory and a processor.
The memory is for storing non-transitory computer readable instructions. In particular, the memory may include one or more computer program products, which may include various forms of computer-readable storage media, such as volatile memory and/or non-volatile memory. The volatile memory may include, for example, random Access Memory (RAM) and/or cache memory (cache), and the like. The non-volatile memory may include, for example, read Only Memory (ROM), hard disk, flash memory, and the like.
The processor may be a Central Processing Unit (CPU) or other form of processing unit having data processing and/or instruction execution capabilities, and may control other components in the electronic device to perform the desired functions. In one embodiment of the present disclosure, the processor is configured to execute the computer readable instructions stored in the memory.
It should be understood by those skilled in the art that, in order to solve the technical problem of how to obtain a good user experience effect, the present embodiment may also include well-known structures such as a communication bus, an interface, and the like, and these well-known structures are also included in the protection scope of the present disclosure.
The detailed description of the present embodiment may refer to the corresponding description in the foregoing embodiments, and will not be repeated herein.
The present disclosure provides a computer-readable storage medium storing a computer program which, when executed by a processor, implements the above-described ant body crack tracking method based on spectral decomposition.
A computer-readable storage medium according to an embodiment of the present disclosure has stored thereon non-transitory computer-readable instructions. When executed by a processor, perform all or part of the steps of the methods of embodiments of the present disclosure described above.
The computer-readable storage medium described above includes, but is not limited to: optical storage media (e.g., CD-ROM and DVD), magneto-optical storage media (e.g., MO), magnetic storage media (e.g., magnetic tape or removable hard disk), media with built-in rewritable non-volatile memory (e.g., memory card), and media with built-in ROM (e.g., ROM cartridge).
The foregoing description of embodiments of the invention has been presented for purposes of illustration and description, and is not intended to be exhaustive or limited to the embodiments disclosed. Many modifications and variations will be apparent to those of ordinary skill in the art without departing from the scope and spirit of the various embodiments described.
Claims (4)
1. The ant body crack tracking method based on frequency spectrum decomposition is characterized by comprising the following steps of:
Performing spectrum decomposition on the full-band post-stack seismic data to obtain a seismic tuning body;
Determining a main frequency band of a seismic data effective signal based on the seismic tuning body;
performing spectrum decomposition on the seismic data in the main frequency band to obtain a discrete single-frequency energy body;
Carrying out ant body tracking calculation on the discrete single-frequency energy bodies of the low frequency band and the discrete single-frequency energy bodies of the high frequency band respectively to obtain a low-frequency-band ant body crack and a high-frequency-band ant body crack;
the performing spectrum decomposition on the full-band post-stack seismic data to obtain a seismic tuning body includes:
Explaining the full-band post-stack seismic data in a time domain to obtain an explained three-dimensional seismic data volume;
Acquiring a data subvolume of a target interval based on the interpreted three-dimensional seismic data volume;
Converting the data subvolumes of the target interval into frequency domains from time domains to obtain an earthquake tuning body, wherein the tuning body is a data body consisting of a plurality of frequency components in a single time window, vertical upper scales represent continuously-changed frequencies, and plane upper scales represent tuning resonance amplitudes corresponding to each single frequency;
Performing discrete Fourier transform on the seismic data in the main frequency band to obtain the discrete single-frequency energy body;
The frequency interval of the discrete single frequency energy bodies is 1Hz.
2. An electronic device, the electronic device comprising:
a memory storing executable instructions;
A processor executing the executable instructions in the memory to implement the spectral decomposition-based ant body crack tracking method of claim 1.
3. A computer-readable storage medium, characterized in that the computer-readable storage medium stores a computer program which, when executed by a processor, implements the ant body crack tracking method based on spectral decomposition according to claim 1.
4. Ant body crack tracking means based on spectral decomposition, characterized by comprising:
The frequency spectrum decomposition module is used for performing frequency spectrum decomposition on the full-band post-stack seismic data to obtain a seismic tuning body;
A main frequency band determining module for determining a main frequency band of the effective signal of the seismic data based on the seismic tuning body;
the energy body acquisition module is used for carrying out frequency spectrum decomposition on the seismic data in the main frequency band to acquire discrete single-frequency energy bodies;
The ant body tracking module is used for respectively carrying out ant body tracking calculation on the discrete single-frequency energy bodies of the low frequency band and the discrete single-frequency energy bodies of the high frequency band to obtain a low-frequency-band ant body crack and a high-frequency-band ant body crack; the performing spectrum decomposition on the full-band post-stack seismic data to obtain a seismic tuning body includes:
Explaining the full-band post-stack seismic data in a time domain to obtain an explained three-dimensional seismic data volume;
Acquiring a data subvolume of a target interval based on the interpreted three-dimensional seismic data volume;
Converting the data subvolumes of the target interval into frequency domains from time domains to obtain an earthquake tuning body, wherein the tuning body is a data body consisting of a plurality of frequency components in a single time window, vertical upper scales represent continuously-changed frequencies, and plane upper scales represent tuning resonance amplitudes corresponding to each single frequency;
Performing discrete Fourier transform on the seismic data in the main frequency band to obtain the discrete single-frequency energy body;
The frequency interval of the discrete single frequency energy bodies is 1Hz.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010968413.4A CN114185091B (en) | 2020-09-15 | 2020-09-15 | Ant body crack tracking method and device based on frequency spectrum decomposition and electronic equipment |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010968413.4A CN114185091B (en) | 2020-09-15 | 2020-09-15 | Ant body crack tracking method and device based on frequency spectrum decomposition and electronic equipment |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN114185091A CN114185091A (en) | 2022-03-15 |
| CN114185091B true CN114185091B (en) | 2024-05-07 |
Family
ID=80539756
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010968413.4A Active CN114185091B (en) | 2020-09-15 | 2020-09-15 | Ant body crack tracking method and device based on frequency spectrum decomposition and electronic equipment |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114185091B (en) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115932967B (en) * | 2023-01-06 | 2023-07-18 | 东北石油大学三亚海洋油气研究院 | Crack identification method based on ant body algorithm |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2956685B1 (en) * | 1998-03-24 | 1999-10-04 | 日本電信電話株式会社 | Observation method of earthquake precursor signals |
| CN103135136A (en) * | 2011-11-25 | 2013-06-05 | 中国石油化工股份有限公司 | Automatic fault interpretation device for three-dimensional seismic data body |
| CN103675913A (en) * | 2012-09-04 | 2014-03-26 | 中国石油天然气集团公司 | Detection method for fractured reservoirs |
| CN105929452A (en) * | 2016-06-02 | 2016-09-07 | 中国石油天然气股份有限公司 | Method and device for predicting underground crack spatial distribution based on seismic data |
| US10139508B1 (en) * | 2016-03-24 | 2018-11-27 | EMC IP Holding Company LLC | Methods and apparatus for automatic identification of faults on noisy seismic data |
| CN109425891A (en) * | 2017-08-28 | 2019-03-05 | 中国石油化工股份有限公司 | A kind of fracture imaging quality detecting method based on geological model |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7248052B2 (en) * | 2003-05-28 | 2007-07-24 | Weaver W Barry | Electric power grid induced geophysical prospecting method and apparatus |
| US9715026B2 (en) * | 2011-03-11 | 2017-07-25 | Schlumberger Technology Corporation | System and method for performing microseismic fracture operations |
-
2020
- 2020-09-15 CN CN202010968413.4A patent/CN114185091B/en active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2956685B1 (en) * | 1998-03-24 | 1999-10-04 | 日本電信電話株式会社 | Observation method of earthquake precursor signals |
| CN103135136A (en) * | 2011-11-25 | 2013-06-05 | 中国石油化工股份有限公司 | Automatic fault interpretation device for three-dimensional seismic data body |
| CN103675913A (en) * | 2012-09-04 | 2014-03-26 | 中国石油天然气集团公司 | Detection method for fractured reservoirs |
| US10139508B1 (en) * | 2016-03-24 | 2018-11-27 | EMC IP Holding Company LLC | Methods and apparatus for automatic identification of faults on noisy seismic data |
| CN105929452A (en) * | 2016-06-02 | 2016-09-07 | 中国石油天然气股份有限公司 | Method and device for predicting underground crack spatial distribution based on seismic data |
| CN109425891A (en) * | 2017-08-28 | 2019-03-05 | 中国石油化工股份有限公司 | A kind of fracture imaging quality detecting method based on geological model |
Non-Patent Citations (1)
| Title |
|---|
| 《分频蚂蚁追踪在识别深层小断层中的应用》;张瑞等;地球物理学进展;第32卷(第1期);第350-355页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN114185091A (en) | 2022-03-15 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| EP2864817B1 (en) | Seismic orthogonal decomposition attribute | |
| US7391675B2 (en) | Microseismic event detection and location by continuous map migration | |
| US10324211B2 (en) | Seismic spectral balancing | |
| JP3390445B2 (en) | Seismic trace analysis method using statistical calibration techniques to estimate geological properties | |
| US20180203144A1 (en) | Interferometric Microseismic Imaging Methods and Apparatus | |
| US20100036614A1 (en) | Locating oil or gas passively by observing a porous oil and gas saturated system giving off its characteristic resonance response to ambient background noise, including optional differentiation of oil, locatinggas and water | |
| CN109425896A (en) | Dolomite oil and gas reservoir distribution forecasting method and device | |
| US20150168574A1 (en) | Seismic trace attribute | |
| CN104360382A (en) | Method for detecting oil and gas by aid of stacked seismic data | |
| US20170254909A1 (en) | Method for determining maximum horizontal stress magnitude and direction using microseismic derived fracture attributes and its application to evaluating hydraulic fracture stimulation induced stress changes | |
| CA2961461C (en) | Method for determining formation stress field using microseismic focal mechanisms | |
| CN114185091B (en) | Ant body crack tracking method and device based on frequency spectrum decomposition and electronic equipment | |
| CN105954800B (en) | Utilize the method for earthquake texture features identification craven fault | |
| CN113970785A (en) | Method and system for predicting development of underground river crack, storage medium and electronic equipment | |
| Raef et al. | Multichannel analysis of surface-waves and integration of downhole acoustic televiewer imaging, ultrasonic Vs and Vp, and vertical seismic profiling in an NEHRP-standard classification, South of Concordia, Kansas, USA | |
| CN113640876A (en) | Method for finely identifying trapping column by using chaotic body attribute | |
| CN112764100B (en) | Reservoir oil gas range detection method and device | |
| RU2789759C1 (en) | Method for determining the boundaries of the fractured zone | |
| CN110941029B (en) | Speed modeling method related to geological capping | |
| Tian et al. | Geology-geophysics-data mining solution enhances the identification of deep fault-karst carbonate reservoirs | |
| Bao | Intra-reservoir sand architecture delineation via depth domain modified cepstral decomposition “seismic resonance” | |
| Lv et al. | Integrated characterization of deep karsted carbonates in the Tahe Oilfield, Tarim Basin | |
| Gabry et al. | Advanced Deep Learning for microseismic events prediction for hydraulic fracture treatment via Continuous Wavelet Transform | |
| CN114879255A (en) | Granite cave body prediction method and device for submarine cave | |
| Sardinal et al. | Site classification based on spectral amplification patterns for microtremor H/V ratios |
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 | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |