CN114113537B - Method and device for quantitatively identifying diagenetic effect of compact sandstone - Google Patents
Method and device for quantitatively identifying diagenetic effect of compact sandstone Download PDFInfo
- Publication number
- CN114113537B CN114113537B CN202010902086.2A CN202010902086A CN114113537B CN 114113537 B CN114113537 B CN 114113537B CN 202010902086 A CN202010902086 A CN 202010902086A CN 114113537 B CN114113537 B CN 114113537B
- Authority
- CN
- China
- Prior art keywords
- rock
- diagenetic
- parameters
- elastic
- test sample
- 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
- 230000000694 effects Effects 0.000 title claims abstract description 97
- 238000000034 method Methods 0.000 title claims abstract description 37
- 239000011435 rock Substances 0.000 claims abstract description 190
- 238000012360 testing method Methods 0.000 claims abstract description 92
- 239000000203 mixture Substances 0.000 claims description 43
- 239000004927 clay Substances 0.000 claims description 34
- 230000000704 physical effect Effects 0.000 claims description 30
- 230000035699 permeability Effects 0.000 claims description 25
- 238000004590 computer program Methods 0.000 claims description 16
- 238000004458 analytical method Methods 0.000 claims description 15
- 238000003860 storage Methods 0.000 claims description 12
- 238000001514 detection method Methods 0.000 claims description 7
- 238000004364 calculation method Methods 0.000 claims description 6
- 238000012163 sequencing technique Methods 0.000 claims description 6
- 230000009471 action Effects 0.000 claims description 5
- 238000005056 compaction Methods 0.000 claims description 5
- 238000000691 measurement method Methods 0.000 claims description 4
- 230000003628 erosive effect Effects 0.000 claims description 3
- 238000010586 diagram Methods 0.000 description 13
- 239000011148 porous material Substances 0.000 description 13
- OYPRJOBELJOOCE-UHFFFAOYSA-N Calcium Chemical compound [Ca] OYPRJOBELJOOCE-UHFFFAOYSA-N 0.000 description 11
- 239000011575 calcium Substances 0.000 description 11
- 229910052791 calcium Inorganic materials 0.000 description 11
- 230000006870 function Effects 0.000 description 10
- 238000011160 research Methods 0.000 description 10
- 239000004568 cement Substances 0.000 description 8
- 238000004891 communication Methods 0.000 description 7
- 238000005516 engineering process Methods 0.000 description 7
- 238000011161 development Methods 0.000 description 6
- 239000012530 fluid Substances 0.000 description 6
- 239000002245 particle Substances 0.000 description 6
- 238000012545 processing Methods 0.000 description 6
- 239000000872 buffer Substances 0.000 description 4
- 238000005553 drilling Methods 0.000 description 4
- 239000004576 sand Substances 0.000 description 4
- 238000005260 corrosion Methods 0.000 description 3
- 230000007797 corrosion Effects 0.000 description 3
- 230000003247 decreasing effect Effects 0.000 description 3
- 238000004445 quantitative analysis Methods 0.000 description 3
- 239000004215 Carbon black (E152) Substances 0.000 description 2
- 206010013883 Dwarfism Diseases 0.000 description 2
- 230000015572 biosynthetic process Effects 0.000 description 2
- 238000005136 cathodoluminescence Methods 0.000 description 2
- 230000006835 compression Effects 0.000 description 2
- 238000007906 compression Methods 0.000 description 2
- 238000004090 dissolution Methods 0.000 description 2
- 238000009826 distribution Methods 0.000 description 2
- 230000002349 favourable effect Effects 0.000 description 2
- 229930195733 hydrocarbon Natural products 0.000 description 2
- 150000002430 hydrocarbons Chemical class 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 239000010453 quartz Substances 0.000 description 2
- VYPSYNLAJGMNEJ-UHFFFAOYSA-N silicon dioxide Inorganic materials O=[Si]=O VYPSYNLAJGMNEJ-UHFFFAOYSA-N 0.000 description 2
- 229910021532 Calcite Inorganic materials 0.000 description 1
- 238000002441 X-ray diffraction Methods 0.000 description 1
- 230000001413 cellular effect Effects 0.000 description 1
- 239000003795 chemical substances by application Substances 0.000 description 1
- 239000002734 clay mineral Substances 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000008021 deposition Effects 0.000 description 1
- 230000004069 differentiation Effects 0.000 description 1
- 238000002637 fluid replacement therapy Methods 0.000 description 1
- 238000002347 injection Methods 0.000 description 1
- 239000007924 injection Substances 0.000 description 1
- 229910052500 inorganic mineral Inorganic materials 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 239000011707 mineral Substances 0.000 description 1
- 238000010295 mobile communication Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 239000000243 solution Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N33/00—Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
- G01N33/24—Earth materials
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume, or surface-area of porous materials
- G01N15/08—Investigating permeability, pore-volume, or surface area of porous materials
- G01N15/082—Investigating permeability by forcing a fluid through a sample
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume, or surface-area of porous materials
- G01N15/08—Investigating permeability, pore-volume, or surface area of porous materials
- G01N15/088—Investigating volume, surface area, size or distribution of pores; Porosimetry
-
- 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. analysis, for interpretation, for correction
- 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
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/50—Corrections or adjustments related to wave propagation
- G01V2210/51—Migration
- G01V2210/512—Pre-stack
-
- 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
- G01V2210/6242—Elastic parameters, e.g. Young, Lamé or Poisson
-
- 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
- G01V2210/6244—Porosity
-
- 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
- G01V2210/6246—Permeability
Abstract
The invention provides a method and a device for quantitatively identifying tight sandstone diagenetic effects, wherein the method comprises the following steps: acquiring a plurality of rock test samples under different diagenetic effects, and respectively detecting physical parameters and rock component factors of the rock test samples; calculating elastic parameters of the rock test samples, and establishing fitting relation equations among the rock component factors, the physical parameters and the elastic parameters under different diagenetic effects through the rock component factors, the physical parameters and the elastic parameters of each rock test sample; obtaining elastic parameters of a region to be tested through pre-stack seismic inversion, and obtaining physical parameters and rock component factors of the region to be tested through inversion of the elastic parameters before stack; and obtaining the diagenetic effect of the region to be tested through the fitting relation equation according to the rock component factors, the elastic parameters and the physical parameters of the region to be tested.
Description
Technical Field
The invention relates to a high-quality reservoir prediction technology of a tight sandstone reservoir in the petroleum industry, in particular to a technology for identifying the differential diagenetic effect of the tight sandstone reservoir by combining geological and earthquake, and particularly relates to a method and a device for quantitatively identifying the diagenetic effect of the tight sandstone.
Background
Tight sandstone generally refers to sandstone with a porosity of <10%, and a permeability of <1.0 mD. At present, dense sandstone gas is found in more than 70 basins worldwide, and is mainly distributed in North America, europe and asia-Tai areas, and dense sandstone gas reservoirs are developed in areas of Sichuan, erdos, songyang and the like in China. Because of the factors of long ground history time, multiple construction period, complex diagenetic effect and the like, the research on the mechanism of the formation of the tight sandstone hydrocarbon reservoir is difficult, and the finding of the distribution of high-quality reservoirs in the tight sandstone is an urgent hope for solving the problems of hydrocarbon exploration and development.
In seismic exploration, seismic waves bring about information about subsurface rock and fluids in the form of travel time, reflected amplitude, and phase changes. However, this information is affected by a combination of factors such as the degree of cementing of the rock particles, the contact conditions of the particles, pressure, temperature, porosity, pore type, fluid type, saturation, seismic frequency, etc. Therefore, the earthquake rock physics builds a bridge with the connection of the earthquake information and the rock most basic parameters, and is a physical foundation for pre-stack reservoir prediction and a tie for connecting the earthquake and the oil reservoir engineering. Seismic petrophysical research can be divided into two major parts, theoretical research and application. Theoretical studies began in the 50 s of the 20 th century, principally mathematically describing the properties of fluid-containing rocks. The application research mainly comprises the following steps: 1) Petrophysical parameter testing and in-reservoir fluid property testing under different fluid saturation conditions in a laboratory simulated field reservoir conditions; 2) Analyzing sample test data, and establishing an empirical relationship model between seismic rock physical parameters and seismic attribute parameters; 3) Fluid replacement and transverse wave data are obtained during forward seismic modeling before the stack, and physical characteristics of reservoir rock are represented and the reservoir oil-gas-containing property is predicted; 4) Development and research of petrophysical test systems, and the like. Since the 90 s, many theories and experimental results of earthquake rock physics have greatly promoted the great development of bright spot (dark spot) technology, AVO technology, DHI technology, 4D earthquake and injection and production monitoring technology. Seismic petrophysical has been widely studied in terms of petrophysical models, petrophysical determination, reservoir parameter analysis, anisotropy, attenuation, etc. 5. Through retrieval, no method for quantitatively predicting diagenetic effect in a tight sandstone reservoir by utilizing seismic data exists at present.
Disclosure of Invention
The invention aims to provide a method and a device for quantitatively identifying the diagenetic effect of compact sandstone, finely describing the influence of the diagenetic effect on the reservoir property of a compact sandstone reservoir, improving the success rate of drilling and saving the exploration investment.
In order to achieve the above purpose, the method for quantitatively identifying the diagenetic effect of the tight sandstone provided by the invention specifically comprises the following steps: acquiring a plurality of rock test samples under different diagenetic effects, and respectively detecting physical parameters and rock component factors of the rock test samples; calculating elastic parameters of the rock test samples, and establishing fitting relation equations among the rock component factors, the physical parameters and the elastic parameters under different diagenetic effects through the rock component factors, the physical parameters and the elastic parameters of each rock test sample; obtaining elastic parameters of a region to be tested through pre-stack seismic inversion, and obtaining physical parameters and rock component factors of the region to be tested through inversion of the elastic parameters before stack; and obtaining the diagenetic effect of the region to be tested through the fitting relation equation according to the rock component factors, the elastic parameters and the physical parameters of the region to be tested.
In the above method for quantitatively identifying tight sandstone diagenesis, it is preferable that the detecting physical property parameters and rock composition factors of the rock test sample, respectively, comprises: retrieving the permeability and the porosity of the rock test sample by a standard gas measurement method, and obtaining physical parameters of the rock test sample according to the permeability and the porosity; and obtaining rock component factors of the rock test sample according to the influence degree of each component in the rock test sample on the permeability and the porosity of the rock test sample under different diagenetic effects.
In the above method for quantitatively identifying tight sandstone diagenetic effects, preferably, obtaining rock composition factors of the rock test sample based on the extent to which each component in the rock test sample affects the permeability and porosity of the rock test sample under different diagenetic effects comprises: detecting the influence degree of each component in the rock test sample on the permeability and the porosity of the rock test sample under the action of different diagenetic rocks through a rock physical test, and obtaining the influence value corresponding to each component; and sequencing the influence values, and obtaining rock component factors corresponding to the influence values according to sequencing results.
In the above method for quantitatively identifying tight sandstone diagenetic events, it is preferable that establishing a fit relationship equation between the rock composition factor, the physical property parameter, and the elastic parameter under different diagenetic events by the rock composition factor, the physical property parameter, and the elastic parameter of each rock test sample comprises: constructing an association equation of the physical property parameter and the elastic parameter of the rock test sample under different diagenetic effects according to the physical property parameter and the elastic parameter; and establishing a fitting relation equation among the rock composition factors, the physical parameters and the elastic parameters under different diagenetic effects according to the rock composition factors and the association equation corresponding to the elastic parameters.
In the above method for quantitatively identifying tight sandstone diagenetic events, preferably, constructing correlation equations for the physical property parameters and the elastic parameters of the rock test sample under different diagenetic events from the physical property parameters and the elastic parameters comprises: when diagenetic is cemented and the rock composition factors are mud siliceous and calcareous, the correlation equation includes: when the rock component factors are mud siliceous:
when the rock composition factor is calcareous:
in the above formula, V p For longitudinal wave velocity, V s In order to be a transverse wave velocity,is porosity.
In the above method for quantitatively identifying tight sandstone diagenetic effects, preferably, establishing a fitting relation equation among the rock composition factor, the physical property parameter and the elastic parameter under different diagenetic effects according to the rock composition factor and the correlation equation corresponding to the elastic parameter includes: when diagenetic is cemented and the rock composition factor is mud siliceous, the fitted relation equation includes:
V p ＝5.99-0.039·V clay -0.064·φ；
V s ＝3.71-0.037·V clay -0.037·φ；
in the above formula, V p For longitudinal wave velocity, V s For transverse wave velocity, V clay For the content of mud-siliceous in the rock test sample,is porosity.
In the above method of quantitatively identifying tight sandstone diagenetic, preferably, the diagenetic comprises any of compaction, cementing, cross-talk, erosion.
The invention also provides a device for quantitatively identifying the diagenetic effect of the compact sandstone, which comprises a detection module, a calculation module, an analysis module and an identification module; the detection module is used for acquiring a plurality of rock test samples under different diagenetic effects and respectively detecting physical parameters and rock component factors of the rock test samples; the calculation module is used for calculating the elasticity parameters of the rock test samples, and establishing fitting relation equations among the rock component factors, the physical parameters and the elasticity parameters under different diagenetic effects through the rock component factors, the physical parameters and the elasticity parameters of each rock test sample; the analysis module is used for obtaining elastic parameters of the region to be tested through pre-stack seismic inversion, and obtaining physical parameters and rock component factors of the region to be tested through pre-stack elastic parameter inversion; the identification module is used for obtaining the diagenetic effect of the region to be tested through the fitting relation equation according to the rock composition factors, the elastic parameters and the physical parameters of the region to be tested.
The invention also provides an electronic device comprising a memory, a processor and a computer program stored on the memory and executable on the processor, the processor implementing the above method when executing the computer program.
The present invention also provides a computer readable storage medium storing a computer program for executing the above method.
By the method and the device for quantitatively identifying the diagenetic effect of the compact sandstone, provided by the invention, the influence of the diagenetic effect on the reservoir property of the compact sandstone reservoir is accurately described, the success rate of drilling is effectively improved, the exploration investment is saved, and the method and the device have good application effects.
Drawings
The accompanying drawings, which are included to provide a further understanding of the invention and are incorporated in and constitute a part of this application, illustrate and together with the description serve to explain the invention. In the drawings:
FIG. 1 is a flow chart of a method for quantitatively identifying tight sandstone diagenetic events according to an embodiment of the present invention;
FIG. 2 is a schematic diagram of a rock test sample detection flow provided by an embodiment of the present invention;
FIG. 3 is a schematic illustration of a representative sheet with varying clay content samples affecting porosity in accordance with one embodiment of the present invention;
FIG. 4 is a schematic diagram of the relationship between longitudinal wave velocity and porosity of a tight sandstone reservoir according to an embodiment of the present invention;
FIG. 5 is a schematic diagram of shear wave velocity versus porosity for different diagenetic effects of a tight sandstone reservoir according to an embodiment of the present invention;
FIG. 6 is a schematic diagram of the relationship between the reservoir sandstone speed and the porosity and clay content according to an embodiment of the present invention;
FIG. 7 is a schematic structural diagram of an apparatus for quantitatively identifying tight sandstone diagenetic effects according to an embodiment of the present invention;
FIG. 8 is a flow chart of a differential diagenetic quantitative prediction according to an embodiment of the present invention;
fig. 9 is a schematic structural diagram of an electronic device according to an embodiment of the invention.
Detailed Description
The following will describe embodiments of the present invention in detail with reference to the drawings and examples, thereby solving the technical problems by applying technical means to the present invention, and realizing the technical effects can be fully understood and implemented accordingly. It should be noted that, as long as no conflict is formed, each embodiment of the present invention and each feature of each embodiment may be combined with each other, and the formed technical solutions are all within the protection scope of the present invention.
Additionally, the steps illustrated in the flowcharts of the figures may be performed in a computer system such as a set of computer executable instructions, and although a logical order is illustrated in the flowcharts, in some cases the steps illustrated or described may be performed in an order other than that herein.
Referring to fig. 1, the method for quantitatively identifying the diagenetic effect of tight sandstone provided by the present invention specifically includes:
s101, acquiring a plurality of rock test samples under different diagenetic effects, and respectively detecting physical parameters and rock component factors of the rock test samples;
s102, calculating elastic parameters of the rock test samples, and establishing fitting relation equations among the rock component factors, the physical parameters and the elastic parameters under different diagenetic effects through the rock component factors, the physical parameters and the elastic parameters of each rock test sample;
s103, obtaining elastic parameters of a region to be tested through pre-stack seismic inversion, and obtaining physical parameters and rock component factors of the region to be tested through inversion of the elastic parameters before stack;
s104, obtaining the diagenetic effect of the region to be tested through the fitting relation equation according to the rock composition factors, the elastic parameters and the physical parameters of the region to be tested.
In the above embodiments, the diagenetic effect may comprise any of compaction, cementing, cross-over, erosion; therefore, the influence of the diagenetic effect on the reservoir performance of the tight sandstone reservoir is accurately described by utilizing the association relation of the rock composition factors, the physical parameters and the elastic parameters under different diagenetic effects, the success rate of drilling is effectively improved, the exploration investment is saved, and the application effect is good.
Referring to fig. 2, in an embodiment of the present invention, detecting physical parameters and rock composition factors of the rock test sample respectively includes:
s201, retrieving the permeability and the porosity of the rock test sample by a standard gas measurement method, and obtaining physical parameters of the rock test sample according to the permeability and the porosity;
s202, obtaining rock component factors of the rock test sample according to the influence degree of each component in the rock test sample on the permeability and the porosity of the rock test sample under the action of different diagenesis.
Wherein, according to the influence degree of each component in the rock test sample on the permeability and the porosity of the rock test sample under different diagenetic effects, the rock component factors of the rock test sample are obtained and comprise: detecting the influence degree of each component in the rock test sample on the permeability and the porosity of the rock test sample under the action of different diagenetic rocks through a rock physical test, and obtaining the influence value corresponding to each component; and sequencing the influence values, and obtaining rock component factors corresponding to the influence values according to sequencing results.
In one embodiment of the present invention, establishing a fit relation equation between the rock composition factor, the physical property parameter, and the elastic parameter under different diagenetic effects by the rock composition factor, the physical property parameter, and the elastic parameter of each rock test sample comprises: constructing an association equation of the physical property parameter and the elastic parameter of the rock test sample under different diagenetic effects according to the physical property parameter and the elastic parameter; and establishing a fitting relation equation among the rock composition factors, the physical parameters and the elastic parameters under different diagenetic effects according to the rock composition factors and the association equation corresponding to the elastic parameters.
In the above embodiment, constructing the correlation equation of the physical property parameter and the elastic parameter of the rock test sample under different diagenetic effects according to the physical property parameter and the elastic parameter includes: when diagenetic is cemented and the rock composition factors are mud siliceous and calcareous, the correlation equation includes: when the rock component factors are mud siliceous:
when the rock composition factor is calcareous:
in the above formula, V p For longitudinal wave velocity, V s In order to be a transverse wave velocity,is porosity.
In the above embodiment, establishing the fitting relation equation among the rock composition factor, the physical property parameter and the elastic parameter under different diagenetic effects according to the rock composition factor and the correlation equation corresponding to the elastic parameter includes: when diagenetic is cemented and the rock composition factor is mud siliceous, the fitted relation equation includes:
V p ＝5.99-0.039·V clay -0.064·φ；
V s ＝3.71-0.037·V clay -0.037·φ；
in the above formula, V p For longitudinal wave velocity, V s For transverse wave velocity, V clay For the content of mud-siliceous in the rock test sample,is porosity.
In order to more clearly understand the method for quantitatively recognizing the diagenetic effect of the tight sandstone provided by the present invention, the following description will be given in detail with reference to practical working examples, which are only application modes for facilitating understanding the method for quantitatively recognizing the diagenetic effect of the tight sandstone provided by the present invention, and are not limited in any way.
Overall, the method for quantitatively identifying tight sandstone diagenetic effects provided by the invention mainly comprises two parts: quantitative analysis of reservoir properties by differential diagenetic effects and establishment of quantitative relationships between differential diagenetic effects and elastic parameters based on petrophysical testing and analysis, specifically:
1. quantitative analysis of reservoir properties by differential diagenetic effects:
the compaction (compression dissolution), cementing, cross-generation, corrosion and other diagenetic effects have important influence on the dwarf-system cold mountain group and the sand-section sandstone reservoir of a certain basin through microscopic thin plates, cathodoluminescence, scanning electron microscope and other analysis.
The sample permeability and porosity are all obtained by adopting a standard gas measurement method. Statistics of the research results of the sampled rock core shows that the porosity of the sandstone sample is obviously controlled by the clay content, and when the clay content is less than 8%, the porosity and the permeability are increased along with the increase of the clay content. The main reason for this is that when the sample contains no or little clay, the initial permeability of the rock is high, so that the fluid containing certain minerals tends to flow in the pores, so that early calcium cementation and later siliceous cementation tend to occur, and pure sandstone tends to occur more easily, so that the effective porosity of the rock is rapidly reduced and densified, and the permeability of the rock is reduced in the course of the reduction of the porosity. As the clay content increases, clay plugs the pore throats causing a gradual decrease in the native permeability of the rock, early cementing is less likely to occur while retaining higher porosity and permeability. When the clay content is more than 8%, the porosity is decreased instead with the increase of the clay content, mainly because the increase of the clay content causes more clay to be filled in the primary pores, so that the porosity is decreased, and at the same time the permeability of the rock sample is gradually decreased, as shown in fig. 3. The sample permeability is obviously affected by calcium cementation, and the porosity and the permeability are reduced along with the increase of the calcium content.
At the same porosity, the speed difference between the calcareous cementing sample and the clay-siliceous cementing sample reaches more than 1000 m/s. As can be seen from analysis of the corresponding cast sheet of the sample, the low velocity mud, siliceous cement sample is primarily clay cement, with clay particles filled between rigid quartz particles, with typical argillaceous contact cement and pore cement, as shown in fig. 4; for calcareous cement samples, very low porosity sample quartz particles are suspended in the calcareous cement, typically basal cementation, as shown in fig. 5. Thus, cement differentiation is a significant cause of differences in rock velocities in the two regions. The pore structures of the two samples are obviously different, the calcium cement sample has strong diagenetic effect, the pores are mainly corrosion pores, and the rigidity is high due to the large aspect ratio of the pores; mud and siliceous diagenetic effects are weak, pores are mainly inter-particle pores, and the aspect ratio of the pores is small so that the rigidity is small. The relation between the elasticity parameters and the physical parameters of the sandstone sample under different cementing conditions can be established by combining the petrophysical test results:
longitudinal wave velocity-porosity relationship of mud and siliceous cementitious sandstone samples:
transverse wave velocity-porosity relationship of mud and siliceous cementitious sandstone samples:
longitudinal wave velocity-porosity relationship of calcareous cemented sandstone sample:
the transverse wave velocity-porosity relationship of the calcareous cemented sandstone sample:
wherein V is p Representing longitudinal wave velocity, V s Representing the velocity of the transverse wave,representing sample porosity.
2. Based on petrophysical testing and analysis, establishing a quantitative relationship between differential diagenetic and elastic parameters:
establishing a quantitative relation formula for researching the calcium content and clay content of a sample to the longitudinal wave speed and transverse wave speed ratio by using a petrophysical test; there is a trend of increasing longitudinal and transverse wave velocities with increasing calcium content at the same porosity, and rock samples exhibit the highest velocity if they are cemented by early calcite. Also the sample longitudinal and transverse wave velocities decrease with increasing clay content. In addition, from the experimental results, the effect of the cementing agent on the speed is not single, and part of samples have lower speed at the conditions of higher calcium content, lower porosity and clay content, and the phenomenon reflects the effect of the pore structure on the rock earthquake elasticity. The longitudinal wave impedance and the transverse wave impedance and the velocity ratio all show increasing trend along with the increase of the calcium content, and the longitudinal wave impedance and the transverse wave impedance gradually decrease along with the increase of the clay content and the velocity ratio shows increasing trend, so that the rock sample shows the characteristic similar to V shape in the intersection diagram of the impedance ratio and the velocity ratio. The purer sandstone shows the lowest longitudinal wave speed ratio and the lowest transverse wave speed ratio, and the speed ratio is increased along with the increase of the clay content and the calcium content, so that the speed ratio has a certain polynomials as a main index for dividing mud and sandstone, and the lithology needs to be comprehensively judged by combining the impedance. The effect of the clay content on the speed ratio of the high-grinding sand temple group sample is more obvious, and the speed ratio is increased from 1.5 to about 1.8 when the clay content is small. The speed ratio is increased from 1.5 to about 1.7 by affecting the speed with more calcium content. Considering the overall effect of porosity and argillaceous content of the reservoir samples of the dwarfism temple group in the high-abrasive region on rock velocity, the velocity was fitted to the relationship between porosity and argillaceous content as shown in fig. 6:
V p ＝5.99-0.039·V clay -0.064·φ (5)
V s ＝3.71-0.037·V clay -0.037·φ (6)
referring to fig. 7, the invention further provides a device for quantitatively identifying the diagenetic effect of the tight sandstone, which comprises a detection module, a calculation module, an analysis module and an identification module; the detection module is used for acquiring a plurality of rock test samples under different diagenetic effects and respectively detecting physical parameters and rock component factors of the rock test samples; the calculation module is used for calculating the elasticity parameters of the rock test samples, and establishing fitting relation equations among the rock component factors, the physical parameters and the elasticity parameters under different diagenetic effects through the rock component factors, the physical parameters and the elasticity parameters of each rock test sample; the analysis module is used for obtaining elastic parameters of the region to be tested through pre-stack seismic inversion, and obtaining physical parameters and rock component factors of the region to be tested through pre-stack elastic parameter inversion; the identification module is used for obtaining the diagenetic effect of the region to be tested through the fitting relation equation according to the rock composition factors, the elastic parameters and the physical parameters of the region to be tested.
In order to better understand the application flow of the device for quantitatively identifying the diagenetic effect of the tight sandstone provided by the present invention, please refer to fig. 8, the following detailed description is given for the flow:
1. selecting petrophysical test samples for different sedimentary diagenetic effects in a study area;
the difference of the formation action of the dwarfism compact sandstone reservoir in a certain basin causes great difference of reservoir storage performance. Therefore, in order to predict the distribution of high-quality reservoirs in a research area, the difference of the diagenetic effects of the reservoirs in the research area needs to be predicted by utilizing the seismic data, so the invention adopts the thought of seismic rock physics to establish the relationship between the differentiated diagenetic effects and the elastic parameters. The above objective is accomplished by first selecting petrophysical test samples for petrophysical testing and analysis for different diagenetic effects in a research area. The embodiment researches that the selected compact sandstone whole core sample is taken from compact sandstone of dwarf system sand temple group in a basin high-grinding area, 42 sandstone samples are taken from 2 sections of the target reservoir sand temple, and 1 section of the sand temple is considered.
2. Analysis of petrophysical characteristics and diagenetic effects of rock samples based on petrophysical test results, X-ray diffraction clay minerals of the samples and whole rock analysis results
Factors influencing reservoir development include the deposition environment, diagenetic effects, and buried history, among others, where diagenetic effects are one of the key factors influencing reservoir physical properties. The compaction (compression dissolution), cementing, cross-generation, corrosion and other diagenetic effects have important influence on the dwarf-system cold mountain group and the sand-section sandstone reservoir of a certain basin through microscopic thin plates, cathodoluminescence, scanning electron microscope and other analysis.
3. Quantitative analysis of reservoir properties by differential diagenetic effects has been described in detail in the foregoing examples and will not be described in detail herein.
4. The fitting relationship between the differential diagenetic effect and the elastic parameter is established and described in detail in the foregoing embodiments, and will not be described in detail herein.
5. Optimized sensitive parameter development pre-stack parameter inversion
And calculating the longitudinal wave speed and the transverse wave speed and the density by utilizing the pre-stack earthquake simultaneous inversion result, and carrying out the prediction of the clay content and the porosity by utilizing the pre-stack inversion result to obtain the clay content and the porosity.
6. Quantitative prediction of differential diagenetic effects
The diagenetic effect of the tight sandstone reservoir can be quantitatively reflected by using the argillaceous content and porosity combined fitting relation equation obtained by inversion of the prestack elastic parameters.
7. Development of favorable zone prediction using differential diagenetic prediction results
The difference of physical properties of the reservoir can be quantitatively established for the mud content and the calcium content caused by different diagenetic effects, and then the favorable region is preferred to guide exploration.
By the method and the device for quantitatively identifying the diagenetic effect of the compact sandstone, provided by the invention, the influence of the diagenetic effect on the reservoir property of the compact sandstone reservoir is accurately described, the success rate of drilling is effectively improved, the exploration investment is saved, and the method and the device have good application effects.
The invention also provides an electronic device comprising a memory, a processor and a computer program stored on the memory and executable on the processor, the processor implementing the above method when executing the computer program.
The present invention also provides a computer readable storage medium storing a computer program for executing the above method.
As shown in fig. 9, the electronic device 600 may further include: a communication module 110, an input unit 120, an audio processing unit 130, a display 160, a power supply 170. It is noted that the electronic device 600 need not include all of the components shown in fig. 9; in addition, the electronic device 600 may further include components not shown in fig. 9, to which reference is made to the related art.
As shown in fig. 9, the central processor 100, sometimes also referred to as a controller or operational control, may include a microprocessor or other processor device and/or logic device, which central processor 100 receives inputs and controls the operation of the various components of the electronic device 600.
The memory 140 may be, for example, one or more of a buffer, a flash memory, a hard drive, a removable media, a volatile memory, a non-volatile memory, or other suitable device. The information about failure may be stored, and a program for executing the information may be stored. And the central processor 100 can execute the program stored in the memory 140 to realize information storage or processing, etc.
The input unit 120 provides an input to the central processor 100. The input unit 120 is, for example, a key or a touch input device. The power supply 170 is used to provide power to the electronic device 600. The display 160 is used for displaying display objects such as images and characters. The display may be, for example, but not limited to, an LCD display.
The memory 140 may be a solid state memory such as Read Only Memory (ROM), random Access Memory (RAM), SIM card, or the like. But also a memory which holds information even when powered down, can be selectively erased and provided with further data, an example of which is sometimes referred to as EPROM or the like. Memory 140 may also be some other type of device. Memory 140 includes a buffer memory 141 (sometimes referred to as a buffer). The memory 140 may include an application/function storage 142, the application/function storage 142 for storing application programs and function programs or a flow for executing operations of the electronic device 600 by the central processor 100.
The memory 140 may also include a data store 143, the data store 143 for storing data, such as contacts, digital data, pictures, sounds, and/or any other data used by the electronic device. The driver storage 144 of the memory 140 may include various drivers of the electronic device for communication functions and/or for performing other functions of the electronic device (e.g., messaging applications, address book applications, etc.).
The communication module 110 is a transmitter/receiver 110 that transmits and receives signals via an antenna 111. A communication module (transmitter/receiver) 110 is coupled to the central processor 100 to provide an input signal and receive an output signal, which may be the same as in the case of a conventional mobile communication terminal.
Based on different communication technologies, a plurality of communication modules 110, such as a cellular network module, a bluetooth module, and/or a wireless local area network module, etc., may be provided in the same electronic device. The communication module (transmitter/receiver) 110 is also coupled to a speaker 131 and a microphone 132 via an audio processor 130 to provide audio output via the speaker 131 and to receive audio input from the microphone 132 to implement usual telecommunication functions. The audio processor 130 may include any suitable buffers, decoders, amplifiers and so forth. In addition, the audio processor 130 is also coupled to the central processor 100 so that sound can be recorded locally through the microphone 132 and so that sound stored locally can be played through the speaker 131.
It will be appreciated by those skilled in the art that embodiments of the present invention may be provided as a method, system, or computer program product. Accordingly, the present invention may take the form of an entirely hardware embodiment, an entirely software embodiment or an embodiment combining software and hardware aspects. Furthermore, the present invention may take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, and the like) having computer-usable program code embodied therein.
The present invention is described with reference to flowchart illustrations and/or block diagrams of methods, apparatus (systems) and computer program products according to embodiments of the invention. It will be understood that each flow and/or block of the flowchart illustrations and/or block diagrams, and combinations of flows and/or blocks in the flowchart illustrations and/or block diagrams, can be implemented by computer program instructions. These computer program instructions may be provided to a processor of a general purpose computer, special purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be stored in a computer-readable memory that can direct a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer-readable memory produce an article of manufacture including instruction means which implement the function specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
The foregoing description of the embodiments has been provided for the purpose of illustrating the general principles of the invention, and is not meant to limit the scope of the invention, but to limit the invention to the particular embodiments, and any modifications, equivalents, improvements, etc. that fall within the spirit and principles of the invention are intended to be included within the scope of the invention.
Claims (8)
1. A method for quantitatively identifying tight sandstone diagenesis, comprising:
acquiring a plurality of rock test samples under different diagenetic effects, and respectively detecting physical parameters and rock component factors of the rock test samples;
calculating elastic parameters of the rock test samples, and establishing fitting relation equations among the rock component factors, the physical parameters and the elastic parameters under different diagenetic effects through the rock component factors, the physical parameters and the elastic parameters of each rock test sample;
obtaining elastic parameters of a region to be tested through pre-stack seismic inversion, and obtaining physical parameters and rock component factors of the region to be tested through inversion of the elastic parameters before stack;
obtaining the diagenetic effect of the region to be tested through the fitting relation equation according to the rock component factors, the elastic parameters and the physical parameters of the region to be tested;
establishing a fitting relation equation among the rock composition factor, the physical property parameter and the elastic parameter under different diagenetic effects through the rock composition factor, the physical property parameter and the elastic parameter of each rock test sample comprises: constructing an association equation of the physical property parameter and the elastic parameter of the rock test sample under different diagenetic effects according to the physical property parameter and the elastic parameter; establishing a fitting relation equation among the rock composition factors, the physical parameters and the elastic parameters under different diagenetic effects according to the rock composition factors and the association equation corresponding to the elastic parameters;
constructing an association equation of the physical property parameter and the elastic parameter of the rock test sample under different diagenetic effects according to the physical property parameter and the elastic parameter comprises: when diagenetic is cemented and the rock composition factors are mud siliceous and calcareous, the correlation equation includes:
when the rock component factors are mud siliceous:
when the rock composition factor is calcareous:
in the above formula, V p For longitudinal wave velocity, V s In order to be a transverse wave velocity,is porosity.
2. The method of quantitatively identifying tight sandstone diagenesis of claim 1, wherein separately detecting physical parameters and rock composition factors of said rock test sample comprises:
retrieving the permeability and the porosity of the rock test sample by a standard gas measurement method, and obtaining physical parameters of the rock test sample according to the permeability and the porosity;
and obtaining rock component factors of the rock test sample according to the influence degree of each component in the rock test sample on the permeability and the porosity of the rock test sample under different diagenetic effects.
3. The method of quantitatively identifying tight sandstone diagenetic according to claim 1, wherein obtaining rock composition factors for the rock test sample based on the extent to which each component in the rock test sample affects the permeability and porosity of the rock test sample for different diagenetic events comprises:
detecting the influence degree of each component in the rock test sample on the permeability and the porosity of the rock test sample under the action of different diagenetic rocks through a rock physical test, and obtaining the influence value corresponding to each component;
and sequencing the influence values, and obtaining rock component factors corresponding to the influence values according to sequencing results.
4. The method of quantitatively identifying tight sandstone diagenetic according to claim 1, wherein establishing a fit relationship equation between the rock composition factor, the physical property parameter, and the elastic parameter for different diagenetic events according to the rock composition factor and the correlation equation corresponding to the elastic parameter comprises:
when diagenetic is cemented and the rock composition factor is mud siliceous, the fitted relation equation includes:
V p ＝5.99-0.039·V clay -0.064·φ；
V s ＝3.71-0.037·V clay -0.037·φ；
in the above formula, V p For longitudinal wave velocity, V s For transverse wave velocity, V clay For the content of mud-siliceous in the rock test sample,is porosity.
5. The method of quantitatively identifying tight sandstone diagenetic according to any of claims 1 to 4, wherein said diagenetic comprises any of compaction, cementing, cross-talk, erosion.
6. The device for quantitatively identifying the diagenetic effect of the compact sandstone is characterized by comprising a detection module, a calculation module, an analysis module and an identification module;
the detection module is used for acquiring a plurality of rock test samples under different diagenetic effects and respectively detecting physical parameters and rock component factors of the rock test samples;
the calculation module is used for calculating the elasticity parameters of the rock test samples, and establishing fitting relation equations among the rock component factors, the physical parameters and the elasticity parameters under different diagenetic effects through the rock component factors, the physical parameters and the elasticity parameters of each rock test sample;
the analysis module is used for obtaining elastic parameters of the region to be tested through pre-stack seismic inversion, and obtaining physical parameters and rock component factors of the region to be tested through pre-stack elastic parameter inversion;
the identification module is used for obtaining the diagenetic effect of the region to be tested through the fitting relation equation according to the rock composition factors, the elastic parameters and the physical parameters of the region to be tested;
establishing a fitting relation equation among the rock composition factor, the physical property parameter and the elastic parameter under different diagenetic effects through the rock composition factor, the physical property parameter and the elastic parameter of each rock test sample comprises: constructing an association equation of the physical property parameter and the elastic parameter of the rock test sample under different diagenetic effects according to the physical property parameter and the elastic parameter; establishing a fitting relation equation among the rock composition factors, the physical parameters and the elastic parameters under different diagenetic effects according to the rock composition factors and the association equation corresponding to the elastic parameters;
constructing an association equation of the physical property parameter and the elastic parameter of the rock test sample under different diagenetic effects according to the physical property parameter and the elastic parameter comprises: when diagenetic is cemented and the rock composition factors are mud siliceous and calcareous, the correlation equation includes:
when the rock component factors are mud siliceous:
when the rock composition factor is calcareous:
in the above formula, V p For longitudinal wave velocity, V s In order to be a transverse wave velocity,is porosity.
7. An electronic 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 method of any of claims 1 to 5 when executing the computer program.
8. A computer readable storage medium, characterized in that the computer readable storage medium stores a computer program for executing the method of any one of claims 1 to 5.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010902086.2A CN114113537B (en) | 2020-09-01 | 2020-09-01 | Method and device for quantitatively identifying diagenetic effect of compact sandstone |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010902086.2A CN114113537B (en) | 2020-09-01 | 2020-09-01 | Method and device for quantitatively identifying diagenetic effect of compact sandstone |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN114113537A CN114113537A (en) | 2022-03-01 |
| CN114113537B true CN114113537B (en) | 2024-03-26 |
Family
ID=80360273
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010902086.2A Active CN114113537B (en) | 2020-09-01 | 2020-09-01 | Method and device for quantitatively identifying diagenetic effect of compact sandstone |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114113537B (en) |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104516017A (en) * | 2013-09-29 | 2015-04-15 | 中国石油化工股份有限公司 | Carbonate rock physical parameter seismic inversion method |
| CN105467466A (en) * | 2015-12-07 | 2016-04-06 | 中国石油大学(华东) | A tight reservoir lithogenous phase prediction method based on multi-scale information |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP3256885B1 (en) * | 2015-02-13 | 2020-12-30 | Services Petroliers Schlumberger | Diagenetic and depositional rock analysis |
-
2020
- 2020-09-01 CN CN202010902086.2A patent/CN114113537B/en active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104516017A (en) * | 2013-09-29 | 2015-04-15 | 中国石油化工股份有限公司 | Carbonate rock physical parameter seismic inversion method |
| CN105467466A (en) * | 2015-12-07 | 2016-04-06 | 中国石油大学(华东) | A tight reservoir lithogenous phase prediction method based on multi-scale information |
Non-Patent Citations (3)
| Title |
|---|
| 成岩作用和岩石结构对砂岩弹性速度的影响;马中高;《石油学报》;第29卷(第1期);第58-63页 * |
| 渤海湾盆地东营凹陷砂岩储层地震特征及对岩石结构的指示;邓继新;王尚旭;伍向阳;;石油与天然气地质(第03期);第294-299页转第309页 * |
| 邓继新 ; 王尚旭 ; 伍向阳 ; .渤海湾盆地东营凹陷砂岩储层地震特征及对岩石结构的指示.石油与天然气地质.2009,(第03期),第294-299页转第309页. * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN114113537A (en) | 2022-03-01 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US7516016B2 (en) | Method for improving prediction of the viability of potential petroleum reservoirs | |
| Xu et al. | An experimental study on porosity and permeability stress-sensitive behavior of sandstone under hydrostatic compression: characteristics, mechanisms and controlling factors | |
| CN104698500A (en) | Method for predicting reservoir lithogenous phase through geology and logging information | |
| Tutuncu | Anisotropy, compaction and dispersion characteristics of reservoir and seal shales | |
| CA2580570A1 (en) | Integrated anisotropic rock physics model | |
| US9892321B2 (en) | Using maximal inscribed spheres for image-based rock property estimation | |
| Martyushev et al. | The Influence of the rock stress state on matrix and fracture permeability under conditions of various lithofacial zones of the Tournaisian–Fammenian oil fields in the Upper Kama Region | |
| Khan et al. | The effect of geological and geomechanical parameters on reservoir stress path and its importance in studying permeability anisotropy | |
| Zhang et al. | Experimental study of fracturing fluid retention in rough fractures | |
| Spain et al. | Beyond volumetrics: unconventional petrophysics for efficient resource appraisal (example from the Khazzan field, Sultanate of Oman) | |
| Ariadji | Effect of vibration on rock and fluid properties: on seeking the vibroseismic technology mechanisms | |
| Cook et al. | A study of the physical mechanisms of sanding and application to sand production prediction | |
| Holt et al. | Core damage effects on compaction behaviour | |
| CN114113537B (en) | Method and device for quantitatively identifying diagenetic effect of compact sandstone | |
| Vega Navarro | Closure of natural fractures caused by increased effective stress, a case study: reservoir Robore III, Bulo Bulo field, Bolivia | |
| Hamid et al. | Stress dependent permeability of carbonate rock | |
| CN105301642B (en) | Non-homogeneous porous rocks and its solid organic matter volume content determines method and device | |
| Grützner et al. | Physical property changes as a monitor of pelagic carbonate diagenesis: an empirically derived diagenetic model for Atlantic Ocean basins | |
| Han et al. | Velocities of deepwater reservoir sands | |
| Li et al. | A new experimental approach for hydraulic fracturing fluid optimization: Especially focus on the ultra-deep tight gas formation | |
| Al-Ameri et al. | The impact of the fracturing fluid spontaneous imbibition on the petrophysical properties of low permeability reservoirs | |
| CN113622901B (en) | Shale oil reservoir quality analysis method, device, equipment and medium | |
| Zhang et al. | Experimental study on deterioration of bedrock strength and P-wave velocity by pressurized seawater | |
| Fang et al. | Study on origin, classification, and identification of complex porous carbonate reservoir: a case study of the middle cretaceous Mishrif Formation in W Oil field | |
| CN114859407A (en) | Method and device for determining volcanic reservoir acoustic characteristic parameters |
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 |