CN113552626B - Method and device for predicting seismic transverse wave velocity of tight sandstone reservoir - Google Patents

Method and device for predicting seismic transverse wave velocity of tight sandstone reservoir Download PDF

Info

Publication number
CN113552626B
CN113552626B CN202110754161.XA CN202110754161A CN113552626B CN 113552626 B CN113552626 B CN 113552626B CN 202110754161 A CN202110754161 A CN 202110754161A CN 113552626 B CN113552626 B CN 113552626B
Authority
CN
China
Prior art keywords
wave velocity
modulus
seismic
elastic modulus
rock
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202110754161.XA
Other languages
Chinese (zh)
Other versions
CN113552626A (en
Inventor
刘俊州
刘致水
包乾宗
时磊
王震宇
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China Petroleum and Chemical Corp
Changan University
Original Assignee
China Petroleum and Chemical Corp
Changan University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by China Petroleum and Chemical Corp, Changan University filed Critical China Petroleum and Chemical Corp
Priority to CN202110754161.XA priority Critical patent/CN113552626B/en
Publication of CN113552626A publication Critical patent/CN113552626A/en
Application granted granted Critical
Publication of CN113552626B publication Critical patent/CN113552626B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28—Processing seismic data, e.g. analysis, for interpretation, for correction
    • G01V1/30—Analysis
    • G01V1/303—Analysis for determining velocity profiles or travel times
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28—Processing seismic data, e.g. analysis, for interpretation, for correction
    • G01V1/30—Analysis
    • G01V1/306—Analysis for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00—Details of seismic processing or analysis
    • G01V2210/60—Analysis
    • G01V2210/62—Physical property of subsurface
    • G01V2210/622—Velocity, density or impedance
    • G01V2210/6222—Velocity; travel time

Abstract

The invention discloses a method and a device for predicting seismic transverse wave velocity of a tight sandstone reservoir, wherein the method comprises the following steps: establishing a differential effective medium rock physical model and a porous rock physical model; obtaining a second elastic modulus according to the differential effective medium rock physical model and the first elastic modulus of the rock matrix; obtaining a third elastic modulus and a fourth elastic modulus according to the porous rock physical model and the first elastic modulus; obtaining the fluid volume modulus of the fluid, and respectively obtaining the theoretical seismic longitudinal wave velocity and the theoretical seismic transverse wave velocity of the saturated fluid rock according to the second elastic module, the third elastic modulus, the fourth elastic modulus and the fluid volume module; determining a first pore volume ratio and a second pore volume ratio according to the seismic longitudinal wave velocity; and predicting the seismic transverse wave velocity according to the first pore volume ratio, the second pore volume ratio and the theoretical seismic transverse wave velocity. The invention improves the accuracy of the predicted seismic transverse wave velocity.

Description

Method and device for predicting seismic transverse wave velocity of tight sandstone reservoir
Technical Field
The invention relates to the technical field of oil and gas geophysical, in particular to a method and a device for predicting seismic transverse wave velocity of a compact sandstone reservoir.
Background
The seismic transverse wave velocity is important data for pre-stack elastic inversion and rock physical analysis, and the calculation of the seismic transverse wave velocity through a rock physical model is an important method for acquiring transverse wave data. The accurate rock physical model can effectively describe the relationship between the rock microscopic characteristics and the elastic modulus, so that the accurate seismic transverse wave velocity is obtained. Different from the conventional sandstone, the compact sandstone undergoes strong mechanical compaction and cementation in the burying process, and the reservoir has the characteristics of low porosity and low permeability. Due to the low porosity of tight sandstone, the pore structure is the main influence factor of the elastic modulus of rock.
At present, scholars apply elliptical pore models with good application effect in conventional sandstone and carbonate rock to tight sandstone reservoirs, the models enable pores in the tight sandstone to be equivalent to 2D ellipses, and the difference between the 2D ellipses and polygonal pores in real rocks is large, so that the problem that the seismic transverse wave velocity prediction of the tight sandstone reservoirs is inaccurate is caused.
Disclosure of Invention
The invention aims to overcome the technical defects and provides a method and a device for predicting the seismic transverse wave velocity of a compact sandstone reservoir, aiming at solving the technical problem that the porous rock physical model in the prior art enables the pores in the compact sandstone to be equivalent to 2D ellipses, and the difference between the pores and the polygonal pores of real rocks is larger, so that the seismic transverse wave velocity predicted by the porous rock physical model is inaccurate.
In one aspect, the invention provides a method for predicting seismic shear wave velocity of a tight sandstone reservoir, which comprises the following steps:
establishing a differential effective medium rock physical model and a porous rock physical model;
obtaining a second elastic modulus of the dry rock containing the cracks according to the differential effective medium rock physical model and the first elastic modulus of the rock matrix;
obtaining a third elastic modulus of the convex-edge-containing triangular dry rock and a fourth elastic modulus of the concave-edge-containing quadrilateral dry rock according to the porous rock physical model and the first elastic modulus;
obtaining the fluid volume modulus of the fluid, and respectively obtaining a first seismic longitudinal wave velocity, a first seismic transverse wave velocity, a second seismic longitudinal wave velocity, a second seismic transverse wave velocity, a third seismic longitudinal wave velocity and a third seismic transverse wave velocity of the saturated fluid rock according to the second elastic module, the third elastic modulus, the fourth elastic modulus and the fluid volume module;
determining a first pore volume ratio and a second pore volume ratio according to the first seismic longitudinal wave velocity, the second seismic longitudinal wave velocity and the third seismic longitudinal wave velocity;
and predicting the seismic transverse wave velocity according to the first pore volume ratio, the second pore volume ratio, the first seismic transverse wave velocity, the second seismic transverse wave velocity and the third seismic transverse wave velocity.
In one possible implementation manner of the present invention, the first elastic modulus includes a first bulk modulus and a first shear modulus, and the second elastic modulus includes a second bulk modulus and a second shear modulus; the differential effective medium rock physical model is as follows:
Figure RE-GDA0003221563280000031
Figure RE-GDA0003221563280000032
Figure RE-GDA0003221563280000033
Figure RE-GDA0003221563280000034
Figure RE-GDA0003221563280000035
in the formula, K * DEM (y) is the second bulk modulus; g * DEM (y) is the second shear modulus; k is i Is the first bulk modulus; g i Is the first shear modulus; y is the volume fraction of cracks in the dry rock; p *i (y)、Q *i (y) is an impact factor associated with a fracture in dry rock; alpha is the first equivalent shape factor.
In a possible implementation manner of the present invention, the third elastic modulus includes a third bulk modulus and a third shear modulus, and the fourth elastic modulus includes a fourth bulk modulus and a fourth shear modulus; the physical model of the porous rock is as follows:
Figure RE-GDA0003221563280000036
Figure RE-GDA0003221563280000037
wherein g is a second equivalent shape factor;
Figure RE-GDA0003221563280000038
is porosity; k is dry Is the third bulk modulus or the fourth bulk modulus; g dry Is the third shear modulus or the fourth shear modulus.
In a possible implementation manner of the present invention, the obtaining, according to the second elastic module, the third elastic modulus, the fourth elastic modulus, and the fluid volume module, a first seismic longitudinal wave velocity, a first seismic transverse wave velocity, a second seismic longitudinal wave velocity, a second seismic transverse wave velocity, a third seismic longitudinal wave velocity, and a third seismic transverse wave velocity of the saturated fluid rock respectively includes:
obtaining a fifth elastic modulus of the saturated fluid rock from the second elastic module and the fluid volume module;
obtaining a sixth elastic modulus of the saturated fluid rock according to the third elastic module and the fluid volume module;
obtaining a seventh elastic modulus of the saturated fluid rock according to the fourth elastic module and the fluid volume module;
and sequentially obtaining the first seismic longitudinal wave velocity, the first seismic transverse wave velocity, the second seismic longitudinal wave velocity, the second seismic transverse wave velocity, the third seismic longitudinal wave velocity and the third seismic transverse wave velocity according to the fifth elastic modulus, the sixth elastic modulus and the seventh elastic modulus.
In a possible implementation manner of the present invention, the fifth elastic modulus includes a fifth bulk modulus and a fifth shear modulus, and the fifth elastic modulus is:
Figure RE-GDA0003221563280000041
Figure RE-GDA0003221563280000042
in the formula, K sat1 Is the fourth bulk modulus; g sat1 Is the fourth shear modulus;
the sixth elastic modulus comprises a sixth bulk modulus and a sixth shear modulus, and the seventh elastic modulus comprises a seventh bulk modulus and a seventh shear modulus; the sixth elastic modulus is:
Figure RE-GDA0003221563280000051
G sat2 ＝G dry1
the seventh elastic modulus is:
Figure RE-GDA0003221563280000052
G sat3 ＝G dry2
in the formula, K sat2 Is the sixth bulk modulus; k sat3 Is the seventh bulk modulus; g sat2 Is the sixth shear modulus; g sat3 Is the seventh bulk modulus; k is dry1 Is the third bulk modulus; g dry1 Is the third shear modulus; k is dry2 Is the fourth bulk modulus; g dry2 Is the fourth shear modulus.
In a possible implementation manner of the present invention, the first seismic longitudinal wave velocity, the first seismic transverse wave velocity, the second seismic longitudinal wave velocity, the second seismic transverse wave velocity, the third seismic longitudinal wave velocity, and the third seismic transverse wave velocity are respectively:
Figure RE-GDA0003221563280000053
Figure RE-GDA0003221563280000054
Figure RE-GDA0003221563280000055
Figure RE-GDA0003221563280000061
Figure RE-GDA0003221563280000062
Figure RE-GDA0003221563280000063
in the formula, V pc_crack Is the first seismic longitudinal wave velocity; v pc_convex The second seismic longitudinal wave velocity; V pc_concave is the third seismic longitudinal wave velocity; v sc_crack Is the first seismic shear velocity; v sc_convex The second seismic transverse wave velocity; v sc_concave Is the third seismic shear velocity; ρ is a unit of a gradient sat Is the density of the saturated fluid rock.
In one possible implementation manner of the present invention, the first pore volume ratio is:
Figure RE-GDA0003221563280000064
the second pore volume ratio is:
Figure RE-GDA0003221563280000065
in the formula, w 1 Is the first pore volume ratio; w is a 2 Is the second pore volume ratio; v pm The measured seismic longitudinal wave velocity is obtained.
In a possible implementation manner of the present invention, the seismic shear wave velocity is:
V sc ＝V sc_convex -w 1 (V sc_convex -V sc_concave )
V sc ＝V sc_concave -w 1 (V sc_concave -V sc_crack )
in the formula, V sc Is the seismic shear wave velocity.
In a possible implementation manner of the present invention, before the obtaining a second elastic modulus of the dry rock containing the fracture according to the differential effective medium rock physical model and the first elastic modulus of the rock matrix, the method further includes:
obtaining a first elastic modulus of the rock matrix;
wherein said obtaining a first elastic modulus of the rock matrix comprises:
determining whether there is more than one mineral species in the rock matrix;
if not, acquiring the first elastic modulus according to a preset parameter table;
if so, acquiring the initial elastic modulus of each mineral in the rock matrix according to a preset parameter table, and acquiring the first elastic modulus according to a Voigt-Reuss-Hill formula and the initial elastic modulus of each mineral.
In another aspect, the present invention provides a device for predicting seismic shear wave velocity of a tight sandstone reservoir, comprising:
the model establishing module is used for establishing a differential effective medium rock physical model and a porous rock physical model;
the first model calculation module is used for obtaining a second elastic modulus of the crack-containing dry rock according to the differential effective medium rock physical model and the first elastic modulus of the rock matrix;
the second model calculation module is used for obtaining a third elastic modulus of the convex-edge-containing triangular dry rock and a fourth elastic modulus of the concave-edge-containing quadrilateral dry rock according to the porous rock physical model and the first elastic modulus of the rock matrix;
the theoretical velocity calculation module is used for obtaining the fluid volume modulus of the fluid and respectively obtaining a first seismic longitudinal wave velocity, a first seismic transverse wave velocity, a second seismic longitudinal wave velocity, a second seismic transverse wave velocity, a third seismic longitudinal wave velocity and a third seismic transverse wave velocity of the saturated fluid rock according to the second elastic module, the third elastic modulus, the fourth elastic modulus and the fluid volume module;
the parameter determining module is used for determining a first pore volume ratio and a second pore volume ratio according to the first seismic longitudinal wave velocity, the second seismic longitudinal wave velocity and the third seismic longitudinal wave velocity;
and the seismic transverse wave velocity prediction module is used for predicting the seismic transverse wave velocity according to the first pore volume ratio, the second pore volume ratio, the first seismic transverse wave velocity, the second seismic transverse wave velocity and the third seismic transverse wave velocity.
In another aspect, the present invention also provides a computer device, including:
one or more processors;
a memory; and
one or more applications, wherein the one or more applications are stored in the memory and configured to be executed by the processor to implement the method of predicting seismic shear wave velocity for tight sand reservoirs of any of the above.
In another aspect, the present invention also provides a computer readable storage medium, on which a computer program is stored, the computer program being loaded by a processor to execute the steps in the method for predicting seismic shear wave velocity of tight sandstone reservoir described in any one of the above.
According to the method, firstly, a differential effective medium rock physical model and a porous rock physical model are constructed, then, a second elastic modulus of a crack-containing dry rock, a third elastic modulus of a convex-edge-triangle-containing dry rock and a fourth elastic modulus of a concave-edge-quadrilateral-containing dry rock are obtained by utilizing the two rock physical models and the first elastic modulus of a rock matrix, then, a first pore volume ratio for representing the pore volume ratio of the convex edge triangle and the concave edge quadrilateral and a second pore volume ratio for representing the pore volume ratio of the concave edge quadrilateral to the crack pore volume ratio are determined according to the third elastic modulus and the fourth elastic modulus, compared with the prior art, the method can describe the pore shape of a sandstone compact reservoir more truly, and finally, the seismic transverse wave velocity is predicted through the first pore volume ratio and the second pore volume ratio, so that the accuracy of the predicted seismic transverse wave velocity is improved.
Drawings
In order to more clearly illustrate the technical solutions in the embodiments of the present invention, the drawings needed to be used in the description of the embodiments will be briefly introduced below, and it is obvious that the drawings in the following description are only some embodiments of the present invention, and it is obvious for those skilled in the art to obtain other drawings based on these drawings without creative efforts.
FIG. 1 is a schematic flow chart of an embodiment of a method for predicting seismic shear wave velocity of a tight sandstone reservoir, provided by an embodiment of the invention;
fig. 2 is a schematic structural diagram of an embodiment of S104 according to the present invention;
FIG. 3 is a schematic flow chart of one embodiment of obtaining a first bulk modulus and a first shear modulus provided by an embodiment of the present invention;
fig. 4 is a schematic structural diagram of an embodiment of a device for predicting seismic shear wave velocity of a tight sandstone reservoir, provided by an embodiment of the present invention;
fig. 5 is a schematic structural diagram of an embodiment of a computer device provided in an embodiment of the present invention.
Detailed Description
The technical solution in the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings in the embodiments of the present invention. It is to be understood that the described embodiments are merely exemplary of the invention, and not restrictive of the full scope of the invention. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
Reference herein to "an embodiment" means that a particular feature, structure, or characteristic described in connection with the embodiment can be included in at least one embodiment of the invention. The appearances of the phrase in various places in the specification are not necessarily all referring to the same embodiment, nor are separate or alternative embodiments mutually exclusive of other embodiments. It is explicitly and implicitly understood by one skilled in the art that the embodiments described herein may be combined with other embodiments.
The invention provides a method and a device for predicting seismic transverse wave velocity of a compact sandstone reservoir, which are respectively explained in detail below.
As shown in fig. 1, a schematic flow chart of an embodiment of a method for predicting seismic shear wave velocity of a tight sandstone reservoir, provided by the embodiment of the present invention, includes:
s101, establishing a differential effective medium rock physical model and a porous rock physical model;
s102, obtaining a second elastic modulus of the dry rock containing the cracks according to the differential effective medium rock physical model and the first elastic modulus of the rock matrix;
s103, obtaining a third elastic modulus of the convex-edge-containing triangular dry rock and a fourth elastic modulus of the concave-edge-containing quadrilateral dry rock according to the porous rock physical model and the first elastic modulus;
s104, obtaining the fluid volume modulus of the fluid, and respectively obtaining a first seismic longitudinal wave velocity, a first seismic transverse wave velocity, a second seismic longitudinal wave velocity, a second seismic transverse wave velocity, a third seismic longitudinal wave velocity and a third seismic transverse wave velocity of the saturated fluid rock according to the second elastic module, the third elastic modulus, the fourth elastic modulus and the fluid volume module;
s105, determining a first pore volume ratio and a second pore volume ratio according to the first seismic longitudinal wave velocity, the second seismic longitudinal wave velocity and the third seismic longitudinal wave velocity;
s106, predicting the seismic transverse wave velocity according to the first pore volume ratio, the second pore volume ratio, the first seismic transverse wave velocity, the second seismic transverse wave velocity and the third seismic transverse wave velocity.
According to the method for predicting the seismic shear wave velocity of the tight sandstone reservoir, firstly, a differential effective medium rock physical model and a porous rock physical model are built, then, a second elastic modulus of a dry rock containing cracks, a third elastic modulus of a dry rock containing a convex triangle and a fourth elastic modulus of a dry rock containing a concave quadrilateral are obtained by utilizing the two rock physical models and a first elastic modulus of a rock matrix, then, a first pore volume ratio for representing the pore volume ratio of the convex triangle and the concave quadrilateral and a second pore volume ratio for representing the pore volume ratio of the concave quadrilateral and the crack are determined according to the third elastic modulus and the fourth elastic modulus, compared with the prior art, the method for predicting the seismic shear wave velocity of the tight sandstone reservoir can describe the pore shape of the tight sandstone reservoir more truly, and finally, the accuracy of the predicted seismic shear wave velocity is improved through the first pore volume ratio and the second pore volume ratio.
It should be noted that: the order of S102 and S103 is not limited, and S103 may be performed simultaneously, or S103 may be performed first and then S102.
Specifically, the first elastic modulus includes a first bulk modulus and a first shear modulus, and the second elastic modulus includes a second bulk modulus and a second shear modulus; the Differential Effective Medium (DEM) rock physical model is:
Figure RE-GDA0003221563280000121
Figure RE-GDA0003221563280000122
Figure RE-GDA0003221563280000123
Figure RE-GDA0003221563280000124
Figure RE-GDA0003221563280000125
in the formula, K * DEM (y) is the second bulk modulus; g * DEM (y) is the second shear modulus; k i Is the first bulk modulus; g i A first shear modulus; y is the volume fraction of cracks in the dry rock; p is *i (y)、Q *i (y) is an impact factor associated with a fracture in dry rock; alpha is the first equivalent shape factor.
It is noted that, in some embodiments of the present invention, α is 0.01.
Further, the third elastic modulus comprises a third bulk modulus and a third shear modulus, and the fourth elastic modulus comprises a fourth bulk modulus and a fourth shear modulus; the physical model of the porous rock is as follows:
Figure RE-GDA0003221563280000131
Figure RE-GDA0003221563280000132
wherein g is a second equivalent shape factor;
Figure RE-GDA0003221563280000133
is porosity; k dry Is a third bulk modulus or the fourth bulk modulus; g dry Either the third shear modulus or the fourth shear modulus.
It should be noted that, in some embodiments of the present invention, g is 3.142 when solving for the third elastic modulus, and g is 7.75 when solving for the fourth elastic modulus.
Further, in some embodiments of the present invention, as shown in fig. 2, S104 includes:
s201, obtaining a fifth elastic modulus of the saturated fluid rock according to the second elastic module and the fluid volume module;
s202, obtaining a sixth elastic modulus of the saturated fluid rock according to the third elastic module and the fluid volume module;
s203, obtaining a seventh elastic modulus of the saturated fluid rock according to the fourth elastic module and the fluid volume module;
and S204, sequentially obtaining a first seismic longitudinal wave velocity, a first seismic transverse wave velocity, a second seismic longitudinal wave velocity, a second seismic transverse wave velocity, a third seismic longitudinal wave velocity and a third seismic transverse wave velocity according to the fifth elastic modulus, the sixth elastic modulus and the seventh elastic modulus.
Through the steps, the first seismic longitudinal wave velocity and the first seismic transverse wave velocity of the saturated fluid rock containing the fracture can be obtained; the second seismic longitudinal wave velocity and the second seismic transverse wave velocity of the saturated fluid rock containing the convex edge triangle; compared with the prior art that the pore in the saturated fluid rock is equivalent to an ellipse, the third seismic longitudinal wave velocity and the third seismic transverse wave velocity of the saturated fluid rock containing the concave side quadrangle are closer to the actual pore, and the reliability and the accuracy of the calculated first seismic longitudinal wave velocity, the calculated first seismic transverse wave velocity, the calculated second seismic longitudinal wave velocity, the calculated third seismic transverse wave velocity and the calculated third seismic transverse wave velocity are improved.
Further, the fifth elastic modulus includes a fifth bulk modulus and a fifth shear modulus, and the fifth elastic modulus is:
Figure RE-GDA0003221563280000141
Figure RE-GDA0003221563280000142
in the formula, K sat1 A fourth bulk modulus; g sat1 Is the fourth shear modulus;
the sixth modulus of elasticity comprises a sixth bulk modulus and a sixth shear modulus, and the seventh modulus of elasticity comprises a seventh bulk modulus and a seventh shear modulus; the sixth modulus of elasticity is:
Figure RE-GDA0003221563280000143
G sat2 ＝G dry1
the seventh elastic modulus is:
Figure RE-GDA0003221563280000144
G sat3 ＝G dry2
in the formula, K sat2 A sixth bulk modulus; k sat3 Is a seventh bulk modulus; g sat2 A sixth shear modulus; g sat3 Is a seventh bulk modulus; k dry1 A third bulk modulus; g dry1 Is the third shear modulus; k is dry2 A fourth bulk modulus; g dry2 Is the fourth shear modulus.
Further, the first seismic longitudinal wave velocity, the first seismic transverse wave velocity, the second seismic longitudinal wave velocity, the second seismic transverse wave velocity, the third seismic longitudinal wave velocity, and the third seismic transverse wave velocity are respectively:
Figure RE-GDA0003221563280000151
Figure RE-GDA0003221563280000152
Figure RE-GDA0003221563280000153
Figure RE-GDA0003221563280000154
Figure RE-GDA0003221563280000155
Figure RE-GDA0003221563280000156
in the formula, V pc_crack Is a first seismic longitudinal wave velocity; v pc_convex The second seismic longitudinal wave velocity; v pc_concave Is the third seismic longitudinal wave velocity; v sc_crack Is the first seismic shear wave velocity; v sc_convex The second seismic transverse wave velocity; v sc_concave Is the third seismic shear velocity; rho sat Is the density of the saturated fluid rock.
In particular, the amount of the solvent to be used,
Figure RE-GDA0003221563280000157
ρ i is the density of the rock matrix; rho f Is the density of the fluid.
Further, the first pore volume ratio is:
Figure RE-GDA0003221563280000161
the second pore volume ratio is:
Figure RE-GDA0003221563280000162
in the formula, w 1 A first pore volume ratio; w is a 2 A second pore volume ratio; v pm The measured seismic longitudinal wave velocity is obtained.
Further, the seismic shear wave velocity is:
V sc ＝V sc_convex -w 1 (V sc_convex -V sc_concave )
V sc ＝V sc_concave -w 1 (V sc_concave -V sc_crack )
in the formula, V sc Is the seismic shear wave velocity.
According to the embodiment of the invention, the first pore volume ratio and the second pore volume ratio are obtained by calculation according to the actually measured seismic longitudinal wave velocity and the first seismic longitudinal wave velocity, the second seismic longitudinal wave velocity and the third seismic longitudinal wave velocity, and the seismic transverse wave velocity is predicted according to the first pore volume ratio and the second pore volume ratio, so that the accuracy of seismic transverse wave velocity prediction can be further improved.
Further, in some embodiments of the present invention, obtaining a fluid volume modulus of the fluid is specifically: obtaining the fluid bulk modulus of the fluid through a British equation; wherein the Brewer's equation is:
Figure RE-GDA0003221563280000163
in the formula, K f Is the fluid bulk modulus of the fluid; k is w Is the bulk modulus of water; k g Is the bulk modulus of air; s w e Is the water saturation; e is the breve coefficient.
Specifically, in some embodiments of the present invention, e is 3.
Further, in some embodiments of the present invention, before S102, the method further includes:
obtaining a first elastic modulus of the rock matrix;
wherein, in some embodiments of the invention, as shown in fig. 3, obtaining the first elastic modulus of the rock matrix comprises:
s301, judging whether the mineral types in the rock matrix are more than one;
s303, if not, acquiring the first elastic modulus according to a preset parameter table;
and S303, if so, acquiring the initial elastic modulus of each mineral in the rock matrix according to a preset parameter table, and acquiring the first elastic modulus according to a Voigt-reus-Hill formula and the initial elastic modulus of each mineral.
Specifically, the preset parameter table is shown in table 1:
TABLE 1
Figure RE-GDA0003221563280000171
Namely: the initial bulk modulus and the initial shear modulus of various minerals can be obtained from table 1.
It should be understood that: the rock matrix may not only comprise mineral species as described in table 1, which is not described in detail herein, but the initial bulk modulus and the initial shear modulus of different mineral species may be obtained from laboratory measurements or well log interpretation.
Further, the Voigt-reus-Hill formula is:
Figure RE-GDA0003221563280000181
Figure RE-GDA0003221563280000182
Figure RE-GDA0003221563280000183
in the formula, M VRH A first elastic modulus of the rock matrix; f. of i Is the volume fraction of the ith mineral in the rock; m i Is the elastic modulus of the ith mineral in rock; m V Is the Voigt upper limit; m R The lower limit of reus.
Further, taking the three-fold system tight sandstone reservoir in the basin of alcdos, china as an example, the method for predicting the seismic transverse wave velocity of the invention and the PSP model and the DEM model in the prior art are used for respectively predicting the seismic transverse wave velocity, and the prediction results are shown in table 2:
TABLE 2 statistical table of errors between measured seismic shear wave velocities and predicted seismic shear wave velocities
Figure RE-GDA0003221563280000184
Figure RE-GDA0003221563280000191
In the table, M is the average value of relative errors, R mse Root mean square error, R 2 Is the correlation coefficient.
As can be seen from the above table, the result predicted by the method for predicting the seismic transverse wave velocity provided by the embodiment of the invention is more accurate than that in the prior art.
In order to better implement the method for predicting the seismic shear wave velocity of the tight sandstone reservoir in the embodiment of the present invention, on the basis of the method for predicting the seismic shear wave velocity of the tight sandstone reservoir, as shown in fig. 4, correspondingly, an apparatus for predicting the seismic shear wave velocity of the tight sandstone reservoir is further provided in the embodiment of the present invention, and the apparatus 400 for predicting the seismic shear wave velocity of the tight sandstone reservoir includes:
the model establishing module 401 is used for establishing a differential effective medium rock physical model and a porous rock physical model;
the first model calculation module 402 is used for obtaining a second elastic modulus of the dry rock containing the cracks according to the differential effective medium rock physical model and the first elastic modulus of the rock matrix;
the second model calculation module 403 is used for obtaining a third elastic modulus of the convex-edge-containing triangular dry rock and a fourth elastic modulus of the concave-edge-containing quadrilateral dry rock according to the porous rock physical model and the first elastic modulus of the rock matrix;
a theoretical velocity calculation module 404, configured to obtain a fluid volume modulus of the fluid, and obtain a first seismic longitudinal wave velocity, a first seismic transverse wave velocity, a second seismic longitudinal wave velocity, a second seismic transverse wave velocity, a third seismic longitudinal wave velocity, and a third seismic transverse wave velocity of the saturated fluid rock according to the second elastic module, the third elastic modulus, the fourth elastic modulus, and the fluid volume module, respectively;
a parameter determination module 405, configured to determine a first pore volume ratio and a second pore volume ratio according to the first seismic longitudinal wave velocity, the second seismic longitudinal wave velocity, and the third seismic longitudinal wave velocity;
and the seismic transverse wave velocity prediction module 406 is configured to predict the seismic transverse wave velocity according to the first pore volume ratio, the second pore volume ratio, the first seismic transverse wave velocity, the second seismic transverse wave velocity, and the third seismic transverse wave velocity.
According to the earthquake transverse wave velocity prediction device 400 of the tight sandstone reservoir, firstly, a differential effective medium rock physical model and a porous rock physical model are established through a model establishing module 401, then, a first model calculating module 402 and a second model calculating module 403 are used for obtaining a second elastic modulus of a dry rock containing cracks, a third elastic modulus of a dry rock containing a convex triangle and a fourth elastic modulus of a dry rock containing a concave quadrilateral by utilizing two rock physical models and a first elastic modulus of a rock matrix, then, a first pore volume ratio for representing pore volume ratios of the convex triangle and the concave quadrilateral and a second pore volume ratio for representing pore volume ratios of the concave triangle and the concave quadrilateral are determined through a parameter determining module 405 according to the third elastic modulus and the fourth elastic modulus, and finally, an earthquake transverse wave velocity is predicted through an earthquake transverse wave velocity prediction module 406 according to the first pore volume ratio and the second pore volume ratio. Compared with the prior art, the embodiment of the invention determines the first pore volume ratio for representing the pore volume ratio of the convex triangle and the concave quadrilateral and the second pore volume ratio for representing the pore volume ratio of the concave quadrilateral to the fracture through the two rock physical models, can describe the pore shape of the tight sandstone reservoir more truly, and finally predicts the seismic shear wave velocity through the first pore volume ratio and the second pore volume ratio, thereby improving the accuracy of the predicted seismic shear wave velocity.
The embodiment of the invention also provides computer equipment which integrates the device for predicting the seismic transverse wave velocity of any compact sandstone reservoir provided by the embodiment of the invention. Fig. 5 is a schematic diagram showing a structure of a computer device according to an embodiment of the present invention, specifically:
the computer device may include components such as a processor 501 of one or more processing cores, memory 502 of one or more computer-readable storage media, a power supply 503, and an input unit 504. Those skilled in the art will appreciate that the computer device configuration illustrated in FIG. 5 does not constitute a limitation of computer devices, and may include more or fewer components than those illustrated, or some components may be combined, or a different arrangement of components. Wherein:
the processor 501 is a control center of the computer device, connects various parts of the entire computer device by using various interfaces and lines, and performs various functions of the computer device and processes data by running or executing software programs and/or modules stored in the memory 502 and calling data stored in the memory 502, thereby monitoring the computer device as a whole. Optionally, processor 501 may include one or more processing cores; preferably, the processor 501 may integrate an application processor, which mainly handles an operating system, an operation user interface, an application program, and the like, and a modem processor, which mainly handles wireless communication. It will be appreciated that the modem processor described above may not be integrated into the processor 501.
The memory 502 may be used to store software programs and modules, and the processor 501 executes various functional applications and data processing by operating the software programs and modules stored in the memory 502. The memory 502 may mainly include a program storage area and a data storage area, wherein the program storage area may store an operating system, an application program required for at least one function (such as a sound playing function, an image playing function, etc.), and the like; the storage data area may store data created according to use of the computer device, and the like. Further, the memory 502 may include high speed random access memory, and may also include non-volatile memory, such as at least one magnetic disk storage device, flash memory device, or other volatile solid state storage device. Accordingly, the memory 502 may also include a memory controller to provide the processor 501 with access to the memory 502.
The computer device further comprises a power supply 503 for supplying power to the various components, and preferably, the power supply 503 is logically connected to the processor 501 through a power management system, so that the functions of charging, discharging, power consumption management and the like are managed through the power management system. The power supply 503 may also include any component of one or more dc or ac power sources, recharging systems, power failure detection circuitry, power converters or inverters, power status indicators, and the like.
The computer device may also include an input unit 504, where the input unit 504 may be used to receive entered numeric or character information and generate keyboard, mouse, joystick, optical or trackball signal inputs related to operating user settings and function controls.
Although not shown, the computer device may further include a display unit and the like, which are not described in detail herein. Specifically, in this embodiment, the processor 501 in the computer device loads the executable file corresponding to the process of one or more application programs into the memory 502 according to the following instructions, and the processor 501 runs the application programs stored in the memory 502, thereby implementing various functions as follows:
establishing a differential effective medium rock physical model and a porous rock physical model;
obtaining a second elastic modulus of the dry rock containing the cracks according to the differential effective medium rock physical model and the first elastic modulus of the rock matrix;
obtaining a third elastic modulus of the convex-edge-containing triangular dry rock and a fourth elastic modulus of the concave-edge-containing quadrilateral dry rock according to the porous rock physical model and the first elastic modulus;
obtaining the fluid volume modulus of the fluid, and respectively obtaining a first seismic longitudinal wave velocity, a first seismic transverse wave velocity, a second seismic longitudinal wave velocity, a second seismic transverse wave velocity, a third seismic longitudinal wave velocity and a third seismic transverse wave velocity of the saturated fluid rock according to the second elastic module, the third elastic modulus, the fourth elastic modulus and the fluid volume module;
determining a first pore volume ratio and a second pore volume ratio according to the first seismic longitudinal wave velocity, the second seismic longitudinal wave velocity and the third seismic longitudinal wave velocity;
and predicting the seismic transverse wave velocity according to the first pore volume ratio, the second pore volume ratio, the first seismic transverse wave velocity, the second seismic transverse wave velocity and the third seismic transverse wave velocity.
It will be understood by those skilled in the art that all or part of the steps of the methods of the above embodiments may be performed by instructions or by associated hardware controlled by the instructions, which may be stored in a computer readable storage medium and loaded and executed by a processor.
To this end, an embodiment of the present invention provides a computer-readable storage medium, which may include: read Only Memory (ROM), random Access Memory (RAM), magnetic or optical disks, and the like. The seismic shear wave velocity prediction method for the tight sandstone reservoir provided by the embodiment of the invention comprises a step of carrying out seismic shear wave velocity prediction on the tight sandstone reservoir, wherein the step is carried out by a processor. For example, the computer program may be loaded by a processor to perform the steps of:
establishing a differential effective medium rock physical model and a porous rock physical model;
obtaining a second elastic modulus of the dry rock containing the cracks according to the differential effective medium rock physical model and the first elastic modulus of the rock matrix;
obtaining a third elastic modulus of the convex-edge-containing triangular dry rock and a fourth elastic modulus of the concave-edge-containing quadrilateral dry rock according to the porous rock physical model and the first elastic modulus;
obtaining the fluid volume modulus of the fluid, and respectively obtaining a first seismic longitudinal wave velocity, a first seismic transverse wave velocity, a second seismic longitudinal wave velocity, a second seismic transverse wave velocity, a third seismic longitudinal wave velocity and a third seismic transverse wave velocity of the saturated fluid rock according to the second elastic module, the third elastic modulus, the fourth elastic modulus and the fluid volume module;
determining a first pore volume ratio and a second pore volume ratio according to the first seismic longitudinal wave velocity, the second seismic longitudinal wave velocity and the third seismic longitudinal wave velocity;
and predicting the seismic transverse wave velocity according to the first pore volume ratio, the second pore volume ratio, the first seismic transverse wave velocity, the second seismic transverse wave velocity and the third seismic transverse wave velocity.
The above-described embodiments of the present invention should not be construed as limiting the scope of the present invention. Any other corresponding changes and modifications made according to the technical idea of the present invention should be included in the protection scope of the claims of the present invention.

Claims (3)

1. A method for predicting seismic shear wave velocity of a tight sandstone reservoir is characterized by comprising the following steps:
establishing a differential effective medium rock physical model and a porous rock physical model;
obtaining a second elastic modulus of the dry rock containing the cracks according to the differential effective medium rock physical model and the first elastic modulus of the rock matrix;
obtaining a third elastic modulus of the convex-edge-containing triangular dry rock and a fourth elastic modulus of the concave-edge-containing quadrilateral dry rock according to the porous rock physical model and the first elastic modulus;
obtaining the fluid volume modulus of the fluid, and respectively obtaining a first seismic longitudinal wave velocity, a first seismic transverse wave velocity, a second seismic longitudinal wave velocity, a second seismic transverse wave velocity, a third seismic longitudinal wave velocity and a third seismic transverse wave velocity of the saturated fluid rock according to the second elastic modulus, the third elastic modulus, the fourth elastic modulus and the fluid volume modulus;
determining a first pore volume ratio and a second pore volume ratio according to the first seismic longitudinal wave velocity, the second seismic longitudinal wave velocity and the third seismic longitudinal wave velocity;
predicting the seismic transverse wave velocity according to the first pore volume ratio, the second pore volume ratio, the first seismic transverse wave velocity, the second seismic transverse wave velocity and the third seismic transverse wave velocity;
the first elastic modulus comprises a first bulk modulus and a first shear modulus, and the second elastic modulus comprises a second bulk modulus and a second shear modulus; the differential effective medium rock physical model is as follows:
Figure FDA0003892748820000011
Figure FDA0003892748820000021
Figure FDA0003892748820000022
Figure FDA0003892748820000023
Figure FDA0003892748820000024
in the formula, K * DEM (y) is the second bulk modulus; g * DEM (y) is the second shear modulus; k i Is the first bulk modulus; g i Is the first shear modulus; y is the volume fraction of cracks in the dry rock;
Figure FDA0003892748820000028
is an influencing factor related to a crack in the dry rock; α is a first equivalent shape factor;
the third modulus of elasticity comprises a third bulk modulus and a third shear modulus, and the fourth modulus of elasticity comprises a fourth bulk modulus and a fourth shear modulus; the physical model of the porous rock is as follows:
Figure FDA0003892748820000025
Figure FDA0003892748820000026
wherein g is a second equivalent shape factor;
Figure FDA0003892748820000027
is porosity; k is dry Is the bulk modulus; g dry Is the shear modulus;
the step of respectively obtaining a first seismic longitudinal wave velocity, a first seismic transverse wave velocity, a second seismic longitudinal wave velocity, a second seismic transverse wave velocity, a third seismic longitudinal wave velocity and a third seismic transverse wave velocity of the saturated fluid rock according to the second elastic modulus, the third elastic modulus, the fourth elastic modulus and the fluid bulk modulus comprises:
obtaining a fifth elastic modulus of the saturated fluid rock according to the second elastic modulus and the fluid volume modulus;
obtaining a sixth elastic modulus of the saturated fluid rock according to the third elastic modulus and the fluid volume modulus;
obtaining a seventh elastic modulus of the saturated fluid rock according to the fourth elastic modulus and the fluid volume modulus;
sequentially obtaining the first seismic longitudinal wave velocity, the first seismic transverse wave velocity, the second seismic longitudinal wave velocity, the second seismic transverse wave velocity, the third seismic longitudinal wave velocity and the third seismic transverse wave velocity according to the fifth elastic modulus, the sixth elastic modulus and the seventh elastic modulus;
the fifth elastic modulus comprises a fifth bulk modulus and a fifth shear modulus, and the fifth elastic modulus is:
Figure FDA0003892748820000031
Figure FDA0003892748820000032
in the formula, K sat1 Is the fifth bulk modulus; g sat1 Is the fifth shear modulus;
the sixth modulus of elasticity comprises a sixth bulk modulus and a sixth shear modulus, and the seventh modulus of elasticity comprises a seventh bulk modulus and a seventh shear modulus; the sixth elastic modulus is:
Figure FDA0003892748820000033
G sat2 ＝G dry1
the seventh elastic modulus is:
Figure FDA0003892748820000041
G sat3 ＝G dry2
in the formula, K sat2 Is the sixth bulk modulus; k sat3 Is the seventh bulk modulus; g sat2 Is the sixth shear modulus; g sat3 Is the seventh shear modulus; k is dry1 Is the third bulk modulus; g dry1 Is the third shear modulus; k dry2 Is the fourth bulk modulus; g dry2 Is the fourth shear modulus;
the first seismic longitudinal wave velocity, the first seismic transverse wave velocity, the second seismic longitudinal wave velocity, the second seismic transverse wave velocity, the third seismic longitudinal wave velocity and the third seismic transverse wave velocity are respectively as follows:
Figure FDA0003892748820000042
Figure FDA0003892748820000043
Figure FDA0003892748820000044
Figure FDA0003892748820000045
Figure FDA0003892748820000046
Figure FDA0003892748820000051
in the formula, V pc_crack Is the first seismic longitudinal wave velocity; v pc_convex The second seismic longitudinal wave velocity; v pc_concave Is the third seismic longitudinal wave velocity; v sc_crack Is the first seismic shear velocity; v sc_convex The second seismic transverse wave velocity; v sc_concave Is the third seismic shear wave velocity; ρ is a unit of a gradient sat Density of saturated fluid rock;
the first pore volume ratio is:
Figure FDA0003892748820000052
the second pore volume ratio is:
Figure FDA0003892748820000053
in the formula, w 1 Is the first pore volume ratio; w is a 2 Is the second pore volume ratio; v pm The measured seismic longitudinal wave velocity is obtained;
the seismic transverse wave velocity is as follows:
V sc ＝V sc_convex -w 1 (V sc_convex -V sc_concave )
V sc ＝V sc_concave -w 2 (V sc_concave -V sc_crack )
in the formula, V sc Is the seismic shear wave velocity.
2. The method for predicting seismic shear wave velocity of tight sandstone reservoir of claim 1, wherein before the step of obtaining the second elastic modulus of the dry rock containing fractures according to the differential effective medium petrophysical model and the first elastic modulus of the rock matrix, the method further comprises:
obtaining a first elastic modulus of the rock matrix;
wherein said obtaining a first modulus of elasticity of the rock matrix comprises:
determining whether there is more than one mineral species in the rock matrix;
if not, acquiring the first elastic modulus according to a preset parameter table;
if so, acquiring the initial elastic modulus of each mineral in the rock matrix according to a preset parameter table, and acquiring the first elastic modulus according to a Voigt-Reuss-Hill formula and the initial elastic modulus of each mineral.
3. The device for predicting the seismic transverse wave velocity of the tight sandstone reservoir is characterized by comprising the following components:
the model building module is used for building a differential effective medium rock physical model and a porous rock physical model;
the first model calculation module is used for obtaining a second elastic modulus of the crack-containing dry rock according to the differential effective medium rock physical model and the first elastic modulus of the rock matrix;
the second model calculation module is used for obtaining a third elastic modulus of the convex-edge-containing triangular dry rock and a fourth elastic modulus of the concave-edge-containing quadrilateral dry rock according to the porous rock physical model and the first elastic modulus of the rock matrix;
the theoretical velocity calculation module is used for obtaining the fluid volume modulus of the fluid and respectively obtaining a first seismic longitudinal wave velocity, a first seismic transverse wave velocity, a second seismic longitudinal wave velocity, a second seismic transverse wave velocity, a third seismic longitudinal wave velocity and a third seismic transverse wave velocity of the saturated fluid rock according to the second elastic modulus, the third elastic modulus, the fourth elastic modulus and the fluid volume modulus;
the parameter determining module is used for determining a first pore volume ratio and a second pore volume ratio according to the first seismic longitudinal wave velocity, the second seismic longitudinal wave velocity and the third seismic longitudinal wave velocity;
the earthquake transverse wave velocity prediction module is used for predicting the earthquake transverse wave velocity according to the first pore volume ratio, the second pore volume ratio, the first earthquake transverse wave velocity, the second earthquake transverse wave velocity and the third earthquake transverse wave velocity;
the first elastic modulus comprises a first bulk modulus and a first shear modulus, and the second elastic modulus comprises a second bulk modulus and a second shear modulus; the differential effective medium rock physical model is as follows:
Figure FDA0003892748820000071
Figure FDA0003892748820000075
Figure FDA0003892748820000072
Figure FDA0003892748820000073
Figure FDA0003892748820000074
in the formula, K * DEM (y) is the second bulk modulus; g * DEM (y) is the second shear modulus; k is i Is the first bulk modulus; g i Is the first shear modulus; y is the volume fraction of the fractures in the dry rock;
Figure FDA0003892748820000084
is an influencing factor related to a crack in the dry rock; α is a first equivalent shape factor;
the third modulus of elasticity comprises a third bulk modulus and a third shear modulus, and the fourth modulus of elasticity comprises a fourth bulk modulus and a fourth shear modulus; the physical model of the porous rock is as follows:
Figure FDA0003892748820000081
Figure FDA0003892748820000082
wherein g is a second equivalent shape factor;
Figure FDA0003892748820000083
is porosity; k dry Is the bulk modulus; g dry Is the shear modulus;
the obtaining of the first seismic longitudinal wave velocity, the first seismic transverse wave velocity, the second seismic longitudinal wave velocity, the second seismic transverse wave velocity, the third seismic longitudinal wave velocity, and the third seismic transverse wave velocity of the saturated fluid rock according to the second elastic modulus, the third elastic modulus, the fourth elastic modulus, and the fluid bulk modulus includes:
obtaining a fifth elastic modulus of the saturated fluid rock according to the second elastic modulus and the fluid volume modulus;
obtaining a sixth elastic modulus of the saturated fluid rock according to the third elastic modulus and the fluid volume modulus;
obtaining a seventh elastic modulus of the saturated fluid rock according to the fourth elastic modulus and the fluid volume modulus;
sequentially obtaining the first seismic longitudinal wave velocity, the first seismic transverse wave velocity, the second seismic longitudinal wave velocity, the second seismic transverse wave velocity, the third seismic longitudinal wave velocity and the third seismic transverse wave velocity according to the fifth elastic modulus, the sixth elastic modulus and the seventh elastic modulus;
the fifth elastic modulus comprises a fifth bulk modulus and a fifth shear modulus, and the fifth elastic modulus is:
Figure FDA0003892748820000091
Figure FDA0003892748820000092
in the formula, K sat1 Is the fifth bulk modulus; g sat1 Is the fifth shear modulus;
the sixth modulus of elasticity comprises a sixth bulk modulus and a sixth shear modulus, and the seventh modulus of elasticity comprises a seventh bulk modulus and a seventh shear modulus; the sixth elastic modulus is:
Figure FDA0003892748820000093
G sat2 ＝G dry1
the seventh elastic modulus is:
Figure FDA0003892748820000094
G sat3 ＝G dry2
in the formula, K sat2 Is the sixth bulk modulus; k sat3 Is the seventh bulk modulus; g sat2 Is the sixth shear modulus; g sat3 Is the seventh shear modulus; k dry1 Is the third bulk modulus; g dry1 Is the third shear modulus; k dry2 Is the fourth bulk modulus; g dry2 Is the fourth shear modulus;
the first seismic longitudinal wave velocity, the first seismic transverse wave velocity, the second seismic longitudinal wave velocity, the second seismic transverse wave velocity, the third seismic longitudinal wave velocity, and the third seismic transverse wave velocity are respectively:
Figure FDA0003892748820000101
Figure FDA0003892748820000102
Figure FDA0003892748820000103
Figure FDA0003892748820000104
Figure FDA0003892748820000105
Figure FDA0003892748820000106
in the formula, V pc_crack Is the first seismic longitudinal wave velocity; v pc_convex The second seismic longitudinal wave velocity; v pc_concave Is the third seismic longitudinal wave velocity; v sc_crack Is the first seismic shear wave velocity; v sc_convex The second seismic transverse wave velocity; v sc_concave Is the third seismic shear velocity; rho sat Is the density of the saturated fluid rock;
the first pore volume ratio is:
Figure FDA0003892748820000107
the second pore volume ratio is:
Figure FDA0003892748820000111
in the formula, w 1 Is the first pore volume ratio; w is a 2 Is the second pore volume ratio; v pm The measured seismic longitudinal wave velocity is obtained;
the seismic transverse wave velocity is as follows:
V sc ＝V sc_convex -w 1 (V sc_convex -V sc_concave )
V sc ＝V sc_concave -w 2 (V sc_concave -V sc_crack )
in the formula, V sc And the seismic transverse wave velocity is used as the seismic transverse wave velocity.
CN202110754161.XA 2021-07-01 2021-07-01 Method and device for predicting seismic transverse wave velocity of tight sandstone reservoir Active CN113552626B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110754161.XA CN113552626B (en) 2021-07-01 2021-07-01 Method and device for predicting seismic transverse wave velocity of tight sandstone reservoir

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110754161.XA CN113552626B (en) 2021-07-01 2021-07-01 Method and device for predicting seismic transverse wave velocity of tight sandstone reservoir

Publications (2)

Publication Number Publication Date
CN113552626A CN113552626A (en) 2021-10-26
CN113552626B true CN113552626B (en) 2023-01-17

Family

ID=78131269

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110754161.XA Active CN113552626B (en) 2021-07-01 2021-07-01 Method and device for predicting seismic transverse wave velocity of tight sandstone reservoir

Country Status (1)

Country Link
CN (1) CN113552626B (en)

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6226390B1 (en) * 1996-11-18 2001-05-01 Institute Francais Du Petrole Method for modeling the pore-size distribution of a porous sample of variable porosity
CN103760081A (en) * 2013-12-31 2014-04-30 中国石油天然气股份有限公司 Gas pool prediction method and system of carbonate reservoir based on pore structure characteristics
EP2816377A1 (en) * 2013-06-20 2014-12-24 IFP Energies nouvelles Method for predicting the amount and composition of fluids produced by mineral reactions operating in a sedimentary basin
CN105974474A (en) * 2016-07-07 2016-09-28 中国石油大学(华东) Method for evaluating reservoir pore type based on pore replacement
WO2018010628A1 (en) * 2016-07-15 2018-01-18 河海大学 Seismic rock physics inversion method based on a large area tight reservoir
CN109581490A (en) * 2017-09-28 2019-04-05 中国石油化工股份有限公司 A kind of Wave velocity estimation method and system of compact reservoir
CN112285801A (en) * 2019-07-24 2021-01-29 中国石油天然气股份有限公司 Method, device and system for predicting permeability of carbonate reservoir

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8731889B2 (en) * 2010-03-05 2014-05-20 Schlumberger Technology Corporation Modeling hydraulic fracturing induced fracture networks as a dual porosity system
CN110824556B (en) * 2019-10-22 2022-03-01 中国石油天然气股份有限公司 Rock physical model building method and application of unconventional tight sandstone reservoir

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6226390B1 (en) * 1996-11-18 2001-05-01 Institute Francais Du Petrole Method for modeling the pore-size distribution of a porous sample of variable porosity
EP2816377A1 (en) * 2013-06-20 2014-12-24 IFP Energies nouvelles Method for predicting the amount and composition of fluids produced by mineral reactions operating in a sedimentary basin
CN103760081A (en) * 2013-12-31 2014-04-30 中国石油天然气股份有限公司 Gas pool prediction method and system of carbonate reservoir based on pore structure characteristics
CN105974474A (en) * 2016-07-07 2016-09-28 中国石油大学(华东) Method for evaluating reservoir pore type based on pore replacement
WO2018010628A1 (en) * 2016-07-15 2018-01-18 河海大学 Seismic rock physics inversion method based on a large area tight reservoir
CN109581490A (en) * 2017-09-28 2019-04-05 中国石油化工股份有限公司 A kind of Wave velocity estimation method and system of compact reservoir
CN112285801A (en) * 2019-07-24 2021-01-29 中国石油天然气股份有限公司 Method, device and system for predicting permeability of carbonate reservoir

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
A shale rock physics model for analysis of brittleness index, mineralogy and porosity in the Barnett Shale;Zhiqi Guo et al.;《J. Geophys. Eng.》;20131231;第10卷;全文 *
Zhishui Liu et al..The differential Kuster–Toksöz rock physics model for predicting S-wave velocity.《J. Geophys. Eng.》.2015,第12卷第839-848页. *
富有机质岩石横波速度预测方法;刘致水等;《石油地球物理勘探》;20210228;第56卷(第1期);第127-136页 *
碳酸盐岩孔洞储层地震岩石物理建模及应用;曹晓初等;《地球物理学进展》;20191231;第34卷(第6期);第2239-2246页 *

Also Published As

Publication number Publication date
CN113552626A (en) 2021-10-26

Similar Documents

Publication Publication Date Title
Zängl et al. The ICON (ICOsahedral Non‐hydrostatic) modelling framework of DWD and MPI‐M: Description of the non‐hydrostatic dynamical core
Pan et al. Improved automatic operational modal analysis method and application to large-scale bridges
CN105549088B (en) The recognition methods of Fractured tight sand Mesosphere and device
CN111695308B (en) Method for predicting change characteristics of longitudinal wave speed along with temperature change
CN113552626B (en) Method and device for predicting seismic transverse wave velocity of tight sandstone reservoir
CN104316962B (en) Rock pore structure parameter inversion method and its system
CN108008444A (en) A kind of method and device of definite petrophysical model
EP2678718A2 (en) Computerized method for the estimation of a value for at least a parameter of a hydrocarbon-producing region, for planning the operation and operating the region
CN109211750A (en) Based on micro-, microscopic structural parameters Chloride Diffusion Coefficient in Concrete prediction technique
CN116738794A (en) Rock physical numerical simulation method, device, equipment and medium for pore fracture medium
CN113552627B (en) Self-adaptive earthquake transverse wave velocity prediction method and device
CN109799164B (en) Method and device for determining adsorption gas amount
CN112347707B (en) Method and device for characterizing flow heterogeneity of oil reservoir
CN112285801B (en) Method, device and system for predicting permeability of carbonate reservoir
CN107015289B (en) Trap evaluation stock number determines method and device
Garitte et al. Modelling benchmark of a laboratory drying test in Opalinus Clay
CN112129482B (en) Shale gas flow capacity analysis method and device in shale reservoir matrix
CN114675346B (en) Method, device and equipment for analyzing seismic logging performance and readable storage medium
CN115079261B (en) Compact sandstone gas reservoir evaluation method based on multiparameter quantitative interpretation template
LU503609B1 (en) Simulation method, system and processing terminal for simulating the numerical model of interbedded rock mass
CN115392090B (en) Prediction method of seismic wave frequency dispersion and attenuation characteristics based on three-dimensional digital rock core
CN114065462B (en) Rock pore structure evolution analysis method and system
CN117421939B (en) Shale oil fracture system simulation agent method based on track piecewise linearization
CN109781039B (en) Method for quantitatively predicting longitudinal dimension of karst cave
CN114458305A (en) Method and device for determining pore structure coefficient

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
CB03 Change of inventor or designer information
CB03 Change of inventor or designer information

Inventor after: Liu Junzhou

Inventor after: Liu Zhishui

Inventor after: Bao Qianzong

Inventor after: Shi Lei

Inventor after: Wang Zhenyu

Inventor before: Liu Zhishui

Inventor before: Bao Qianzong

Inventor before: Liu Junzhou

Inventor before: Shi Lei

Inventor before: Wang Zhenyu

TA01 Transfer of patent application right
TA01 Transfer of patent application right

Effective date of registration: 20221123

Address after: 100000 22 Chaoyangmen North Street, Chaoyang District, Beijing.

Applicant after: CHINA PETROLEUM & CHEMICAL Corp.

Applicant after: CHANG'AN University

Address before: No. 126, south section of Yanta Road, Yanta District, Xi'an, Shaanxi 710000

Applicant before: CHANG'AN University

GR01 Patent grant
GR01 Patent grant