Method and device for generating elastic parameters of organic nanoporous shale reservoir

Abstract

The present invention provides a method and device for generating elastic parameters of a shale reservoir containing organic nanopores, which relates to the field of geophysical technology. The method comprises: obtaining target measurement parameters and physical morphology parameters of organic matter in a shale reservoir containing organic nanopores; generating organic elastic property parameters containing nanopores according to the physical morphology parameters of organic matter; generating inorganic elastic stiffness coefficients according to target measurement parameters; generating coupled elastic property parameters of shale oil and gas reservoirs according to the organic elastic property parameters and the inorganic elastic stiffness coefficients; generating elastic parameters of shale oil and gas reservoirs according to the coupled elastic property parameters of shale oil and gas reservoirs. The present invention can generate elastic parameters that take into account the influence of organic nanopores, and the obtained elastic parameters can be used to improve the accuracy of prediction of porosity and related physical property parameters of shale oil and gas reservoirs.

Description

Technical Field

[0001] The present invention relates to the field of geophysical technology, and in particular to a method and device for generating elastic parameters of a shale reservoir containing organic nanopores. Background Art

[0002] The mineral composition of shale oil and gas reservoirs is complex, the pore structure is complex, and a large number of micro-nano pores are developed. Nano-scale pores are the main storage space for shale oil and gas. Since the high specific surface area of ​​nano-scale pores can cause a strong surface energy effect, the existing methods of predicting the porosity of shale reservoirs using logging and seismic elastic characteristics are often based on the equivalent medium theory of classical elastic mechanics, without considering the influence of the difference in surface and interface effects caused by different nanopore sizes on the elastic characteristics, which leads to the existing prediction of the porosity and related physical properties of shale oil and gas reservoirs is not accurate enough. Summary of the invention

[0003] The present invention provides a method and device for generating elastic parameters of a shale reservoir containing organic nanopores, which can generate elastic parameters that take into account the influence of organic nanopores. Therefore, the obtained elastic parameters can be used to improve the accuracy of predicting the porosity and related physical property parameters of shale oil and gas reservoirs.

[0004] In the first aspect, an embodiment of the present invention provides a method for generating elastic parameters of a shale reservoir containing organic nanopores, the method comprising: obtaining target measurement parameters and physical morphology parameters of organic matter of the shale reservoir containing organic nanopores; generating organic elastic property parameters containing nanopores according to the physical morphology parameters of organic matter; generating inorganic elastic stiffness coefficients according to the target measurement parameters; generating coupled elastic property parameters of shale oil and gas reservoirs according to the organic elastic property parameters and the inorganic elastic stiffness coefficients; generating shale oil and gas reservoir elastic parameters according to the coupled elastic property parameters of the shale oil and gas reservoirs.

[0005] In a second aspect, an embodiment of the present invention further provides a device for generating elastic parameters of a shale reservoir containing organic nanopores, the device comprising: an acquisition module for acquiring target measurement parameters and physical morphology parameters of the shale reservoir containing organic nanopores; an organic module for generating organic elastic property parameters containing nanopores according to the physical morphology parameters of the organic matter; an inorganic module for generating inorganic elastic stiffness coefficients according to the target measurement parameters; a synthesis module for generating coupled elastic property parameters of a shale oil and gas reservoir according to the organic elastic property parameters and the inorganic elastic stiffness coefficients; and a calculation module for generating elastic parameters of a shale oil and gas reservoir according to the coupled elastic property parameters of the shale oil and gas reservoir.

[0006] In a third aspect, an embodiment of the present invention further provides a computer device, comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor implements the above-mentioned method for generating elastic parameters of organic nanoporous shale reservoirs when executing the computer program.

[0007] In a fourth aspect, an embodiment of the present invention further provides a computer-readable storage medium, wherein the computer-readable storage medium stores a computer program for executing the above-mentioned method for generating elastic parameters of an organic-containing nanoporous shale reservoir.

[0008] The embodiments of the present invention bring the following beneficial effects: the embodiments of the present invention provide a scheme for generating elastic parameters of shale reservoirs containing organic nanopores, which obtains target measurement parameters and physical morphology parameters of organic matter in shale reservoirs containing organic nanopores; generates organic elastic property parameters containing nanopores according to the physical morphology parameters of organic matter; generates inorganic elastic stiffness coefficients according to target measurement parameters; generates coupled elastic property parameters of shale oil and gas reservoirs according to organic elastic property parameters and inorganic elastic stiffness coefficients; generates elastic parameters of shale oil and gas reservoirs according to coupled elastic property parameters of shale oil and gas reservoirs. The embodiments of the present invention can generate elastic parameters that take into account the influence of organic nanopores, and the obtained elastic parameters can be used to improve the accuracy of prediction of porosity and related physical property parameters of shale oil and gas reservoirs.

[0009] Other features and advantages of the present invention will be described in the following description, and partly become apparent from the description, or understood by practicing the present invention. The purpose and other advantages of the present invention are realized and obtained by the structures particularly pointed out in the description, claims and drawings.

[0010] In order to make the above-mentioned objects, features and advantages of the present invention more obvious and easy to understand, preferred embodiments are given below and described in detail with reference to the accompanying drawings. BRIEF DESCRIPTION OF THE DRAWINGS

[0011] In order to more clearly illustrate the specific implementation methods of the present invention or the technical solutions in the prior art, the drawings required for use in the specific implementation methods or the description of the prior art will be briefly introduced below. Obviously, the drawings described below are some implementation methods of the present invention. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying creative work.

[0012] Figure 1 A flow chart of a method for generating elastic parameters of a shale reservoir containing organic nanopores provided in an embodiment of the present invention;

[0013] Figure 2 A workflow diagram for rock physics modeling of shale oil and gas reservoirs taking into account the influence of nanopores provided in an embodiment of the present invention;

[0014] Figure 3 A shale oil reservoir nanopore size-total porosity seismic rock physics template provided by an embodiment of the present invention;

[0015] Figure 4 A structural block diagram of a device for generating elastic parameters of organic nanoporous shale reservoirs provided in an embodiment of the present invention;

[0016] Figure 5 A structural block diagram of a computer device provided in an embodiment of the present invention. DETAILED DESCRIPTION

[0017] In order to make the purpose, technical solution and advantages of the embodiments of the present invention clearer, the technical solution of the present invention will be clearly and completely described below in conjunction with the accompanying drawings. Obviously, the described embodiments are part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of the present invention.

[0018] For unconventional shale oil reservoirs, the extensively developed nanopore system will have a very important impact on the elastic characteristics of the rock. The surface effect of the nanopore size can cause the impact of the longitudinal wave impedance and the longitudinal and shear wave velocity ratio exceeding 10%.

[0019] Based on this, an embodiment of the present invention provides a method and device for generating elastic parameters of a shale reservoir containing organic nanopores. This method targets the phenomenon of widespread nanoscale pores in shale oil and gas reservoirs, and establishes a seismic rock physics model of shale oil and gas reservoirs that takes into account the influence of organic nanopores, thereby greatly improving the accuracy of the existing predictions of porosity and related physical parameters of shale oil and gas reservoirs. This method introduces the microscopic mechanical equations of nanoporous media into the elastic properties of organic matter in shale oil and gas reservoirs, and then establishes a seismic rock physics model of shale oil and gas reservoirs that takes into account the influence of organic nanopores. This rock physics model can be used to guide quantitative seismic interpretation of shale oil and gas reservoirs based on logging and seismic data.

[0020] To facilitate understanding of this embodiment, a method for generating elastic parameters of an organic nanoporous shale reservoir disclosed in an embodiment of the present invention is first introduced in detail.

[0021] The embodiment of the present invention provides a method for generating elastic parameters of a shale reservoir containing organic nanopores, see Figure 1 A flow chart of a method for generating elastic parameters of an organic nanoporous shale reservoir is shown, the method comprising the following steps:

[0022] Step S102, obtaining target measurement parameters and physical morphology parameters of organic matter containing nanoporous shale reservoir.

[0023] In the embodiment of the present invention, for a specific shale oil and gas reservoir, the main physical parameters affecting the adsorption performance and elastic characteristics of the organic shale are extracted, mainly including target measurement parameters and organic matter physical morphology parameters.

[0024] Based on scanning electron microscopy and CT (Computed Tomography) analysis, the relevant parameters characterizing the physical form of organic matter are determined, including the aspect ratio of the organic phase, the radius of the organic pore, the aspect ratio of the organic pore and the inorganic pore, and the relative percentage of solid asphalt and kerogen. Based on XRD (Diffraction of X-Rays) mineral analysis, geochemical analysis, and porosity measurement, the content of inorganic minerals, organic carbon content, inorganic porosity, organic porosity, etc. are determined to obtain the target measurement parameters.

[0025] For example, see Table 1, which shows the petrophysical parameters required for shale oil reservoir modeling considering nanopore surface effects.

[0026]

[0027] Table 1

[0028] Step S104, generating elastic property parameters of organic materials containing nanopores according to physical morphology parameters of organic matter.

[0029] In an embodiment of the present invention, after obtaining the physical morphology parameters of organic matter, nanopore size information, organic porosity information, pore morphology information, etc. can be obtained, and this information can be used to calculate the organic elastic property parameters of shale oil and gas reservoirs taking into account the influence of organic nanopores.

[0030] Step S106, generating an inorganic elastic stiffness coefficient according to the target measurement parameter.

[0031] In the embodiment of the present invention, the composition of the inorganic phase of shale generally includes a clay portion containing bound water and a non-clay portion. The non-clay portion includes silicate minerals, carbonate minerals and inorganic pores containing free water.

[0032] The equivalent elastic properties of silicate minerals, carbonate minerals and inorganic pores containing free water are calculated according to the target measurement parameters to obtain the inorganic elastic stiffness coefficient.

[0033] Step S108, generating coupled shale oil and gas reservoir elastic property parameters according to the organic elastic property parameters and the inorganic elastic stiffness coefficients.

[0034] In the embodiment of the present invention, the organic elastic property parameters and the inorganic elastic stiffness coefficients are coupled to obtain the elastic property parameters of the shale oil and gas reservoir coupled with the organic phase and the inorganic phase.

[0035] Step S110, generating shale oil and gas reservoir elastic parameters according to the coupled shale oil and gas reservoir elastic property parameters.

[0036] In the embodiment of the present invention, the elastic characteristics of the shale oil and gas reservoir taking into account the coupling of adsorption performance and pore pressure, that is, the stiffness coefficient of the entire shale, are obtained. For example, it can be expressed in a matrix as follows:

[0037]

[0038] in, Then the density of the entire shale is obtained by taking the weighted average of the volume fractions of each phase, combined with the following formula: and In the formula, ρ represents density, V s (0) and V p (0) represents velocity, and the longitudinal wave velocity and transverse wave velocity of the shale oil and gas reservoir can be further calculated to obtain the elastic parameters of the shale oil and gas reservoir.

[0039] It should be noted that the elastic parameters of shale oil and gas reservoirs can be expressed as the ratio of longitudinal and transverse wave velocities or longitudinal wave impedance, which can be set according to actual needs and are not specifically limited in the embodiments of the present invention.

[0040] The embodiment of the present invention provides a scheme for generating elastic parameters of shale reservoirs containing organic nanopores, which obtains target measurement parameters and physical morphology parameters of organic matter in shale reservoirs containing organic nanopores; generates organic elastic property parameters containing nanopores according to the physical morphology parameters of organic matter; generates inorganic elastic stiffness coefficients according to target measurement parameters; generates coupled elastic property parameters of shale oil and gas reservoirs according to the organic elastic property parameters and the inorganic elastic stiffness coefficients; and generates elastic parameters of shale oil and gas reservoirs according to the coupled elastic property parameters of shale oil and gas reservoirs. The embodiment of the present invention can generate elastic parameters that take into account the influence of organic nanopores, and the obtained elastic parameters can be used to improve the accuracy of prediction of porosity and related physical property parameters of shale oil and gas reservoirs.

[0041] Considering the surface energy effect caused by the higher specific surface area of ​​nano-scale pores, the micromechanical equations of nanoporous media are introduced into the characterization of the elastic properties of organic matter in shale oil and gas reservoirs. The elastic property parameters of organic matter containing nanopores are generated according to the physical morphology parameters of organic matter, which can be performed as follows:

[0042] The bulk modulus of the organic phase containing nanopores or the shear modulus of the organic phase containing nanopores is generated according to the physical morphology parameters of the organic matter.

[0043] In an embodiment of the present invention, the bulk modulus and shear modulus of the organic phase portion of the shale containing solid kerogen and organic nanopore fluid are calculated based on the physical morphology parameters of the organic matter.

[0044] In one embodiment, the following formula can be used to generate the bulk modulus of the nanopore-containing organic phase according to the physical morphology parameters of the organic matter:

[0045]

[0046]

[0047] Among them, K o is the bulk modulus of the organic phase containing nanopores, K s is the solid kerogen bulk modulus, μ s is the solid kerogen shear modulus, K f is the bulk modulus of the fluid filling the organic nanopores, f is the volume fraction of the organic nanopores in the organic phase, r is the radius of the organic nanopores, K surface It is the interfacial bulk modulus that characterizes the interfacial effect of the nanopore surface, and its unit is N / m.

[0048] It should be noted that, in the embodiment of the present invention, if the surface energy effect is not considered, the above equation will degenerate into the classic Hashin-Shtrikman equation in elastic mechanics.

[0049] In one embodiment, the shear modulus of the nanopore-containing organic phase can be generated according to the physical morphology parameters of the organic matter using the following formula:

[0050]

[0051]

[0052]

[0053]

[0054]

[0055]

[0056]

[0057]

[0058] Among them, μ o is the shear modulus of the organic phase containing nanopores, μ s is the solid kerogen shear modulus, ν s is the Poisson's ratio of solid kerogen, μ fis the shear modulus of the fluid filling the organic nanopores, f is the volume fraction of the organic nanopores in the organic phase, ν f is the Poisson's ratio of the fluid filling the organic nanopore, r is the radius of the organic nanopore, μ surface is the interfacial shear modulus that characterizes the interfacial effect of the nanopore surface, K surface It is the interfacial bulk modulus that characterizes the interfacial effect of the nanopore surface, and its unit is N / m.

[0059] It should be noted that, in the embodiment of the present invention, if the surface energy effect is not considered, the above equation will degenerate into the classic Hashin-Shtrikman equation in elastic mechanics.

[0060] In one embodiment, generating an inorganic elastic stiffness coefficient according to a target measurement parameter may be performed according to the following steps:

[0061] The comprehensive elastic modulus of silicate minerals and carbonate minerals is calculated according to the target measurement parameters; the equivalent modulus of the non-clay part is determined according to the inorganic pore parameter and the comprehensive elastic modulus; the inorganic elastic stiffness coefficient is calculated according to the equivalent modulus of the non-clay part and the target measurement parameters.

[0062] In the embodiment of the present invention, the non-clay part includes silicate minerals (quartz and feldspar), carbonate minerals (calcite and dolomite) and inorganic pores. Then, the comprehensive results of silicate minerals and carbonate minerals are calculated using Hill average, and then the inorganic phase is constructed by integrating clay and non-clay. The anisotropic SCA model is used to integrate the clay part and the non-clay part to obtain the equivalent elastic modulus of the inorganic phase part of the shale oil and gas reservoir, which is used as the inorganic elastic stiffness coefficient to complete the construction of the equivalent elastic properties of the inorganic phase part.

[0063] In one embodiment, the comprehensive elastic modulus of silicate minerals and carbonate minerals is calculated according to the target measurement parameters according to the following formula: Among them, M VRH Represents the comprehensive elastic modulus of silicate minerals and carbonate minerals, f i represents the volume fraction of the i-th phase, M i represents the elastic modulus of the ith phase.

[0064] In the embodiment of the present invention, M i It can be the bulk modulus or the shear modulus. According to the above formula, when N is 2, the calculation can be performed according to the following steps:

[0065] Elastic modulus of silicate minerals: M1 = f1 × M feldspar +(1-f1)×M quartz ;

[0066] Elastic modulus of carbonate minerals: M2 = f2 × M calcite +(1-f2)×M dolomite ;

[0067] in,

[0068] The comprehensive results obtained according to Hill average are as follows:

[0069] M V =(f feldspar +f quartz )M1+(f calcite +f dolomite )M2

[0070]

[0071] Among them, f feldspare 、f quartz 、f calcite 、f dolomite Represent the volume fractions of feldspar, quartz, calcite and dolomite respectively; M feldspare 、M quartz 、M calcite 、M dolomite represents the corresponding elastic modulus, and M can be either the bulk modulus or the shear modulus.

[0072] In one embodiment, the equivalent modulus of the non-clay part is determined according to the inorganic pore parameter and the comprehensive elastic modulus according to the following formula: Where i represents the i-th material, x i represents the volume content of the ith material, represents the equivalent bulk modulus of the non-clay part, It represents the equivalent shear modulus of the non-clay part. P and Q are geometric factors related to the aspect ratio of each phase. The superscript *i of P and Q means that this geometric factor is for the self-compatible equivalent modulus. and The inclusion material i in the background medium, K i represents the bulk modulus of the non-clay part, μ i Represents the shear modulus of the non-clay part.

[0073] In the embodiment of the present invention, after the comprehensive elastic modulus of silicate minerals and carbonate minerals is calculated, the inorganic pores are added using the SCA model to obtain the equivalent elastic modulus of the non-clay part.

[0074] In one embodiment, the inorganic elastic stiffness coefficient is calculated according to the equivalent modulus of the non-clay part and the target measurement parameter according to the following formula: Among them, Csca represents the inorganic elastic stiffness coefficient, I is the unit tensor, v n Indicates the volume percentage of clay and non-clay parts, C n is the stiffness matrix of the n-th phase material, which is converted from the equivalent modulus of the non-clay part, and the tensor G is the Eshelby tensor.

[0075] In an embodiment of the present invention, the value of N can be 2. When n is 1, it can correspond to the calculation of parameters related to the clay part. When n is 2, it can correspond to the calculation of parameters related to the non-clay part. The above formula is used to comprehensively calculate the clay part and the non-clay part to obtain the inorganic elastic stiffness coefficient.

[0076] The formula of the anisotropic SCA model is as follows:

[0077]

[0078] The clay part is given a modulus setting based on the VTI anisotropy assumption (Sayers et al., 2013): c 11 =44.9GPa,c 33 =24.2Gpa,c 44 =3.7Gpa; c 66 =11.6Gpa,c 13 =18.1GPa c 12 =c 11 -2×c 66 ;

[0079]

[0080] The equivalent elastic modulus of the non-clay part is or Its formula is as follows:

[0081] in

[0082]

[0083]

[0084]

[0085] Here, v n (n=1,2……N) represents the volume percentage of each component, C n The stiffness coefficients representing the various components are given by the equivalent modulus of the non-clay part and Convert to C SCA Represents the hybrid stiffness coefficient of the self-consistent overall inorganic phase part, which is a second-order tensor. It is a fourth-order tensor that uses Green's function to characterize the geometry of inclusion compounds.

[0086] In one embodiment, generating coupled shale oil and gas reservoir elastic property parameters according to organic elastic property parameters and inorganic elastic stiffness coefficients may be performed according to the following steps:

[0087] The elastic flexibility coefficient is generated according to the organic elastic property parameters and the inorganic elastic stiffness coefficient; and the coupled shale oil and gas reservoir elastic property parameters are generated according to the elastic flexibility coefficient.

[0088] In the embodiment of the present invention, the anisotropic SCA model is used to fill the kerogen containing organic pores into the inorganic skeleton to obtain a solid background replaced by a solid. The specific calculation can also refer to the formula: The aspect ratio of the inclusions is consistent with that of porous kerogen. Since solid asphalt does not play a skeletal supporting role in rocks, the anisotropic SCA model is used to fill the empty pores with an elastic modulus of 0 into the solid replacement background to obtain the elastic stiffness coefficient of organic shale without asphalt. The aspect ratio of the empty pores is consistent with that of organic pore asphalt, and the volume of the empty pores is the same as that of organic pore asphalt. Anisotropic Gassmann solid replacement is used to characterize the effect of filling with organic pore asphalt on the elastic parameters of the overall shale oil and gas reservoir. The solid replacement formula of anisotropic media can be expressed by the flexibility coefficient.

[0089] In one embodiment, the coupled shale oil and gas reservoir elastic property parameters are generated according to the elastic compliance coefficient according to the following formula: in, represents the elastic compliance coefficient, Represents the coupled shale oil and gas reservoir elastic property parameters, represents the effective flexibility tensor of the mineral, β fl represents the compressibility of pore fluid, β0 represents the compressibility of mineral, and φ represents the porosity.

[0090] In the embodiment of the present invention, the solid replacement formula of the anisotropic medium can be expressed by the flexibility coefficient. The flexibility tensor is obtained by the inverse matrix of the elastic stiffness coefficient. The fourth-order tensor of the stiffness coefficient is required. The conversion relationship is shown in Table 2 below:

[0091] ij(kl) I(J) 11 1 22 2 33 3 23,32 4 13,31 5 12,21 6

[0092] Table 2

[0093]

[0094]

[0095] are the effective elastic stiffness coefficients of dry rock (in the above formula, C SCA The fourth-order tensor of the mineral effective elastic stiffness coefficient (in the above formula ), They are the effective flexibility tensor of rock with fillings, the effective flexibility tensor of dry rock, the effective flexibility tensor of minerals, and the flexibility tensor is the inverse matrix of the stiffness coefficient; K fl is the elastic modulus of the fluid, K0 is the bulk modulus of the solid background, β fl , β0 are the compressibility of pore fluid and mineral, respectively K fl and the reciprocal of K0. This is the final elastic property of shale oil and gas that reflects the inorganic and organic coupling parts. The inverse conversion of the flexibility matrix into a stiffness matrix can obtain the longitudinal and transverse wave velocities of the shale oil and gas reservoir.

[0096] See also Figure 2 The flowchart of rock physics modeling of shale oil and gas reservoirs considering the influence of nanopores is shown in FIG. 1 . The embodiment of the present invention provides a method and device for generating elastic parameters of shale reservoirs containing organic nanopores, see Figure 3 The shale oil reservoir nanopore size-total porosity seismic rock physics template shown in the figure shows the shale oil reservoir nanopore size-porosity seismic rock physics template considering the nanopore surface interface effect. It can be seen that the nanopores under small size (1nm) and large size conditions, for shale oil reservoirs with a porosity of 10%, have an influence of more than 10% on the longitudinal wave impedance and longitudinal and transverse wave velocity ratio. As the nanopore size increases, the longitudinal wave impedance and aspect ratio velocity ratio show a downward trend. As the nanopore size increases, the organic pore surface energy also decreases accordingly. It can be seen that this method can construct different seismic rock physics templates to characterize the influence of different reservoir physical parameters on the elastic characteristics of shale oil and gas reservoirs under different nanopore size conditions.

[0097] The present invention also provides an organic nanoporous shale reservoir elastic parameter generation device, as described in the following embodiments. Since the principle of the device to solve the problem is similar to the method for generating elastic parameters of organic nanoporous shale reservoirs, the implementation of the device can refer to the implementation of the method for generating elastic parameters of organic nanoporous shale reservoirs, and the repeated parts will not be repeated. Figure 4 The structural block diagram of the device for generating elastic parameters of organic nanoporous shale reservoirs shown in FIG. 1 includes:

[0098] An acquisition module 71 is used to acquire target measurement parameters and physical morphology parameters of organic matter containing nanopores in a shale reservoir; an organic module 72 is used to generate organic elastic property parameters containing nanopores according to the physical morphology parameters of organic matter; an inorganic module 73 is used to generate inorganic elastic stiffness coefficients according to target measurement parameters; a synthesis module 74 is used to generate coupled elastic property parameters of a shale oil and gas reservoir according to the organic elastic property parameters and the inorganic elastic stiffness coefficients; and a calculation module 75 is used to generate elastic parameters of a shale oil and gas reservoir according to the coupled elastic property parameters of a shale oil and gas reservoir.

[0099] In one embodiment, the organic module is specifically used to generate a bulk modulus of an organic phase containing nanopores or a shear modulus of an organic phase containing nanopores according to physical morphology parameters of organic matter.

[0100] In one embodiment, the organic module is specifically used to generate the bulk modulus of the organic phase containing nanopores according to the physical morphology parameters of the organic matter using the following formula:

[0101]

[0102]

[0103] Among them, K o is the bulk modulus of the organic phase containing nanopores, K s is the solid kerogen bulk modulus, μ s is the solid kerogen shear modulus, K f is the bulk modulus of the fluid filling the organic nanopores, f is the volume fraction of the organic nanopores in the organic phase, r is the radius of the organic nanopores, K surface It is the interfacial bulk modulus that characterizes the interfacial effect of the nanopore surface, and its unit is N / m.

[0104] In one embodiment, the organic module is specifically used to generate the shear modulus of the nanopore-containing organic phase according to the physical morphology parameters of the organic matter using the following formula:

[0105]

[0106]

[0107]

[0108]

[0109]

[0110]

[0111]

[0112]

[0113] Among them, μ o is the shear modulus of the organic phase containing nanopores, μ s is the solid kerogen shear modulus, ν s is the Poisson's ratio of solid kerogen, μ f is the shear modulus of the fluid filling the organic nanopores, f is the volume fraction of the organic nanopores in the organic phase, ν f is the Poisson's ratio of the fluid filling the organic nanopore, r is the radius of the organic nanopore, μ surface is the interfacial shear modulus that characterizes the interfacial effect of the nanopore surface, K surface It is the interfacial bulk modulus that characterizes the interfacial effect of the nanopore surface, and its unit is N / m.

[0114] In one embodiment, the inorganic module is specifically used to: calculate the comprehensive elastic modulus of silicate minerals and carbonate minerals according to the target measurement parameters; determine the equivalent modulus of the non-clay part according to the inorganic pore parameters and the comprehensive elastic modulus; calculate the inorganic elastic stiffness coefficient according to the equivalent modulus of the non-clay part and the target measurement parameters.

[0115] In one embodiment, the inorganic module is specifically used to calculate the comprehensive elastic modulus of silicate minerals and carbonate minerals according to the target measurement parameters according to the following formula: Among them, M VRH Represents the comprehensive elastic modulus of silicate minerals and carbonate minerals, f i represents the volume fraction of the i-th phase, M i represents the elastic modulus of the i-th phase.

[0116] In one embodiment, the inorganic module is specifically used to determine the equivalent modulus of the non-clay part according to the inorganic pore parameter and the comprehensive elastic modulus according to the following formula: Where i represents the i-th material, x i represents the volume content of the ith material, represents the equivalent bulk modulus of the non-clay part, It represents the equivalent shear modulus of the non-clay part. P and Q are geometric factors related to the aspect ratio of each phase. The superscript *i of P and Q means that this geometric factor is for the self-compatible equivalent modulus. and The inclusion material i in the background medium, K i represents the bulk modulus of the non-clay part, μ i represents the shear modulus of the non-clay part

[0117] In one embodiment, the inorganic module is specifically used to calculate the inorganic elastic stiffness coefficient according to the equivalent modulus of the non-clay part and the target measurement parameter according to the following formula: Among them, C sca represents the inorganic elastic stiffness coefficient, I is the unit tensor, v n Indicates the volume percentage of clay and non-clay parts, C n is the stiffness matrix of the n-th phase material, which is converted from the equivalent modulus of the non-clay part, and the tensor G is the Eshelby tensor.

[0118] In one embodiment, the synthesis module is specifically used to: generate an elastic compliance coefficient according to organic elastic property parameters and inorganic elastic stiffness coefficients; and generate coupled shale oil and gas reservoir elastic property parameters according to the elastic compliance coefficients.

[0119] In one embodiment, the synthesis module is specifically used to generate coupled shale oil and gas reservoir elastic property parameters according to the elastic compliance coefficient according to the following formula: in, represents the elastic compliance coefficient, Represents the coupled shale oil and gas reservoir elastic property parameters, represents the effective flexibility tensor of the mineral, β fl represents the compressibility of pore fluid, β0 represents the compressibility of mineral, and φ represents the porosity.

[0120] The embodiment of the present invention also provides a computer device, see Figure 5 The computer device structure schematic block diagram shown includes a memory 81, a processor 82, and a computer program stored in the memory and executable on the processor. When the processor executes the computer program, the steps of any one of the above-mentioned methods for generating elastic parameters of organic nanoporous shale reservoirs are implemented.

[0121] Those skilled in the art can clearly understand that, for the convenience and brevity of description, the specific working process of the computer device described above can refer to the corresponding process in the aforementioned method embodiment, and will not be repeated here.

[0122] An embodiment of the present invention further provides a computer-readable storage medium storing a computer program for executing any of the above-mentioned methods for generating elastic parameters of an organic-containing nanoporous shale reservoir.

[0123] Those skilled in the art will appreciate that embodiments of the present invention may be provided as methods, systems, or computer program products. Therefore, the present invention may take the form of a complete hardware embodiment, a complete software embodiment, or an embodiment combining software and hardware. Moreover, the present invention may take the form of a computer program product implemented on one or more computer-usable storage media (including but not limited to disk storage, CD-ROM, optical storage, etc.) containing computer-usable program code.

[0124] The present invention is described with reference to flowcharts and / or block diagrams of methods, devices (systems), and computer program products according to embodiments of the present invention. It should be understood that each process and / or block in the flowchart and / or block diagram, as well as the combination of processes and / or blocks in the flowchart and / or block diagram, can be implemented by computer program instructions. These computer program instructions can be provided to a processor of a general-purpose computer, a special-purpose computer, an embedded processor, or other programmable data processing device to produce a machine, so that the instructions executed by the processor of the computer or other programmable data processing device generate instructions for implementing the processes in the flowchart and / or block diagram. Figure 1 A process or multiple processes and / or boxes Figure 1 A device that provides the functions specified in a block or multiple blocks.

[0125] These computer program instructions may also be stored in a computer-readable memory capable of directing a computer or other programmable data processing device to operate in a specific manner, so that the instructions stored in the computer-readable memory produce an article of manufacture comprising an instruction device, which implements the process Figure 1 A process or multiple processes and / or boxes Figure 1 A function specified in one or more boxes.

[0126] These computer program instructions can also be loaded onto a computer or other programmable data processing device so that a series of operating steps are executed on the computer or other programmable device to produce a computer-implemented process, thereby providing instructions for implementing the process. Figure 1 A process or multiple processes and / or boxes Figure 1 The steps for the functions specified in one or more boxes.

[0127] Finally, it should be noted that the above-described embodiments are only specific implementations of the present invention, which are used to illustrate the technical solutions of the present invention, rather than to limit them. The protection scope of the present invention is not limited thereto. Although the present invention is described in detail with reference to the above-described embodiments, ordinary technicians in the field should understand that any technician familiar with the technical field can still modify the technical solutions recorded in the above-described embodiments within the technical scope disclosed by the present invention, or can easily think of changes, or make equivalent replacements for some of the technical features therein; and these modifications, changes or replacements do not make the essence of the corresponding technical solutions deviate from the spirit and scope of the technical solutions of the embodiments of the present invention, and should be included in the protection scope of the present invention. Therefore, the protection scope of the present invention shall be subject to the protection scope of the claims.

Claims

1. A method for generating elastic parameters of organic nanoporous shale reservoirs, characterized in that: include: Obtaining target measurement parameters and physical morphology parameters of the organic-containing nanoporous shale reservoir; Generating elastic property parameters of organic matter containing nanopores according to the physical morphology parameters of the organic matter; generating an inorganic elastic stiffness coefficient according to the target measurement parameter; generating coupled shale oil and gas reservoir elastic property parameters according to the organic elastic property parameters and the inorganic elastic stiffness coefficient; Generating shale oil and gas reservoir elasticity parameters according to the coupled shale oil and gas reservoir elasticity property parameters; Generating elastic property parameters of organic matter containing nanopores according to the physical morphology parameters of the organic matter includes: Generating a bulk modulus of an organic phase containing nanopores or a shear modulus of an organic phase containing nanopores according to the physical morphology parameters of the organic matter; The bulk modulus of the nanopore-containing organic phase is generated according to the physical morphology parameters of the organic matter using the following formula: Among them, K o is the bulk modulus of the organic phase containing nanopores, K s is the solid kerogen bulk modulus, μ s is the solid kerogen shear modulus, K f is the bulk modulus of the fluid filling the organic nanopores, f is the volume fraction of the organic nanopores in the organic phase, r is the radius of the organic nanopores, K surface The interfacial bulk modulus characterizes the interfacial effect of the nanopore surface, in N / m; The shear modulus of the nanopore-containing organic phase is generated according to the physical morphology parameters of the organic matter using the following formula: Among them, μ o is the shear modulus of the organic phase containing nanopores, μ s is the solid kerogen shear modulus, ν s is the Poisson's ratio of solid kerogen, μ f is the shear modulus of the fluid filling the organic nanopores, f is the volume fraction of the organic nanopores in the organic phase, ν f is the Poisson's ratio of the fluid filling the organic nanopore, r is the radius of the organic nanopore, μ surface is the interfacial shear modulus that characterizes the interfacial effect of the nanopore surface, K surface It is the interfacial bulk modulus that characterizes the interfacial effect of the nanopore surface, and its unit is N / m.

2. The method according to claim 1, characterized in that Generating an inorganic elastic stiffness coefficient according to the target measurement parameter includes: Calculating the comprehensive elastic modulus of silicate minerals and carbonate minerals according to the target measurement parameters; Determining the equivalent modulus of the non-clay part according to the inorganic pore parameter and the comprehensive elastic modulus; The inorganic elastic stiffness coefficient is calculated according to the equivalent modulus of the non-clay part and the target measurement parameter.

3. The method according to claim 2, characterized in that The method includes calculating the comprehensive elastic modulus of silicate minerals and carbonate minerals according to the target measurement parameters according to the following formula: Among them, M VRH Represents the comprehensive elastic modulus of silicate minerals and carbonate minerals, f i represents the volume fraction of the i-th phase, M i represents the elastic modulus of the ith phase.

4. The method according to claim 2, characterized in that: The method includes determining the equivalent modulus of the non-clay part according to the inorganic pore parameter and the comprehensive elastic modulus according to the following formula: Where i represents the i-th material, x i represents the volume content of the ith material, represents the equivalent bulk modulus of the non-clay part, It represents the equivalent shear modulus of the non-clay part. P and Q are geometric factors related to the aspect ratio of each phase. The superscript *i of P and Q means that this geometric factor is for the self-compatible equivalent modulus. and The inclusion material i in the background medium, K i represents the bulk modulus of the non-clay part, μ i Represents the shear modulus of the non-clay part.

5. The method according to claim 2, characterized in that: The method includes calculating the inorganic elastic stiffness coefficient according to the equivalent modulus of the non-clay part and the target measurement parameter according to the following formula: Among them, C sca represents the inorganic elastic stiffness coefficient, I is the unit tensor, v n Indicates the volume percentage of clay and non-clay parts, C n is the stiffness matrix of the n-phase material, which is converted from the equivalent modulus of the non-clay part, and the tensor G is the Eshelby tensor; is a fourth-order tensor that uses Green's function to characterize the geometry of the inclusion compound; Represent the fourth-order tensor The index arrangements are in the form of ijkl, ikjl, and jkil, which are used to describe the tensor characteristics under different index arrangements.

6. The method according to claim 1, characterized in that Generate coupled shale oil and gas reservoir elastic property parameters according to the organic elastic property parameters and the inorganic elastic stiffness coefficient, including: generating an elastic compliance coefficient according to the organic elastic property parameter and the inorganic elastic stiffness coefficient; The coupled shale oil and gas reservoir elastic property parameters are generated according to the elastic flexibility coefficient.

7. The method according to claim 6, characterized in that The method includes generating coupled shale oil and gas reservoir elastic property parameters according to the elastic flexibility coefficient according to the following formula: in, and They represent the index arrangement of the elastic flexibility tensor in the form of ijkl, ijaa, klaa and ααββ respectively; Represents the coupled shale oil and gas reservoir elastic property parameters, as well as They represent the index arrangement of the mineral effective flexibility tensor as ijaa, klaa and ααββ, β fl represents the compressibility of pore fluid, β0 represents the compressibility of mineral, and φ represents the porosity.

8. A device for generating elastic parameters of organic nanoporous shale reservoirs, characterized in that: include: An acquisition module, used to acquire target measurement parameters and physical form parameters of organic matter in the organic nanoporous shale reservoir; An organic module, used for generating elastic property parameters of organic materials containing nanopores according to the physical morphology parameters of the organic matter; An inorganic module, used for generating an inorganic elastic stiffness coefficient according to the target measurement parameter; A synthesis module, used for generating coupled shale oil and gas reservoir elastic property parameters according to the organic elastic property parameters and the inorganic elastic stiffness coefficient; A calculation module, used for generating shale oil and gas reservoir elastic parameters according to the coupled shale oil and gas reservoir elastic property parameters; The organic module is specifically used for: Generating a bulk modulus of an organic phase containing nanopores or a shear modulus of an organic phase containing nanopores according to the physical morphology parameters of the organic matter; The organic module is specifically used to generate the bulk modulus of the nanopore-containing organic phase according to the physical morphology parameters of the organic matter using the following formula: Among them, K o is the bulk modulus of the organic phase containing nanopores, K s is the solid kerogen bulk modulus, μ s is the solid kerogen shear modulus, K f is the bulk modulus of the fluid filling the organic nanopores, f is the volume fraction of the organic nanopores in the organic phase, r is the radius of the organic nanopores, K surface The interfacial bulk modulus characterizes the interfacial effect of the nanopore surface, in N / m; The organic module is specifically used to generate the shear modulus of the nanopore-containing organic phase according to the physical morphology parameters of the organic matter using the following formula: Among them, μ o is the shear modulus of the organic phase containing nanopores, μ s is the solid kerogen shear modulus, ν s is the Poisson's ratio of solid kerogen, μ f is the shear modulus of the fluid filling the organic nanopores, f is the volume fraction of the organic nanopores in the organic phase, ν f is the Poisson's ratio of the fluid filling the organic nanopore, r is the radius of the organic nanopore, μ surface is the interfacial shear modulus that characterizes the interfacial effect of the nanopore surface, K surface It is the interfacial bulk modulus that characterizes the interfacial effect of the nanopore surface, and its unit is N / m.

9. The device according to claim 8, characterized in that The inorganic module is specifically used for: Calculating the comprehensive elastic modulus of silicate minerals and carbonate minerals according to the target measurement parameters; Determining the equivalent modulus of the non-clay part according to the inorganic pore parameter and the comprehensive elastic modulus; The inorganic elastic stiffness coefficient is calculated according to the equivalent modulus of the non-clay part and the target measurement parameter.

10. The device according to claim 9, characterized in that The inorganic module is specifically used to calculate the comprehensive elastic modulus of silicate minerals and carbonate minerals according to the target measurement parameters according to the following formula: Among them, M VRH Represents the comprehensive elastic modulus of silicate minerals and carbonate minerals, f i represents the volume fraction of the i-th phase, M i represents the elastic modulus of the ith phase.

11. The device according to claim 9, characterized in that The inorganic module is specifically used to determine the equivalent modulus of the non-clay part according to the inorganic pore parameter and the comprehensive elastic modulus according to the following formula: Where i represents the i-th material, x i represents the volume content of the ith material, represents the equivalent bulk modulus of the non-clay part, It represents the equivalent shear modulus of the non-clay part. P and Q are geometric factors related to the aspect ratio of each phase. The superscript *i of P and Q means that this geometric factor is for the self-compatible equivalent modulus. and The inclusion material i in the background medium, K i represents the bulk modulus of the non-clay part, μ i Represents the shear modulus of the non-clay part.

12. The device according to claim 9, characterized in that The inorganic module is specifically used to calculate the inorganic elastic stiffness coefficient according to the equivalent modulus of the non-clay part and the target measurement parameter according to the following formula: Among them, C sca represents the inorganic elastic stiffness coefficient, I is the unit tensor, v n Indicates the volume percentage of clay and non-clay parts, C n is the stiffness matrix of the n-phase material, which is converted from the equivalent modulus of the non-clay part, and the tensor G is the Eshelby tensor; is a fourth-order tensor that uses Green's function to characterize the geometry of the inclusion compound; Represent the fourth-order tensor The index arrangements are in the form of ijkl, ikjl, and jkil, which are used to describe the tensor characteristics under different index arrangements.

13. The device according to claim 8, characterized in that The synthesis module is specifically used for: generating an elastic compliance coefficient according to the organic elastic property parameter and the inorganic elastic stiffness coefficient; The coupled shale oil and gas reservoir elastic property parameters are generated according to the elastic flexibility coefficient.

14. The device according to claim 13, characterized in that The synthesis module is specifically used to generate coupled shale oil and gas reservoir elastic property parameters according to the elastic flexibility coefficient according to the following formula: in, and They represent the index arrangement of the elastic flexibility tensor in the form of ijkl, ijaa, klaa and ααββ respectively; Represents the coupled shale oil and gas reservoir elastic property parameters, as well as They represent the index arrangement of the mineral effective flexibility tensor as ijaa, klaa and ααββ, β fl represents the compressibility of pore fluid, β0 represents the compressibility of mineral, and φ represents the porosity.

15. A computer device comprising a memory, a processor and a computer program stored in the memory and executable on the processor, characterized in that: When the processor executes the computer program, the method for generating elastic parameters of an organic-containing nanoporous shale reservoir as claimed in any one of claims 1 to 7 is implemented.

16. A computer-readable storage medium, characterized in that: The computer-readable storage medium stores a computer program, and when the computer program is executed by a processor, the method for generating elastic parameters of an organic-containing nanoporous shale reservoir according to any one of claims 1 to 7 is implemented.

Images