CN112505761A - Reservoir gas-containing property detection method and device - Google Patents
Reservoir gas-containing property detection method and device Download PDFInfo
- Publication number
- CN112505761A CN112505761A CN202011175054.3A CN202011175054A CN112505761A CN 112505761 A CN112505761 A CN 112505761A CN 202011175054 A CN202011175054 A CN 202011175054A CN 112505761 A CN112505761 A CN 112505761A
- Authority
- CN
- China
- Prior art keywords
- gas
- containing sandstone
- prediction factor
- reservoir
- sandstone
- 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.)
- Granted
Links
- 238000001514 detection method Methods 0.000 title claims abstract description 32
- 238000004458 analytical method Methods 0.000 claims abstract description 53
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 53
- 238000000034 method Methods 0.000 claims abstract description 41
- 239000011435 rock Substances 0.000 claims abstract description 37
- 238000004364 calculation method Methods 0.000 claims abstract description 25
- 239000004576 sand Substances 0.000 claims description 23
- 238000010276 construction Methods 0.000 claims description 6
- 238000011161 development Methods 0.000 abstract description 3
- 238000011160 research Methods 0.000 description 16
- 238000010586 diagram Methods 0.000 description 5
- 238000010206 sensitivity analysis Methods 0.000 description 4
- 239000012530 fluid Substances 0.000 description 3
- 238000012360 testing method Methods 0.000 description 3
- 230000035945 sensitivity Effects 0.000 description 2
- 239000004215 Carbon black (E152) Substances 0.000 description 1
- 238000007418 data mining Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 229930195733 hydrocarbon Natural products 0.000 description 1
- 150000002430 hydrocarbons Chemical class 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000002194 synthesizing effect Effects 0.000 description 1
Images
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/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/40—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging
-
- 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/61—Analysis by combining or comparing a seismic data set with other data
- G01V2210/616—Data from specific type of measurement
- G01V2210/6161—Seismic or acoustic, e.g. land or sea measurements
-
- 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/61—Analysis by combining or comparing a seismic data set with other data
- G01V2210/616—Data from specific type of measurement
- G01V2210/6169—Data from specific type of measurement using well-logging
-
- 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/62—Physical property of subsurface
- G01V2210/624—Reservoir parameters
Landscapes
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Remote Sensing (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 provides a reservoir gas-containing property detection method and device. The reservoir gas-containing detection method comprises the following steps: (1) acquiring logging data and seismic data; (2) constructing a rock physical analysis quantity version based on the logging data; (3) fitting an aqueous sandstone linear trend line based on the petrophysical analysis quantity plate; (4) constructing a gas-containing sandstone prediction factor based on the rock physical analysis quantity version and the water-containing sandstone linear trend line; (5) substituting the seismic data into a calculation formula of the gas-containing sandstone prediction factor to obtain a gas-containing sandstone prediction factor seismic data volume; (6) and performing reservoir gas content detection based on the gas-containing sandstone prediction factor seismic data volume. The method can reduce the multi-solution of reservoir prediction and improve the precision of exploration and development.
Description
Technical Field
The invention relates to the field of geophysical exploration of petroleum, in particular to a method and a device for predicting a high-porosity hydrocarbon-containing sandstone reservoir.
Background
In the oil exploration process, prestack inversion and rock physics analysis technology are applied more and more, good application effects are obtained in different research blocks, the reservoir identification factor constructed on the basis of attributes such as longitudinal wave impedance, transverse wave impedance and density can realize the function of distinguishing sandstone reservoirs and mudstone cover layers to a certain extent, and the classical reservoir identification factor comprises: the longitudinal and transverse wave velocity ratio, the Poisson impedance, the Russell fluid factors and the like have different applicability to different identification factors in different work areas and basins, and sensitivity analysis is required in the practical application process. Under the general condition, identification factors with strong sensitivity to reservoir lithology and fluid need to be introduced into undetermined coefficients in the construction process, and the introduction of the undetermined coefficients enhances the applicability of attribute factors on one hand, and increases the artificial interference of reservoir prediction and the multi-solution of prediction results on the other hand, and further aggravates the uncontrollable property of the prediction results for the lack of rock physical basic data of a research block with low exploration degree.
An effective solution to the above problems is to determine reservoir identification factors in a data-driven manner, construct attribute factors with work area applicability based on collected data, fully utilize empirical data fitting and trend analysis, perform data mining deeply, and summarize local laws. However, the existing reservoir gas-containing property detection method still has the defects of low precision and uncontrollable prediction result. Therefore, a method for detecting the gas content of the reservoir with higher precision is needed to be found.
Disclosure of Invention
In order to solve the above problems, an object of the present invention is to provide a reservoir gas-bearing property detection method, which can reduce the ambiguity of reservoir prediction and improve the accuracy of exploration and development.
Another object of the present invention is to provide a reservoir gas-containing property detection device.
In order to achieve the above object, in one aspect, the present invention provides a method for detecting reservoir gas-bearing property, wherein the method comprises: (1) acquiring logging data and seismic data; (2) constructing a rock physical analysis quantity version based on the logging data; (3) fitting an aqueous sandstone linear trend line based on the petrophysical analysis quantity plate; (4) constructing a gas-containing sandstone prediction factor based on the rock physical analysis quantity version and the water-containing sandstone linear trend line; (5) substituting the seismic data into a calculation formula of the gas-containing sandstone prediction factor to obtain a gas-containing sandstone prediction factor seismic data volume; (6) and performing reservoir gas content detection based on the gas-containing sandstone prediction factor seismic data volume.
According to some embodiments of the invention, the logging data in step (1) comprises: longitudinal wave velocity curve, transverse wave velocity curve, density curve and well logging oil gas water interpretation result.
According to some embodiments of the invention, the seismic data in step (1) comprises: and obtaining longitudinal wave velocity, transverse wave velocity and density seismic data volume by pre-stack seismic inversion.
According to some embodiments of the present invention, in step (2), the method for constructing the petrophysical analysis quantitative plate based on the log data comprises: constructing a rock physical analysis quantity version based on a longitudinal wave velocity curve, a transverse wave velocity curve, a density curve and a logging oil-gas-water interpretation result in logging data, wherein rocks areThe abscissa of the physical analysis measuring plate is X, the ordinate is Y, and X is rho/vs，Y＝ρvpRho is density in g/cm3，vsIs the transverse wave velocity in Km/s, vpThe unit is the velocity of the longitudinal wave, Km/s.
According to some specific embodiments of the present invention, in the step (3), the specific method for fitting the linear trend line of the aqueous sandstone based on the petrophysical analysis quantitative plate comprises: in the rock physical analysis quantity version, selecting data points of a water-containing sandstone storage area of a target interval to perform linear fitting to obtain a water-containing sandstone linear trend line, wherein the trend line formula obtained by fitting is as follows:
Y＝kX+d
wherein X is rho/vs，Y＝ρvp，
k is the slope, d is the intercept, ρ is the density in g/cm3，vsIs the transverse wave velocity in Km/s, vpThe unit is the velocity of the longitudinal wave, Km/s.
According to some specific embodiments of the present invention, in the step (4), the specific method for constructing the gas sand predicting factor based on the petrophysical analysis quantitative plate and the water-containing sandstone linear trend line comprises: constructing a calculation formula of the gas-containing sandstone prediction factor based on the compressional wave velocity, the shear wave velocity, the density curve and the linear trend line of the water-containing sandstone obtained by fitting, wherein the calculation formula is as follows:
wherein F is a gas-containing sandstone prediction factor, k is a slope, rho is a density, and the unit is g/cm3，vsIs the transverse wave velocity in Km/s, vpThe unit is the velocity of the longitudinal wave, Km/s.
According to some specific embodiments of the present invention, the specific method for reservoir gas bearing detection based on the gas sand predictor seismic data volume comprises: and predicting and researching the distribution rule of the gas reservoir in the work area based on the gas-containing sandstone prediction factor seismic data volume, wherein the distribution rule corresponds to the distribution range of the gas-containing sandstone reservoir when the gas-containing sandstone prediction factor is smaller than a gas-containing sandstone prediction threshold value, corresponds to the distribution range of mudstone when the gas-containing sandstone prediction factor is larger than the mudstone prediction threshold value, and corresponds to the distribution range of the water-containing sandstone reservoir when the gas-containing sandstone prediction factor is positioned between the gas-containing sandstone prediction factor and the mudstone.
In another aspect, the present invention further provides a device for detecting reservoir gas content, wherein the device is used to implement the method for detecting reservoir gas content as claimed above, and the device includes: the input unit is used for inputting logging data and seismic data; the analysis unit is used for constructing a rock physical analysis quantity version based on the logging data; the fitting unit is used for fitting the linear trend line of the water-containing sandstone based on the rock physical analysis quantity plate; the device comprises a construction unit, a calculation unit and a calculation unit, wherein the construction unit is used for constructing a gas-containing sandstone prediction factor based on a rock physical analysis quantity version and a water-containing sandstone linear trend line; the calculation unit is used for substituting the seismic data into a gas-containing sandstone prediction factor calculation formula to obtain a gas-containing sandstone prediction factor seismic data volume; and the detection unit is used for detecting the gas content of the reservoir based on the gas-containing sandstone prediction factor seismic data body.
According to some embodiments of the invention, the well log data comprises: longitudinal wave velocity curve, transverse wave velocity curve, density curve and well logging oil gas water interpretation result.
According to some embodiments of the invention, the seismic data comprises: and obtaining longitudinal wave velocity, transverse wave velocity and density seismic data volume by pre-stack seismic inversion.
According to some embodiments of the present invention, the method for constructing the petrophysical analysis quantitative plate based on the well log data comprises: constructing a rock physical analysis quantity version based on a longitudinal wave velocity curve, a transverse wave velocity curve, a density curve and a logging oil-gas-water interpretation result in logging data, wherein the abscissa of the rock physical analysis quantity version is X, the ordinate is Y, and X is rho/vs，Y＝ρvpRho is density in g/cm3，vsIs the transverse wave velocity in Km/s, vpThe unit is the velocity of the longitudinal wave, Km/s.
According to some embodiments of the invention, the petrophysical analysis metric-based simulationThe specific method for synthesizing the linear trend line of the aqueous sandstone comprises the following steps: in the rock physical analysis quantity version, selecting data points of a water-containing sandstone storage area of a target interval to perform linear fitting to obtain a water-containing sandstone linear trend line, wherein the trend line formula obtained by fitting is as follows: y ═ kX + d, where X ═ ρ/vs，Y＝ρvpK is the slope, d is the intercept, ρ is the density in g/cm3，vsIs the transverse wave velocity in Km/s, vpThe unit is the velocity of the longitudinal wave, Km/s.
According to some specific embodiments of the present invention, the specific method for constructing the gas sand prediction factor based on the petrophysical analysis quantitative plate and the water-containing sand linear trend line comprises the following steps: constructing a calculation formula of the gas-containing sandstone prediction factor based on the compressional wave velocity, the shear wave velocity, the density curve and the linear trend line of the water-containing sandstone obtained by fitting, wherein the calculation formula is as follows:
wherein F is a gas-containing sandstone prediction factor, k is a slope, rho is a density, and the unit is g/cm3，vsIs the transverse wave velocity in Km/s, vpThe unit is the velocity of the longitudinal wave, Km/s.
According to some specific embodiments of the present invention, the specific method for reservoir gas bearing detection based on the gas sand predictor seismic data volume comprises: and predicting and researching the distribution rule of the gas reservoir in the work area based on the gas-containing sandstone prediction factor seismic data volume, wherein the distribution rule corresponds to the distribution range of the gas-containing sandstone reservoir when the gas-containing sandstone prediction factor is smaller than a gas-containing sandstone prediction threshold value, corresponds to the distribution range of mudstone when the gas-containing sandstone prediction factor is larger than the mudstone prediction threshold value, and corresponds to the distribution range of the water-containing sandstone reservoir when the gas-containing sandstone prediction factor is positioned between the gas-containing sandstone prediction factor and the mudstone.
The reservoir gas-containing detection method provided by the invention constructs a self-adaptive reservoir gas-containing detection attribute factor based on a rock physical analysis quantity version and elastic parameters of pre-stack inversion, in the implementation process, a linear trend line of a water-containing sandstone reservoir is fitted based on the rock physical quantity version, then a gas-containing sandstone prediction factor calculation formula is constructed by combining reservoir sensitivity analysis, and finally the gas-containing sandstone prediction factor seismic data is substituted into the pre-stack inversion elastic parameters to obtain a gas-containing sandstone prediction factor seismic data body of a research work area and perform reservoir gas-containing detection. The reservoir gas-containing detection method is based on a data-driven research idea, reduces interference of human factors, makes full use of effective information of the existing data as far as possible, reduces the multi-solution of reservoir prediction, and improves the precision of exploration and development.
Drawings
FIG. 1 is a flow chart of a reservoir gas bearing detection method of the present invention;
FIG. 2 is a block diagram of a reservoir gas bearing detection apparatus of the present invention;
FIG. 3 is a fitting analysis diagram of a rock physical quantity plate and a water-containing sandstone constructed in a target interval of a research work area in example 1 of the present invention;
FIG. 4 is a histogram of analysis of the sensitivity of the gas-bearing sandstone prediction factor in the target interval of the research work area in example 1 of the present invention;
FIG. 5 is a seismic profile of a gas-sand predicting factor through a well for a target interval of a research work area in example 1 of the present invention;
fig. 6 is a planar distribution diagram of gas sand sandstone prediction factor in the target interval of the research work area in the embodiment 1 of the invention.
Detailed Description
The technical solution in the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings in the embodiments of the present invention. It is to be understood that the described embodiments are merely exemplary of the invention, and not restrictive of the full scope of the invention. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
Example 1
The embodiment provides a reservoir gas-bearing property detection method (the specific flow is shown in figure 1 in detail), which comprises the following steps:
step 101: inputting well log data and seismic data
In practice, the input logging data includes: the longitudinal wave velocity curve, the transverse wave velocity curve, the density curve and the well logging oil-gas-water interpretation result, the input seismic data comprise: and obtaining longitudinal wave velocity, transverse wave velocity and density seismic data volume by pre-stack seismic inversion.
Step 102: method for constructing rock physical analysis quantity version based on logging data
In the implementation process, a rock physical analysis quantity plate is constructed based on an input logging curve and a logging oil-gas-water interpretation result, the abscissa of the quantity plate is X, the ordinate is Y, and X is rho/vs，Y＝ρvpRho is density in g/cm3,vsIs the transverse wave velocity in Km/s, vpThe unit is the velocity of the longitudinal wave, Km/s. Fig. 3 shows a rock physical quantity plate of a target interval of a research work area constructed in the implementation, and it can be seen from the figure that there are obvious differences in distribution areas of different lithologies and different fluid-containing reservoirs in the quantity plate, wherein mudstone (dark color dots) is mainly distributed in the upper right corner area of the quantity plate, gas sand (light color dots) is mainly distributed in the lower left corner area of the quantity plate, and water-containing sand (black dots, oval circles) is distributed in the narrow area between the two areas. The distribution differences of different types of reservoirs provide theoretical basis for constructing the gas-containing sandstone prediction factor, and how to specifically describe the distribution range of the gas-containing sandstone reservoir is the key point of the next research.
Step 103: fitting water-containing sandstone linear trend line based on rock physical analysis quantity plate
In the implementation process, in the rock physical quantity version, selecting data points of a water-containing sandstone storage zone of a target interval to perform linear fitting to obtain a water-containing sandstone linear trend line, wherein the trend line formula obtained by fitting is as follows:
Y＝kX+d
wherein X is rho/vs，Y＝ρvp(ii) a k is the slope, d is the intercept, ρ is the density in g/cm3，vsIs the transverse wave velocity in Km/s, vpThe unit is the velocity of the longitudinal wave, Km/s.
In this embodiment, as shown in fig. 3, a linear fitting is performed on the data points (in the elliptical circles) of the water-containing sandstone reservoir, the black solid line in the figure is a linear trend line obtained by the fitting, and a trend line formula obtained by the fitting is as follows:
Y＝-8.05*X+18.9
that is, k is-8.05, d is 18.9, and the fitting correlation coefficient is 0.84, indicating that the linear fitting relationship is good.
Step 104: and constructing a gas-containing sandstone prediction factor based on the rock physical analysis quantity version and the water-containing sandstone linear trend line.
In the implementation process, a calculation formula of the gas-containing sandstone prediction factor is constructed based on the input compressional wave velocity, shear wave velocity and density curve and the linear trend line of the water-containing sandstone obtained by fitting, and is as follows:
wherein F is a gas-containing sandstone prediction factor, k is a slope, rho is a density, and the unit is g/cm3，vsIs the transverse wave velocity in Km/s, vpThe unit is the velocity of the longitudinal wave, Km/s.
In this embodiment, the slope k obtained in step 103 is substituted into-8.05 to obtain a prediction factor applicable to the target interval of the study area, where the calculation formula is:
fig. 4 shows a gas sand prediction factor sensitivity analysis histogram, in which the abscissa is a gas sand prediction factor F, the ordinate is a percentage, and two black vertical solid lines in the histogram divide the entire region into three parts according to the numerical value of the gas sand prediction factor, wherein the left bright color region (i.e., F <0.245) corresponds to a gas sand reservoir, the middle black region (i.e., F0.245 ≦ 0.315) corresponds to a water sand reservoir, and the right dark color region (i.e., F >0.315) corresponds to mudstone.
Step 105: substituting the seismic data into the calculation formula of the gas-containing sandstone prediction factor to obtain a gas-containing sandstone prediction factor seismic data volume
In the implementation process, the input prestack inversion longitudinal wave velocity, transverse wave velocity and density seismic data volume is substituted into the gas-containing sandstone prediction factor calculation formula in the step 104, so that a gas-containing sandstone prediction factor seismic data volume in the range of the target zone of the research work area can be obtained, and a data basis is provided for the subsequent prediction of the gas-containing reservoir distribution range.
Step 106: reservoir gas-bearing property detection based on gas-bearing sandstone prediction factor seismic data volume
In the implementation process, the distribution rule of the gas reservoir in the research work area is predicted by using the gas sandstone prediction factor seismic data volume based on the prediction factor sensitivity analysis result in the step 104, and the prediction result is shown in fig. 5 and 6. Fig. 5 shows a seismic profile of the gas sand prediction factor through the well at the target interval of the research work area, wherein the dark area represents the gas sand prediction factor with low value F <0.245, namely the gas sand area, and the light area represents the gas sand prediction factor with high value F >0.315, namely the mudstone area. According to the diagram, two sets of gas-containing sandstone reservoirs are obtained by predicting that a W-2 well and a W-3 well correspond to each other in the solid line ellipse range, a W-1 well corresponds to each other in the dotted line ellipse range, the prediction result is completely consistent with the oil testing result of three wells in a target interval, the oil testing result shows that the W-1 well is a water well, the W-2 well and the W-3 well are high-yield gas wells, and the geological research shows that the W-1 well is positioned below an oil-water interface, so that the target that the trap structure high point is higher than the oil-water interface is preferentially considered when the structure zone is subsequently deployed. Fig. 6 shows a planar distribution diagram of the gas sand prediction factor of the target interval of the research work area, wherein the dark color area represents the gas sand prediction factor with low value F <0.245, namely the gas sand area, and the light color area represents the gas sand prediction factor with high value F >0.315, namely the mudstone area. The distribution range and the distribution rule of the gas-containing sandstone in the research area can be seen from the figure, the prediction results of the three wells in the figure are consistent with the oil testing conclusion, and meanwhile, the W-2 well and the W-3 well belong to two different structural units, and a mudstone separation zone exists between the two structural units, so that the two structural units are treated differently when oil-water system analysis is carried out. The effectiveness of the method is verified through analysis of the embodiment, and effective technical support is provided for reservoir gas-bearing detection.
Example 2
The embodiment also provides a reservoir gas content detection device for implementing the reservoir gas content detection method, and in detail referring to fig. 2, the device includes:
an input unit 201 for inputting logging data and seismic data; the logging data comprise a longitudinal wave velocity curve, a transverse wave velocity curve, a density curve and a logging oil-gas-water interpretation result; the seismic data includes: and (3) obtaining the longitudinal wave velocity, the transverse wave velocity and the density by pre-stack seismic inversion.
The analysis unit 202 is used for constructing a rock physical analysis quantity version based on the logging data; the specific method comprises the following steps: constructing a rock physical analysis quantity version based on a longitudinal wave velocity curve, a transverse wave velocity curve, a density curve and a logging oil-gas-water interpretation result in logging data, wherein the abscissa of the rock physical analysis quantity version is X, the ordinate is Y, and X is rho/vs，Y＝ρvpRho is density in g/cm3,vsIs the transverse wave velocity in Km/s, vpThe unit is the velocity of the longitudinal wave, Km/s.
The fitting unit 203 is used for fitting the linear trend line of the water-containing sandstone based on the rock physical analysis quantity plate; the specific method comprises the following steps: in the rock physical analysis quantity version, selecting data points of a water-containing sandstone storage area of a target interval to perform linear fitting to obtain a water-containing sandstone linear trend line, wherein the trend line formula obtained by fitting is as follows:
Y＝kX+d
wherein X is rho/vs，Y＝ρvp(ii) a k is the slope, d is the intercept, ρ is the density in g/cm3，vsIs the transverse wave velocity in Km/s, vpThe unit is the velocity of the longitudinal wave, Km/s.
The construction unit 204 is used for constructing the gas-containing sandstone prediction factor based on the petrophysical analysis quantity version and the water-containing sandstone linear trend line; the specific method comprises the following steps:
constructing a calculation formula of the gas-containing sandstone prediction factor based on the compressional wave velocity, the shear wave velocity, the density curve and the linear trend line of the water-containing sandstone obtained by fitting, wherein the calculation formula is as follows:
wherein F is a gas-containing sandstone prediction factor, k is a slope, rho is a density, and the unit is g/cm3，vsIs the transverse wave velocity in Km/s, vpThe unit is the velocity of the longitudinal wave, Km/s.
And the calculating unit 205 is used for substituting the seismic data into the gas-containing sandstone prediction factor calculation formula to obtain a gas-containing sandstone prediction factor seismic data volume.
The detection unit 206 is used for detecting the reservoir gas content based on the gas-containing sandstone prediction factor seismic data volume; the specific method comprises the following steps: and predicting and researching the distribution rule of the gas reservoir in the work area based on the gas-containing sandstone prediction factor seismic data volume, wherein the distribution rule corresponds to the distribution range of the gas-containing sandstone reservoir when the gas-containing sandstone prediction factor is smaller than a gas-containing sandstone prediction threshold value, corresponds to the distribution range of mudstone when the gas-containing sandstone prediction factor is larger than the mudstone prediction threshold value, and corresponds to the distribution range of the water-containing sandstone reservoir when the gas-containing sandstone prediction factor is positioned between the gas-containing sandstone prediction factor and the mudstone prediction threshold value.
The above-mentioned embodiments, objects, technical solutions and advantages of the present invention are further described in detail, it should be understood that the above-mentioned embodiments are only examples of the present invention, and are not intended to limit the scope of the present invention, and any modifications or equivalent substitutions made within the spirit and principle of the present invention should be covered by the claims of the present invention.
Claims (10)
1. A method of reservoir gas bearing detection, wherein the method comprises:
(1) acquiring logging data and seismic data;
(2) constructing a rock physical analysis quantity version based on the logging data;
(3) fitting an aqueous sandstone linear trend line based on the petrophysical analysis quantity plate;
(4) constructing a gas-containing sandstone prediction factor based on the rock physical analysis quantity version and the water-containing sandstone linear trend line;
(5) substituting the seismic data into a calculation formula of the gas-containing sandstone prediction factor to obtain a gas-containing sandstone prediction factor seismic data volume;
(6) and performing reservoir gas content detection based on the gas-containing sandstone prediction factor seismic data volume.
2. The method of claim 1, wherein the logging data in step (1) comprises: longitudinal wave velocity curve, transverse wave velocity curve, density curve and well logging oil gas water interpretation result.
3. The method of claim 1 or 2, wherein the seismic data in step (1) comprises: and (3) obtaining the longitudinal wave velocity, the transverse wave velocity and the density by pre-stack seismic inversion.
4. The method according to any one of claims 1-3, wherein in step (2), the specific method for constructing the petrophysical analysis quantitative plate based on the well log data comprises the following steps:
constructing a rock physical analysis quantity version based on a longitudinal wave velocity curve, a transverse wave velocity curve, a density curve and a logging oil-gas-water interpretation result in logging data, wherein the abscissa of the rock physical analysis quantity version is X, the ordinate is Y, and X is rho/vs，Y＝ρvpRho is density in g/cm3，vsIs the transverse wave velocity in Km/s, vpThe unit is the velocity of the longitudinal wave, Km/s.
5. The method of any one of claims 1 to 4, wherein in the step (3), the specific method for fitting the linear trend line of the aqueous sandstone based on the petrophysical analysis quantitative plate comprises the following steps:
in the rock physical analysis quantity version, selecting data points of a water-containing sandstone storage area of a target interval to perform linear fitting to obtain a water-containing sandstone linear trend line, wherein the trend line formula obtained by fitting is as follows:
Y＝kX+d
wherein X is rho/vs，Y＝ρvp；
k is the slope, d is the intercept, ρ is the density in g/cm3，vsIs the transverse wave velocity in Km/s, vpThe unit is the velocity of the longitudinal wave, Km/s.
6. The method of any one of claims 1 to 5, wherein in the step (4), the specific method for constructing the gas sandstone prediction factor based on the petrophysical analysis quantitative plate and the water-containing sandstone linear trend line comprises the following steps:
constructing a calculation formula of the gas-containing sandstone prediction factor based on the compressional wave velocity, the shear wave velocity, the density curve and the linear trend line of the water-containing sandstone obtained by fitting, wherein the calculation formula is as follows:
wherein F is a gas-containing sandstone prediction factor, k is a slope, rho is a density, and the unit is g/cm3，vsIs the transverse wave velocity in Km/s, vpThe unit is the velocity of the longitudinal wave, Km/s.
7. The method of any one of claims 1-6, wherein in step (6), the specific method for reservoir gas bearing detection based on the gas sand predictor seismic data volume comprises:
and predicting and researching the distribution rule of the gas reservoir in the work area based on the gas-containing sandstone prediction factor seismic data volume, wherein the distribution rule corresponds to the distribution range of the gas-containing sandstone reservoir when the gas-containing sandstone prediction factor is smaller than a gas-containing sandstone prediction threshold value, corresponds to the distribution range of mudstone when the gas-containing sandstone prediction factor is larger than the mudstone prediction threshold value, and corresponds to the distribution range of the water-containing sandstone reservoir when the gas-containing sandstone prediction factor is positioned between the gas-containing sandstone prediction factor and the mudstone.
8. A reservoir gas bearing detection apparatus for carrying out the reservoir gas bearing detection method according to any one of claims 1 to 7, the apparatus comprising:
the input unit is used for inputting logging data and seismic data;
the analysis unit is used for constructing a rock physical analysis quantity version based on the logging data;
the fitting unit is used for fitting the linear trend line of the water-containing sandstone based on the rock physical analysis quantity plate;
the device comprises a construction unit, a calculation unit and a calculation unit, wherein the construction unit is used for constructing a gas-containing sandstone prediction factor based on a rock physical analysis quantity version and a water-containing sandstone linear trend line;
the calculation unit is used for substituting the seismic data into a gas-containing sandstone prediction factor calculation formula to obtain a gas-containing sandstone prediction factor seismic data volume;
and the detection unit is used for detecting the gas content of the reservoir based on the gas-containing sandstone prediction factor seismic data body.
9. A reservoir gas bearing detection apparatus as defined in claim 8, wherein the input unit is configured to input well log data comprising: longitudinal wave velocity curve, transverse wave velocity curve, density curve and well logging oil gas water interpretation result.
10. A reservoir gas bearing detection apparatus as claimed in claim 8 or 9, wherein in the input unit, the seismic data comprises: and obtaining longitudinal wave velocity, transverse wave velocity and density seismic data volume by pre-stack seismic inversion.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011175054.3A CN112505761B (en) | 2020-10-28 | 2020-10-28 | Reservoir gas content detection method and device |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011175054.3A CN112505761B (en) | 2020-10-28 | 2020-10-28 | Reservoir gas content detection method and device |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112505761A true CN112505761A (en) | 2021-03-16 |
| CN112505761B CN112505761B (en) | 2022-10-04 |
Family
ID=74954392
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202011175054.3A Active CN112505761B (en) | 2020-10-28 | 2020-10-28 | Reservoir gas content detection method and device |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112505761B (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115822580A (en) * | 2022-12-15 | 2023-03-21 | 吉林大学 | Method for quantitatively predicting spatial distribution of deep volcanic gas layer |
Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2012006604A2 (en) * | 2010-07-08 | 2012-01-12 | Geco Technology B.V. | Methods and Devices for Transformation of Collected Data for Improved Visualization Capability |
| CN102854531A (en) * | 2012-09-11 | 2013-01-02 | 中国石油天然气股份有限公司 | Multi-scale rock physical plate method and device for reservoir hydrocarbon detection |
| CN103527184A (en) * | 2013-10-28 | 2014-01-22 | 北京大学 | Method and system for predicting dolomite reservoir |
| US20150362623A1 (en) * | 2014-06-12 | 2015-12-17 | Westerngeco, Llc | Joint inversion of attributes |
| WO2016161914A1 (en) * | 2015-04-07 | 2016-10-13 | 四川行之智汇知识产权运营有限公司 | Method for predicting reservoir lithogenous phase using geology and logging information |
| CN106597544A (en) * | 2016-11-25 | 2017-04-26 | 中国石油天然气股份有限公司 | Method and device for predicting brittleness of compact oil and gas reservoir |
| CN109115987A (en) * | 2018-07-20 | 2019-01-01 | 中国石油天然气股份有限公司 | Fluid factor evaluation method and device based on rock physical model |
| US20190179049A1 (en) * | 2017-12-07 | 2019-06-13 | Saudi Arabian Oil Company | Mapping chemostratigraphic signatures of a reservoir with rock physics and seismic inversion |
-
2020
- 2020-10-28 CN CN202011175054.3A patent/CN112505761B/en active Active
Patent Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2012006604A2 (en) * | 2010-07-08 | 2012-01-12 | Geco Technology B.V. | Methods and Devices for Transformation of Collected Data for Improved Visualization Capability |
| CN102854531A (en) * | 2012-09-11 | 2013-01-02 | 中国石油天然气股份有限公司 | Multi-scale rock physical plate method and device for reservoir hydrocarbon detection |
| CN103527184A (en) * | 2013-10-28 | 2014-01-22 | 北京大学 | Method and system for predicting dolomite reservoir |
| US20150362623A1 (en) * | 2014-06-12 | 2015-12-17 | Westerngeco, Llc | Joint inversion of attributes |
| WO2016161914A1 (en) * | 2015-04-07 | 2016-10-13 | 四川行之智汇知识产权运营有限公司 | Method for predicting reservoir lithogenous phase using geology and logging information |
| CN106597544A (en) * | 2016-11-25 | 2017-04-26 | 中国石油天然气股份有限公司 | Method and device for predicting brittleness of compact oil and gas reservoir |
| US20190179049A1 (en) * | 2017-12-07 | 2019-06-13 | Saudi Arabian Oil Company | Mapping chemostratigraphic signatures of a reservoir with rock physics and seismic inversion |
| CN109115987A (en) * | 2018-07-20 | 2019-01-01 | 中国石油天然气股份有限公司 | Fluid factor evaluation method and device based on rock physical model |
Non-Patent Citations (5)
| Title |
|---|
| LEI WANG ET AL.: "Reservoir fluid discrimination with an alternative fluid indicator based on rock physics templates", 《2016 SEG INTERNATIONAL EXPOSITION AND 86TH ANNUAL MEETING》 * |
| PER AVSETH ET AL.: "Seismic screening for hydrocarbon prospects using rock-physics attributes", 《THE LEADING EDGE》 * |
| 印兴耀等: "岩石物理驱动下地震流体识别研究", 《中国科学: 地球科学》 * |
| 王永刚等: "岩石物理模型分析在致密砂岩含气储层预测中的应用", 《2015 年物探技术研讨会》 * |
| 王磊等: "基于广义泊松阻尼因子的流体检测技术与应用", 《中国石油学会 2019 年物探技术研讨会》 * |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115822580A (en) * | 2022-12-15 | 2023-03-21 | 吉林大学 | Method for quantitatively predicting spatial distribution of deep volcanic gas layer |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112505761B (en) | 2022-10-04 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Kadkhodaie-Ilkhchi et al. | A review of reservoir rock typing methods in carbonate reservoirs: Relation between geological, seismic, and reservoir rock types | |
| US4646240A (en) | Method and apparatus for determining geological facies | |
| Li et al. | Prediction of fracture density using genetic algorithm support vector machine based on acoustic logging data | |
| CN111025384B (en) | Reservoir stratum prediction method and device based on waveform classification intersection fusion | |
| CN104834006B (en) | The treating method and apparatus of structure origin cause of formation deposition parameter figure | |
| US6718265B2 (en) | Petrophysical property estimation using an acoustic calibration relationship | |
| CN110632652A (en) | Method and device for determining oil-gas migration mode of fault block oil-gas reservoir | |
| CN112684503B (en) | Reservoir fluid identification method and device based on post-stack seismic attributes | |
| Bourges et al. | A practical review of geostatistical processing applied to geophysical data: methods and applications | |
| CN105240006A (en) | Oil and water layer recognition method suitable for volcanic reservoir | |
| CN116168224A (en) | Machine learning lithology automatic identification method based on imaging gravel content | |
| Oyeyemi et al. | Hydrocarbon resource evaluation using combined petrophysical analysis and seismically derived reservoir characterization, offshore Niger Delta | |
| Bagheri et al. | Fracture permeability estimation utilizing conventional well logs and flow zone indicator | |
| CN112505761B (en) | Reservoir gas content detection method and device | |
| RU2619803C2 (en) | Determining fluid motion path | |
| CN108073744B (en) | Turbidimetric fan propulsion distance prediction method based on geology statistics | |
| Zhang et al. | Data-processing and recognition of seepage and microseepage anomalies of acid-extractable hydrocarbons in the south slope of the Dongying depression, eastern China | |
| George et al. | Estimation of aquifer hydraulic parameters via complementing surfacial geophysical measurement by laboratory measurements on the aquifer core samples | |
| CN114509812B (en) | Clastic rock reservoir dessert area prediction method, device, medium and electronic equipment | |
| Ge et al. | Determination of total organic carbon (TOC) in tight reservoir using empirical mode decomposition-support vector regression (EMD-SVR): a case study from XX-1 Basin, Western China | |
| Jia et al. | Seismic response analysis and distribution prediction of source rocks in a survey of the South China Sea | |
| CN113806674B (en) | Quantification method and device for longitudinal scale of ancient river channel, electronic equipment and storage medium | |
| Méndez et al. | Predicting and 3D modeling of karst zones using seismic facies analysis in Ordovician carbonates of the Tahe oilfield, China | |
| Zhang et al. | A method for eliminating caprock thickness influence on anomaly intensities in geochemical surface survey for hydrocarbons | |
| CN112505764A (en) | High-porosity hydrocarbon-containing sandstone reservoir prediction method and device |
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 |