阅读说明：本技术 盐胶结砂岩岩石物理模型的建立方法 (Method for establishing salt-cemented sandstone rock physical model ) 是由 韩利 韩文明 阳怀忠 李欣 袁野 刘志国 程涛 黄健良 陈全红 刘新颖 于 2020-07-08 设计创作，主要内容包括：本发明涉及一种盐胶结砂岩岩石物理模型的建立方法,其特征在于,包括如下步骤：1)首先根据测井、壁心分析资料获取岩石各组分参数,并根据岩石物理手册获取各组分纵、横波速度和密度值；然后设置岩石总孔隙度为石盐体积分量与岩石盐胶结后的有效孔隙度之和；2)根据步骤1)获取的参数,依次使用Voigt-Reuss-Hill模型、修改的Xu-White模型、Wood模型和Gassmann方程四个等效模型,建立盐胶结砂岩岩石物理等效模型。本发明能够基于该模型实现流体/固体替代。(The invention relates to a method for establishing a salt cemented sandstone rock physical model, which is characterized by comprising the following steps of: 1) firstly, acquiring parameters of each component of rock according to well logging and wall center analysis data, and acquiring longitudinal and transverse wave velocity and density values of each component according to a rock physics manual; then setting the total porosity of the rock to be the sum of the volume fraction of rock salt and the effective porosity of the cemented rock salt; 2) and (2) according to the parameters obtained in the step 1), sequentially using a Voigt-reus-Hill model, a modified Xu-White model, a Wood model and a Gassmann equation to establish a salt cemented sandstone rock physical equivalent model. The present invention enables fluid/solid substitution based on this model.)

1. The method for establishing the salt cementation sandstone rock physical model is characterized by comprising the following steps of:

1) firstly, acquiring parameters of each component of rock according to well logging and wall center analysis data, and acquiring longitudinal and transverse wave velocity and density values of each component according to a rock physics manual; then setting the total porosity of the rock to be the sum of the volume fraction of rock salt and the effective porosity of the cemented rock salt;

2) and (2) according to the parameters obtained in the step 1), sequentially using a Voigt-reus-Hill model, a modified Xu-White model, a Wood model and a Gassmann equation to establish a salt cemented sandstone rock physical equivalent model.

2. The method for establishing the salt-cemented sandstone petrophysical model according to claim 1, wherein in the step 1), the sum of the rock total porosity, which is the rock salt volume fraction, and the rock salt-cemented effective porosity is:

wherein the content of the first and second substances,the total porosity of the rock before salt cementation; f. of_haliteThe volume fraction of the rock salt is;is the effective porosity after rock salt cementation.

3. The method for establishing the salt consolidated sandstone petrophysical model according to claim 1, wherein the step 2) comprises the following substeps:

calculating the equivalent modulus and density of the rock matrix;

calculating the equivalent modulus and density of the rock skeleton;

thirdly, calculating equivalent modulus and density of the pore mixed fluid;

fourthly, calculating the equivalent modulus and density of the saturated rock;

calculating the longitudinal wave speed, the transverse wave speed and the density of the rock;

when the rock is dry rock, namely the rock pore space is not filled with any fluid, skipping the step (c) and the step (c); and (c) when the rock pore space is filled with fluid, implementing the step (c) and the step (c).

4. The method for building the salt consolidated sandstone petrophysical model of claim 3, wherein in the step ①, the Voigt-reus-Hill model is used for calculating the equivalent modulus of the rock matrix, and in the Voigt-reus-Hill model, the Voigt strain average M is_VIs the upper limit of the equivalent elastic modulus, Reuss stress mean M, of the mixture of N minerals_RIs the lower limit of the equivalent modulus of elasticity, Hill mean M_VRHIs the arithmetic mean of the upper Voigt limit and the lower reus limit, i.e.:

M_VRH＝(M_V+M_R)/2 (4)

wherein f is_i，M_iRespectively, the volume fraction and the elastic modulus of the ith composition component of the rock matrix; n is the type of the composition of the rock matrix;

when M in the formula (2-4) is rock bulk modulus and shear modulus, calculating the equivalent bulk modulus K of the rock matrix by the formula (5) and the formula (6) respectively₀And equivalent shear modulus μ₀：

Wherein, K_i，μ_iThe bulk modulus and shear modulus of the ith component of the rock matrix;

meanwhile, the average density ρ of the rock matrix is calculated by equation (7)₀：

Where ρ is_iIs the density of the ith constituent of the rock matrix; n is the type of rock matrix constituent.

5. The method for establishing the salt bonded sandstone petrophysical model according to claim 4, wherein in the step (II), the modified Xu-White model is used for calculating the equivalent modulus of the rock skeleton, namely the equivalent bulk modulus and the equivalent shear modulus of the rock skeleton are respectively calculated by solving linear ordinary differential equations (8) and (9):

wherein the content of the first and second substances,andrespectively effective porosity ofThe equivalent bulk modulus and the equivalent shear modulus of the rock skeleton; the coefficients p and q are defined as:

wherein the rock skeleton is dry rock, v_sAnd v_cVolume fractions of sandstone and mudstone in the rock skeleton, respectively, sandstone being other rock matrix except clay, mudstone being clay, and v_s＝1-v_cα is the pore aspect ratio, α_sAnd α_cThe aspect ratios of the sandstone and mudstone pores are respectively; t is_iijj(α) and F (α) are functions of pore aspect ratio,

meanwhile, the average density ρ of the rock skeleton is calculated by equation (12)_dry：

6. The method of establishing the salt consolidated sandstone petrophysical model of claim 5, wherein in the step ③, the equivalent modulus and density of the pore-mixing fluid are calculated by using the Wood model, that is, the equivalent bulk modulus K of the pore-mixing fluid is calculated by the formula (13) and the formula (14)_flAnd average density ρ_fl：

7. The method for building the salt consolidated sandstone petrophysical model of claim 6, wherein in the step ④, the Gassmann equation is used to add the pore mixture fluid into the rock skeleton model, and the equivalent bulk modulus K of the saturated rock is calculated by the formula (15) and the formula (16) respectively_satAnd equivalent shear modulus μ_sat：

μ_sat＝μ_dry(16)

Meanwhile, the average density ρ of the saturated rock is calculated by equation (17)_sat：

8. The method for building the physical model of salt-cemented sandstone rock of claim 7, wherein in the step ⑤, when the rock is dry rock, i.e. the rock pores are not filled with any fluid, the fluid comprises oil, gas or water, and the rock longitudinal wave velocity V is calculated by the equations (18), (19) and (20) respectively by using the equivalent bulk modulus, equivalent shear modulus and average density of the rock skeleton output in the step ②_pTransverse wave velocity V_sAnd density ρ:

ρ＝ρ_dry(20)

when the rock pores are filled with fluid, which includes oil, gas or water, the longitudinal wave velocity V of the rock is calculated by equations (21), (22) and (23) using the equivalent bulk modulus, equivalent shear modulus and average density of the saturated rock output in step ④, respectively_pTransverse wave velocity V_sAnd density ρ:

ρ＝ρ_sat(23)。

9. the method of building a salt consolidated sandstone petrophysical model according to any of claims 1 to 8, characterized in that the building method further comprises the following steps:

3) and calibrating and verifying the salt cemented sandstone rock physical equivalent model by using the logging and wall center analysis data, and adjusting and determining the aspect ratio parameter of the pore.

10. The method of constructing a salt-cemented sandstone petrophysical model according to claim 9, further comprising the steps of:

4) and (3) realizing fluid/solid substitution by using the calibrated salt cemented sandstone rock physical equivalent model, and obtaining the elastic parameters of the condition that the fluid is substituted by the rock salt or the fluid is substituted by the rock salt.

Technical Field

The invention relates to the field of oil and gas geophysical exploration, in particular to a method for establishing a salt cementing sandstone rock physical model.

Background

West non salt down-hole drilling revealed that some sandstone was densified due to salt cementation. The salt cementation development range influences the reserve distribution of oil and gas fields and exploration and development strategies. The mural sandstone sheet shows that salt-free cemented sandstone and salt-cemented sandstone are distinguished in that the former pore filling is a fluid (including a gas), while the latter pore filling contains solid rock salt. Different pore fillers and different filler contents give rise to complex seismic response characteristics. Therefore, the development of the seismic response characteristic research of the salt-cemented sandstone and the salt-free cemented hydrocarbon-bearing sandstone in the potential salt-cemented development area is particularly important.

Fluid substitution is a common means for studying the change of seismic response with the change of pore fillers, but the traditional fluid substitution technology based on the Gassmann theory at present requires pore communication, pore pressure balance and zero pore filler shear modulus, and is only suitable for fluid/fluid substitution. And the solid rock salt in the salt cemented sandstone fills pores, so that the pores are not communicated any more, the shear modulus of the pore filler is not zero any more, the pore pressure is unbalanced, and the traditional fluid substitution technology is not suitable for fluid/solid substitution in the salt cemented sandstone seismic response characteristic research.

Ciz and Shapiro (2007) develop on the basis of Gassman equation aiming at the condition that the pore filling is heavy oil, and provide a substitution equation suitable for the pore filling with a non-zero shear modulus, but the equation still assumes that pores are communicated and pore pressure is consistent, and is still not suitable for the condition of salt cemented sandstone, and other scholars (Sun et al, 2018) prove that the equation has errors in fluid/solid substitution prediction through physical experiments and numerical simulation. Before no new effective theory has emerged, using petrophysical modeling is the most feasible way to achieve fluid/solid substitution, but at present there is no petrophysical model for salt-cemented sandstone.

Disclosure of Invention

In view of the above problems, the present invention aims to provide a method for establishing a physical rock model of salt-cemented sandstone, and to enable fluid/solid substitution based on the model.

In order to achieve the purpose, the invention adopts the following technical scheme:

the method for establishing the salt bonded sandstone rock physical model comprises the following steps of:

1) firstly, acquiring parameters of each component of rock according to well logging and wall center analysis data, and acquiring longitudinal and transverse wave velocity and density values of each component according to a rock physics manual; then setting the total porosity of the rock to be the sum of the volume fraction of rock salt and the effective porosity of the cemented rock salt;

2) and (2) according to the parameters obtained in the step 1), sequentially using a Voigt-reus-Hill model, a modified Xu-White model, a Wood model and a Gassmann equation to establish a salt cemented sandstone rock physical equivalent model.

Preferably, in the step 1), the total rock porosity is the sum of the rock salt volume fraction and the effective porosity after rock salt cementation, and is as follows:

wherein the content of the first and second substances,the total porosity of the rock before salt cementation; f. of_haliteThe volume fraction of the rock salt is;is the effective porosity after rock salt cementation.

The method for establishing the salt-cemented sandstone rock physical model preferably comprises the following substeps in the step 2):

calculating the equivalent modulus and density of the rock matrix;

calculating the equivalent modulus and density of the rock skeleton;

thirdly, calculating equivalent modulus and density of the pore mixed fluid;

fourthly, calculating the equivalent modulus and density of the saturated rock;

calculating the longitudinal wave speed, the transverse wave speed and the density of the rock;

when the rock is dry rock, namely the rock pore space is not filled with any fluid, skipping the step (c) and the step (c); and (c) when the rock pore space is filled with fluid, implementing the step (c) and the step (c).

Preferably, in the step ①, the method for establishing the salt-cemented sandstone rock physical model calculates the equivalent modulus of the rock matrix by using a Voigt-reus-Hill model, and in the Voigt-reus-Hill model, the Voigt strain average M is_VIs the upper limit of the equivalent elastic modulus, Reuss stress mean M, of the mixture of N minerals_RIs the lower limit of the equivalent modulus of elasticity, Hill mean M_VRHIs the arithmetic mean of the upper Voigt limit and the lower reus limit, i.e.:

M_VRH＝(M_V+M_R)/2 (4)

wherein f is_i，M_iRespectively, the volume fraction and the elastic modulus of the ith composition component of the rock matrix; n is the type of rock matrix constituent.

When M in the formula (2-4) is rock bulk modulus and shear modulus, respectively represented by the formula (5) and the formula (6)

Calculating the equivalent bulk modulus K of the rock matrix₀And equivalent shear modulus μ₀：

Wherein K, mu are the bulk modulus and shear modulus of the ith component of the rock matrix respectively;

meanwhile, the average density ρ of the rock matrix is calculated by equation (7)₀：

Where ρ is_iIs the density of the ith constituent of the rock matrix; n is the type of rock matrix constituent.

Preferably, in the second step, the modified Xu-White model is used to calculate the equivalent modulus of the rock skeleton, that is, the equivalent bulk modulus and the equivalent shear modulus of the rock skeleton are calculated by solving linear ordinary differential equations (8) and (9):

wherein the content of the first and second substances,and

respectively effective porosity ofEquivalent bulk modulus and equivalent shear modulus of the rock skeleton (dry rock) in time; the coefficients p and q are defined as:

wherein v is_sAnd v_cVolume fractions of sandstone (other rock matrix than clay) and mudstone (clay) in the rock skeleton, respectively, and v_s＝1-v_cα is the pore aspect ratio, α_sAnd α_cThe aspect ratios of the sandstone and mudstone pores are respectively; t is_iijj(α) and F (α) are functions of pore aspect ratio,

meanwhile, the average density ρ of the rock skeleton is calculated by equation (12)_dry：

Preferably, in the step ③, the Wood model is used to calculate the equivalent modulus and density of the pore mixture fluid, that is, the equivalent bulk modulus K of the pore mixture fluid is calculated by the equations (13) and (14) respectively_flAnd average density ρ_fl：

The method for establishing the physical model of the salt-cemented sandstone rock preferably utilizes the method in step ④Adding pore mixed fluid into a rock skeleton model by a Gassmann equation, and respectively calculating the equivalent bulk modulus K of the saturated rock by the formulas (15) and (16)_satAnd equivalent shear modulus μ_sat：

μ_sat＝μ_dry(16)

Meanwhile, the average density ρ of the saturated rock is calculated by equation (17)_sat：

Preferably, in the fifth step, when the rock is dry rock, that is, when the rock pores are not filled with any fluid (including oil, gas or water), the second step is used to establish the physical model of the salt-cemented sandstone rock

The equivalent bulk modulus, the equivalent shear modulus and the average density of the rock skeleton (dry rock) are respectively output, and the longitudinal wave velocity V of the rock is calculated by the formulas (18), (19) and (20)_pTransverse wave velocity V_sAnd density ρ:

ρ＝ρ_dry(20)

when the rock pores are filled with fluid (including oil, gas or water), the equivalent bulk modulus, the equivalent shear modulus and the average density of the saturated rock output in the step (r) are respectively expressed by the formulas (21), (22) and (23)

Calculating the longitudinal wave velocity V of the rock_pTransverse wave velocity V_sAnd density ρ:

ρ＝ρ_sat(23)

the method for establishing the salt cementation sandstone rock physical model preferably further comprises the following steps:

3) and calibrating and verifying the salt cemented sandstone rock physical equivalent model by using the logging and wall center analysis data, and adjusting and determining the aspect ratio parameter of the pore.

The method for establishing the salt cementation sandstone rock physical model preferably further comprises the following steps:

4) and (3) realizing fluid/solid substitution by using the calibrated salt cemented sandstone rock physical equivalent model, and obtaining the elastic parameters of the condition that the fluid is substituted by the rock salt or the fluid is substituted by the rock salt.

Due to the adoption of the technical scheme, the invention has the following advantages: 1. the salt-cemented petrophysical models of the present invention are applicable not only to conventional fluid/fluid substitution, but also to fluid/solid substitution. 2. In the present invention, the halite component serves both as a pore filler to occupy the total pore space prior to salt cementation and as a rock matrix to support the elastic modulus of the rock. The salt-cemented sandstone rock physical model is suitable for the condition of filling a pore with solid or filling a pore with solid-liquid mixture, and the model has high coincidence degree of the calculation result and the actual logging result and is reliable.

Drawings

FIG. 1 is a schematic illustration of a wall core sample sandstone rock slice; wherein figure 1(a) shows a salt-free cemented sandstone sheet and figure 1(b) shows a salt-cemented sandstone sheet;

FIG. 2 is a schematic of the construction of sandstone rock; wherein fig. 2(a) is a schematic view of the framework of salt-free cemented sandstone rock, and fig. 2(b) is a schematic view of the framework of localized salt-cemented sandstone rock;

FIG. 3 is a schematic illustration of the composition of a salt cemented sandstone rock component;

FIG. 4 is a flow chart of a method for establishing a physical model of salt cemented sandstone rock;

FIG. 5 is a flow chart of a method for establishing a physical equivalent model of salt cemented sandstone rock;

FIG. 6 is a verification diagram of a salt cemented sandstone petrophysical model;

FIG. 7 is a graph of elastic response and amplitude response for different rock salt content models; wherein, FIG. 7(a) is a geological model map of sandstone with different salt content, FIGS. 7(b-d) are velocity, density and compressional impedance curves, respectively, and FIGS. 7(e-f) are seismic record and maximum amplitude curve, respectively;

FIG. 8 is a schematic diagram of the relationship between salt content, total porosity and longitudinal wave impedance; wherein, fig. 8(a) is a relation model diagram of total porosity and longitudinal wave impedance under the condition of different salt contents, and fig. 8(b) is a relation model diagram of salt content and longitudinal wave impedance under the condition of different total porosity;

FIG. 9 is a graph of fluid/solid (halite) alternative seismic amplitude response; wherein, fig. 9(a) is a geological model diagram, fig. 9(b-e) are a salt-cemented dry sandstone seismic record diagram, a petrophysical model simulation salt-cemented dry sandstone seismic record diagram, a petrophysical model seismic record diagram of replacing rock salt with water condition and a petrophysical model seismic record diagram of replacing rock salt with gas condition respectively according to in-situ logging parameters;

figure 10 is a graph of AVO response for this zone based on fluid/solid (halite) substitution.

Detailed Description

The preferred embodiments of the present invention will be described in detail below with reference to the accompanying drawings so that the objects, features and advantages of the invention can be more clearly understood. It should be understood that the embodiments shown in the drawings are not intended to limit the scope of the present invention, but are merely intended to illustrate the spirit of the technical solution of the present invention.

Before describing the technical solution of the present invention, it is necessary to first make a brief description of the characteristics of salt cemented sandstone in order that those skilled in the art can better understand the essential spirit of the technical solution of the present invention. Figure 1(a) shows a salt-free cemented sandstone sheet with intergranular pore development, pore communication, fluid (including oil, gas, water) fill, pore communication, pore pressure equalization, and zero pore fill shear modulus; while figure 1(b) shows a salt-cemented sandstone sheet in which the intergranular pores are filled with solid rock salt, the pores are no longer connected and the pore filling shear modulus is no longer zero.

As shown in fig. 2, salt free cemented sandstone and salt cemented sandstone rock frameworks were constructed from the rock slices. Sandstone mainly comprises solid parts such as rock minerals (quartz, feldspar and the like), rock debris, cement and the like, and pores and pore fillers. The non-salt cemented sandstone and the salt cemented sandstone are different in that the filler in the pores of the former is fluid (including gas) (see fig. 1(a) and 2(a)), and the filler in the pores of the latter contains solid rock salt (see fig. 1(b) and 2 (b)).

As shown in fig. 3, according to the wall core analysis, the sandstone rock before salt cementation is mainly composed of quartz, feldspar, rock debris, cement, pore fluid and other components, wherein the cement comprises quartz cementation (quartz), carbonate cementation (dolomite), clay cementation (illite, chlorite and kaolinite), and the pore fluid is oil, gas or water. After salt cementation, most of the pore space is filled with the stone salt, so that the effective porosity is reduced. The bonding of the rock salt as a solid with the mineral particles as part of the rock matrix not only affects the bulk modulus of the rock but also increases the rock shear modulus.

As shown in fig. 4, the method for establishing the physical rock model of the salt cemented sandstone, provided by the invention, aiming at the characteristics of the salt cemented sandstone, comprises the following steps:

the first step is as follows: firstly, obtaining parameters of each component of the rock through well logging and wall center data, wherein the parameters comprise integral quantity, speed, density and the like of each component of the rock, and obtaining longitudinal and transverse wave speed and density values of each component according to a rock physics manual (Mavko et al, 1998); then setting the total rock porosity as the sum of the rock salt volume fraction and the effective porosity after rock salt cementation, namely:

wherein the content of the first and second substances,

the total porosity of the rock before salt cementation; f. of_haliteThe volume fraction of the rock salt is;is the effective porosity after rock salt cementation.

The second step is that: as shown in fig. 5, according to the parameters obtained in the first step, four equivalent models, namely a Voigt-Reuss-Hill model, a modified Xu-White model, a Wood model and a Gassmann equation, are sequentially used to establish a physical equivalent model of the salt cemented sandstone rock, which specifically comprises the following steps:

1) calculating the equivalent modulus and density of the rock matrix:

in the invention, the equivalent modulus of the rock matrix is calculated by utilizing a Voigt-reus-Hill model, and in the Voigt-reus-Hill model, the Voigt strain is averaged to be M_VIs the upper limit of the equivalent elastic modulus, Reuss stress mean M, of the mixture of N minerals_RIs the lower limit of the equivalent modulus of elasticity, Hill mean M_VRHIs the arithmetic mean of the upper Voigt limit and the lower reus limit, i.e.:

M_VRH＝(M_V+M_R)/2 (4)

wherein f is_i，M_iRespectively, the volume fraction and the elastic modulus of the ith composition component of the rock matrix; n is the type of rock matrix constituent.

When M in the formula (2-4) is rock bulk modulus and shear modulus, calculating the equivalent bulk modulus K of the rock matrix by the formula (5) and the formula (6) respectively₀And equivalent shear modulus μ₀：

Wherein, K_i，μ_iRespectively, the bulk modulus and shear modulus of the ith constituent component of the rock matrix. The rock matrix comprises quartz, feldspar, rock debris, conventional cement, rock salt and other solid components.

Meanwhile, the average density ρ of the rock matrix is calculated by equation (7)₀：

Where ρ is_iIs the density of the ith constituent of the rock matrix; n is the type of rock matrix constituent.

2) Calculating equivalent modulus and density of the rock skeleton:

in the invention, the equivalent modulus of the rock skeleton is calculated by using the modified Xu-White model, namely the equivalent bulk modulus and the equivalent shear modulus of the rock skeleton are respectively calculated by solving linear ordinary differential equations (8) and (9):

wherein the content of the first and second substances,andrespectively effective porosity ofEquivalent bulk modulus and equivalence of rock skeleton (dry rock)A shear modulus; the coefficients p and q are defined as:

wherein v is_sAnd v_cVolume fractions of sandstone (other rock matrix than clay) and mudstone (clay) in the rock framework, respectively, and v_s＝1-v_cα is the pore aspect ratio, α_sAnd α_cThe aspect ratios of the sandstone and mudstone pores are respectively; t is_iijj(α) and F (α) are functions of pore aspect ratio, as defined by Kuster and(1974) and Xu and White (1995).

Meanwhile, the average density ρ of the rock skeleton is calculated by equation (12)_dry：

3) Calculating equivalent modulus and density of pore mixed fluid:

in the invention, the equivalent modulus and density of the pore mixed fluid are calculated by using a Wood model, namely the equivalent bulk modulus (Reuss average) K of the pore mixed fluid is calculated by respectively using an equation (13) and an equation (14)_flAnd average density ρ_fl：

4) Calculating equivalent modulus and density of saturated rock:

in the present invention, the Gassmann equation is utilized toAdding the pore mixed fluid into a rock skeleton model, and respectively calculating the equivalent bulk modulus K of the saturated rock by using a formula (15) and a formula (16)_satAnd equivalent shear modulus μ_sat：

μ_sat＝μ_dry(16)

Meanwhile, the average density ρ of the saturated rock is calculated by equation (17)_sat：

5) Calculating the longitudinal wave velocity, the transverse wave velocity and the density of the rock:

in the present invention, when the rock is dry rock, that is, when the rock pores are not filled with any fluid (including oil, gas or water), the longitudinal wave velocity V of the rock is calculated by equations (18), (19) and (20) using the equivalent bulk modulus, equivalent shear modulus and average density of the rock skeleton (dry rock) output in said step 2), respectively_pTransverse wave velocity V_sAnd density ρ:

ρ＝ρ_dry(20)

in the present invention, when the rock pores are filled with fluid (including oil, gas or water), said step 4) is utilized

The equivalent volume modulus, equivalent shear modulus and average density of the saturated rock output in the process are respectively represented by the formula (21),

(22) And (23) calculating the longitudinal wave velocity V of the rock_pTransverse wave velocity V_sAnd density ρ:

ρ＝ρ_sat(23)

in the above embodiment, when the rock is dry rock, i.e. the rock pores are not filled with any fluid (including oil, gas or water), step 3) and step 4) are skipped; when the rock pores are filled with fluid, step 3) and step 4) are performed.

In the above embodiment, preferably, a third step may be further included: and (3) calibrating and verifying the salt bonded sandstone rock physical equivalent model by using the drilled well measured data (comprising the well logging data and the wall center analysis data). The parameters for determining the aspect ratio of the pores can be adjusted, and specifically comprise: the rock pore aspect ratio parameter is a parameter of the pore morphology of the reaction salt cemented sandstone after compaction, influences the calculation of the rock skeleton equivalent model, and is expressed as a value range which can influence the calculation of the elastic parameter of the model integrally. The volume content of each component of the wall center analysis data is used as model input, the elastic parameters of the salt bonded sandstone rock simulated by the model are compared with well measurement data, and the aspect ratio parameters of the salt bonded sandstone rock pore can be calibrated. Within a zone of interest, the parameters are usually fixed because diagenesis is similar. After the parameters are determined, the correlation degree of simulation results of all sample models and well measured data is calculated, and when the correlation degree is higher, the model is proved to be effective, so that the sandstone elastic parameters in the area can be well simulated.

In the above embodiment, preferably, the method may further include the fourth step of: the model is used for fluid/solid (halite) replacement, and the elastic parameters of the halite replacement fluid or the fluid replacement halite are obtained to research the elastic parameters and the seismic amplitude response characteristics of the salt cemented sandstone rock under the conditions of different pore filling types and different halite contents.

The physical modeling of the salt bonded sandstone rock and the fluid/solid substitution based on the model of the invention are further described below by way of example.