CN117807900A - Method and device for calculating dry rock modulus of deep carbonate stratum - Google Patents
Method and device for calculating dry rock modulus of deep carbonate stratum Download PDFInfo
- Publication number
- CN117807900A CN117807900A CN202211211549.6A CN202211211549A CN117807900A CN 117807900 A CN117807900 A CN 117807900A CN 202211211549 A CN202211211549 A CN 202211211549A CN 117807900 A CN117807900 A CN 117807900A
- Authority
- CN
- China
- Prior art keywords
- rock
- modulus
- density
- dry rock
- pressure
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 239000011435 rock Substances 0.000 title claims abstract description 269
- BVKZGUZCCUSVTD-UHFFFAOYSA-L Carbonate Chemical compound [O-]C([O-])=O BVKZGUZCCUSVTD-UHFFFAOYSA-L 0.000 title claims abstract description 70
- 238000000034 method Methods 0.000 title claims abstract description 48
- 239000011159 matrix material Substances 0.000 claims abstract description 65
- 229910052500 inorganic mineral Inorganic materials 0.000 claims abstract description 60
- 239000011707 mineral Substances 0.000 claims abstract description 60
- 239000011148 porous material Substances 0.000 claims abstract description 17
- 238000004364 calculation method Methods 0.000 claims abstract description 10
- 230000015572 biosynthetic process Effects 0.000 claims description 29
- 238000004590 computer program Methods 0.000 claims description 13
- 238000004458 analytical method Methods 0.000 claims description 10
- 238000003860 storage Methods 0.000 claims description 9
- 238000012937 correction Methods 0.000 claims description 5
- 238000002441 X-ray diffraction Methods 0.000 claims description 4
- 238000002591 computed tomography Methods 0.000 claims description 3
- 238000007405 data analysis Methods 0.000 claims description 3
- 239000012530 fluid Substances 0.000 description 18
- 238000012545 processing Methods 0.000 description 10
- 238000010586 diagram Methods 0.000 description 7
- 230000006870 function Effects 0.000 description 7
- 230000008569 process Effects 0.000 description 6
- 230000009471 action Effects 0.000 description 4
- 230000008859 change Effects 0.000 description 4
- 238000004891 communication Methods 0.000 description 4
- 238000002474 experimental method Methods 0.000 description 4
- 238000002637 fluid replacement therapy Methods 0.000 description 4
- 230000002093 peripheral effect Effects 0.000 description 4
- 238000007596 consolidation process Methods 0.000 description 3
- 238000005553 drilling Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 238000012360 testing method Methods 0.000 description 3
- 229910021532 Calcite Inorganic materials 0.000 description 2
- 239000004927 clay Substances 0.000 description 2
- 229910052570 clay Inorganic materials 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 239000010459 dolomite Substances 0.000 description 2
- 229910000514 dolomite Inorganic materials 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000003287 optical effect Effects 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
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 2
- 238000012935 Averaging Methods 0.000 description 1
- VTYYLEPIZMXCLO-UHFFFAOYSA-L Calcium carbonate Chemical compound [Ca+2].[O-]C([O-])=O VTYYLEPIZMXCLO-UHFFFAOYSA-L 0.000 description 1
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 description 1
- 238000003491 array Methods 0.000 description 1
- 238000003556 assay Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 230000001413 cellular effect Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 239000007795 chemical reaction product Substances 0.000 description 1
- 150000001875 compounds Chemical class 0.000 description 1
- 230000006835 compression Effects 0.000 description 1
- 238000007906 compression Methods 0.000 description 1
- 239000000470 constituent Substances 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 238000004836 empirical method Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000013213 extrapolation Methods 0.000 description 1
- 229910021389 graphene Inorganic materials 0.000 description 1
- 239000010440 gypsum Substances 0.000 description 1
- 229910052602 gypsum Inorganic materials 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 238000011545 laboratory measurement Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000011022 operating instruction Methods 0.000 description 1
- 230000000704 physical effect Effects 0.000 description 1
- 239000000047 product Substances 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- MWNQXXOSWHCCOZ-UHFFFAOYSA-L sodium;oxido carbonate Chemical compound [Na+].[O-]OC([O-])=O MWNQXXOSWHCCOZ-UHFFFAOYSA-L 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 238000000547 structure data Methods 0.000 description 1
- 239000000758 substrate Substances 0.000 description 1
- 239000011800 void material Substances 0.000 description 1
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention provides a calculation method and a device for the dry rock modulus of a deep carbonate stratum, wherein the method comprises the steps of obtaining the matrix mineral component of a target stratum rock, the volume content, the porosity and the pore structure parameters, the depth of the target stratum and the corresponding pressure; determining the bulk modulus, shear modulus and density of a rock matrix skeleton according to the mineral components and the bulk modulus, shear modulus and density of the mineral of the rock matrix of the target layer; selecting a theoretical model of carbonate rock suitable for the depth of the target layer based on the rock characteristics, and calculating the intrinsic bulk modulus, shear modulus and dry rock density of the dry rock; according to rock physical experimental data of the target layer rock sample, establishing a relation between dry rock modulus and pressure and a relation between dry rock density and pressure; substituting the intrinsic bulk modulus, the shear modulus and the dry rock density of the dry rock into a relation between the dry rock modulus and the pressure, correcting the relation between the dry rock density and the pressure, and combining the target layer depth and the corresponding pressure to obtain the carbonate stratum dry rock modulus and the carbonate stratum density of the target layer depth.
Description
Technical Field
The invention relates to the technical field of rock physics, in particular to a method and a device for calculating the dry rock modulus of a deep carbonate stratum, a computer readable storage medium and electronic equipment.
Background
In oil and gas exploration and development, the type, property and distribution of fluid are predicted by using seismic data, and the prediction basis is that the elastic parameters and the seismic attributes caused by fluid change have obvious differences and can be identified. After the reservoir layer obtains the parameters of physical property parameters, fluid (oil, gas and water) saturation, fluid density, stratum lithology and the like which are actually measured, the reservoir layer porosity and other parameters such as rock skeleton and the like are assumed to be unchanged, and the elasticity parameters are calculated by utilizing a rock physical model through changing the water and oil (gas) saturation, so that the corresponding elasticity parameters and the change of earthquake attributes caused by different fluids or fluid saturation in the same stratum are simulated, and a basis is established for earthquake reservoir layer prediction and fluid identification.
Extrapolation from known well data analyzes seismic signature changes caused by fluid and pressure changes, and many times the predicted results are far from actual because there is a very important set of parameters in the petrophysical model used for flow replacement analysis-dry rock modulus, the exact calculation of which is the difficulty and key to fluid replacement analysis. The dry rock modulus can be obtained indirectly by measuring the longitudinal and transverse wave velocity and density of the dry rock sample in a laboratory, but is more calculated by using an empirical formula or an equivalent medium model.
Laboratory measurement methods require the collection of underground core samples, which are not only costly, but also difficult to do in practice with dense sampling. At present, empirical methods including a Biot coefficient method, a prism consolidation parameter method and a Nur critical porosity method are commonly adopted in fluid replacement analysis, and the size of the dry rock modulus is determined by one parameter. The equivalent medium model method comprises the steps of obtaining the modulus of dry rock by a K-T model, a DEM model, an SC model and the like, wherein thinking is basically similar, and can be summarized in that a substrate is taken as a background medium, inclusions are all based on idealized ellipsoidal inclusion shapes, the inclusions are isolated from each other, and the inclusions are added into the background medium to form the equivalent medium. The essence is mainly calculated by relation to the porosity of the rock, the pore shape or the consolidation parameters of the rock. It is well known that a subterranean reservoir is under a pressure environment, and whether the formation porosity, pore shape or consolidated rock changes, the dry rock bulk modulus is affected by the pressure, but the existing dry rock modulus and density calculation methods do not take the pressure effect into account.
Disclosure of Invention
In view of the foregoing, embodiments of the present invention provide a method, an apparatus, a computer-readable storage medium, and an electronic device for calculating a dry rock modulus of a deep carbonate formation.
In a first aspect, an embodiment of the present invention provides a method for calculating a dry rock modulus of a deep carbonate formation, including:
s100, obtaining a target layer rock matrix mineral component, volume content, porosity and pore structure parameters, a target layer depth and corresponding pressure according to target layer rock data;
s200, determining the bulk modulus, shear modulus and density of a rock matrix skeleton according to the mineral components and the volume content of the rock matrix of the target layer;
s300, establishing a carbonate rock dry rock model based on rock characteristics, selecting a theoretical model of carbonate rock suitable for a target depth of layer, and calculating intrinsic bulk modulus, shear modulus and dry rock density of the dry rock;
s400, establishing a relation between the modulus of dry rock and the pressure and a relation between the density of dry rock and the pressure according to rock physical experimental data of a target layer rock sample;
s500, substituting the intrinsic bulk modulus, the shear modulus and the dry rock density of the dry rock into a relation between the dry rock modulus and the pressure and a relation between the dry rock density and the pressure, correcting, and combining the target depth of layer and the corresponding pressure to obtain the dry rock modulus and the density of the carbonate rock stratum of the target depth of layer.
According to an embodiment of the present invention, step S100 includes: carrying out geological slice identification, X-ray diffraction analysis and CT scanning according to a carbonate reservoir core sample, and determining the mineral components, the volume content and the pore structure of a target reservoir rock matrix; the volume content of the mineral components of the target layer rock matrix, the rock porosity parameters and the pore structure parameters are obtained through logging data, and the depth of the target layer and the corresponding overburden formation pressure are obtained by combining well completion data.
According to an embodiment of the present invention, step S200 includes: according to the mineral components of the rock matrix of the target layer, the volume content of the mineral components, the volume modulus and the shear modulus of the mineral, calculating the volume modulus and the shear modulus of the skeleton of the rock matrix of the target layer by utilizing a V-R-H average formula; the density of the target layer rock matrix is calculated according to the weighted average of the mineral components according to the mineral components of the target layer rock matrix, the volume content and the mineral density.
According to an embodiment of the present invention, step S300 includes: taking a rock matrix skeleton and the porosity as a whole, and establishing a carbonate rock dry rock model; a theoretical model of carbonate rock suitable for the depth of layer of interest is selected and the intrinsic bulk modulus, shear modulus and dry rock density of the dry rock are calculated.
According to the embodiment of the invention, the theoretical model of the carbonate rock suitable for the target depth of layer comprises a K-T model or a self-consistent theoretical model based on a first-order scattering theory.
According to an embodiment of the present invention, step S400 includes: and establishing a relation between the low-pressure dry rock modulus and the pressure and a relation between the high-pressure dry rock modulus and the pressure according to petrophysical experimental data of the carbonate dry rock sample with the target depth of layer under the pressure.
According to an embodiment of the invention, the low pressure and the high pressure are distinguished by a critical pressure, which is determined according to the formation conditions of the layer of interest.
In a second aspect, the present invention also provides a device for calculating the dry rock modulus of a deep carbonate formation, which is characterized by comprising:
the rock data analysis module is used for obtaining the mineral components of the target layer rock matrix, the volume content, the porosity and the pore structure parameters, the target layer depth and the corresponding pressure according to the target layer rock data;
the modulus density determining module is used for determining the bulk modulus, the shear modulus and the density of the rock matrix skeleton according to the mineral components and the volume content of the rock matrix of the target layer;
the intrinsic modulus determining module is used for establishing a carbonate rock dry rock model based on a rock matrix, selecting a theoretical model of carbonate rock suitable for the depth of a target layer, and calculating the intrinsic bulk modulus, the shear modulus and the dry rock density of the dry rock;
the pressure relation analysis module is used for establishing a relation between the dry rock modulus and the pressure and a relation between the dry rock density and the pressure according to rock physical experimental data of the target layer rock sample;
and the modulus density correction module is used for substituting the intrinsic bulk modulus, the shear modulus and the dry rock density of the dry rock into a relation between the dry rock modulus and the pressure and a relation between the dry rock density and the pressure, and correcting the relation between the dry rock density and the pressure, and combining the target layer depth and the corresponding pressure to obtain the dry rock modulus and the density of the carbonate rock stratum with the target layer depth.
In a third aspect, an embodiment of the present invention provides a computer readable storage medium having stored thereon a computer program which, when executed by a processor, implements a method for calculating a dry rock modulus of a deep carbonate formation as described in the first aspect.
In a fourth aspect, an embodiment of the present invention provides an electronic device, including:
a processor;
a memory for storing the processor-executable instructions;
wherein the processor is configured to execute the instructions to implement a method of calculating a dry rock modulus of a deep carbonate formation as described in the first aspect above.
Compared with the prior art, the technical scheme of the invention has the following beneficial effects:
compared with the prior art, the method has the innovation point that the relation between the rock skeleton density and the modulus and the stratum pressure is introduced, so that the calculation of the carbonate rock dry rock modulus under the action of the pressure in a deep environment is more in line with the actual rock, and the dry modulus is applied to fluid replacement analysis of different fluids or fluid saturation and analysis of response of the fluids to earthquake attributes, so that the accuracy of simulating corresponding elastic parameters and earthquake attribute variation caused by different fluids or fluid saturation in the same stratum can be improved, and the risk of earthquake fluid identification is reduced.
Drawings
In order to more clearly illustrate the embodiments of the invention or the technical solutions in the prior art, the drawings that are required in the embodiments or the description of the prior art will be briefly described, it being obvious that the drawings in the following description are only some embodiments of the invention, and that other drawings may be obtained according to these drawings without inventive effort for a person skilled in the art.
FIG. 1 is a flow chart of steps of a method for calculating the dry rock modulus of a deep carbonate formation provided by an embodiment of the present invention;
FIG. 2 is a graph of the porosity, density, bulk modulus of dry rock and shear modulus of a carbonate rock of the destination layer provided by an embodiment of the present invention;
fig. 3 is a schematic diagram of an electronic device according to an embodiment of the present invention.
Detailed Description
For the purpose of making the objects, technical solutions and advantages of the embodiments of the present invention more apparent, the technical solutions of 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, and it is apparent that the described embodiments are some embodiments of the present invention, but not all embodiments of the present invention. All other embodiments, which can be made by those skilled in the art based on the embodiments of the invention without making any inventive effort, are intended to be within the scope of the invention.
Example 1
The invention aims to solve the problems that the influence of formation pressure on elastic parameters of rock is not considered and the calculation accuracy of dry rock modulus is poor in the prior art, and provides a novel calculation method for the dry rock modulus and the formation pressure based on a high-pressure rock physical experiment so as to improve the accuracy of simulating corresponding elastic parameters and seismic attribute variation caused by different fluids or fluid saturation in the same stratum. In order to achieve the above object, the present embodiment provides the following method for calculating the dry rock modulus of a deep carbonate formation. The method mainly comprises the following steps:
and step 1, acquiring basic data of a deep carbonate reservoir.
(1) And collecting a target layer rock/wall core sample, carrying out rock physical experiment tests such as sheet identification, X-ray diffraction analysis and the like, obtaining rock matrix mineral composition and pore structure data, and calibrating logging by using core data.
(2) And obtaining basic parameters of the rock through logging data. The basic rock parameters include matrix mineral composition and its volume content, porosity and pore structure parameters, depth and pressure. Wherein the depth H of the formation of interest and its corresponding pressure P can be obtained in connection with the drilling and completion data.
Carbonate matrix minerals may include calcite, dolomite, clay, aragonite, quartz, gypsum, and the like.
In particular, porosity refers to the sum of all void space volumes in a unit reservoir rock sample, expressed as a percentage. The rock matrix mineral component refers to the primary mineral in the carbonate reservoir rock sample, such as the relatively more abundant mineral. Rock matrix mineral volume content refers to the volume percent of matrix mineral per carbonate reservoir rock.
And 2, determining the skeleton modulus and the density of the rock matrix corresponding to the carbonate reservoir according to the matrix mineral components and the volume content.
Based on matrix mineral composition and volume content in the base data and its volume modelThe volume modulus K of the matrix skeleton of the carbonate rock is calculated and determined by the amount, the shear modulus and the density m Shear modulus mu m And density ρ m 。
Specifically, the bulk modulus K of the rock matrix framework is calculated by utilizing a V-R-H average formula m Shear modulus mu m . Calculating the density ρ of the rock matrix skeleton according to the weighted average of the mineral components m Wherein the volume content of the various minerals has been obtained in the basic data.
V-R-H average:
wherein f i 、K i Sum mu i The i-th mineral component accounts for the volume content, the volume modulus and the shear modulus of the total mineral, and N is the number of the mineral components.
Density ρ of matrix backbone m Calculated on the basis of a weighted average of the mineral constituents of the matrix, i.e
ρ i For the density of the i-th mineral component, f i The i-th mineral component accounts for the volume content of the total mineral.
And 3, building a carbonate rock dry rock model on the basis of the rock matrix. According to the rock characteristics, selecting a theoretical model suitable for deep carbonate rock, and calculating the intrinsic bulk modulus of the dry rockShear modulus>And dry rock Density->
Specifically, after the rock matrix model is established, a dry rock model of the carbonate reservoir is established with the rock matrix skeleton and the porosity as a whole. The K-T model based on the first-order scattering theory or the self-consistent theoretical model can be used for calculating the intrinsic bulk modulus of the dry rockShear modulus>Dry rock density->Is the average of matrix backbone and porosity.
Wherein i represents the ith inclusion compound and the corresponding volume content is x i And:
coefficient P *i And Q *i Is an influencing factor describing the effect after the inclusion material i has been added to the background medium m.For dry rock bulk modulus>Is dryShear modulus of rock, K m Is the bulk modulus, mu, of the rock matrix skeleton m Is the shear modulus of the rock matrix skeleton.
Dry rock densityObtained by a weighted average of rock matrix skeletal density and porosity:
ρ m for the matrix skeleton density, the matrix skeleton density is equal to the matrix skeleton density,is the porosity ρ p The dry cell density may be approximately 0.
Step 4, establishing low-pressure and high-pressure dry rock modulus K according to rock physical experiment test data of the dry rock sample of the deep carbonate reservoir under the pressure change d 、μ d Relationship to pressure P:
K d (P)＝b+cln(P) P≤P c (6)
μ d (P)＝x+yln(P) P≤P c
wherein,the intrinsic modulus of rock at 1 atmosphere pressure, a, b, c, d, x, y, is the fitting coefficient. P is the pressure, P c For critical pressure, 55-80Mpa is generally selected, as determined by the formation conditions.
Establishing a relation between dry rock density and pressure:
ρ d (P) is the density of dry rock in the deep layer of the underground, in g/cm 3 ，Is the intrinsic density of the rock framework under the normal pressure of the earth surface, and the unit g/cm 3 B is the average body compression coefficient of the rock framework, and the unit is MPa -1 Obtained from experimental assays. P is the static pressure of the stratum at the underground depth (H meters) and is in MPa.
Step 5, utilizing the established relation between the dry rock modulus, density and differential pressure to calculate the dry rock modulusAnd Density->Correcting to obtain the depth (H) carbonate reservoir dry rock modulus K d (H)、μ d (H) And density ρ d (H)。
Wherein the pressure P is converted to a depth H representation according to the formation depth and pressure relationship. The formulas (6) - (8) are rewritten as dry rock modulus and density as a function of depth H:
dry rock modulus correction formula:
K d(H)＝bln(P)+c P≤Pc
μ d (H)＝x+yln(P) P≤Pc
the correction formula of the rock skeleton density:
where depth H is a variable.
By adopting the method, the equivalent effect of the influence of the pressure of the overburden stratum on the porosity, different pore shapes and rock consolidation degree of the deep carbonate rock can be calculated, the condition that the bulk modulus parameter of the dry rock is more close to the condition under the deep environment is obtained, and the accuracy of fluid replacement analysis and the influence of the fluid on the seismic attribute calculated based on the rock physical model is improved. The method is widely applied to the field of oil and gas exploration and development, comprises reservoir quantitative description, fluid prediction, analysis of reservoir description and monitoring processes, time-lapse seismic feasibility analysis and the like, and can also be applied to research and production of concealed oil and gas reservoirs, unconventional oil and gas and the like.
Example two
As shown in fig. 1, an embodiment of the present invention provides a method for calculating a dry rock modulus of a deep carbonate reservoir, including:
and step 1, performing geological slice identification, X-ray diffraction analysis and CT scanning according to a carbonate reservoir rock/drilling wall core sample, and determining the mineral components and pore structures of a target layer rock matrix. The main mineral components include dolomite, calcite, clay, quartz, etc.
Obtaining the volume content f of main mineral components at each depth point of a target layer through logging data i Porosity parameter of rockAnd combining with the pore structure parameters and the drilling and completion data to obtain the depth H of the target layer of 6440-6539 m and the corresponding overburden stratum pressure P of 142.1-144.3Mpa.
Step 2, obtaining the elastic modulus K of each mineral of the reservoir according to the matrix mineral components and the volume content of each depth of the target layer i 、μ i Mineral density ρ i 。
According to VRH method, i.e. formulaAnd->Calculating K of matrix skeleton m ,μ m 。
According to the volume content f of each mineral in the reservoir i And the respective mineral density ρ i By volume weighted averagingCalculating the skeleton density rho of rock matrix m 。
Step 3, after the rock matrix model is built, taking the rock matrix skeleton and the porosity as a whole, building a dry rock model of a carbonate reservoir, and calculating the intrinsic bulk modulus of the dry rock by adopting a K-T model based on a first-order scattering theoryShear modulus>And dry rock Density->
Step 4, determining low-pressure and high-pressure dry rock modulus K according to rock physical experiment test data of the deep carbonate reservoir dry rock sample under the pressure change d 、μ d Coefficient a, b, c, d, x, y in relation to pressure P and dry rock density ρ d And coefficient B in the pressure relationship.
Step 5, bulk modulus of the dry rockAnd shear modulus->Substituting the formulas (9) - (11) to correct the pressure to obtain the dry rock modulus K of the deep carbonate reservoir d (H)、μ d (H)And density ρ d (H) The conversion of depth H and pressure P (as shown in fig. 2) is required during implementation.
The foregoing description of the embodiments of the present disclosure has been presented for purposes of illustration and description, and is not intended to be exhaustive or limited to the embodiments disclosed. Many modifications and variations will be apparent to those of ordinary skill in the art without departing from the scope and spirit of the various embodiments described. The terminology used herein was chosen in order to best explain the principles of the embodiments, the practical application, or the technical improvements in the marketplace, or to enable others of ordinary skill in the art to understand the embodiments disclosed herein.
Example III
The following are examples of the apparatus of the present invention that may be used to perform the method embodiments of the present invention. For details not disclosed in the embodiments of the apparatus of the present invention, please refer to the embodiments of the method of the present invention.
The embodiment provides a deep carbonate formation dry rock modulus calculation device, which comprises:
the rock data analysis module is used for obtaining the mineral components of the target layer rock matrix, the volume content, the porosity and the pore structure parameters, the target layer depth and the corresponding pressure according to the target layer rock data;
the modulus density determining module is used for determining the bulk modulus, the shear modulus and the density of the rock matrix skeleton according to the mineral components and the volume contents of the mineral components of the rock matrix of the target layer and the bulk modulus, the shear modulus and the density of the mineral;
the intrinsic modulus determining module is used for establishing a carbonate rock dry rock model based on rock characteristics, selecting a theoretical model of carbonate rock suitable for the depth of a target layer, and calculating intrinsic bulk modulus, shear modulus and dry rock density of the dry rock;
the pressure relation analysis module is used for establishing a relation between the dry rock modulus and the pressure and a relation between the dry rock density and the pressure according to rock physical experimental data of the target layer rock sample;
and the modulus density correction module is used for substituting the intrinsic bulk modulus, the shear modulus and the dry rock density of the dry rock into a relation between the dry rock modulus and the pressure and a relation between the dry rock density and the pressure, and correcting the relation between the dry rock density and the pressure, and combining the target layer depth and the corresponding pressure to obtain the dry rock modulus and the density of the carbonate rock stratum with the target layer depth.
Example IV
The present embodiment provides a computer readable medium having stored thereon a computer program which, when executed by a processor, implements the steps of a method for calculating the dry rock modulus of a deep carbonate formation as described in the above embodiments.
It should be noted that, all or part of the flow of the method of the above embodiment may be implemented by a computer program, which may be stored in a computer readable storage medium and which, when executed by a processor, implements the steps of the method embodiments described above. Wherein the computer program comprises computer program code which may be in source code form, object code form, executable file or some intermediate form etc. The computer readable medium may include: any entity or device capable of carrying the computer program code, a recording medium, a U disk, a removable hard disk, a magnetic disk, an optical disk, a computer Memory, a Read-Only Memory (ROM), a random access Memory (RAM, random Access Memory), an electrical carrier signal, a telecommunications signal, a software distribution medium, and so forth. Of course, there are other ways of readable storage medium, such as quantum memory, graphene memory, etc. It should be noted that the computer readable medium contains content that can be appropriately scaled according to the requirements of jurisdictions in which such content is subject to legislation and patent practice, such as in certain jurisdictions in which such content is subject to legislation and patent practice, the computer readable medium does not include electrical carrier signals and telecommunication signals.
Example five
Fig. 3 is a schematic structural view of an electronic device according to an embodiment of the present invention. As shown in fig. 3, at the hardware level, the electronic device comprises a processor, optionally together with an internal bus, a network interface, a memory. The Memory may include a Memory, such as a Random-Access Memory (RAM), and may further include a non-volatile Memory (non-volatile Memory), such as at least 1 disk Memory. Of course, the electronic device may also include hardware required for other services.
The processor, network interface, and memory may be interconnected by an internal bus, which may be an ISA (Industry Standard Architecture ) bus, a PCI (PeripheralComponent Interconnect, peripheral component interconnect standard) bus, or EISA (Extended Industry StandardArchitecture ) bus, among others. The buses may be classified as address buses, data buses, control buses, etc. For ease of illustration, the figures are shown with only line segments, but not with only one bus or one type of bus.
And the memory is used for storing programs. In particular, the program may include program code including computer-operating instructions. The memory may include memory and non-volatile storage and provide instructions and data to the processor. The processor reads the corresponding computer program from the non-volatile memory into the memory and then runs. The processor executes the program stored in the memory to perform all the steps in the above-mentioned method for calculating the dry rock modulus of the deep carbonate formation.
The communication bus mentioned by the above devices may be a peripheral component interconnect standard (Peripheral Component Interconnect, PCI) bus or an extended industry standard architecture (Extended Industry Standard Architecture, EISA) bus, etc. The communication bus may be classified as an address bus, a data bus, a control bus, or the like. For ease of illustration, the figures are shown with only one bold line, but not with only one bus or one type of bus. The communication interface is used for communication between the electronic device and other devices.
The bus includes hardware, software, or both for coupling the above components to each other. For example, the bus may include an Accelerated Graphics Port (AGP) or other graphics bus, an Enhanced Industry Standard Architecture (EISA) bus, a Front Side Bus (FSB), a HyperTransport (HT) interconnect, an Industry Standard Architecture (ISA) bus, an infiniband interconnect, a Low Pin Count (LPC) bus, a memory bus, a micro channel architecture (MCa) bus, a Peripheral Component Interconnect (PCI) bus, a PCI-Express (PCI-X) bus, a Serial Advanced Technology Attachment (SATA) bus, a video electronics standards association local (VLB) bus, or other suitable bus, or a combination of two or more of the above. The bus may include one or more buses, where appropriate. Although embodiments of the invention have been described and illustrated with respect to a particular bus, the invention contemplates any suitable bus or interconnect.
The Memory may include random access Memory (Random Access Memory, RAM) or may include Non-Volatile Memory (NVM), such as at least one disk Memory. Optionally, the memory may also be at least one memory device located remotely from the aforementioned processor.
The memory may include mass storage for data or instructions. By way of example, and not limitation, the memory may comprise a Hard Disk Drive (HDD), floppy Disk Drive, flash memory, optical Disk, magneto-optical Disk, magnetic tape, or universal serial bus (Universal Serial Bus, USB) Drive, or a combination of two or more of the foregoing. The memory may include removable or non-removable (or fixed) media, where appropriate. In a particular embodiment, the memory is a non-volatile solid state memory. In a particular embodiment, the memory includes Read Only Memory (ROM). The ROM may be mask programmed ROM, programmable ROM (PROM), erasable PROM (EPROM), electrically Erasable PROM (EEPROM), electrically rewritable ROM (EAROM), or flash memory, or a combination of two or more of these, where appropriate.
The processor may be a general-purpose processor, including a central processing unit (Central Processing Unit, CPU), a network processor (Network Processor, NP), etc.; but also digital signal processors (Digital Signal Processing, DSP), application specific integrated circuits (Application Specific Integrated Circuit, ASIC), field programmable gate arrays (Field-Programmable Gate Array, FPGA) or other programmable logic devices, discrete gate or transistor logic devices, discrete hardware components.
It should be noted that, for convenience and brevity of description, only the above-described division of the functional units and modules is illustrated, and in practical application, the above-described functional allocation may be performed by different functional units and modules according to needs, that is, the internal structure of the apparatus is divided into different functional units or modules to perform all or part of the functions described above. The functional units and modules in the embodiment may be integrated in one processing unit, or each unit may exist alone physically, or two or more units may be integrated in one unit, where the integrated units may be implemented in a form of hardware or a form of a software functional unit. In addition, the specific names of the functional units and modules are only for distinguishing from each other, and are not used for limiting the protection scope of the present invention. The specific working process of the units and modules in the above system may refer to the corresponding process in the foregoing method embodiment, which is not described herein again.
The apparatus, device, system, module or unit described in the above embodiments may be implemented in particular by a computer chip or entity or by a product having a certain function. One typical implementation is a computer. In particular, the computer may be, for example, a personal computer, a laptop computer, a car-mounted human-computer interaction device, a cellular telephone, a camera phone, a smart phone, a personal digital assistant, a media player, a navigation device, an email device, a game console, a tablet computer, a wearable device, or a combination of any of these devices.
Although the invention provides method operational steps as described in the examples or flowcharts, more or fewer operational steps may be included based on conventional or non-inventive means. The order of steps recited in the embodiments is merely one way of performing the order of steps and does not represent a unique order of execution. When implemented in an actual device or end product, the instructions may be executed sequentially or in parallel (e.g., in a parallel processor or multi-threaded processing environment, or even in a distributed data processing environment) as illustrated by the embodiments or by the figures.
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.
It should be noted that in this document relational terms such as first and second, and the like are used solely to distinguish one entity or action from another entity or action without necessarily requiring or implying any actual such relationship or order between such entities or actions. Moreover, the terms "comprises," "comprising," or any other variation thereof, are intended to cover a non-exclusive inclusion, such that a process, method, article, or apparatus that comprises a list of elements does not include only those elements but may include other elements not expressly listed or inherent to such process, method, article, or apparatus. Without further limitation, an element defined by the phrase "comprising one … …" does not exclude the presence of other like elements in a process, method, article, or apparatus that comprises the element.
In this specification, each embodiment is described in a related manner, and identical and similar parts of each embodiment are all referred to each other, and each embodiment mainly describes differences from other embodiments. In particular, for apparatus, electronic devices, and readable storage medium embodiments, since they are substantially similar to method embodiments, the description is relatively simple, and references to parts of the description of method embodiments are only required.
The foregoing description is only of the preferred embodiments of the present invention and is not intended to limit the scope of the present invention. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the present invention are included in the protection scope of the present invention.
Claims (10)
1. The calculation method of the dry rock modulus of the deep carbonate stratum is characterized by comprising the following steps of:
s100, obtaining a target layer rock matrix mineral component, volume content, porosity and pore structure parameters, a target layer depth and corresponding pressure according to target layer rock data;
s200, determining the bulk modulus, shear modulus and density of a rock matrix skeleton according to the mineral components and the bulk modulus, shear modulus and density of the mineral of the rock matrix of the target layer;
s300, based on rock characteristics, establishing a carbonate rock dry rock model, selecting a theoretical model of carbonate rock suitable for the depth of a target layer, and calculating the intrinsic bulk modulus, shear modulus and dry rock density of the dry rock;
s400, establishing a relation between the modulus of dry rock and the pressure and a relation between the density of dry rock and the pressure according to rock physical experimental data of a target layer rock sample;
s500, substituting the intrinsic bulk modulus, the shear modulus and the dry rock density of the dry rock into a relation between the dry rock modulus and the pressure and a relation between the dry rock density and the pressure, correcting, and combining the target depth of layer and the corresponding pressure to obtain the dry rock modulus and the density of the carbonate rock stratum of the target depth of layer.
2. The method of calculating the dry rock modulus of a deep carbonate formation of claim 1, wherein step S100 comprises:
carrying out geological slice identification, X-ray diffraction analysis and CT scanning according to a carbonate reservoir core sample, and determining a target reservoir rock matrix mineral component and a pore structure;
the volume content of the mineral components of the target layer rock matrix, the rock porosity parameters and the pore structure parameters are obtained through logging data, and the depth of the target layer and the corresponding overburden formation pressure are obtained by combining well completion data.
3. The method of calculating the dry rock modulus of a deep carbonate formation of claim 1, wherein step S200 comprises:
calculating the bulk modulus and the shear modulus of the matrix skeleton of the rock of the target layer by utilizing a V-R-H average formula according to the mineral components and the bulk modulus and the shear modulus of the mineral of the matrix of the rock of the target layer;
and calculating the density of the matrix skeleton of the rock matrix of the target layer according to the weighted average of the mineral components according to the mineral components of the matrix of the rock of the target layer, the volume content and the mineral density of the matrix of the rock of the target layer.
4. The method of calculating the dry rock modulus of a deep carbonate formation of claim 1, wherein step S300 comprises:
taking a rock matrix skeleton and the porosity as a whole, and establishing a carbonate rock dry rock model;
a theoretical model of carbonate rock suitable for the depth of layer of interest is selected and the intrinsic bulk modulus, shear modulus and dry rock density of the dry rock are calculated.
5. The method for calculating the dry rock modulus of a deep carbonate formation of claim 4,
the theoretical model of the carbonate rock suitable for the target depth of layer comprises a K-T model or a self-consistent theoretical model based on a first-order scattering theory.
6. The method of calculating the dry rock modulus of a deep carbonate formation of claim 1, wherein step S400 comprises:
and establishing a relation between the low-pressure dry rock modulus and the pressure and a relation between the high-pressure dry rock modulus and the pressure according to petrophysical experimental data of the carbonate dry rock sample with the target depth of layer under the pressure.
7. The method for calculating the dry rock modulus of a deep carbonate formation of claim 6,
the low pressure and the high pressure are distinguished by a critical pressure, which is determined based on the formation conditions of the layer of interest.
8. A deep carbonate formation dry rock modulus calculation device, comprising:
the rock data analysis module is used for obtaining the mineral components of the target layer rock matrix, the volume content, the porosity and the pore structure parameters, the target layer depth and the corresponding pressure according to the target layer rock data;
the modulus density determining module is used for determining the bulk modulus, the shear modulus and the density of the rock matrix skeleton according to the mineral components and the volume contents of the mineral components of the rock matrix of the target layer and the bulk modulus, the shear modulus and the density of the mineral;
the intrinsic modulus determining module is used for establishing a carbonate dry rock model based on rock characteristics, selecting a theoretical model of carbonate rock suitable for the depth of a target layer, and calculating intrinsic bulk modulus, shear modulus and dry rock density of the dry rock;
the pressure relation analysis module is used for establishing a relation between the dry rock modulus and the pressure and a relation between the dry rock density and the pressure according to rock physical experimental data of the target layer rock sample;
and the modulus density correction module is used for substituting the intrinsic bulk modulus, the shear modulus and the dry rock density of the dry rock into a relation between the dry rock modulus and the pressure and a relation between the dry rock density and the pressure, and correcting the relation between the dry rock density and the pressure, and combining the target layer depth and the corresponding pressure to obtain the dry rock modulus and the density of the carbonate rock stratum with the target layer depth.
9. A computer readable storage medium, characterized in that it has stored thereon a computer program which, when executed by a processor, implements a method for calculating the dry rock modulus of a deep carbonate formation according to any one of claims 1 to 7.
10. An electronic device, comprising:
a processor;
a memory for storing the processor-executable instructions;
wherein the processor is configured to execute the instructions to implement a deep carbonate formation dry rock modulus calculation method as claimed in any one of claims 1 to 7.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202211211549.6A CN117807900A (en) | 2022-09-30 | 2022-09-30 | Method and device for calculating dry rock modulus of deep carbonate stratum |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202211211549.6A CN117807900A (en) | 2022-09-30 | 2022-09-30 | Method and device for calculating dry rock modulus of deep carbonate stratum |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN117807900A true CN117807900A (en) | 2024-04-02 |
Family
ID=90420667
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202211211549.6A Pending CN117807900A (en) | 2022-09-30 | 2022-09-30 | Method and device for calculating dry rock modulus of deep carbonate stratum |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN117807900A (en) |
-
2022
- 2022-09-30 CN CN202211211549.6A patent/CN117807900A/en active Pending
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Pinheiro et al. | Geostatistical simulation to map the spatial heterogeneity of geomechanical parameters: A case study with rock mass rating | |
| US8583410B2 (en) | Method for obtaining consistent and integrated physical properties of porous media | |
| CN110318745B (en) | Particle size lithology logging evaluation method under deposition microphase constraint | |
| CN111077568B (en) | Method and equipment for detecting oil and gas reservoir by fluid factor of tight oil and gas reservoir | |
| US8359184B2 (en) | Method, program and computer system for scaling hydrocarbon reservoir model data | |
| Chehrazi et al. | Pore-facies as a tool for incorporation of small-scale dynamic information in integrated reservoir studies | |
| RU2731842C1 (en) | Methods and systems for determining bulk density, porosity and pore size distribution of a subsurface formation | |
| CN109709301B (en) | Method, device and system for determining distribution of fracture pore type tight sandstone reservoir | |
| Rezaei et al. | Studying the correlation of rock properties with P-wave velocity index in dry and saturated conditions | |
| US9892321B2 (en) | Using maximal inscribed spheres for image-based rock property estimation | |
| Karaman et al. | A comparative assessment of rock mass deformation modulus | |
| Aljuboori et al. | The impact of diagenesis precipitation on fracture permeability in naturally fractured carbonate reservoirs | |
| CN104316962A (en) | Rock pore structure parameter inversion method and system | |
| CN109323954A (en) | A kind of predicting method of formation pore pressure for car-bonate rock | |
| Gong et al. | Quantification of fracture surface roughness and its insights to mechanical rock properties determination using image analysis techniques | |
| CN115522912A (en) | Cable adsorption card risk quantitative evaluation method, system and medium based on big data | |
| Shchipanov et al. | A new approach to deformable fractured reservoir characterization: case study of the Ekofisk field | |
| CN117807900A (en) | Method and device for calculating dry rock modulus of deep carbonate stratum | |
| Danquigny et al. | Intra-and inter-facies variability of multi-physics data in carbonates. New insights from database of ALBION R&D project | |
| CN109580906B (en) | Method and system for manufacturing shale brittleness identification chart based on rock physics | |
| CN116930023A (en) | Fine interpretation method and device for dense sandstone phase-control classified porosity logging | |
| Chaudhary et al. | Physico-mechanical characteristics of vindhyan sandstone, India | |
| CN114153002A (en) | Three-dimensional geological modeling method and device for natural fracture of reservoir, electronic equipment and medium | |
| Zhang et al. | Characteristics, Main Controlling Factors and Prediction of Complex Fracture–Cavity Reservoirs of Buried Carbonate Hills in the Weixinan Depression, Western South China Sea | |
| CN117077446B (en) | Unconventional reservoir comprehensive brittleness index evaluation method and application system |
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 |