CN113700466B - Method, device, equipment and medium for detecting oil gas in deep carbonate reservoir - Google Patents
Method, device, equipment and medium for detecting oil gas in deep carbonate reservoir Download PDFInfo
- Publication number
- CN113700466B CN113700466B CN202010440587.3A CN202010440587A CN113700466B CN 113700466 B CN113700466 B CN 113700466B CN 202010440587 A CN202010440587 A CN 202010440587A CN 113700466 B CN113700466 B CN 113700466B
- Authority
- CN
- China
- Prior art keywords
- azimuth
- stack
- post
- data
- fault
- 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
- BVKZGUZCCUSVTD-UHFFFAOYSA-L Carbonate Chemical compound [O-]C([O-])=O BVKZGUZCCUSVTD-UHFFFAOYSA-L 0.000 title claims abstract description 42
- 238000000034 method Methods 0.000 title claims abstract description 18
- 238000001514 detection method Methods 0.000 claims abstract description 93
- 238000004458 analytical method Methods 0.000 claims abstract description 45
- 238000009826 distribution Methods 0.000 claims abstract description 23
- 239000004215 Carbon black (E152) Substances 0.000 claims description 30
- 229930195733 hydrocarbon Natural products 0.000 claims description 30
- 150000002430 hydrocarbons Chemical class 0.000 claims description 24
- 238000003860 storage Methods 0.000 claims description 13
- 230000004044 response Effects 0.000 claims description 12
- 238000012512 characterization method Methods 0.000 claims description 9
- 238000004590 computer program Methods 0.000 claims description 9
- 238000012216 screening Methods 0.000 claims description 9
- -1 carbonate hydrocarbon Chemical class 0.000 claims description 6
- 238000007405 data analysis Methods 0.000 claims description 5
- 238000012360 testing method Methods 0.000 claims 1
- 238000004364 calculation method Methods 0.000 abstract description 3
- 238000010586 diagram Methods 0.000 description 14
- 238000010183 spectrum analysis Methods 0.000 description 7
- 238000005516 engineering process Methods 0.000 description 6
- 238000012545 processing Methods 0.000 description 6
- 238000010276 construction Methods 0.000 description 3
- 238000001228 spectrum Methods 0.000 description 3
- 230000005856 abnormality Effects 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 230000002349 favourable effect Effects 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 238000005457 optimization Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 230000002159 abnormal effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000005012 migration Effects 0.000 description 1
- 238000013508 migration Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 238000004321 preservation Methods 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 239000011435 rock Substances 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
Classifications
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B47/00—Survey of boreholes or wells
-
- 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
-
- 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/306—Analysis for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles
-
- 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
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/30—Assessment of water resources
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Remote Sensing (AREA)
- Geology (AREA)
- Environmental & Geological Engineering (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geophysics (AREA)
- Acoustics & Sound (AREA)
- General Physics & Mathematics (AREA)
- Mining & Mineral Resources (AREA)
- Fluid Mechanics (AREA)
- Geochemistry & Mineralogy (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The embodiment of the application provides a method, a device, equipment and a medium for detecting oil gas in a deep carbonate reservoir, wherein the method comprises the following steps: determining the number of azimuth according to the distribution characteristics of the main fault and the seismic acquisition parameters of the region to be detected; dividing the orientation of the seismic pre-stack gather data of the area to be detected according to the number of the orientations, and acquiring post-stack orientation data bodies corresponding to the orientations according to the seismic pre-stack gather data corresponding to the orientations, wherein the orientations comprise orientations of parallel main fault strike and orientations of vertical main fault strike; determining the post-stack azimuth data body with the minimum anisotropism in the post-stack azimuth data bodies corresponding to the plurality of azimuths; and carrying out oil gas detection based on the determined post-stack azimuth data body. According to the scheme, the data volume of oil gas detection is reduced, time consumption of oil gas detection calculation and analysis is reduced, detection efficiency is improved, influence of anisotropy on oil gas detection is reduced, and detection accuracy is improved.
Description
Technical Field
The application relates to the technical field of petroleum geophysical exploration, in particular to a method, a device, equipment and a medium for detecting oil gas in a deep carbonate reservoir.
Background
In recent years, with the great breakthrough of a plurality of basins in China and western China in deep carbonate oil and gas exploration, deep carbonate oil and gas reservoirs have become one of important exploration fields for oil and gas production and up-storage in the future of China. Foreign scholars define deep oil gas as oil gas resources buried below 4000 meters, and domestic scholars mostly define oil gas reservoirs larger than 4500 meters as deep oil gas reservoirs. The deep carbonate hydrocarbon reservoir has the characteristics of large reserve scale and high yield, but the deep hydrocarbon resource exploration rate in China is far lower than that of the middle and shallow hydrocarbon resource exploration rate, and the effective exploration of the deep carbonate hydrocarbon reservoir is important to our attention. Seismic technology has played an important role in oil and gas exploration, but the following research difficulties still exist in the deep carbonate oil and gas detection: (1) because the deep carbonate reservoir is buried deeply, the earthquake wave attenuation is serious, the signal-to-noise ratio and the frequency of the earthquake data are low, and the oil gas response characteristics are not obvious; (2) the hydrocarbon response characteristics are complicated due to the rapid and heterogeneous lateral lithology changes of the carbonate reservoir due to the development of cracks and holes.
Hydrocarbon detection is an important parameter for the optimization of a favorable target in hydrocarbon exploration, and the use of seismic data for hydrocarbon detection has long been an important research content in hydrocarbon exploration. Currently, oil and gas detection techniques based on seismic data can be divided into two major categories, pre-stack and post-stack. The pre-stack oil gas detection technology mainly utilizes pre-stack seismic data to carry out pre-stack inversion, fluid identification, anisotropic analysis and the like; the post-stack seismic oil-gas detection technology mainly utilizes abnormal response of post-stack seismic data caused by an oil-gas reservoir to carry out oil-gas detection, and the current common technology mainly comprises amplitude abnormality, frequency abnormality and the like, and is most common in oil-gas exploration due to high overall signal-to-noise ratio and small data volume of the post-stack seismic data.
Although the two methods play an important role in oil and gas detection, the following defects still exist for deep oil and gas detection, 1) because pre-stack seismic data are large, time consumption for calculation and analysis of the pre-stack oil and gas detection is long, and defects of weak energy, tensile distortion and the like exist in a near offset distance and a large offset distance of a pre-stack gather respectively, so that an oil and gas detection result has certain uncertainty. 2) The post-stack oil gas detection is not enough in accuracy on deep carbonate rock oil gas detection under the influence of seismic data, prediction errors can be obviously increased, and deep oil gas exploration cannot be effectively guided.
Disclosure of Invention
The embodiment of the application provides a method for detecting oil gas in a deep carbonate oil gas reservoir, which aims to solve the technical problems of time consumption and low precision in deep carbonate oil gas detection in the prior art. The method comprises the following steps:
determining the number of azimuth according to the distribution characteristics of the main fault and the seismic acquisition parameters of the region to be detected;
dividing the data of the seismic pre-stack gathers of the area to be detected according to the number of the orientations, and acquiring post-stack orientation data bodies corresponding to the orientations according to the data of the seismic pre-stack gathers corresponding to the orientations, wherein the orientations comprise orientations of parallel main fault strike and orientations of vertical main fault strike;
determining the post-stack azimuth data body with the minimum anisotropism in the post-stack azimuth data bodies corresponding to the plurality of azimuths;
oil gas detection is carried out based on the determined post-stack azimuth data body;
determining the post-stack azimuth data body with the minimum anisotropism in the post-stack azimuth data bodies corresponding to the plurality of azimuths, wherein the post-stack azimuth data body comprises the following steps:
performing amplitude analysis and fault characterization on the post-stack azimuth data body corresponding to each azimuth;
AVO analysis is carried out on the seismic prestack gather data of the area to be detected;
determining a post-stack orientation data volume satisfying the following condition among the plurality of orientations as a post-stack orientation data volume having minimum anisotropy:
the amplitude is maximum, the fault depiction is weakest, and the AVO analysis oil gas response characteristics meet the requirements.
The embodiment of the application also provides a device for detecting the oil gas in the deep carbonate oil gas reservoir, which aims to solve the technical problems of time consumption and low precision in the detection of the deep carbonate oil gas in the prior art. The device comprises:
the azimuth number determining module is used for determining the azimuth number according to the main fault distribution characteristics and the seismic acquisition parameters of the region to be detected;
the azimuth dividing module is used for dividing the earthquake pre-stack gather data of the area to be detected according to the azimuth number, and acquiring post-stack azimuth dividing data bodies corresponding to a plurality of azimuths according to the earthquake pre-stack gather data corresponding to the azimuths, wherein the azimuths comprise azimuth of parallel main fault strike and azimuth of vertical main fault strike;
the azimuth screening module is used for determining the post-stack azimuth data body with the minimum anisotropism from the post-stack azimuth data bodies corresponding to the plurality of azimuths;
the detection module is used for carrying out oil gas detection based on the determined post-stack azimuth data body;
the azimuth screening module comprises:
the data analysis unit is used for carrying out amplitude analysis and fault depiction on the post-stack azimuth data body corresponding to each azimuth;
the AVO analysis unit is used for carrying out AVO analysis on the seismic prestack gather data of the area to be detected;
an azimuth screening unit configured to determine, as a post-stack azimuth data volume having minimum anisotropy, a post-stack azimuth data volume satisfying the following condition among a plurality of azimuth:
the amplitude is maximum, the fault depiction is weakest, and the AVO analysis oil gas response characteristics meet the requirements.
The embodiment of the application also provides computer equipment, which comprises a memory, a processor and a computer program stored on the memory and capable of running on the processor, wherein the processor realizes the arbitrary deep carbonate reservoir oil gas detection method when executing the computer program so as to solve the technical problems of time consumption and low precision in deep carbonate reservoir oil gas detection in the prior art.
The embodiment of the application also provides a computer readable storage medium which stores a computer program for executing the arbitrary deep carbonate hydrocarbon reservoir hydrocarbon detection method, so as to solve the technical problems of time consumption and low precision in deep carbonate hydrocarbon detection in the prior art.
In the embodiment of the application, the number of the azimuth is determined according to the main fault distribution characteristics and the earthquake acquisition parameters of the area to be detected, the azimuth is further divided into the earthquake pre-stack gather data of the area to be detected according to the number of the azimuth, and the post-stack azimuth data body corresponding to a plurality of azimuth is obtained according to the earthquake pre-stack gather data corresponding to a plurality of azimuth, wherein the plurality of azimuth comprises azimuth parallel to the trend of the main fault and azimuth perpendicular to the trend of the main fault, so that the division of the earthquake pre-stack gather data of the area to be detected into a plurality of azimuth based on the main fault distribution characteristics and the earthquake acquisition parameters of the area to be detected is realized, and the earthquake acquisition parameters and the area construction characteristics of the pre-stack gather are fully considered, so that the division of azimuth has scientificity; meanwhile, in the post-stack azimuth data bodies corresponding to a plurality of azimuth, the post-stack azimuth data body with the minimum anisotropism is determined, and the oil gas detection is carried out based on the determined post-stack azimuth data body.
Drawings
The accompanying drawings, which are included to provide a further understanding of the application and are incorporated in and constitute a part of this specification, illustrate and together with the description serve to explain the application. In the drawings:
FIG. 1 is a flow chart of a method for detecting oil and gas in a deep carbonate reservoir according to an embodiment of the present application;
FIG. 2 is a schematic diagram of an omnidirectional seismic data fault detection provided by an embodiment of the application;
FIG. 3 is a schematic diagram of pre-stack angle gather data after OVT domain migration processing according to an embodiment of the present application;
FIG. 4 is a schematic diagram of amplitude analysis versus ellipse for a split azimuth data volume according to an embodiment of the present application;
FIG. 5 is a schematic diagram of a split azimuth data volume fault detection contrast provided by an embodiment of the present application;
FIG. 6 (a) is a schematic diagram of AVO analysis of a certain azimuth data volume according to an embodiment of the present application;
FIG. 6 (b) is a schematic diagram of an AVO analysis of another azimuth data volume provided by an embodiment of the present application;
FIG. 7 (a) is a schematic diagram of an omnidirectional seismic spectrum analysis provided by an embodiment of the application;
FIG. 7 (b) is a schematic illustration of a dominant azimuthal seismic spectral analysis provided by an embodiment of the present application;
FIG. 8 is a schematic diagram of an omni-directional and dominant directional hydrocarbon detection contrast provided by embodiments of the present application;
FIG. 9 is a schematic diagram of a result of detecting oil and gas by the deep carbonate reservoir oil and gas detection method according to an embodiment of the present application;
FIG. 10 is a block diagram of a computer device according to an embodiment of the present application;
fig. 11 is a block diagram of a deep carbonate reservoir oil-gas detection device according to an embodiment of the present application.
Detailed Description
The present application will be described in further detail with reference to the following embodiments and the accompanying drawings, in order to make the objects, technical solutions and advantages of the present application more apparent. The exemplary embodiments of the present application and the descriptions thereof are used herein to explain the present application, but are not intended to limit the application.
In an embodiment of the present application, a method for detecting oil and gas in a deep carbonate reservoir is provided, as shown in fig. 1, where the method includes:
step 102: determining the number of azimuth according to the distribution characteristics of the main fault and the seismic acquisition parameters of the region to be detected;
step 104: dividing the data of the seismic pre-stack gathers of the area to be detected according to the number of the orientations, and acquiring post-stack orientation data bodies corresponding to the orientations according to the data of the seismic pre-stack gathers corresponding to the orientations, wherein the orientations comprise orientations of parallel main fault strike and orientations of vertical main fault strike;
step 106: determining the post-stack azimuth data body with the minimum anisotropism in the post-stack azimuth data bodies corresponding to the plurality of azimuths;
step 108: and carrying out oil gas detection based on the determined post-stack azimuth data body.
As can be seen from the flow shown in fig. 1, in the embodiment of the present application, the number of azimuth is determined according to the main fault distribution characteristics and the seismic acquisition parameters of the area to be measured, and then the data of the pre-stack trace set of the area to be measured is azimuth divided according to the number of azimuth, and the post-stack azimuth data volume corresponding to the plurality of azimuth is obtained according to the data of the pre-stack trace set of the area to be measured, where the plurality of azimuth includes the azimuth parallel to the main fault trend and the azimuth perpendicular to the main fault trend, so that the data of the pre-stack trace set of the area to be measured is divided into a plurality of azimuth based on the main fault distribution characteristics and the seismic acquisition parameters of the area to be measured, and the pre-stack trace set seismic acquisition parameters and the area construction characteristics are fully considered, so that the azimuth division is scientific; meanwhile, in the post-stack azimuth data bodies corresponding to a plurality of azimuth, the post-stack azimuth data body with the minimum anisotropism is determined, and the oil gas detection is carried out based on the determined post-stack azimuth data body.
In particular, in order to further improve the detection precision, in this embodiment, OVT (wide azimuth vector offset) domain offset processing may be performed on deep carbonate reservoirs, and signal-to-noise ratio improving processing and high-frequency information compensating processing may be performed on pre-stack gather data, so as to finally obtain seismic pre-stack gather data with high fidelity, high resolution and high signal-to-noise ratio, as shown in fig. 3.
In particular implementation, in order to further improve the scientificity of dividing the data of the pre-stack trace set of the earthquake, in this embodiment, the omnidirectional fault prediction can be performed on the omnidirectional seismic data body of the area to be measured, by extracting the coherence, curvature and other seismic attributes of the main target layer of the omnidirectional seismic data body of the area to be measured, the fracture distribution characteristics of the area to be measured are defined, the fracture distribution characteristics mainly comprise the fracture scale and trend, the main fault of the area to be measured is determined according to the size of the scale, and then the trend of the main fault of the area to be measured is defined, so that a basis is provided for determining the number of the azimuth. For example, as shown in FIG. 2, it can be seen that the dominant fault in the region is a northeast distribution and the individual faults are northeast distributions.
In particular, in order to further improve the scientificity of dividing the data of the pre-stack trace set into azimuth, in this embodiment, as shown in table 1 below, the seismic acquisition parameters may be the bin size, coverage frequency, aspect ratio, and the like of the seismic acquisition, so as to provide a basis for determining the azimuth number.
TABLE 1
In particular, in order to further improve the detection precision, in this embodiment, the number of azimuth can be determined by considering the coverage times of seismic acquisition on the basis of the trend of the main fault, so that the number of azimuth can be determined by referring to the coverage times, for example, the coverage times in each azimuth can be not less than 60 times, so that the number of azimuth data for dividing the seismic pre-stack gather data can be determined, and the coverage times and the regional structural characteristics of the seismic pre-stack gather data are fully considered in the process of dividing the azimuth, so that the azimuth dividing scheme can be determined on the basis of scientific demonstration, the accuracy of the azimuth optimization can be realized according to the demonstration, and the effectiveness and the feasibility of oil gas detection can be improved.
In specific implementation, after determining the number of azimuth, the angle covered by each azimuth may be determined, for example, as shown in fig. 4, the angle covered by a certain azimuth may be 28 degrees, the angle covered by a certain azimuth may be 34 degrees, and the like, and different azimuth may cover different angles.
In specific implementation, after the number of azimuth is determined, the azimuth can be divided into the seismic pre-stack gather data of the area to be detected according to the number of azimuth, the direction parallel to the trend of the main fault can be determined to be 0 degree, then the angle of each azimuth is divided into the seismic pre-stack gather data of the area to be detected according to the number of azimuth, the seismic pre-stack gather data corresponding to each azimuth can be overlapped, the post-stack azimuth data body corresponding to each azimuth is obtained, and then the post-stack azimuth data bodies corresponding to a plurality of azimuth are obtained.
In particular, in order to implement high-precision oil and gas detection, in this embodiment, a dominant azimuth may be optimized by the following steps, so as to perform oil and gas detection based on the dominant azimuth, for example, determining a post-stack split azimuth data body with minimum anisotropy from post-stack split azimuth data bodies corresponding to a plurality of azimuths, including:
performing amplitude analysis and fault characterization on the post-stack azimuth data body corresponding to each azimuth;
AVO (amplitude versus offset) analysis is carried out on seismic prestack gather data of a region to be detected;
determining a post-stack orientation data volume satisfying the following condition among the plurality of orientations as a post-stack orientation data volume having minimum anisotropy:
the amplitude is maximum, the fault depiction is weakest, and the AVO analysis oil gas response characteristics meet the requirements.
In specific implementation, for the post-stack azimuth data volume corresponding to each azimuth, the amplitude intensity analysis is performed, as shown in fig. 4, the amplitude is strongest in the direction parallel to the fault azimuth, and the amplitude is weakest in the direction perpendicular to the fault azimuth, so that the post-stack azimuth data volume corresponding to the fault azimuth is more suitable for oil and gas detection.
In the specific implementation, fault depiction is carried out on the post-stack azimuth data body corresponding to each azimuth, the fault effect is best when the azimuth perpendicular to the fault is distinguished, and the anisotropy is strongest; the orientation fault description parallel to the fault is not obvious, and the anisotropy is weakest, so that the post-stack orientation data volume corresponding to the orientation parallel to the fault is more suitable for carrying out oil gas detection.
In specific implementation, AVO analysis is carried out on seismic pre-stack trace set data of a region to be detected, the method also shows that the azimuth oil gas response characteristics parallel to faults are most obvious, by combining the results of amplitude analysis and fault characterization on post-stack azimuth data bodies corresponding to all azimuths, the post-stack azimuth data bodies corresponding to the azimuth with the maximum amplitude, the weakest fault characterization and the AVO analysis oil gas response characteristics meeting requirements can be determined to be the smallest in anisotropy, and the post-stack azimuth data body with the smallest determined anisotropy is the post-stack azimuth data body corresponding to the dominant azimuth, so that oil gas detection is carried out based on the post-stack azimuth data body corresponding to the dominant azimuth.
In order to further reduce the detected data amount, in this embodiment, the amplitude analysis and tomographic characterization are performed on the post-stack azimuth data volume corresponding to each azimuth, including:
performing horizon interpretation on the post-stack azimuth data body corresponding to each azimuth; specifically, the horizon of the oil gas detection of the deep carbonate reservoir, namely the target horizon, can be determined based on horizon interpretation.
And carrying out amplitude analysis and fault depiction on the target horizon of the post-stack azimuth data body corresponding to each azimuth.
In specific implementation, after the post-stack azimuth data body of the dominant azimuth with the minimum anisotropy is determined, the frequency band of low-frequency resonance and high-frequency attenuation in the post-stack azimuth data body of the dominant azimuth is determined through spectrum analysis based on the Biot biphase medium theory, and then oil gas detection is carried out. The anisotropic influence is reduced in the oil gas detection in the dominant azimuth, the oil gas detection result is more accurate, the effectiveness and the applicability of the oil gas detection method for the deep carbonate reservoir are verified, meanwhile, the method is high in calculation efficiency, the restriction of a known well is not needed, and the calculated oil gas detection is more accurate under the condition of less deep oil gas exploration logging data.
In specific implementation, the effectiveness and applicability of the deep carbonate reservoir oil-gas detection method are verified through the following steps.
1) Omnibearing fault prediction: extracting all-round seismic data bulk property of the region to be detected, and determining fracture distribution characteristics of the region to be detected, wherein the fracture distribution characteristics mainly comprise fracture dimensions and trends, for example, as shown in fig. 2, main faults in the region can be seen to be northeast distribution, and individual faults are respectively north-northeast distribution;
2) And (3) analyzing the seismic acquisition parameters: the coverage times of the seismic acquisition are mainly analyzed, a basis is provided for selecting the number of the sub-azimuth of the seismic data, for example, the coverage times are 300 times in the example;
3) Acquiring seismic trace set data: obtaining post-stack seismic data and pre-stack CRP gathers through seismic data amplitude preservation processing, carrying out corresponding processing on deep carbonate reservoirs to obtain seismic pre-stack gather data with high fidelity, high resolution and high signal to noise ratio, and providing guarantee for obtaining azimuth-dividing superposition data volumes as shown in figure 3;
4) Determining a seismic data azimuth scheme: according to the fracture spread characteristics of the step 1) and the analysis results of the seismic acquisition parameters of the step 2), according to the standard that the coverage times of each azimuth dividing data body is not less than 60 times, the area to be measured is suitable for dividing 5 azimuth, and the azimuth dividing angles are mainly parallel to the northeast faults and perpendicular to the northeast faults,
5) Obtaining azimuth data volume: dividing the seismic pre-stack gather data in the step 3) into 5 pieces according to the azimuth dividing preferred scheme in the step 4), respectively superposing the 5 pieces of data divided into corresponding pieces to obtain 5 post-stack azimuth dividing data bodies, and providing preparation data for the azimuth selection of the oil and gas detection advantages;
6) Seismic horizon interpretation: after the 5 post-stack azimuth data bodies are obtained in the step 5), respectively carrying out main purpose layer-by-layer azimuth interpretation on the 5 post-stack azimuth data bodies, and providing reliable data for subsequent attribute analysis;
7) Azimuth amplitude analysis: respectively carrying out amplitude analysis and comparison along a target layer on the 5 post-stack azimuth data volumes obtained in the step 5), wherein as shown in fig. 4, the amplitude is strongest in the azimuth of the fault parallel to the northeast direction, the amplitude is weakest in the azimuth of the fault perpendicular to the northeast direction, and the post-stack azimuth data volumes parallel to the azimuth of the northeast direction are more suitable for carrying out oil gas detection;
8) Azimuth anisotropy analysis: carrying out fault detection according to coherence attribute along a target layer on the 5 post-stack azimuth data volumes obtained in the step 5), wherein as shown in fig. 5, the fault is best depicted in the azimuth vertical to the northeast fault, and the anisotropy is strongest; the azimuth fault description parallel to the northeast fault is not obvious, the anisotropy is weakest, and the post-stack azimuth data volume of the azimuth is more suitable for carrying out oil gas detection;
9) Split azimuth prestack AVO analysis: carrying out AVO analysis of different orientations by utilizing the pre-stack trace set obtained in the step 3), as shown in fig. 6 (a) and 6 (b), the orientation oil gas response characteristics parallel to the northeast fault are most obvious, and combining the conclusion of the step 7) and the step 8), the orientation parallel to the northeast fault can be determined to be the dominant orientation of oil gas detection, namely the anisotropy of the post-stack split orientation data volume parallel to the northeast fault orientation is minimum;
10 Seismic spectrum analysis: carrying out spectrum analysis of a dry well, a differential gas well and an oil-gas well on an all-dimensional post-stack seismic data body and an advantageous-dimensional post-stack azimuth data body respectively, determining azimuth data bodies with most obvious low-frequency resonance and high-frequency attenuation, as shown in fig. 7 (a) and 7 (b) (the fig. 7 (a) is an all-dimensional seismic spectrum analysis diagram, and the fig. 7 (b) is an advantageous-dimensional seismic spectrum analysis diagram), finally determining azimuth data bodies parallel to a northeast fault as optimal azimuth data bodies for oil-gas detection, and further determining the spectrum distribution range of high-frequency attenuation and low-frequency resonance of an oil-gas-containing interval at the same time, so as to provide basis for oil-gas detection;
11 Omnibearing oil gas detection: carrying out oil gas detection on the omnibearing post-stack seismic data body through the frequency spectrum change range of low-frequency resonance and high-frequency attenuation of the oil-gas layer section of the omnibearing post-stack seismic data body determined in the step 10), and comparing the subsequent oil gas detection result with the dominant azimuth oil gas detection result;
12 Dominant azimuth oil gas detection: developing dominant azimuth oil gas detection through the spectrum variation range of low-frequency resonance and high-frequency attenuation of the oil-gas reservoir section of the dominant azimuth data body determined in the step 10), as shown in fig. 9;
13 Oil gas detection result analysis: by comparing the omni-directional hydrocarbon detection result in step 11) with the optimal directional hydrocarbon detection result in step 12), as shown in fig. 8 ((a) is an omni-directional hydrocarbon detection plan view, (b) is a dominant directional hydrocarbon detection plan view), the advantages of the sub-directional hydrocarbon detection result based on the above-mentioned deep carbonate hydrocarbon reservoir hydrocarbon detection method are further demonstrated, and the prediction of the favorable hydrocarbon zone of deep carbonate is finally completed. From the known gas producing well, the influence of anisotropy is reduced by the dominant azimuth oil gas detection, the oil gas detection result is more accurate, and the effectiveness and applicability of the deep carbonate oil gas reservoir oil gas detection method are verified.
In this embodiment, a computer device is provided, as shown in fig. 10, including a memory 1002, a processor 1004, and a computer program stored on the memory and capable of running on the processor, where the processor implements any of the above deep carbonate reservoir hydrocarbon detection methods when executing the computer program.
In particular, the computer device may be a computer terminal, a server or similar computing means.
In this embodiment, a computer-readable storage medium storing a computer program for executing any of the above deep carbonate reservoir hydrocarbon detection methods is provided.
In particular, computer-readable storage media, including both permanent and non-permanent, removable and non-removable media, may be used to implement information storage by any method or technology. The information may be computer readable instructions, data structures, modules of a program, or other data. Examples of computer-readable storage media include, but are not limited to, phase-change memory (PRAM), static Random Access Memory (SRAM), dynamic Random Access Memory (DRAM), other types of Random Access Memory (RAM), read Only Memory (ROM), electrically Erasable Programmable Read Only Memory (EEPROM), flash memory or other memory technology, compact disc read only memory (CD-ROM), digital Versatile Disks (DVD) or other optical storage, magnetic cassettes, magnetic tape magnetic disk storage or other magnetic storage devices, or any other non-transmission medium, which can be used to store information that can be accessed by a computing device. Computer-readable storage media, as defined herein, does not include transitory computer-readable media (transmission media), such as modulated data signals and carrier waves.
Based on the same inventive concept, the embodiment of the application also provides a deep carbonate reservoir oil-gas detection device, as described in the following embodiment. Because the principle of the deep carbonate reservoir oil gas detection device for solving the problem is similar to that of the deep carbonate reservoir oil gas detection method, the implementation of the deep carbonate reservoir oil gas detection device can be referred to the implementation of the deep carbonate reservoir oil gas detection method, and the repetition is omitted. As used below, the term "unit" or "module" may be a combination of software and/or hardware that implements the intended function. While the means described in the following embodiments are preferably implemented in software, implementation in hardware, or a combination of software and hardware, is also possible and contemplated.
Fig. 11 is a block diagram of a deep carbonate reservoir oil and gas detection apparatus according to an embodiment of the present application, as shown in fig. 11, the apparatus includes:
the azimuth number determining module 1102 is configured to determine the azimuth number according to the main fault distribution characteristics and the seismic acquisition parameters of the area to be detected;
the azimuth dividing module 1104 is configured to divide azimuth for the seismic pre-stack gather data of the area to be measured according to the number of azimuth, and obtain post-stack azimuth data volumes corresponding to a plurality of azimuth according to the seismic pre-stack gather data corresponding to a plurality of azimuth, where the plurality of azimuth includes azimuth of parallel main fault strike and azimuth of perpendicular main fault strike;
the azimuth screening module 1106 is configured to determine, from the post-stack azimuth data volumes corresponding to the plurality of azimuth, a post-stack azimuth data volume with minimum anisotropy;
the detection module 1108 is configured to perform oil and gas detection based on the determined post-stack azimuth data body.
In one embodiment, the azimuth number determining module is specifically configured to determine the azimuth number according to the trend of the main fault of the area to be detected and the coverage frequency in each azimuth, where the coverage frequency in each azimuth makes the data corresponding to each azimuth meet the preset signal-to-noise ratio requirement.
In one embodiment, the position screening module includes:
the data analysis unit is used for carrying out amplitude analysis and fault depiction on the post-stack azimuth data body corresponding to each azimuth;
the AVO analysis unit is used for carrying out AVO analysis on the seismic prestack gather data of the area to be detected;
an azimuth screening unit configured to determine, as a post-stack azimuth data volume having minimum anisotropy, a post-stack azimuth data volume satisfying the following condition among a plurality of azimuth:
the amplitude is maximum, the fault depiction is weakest, and the AVO analysis oil gas response characteristics meet the requirements.
In one embodiment, the data analysis unit is specifically configured to perform horizon interpretation for the post-stack sub-azimuth data body corresponding to each azimuth, and perform amplitude analysis and fault characterization for the target horizon of the post-stack sub-azimuth data body corresponding to each azimuth.
The embodiment of the application realizes the following technical effects: the method comprises the steps of determining the number of azimuth according to the main fault distribution characteristics and the earthquake acquisition parameters of a region to be detected, dividing the earthquake pre-stack gather data of the region to be detected according to the number of azimuth, and obtaining post-stack azimuth data bodies corresponding to the azimuth according to the earthquake pre-stack gather data corresponding to the azimuth, wherein the azimuth comprises azimuth parallel to the main fault trend and azimuth perpendicular to the main fault trend, dividing the earthquake pre-stack gather data of the region to be detected into azimuth based on the main fault distribution characteristics and the earthquake acquisition parameters of the region to be detected, fully considering the earthquake acquisition parameters and the region construction characteristics of the pre-stack gather, and enabling the division of azimuth to be scientific; meanwhile, in the post-stack azimuth data bodies corresponding to a plurality of azimuth, the post-stack azimuth data body with the minimum anisotropism is determined, and the oil gas detection is carried out based on the determined post-stack azimuth data body.
It will be apparent to those skilled in the art that the modules or steps of the embodiments of the application described above may be implemented in a general purpose computing device, they may be concentrated on a single computing device, or distributed across a network of computing devices, they may alternatively be implemented in program code executable by computing devices, so that they may be stored in a storage device for execution by computing devices, and in some cases, the steps shown or described may be performed in a different order than what is shown or described, or they may be separately fabricated into individual integrated circuit modules, or a plurality of modules or steps in them may be fabricated into a single integrated circuit module. Thus, embodiments of the application are not limited to any specific combination of hardware and software.
The above description is only of the preferred embodiments of the present application and is not intended to limit the present application, and various modifications and variations can be made to the embodiments of the present application by those skilled in the art. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the present application should be included in the protection scope of the present application.
Claims (8)
1. A method for detecting oil and gas in a deep carbonate reservoir, comprising the steps of:
determining the number of azimuth according to the distribution characteristics of the main fault and the seismic acquisition parameters of the region to be detected;
dividing the data of the seismic pre-stack gathers of the area to be detected according to the number of the orientations, and acquiring post-stack orientation data bodies corresponding to the orientations according to the data of the seismic pre-stack gathers corresponding to the orientations, wherein the orientations comprise orientations of parallel main fault strike and orientations of vertical main fault strike;
determining the post-stack azimuth data body with the minimum anisotropism in the post-stack azimuth data bodies corresponding to the plurality of azimuths;
oil gas detection is carried out based on the determined post-stack azimuth data body;
determining the post-stack azimuth data body with the minimum anisotropism in the post-stack azimuth data bodies corresponding to the plurality of azimuths, wherein the post-stack azimuth data body comprises the following steps:
performing amplitude analysis and fault characterization on the post-stack azimuth data body corresponding to each azimuth;
AVO analysis is carried out on the seismic prestack gather data of the area to be detected;
determining a post-stack orientation data volume satisfying the following condition among the plurality of orientations as a post-stack orientation data volume having minimum anisotropy:
the amplitude is maximum, the fault depiction is weakest, and the AVO analysis oil gas response characteristics meet the requirements.
2. The method for detecting hydrocarbon in a deep carbonate hydrocarbon reservoir according to claim 1, wherein determining the number of azimuth according to the main fault distribution characteristics and the seismic acquisition parameters of the area to be detected comprises:
and determining the number of the azimuth according to the trend of the main fault of the area to be detected and the coverage times in each azimuth, wherein the coverage times in each azimuth enable the data corresponding to each azimuth to meet the preset signal-to-noise ratio requirement.
3. The deep carbonate reservoir hydrocarbon testing method of claim 1, wherein for each azimuth corresponding post stack azimuth data volume, performing amplitude analysis and fault characterization comprises:
performing horizon interpretation on the post-stack azimuth data body corresponding to each azimuth;
and carrying out amplitude analysis and fault depiction on the target horizon of the post-stack azimuth data body corresponding to each azimuth.
4. A deep carbonate reservoir oil gas detection device, comprising:
the azimuth number determining module is used for determining the azimuth number according to the main fault distribution characteristics and the seismic acquisition parameters of the region to be detected;
the azimuth dividing module is used for dividing the earthquake pre-stack gather data of the area to be detected according to the azimuth number, and acquiring post-stack azimuth dividing data bodies corresponding to a plurality of azimuths according to the earthquake pre-stack gather data corresponding to the azimuths, wherein the azimuths comprise azimuth of parallel main fault strike and azimuth of vertical main fault strike;
the azimuth screening module is used for determining the post-stack azimuth data body with the minimum anisotropism from the post-stack azimuth data bodies corresponding to the plurality of azimuths;
the detection module is used for carrying out oil gas detection based on the determined post-stack azimuth data body;
the azimuth screening module comprises:
the data analysis unit is used for carrying out amplitude analysis and fault depiction on the post-stack azimuth data body corresponding to each azimuth;
the AVO analysis unit is used for carrying out AVO analysis on the seismic prestack gather data of the area to be detected;
an azimuth screening unit configured to determine, as a post-stack azimuth data volume having minimum anisotropy, a post-stack azimuth data volume satisfying the following condition among a plurality of azimuth:
the amplitude is maximum, the fault depiction is weakest, and the AVO analysis oil gas response characteristics meet the requirements.
5. The deep carbonate reservoir hydrocarbon detection device of claim 4, wherein the azimuth number determination module is specifically configured to determine the azimuth number according to the trend of the main fault of the area to be detected and the coverage number in each azimuth, where the coverage number in each azimuth makes the data corresponding to each azimuth meet the preset signal-to-noise ratio requirement.
6. The deep carbonate reservoir hydrocarbon detection device of claim 4, wherein the data analysis unit is specifically configured to perform horizon interpretation for a post-stack sub-azimuth data volume corresponding to each azimuth, and perform amplitude analysis and fault characterization for a target horizon of the post-stack sub-azimuth data volume corresponding to each azimuth.
7. A computer device comprising a memory, a processor and a computer program stored on the memory and executable on the processor, characterized in that the processor implements the deep carbonate reservoir hydrocarbon detection method of any one of claims 1 to 3 when the computer program is executed.
8. A computer-readable storage medium, characterized in that the computer-readable storage medium stores a computer program for executing the deep carbonate reservoir hydrocarbon detection method of any one of claims 1 to 3.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010440587.3A CN113700466B (en) | 2020-05-22 | 2020-05-22 | Method, device, equipment and medium for detecting oil gas in deep carbonate reservoir |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010440587.3A CN113700466B (en) | 2020-05-22 | 2020-05-22 | Method, device, equipment and medium for detecting oil gas in deep carbonate reservoir |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113700466A CN113700466A (en) | 2021-11-26 |
| CN113700466B true CN113700466B (en) | 2023-10-31 |
Family
ID=78646142
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010440587.3A Active CN113700466B (en) | 2020-05-22 | 2020-05-22 | Method, device, equipment and medium for detecting oil gas in deep carbonate reservoir |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113700466B (en) |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103527184A (en) * | 2013-10-28 | 2014-01-22 | 北京大学 | Method and system for predicting dolomite reservoir |
| CN104977618A (en) * | 2014-04-09 | 2015-10-14 | 中国石油集团东方地球物理勘探有限责任公司 | Method for evaluating shale gas reservoir and finding dessert area |
| CN105425292A (en) * | 2015-11-20 | 2016-03-23 | 中国石油天然气集团公司 | Oil and gas prediction method and oil and gas prediction device |
| WO2016041189A1 (en) * | 2014-09-19 | 2016-03-24 | 杨顺伟 | Method for evaluating shale gas reservoir and seeking desert area |
| CN106842313A (en) * | 2015-12-04 | 2017-06-13 | 中国石油化工股份有限公司 | Anisotropic parameters inversion method based on orientation earthquake data before superposition |
-
2020
- 2020-05-22 CN CN202010440587.3A patent/CN113700466B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103527184A (en) * | 2013-10-28 | 2014-01-22 | 北京大学 | Method and system for predicting dolomite reservoir |
| CN104977618A (en) * | 2014-04-09 | 2015-10-14 | 中国石油集团东方地球物理勘探有限责任公司 | Method for evaluating shale gas reservoir and finding dessert area |
| WO2016041189A1 (en) * | 2014-09-19 | 2016-03-24 | 杨顺伟 | Method for evaluating shale gas reservoir and seeking desert area |
| CN105425292A (en) * | 2015-11-20 | 2016-03-23 | 中国石油天然气集团公司 | Oil and gas prediction method and oil and gas prediction device |
| CN106842313A (en) * | 2015-12-04 | 2017-06-13 | 中国石油化工股份有限公司 | Anisotropic parameters inversion method based on orientation earthquake data before superposition |
Non-Patent Citations (1)
| Title |
|---|
| 海拉尔盆地乌东地区裂缝储层综合分析;刘俊峰;孟小红;王建民;陈守田;安朝晖;;石油地球物理勘探(S1);第181-184页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113700466A (en) | 2021-11-26 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Wang et al. | Seismic sensor misorientation measurement using P‐wave particle motion: An application to the NECsaids array | |
| US7539578B2 (en) | VH signal integration measure for seismic data | |
| D’Alessandro et al. | Seismic network evaluation through simulation: An application to the Italian National Seismic Network | |
| CN112883564B (en) | Water body temperature prediction method and prediction system based on random forest | |
| CN106556861B (en) | A kind of azimuthal AVO inversion method based on Omnibearing earthquake auto data | |
| Zhu et al. | Determining the orientation of ocean‐bottom seismometers on the seafloor and correcting for polarity flipping via polarization analysis and waveform modeling | |
| CN108646292B (en) | Crack density prediction method, crack density prediction device and computer storage medium | |
| CN105445814B (en) | The azimuthal anisotropy attribute difference crack prediction method of marine narrow azimuth data | |
| Yang et al. | Determination of the local magnitudes of small earthquakes using a dense seismic array in the Changning–Zhaotong Shale Gas Field, Southern Sichuan Basin | |
| Malagnini et al. | Estimating absolute site effects | |
| Sun et al. | P-wave fracture prediction algorithm using data with limited azimuthal distribution | |
| Wang et al. | Distribution of Rayleigh wave microseisms constrained by multiple seismic arrays | |
| O’Rourke et al. | A spectrogram‐based method of Rg detection for explosion monitoring | |
| Bayer et al. | Tracking unilateral earthquake rupture by P-wave polarization analysis | |
| CN105445787A (en) | Crack prediction method for preferred orientation daughter coherence | |
| EP2703847A1 (en) | Method of dynamic seismic prospecting | |
| CN113700466B (en) | Method, device, equipment and medium for detecting oil gas in deep carbonate reservoir | |
| CN108562936A (en) | A kind of crack prediction method, system, storage medium and terminal | |
| Xie et al. | Evaluating Global Tomography Models With Antipodal Ambient Noise Cross‐Correlation Functions | |
| Ikawa et al. | A case study of seismic anisotropy analysis on multi-azimuth OBC seismic data in an offshore carbonate oil field, UAE | |
| Luo et al. | Abnormal transmission attenuation and its impact on seismic-fracture prediction—A physical modeling study | |
| Cristiano et al. | Teleseismic P-wave polarization analysis at the Gräfenberg array | |
| Karamzadeh et al. | Small-aperture array as a tool to monitor fluid injection-and extraction-induced microseismicity: applications and recommendations | |
| Rentsch et al. | Migration-based location of seismicity recorded with an array installed in the main hole of the San Andreas Fault Observatory at Depth (SAFOD) | |
| Albuquerque et al. | Estimating seismometer component orientation of the Brazilian seismographic network using teleseismic P-wave particle motion analysis and directional statistics |
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 |