Rock physical chart establishing method and device, electronic equipment and storage medium

Abstract

The embodiment of the invention discloses a rock physical plate establishing method and device, electronic equipment and a storage medium. The method comprises the following steps: for a shale gas reservoir to be subjected to rock physical chart establishment, determining the elastic modulus of a rock matrix corresponding to the shale gas reservoir; for the dry pores to be added into the rock matrix, obtaining the elastic modulus of the first rock skeleton according to the elastic modulus of the dry pores and the elastic modulus of the rock matrix; for a directional crack to be added into the first rock skeleton, obtaining an elastic stiffness matrix of a second rock skeleton according to the directional weakness of the directional crack and the elastic modulus of the first rock skeleton; according to the elastic stiffness matrix of the second rock skeleton, a rock physical model is established, and according to the rock physical model, a rock physical chart of the shale gas reservoir is established. According to the technical scheme provided by the embodiment of the invention, the anisotropic rock physical plate can be established, so that accurate characterization of the shale gas reservoir is realized.

Description

Technical Field

[0001] The embodiments of the present invention relate to the field of geophysical exploration technology, and in particular to a method, device, electronic device and storage medium for establishing a rock physics plate. Background Art

[0002] Shale oil and gas reservoirs have become an important area for increasing onshore oil and gas reserves, and tight shale gas reservoirs are gradually becoming a new research hotspot for seismic reservoir prediction. On this basis, various scholars have proved that rock physics maps can be used to predict the physical parameters of shale gas reservoirs intuitively, efficiently and conveniently.

[0003] However, the rock physics maps currently established cannot accurately characterize shale gas reservoirs and need to be improved. Summary of the invention

[0004] The embodiments of the present invention provide a method, device, electronic device and storage medium for establishing a rock physics plate, so as to establish an anisotropic rock physics plate, thereby achieving accurate characterization of shale gas reservoirs.

[0005] According to one aspect of the present invention, a method for establishing a rock physics plate is provided, which may include:

[0006] For the shale gas reservoir to be subjected to rock physics plate establishment, determine the elastic modulus of the rock matrix corresponding to the shale gas reservoir;

[0007] For the dry pores to be added to the rock matrix, the elastic modulus of the first rock skeleton is obtained according to the elastic modulus of the dry pores and the elastic modulus of the rock matrix, wherein the first rock skeleton is obtained by adding the dry pores to the rock matrix;

[0008] For the directional cracks to be added to the first rock skeleton, an elastic stiffness matrix of the second rock skeleton is obtained according to the directional weakness of the directional cracks and the elastic modulus of the first rock skeleton, wherein the second rock skeleton is obtained by adding the directional cracks to the first rock skeleton;

[0009] A rock physics model is established based on the elastic stiffness matrix of the second rock skeleton, and a rock physics map of the shale gas reservoir is established based on the rock physics model.

[0010] According to another aspect of the present invention, a rock physics plate creation device is provided, which may include:

[0011] A rock matrix elastic modulus obtaining module is used to determine the elastic modulus of the rock matrix corresponding to the shale gas reservoir for which the rock physics map is to be established;

[0012] A rock skeleton elastic modulus obtaining module is used to obtain the elastic modulus of a first rock skeleton according to the elastic modulus of the dry pores and the elastic modulus of the rock matrix for the dry pores to be added to the rock matrix, wherein the first rock skeleton is obtained by adding the dry pores to the rock matrix;

[0013] An elastic stiffness matrix obtaining module is used to obtain an elastic stiffness matrix of a second rock skeleton according to the directional weakness of the directional cracks to be added to the first rock skeleton and the elastic modulus of the first rock skeleton, wherein the second rock skeleton is obtained by adding the directional cracks to the first rock skeleton;

[0014] The rock physics plate establishing module is used to establish a rock physics model according to the elastic stiffness matrix of the second rock skeleton, and to establish a rock physics plate of the shale gas reservoir according to the rock physics model.

[0015] According to another aspect of the present invention, there is provided an electronic device, which may include:

[0016] at least one processor; and

[0017] a memory communicatively connected to at least one processor; wherein,

[0018] The memory stores a computer program that can be executed by at least one processor. The computer program is executed by at least one processor so that the at least one processor implements the rock physical plate establishment method provided by any embodiment of the present invention when executing the computer program.

[0019] According to another aspect of the present invention, a computer-readable storage medium is provided, on which computer instructions are stored, and the computer instructions are used to implement the rock physical plate creation method provided by any embodiment of the present invention when executed by a processor.

[0020] The technical solution of the embodiment of the present invention is to determine the elastic modulus of the rock matrix corresponding to the shale gas reservoir for which a rock physics map is to be established; for the dry pores to be added to the rock matrix, the elastic modulus of the first rock skeleton containing the dry pores is obtained according to the elastic modulus of the dry pores and the elastic modulus of the rock matrix; for the directional cracks to be added to the first rock skeleton, the elastic stiffness matrix of the second rock skeleton containing the dry pores and the directional cracks is obtained according to the directional weakness of the directional cracks and the elastic modulus of the first rock skeleton; a rock physics model is established according to the elastic stiffness matrix of the second rock skeleton, and a rock physics map of the shale gas reservoir is established according to the rock physics model. The above technical scheme, through the anisotropic equivalent medium theory and the poroelasticity theory, establishes a rock physics model of shale gas reservoirs that can reasonably characterize the fracture-pore type (i.e., porous structure). Then, based on the rock physics model, a rock physics plate can be established that can accurately characterize shale gas reservoirs with strong heterogeneity and developed anisotropic fractures. The rock physics plate can be used to accurately characterize the shale gas reservoirs.

[0021] It should be understood that the contents described in this section are not intended to identify the key or important features of the embodiments of the present invention, nor are they intended to limit the scope of the present invention. Other features of the present invention will become easily understood through the following description. BRIEF DESCRIPTION OF THE DRAWINGS

[0022] In order to more clearly illustrate the technical solutions in the embodiments of the present invention, the following briefly introduces the drawings required for use in the description of the embodiments. Obviously, the drawings described below are only some embodiments of the present invention. For ordinary technicians in this field, other drawings can be obtained based on these drawings without creative work.

[0023] Figure 1 is a flow chart of a method for establishing a rock physics plate according to an embodiment of the present invention;

[0024] Figure 2 is a flow chart of another method for establishing a rock physics plate provided according to an embodiment of the present invention;

[0025] Figure 3 is a flow chart of another method for establishing a rock physics plate provided according to an embodiment of the present invention;

[0026] Figure 4 It is a flowchart of an example of inversion of physical property parameters of shale gas reservoirs based on anisotropic three-dimensional rock physics plates in another rock physics plate establishment method provided in an embodiment of the present invention;

[0027] Figure 5aIt is a schematic diagram of the comparison result between anisotropic three-dimensional rock physical plate and logging data (TOC content) in another rock physical plate establishment method provided according to an embodiment of the present invention;

[0028] Figure 5b is a schematic diagram of a comparison result between an anisotropic three-dimensional rock physics plate and well logging data (POR) in another rock physics plate establishment method provided according to an embodiment of the present invention;

[0029] Figure 5c is a schematic diagram of the comparison result between anisotropic three-dimensional rock physical plate and well logging data (VSH) in another rock physical plate establishment method provided according to an embodiment of the present invention;

[0030] Figure 6 It is a schematic diagram of TOC content, porosity and shale content of a target layer section of a shale gas well based on anisotropic three-dimensional rock physics plate inversion in another rock physics plate establishment method provided in an embodiment of the present invention;

[0031] Figure 7 It is a schematic diagram of longitudinal wave impedance, shear wave impedance and density obtained by AVO inversion based on HTI medium reflection coefficient equation in another rock physics plate establishment method provided in an embodiment of the present invention;

[0032] Figure 8 It is a schematic diagram of comparison results between TOC content, porosity and shale content obtained by inversion of anisotropic three-dimensional rock physical plate and logging interpretation in another rock physical plate establishment method provided in an embodiment of the present invention;

[0033] Fig. 9 is a structural block diagram of a rock physics plate establishing device provided according to an embodiment of the present invention;

[0034] Fig.10 It is a schematic diagram of the structure of an electronic device for implementing the rock physics plate establishment method in an embodiment of the present invention. DETAILED DESCRIPTION

[0035] In order to enable those skilled in the art to better understand the scheme of the present invention, the technical scheme in the embodiments of the present invention will be clearly and completely described below in conjunction with the drawings in the embodiments of the present invention. Obviously, the described embodiments are only 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 should fall within the scope of protection of the present invention.

[0036] It should be noted that the terms "first", "second", etc. in the specification and claims of the present invention and the above-mentioned drawings are used to distinguish similar objects, and are not necessarily used to describe a specific order or sequence. It should be understood that the data used in this way can be interchangeable where appropriate, so that the embodiments of the present invention described herein can be implemented in an order other than those illustrated or described herein. The situations of "target", "original", etc. are similar and will not be repeated here. In addition, the terms "including" and "having" and any of their variations are intended to cover non-exclusive inclusions, for example, a process, method, system, product or device that includes a series of steps or units is not necessarily limited to those steps or units that are clearly listed, but may include other steps or units that are not clearly listed or inherent to these processes, methods, products or devices.

[0037] Figure 1 1 is a flow chart of a method for establishing a rock physics plate provided in an embodiment of the present invention. This embodiment is applicable to the case of establishing anisotropic rock physics plates. The method can be executed by a rock physics plate establishing device provided in an embodiment of the present invention, which can be implemented by software and / or hardware, and the device can be integrated in an electronic device, which can be various user terminals or servers.

[0038] See also Figure 1 The method of the embodiment of the present invention specifically includes the following steps:

[0039] S110. For a shale gas reservoir for which a rock physics map is to be established, determine the elastic modulus of a rock matrix corresponding to the shale gas reservoir.

[0040] The shale gas reservoir can be understood as a shale gas reservoir for which a rock physics map is to be established. The rock matrix can be understood as a rock matrix corresponding to the shale gas reservoir. In practical applications, the rock matrix can also be referred to as a solid mineral matrix. The elastic modulus of the rock matrix can be, for example, the bulk modulus and / or shear modulus of the rock matrix.

[0041] S120. For the dry pores to be added into the rock matrix, the elastic modulus of the first rock skeleton is obtained according to the elastic modulus of the dry pores and the elastic modulus of the rock matrix, wherein the first rock skeleton is obtained by adding the dry pores into the rock matrix.

[0042] Among them, dry pores can be understood as dry (i.e., containing no pore fluid) pores to be added to the rock matrix. In practical applications, optionally, the dry pores can be randomly distributed dry ellipsoidal pores with a single pore aspect ratio; and optionally, the dry pores can be added to the rock matrix through a DEM model. The first rock skeleton can be understood as the rock skeleton obtained after adding the dry pores to the rock matrix, that is, the first rock skeleton is a rock skeleton containing dry pores.

[0043] The elastic modulus of the dry pores is obtained, and the elastic modulus of the first rock skeleton is obtained according to the elastic modulus of the dry pores and the elastic modulus of the rock matrix. In combination with the above example, the elastic modulus of the first rock skeleton can be obtained by combining the elastic modulus of the dry pores and the elastic modulus of the rock matrix through a DEM model.

[0044] S130. For the directional cracks to be added to the first rock skeleton, obtain the elastic stiffness matrix of the second rock skeleton according to the directional weakness of the directional cracks and the elastic modulus of the first rock skeleton, wherein the second rock skeleton is obtained by adding the directional cracks to the first rock skeleton.

[0045] The directional fractures may be understood as the directional fractures to be added to the first rock skeleton. The directional weakness may be understood as the weakness of the directional fracture in a certain direction, for example, normal weakness and / or tangential weakness.

[0046] The second rock skeleton can be understood as a rock skeleton obtained after adding directional fractures to the first rock skeleton, that is, the second rock skeleton is a rock skeleton containing dry pores and directional fractures. Here, by adding directional fractures to the first rock skeleton, the anisotropy induced by directional fractures can be considered in the rock physics modeling process. In practical applications, optionally, directional fractures can be added to the shale gas reservoir, specifically to the first rock skeleton, through a linear sliding model.

[0047] On this basis, the elastic stiffness matrix of the second rock skeleton can be obtained according to the directional weakness and the elastic modulus of the first rock skeleton. In practical applications, optionally, the Lame constant can be obtained according to the elastic modulus of the first rock skeleton, and then the elastic stiffness matrix of the second rock skeleton can be obtained according to the Lame constant and the directional weakness.

[0048] S140. A rock physics model is established based on the elastic stiffness matrix of the second rock skeleton, and a rock physics map of the shale gas reservoir is established based on the rock physics model.

[0049] Among them, since the elastic stiffness matrix of the second rock skeleton is obtained on the basis of considering directional fractures and dry pores, that is, it is an elastic stiffness matrix that matches the porous structure, and the shale gas reservoir is a porous structure, therefore, according to the elastic stiffness matrix of the second rock skeleton, a rock physics model that can accurately characterize the shale gas reservoir can be established, that is, an anisotropic rock physics model of the shale gas reservoir can be established. On this basis, further, according to the rock physics model, a rock physics plate of the shale gas reservoir is established, and the rock physics plate is an anisotropic three-dimensional rock physics plate.

[0050] The technical solution of the embodiment of the present invention is to determine the elastic modulus of the rock matrix corresponding to the shale gas reservoir for which a rock physics map is to be established; for the dry pores to be added to the rock matrix, the elastic modulus of the first rock skeleton containing the dry pores is obtained according to the elastic modulus of the dry pores and the elastic modulus of the rock matrix; for the directional cracks to be added to the first rock skeleton, the elastic stiffness matrix of the second rock skeleton containing the dry pores and the directional cracks is obtained according to the directional weakness of the directional cracks and the elastic modulus of the first rock skeleton; a rock physics model is established according to the elastic stiffness matrix of the second rock skeleton, and a rock physics map of the shale gas reservoir is established according to the rock physics model. The above technical scheme, through the anisotropic equivalent medium theory and the poroelasticity theory, establishes a rock physics model of shale gas reservoirs that can reasonably characterize the fracture-pore type (i.e., porous structure). Then, based on the rock physics model, a rock physics plate can be established that can accurately characterize shale gas reservoirs with strong heterogeneity and developed anisotropic fractures. The rock physics plate can be used to accurately characterize the shale gas reservoirs.

[0051] On this basis, an optional technical solution is to obtain the elastic modulus of the first rock skeleton according to the elastic modulus of the dry pores and the elastic modulus of the rock matrix, including: determining the porosity of the dry pores; solving the elastic modulus of the first rock skeleton according to the elastic modulus and porosity of the dry pores and taking the elastic modulus of the rock matrix as the initial value.

[0052] The porosity can be understood as the porosity of dry pores. In practical applications, the porosity can be determined based on the total porosity and fracture volume corresponding to the shale gas reservoir. The elastic modulus of the rock matrix is ​​used as the initial value for the solution, and the elastic modulus of the dry pores and the porosity of the dry pores are combined to solve the elastic modulus of the first rock skeleton.

[0053] The above technical solution realizes the accurate solution of the elastic modulus of the first rock skeleton.

[0054] Another optional technical solution is to establish a rock physics model based on the elastic stiffness matrix of the second rock skeleton, which may include: obtaining the bulk modulus of the pore fluid in the dry pores and oriented fractures based on the bulk modulus of water, the bulk modulus of gas and the water saturation corresponding to the shale gas reservoir; obtaining the elastic stiffness matrix of saturated rock based on the elastic modulus of the rock matrix, the elastic stiffness matrix of the second rock skeleton and the bulk modulus of the pore fluid, wherein the dry pores and oriented fractures in the saturated rock are filled with saturated pore fluid; and establishing a rock physics model based on the elastic stiffness matrix of the saturated rock.

[0055] Among them, pore fluid can be understood as the fluid in dry pores and oriented fractures. The bulk modulus of pore fluid is obtained according to the bulk modulus of water, bulk modulus of gas and water saturation corresponding to the shale gas reservoir. In practical applications, the bulk modulus of pore fluid can be calculated using the Wood formula for the non-uniform gas and water distribution characteristics of shale gas reservoirs.

[0056] Furthermore, saturated rock can be understood as a rock in which dry pores and directional fractures are filled with saturated pore fluid. The elastic stiffness matrix of saturated rock is obtained according to the bulk modulus, the elastic modulus of the rock matrix and the elastic stiffness matrix of the second rock skeleton. For example, the elastic stiffness matrix of saturated rock can be calculated by the anisotropic Gassmann equation.

[0057] The above technical solution fully considers the influence of water saturation on elastic parameters during the rock physics modeling process, which can greatly improve the subsequent prediction accuracy of physical property parameters.

[0058] Another optional technical solution, determining the elastic modulus of the rock matrix corresponding to the shale gas reservoir, may include: determining the elastic modulus of the brittle mineral mixture corresponding to the shale gas reservoir; determining the elastic modulus of the organic-rich clay minerals corresponding to the shale gas reservoir, wherein the organic-rich clay minerals are obtained by adding organic matter to the clay minerals; for the rock matrix composed of the brittle mineral mixture and the organic-rich clay minerals, obtaining the elastic modulus of the rock matrix according to the elastic modulus of the brittle mineral mixture and the elastic modulus of the organic-rich clay minerals.

[0059] The brittle mineral mixture can be understood as a mixture of minerals such as quartz, feldspar and calcite in a shale gas reservoir, specifically, a brittle mineral mixture without clay. The elastic modulus of the brittle mineral mixture is determined. In practical applications, the elastic modulus can be calculated by a Voigt-Reuss-Hill (VRH) boundary average model.

[0060] Organic-rich clay minerals can be obtained by adding organic matter to clay minerals. The elastic modulus of organic-rich clay minerals is determined. In practical applications, optionally, the clay mineral can also be called clay mineral; further optionally, the organic matter can be kerogen, and on this basis, optionally, when the kerogen is evenly distributed in the clay mineral and is spherical in shape, the KT model can be used to add kerogen to the clay mineral to obtain organic-rich clay minerals.

[0061] For a rock matrix composed of a brittle mineral mixture and organic-rich clay minerals, the elastic modulus of the rock matrix can be obtained by using the elastic modulus of the brittle mineral mixture and the elastic modulus of the organic-rich clay minerals. In practical applications, the elastic modulus of the rock matrix can be calculated by the SCA model.

[0062] The above technical solution realizes the accurate calculation of the elastic modulus of the rock matrix.

[0063] Figure 2 It is a flow chart of another method for establishing a rock physical plate provided in an embodiment of the present invention. This embodiment is optimized based on the above-mentioned technical solutions. In this embodiment, optionally, according to the rock physical model, a rock physical plate of a shale gas reservoir is established, including: according to the logging data corresponding to the shale gas reservoir, a range of variation of physical property parameters is obtained, and grid division is performed within the range of variation of the physical property parameters to obtain multiple grid nodes: for each grid node in the multiple grid nodes, according to the rock physical model and the physical property parameters represented by the grid node, the elastic parameters represented by the grid node are obtained; according to the elastic parameters respectively represented by the multiple grid nodes, a rock physical plate of the shale gas reservoir is constructed. Among them, the explanations of the terms that are the same or corresponding to the above-mentioned embodiments are not repeated here.

[0064] See also Figure 2 The method of this embodiment may specifically include the following steps:

[0065] S210: For the shale gas reservoir for which a rock physics map is to be established, determine the elastic modulus of the rock matrix corresponding to the shale gas reservoir.

[0066] S220. For the dry pores to be added to the rock matrix, the elastic modulus of the first rock skeleton is obtained according to the elastic modulus of the dry pores and the elastic modulus of the rock matrix, wherein the first rock skeleton is obtained by adding the dry pores to the rock matrix.

[0067] S230. For the directional cracks to be added to the first rock skeleton, an elastic stiffness matrix of the second rock skeleton is obtained according to the directional weakness of the directional cracks and the elastic modulus of the first rock skeleton, wherein the second rock skeleton is obtained by adding the directional cracks to the first rock skeleton.

[0068] S240. A rock physics model is established based on the elastic stiffness matrix of the second rock skeleton.

[0069] S250, obtaining a variation range of physical property parameters according to the well logging data corresponding to the shale gas reservoir, and performing grid division within the variation range of the physical property parameters to obtain a plurality of grid nodes.

[0070] Among them, according to the well logging data corresponding to the shale gas reservoir, physical property parameters such as TOC content, porosity and shale content can be obtained, and then according to these physical property parameters, the range of physical property parameter variation (i.e., physical property parameter domain) can be obtained. Further, grid division is performed within the range of physical property parameter variation to obtain multiple grid nodes. It can be understood that each of the multiple grid nodes corresponds to its own physical property parameter.

[0071] S260. For each grid node among the multiple grid nodes, obtain an elastic parameter represented by the grid node according to the rock physics model and the physical property parameters represented by the grid node.

[0072] Among them, the rock physics model established through the above steps can accurately reflect the mapping relationship or quantitative relationship between the elastic parameters and physical parameters of the shale gas reservoir. Therefore, for each grid node in the multiple grid nodes, the elastic parameters represented by the grid node can be calculated according to the physical parameters represented by the grid node and the rock physics model. In practical applications, the elastic parameters can be, optionally, longitudinal and transverse wave impedance and / or density.

[0073] S270. Establish a rock physics map of the shale gas reservoir based on the elastic parameters respectively represented by the plurality of grid nodes.

[0074] A rock physics map is established according to the elastic parameters represented by each of the multiple grid nodes. For example, the rock physics map can be generated by connecting contour lines, that is, connecting elastic parameters with the same values.

[0075] The technical solution of the embodiment of the present invention obtains the range of variation of physical property parameters through logging data, and then performs grid division within the range of variation of physical property parameters to obtain multiple grid nodes; further, the physical property parameters respectively represented by each grid node in the multiple grid nodes and the rock physics model can be used to calculate the elastic parameters respectively represented by these grid nodes, and then a rock physics plate can be established according to these elastic parameters, thereby realizing the effective establishment of the rock physics plate.

[0076] On this basis, an optional technical solution is that after establishing a rock physical map of the shale gas reservoir according to the elastic parameters respectively represented by multiple grid nodes, the above-mentioned rock physical map establishment method also includes: projecting the logging data into the rock physical map to adjust the elastic parameters respectively represented by the multiple grid nodes based on the obtained projection results; and correcting the rock physical map according to the adjusted elastic parameters respectively represented by the multiple grid nodes.

[0077] The above technical solution adjusts the elastic parameters respectively represented by multiple grid nodes by projecting logging data, and then corrects the rock physics map according to the adjusted elastic parameters, thereby further improving the accuracy of shale gas reservoir characterization.

[0078] Figure 3 It is a flow chart of another method for establishing a rock physics plate provided in an embodiment of the present invention. This embodiment is optimized based on the above-mentioned technical solutions. In this embodiment, optionally, after establishing a rock physics plate of a shale gas reservoir according to the elastic parameters respectively represented by a plurality of grid nodes, the above-mentioned method for establishing a rock physics plate may also include: obtaining target data, wherein the target data is logging data, or elastic parameters obtained by inverting seismic data corresponding to the shale gas reservoir; projecting the target data into the rock physics plate to determine the target node matching the target data from each grid node in the rock physics plate; and using the physical property parameters represented by the target node as the physical property parameters corresponding to the target data. Among them, the explanations of the terms that are the same or corresponding to the above-mentioned embodiments are not repeated here.

[0079] See also Figure 3 The method of this embodiment may specifically include the following steps:

[0080] S3010. For a shale gas reservoir for which a rock physics map is to be established, determine the elastic modulus of the rock matrix corresponding to the shale gas reservoir.

[0081] S3020. For the dry pores to be added to the rock matrix, the elastic modulus of the first rock skeleton is obtained according to the elastic modulus of the dry pores and the elastic modulus of the rock matrix, wherein the first rock skeleton is obtained by adding the dry pores to the rock matrix.

[0082] S3030. For the directional cracks to be added to the first rock skeleton, obtain the elastic stiffness matrix of the second rock skeleton according to the directional weakness of the directional cracks and the elastic modulus of the first rock skeleton, wherein the second rock skeleton is obtained by adding the directional cracks to the first rock skeleton.

[0083] S3040. A rock physics model is established based on the elastic stiffness matrix of the second rock skeleton.

[0084] S3050. Obtain a range of physical property parameters according to the well logging data corresponding to the shale gas reservoir, and perform grid division within the range of physical property parameters to obtain a plurality of grid nodes.

[0085] S3060. For each grid node among the multiple grid nodes, obtain the elastic parameters represented by the grid node according to the rock physics model and the physical property parameters represented by the grid node.

[0086] S3070. Establish a rock physics map of the shale gas reservoir based on the elastic parameters respectively represented by the multiple grid nodes.

[0087] S3080, obtaining target data, wherein the target data is well logging data, or elastic parameters obtained by inverting seismic data corresponding to the shale gas reservoir.

[0088] S3090, projecting the target data into the rock physics plate to determine the target node matching the target data from each grid node in the rock physics plate.

[0089] The target data is projected into the rock physics map, so that a target node matching the target data can be found in the rock physics map. For example, an optimization algorithm can be used to find the target node closest to the target data in the rock physics map.

[0090] S3100. Use the physical property parameters represented by the target node as the physical property parameters corresponding to the target data.

[0091] The technical solution of the embodiment of the present invention uses a grid search method to carry out the inversion process of the physical property parameters of the shale gas reservoir using rock physics plates, thereby achieving a high-precision and stable quantitative characterization of the physical property parameters.

[0092] In order to better understand the above-mentioned technical solutions as a whole, the following is an exemplary description with reference to specific examples. Figure 4 It can be mainly divided into three parts: establishment of rock physics model of fracture-pore type shale gas reservoir, establishment of anisotropic three-dimensional rock physics map and inversion of anisotropic three-dimensional rock physics map.

[0093] 1) Establishment of rock physics model of fracture-pore shale gas reservoir

[0094] VRH boundary average model, KT model, SCA model, DEM model, linear sliding model, Wood formula and anisotropic Gassmann equation are comprehensively used to gradually couple the mineral components, pore fluid and pore structure to establish a fracture-pore shale gas reservoir rock physics model. The specific implementation process is as follows: Steps (1) to (6):

[0095] (1) Using the VRH boundary average model, the elastic modulus (i.e., bulk modulus or shear modulus) of a brittle mineral mixture composed of quartz, feldspar, and calcite is calculated as:

[0096]

[0097] Among them, f i represents the volume content of the i-th component, M i represents the bulk modulus or shear modulus of the i-th component, M VRH represents the elastic modulus of a brittle mineral mixture without clay, and M VRH The bulk modulus is specifically denoted as K nck And the shear modulus is specifically denoted as μ nck .

[0098] (2) Assuming that kerogen is uniformly distributed in clay minerals and is spherical in shape, the KT model is used to add organic matter (i.e., kerogen) to clay to obtain the bulk modulus K of organic-rich clay minerals: ck and the shear modulus μ ck :

[0099]

[0100]

[0101]

[0102] Where: K c is the bulk modulus of clay, K k is the bulk modulus of kerogen, μ c is the shear modulus of clay, μ k is the shear modulus of kerogen, v k is the volume fraction of kerogen, Q ck and P ck are all geometric factors.

[0103] (3) Using the SCA model, the bulk modulus K of the rock matrix based on a mixture of brittle minerals and organic-rich clay minerals is calculated. ma and shear modulus μ ma :

[0104] f ck (K ck -K ma ) ck +f nck (K nck -K ma ) nck =0 (5)

[0105] fck (μ ck -μ ma )Q ck +f nck (μ nck -μ ma )Q nck =0 (6)

[0106] Among them, f ck is the volume content of organic-rich clay minerals, f nck is the volume content of the brittle mineral mixture, satisfying f nck +f ck =1,P nck , Q nck , P ck and Q ck are all geometric factors.

[0107] (4) Using the DEM model, randomly distributed ellipsoidal dry pores with a single pore aspect ratio are added to the rock matrix to obtain the bulk modulus of the first rock skeleton containing dry pores (i.e., the dry pore rock skeleton): and shear modulus

[0108]

[0109]

[0110] Among them, K i is the bulk modulus of the dry pore, μ i is the shear modulus of the dry pore, which can be taken as 0; φ p =φ(1-ξ f ), φ p is the porosity of dry pores, φ is the total porosity, ξ f is the crack volume; P *i and Q *i are polarization factors. It is worth noting that the DEM model is very effective in solving and When K ma and μ ma As the initial value, that is

[0111] On this basis, in order to consider the anisotropy induced by directional fractures, directional fractures can be added to the shale gas reservoir through the linear sliding model based on the dry pore rock skeleton calculated above, and the elastic stiffness matrix C of the second rock skeleton containing dry pores and directional fractures can be obtained: dry :

[0112]

[0113] Where λ and μ are both Lame constants, Δ N is the normal weakness of the directional crack, Δ T is the tangential weakness of a directional crack.

[0114] (5) Based on the fluid distribution characteristics of shale gas reservoirs, the Wood formula is used to calculate the bulk modulus K of the pore fluid: f :

[0115]

[0116] Among them, S w is water saturation, K w is the bulk modulus of water, K g is the bulk modulus of gas.

[0117] (6) After determining the bulk modulus K of the rock matrix ma , the elastic stiffness matrix C of the second rock skeleton dry (That is, in the following text ) and the bulk modulus of the pore fluid K f Afterwards, the elastic stiffness matrix (i.e. bulk modulus and shear modulus) of the saturated rock can be calculated using the anisotropic Gassmann equation:

[0118]

[0119] Among them, for m=1,2,3, where a 4 =a 5 =a 6 =0;

[0120]

[0121]

[0122]

[0123]

[0124] Where L is the fluid spatial modulus. So far, the rock physics model has been established.

[0125] (7) According to The longitudinal wave velocity V can be calculated p and shear wave velocity V s :

[0126] N=[(C 11 -C 55 )sin 2θ-(C 33 -C 55 )cos 2 θ] 2 +4(C 31 +C 55 ) 2 sin 2 θcos 2 θ (16)

[0128]

[0129]

[0130] Among them, θ is the angle between the wavefront direction and the Z axis, and ρ is the density of the rock matrix.

[0131] 2) Establishment of anisotropic three-dimensional rock physics map

[0132] Physical parameters such as TOC content, porosity and shale content are obtained by combining well logging data, and then the range of physical parameter variation can be obtained based on these physical parameters. Further, three-dimensional grid division is performed within the range of physical parameter variation, and then the above-established fracture-pore shale gas reservoir rock physics model is used to calculate the elastic parameters of each grid node (for example, longitudinal and transverse wave impedance and density), and the contour lines are connected to generate anisotropic three-dimensional rock physics plates. Further, the well logging data is projected onto the rock physics plate, and the elastic parameters are dynamically adjusted to achieve rock physics plate correction.

[0133] 3) Anisotropic 3D rock physics map inversion

[0134] The deterministic inversion of anisotropic three-dimensional rock physics map is to project the logging data or the elastic parameters of prestack inversion onto the rock physics map, and use the optimization algorithm to find the grid node closest to the logging data or elastic parameters. The TOC content, porosity and shale content corresponding to the grid node are the predicted physical property parameters.

[0135] On this basis, in order to prove the effectiveness of the rock physics map established above, relevant experiments were carried out and the specific experimental results are given below.

[0136] See also Figure 5a-5c, they are the anisotropic three-dimensional rock physics plates of logarithmic P-wave impedance-logarithmic S-wave impedance-density for TOC content, porosity (POR) and shale content (VSH) based on the rock physics model of fracture-pore shale gas reservoirs, and the logging data of the target layer of a shale gas well in southwest China are projected onto the rock physics plates. By comparing the logging data and the rock physics plates, it is found that the projected data points all fall inside the rock physics plates, and the scatter distribution characteristics are consistent with the change law of the rock physics plates.

[0137] See also Figure 6 , which shows the TOC content, porosity and shale content predicted from the P- and S-wave impedance and density logging data by combining anisotropic 3D rock physics charts and grid search method. The dotted line represents the logging interpretation result, and the solid line represents the inversion result, i.e. the prediction result. The application results of logging data show that the predicted TOC content, porosity and shale content are generally consistent with the logging interpretation results.

[0138] See also Figure 7 , which shows the elastic parameter profile of a shale gas well obtained by anisotropic AVO inversion of the longitudinal isotropic (Horizontal Transverse Isotropy, HTI) medium. As can be seen from the figure, the target layer has low longitudinal wave impedance (10-11.5km / s*g / cm3), low density (2.4-2.55g / cm3) and relatively low shear wave impedance (6.5-7.5km / s*g / cm3). Overall, the elastic parameters of anisotropic AVO inversion can effectively characterize the distribution of underground elastic parameters and provide reliable data support for the quantitative prediction of physical parameters of shale gas reservoirs.

[0139] See also Figure 8 It can be seen that the TOC content varies between 0 and 4%, which is basically consistent with the logging interpretation results; the porosity distribution range is between 0 and 8%, and the predicted value is slightly higher than the logging interpretation results, but the overall change trend is consistent with the logging interpretation; the shale content varies between 0 and 60%, which is highly consistent with the logging interpretation results. Overall, it can be seen that the seismic prediction of physical parameters is consistent with the logging interpretation results, which verifies the effectiveness of the anisotropic three-dimensional rock physics plate inversion method.

[0140] Fig. 9 This is a structural block diagram of a rock physics plate establishment device provided in an embodiment of the present invention, and the device is used to execute the rock physics plate establishment method provided in any of the above embodiments. The device and the rock physics plate establishment method of the above embodiments belong to the same inventive concept. For details not described in detail in the embodiments of the rock physics plate establishment device, please refer to the embodiments of the above rock physics plate establishment method. Fig. 9 The device specifically includes: a rock matrix elastic modulus obtaining module 410, a rock skeleton elastic modulus obtaining module 420, an elastic stiffness matrix obtaining module 430 and a rock physics plate establishing module 440. Among them,

[0141] The rock matrix elastic modulus obtaining module 410 can be used to determine the elastic modulus of the rock matrix corresponding to the shale gas reservoir for which the rock physics map is to be established;

[0142] The rock skeleton elastic modulus obtaining module 420 may be used to obtain the elastic modulus of a first rock skeleton according to the elastic modulus of the dry pores and the elastic modulus of the rock matrix for the dry pores to be added to the rock matrix, wherein the first rock skeleton is obtained by adding the dry pores to the rock matrix;

[0143] The elastic stiffness matrix obtaining module 430 can be used to obtain the elastic stiffness matrix of the second rock skeleton according to the directional weakness of the directional cracks to be added to the first rock skeleton and the elastic modulus of the first rock skeleton, wherein the second rock skeleton is obtained by adding the directional cracks to the first rock skeleton;

[0144] The rock physics plate establishing module 440 can be used to establish a rock physics model according to the elastic stiffness matrix of the second rock skeleton, and to establish a rock physics plate of the shale gas reservoir according to the rock physics model.

[0145] Optionally, the rock physics plate building module 440 may include:

[0146] The grid nodes are used to obtain the variation range of physical parameters according to the well logging data corresponding to the shale gas reservoir, and to perform grid division within the variation range of the physical parameters to obtain multiple grid nodes:

[0147] An elastic parameter obtaining unit is used to obtain, for each grid node among the plurality of grid nodes, an elastic parameter represented by the grid node according to a rock physics model and a physical property parameter represented by the grid node;

[0148] The rock physics map establishing unit is used to establish a rock physics map of the shale gas reservoir according to the elastic parameters respectively represented by a plurality of grid nodes.

[0149] On this basis, optionally, the rock physics plate establishing device further includes:

[0150] The elastic parameter adjustment module can be used to project the logging data into the rock physics map after constructing the rock physics map of the shale gas reservoir according to the elastic parameters respectively represented by the multiple grid nodes, so as to adjust the elastic parameters respectively represented by the multiple grid nodes based on the obtained projection results;

[0151] The rock physics plate correction module can be used to correct the rock physics plate according to the elastic parameters respectively represented by the adjusted multiple grid nodes.

[0152] Alternatively, the rock physics plate establishing device may further include:

[0153] The target data acquisition module can be used to acquire target data after constructing a rock physical map of the shale gas reservoir according to the elastic parameters respectively represented by the multiple grid nodes, wherein the target data is well logging data, or elastic parameters obtained by inverting the seismic data corresponding to the shale gas reservoir;

[0154] A target node determination module may be used to project the target data into a rock physics plate to determine a target node matching the target data from each grid node in the rock physics plate;

[0155] The physical property parameter obtaining module can be used to use the physical property parameters represented by the target node as the physical property parameters corresponding to the target data.

[0156] Optionally, the rock skeleton elastic modulus obtaining module 420 includes:

[0157] a porosity determination unit for determining the porosity of the dried pores;

[0158] The rock skeleton elastic modulus solving unit is used to solve the elastic modulus of the first rock skeleton according to the elastic modulus and porosity of the dry pores and the elastic modulus of the rock matrix as an initial value.

[0159] Optionally, the rock physics plate building module 440 may include:

[0160] A bulk modulus obtaining unit is used to obtain the bulk modulus of pore fluids in dry pores and oriented fractures according to the bulk modulus of water, the bulk modulus of gas and the water saturation corresponding to the shale gas reservoir;

[0161] A rock stiffness matrix obtaining unit is used to obtain an elastic stiffness matrix of a saturated rock according to an elastic modulus of a rock matrix, an elastic stiffness matrix of a second rock skeleton, and a bulk modulus of a pore fluid, wherein dry pores and oriented fractures in the saturated rock are filled with saturated pore fluid;

[0162] The rock physics model building unit is used to build a rock physics model according to the elastic stiffness matrix of saturated rock.

[0163] Optionally, the rock matrix elastic modulus obtaining module 410 may include:

[0164] A brittle mineral mixture elastic modulus determination unit, used to determine the elastic modulus of the brittle mineral mixture corresponding to the shale gas reservoir;

[0165] An organic-rich clay mineral elastic modulus determination unit, used to determine the elastic modulus of the organic-rich clay mineral corresponding to the shale gas reservoir, wherein the organic-rich clay mineral is obtained by adding organic matter to the clay mineral;

[0166] The rock matrix elastic modulus obtaining unit is used for obtaining the elastic modulus of the rock matrix based on the brittle mineral mixture and the organic-rich clay mineral composition according to the elastic modulus of the brittle mineral mixture and the elastic modulus of the organic-rich clay mineral.

[0167] The rock physics plate establishing device provided by the embodiment of the present invention determines the elastic modulus of the rock matrix corresponding to the shale gas reservoir for which the rock physics plate is to be established through a rock matrix elastic modulus obtaining module; obtains the elastic modulus of the first rock skeleton containing the dry pores according to the elastic modulus of the dry pores and the elastic modulus of the rock matrix for the dry pores to be added to the rock matrix through a rock skeleton elastic modulus obtaining module; obtains the elastic stiffness matrix of the second rock skeleton containing both the dry pores and the directional fractures according to the directional weakness of the directional fractures and the elastic modulus of the first rock skeleton for the directional fractures to be added to the first rock skeleton through an elastic stiffness matrix obtaining module; establishes a rock physics model according to the elastic stiffness matrix of the second rock skeleton through a rock physics plate establishing module, and establishes a rock physics plate of the shale gas reservoir according to the rock physics model. The above-mentioned device, through the anisotropic equivalent medium theory and the poroelasticity theory, establishes a rock physics model of shale gas reservoirs that can reasonably characterize the fracture-pore type (i.e., porous structure). Then, based on the rock physics model, a rock physics plate can be established that can accurately characterize shale gas reservoirs with strong heterogeneity and developed anisotropic fractures. The rock physics plate can be used to accurately characterize the shale gas reservoirs.

[0168] The rock physics plate establishing device provided in the embodiment of the present invention can execute the rock physics plate establishing method provided in any embodiment of the present invention, and has the corresponding functional modules and beneficial effects of the execution method.

[0169] It is worth noting that in the embodiment of the above-mentioned rock physical plate establishment device, the various units and modules included are only divided according to functional logic, but are not limited to the above-mentioned division, as long as the corresponding functions can be achieved; in addition, the specific names of the functional units are only for the convenience of distinguishing each other, and are not used to limit the scope of protection of the present invention.

[0170] Fig.10 A schematic diagram of the structure of an electronic device 10 that can be used to implement an embodiment of the present invention is shown. The electronic device is intended to represent various forms of digital computers, such as laptop computers, desktop computers, workstations, personal digital assistants, servers, blade servers, mainframe computers, and other suitable computers. The electronic device can also represent various forms of mobile devices, such as personal digital processing, cellular phones, smart phones, wearable devices (such as helmets, glasses, watches, etc.) and other similar computing devices. The components shown herein, their connections and relationships, and their functions are merely examples and are not intended to limit the implementation of the present invention described and / or required herein.

[0171] like Fig.10 As shown, the electronic device 10 includes at least one processor 11, and a memory connected to the at least one processor 11, such as a read-only memory (ROM) 12, a random access memory (RAM) 13, etc., wherein the memory stores a computer program that can be executed by at least one processor, and the processor 11 can perform various appropriate actions and processes according to the computer program stored in the read-only memory (ROM) 12 or the computer program loaded from the storage unit 18 to the random access memory (RAM) 13. In the RAM 13, various programs and data required for the operation of the electronic device 10 can also be stored. The processor 11, the ROM 12, and the RAM 13 are connected to each other through a bus 14. An input / output (I / O) interface 15 is also connected to the bus 14.

[0172] A number of components in the electronic device 10 are connected to the I / O interface 15, including: an input unit 16, such as a keyboard, a mouse, etc.; an output unit 17, such as various types of displays, speakers, etc.; a storage unit 18, such as a disk, an optical disk, etc.; and a communication unit 19, such as a network card, a modem, a wireless communication transceiver, etc. The communication unit 19 allows the electronic device 10 to exchange information / data with other devices through a computer network such as the Internet and / or various telecommunication networks.

[0173] The processor 11 may be a variety of general and / or special processing components with processing and computing capabilities. Some examples of the processor 11 include, but are not limited to, a central processing unit (CPU), a graphics processing unit (GPU), various special artificial intelligence (AI) computing chips, various processors running machine learning model algorithms, a digital signal processor (DSP), and any appropriate processor, controller, microcontroller, etc. The processor 11 executes the various methods and processes described above, such as a rock physics plate establishment method.

[0174] In some embodiments, the rock physical plate establishment method can be implemented as a computer program, which is tangibly contained in a computer-readable storage medium, such as a storage unit 18. In some embodiments, part or all of the computer program can be loaded and / or installed on the electronic device 10 via the ROM 12 and / or the communication unit 19. When the computer program is loaded into the RAM 13 and executed by the processor 11, one or more steps of the rock physical plate establishment method described above can be performed. Alternatively, in other embodiments, the processor 11 can be configured to execute the rock physical plate establishment method in any other appropriate manner (for example, by means of firmware).

[0175] Various implementations of the systems and techniques described above herein can be implemented in digital electronic circuit systems, integrated circuit systems, field programmable gate arrays (FPGAs), application specific integrated circuits (ASICs), application specific standard products (ASSPs), systems on chips (SOCs), load programmable logic devices (CPLDs), computer hardware, firmware, software, and / or combinations thereof. These various implementations can include: being implemented in one or more computer programs that can be executed and / or interpreted on a programmable system including at least one programmable processor, which can be a special purpose or general purpose programmable processor that can receive data and instructions from a storage system, at least one input device, and at least one output device, and transmit data and instructions to the storage system, the at least one input device, and the at least one output device.

[0176] Computer programs for implementing the methods of the present invention may be written in any combination of one or more programming languages. These computer programs may be provided to a processor of a general-purpose computer, a special-purpose computer, or other programmable data processing device, so that when the computer program is executed by the processor, the functions / operations specified in the flow chart and / or block diagram are implemented. The computer program may be executed entirely on the machine, partially on the machine, partially on the machine as a stand-alone software package and partially on a remote machine, or entirely on a remote machine or server.

[0177] In the context of the present invention, a computer-readable storage medium may be a tangible medium that may contain or store a computer program for use by or in combination with an instruction execution system, device or equipment. A computer-readable storage medium may include, but is not limited to, electronic, magnetic, optical, electromagnetic, infrared, or semiconductor systems, devices or equipment, or any suitable combination of the foregoing. Alternatively, a computer-readable storage medium may be a machine-readable signal medium. A more specific example of a machine-readable storage medium may include an electrical connection based on one or more lines, a portable computer disk, a hard disk, a random access memory (RAM), a read-only memory (ROM), an erasable programmable read-only memory (EPROM or flash memory), an optical fiber, a portable compact disk read-only memory (CD-ROM), an optical storage device, a magnetic storage device, or any suitable combination of the foregoing.

[0178] To provide interaction with a user, the systems and techniques described herein may be implemented on an electronic device having: a display device (e.g., a CRT (cathode ray tube) or LCD (liquid crystal display) monitor) for displaying information to the user; and a keyboard and a pointing device (e.g., a mouse or trackball) through which the user can provide input to the electronic device. Other types of devices may also be used to provide interaction with the user; for example, the feedback provided to the user may be any form of sensory feedback (e.g., visual feedback, auditory feedback, or tactile feedback); and input from the user may be received in any form (including acoustic input, voice input, or tactile input).

[0179] The systems and techniques described herein may be implemented in a computing system that includes backend components (e.g., as a data server), or a computing system that includes middleware components (e.g., an application server), or a computing system that includes frontend components (e.g., a user computer with a graphical user interface or a web browser through which a user can interact with implementations of the systems and techniques described herein), or a computing system that includes any combination of such backend components, middleware components, or frontend components. The components of the system may be interconnected by any form or medium of digital data communication (e.g., a communication network). Examples of communication networks include: a local area network (LAN), a wide area network (WAN), a blockchain network, and the Internet.

[0180] A computing system may include a client and a server. The client and the server are generally remote from each other and usually interact through a communication network. The client and server relationship is generated by computer programs running on the corresponding computers and having a client-server relationship with each other. The server may be a cloud server, also known as a cloud computing server or cloud host, which is a host product in the cloud computing service system to solve the defects of difficult management and weak business scalability in traditional physical hosts and VPS services.

[0181] It should be understood that the various forms of processes shown above can be used to reorder, add or delete steps. For example, the steps described in the present invention can be executed in parallel, sequentially or in different orders, as long as the desired results of the technical solution of the present invention can be achieved, and this document does not limit this.

[0182] The above specific implementations do not constitute a limitation on the protection scope of the present invention. It should be understood by those skilled in the art that various modifications, combinations, sub-combinations and substitutions can be made according to design requirements and other factors. Any modification, equivalent substitution and improvement made within the spirit and principle of the present invention should be included in the protection scope of the present invention.

Claims

1. A method for establishing a rock physics plate, characterized in that: include: For a shale gas reservoir for which a rock physics map is to be established, determining the elastic modulus of a rock matrix corresponding to the shale gas reservoir; For the dry pores to be added to the rock matrix, the elastic modulus of a first rock skeleton is obtained according to the elastic modulus of the dry pores and the elastic modulus of the rock matrix, wherein the first rock skeleton is obtained by adding the dry pores to the rock matrix; For the directional cracks to be added to the first rock skeleton, an elastic stiffness matrix of a second rock skeleton is obtained according to the directional weakness of the directional cracks and the elastic modulus of the first rock skeleton, wherein the second rock skeleton is obtained by adding the directional cracks to the first rock skeleton; A rock physics model is established based on the elastic stiffness matrix of the second rock skeleton, and a rock physics plate of the shale gas reservoir is established based on the rock physics model.

2. The method according to claim 1, characterized in that The step of establishing a rock physics plate of the shale gas reservoir according to the rock physics model comprises: According to the well logging data corresponding to the shale gas reservoir, the variation range of the physical property parameters is obtained, and grid division is performed within the variation range of the physical property parameters to obtain a plurality of grid nodes: For each grid node among the plurality of grid nodes, according to the rock physics model and the physical property parameters represented by the grid node, obtaining the elastic parameter represented by the grid node; A rock physics map of the shale gas reservoir is established based on the elastic parameters respectively represented by the multiple grid nodes.

3. The method according to claim 2, characterized in that After establishing the rock physics map of the shale gas reservoir according to the elastic parameters respectively represented by the plurality of grid nodes, the method further includes: Projecting the logging data onto the rock physics plate, so as to adjust the elastic parameters respectively represented by the plurality of grid nodes based on the obtained projection result; The rock physics plate is corrected according to the adjusted elastic parameters respectively represented by the plurality of grid nodes.

4. The method according to claim 2, characterized in that: After establishing the rock physics map of the shale gas reservoir according to the elastic parameters respectively represented by the plurality of grid nodes, the method further includes: Acquiring target data, wherein the target data is the well logging data, or elastic parameters obtained by inverting seismic data corresponding to the shale gas reservoir; Projecting the target data into the rock physics plate to determine a target node matching the target data from each grid node in the rock physics plate; The physical property parameters represented by the target node are used as the physical property parameters corresponding to the target data.

5. The method according to claim 1, characterized in that The step of obtaining the elastic modulus of the first rock skeleton according to the elastic modulus of the dry pores and the elastic modulus of the rock matrix comprises: determining the porosity of the dried pores; According to the elastic modulus of the dry pores and the porosity, and taking the elastic modulus of the rock matrix as an initial value, the elastic modulus of the first rock skeleton is solved.

6. The method according to claim 1, characterized in that The step of establishing a rock physics model according to the elastic stiffness matrix of the second rock skeleton comprises: Obtaining the bulk modulus of the pore fluid in the dry pores and the oriented fractures according to the bulk modulus of water, the bulk modulus of gas and the water saturation corresponding to the shale gas reservoir; Obtaining an elastic stiffness matrix of saturated rock according to the elastic modulus of the rock matrix, the elastic stiffness matrix of the second rock skeleton, and the bulk modulus of the pore fluid, wherein the dry pores and the directional fractures in the saturated rock are filled with saturated pore fluid; A rock physics model is established based on the elastic stiffness matrix of the saturated rock.

7. The method according to claim 1, characterized in that The determining of the elastic modulus of the rock matrix corresponding to the shale gas reservoir comprises: Determining the elastic modulus of the brittle mineral mixture corresponding to the shale gas reservoir; Determining the elastic modulus of an organic-rich clay mineral corresponding to the shale gas reservoir, wherein the organic-rich clay mineral is obtained by adding organic matter to a clay mineral; For a rock matrix composed of the brittle mineral mixture and the organic-rich clay mineral, the elastic modulus of the rock matrix is ​​obtained according to the elastic modulus of the brittle mineral mixture and the elastic modulus of the organic-rich clay mineral.

8. A rock physics plate creation device, characterized in that: include: A rock matrix elastic modulus obtaining module is used to determine the elastic modulus of the rock matrix corresponding to the shale gas reservoir for which the rock physics map is to be established; A rock skeleton elastic modulus obtaining module, for obtaining the elastic modulus of a first rock skeleton according to the elastic modulus of the dry pores to be added to the rock matrix and the elastic modulus of the rock matrix, wherein the first rock skeleton is obtained by adding the dry pores to the rock matrix; An elastic stiffness matrix obtaining module is used to obtain an elastic stiffness matrix of a second rock skeleton according to the directional weakness of the directional cracks to be added to the first rock skeleton and the elastic modulus of the first rock skeleton, wherein the second rock skeleton is obtained by adding the directional cracks to the first rock skeleton; The rock physics plate establishing module is used to establish a rock physics model according to the elastic stiffness matrix of the second rock skeleton, and to establish a rock physics plate of the shale gas reservoir according to the rock physics model.

9. An electronic device, characterized in that: include: at least one processor; as well as a memory communicatively connected to the at least one processor; wherein, The memory stores a computer program executable by the at least one processor, and the computer program is executed by the at least one processor so that the at least one processor performs the rock physical plate creation method according to any one of claims 1 to 7.

10. A computer-readable storage medium, characterized in that: The computer-readable storage medium stores computer instructions, and the computer instructions are used to enable a processor to implement the rock physical plate establishment method as described in any one of claims 1 to 7 when executed.

Images