A method and system for seismic petrophysical modeling of deep high temperature high pressure hot dry rock

Abstract

The application discloses a deep high-temperature high-pressure dry hot rock seismic rock physics modeling method and system, relates to the field of seismic exploration, determines the mineral composition and content of the dry hot rock matrix, and calculates the elastic modulus of the dry hot rock matrix; the differential equivalent medium model is adopted to calculate the elastic modulus of the dry hot rock skeleton containing hard pores; the volume modulus of the rock skeleton of a rock sample under different pressure conditions is tested, and the effective pressure constant is estimated; based on the pore space stiffness theory, the elastic modulus of the dry hot rock skeleton containing cracks affected by pressure is calculated; the volume modulus of the pore fluid affected by temperature and pressure is calculated; and the elastic modulus and density of the dry hot rock containing fluid are calculated. The application quantitatively constructs the analytical relationship between the temperature and pressure conditions and the elastic response and the crack content, provides a theoretical model basis for seismic prediction of the physical property parameters of the dry hot rock reservoir, can significantly improve the seismic prediction precision and resource evaluation efficiency of the dry hot rock reservoir, and provides technical support for the exploration and development of the dry hot rock geothermal resources.

Description

Technical Field

[0001] This invention relates to the field of seismic exploration technology for unconventional geothermal resources, and more specifically to a seismic rock physics modeling method and system for deep, high-temperature, high-pressure dry hot rocks. Background Technology

[0002] Currently, hot dry rock, as an important unconventional geothermal resource, has attracted much attention in recent years due to its abundant reserves and potential for sustainable development and utilization. The exploration and development of hot dry rock geothermal resources requires accurate characterization of reservoir physical parameters (such as porosity, fracture density, and fluid saturation), and seismic rock physics modeling is a crucial technical means to achieve refined characterization of hot dry rock reservoirs. By combining seismic wave propagation characteristics with rock physics models, seismic rock physics modeling can quantitatively predict reservoir physical parameters, providing theoretical support for the evaluation and development of hot dry rock reservoirs.

[0003] However, deep hot dry rocks possess significantly different physical and elastic characteristics compared to conventional geothermal reservoirs, such as complex pore structures and high-temperature, high-pressure environments. These characteristics significantly affect the propagation of seismic waves in hot dry rocks, thereby impacting seismic exploration and prediction of such geothermal reservoirs. Existing seismic rock physics modeling methods are primarily designed for conventional reservoirs and are significantly inadequate in describing the influence of complex fracture and pore structures on seismic response. Furthermore, the impact of high-temperature, high-pressure environments on the mechanical properties and elastic parameters of reservoir rocks has not been fully considered, further limiting the applicability of traditional modeling methods.

[0004] Therefore, for high-temperature and high-pressure dry hot rock reservoirs, it is urgent to carry out research on seismic rock physics modeling methods that take into account the effects of high-temperature and high-pressure conditions and complex pore structures. Summary of the Invention

[0005] In view of this, the present invention provides a seismic rock physics modeling method and system for deep, high-temperature and high-pressure dry hot rocks. By quantitatively constructing the analytical relationship between temperature and pressure conditions, fracture content and elastic response, it provides a theoretical model basis for seismic prediction of the physical parameters of dry hot rock reservoirs, which can significantly improve the accuracy of seismic prediction and resource evaluation efficiency of dry hot rock reservoirs, and provide technical support for the exploration and development of dry hot rock geothermal resources.

[0006] To achieve the above objectives, the present invention adopts the following technical solution:

[0007] A seismic rock physics modeling method for deep, high-temperature, high-pressure dry hot rocks, comprising:

[0008] Step 1: Determine the mineral composition and content of the hot dry rock matrix, and calculate the elastic modulus of the hot dry rock matrix;

[0009] Step 2: Using a differential equivalent medium model, calculate the elastic modulus of the hard-pore dry hot rock skeleton;

[0010] Step 3: Test the bulk modulus of the rock skeleton under different pressure conditions and estimate the effective pressure constant;

[0011] Step 4: Based on the pore space stiffness theory, and on the obtained elastic modulus of the dry hot rock matrix, the elastic modulus of the dry hot rock skeleton with hard pores, and the effective pressure constant, calculate the elastic modulus of the dry hot rock skeleton with fractures affected by pressure.

[0012] Step 5: Calculate the bulk modulus of the pore fluid affected by temperature and pressure;

[0013] Step 6: Based on the bulk modulus of the pore fluid affected by temperature and pressure and the elastic modulus of the fractured dry hot rock skeleton affected by pressure, calculate the elastic modulus and density of the fluid-bearing dry hot rock under high temperature and high pressure conditions.

[0014] Optionally, step one specifically includes:

[0015] Core X-ray spectroscopy analysis was used to determine the mineral composition and percentage content of the hot dry rock. The Voigt-Reuss-Hill model was used to calculate the bulk modulus and shear modulus of the hot dry rock matrix.

[0016]

[0017] The density of the hot dry rock is calculated as follows:

[0018]

[0019] In the formula, K m μ m and ρ m These represent the bulk modulus, shear modulus, and density of the dry hot rock matrix, K. mi μ mi ρ mi and V mi These represent the bulk modulus, shear modulus, density, and percentage content of the i-th mineral, respectively, while N represents the total number of mineral types.

[0020] Optionally, step two specifically includes:

[0021] The pores of hot dry rocks are divided into fractures and hard pores. A differential equivalent medium model is used to add hard pores to the hot dry rock matrix, forming a hot dry rock skeleton with hard pores. The elastic modulus of the hot dry rock skeleton with hard pores is then calculated.

[0022]

[0023] In the formula, K * and μ *These are the bulk modulus and shear modulus of a hard-porous, dry, hot rock skeleton, φ st For hard pores, K2 and μ2 are the bulk modulus and shear modulus of the inclusions, respectively, and P* is the porosity. 2 and Q* 2 All of these are geometric factors of the inclusions.

[0024] Optionally, step three specifically includes:

[0025] Twenty standard dry hot rock samples were prepared, and the total porosity φ of each sample was measured. t|k Where k = 1, 2, 3, ..., 20, the bulk modulus K of each rock sample was tested under effective pressures P = 5, 10, 20, 30, 40, and 50 MPa. d|k,p , where k = 1, 2, 3, ..., 20; p = 5, 10, 20, 30, 40, 50;

[0026] Based on the bulk modulus and porosity measurements of 20 rock samples, and using the bulk modulus measurements under each effective pressure condition, the least squares method was used to fit the porosity and bulk modulus to obtain the pore stiffness ratio k. The relationship used for fitting is as follows:

[0027]

[0028] For the measured values ​​under all effective pressure conditions, porosity and bulk modulus were fitted to obtain the variation of porosity stiffness ratio k with effective pressure. p Where p = 5, 10, 20, 30, 40, 50, the least squares method is further used to fit the pore stiffness ratio and the effective pressure to obtain the effective pressure constants G and M. The relationship used for fitting is as follows:

[0029] k = G + Mln(P).

[0030] Optionally, step four specifically includes:

[0031] Treating the entire hard-pore-bearing hot dry rock skeleton as a new matrix, a stiffness-pore space model was used to incorporate fractures into the new matrix, yielding the bulk modulus and shear modulus of the fractured hot dry rock skeleton under pressure.

[0032]

[0033] In the formula, K dry and μ dry These are the bulk modulus and shear modulus of a fractured, dry, hot rock skeleton, K. * and μ * These are the bulk modulus and shear modulus of a hard-porous, dry, hot rock skeleton, φ t and φ stThese represent total porosity and hard pore porosity, respectively; G and M are effective pressure constants; P is the effective pressure; and V is the effective porosity. ck This represents the fracture content.

[0034] Optionally, step five specifically includes:

[0035] Calculate the bulk modulus and density of pore fluids affected by temperature and pressure.

[0036]

[0037] in,

[0038]

[0039] γ o =0.85+5.6 / (P) pr +2)+27.1 / (T pr +3.5) 2 -8.7exp[-0.65(P pr +1)];

[0040] P pr =P / (4.892-0.405G);

[0041] T pr =T a / (94.72+170.75G);

[0042] In the formula, G is the specific gravity of air, and T is the specific gravity of air. a Where T is the absolute temperature, P is the effective pressure, and T is the effective pressure. pr P is the critical temperature. pr Let R be the critical pressure and R be the gas constant.

[0043] Optionally, step six specifically includes:

[0044] The bulk modulus and shear modulus of fluid-bearing dry hot rock were calculated using the Biot-Gassmann equation under preset temperature and pressure conditions:

[0045]

[0046] μsat = μdry;

[0047] In the formula, K sat and μ sat These are the bulk modulus and shear modulus of hot, dry rock containing fluid, respectively, K. fl The bulk modulus of the pore fluid;

[0048] Calculate the elastic parameters of hot, dry rock containing fluid:

[0049]

[0050] ρ=ρ m (1-φ t )+ρ fl φ t ;

[0051] In the formula, V P V represents the longitudinal wave velocity of hot dry rock. S ρ is the shear wave velocity of hot dry rock, and ρ is the density of hot dry rock. m and ρ fl These represent the densities of the dry, hot rock matrix and the pore fluid, respectively.

[0052] A seismic rock physics modeling system for deep, high-temperature, high-pressure dry hot rocks includes:

[0053] The first calculation module determines the mineral composition and content of the hot dry rock matrix and calculates the elastic modulus of the hot dry rock matrix.

[0054] The second calculation module uses a differential equivalent medium model to calculate the elastic modulus of a hard-pore dry hot rock skeleton.

[0055] The test estimation module measures the bulk modulus of the rock skeleton under different pressure conditions and estimates the effective pressure constant.

[0056] The third calculation module, based on the pore space stiffness theory, calculates the elastic modulus of the fractured dry hot rock skeleton under pressure, based on the obtained elastic modulus of the dry hot rock matrix, the elastic modulus of the dry hot rock skeleton with hard pores, and the effective pressure constant.

[0057] The fourth calculation module calculates the bulk modulus of pore fluids affected by temperature and pressure.

[0058] The fifth calculation module calculates the elastic modulus and density of fluid-bearing dry hot rock under high temperature and high pressure conditions, based on the bulk modulus of pore fluid affected by temperature and pressure and the elastic modulus of fractured dry hot rock skeleton affected by pressure.

[0059] As can be seen from the above technical solution, compared with the prior art, the present invention discloses a seismic rock physics modeling method and system for deep, high-temperature, high-pressure, dry hot rocks, which has the following technical effects:

[0060] 1. A rock physics modeling method considering the complex pore structure of hot dry rocks is proposed. Compared with existing modeling methods, it can effectively describe the elastic characteristics of hot dry rock reservoirs with complex fracture and pore structure development.

[0061] 2. A rock physics modeling method that considers the effects of temperature and pressure is proposed. Compared with existing modeling methods, it can describe the elastic response characteristics of hot dry rocks under the influence of high temperature and high pressure environments, and improve the accuracy of seismic rock physics modeling of hot dry rock reservoirs. Attached Figure Description

[0062] To more clearly illustrate the technical solutions in the embodiments of the present invention or the prior art, the drawings used in the description of the embodiments or the prior art will be briefly introduced below. Obviously, the drawings described below are only embodiments of the present invention. For those skilled in the art, other drawings can be obtained based on the provided drawings without creative effort.

[0063] Figure 1 A flowchart illustrating the method provided by the present invention;

[0064] Figure 2(a) shows the curves of longitudinal wave velocity as a function of porosity under different effective pressure conditions provided by the present invention; Figure 2(b) shows the curves of transverse wave velocity as a function of porosity under different effective pressure conditions provided by the present invention.

[0065] Figure 3(a) shows the curves of longitudinal wave velocity as a function of porosity under different temperature conditions provided by the present invention; Figure 3(b) shows the curves of transverse wave velocity as a function of porosity under different temperature conditions provided by the present invention.

[0066] Figure 4(a) shows the curves of longitudinal wave velocity as a function of porosity under different fracture contents provided by the present invention; Figure 4(b) shows the curves of transverse wave velocity as a function of porosity under different fracture contents provided by the present invention. Detailed Implementation

[0067] The technical solutions of the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some embodiments of the present invention, and not all embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those skilled in the art without creative effort are within the scope of protection of the present invention.

[0068] This invention discloses a seismic rock physics modeling method for deep, high-temperature, high-pressure, dry hot rocks, such as... Figure 1 As shown, it includes:

[0069] Step 1: Determine the mineral composition and content of the hot dry rock matrix, and calculate the elastic modulus of the hot dry rock matrix;

[0070] Step 2: Using a differential equivalent medium model, calculate the elastic modulus of the hard-pore dry hot rock skeleton;

[0071] Step 3: Test the bulk modulus of the rock skeleton under different pressure conditions and estimate the effective pressure constant;

[0072] Step 4: Based on the pore space stiffness theory, and on the obtained elastic modulus of the dry hot rock matrix, the elastic modulus of the dry hot rock skeleton with hard pores, and the effective pressure constant, calculate the elastic modulus of the dry hot rock skeleton with fractures affected by pressure.

[0073] Step 5: Calculate the bulk modulus of the pore fluid affected by temperature and pressure;

[0074] Step 6: Based on the bulk modulus of the pore fluid affected by temperature and pressure and the elastic modulus of the fractured dry hot rock skeleton affected by pressure, calculate the elastic modulus and density of the fluid-bearing dry hot rock under high temperature and high pressure conditions.

[0075] In one specific embodiment, step one specifically includes:

[0076] Core X-ray spectroscopy analysis was used to determine the mineral composition and percentage content of the hot dry rock. The Voigt-Reuss-Hill model was used to calculate the bulk modulus and shear modulus of the hot dry rock matrix.

[0077]

[0078] The density of the hot dry rock is calculated as follows:

[0079]

[0080] In the formula, K m μ m and ρ m These represent the bulk modulus, shear modulus, and density of the dry hot rock matrix, K. mi μ mi ρ mi and V mi These represent the bulk modulus, shear modulus, density, and percentage content of the i-th mineral, respectively, while N represents the total number of mineral types.

[0081] In one specific embodiment, step two specifically includes:

[0082] The pores of hot dry rocks are divided into fractures and hard pores. A differential equivalent medium model is used to add hard pores to the hot dry rock matrix, forming a hot dry rock skeleton with hard pores. The elastic modulus of the hot dry rock skeleton with hard pores is then calculated.

[0083]

[0084] In the formula, K * and μ * These are the bulk modulus and shear modulus of a hard-porous, dry, hot rock skeleton, φ st For hard pores, K2 and μ2 are the bulk modulus and shear modulus of the inclusions, respectively, and P* is the porosity. 2 and Q*2 All of these are geometric factors of the inclusions.

[0085] In this embodiment, the rock physical parameters are set as follows: under standard conditions (temperature T = 20℃, pressure P = 0.1MPa), K m =72 GPa, μ m =43 GPa, K fl =1.0 GPa, ρ m =2.71g / cm 3 , ρ fl =1.0g / cm 3 .

[0086] In one specific embodiment, step three specifically includes:

[0087] Twenty standard dry hot rock samples were prepared, and the total porosity φ of each sample was measured. t|k (k=1,2,3,...,20), the bulk modulus K of each rock sample was measured under effective pressure P=5,10,20,30,40,50MPa. d|k,p (k=1,2,3,...,20; p=5,10,20,30,40,50);

[0088] Based on the bulk modulus and porosity measurements of 20 rock samples, and using the bulk modulus measurements under each effective pressure condition, the least squares method was used to fit the porosity and bulk modulus to obtain the pore stiffness ratio k. The relationship used for fitting is as follows:

[0089]

[0090] For the measured values ​​under all effective pressure conditions, porosity and bulk modulus were fitted to obtain the variation of porosity stiffness ratio k with effective pressure. p Where p = 5, 10, 20, 30, 40, 50, the least squares method is further used to fit the pore stiffness ratio and the effective pressure to obtain the effective pressure constants G and M. The relationship used for fitting is as follows:

[0091] k = G + Mln(P).

[0092] In this embodiment, the effective pressure constants obtained by fitting experimental data are G = 0.062 and M = 0.025, respectively.

[0093] In one specific embodiment, step four specifically includes:

[0094] Treating the entire hard-pore-bearing hot dry rock skeleton as a new matrix, a stiffness-pore space model was used to incorporate fractures into the new matrix, yielding the bulk modulus and shear modulus of the fractured hot dry rock skeleton under pressure.

[0095]

[0096] In the formula, K dry and μ dry These are the bulk modulus and shear modulus of a fractured, dry, hot rock skeleton, K. * and μ * These are the bulk modulus and shear modulus of a hard-porous, dry, hot rock skeleton, φ t and φ st These represent total porosity and hard pore porosity, respectively; G and M are effective pressure constants; P is the effective pressure; and V is the effective porosity. ck This represents the fracture content.

[0097] In one specific embodiment, step five specifically includes:

[0098] Calculate the bulk modulus and density of pore fluids affected by temperature and pressure.

[0099]

[0100] in,

[0101]

[0102]

[0103] γ o =0.85+5.6 / (P) pr +2)+27.1 / (T pr +3.5) 2 -8.7exp[-0.65(P pr +1)];

[0104] P pr =P / (4.892-0.405G);

[0105] T pr =T a / (94.72+170.75G);

[0106] In the formula, G is the specific gravity of air, and T is the specific gravity of air. a Where T is the absolute temperature, P is the effective pressure, and T is the effective pressure. pr P is the critical temperature. pr Let R be the critical pressure and R be the gas constant.

[0107] In this embodiment, the pore fluid parameters are set as follows: under standard conditions (temperature T = 20℃, pressure P = 0.1MPa), K fl =1.0 GPa, ρ fl =1.0g / cm 3 .

[0108] In one specific embodiment, step six specifically includes:

[0109] The bulk modulus and shear modulus of fluid-bearing dry hot rock were calculated using the Biot-Gassmann equation under preset temperature and pressure conditions:

[0110]

[0111] μsat = μdry;

[0112] In the formula, K sat and μ sat These are the bulk modulus and shear modulus of hot, dry rock containing fluid, respectively, K. fl The bulk modulus of the pore fluid;

[0113] Calculate the elastic parameters of hot, dry rock containing fluid:

[0114]

[0115] ρ=ρ m (1-φ t )+ρ fl φ t ;

[0116] In the formula, V P V represents the longitudinal wave velocity of hot dry rock. S ρ is the shear wave velocity of hot dry rock, and ρ is the density of hot dry rock. m and ρ fl These represent the densities of the dry, hot rock matrix and the pore fluid, respectively.

[0117] Figure 2(a) shows the variation curves of P-wave velocity of hot dry rock with porosity under different effective pressure conditions in this embodiment, and Figure 2(b) shows the variation curves of S-wave velocity of hot dry rock with porosity under different effective pressure conditions in this embodiment. The effective pressure P = 0.1, 10, 25, 50 MPa, and the porosity ranges from 0 to 0.2. As can be seen from Figures 2(a) and 2(b), both P-wave and S-wave velocities of the hot dry rock increase with increasing effective pressure. Figure 3(a) shows the variation curves of P-wave velocity of hot dry rock with porosity under different temperature conditions in this embodiment, and Figure 3(b) shows the variation curves of S-wave velocity of hot dry rock with porosity under different temperature conditions in this embodiment. The temperature T = 20, 50, 100, 200 °C, and the porosity ranges from 0 to 0.2. As can be seen from Figures 3(a) and 3(b), the P-wave velocity of the hot dry rock decreases with increasing temperature, but the S-wave velocity increases slightly. Figure 4(a) shows the variation curves of P-wave velocity of dry hot rock with porosity under different fracture contents in this embodiment, and Figure 4(b) shows the variation curves of S-wave velocity of dry hot rock with porosity under different fracture contents in this embodiment, where the fracture content V ck=0.1, 0.2, 0.4, 0.6, with porosity ranging from 0 to 0.2. As shown in Figures 4(a) and 4(b), the P-wave and S-wave velocities of hot dry rock decrease significantly with increasing fracture content. Therefore, the method of this invention can effectively describe the elastic response characteristics of hot dry rock reservoirs as a function of pressure, temperature, and fracture content.

[0118] A seismic rock physics modeling system for deep, high-temperature, high-pressure dry hot rocks includes:

[0119] The first calculation module determines the mineral composition and content of the hot dry rock matrix and calculates the elastic modulus of the hot dry rock matrix.

[0120] The second calculation module uses a differential equivalent medium model to calculate the elastic modulus of a hard-pore dry hot rock skeleton.

[0121] The test estimation module measures the bulk modulus of the rock skeleton under different pressure conditions and estimates the effective pressure constant.

[0122] The third calculation module, based on the pore space stiffness theory, calculates the elastic modulus of the fractured dry hot rock skeleton under pressure, based on the obtained elastic modulus of the dry hot rock matrix, the elastic modulus of the dry hot rock skeleton with hard pores, and the effective pressure constant.

[0123] The fourth calculation module calculates the bulk modulus of pore fluids affected by temperature and pressure.

[0124] The fifth calculation module calculates the elastic modulus and density of fluid-bearing dry hot rock under high temperature and high pressure conditions, based on the bulk modulus of pore fluid affected by temperature and pressure and the elastic modulus of fractured dry hot rock skeleton affected by pressure.

[0125] The various embodiments in this specification are described in a progressive manner, with each embodiment focusing on its differences from other embodiments. Similar or identical parts between embodiments can be referred to interchangeably. For the apparatus disclosed in the embodiments, since they correspond to the methods disclosed in the embodiments, the description is relatively simple; relevant parts can be referred to the method section.

[0126] The above description of the disclosed embodiments enables those skilled in the art to make or use the invention. Various modifications to these embodiments will be readily apparent to those skilled in the art, and the general principles defined herein may be implemented in other embodiments without departing from the spirit or scope of the invention. Therefore, the invention is not to be limited to the embodiments shown herein, but is to be accorded the widest scope consistent with the principles and novel features disclosed herein.

Claims

1. A method of seismic petrophysical modeling of deep high temperature high pressure hot dry rock characterized in that, The method comprises the following steps: Step one, determining the mineral composition and content of the hot dry rock matrix, and calculating the elastic modulus of the hot dry rock matrix; Step two, calculating the elastic modulus of the hot dry rock skeleton containing hard pores by using a differential equivalent medium model; Step three, testing the bulk modulus of the rock skeleton under different pressure conditions, and estimating the effective pressure constant; Step four, calculating the elastic modulus of the hot dry rock skeleton containing fractures under the influence of pressure based on the pore space stiffness theory and on the basis of the obtained elastic modulus of the hot dry rock matrix, the elastic modulus of the hot dry rock skeleton containing hard pores, and the effective pressure constant; specifically comprising: regarding the hot dry rock skeleton containing hard pores as a new matrix, adding fractures to the new matrix by using the stiffness pore space model, and obtaining the bulk modulus and shear modulus of the hot dry rock skeleton containing fractures under the influence of pressure: ； ； ； wherein K dry and Step five, calculating the bulk modulus of the pore fluid under the influence of temperature and pressure; specifically comprising: dry are the bulk and shear moduli of the fractured hot dry rock skeleton, respectively, and are the bulk and shear moduli of the hard-porous hot dry rock skeleton, respectively, t and st are the total and hard-pore porosities, respectively, G and M are the effective stress constants, P is the effective stress, V ck is the fracture content; calculating the bulk modulus and density of the pore fluid under the influence of temperature and pressure: wherein, ； ； Step six, calculating the elastic modulus and density of the hot dry rock containing fluid under high temperature and high pressure conditions based on the bulk modulus of the pore fluid under the influence of temperature and pressure and the elastic modulus of the hot dry rock skeleton containing fractures under the influence of pressure; specifically comprising: ； ； ； ； ； wherein G 1 is the air specific weight, T a T is the absolute temperature, P P is the effective pressure, T pr Tc is the critical temperature, P pr Pc is the critical pressure, R R is the gas constant; calculating the bulk modulus and shear modulus of the hot dry rock containing fluid under the preset temperature and pressure conditions by using the Biot-Gassmann equation: calculating the elastic parameters of the hot dry rock containing fluid: ； ； wherein K dry and Step one specifically comprises: dry Ei and Gj are the bulk and shear moduli of the fractured hot dry rock skeleton, respectively, K m Ei is the bulk modulus of the hot dry rock matrix, K sat and determining the mineral composition and percentage content of the hot dry rock by using core X-ray spectroscopy analysis, and calculating the bulk modulus and shear modulus of the hot dry rock matrix by using the Voigt-Reuss-Hill model: sat Ei and Gj are the bulk and shear moduli of the fractured hot dry rock skeleton, respectively, K fl Ei is the bulk modulus of the hot dry rock matrix, calculating the density of the hot dry rock: ； ； ； wherein V P Vp is the dry hot rock P-wave velocity, V S Vs is the dry hot rock S-wave velocity, Step three specifically comprises: p is the dry hot rock density, applying the deep high temperature and high pressure hot dry rock seismic rock physics modeling method according to any one of claims 1-3, which comprises: m and a first calculation module for determining the mineral composition and content of the hot dry rock matrix, and calculating the elastic modulus of the hot dry rock matrix; fl p and p are the densities of the dry hot rock matrix and pore fluid, respectively.

2. The method of claim 1, wherein, a second calculation module for calculating the elastic modulus of the hot dry rock skeleton containing hard pores by using a differential equivalent medium model; a test and estimation module for testing the bulk modulus of the rock skeleton under different pressure conditions, and estimating the effective pressure constant; ； ； a third calculation module for calculating the elastic modulus of the hot dry rock skeleton containing fractures under the influence of pressure based on the pore space stiffness theory and on the basis of the obtained elastic modulus of the hot dry rock matrix, the elastic modulus of the hot dry rock skeleton containing hard pores, and the effective pressure constant; ； In the formula, K m , a fourth calculation module for calculating the bulk modulus of the pore fluid under the influence of temperature and pressure; m and a fifth calculation module for calculating the elastic modulus and density of the hot dry rock containing fluid under high temperature and high pressure conditions based on the bulk modulus of the pore fluid under the influence of temperature and pressure and the elastic modulus of the hot dry rock skeleton containing fractures under the influence of pressure. m These are the bulk modulus, shear modulus, and density of the dry, hot rock matrix, respectively. K mi , ​ mi , ​ mi and V mi The first i The bulk modulus, shear modulus, density, and percentage content of a mineral. N This represents the total number of mineral types.

3. The method of claim 1, wherein, ​ Twenty dry hot rock standard rock samples are prepared, and the total porosity of each rock sample is measured t|k wherein k = 1, 2, 3,..., 20, each rock sample is tested for bulk modulus at effective pressures P = 5, 10, 20, 30, 40, 50 MPa K d|k,P ; From the bulk modulus and porosity measurements of the 20 rock samples, the porosity-stiffness ratio was obtained by fitting the porosity and bulk modulus measurements at each effective pressure condition using the least squares method k 1, the relationship for fitting is as follows: ； The fitting of porosity and bulk modulus for all the measured values under different effective pressures was performed, and the ratio of pore stiffness was obtained k 1 Value of change with effective pressure k p Further, the least square method was used to fit the ratio of pore stiffness and effective pressure, and the effective pressure constant was obtained G and M The relationship for fitting is as follows: 。 4. A deep high temperature high pressure hot dry rock seismic petrophysical modeling system characterized by, ​ ​ ​ ​ ​ ​ ​

Images