Rock physical model construction method based on SPC lithofacies parameter constraint

Abstract

The invention relates to the technical field of oil and gas field exploration and development, in particular to a rock physical model construction method based on SPC lithofacies parameter constraint. The method comprises the following steps: acquiring S-type parameters, P-type parameters and C-type parameters through core data and logging data; performing quantitative analysis on the SPC type parameters and the physical property parameters, and determining the influence of the complex lithofacies sedimentary structure on the physical property of the rock; constructing an SPC improved solid mineral matrix model; constructing a dry rock skeleton model based on multiple pore aspect ratio factors; constructing a pore fluid model; and according to the SPC improved solid mineral matrix model, the dry rock skeleton model and the pore fluid model, on the basis of a Gassmann theory, constructing a rock physical model based on SPC lithofacies parameter constraint. According to the method, the adaptability and accuracy of the rock physical model in different lithofacies prediction processes can be improved; and a communication bridge is built for further predicting the lithology, the physical property and the oil-gas possibility of the reservoir.

Description

Technical Field

[0001] The present invention relates to the technical field of oil and gas field exploration and development, and particularly relates to a method for constructing a rock physics model based on SPC lithofacies parameter constraints. Background Art

[0002] Seismic rock physics is a science that studies the rock physical properties related to seismic characteristics and the relationships between these physical properties and seismic responses. It builds a bridge between seismic data and oil and gas characteristics and reservoir parameters, greatly expanding the application fields of seismic technology.

[0003] In the past few decades, great progress has been made in studying the rock physical properties related to seismic exploration and natural seismology. During this period, many theories have been developed and many experiments have been carried out. Many rock physical theories and experimental results have played an important guiding role in promoting earth science and exploration technology.

[0004] In recent years, rock physical technology has played an important role in oilfield exploration and development, promoting the development of oil and gas detection technologies such as time-lapse seismic oil and gas monitoring, seismic lithology identification, and reservoir fluid identification. The rock physics model is one of the main methods for rock physical research. It is a model that evaluates the physical properties based on the basic rock information of geological bodies (such as mineral composition, content, and distribution). It is an important bridge connecting geology, logging, and seismology, and is the basis for carrying out research work such as seismic attribute analysis and fluid inversion. Only a suitable rock physics model can ensure the accuracy and reliability of reservoir geophysical prediction results.

[0005] As more and more oil and gas fields enter the mature stage and many methods and technologies are further improved and developed, including the mapping of porosity and fluid flow, the determination of changes in reservoir pressure, temperature, and saturation, and even the inference of the type of oil and gas-bearing. In such a development process, rock physics modeling is a very important technical means for quantitatively interpreting seismic data.

[0006] CN110909486B discloses a method for establishing an orthotropic shale rock physics model, including the following steps: establishing a rock physics model that can predict the shale stiffness tensor through logging, mud logging, laboratory tests, anisotropic SCA and DEM models, isotropic SCA and DEM models, Schoenberg linear slip model, and Brown-Korringa model.

[0007] CN110824556B discloses a method for establishing a petrophysical model of an unconventional tight sandstone reservoir. The method includes the following steps: Step 1: Establish a rock skeleton model: According to the test temperature, pressure, rock mineral component content, elastic modulus, rock porosity, pore shape, density, and fluid property parameters, simulate the temperature and pressure environment of the core sample test. Based on the equivalent medium theory, establish a rock skeleton model and obtain the equivalent elastic modulus data of the rock skeleton. Step 2: Establish a mixed fluid model: According to the temperature and pressure environment of the core sample test simulated in the first step, use the equivalent modulus formula of the fluid mixture to establish a mixed fluid model and obtain the equivalent bulk modulus data of the fluid mixture. Step 3: Establish an isotropic saturated rock model: According to the equivalent elastic modulus data of the rock skeleton obtained in Step 1 and the equivalent bulk modulus data of the fluid mixture obtained in Step 2, establish an isotropic saturated rock model. The forward modeling of the model obtains the longitudinal wave velocity, transverse wave velocity, and Poisson's ratio data. Step 4: Establish a petrophysical model of the unconventional tight sandstone reservoir.

[0008] CN109655940B discloses a method for modeling an anisotropic petrophysical model of shale, including regarding the shale matrix as a mixture composed of brittle minerals, organic matter, and clay; regarding clay particles as anisotropic elements with a fixed elastic stiffness matrix, and introducing a clay particle orientation index to characterize the degree of oriented arrangement of clay particles; dividing the total pores of shale into three types of pores: brittle pores, clay pores, and organic matter pores; obtaining the elastic parameters of the mixture of porous brittle minerals and porous organic matter; obtaining the elastic parameters of the mixture of brittle minerals and organic matter; obtaining the elastic parameters of the pore clay medium containing bound water; obtaining the equivalent elastic parameters of the shale matrix and the dry rock skeleton of shale; and obtaining the equivalent elastic parameters of fluid-saturated shale based on the equivalent elastic parameters of the shale matrix and the dry rock skeleton of shale, thus completing the construction of the anisotropic petrophysical model of shale in the fluid-saturated state.

[0009] How to construct a practical and accurate petrophysical model is the part of constructing a new petrophysical model for reservoir physical property estimation. Currently, it mainly includes laboratory measurement and theoretical model research. The theoretical model idealizes the actual rock through certain assumptions and establishes a general relationship through the internal physical principles. Some models assume that the pores and grains of the rock are arranged in layers, and some consider the rock as an aggregate composed of grains and pores of a single geometric shape. Judging from the current exploration degree and understanding, these theoretical models often have certain limitations and lack targeted differential characteristics for the identification of complex lithofacies, which to a certain extent increases the multi-solution of later research work such as reservoir physical property prediction and fluid inversion. Summary of the Invention

[0010] To solve the above-mentioned problems, the present invention provides a method for constructing a rock physics model based on SPC lithofacies parameter constraints. The method of the present invention can improve the adaptability and accuracy of the rock physics model in the prediction process of different lithofacies.

[0011] To achieve the above object, the present invention adopts the following technical solutions:

[0012] The present invention provides a method for constructing a rock physics model based on SPC lithofacies parameter constraints, which includes the following steps:

[0013] Obtain S-type parameters, P-type parameters, and C-type parameters through core data and logging data;

[0014] Conduct quantitative analysis on SPC lithofacies parameters and physical properties parameters to clarify the influence of complex lithofacies sedimentary fabrics on rock physical properties;

[0015] Construct an improved SPC solid mineral matrix model;

[0016] Use the Self-Consistence model to construct a dry rock skeleton model based on multiple pore aspect ratio factors;

[0017] Use the Batzle and Wang equation to construct a pore fluid model;

[0018] According to the above-mentioned improved SPC solid mineral matrix model, dry rock skeleton model, and pore fluid model, based on the Gassmann theory, construct a rock physics model based on SPC lithofacies parameter constraints.

[0019] Furthermore, the S-type parameter is a fluid-insensitive elastic parameter related to the shear modulus G; the P-type parameter is a parameter that can reflect the coupling of the solid medium and structure and the fluid and is sensitive to the compressibility of the rock; the C-type parameter is a combination of the S-type parameter and the P-type parameter. The present invention conducts SPC characteristic parameter analysis to quantitatively analyze the influence of complex lithofacies sedimentary fabrics on rock physical properties.

[0020] Even further, the S-type parameters include the shear wave velocity Vs, Lame coefficient λ, shear modulus μ, and shear wave impedance Zs, etc. These parameters are insensitive to the presence and content of fluids in the rock pores. From the perspective of rock deformation, these parameters mainly respond to shear, and the shear parameters of different rocks are different. Therefore, analyzing the S parameters helps to distinguish lithologies.

[0021] Furthermore, the P-type parameters include bulk modulus K, P-wave velocity Vp, P-wave impedance Zp, Lame coefficient, etc. In terms of physical nature, bulk modulus K, P-wave velocity Vp, P-wave impedance Zp, Lame coefficient λ, etc. reflect the coupling of solid media and structures, and fluids. From the perspective of rock deformation, these parameters are sensitive to the compressibility of rocks.

[0022] Furthermore, the C-type parameters include Vp / Vs, Poisson's ratio, Zp 2 -cZs 2 etc.

[0023] Furthermore, the SPC improved solid mineral matrix model is expressed by the following formula:

[0024]

[0025] where K m is the bulk modulus of the rock matrix, K SPC represents the Voigt upper bound of the effective bulk modulus of the complex lithofacies, K R represents the Reuss lower bound of the effective bulk modulus of the rock, μ m is the shear modulus of the rock matrix, μ spc represents the Voigt upper bound of the effective shear modulus of the complex lithofacies, μ R represents the Reuss lower bound of the effective shear modulus of the rock.

[0026] The Voigt-Reuss-Hill model is to estimate the Voigt upper bound and Reuss lower bound of the effective elastic modulus of the rock medium under the same strain state and the same stress state respectively, given the relative content and elastic modulus of each phase medium composing the rock, and calculate the effective elastic modulus of the rock by taking the arithmetic mean of the two. The model is shown in the following formula:

[0027]

[0028] where:

[0029] f i —— the volume fraction of the i-th component of the rock, %;

[0030] M i —— the elastic modulus of the i-th medium, MPa;

[0031] —— the Voigt upper bound of the effective elastic modulus of the rock mixture, MPa;

[0032] —— the Reuss lower bound of the effective elastic modulus of the rock mixture, MPa;

[0033] After the components of the rock minerals are determined, the elastic moduli K m and μ m can be obtained by the Voigt-Reuss-Hill (VRH) averaging as follows:

[0034]

[0035] where K m is the bulk modulus of the rock matrix, K V represents the Voigt upper bound of the effective bulk modulus of the rock, K R represents the Reuss lower bound of the effective bulk modulus of the rock, μ m is the shear modulus of the rock matrix, μ V represents the Voigt upper bound of the effective shear modulus of the rock, μ R represents the Reuss lower bound of the effective shear modulus of the rock.

[0036] This averaging has no theoretical basis or physical meaning. This model is more suitable for calculating the effective bulk modulus of mineral components and its possible maximum upper and lower limits, and is not suitable for obtaining the total bulk modulus, shear modulus of the rock, and the case of gas-saturated rock.

[0037] The present invention introduces an S-type parameter that can reflect the shear performance, a P-type parameter that reflects the coupling of solid media, structures, and fluids, and a C-type combined parameter for complex lithofacies types such as glutenite bodies, beach-bar sands, and turbidites, to realize the construction of an SPC improved (V-R-H) solid mineral matrix model.

[0038] Furthermore, the dry rock skeleton model based on multiple pore aspect ratio factors is expressed by the following formula:

[0039]

[0040] where:

[0041]

[0042] K dry represents the bulk modulus of the dry rock, μ dry represents the shear modulus of the dry rock, the subscript m represents the background material, i represents the inclusion material, x i represents the volume content of the inclusion material, Φ and ξ are aspect ratio factors representing pore fillers or dry pore shapes, related to the mineral matrix modulus and pore aspect ratio, is the porosity, and α is the consolidation index.

[0043] When only the elastic moduli and volume fractions of the mineral components of the rock are known, without knowing the geometric details of how the components are combined with each other, most methods first use the theoretical solution of the elastic deformation of a single inclusion of one material added to an infinite background medium of another material, and then use different methods to estimate the equivalent modulus when the inclusions are distributed in a certain way. The Self-Consistence model still applies the mathematical solution of the deformation of the inclusions isolated from each other. The model is shown as follows:

[0044]

[0045] For complex lithofacies types such as glutenite bodies, beach-bar sands, and turbidites, the pore space contains both easily deformable shale pores (low pore aspect ratio) and hard sand pores (high pore aspect ratio). Therefore, the present invention introduces the pore aspect ratio factors Φ and ξ, and constructs a dry rock skeleton model based on multiple pore aspect ratio factors.

[0046] Furthermore, the pore fluid model is:

[0047]

[0048] where K f is the bulk modulus of the fluid, β is the compressibility, P is the pressure, and V is the volume.

[0049] The bulk modulus of the fluid in the reservoir increases with the increase of the formation pressure and decreases with the increase of the temperature. For oil and water, the influence of pressure and temperature can be ignored, but for gas, pressure and temperature will have a great influence on its bulk modulus. Therefore, the influence of pressure and temperature cannot be ignored.

[0050] Through the analysis and test data of oil samples, gas samples, and formation water, the elastic parameters of the fluid can be calculated by the Batzle and Wang equation (1992). The specific calculation formula:

[0051] ① Formation water

[0052] The density of formation water with salinity S is:

[0053] ρ B = ρ W + S{0.668 + 0.44S + 10 -6 [300P - 2400PS + T(80 + 3T - 3300S - 13P + 47PS)]}

[0054] The acoustic velocity of formation water with salinity S is:

[0055] V B = V W + S(1170 - 9.6T + 0.055T2 -8.5x10 -5 T 3 +2.6P - 0.0029TP - 0.0476P 2 ) + S 1.5 (780 - 10P + 0.16P 2 ) - 1820S 2

[0056] ② Gas

[0057]

[0058]

[0059] ③ Oil

[0060]

[0061] ④ Live oil

[0062]

[0063] ⑤ Dead oil

[0064] ρ P = ρ 0 +(0.00277P - 1.71×10 -7 P 3 )(ρ 0 - 1.15) 2 +3.49×10 -4 P

[0065]

[0066] From this, the bulk modulus K of the fluid can be obtained f .

[0067] Furthermore, the rock physics model based on SPC lithofacies parameter constraints is expressed by the following formula:

[0068]

[0069] where V p is the P-wave velocity of the rock, V s is the S-wave velocity of the rock, K sat is the bulk modulus of the saturated rock, μ sat is the shear modulus of the saturated rock, and ρ sat is the density of the saturated rock.

[0070] The Gassmann (1951) theory assumes the macroscopic elastic characteristics of the subsurface medium as a homogeneous isotropic, non-viscous fluid-saturated rock, and establishes a quantitative relationship between the modulus of the saturated rock and the dry rock skeleton, mineral matrix, porosity, and equivalent fluid modulus:

[0071]

[0072] Among them, K dry is the bulk modulus of the dry rock skeleton, K mat is the bulk modulus of the rock matrix, and K f is the bulk modulus of the mixed fluid.

[0073] Combined with the SPC improved (V-R-H) solid mineral matrix model, the dry rock skeleton model based on multiple pore aspect ratio factors, and the pore fluid model, the above-mentioned rock physics model constrained by SPC lithofacies parameters is finally constructed.

[0074] Compared with the prior art, the present invention has the following advantages:

[0075] The method of the present invention comprehensively considers multiple models such as the solid mineral matrix, dry rock skeleton, and pore fluid model, realizes the construction of a rock physics model based on lithofacies constraints, and improves the adaptability and accuracy of the rock physics model in the prediction process of different lithofacies. The rock physics model constructed by the present invention is more targeted and applicable to complex lithofacies such as glutenite bodies, beach-bar sands, and turbidites. BRIEF DESCRIPTION OF THE DRAWINGS

[0076] Figure 1 is a flowchart of the method for constructing a rock physics model constrained by SPC lithofacies parameters according to an embodiment of the present invention.

[0077] Figure 2 is a relationship diagram between SPC parameters and physical properties of different lithologies according to Embodiment 2 of the present invention.

[0078] Figure 3 Schematic diagram of SPC parameter templates for different facies belts in Embodiment 2 of the present invention. DETAILED DESCRIPTION OF THE EMBODIMENTS

[0079] It should be noted that the following detailed description is exemplary and is intended to provide further explanation of the present invention. Unless otherwise specified, all technical and scientific terms used herein have the same meaning as commonly understood by those of ordinary skill in the technical field to which the present invention belongs.

[0080] It should be noted that the terms used herein are only for describing specific embodiments and are not intended to limit the exemplary embodiments according to the present invention. As used herein, unless the context clearly indicates otherwise, the singular form is also intended to include the plural form. In addition, it should also be understood that when the terms "comprising" and "including" are used in this specification, they indicate the presence of features, steps, operations, and combinations thereof.

[0081] In order to enable those skilled in the art to more clearly understand the technical solution of the present invention, the technical solution of the present invention will be described in detail below in conjunction with specific embodiments.

[0082] The SPC lithofacies parameters described in the present invention refer to the general term of S-type parameters, P-type parameters, and C-type parameters.

[0083] Example 1

[0084] A method for constructing a rock physics model based on SPC lithofacies parameter constraints, which includes the following steps:

[0085] Step 1. Obtain S-type parameters, P-type parameters, and C-type parameters through core data and logging data;

[0086] The S-type parameter is fluid-insensitive and an elastic parameter related to the shear modulus G. The S-type parameters include shear wave velocity Vs, Lame coefficient λ, shear modulus μ, and shear wave impedance Zs, etc. These parameters are insensitive to the presence and content of fluids in rock pores. From the perspective of rock deformation, these parameters mainly respond to shear, and the shear parameters of different rocks are different. Therefore, analyzing the S parameter helps to distinguish lithologies.

[0087] The P-type parameter is a parameter that can reflect the coupling of solid media and structures and fluids and is relatively sensitive to the compressibility of rocks. The P-type parameters include bulk modulus K, longitudinal wave velocity Vp, longitudinal wave impedance Zp, and Lame coefficient, etc. In terms of physical essence, the bulk modulus K, longitudinal wave velocity Vp, longitudinal wave impedance Zp, Lame coefficient λ, etc. reflect the coupling of solid media and structures and fluids. From the perspective of rock deformation, this type of parameter is relatively sensitive to the compressibility of rocks.

[0088] The C-type parameter is a combination of the S-type parameter and the P-type parameter. The C-type parameters include Vp / Vs, Poisson's ratio, Zp2 - cZs2, etc.

[0089] Conduct quantitative analysis on SPC lithofacies parameters and physical properties parameters to clarify the influence of complex lithofacies sedimentary fabrics on rock physical properties.

[0090] Step 2. Construct an SPC improved solid mineral matrix model

[0091] The SPC improved solid mineral matrix model is represented by the following formula:

[0092]

[0093] where K m is the bulk modulus of the rock matrix, K SPC represents the Voigt upper bound of the effective bulk modulus of the complex lithofacies, K R represents the Reuss lower bound of the effective bulk modulus of the rock, μ m is the shear modulus of the rock matrix, μ spc represents the Voigt upper bound of the effective shear modulus of the complex lithofacies, μ R represents the Reuss lower bound of the effective shear modulus of the rock.

[0094] Step 3. Construct a dry rock skeleton model based on multiple pore aspect ratio factors using the Self-Consistence model:

[0095] The dry rock skeleton model based on multiple pore aspect ratio factors is expressed by the following formula:

[0096]

[0097] where:

[0098]

[0099] K dry represents the bulk modulus of the dry rock, μ dry represents the shear modulus of the dry rock, the subscript m represents the background material, i represents the inclusion material, x i represents the volume content of the inclusion material, Φ and ξ are aspect ratio factors describing the pore filler or the shape of the dry pores, related to the mineral matrix modulus and the pore aspect ratio, is the porosity, and α is the consolidation index.

[0100] Step 4. Construct a pore fluid model using the Batzle and Wang equation

[0101] The pore fluid model is:

[0102]

[0103] where K f is the bulk modulus of the fluid, β is the compressibility, P is the pressure, and V is the volume.

[0104] Through the analysis and test data of oil samples, gas samples, and formation water, the elastic parameters of the fluid can be calculated by the Batzle and Wang equation (1992). The specific calculation formula:

[0105] ① Formation water

[0106] The density of formation water with salinity s is:

[0107] ρ B = ρ W + S{0.668 + 0.44S + 10 -6 [300P - 2400PS + T(80 + 3T - 3300S - 13P + 47PS)]}

[0108] The acoustic velocity of formation water with salinity S is:

[0109] V B = V W + S(1170 - 9.6T + 0.055T 2 - 8.5×10 -5 T 3 + 2.6P - 0.0029TP - 0.0476P 2 ) + S 1.5 (780 - 10P + 0.16P 2 ) - 1820S 2

[0110] ② Gas

[0111]

[0112]

[0113] ③ Oil

[0114]

[0115] ④ Live oil

[0116]

[0117] ⑤ Dead oil

[0118] ρ P = ρ 0 + (0.00277P - 1.71×10 -7 P 3 )(ρ 0 - 1.15) 2 + 3.49×10 -4 P

[0119]

[0120] From this, the bulk modulus K of the fluid can be obtained f .

[0121] Step 5. Combine the SPC improved (V - R - H) solid mineral matrix model, the dry rock skeleton model based on multiple pore aspect ratio factors, and the pore fluid model to finally construct a rock physics model constrained by SPC lithofacies parameters.

[0122] Gassmann (1951) theory assumes the macroscopic elastic characteristics of the underground medium as a homogeneous isotropic, non - viscous fluid - saturated rock, and establishes a quantitative relationship between the modulus of the saturated rock and the dry rock skeleton, mineral matrix, porosity, and equivalent fluid modulus:

[0123]

[0124] Among them, Kdr y is the bulk modulus of the dry rock skeleton, K mat is the bulk modulus of the rock matrix, K f is the bulk modulus of the mixed fluid.

[0125] Combine the SPC improved (V - R - H) solid mineral matrix model, the dry rock skeleton model based on multiple pore aspect ratio factors, and the pore fluid model to finally construct the above - mentioned rock physics model constrained by SPC lithofacies parameters:

[0126]

[0127] Among them, V p is the P - wave velocity of the rock, V s is the S - wave velocity of the rock, K sat is the bulk modulus of the saturated rock, μ sat is the shear modulus of the saturated rock, ρ sat is the density of the saturated rock.

[0128] Example 2

[0129] Taking glutenite as an example, to quantitatively analyze the influence of gravel type, content, and sand - grade particle debris fabric of glutenite on rock physical properties, and objectively analyze the differences in rock physical properties between glutenite bodies and conventional sandstones, a method for constructing a rock physics model constrained by SPC lithofacies parameters for glutenite is carried out.

[0130] Step 1. Obtain S - type parameters, P - type parameters, and C - type parameters through core data and logging data;

[0131] Step 2. Conduct a quantitative analysis of SPC lithofacies parameters and physical properties parameters to clarify the influence of complex lithofacies sedimentary fabrics on rock physical properties ( Figure 2 ).

[0132] According to the differences in rock structures such as gravel content and physical properties, the sedimentation of subaqueous fan sandy conglomerate bodies near the shore is divided into three lithofacies zones: fan tip, fan middle, and fan root. The fan root facies zone is mainly composed of conglomerate with a relatively high gravel content and poor physical properties. The middle fan subfacies is mainly composed of pebbly sandstone and gravelly sandstone with a medium gravel content and good physical properties. The fan tip is mainly composed of mudstone.

[0133] Step 3: According to the process of Example 1, establish dry rock skeleton models and pore fluid models for sandy conglomerates of different lithofacies, and finally construct the rock physics model based on SPC lithofacies parameter constraints described above.

[0134] The velocity intervals and the resulting elastic parameters of different lithofacies zones are different. The longitudinal wave velocity of the conglomerate facies with the highest gravel content is the highest, and the ratio of longitudinal wave velocity to transverse wave velocity is also the largest. The longitudinal wave velocity of the sandstone in the middle fan - fan tip facies zone with a low gravel content is relatively lower than that of the gravel - containing lithofacies zone, and the ratio of longitudinal wave velocity to transverse wave velocity is the smallest( Figure 3 ). Thus, elastic parameter templates for different facies zones with varying porosities can be established. Based on this, the elastic parameter values varying with gravel content under the corresponding porosity conditions of different facies zones can be directly obtained. Further, by analyzing the changes in their elastic parameters, the distribution law of lithofacies zones can be inferred inversely.

[0135] The above embodiments are preferred embodiments of the present invention, but the embodiments of the present invention are not limited by the above embodiments. Any other changes, modifications, substitutions, combinations, and simplifications made without departing from the spirit and principle of the present invention shall be equivalent replacement methods and are all included in the protection scope of the present invention.

Claims

1. A method for constructing a rock physics model constrained by SPC petrophysical parameters, characterized in that, it includes the following steps: Obtain S-type parameters, P-type parameters, and C-type parameters through core data and logging data; the S-type parameters are fluid-insensitive elastic parameters related to the shear modulus G; the P-type parameters are parameters that can reflect the coupling of solid media and structures and fluids and are sensitive to the compressibility of rocks; the C-type parameters are a combination of S-type parameters and P-type parameters; Conduct quantitative analysis on SPC petrophysical parameters and physical property parameters to clarify the influence of complex lithofacies sedimentary fabrics on rock physical properties; Construct an SPC improved solid mineral matrix model; Construct a dry rock skeleton model based on multiple pore aspect ratio factors using the Self-Consistence model; Construct a pore fluid model using the Batzle and Wang equation; According to the above-mentioned SPC improved solid mineral matrix model, dry rock skeleton model, and pore fluid model, construct a rock physics model constrained by SPC petrophysical parameters based on the Gassmann theory.

2. The method for constructing a rock physics model constrained by SPC petrophysical parameters according to claim 1, characterized in that, the SPC improved solid mineral matrix model is expressed by the following formula: Among them, K m is the bulk modulus of the rock matrix, K SPC represents the Voigt upper bound of the effective bulk modulus of the complex rock facies, K R represents the Reuss lower bound of the effective bulk modulus of the rock, μ m is the shear modulus of the rock matrix, μ spc represents the Voigt upper bound of the effective shear modulus of the complex rock facies, μ R represents the Reuss lower bound of the effective shear modulus of the rock.

3. The method for constructing a rock physics model constrained by SPC petrophysical parameters according to claim 1, characterized in that, the dry rock skeleton model based on multiple pore aspect ratio factors is expressed by the following formula: where: K dry represents the bulk modulus of dry rock, μ dry represents the shear modulus of dry rock, the subscript m represents the background material, i represents the inclusion material, x i represents the volume content of the inclusion material, Φ, ξ are the aspect ratio factors representing the pore filler or the dry pore shape, related to the mineral matrix modulus and the pore aspect ratio, is the porosity, α is the consolidation index.

4. The method for constructing a rock physics model constrained by SPC petrophysical parameters according to claim 1, characterized in that, the pore fluid model is: where K f is the bulk modulus of the fluid, β is the compressibility, P is the pressure, and V is the volume.

5. The method for constructing a rock physics model constrained by SPC petrophysical parameters according to claim 1, characterized in that, the rock physics model constrained by SPC petrophysical parameters is expressed by the following formula: Among them, V p is the longitudinal wave velocity of the rock, V s is the shear wave velocity of the rock, K sat is the bulk modulus of the saturated rock, μ sat is the shear modulus of the saturated rock, ρ sat is the density of the saturated rock.

6. The method for constructing a rock physics model constrained by SPC petrophysical parameters according to claim 1, characterized in that, the S-type parameters include the shear wave velocity Vs, Lame coefficient λ, shear modulus μ, and shear wave impedance Zs.

7. The method for constructing a rock physics model constrained by SPC petrophysical parameters according to claim 1, characterized in that, the P-type parameters include the bulk modulus K, longitudinal wave velocity Vp, longitudinal wave impedance Zp, and shear modulus μ.

8. The method for constructing a rock physics model constrained by SPC petrophysical parameters according to claim 1, characterized in that, The C-type parameters include Vp / Vs, Poisson's ratio, and Zp 2 -cZs 2 .

Images