CN116819607B - Transverse wave prediction method and system based on tunnel curvature - Google Patents
Transverse wave prediction method and system based on tunnel curvature Download PDFInfo
- Publication number
- CN116819607B CN116819607B CN202310548829.4A CN202310548829A CN116819607B CN 116819607 B CN116819607 B CN 116819607B CN 202310548829 A CN202310548829 A CN 202310548829A CN 116819607 B CN116819607 B CN 116819607B
- Authority
- CN
- China
- Prior art keywords
- rock
- pore
- saturated
- fluid
- equivalent
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 36
- 239000011435 rock Substances 0.000 claims abstract description 289
- 239000011148 porous material Substances 0.000 claims abstract description 111
- 239000012530 fluid Substances 0.000 claims abstract description 102
- 229920006395 saturated elastomer Polymers 0.000 claims abstract description 98
- 239000011159 matrix material Substances 0.000 claims abstract description 25
- 238000004364 calculation method Methods 0.000 claims description 21
- 238000005452 bending Methods 0.000 claims description 2
- 229910052500 inorganic mineral Inorganic materials 0.000 abstract description 4
- 239000011707 mineral Substances 0.000 abstract description 4
- 230000008569 process Effects 0.000 abstract description 4
- 238000007796 conventional method Methods 0.000 description 5
- 238000004458 analytical method Methods 0.000 description 4
- 230000008901 benefit Effects 0.000 description 4
- 238000011160 research Methods 0.000 description 4
- 239000004927 clay Substances 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 238000011156 evaluation Methods 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- BVKZGUZCCUSVTD-UHFFFAOYSA-L Carbonate Chemical compound [O-]C([O-])=O BVKZGUZCCUSVTD-UHFFFAOYSA-L 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000013211 curve analysis Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000003631 expected effect Effects 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000004576 sand Substances 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention provides a transverse wave prediction method and a transverse wave prediction system based on tunnel curvature, wherein the method comprises the following steps: constructing a rock equivalent volume model of the reservoir; characterizing the pore flatness of the rock in the reservoir by resistivity, porosity and pore tortuosity of the pore-saturated fluid rock based on the rock equivalent volume model; calculating the elastic modulus of the rock framework and the rock matrix; calculating the equivalent bulk modulus and the equivalent shear modulus of the pore-saturated fluid rock in the reservoir according to the rock skeleton and the rock matrix elastic modulus; and predicting the shear wave velocity of the rock in the reservoir according to the equivalent shear modulus and the equivalent density of the pore-saturated fluid rock. The rock equivalent volume model provided by the invention can be used for predicting the longitudinal wave speed and the transverse wave speed of the rock full of the pore fluid well, reduces the requirement on the input precision of porosity, mineral components and the like in the use process of the traditional Xu-White model, and has strong practicability.
Description
Technical Field
The invention belongs to the technical field of geophysics, and particularly relates to a transverse wave prediction method and a transverse wave prediction system based on pore tortuosity.
Background
Rock physics experiments show that the acoustic wave speed of a reservoir is sensitive to porosity, fluid properties and clay content, so that the longitudinal and transverse wave speed is widely applied to seismic and logging interpretation works such as AVA, AVO attribute analysis, pre-stack seismic inversion, reservoir prediction, fluid identification, reservoir rock mechanical parameter calculation and the like. However, in actual work, the transverse wave speed logging is expensive, and many old well research areas lack transverse wave speed information, so that the development of transverse wave speed prediction research has economic benefit and practical significance.
In the aspect of transverse wave prediction, researchers at home and abroad have conducted a great deal of research work, and the technical routes of the research work are mainly divided into two categories: one is to use a great amount of logging information such as longitudinal and transverse waves, density, gamma and the like to carry out statistical fitting, and construct an empirical formula between the transverse waves and other parameters and logging information, such as: based on known logging data, castagna et al establish an empirical regression formula of carbonate and sand shale longitudinal and transverse wave velocities in a saturated state; han et al further fit the longitudinal and transverse wave relationship taking into account the effects of porosity and clay content. The method is based on petrophysical modeling, numerical simulation is carried out on pore morphology and the like through known information, a petrophysical model is built, and transverse wave speed prediction is carried out by combining with a fluctuation theory, such as: xu and White are combined with Gassmann equation, K-T model and differential equivalent medium theory to provide a method for predicting transverse wave speed by utilizing porosity and clay content; keys et al improve on the Xu-White model and calculate the longitudinal and transverse wave velocities by solving a linear ordinary differential equation to determine the elastic modulus of the rock skeleton.
In practical applications, the above models are not adequate. The method for constructing the well formula by statistical fitting relies on experience relation, has strong regional limitation, and has extremely small application range; the petrophysical modeling method is mostly based on the improvement of an Xu-White model, a large amount of parameter input is needed for speed calculation, most parameters are not easy to directly obtain, parameter extraction is needed by combining a laboratory petrophysical experiment and log curve analysis, and the uncertainty of estimation is increased due to the simplification of a theoretical model to a rock structure, so that a large error exists in a prediction result.
Disclosure of Invention
Aiming at the problems, the invention provides a transverse wave prediction method and a transverse wave prediction system based on tunnel tortuosity, which provide powerful technical support for mechanical evaluation and seismic inversion of stratum rock lacking transverse wave logging data.
A transverse wave prediction method based on tunnel curvature comprises the following steps:
constructing a rock equivalent volume model of the reservoir, wherein the rock is pore-fluid-saturated rock;
Characterizing the pore flatness of the rock in the reservoir by resistivity, porosity and pore tortuosity of the pore-saturated fluid rock based on the rock equivalent volume model;
According to a first expression of the elastic modulus of the rock skeleton and a second expression of the correlation constant of the elastic modulus of the rock skeleton and the porosity, calculating to obtain the elastic modulus of the rock skeleton and the rock matrix;
calculating the equivalent bulk modulus and the equivalent shear modulus of the pore-saturated fluid rock in the reservoir according to the rock skeleton and the rock matrix elastic modulus;
And predicting the shear wave velocity of the rock in the reservoir according to the equivalent shear modulus and the equivalent density of the pore-saturated fluid rock.
Further, constructing an equivalent volume model of the pore-saturated fluid rock of the reservoir, which is specifically as follows:
The cross-sectional area of the saturated pore fluid rock is defined as A, the length of the saturated pore fluid rock is defined as L, the length of the equivalent throat of the saturated pore fluid rock is defined as L w, and the cross-sectional area of the equivalent throat of the saturated pore fluid rock is defined as A w.
Further, based on the rock equivalent volume model, the pore volume of the rock in the reservoir is characterized by the resistivity and pore tortuosity of the pore-saturated fluid rock as follows:
Based on the rock equivalent volume model, constructing a first relation of rock resistivity saturated with pore fluid according to a resistance parallel calculation formula;
determining a second relationship between the saturated water-free rock resistivity and the pore fluid resistivity from the first relationship;
and defining the pore canal curvature of the rock, and determining the porosity and the flat rate of the rock in the reservoir according to the porosity of the rock and the second relation.
Further, the first relation is specifically as follows:
Where r 0、rma and r w are water-saturated pure rock resistance, rock matrix resistance, and pore fluid resistance, respectively.
Further, a second relation between the resistivity of saturated pure rock and the resistivity of pore fluid is determined according to the first relation, and is specifically as follows:
defining the rock skeleton resistance to be infinite, solving a first relation according to the relation between the resistivity and the resistance to obtain a second relation, wherein the specific relation is as follows:
Where R o represents the electrical resistance of the rock, R w represents the pore fluid resistance, A represents the cross-sectional area of the saturated pore fluid rock, L represents the length of the saturated pore fluid rock, L w represents the length of the equivalent throat of the saturated pore fluid rock, and A w represents the cross-sectional area of the equivalent throat of the saturated pore fluid rock.
Further, the curvature of the pore canal of the rock is thatPorosity of rock/>The method comprises the following steps:
where A represents the cross-sectional area of the saturated pore fluid rock, L represents the length of the saturated pore fluid rock, L w represents the length of the equivalent throat of the saturated pore fluid rock, and A w represents the cross-sectional area of the equivalent throat of the saturated pore fluid rock.
Further, according to the porosity of the rock and the second relation, the porosity flat of the rock in the reservoir is determined, specifically as follows:
according to the porosity of the rock and the second relation, the relation between the pore canal curvature of the rock, the rock resistance, the pore fluid resistance and the porosity of the rock is determined, and the relation is specifically as follows:
the logarithm of the curvature of the pore canal is obtained, and the porosity flat rate alpha of the rock is obtained, specifically as follows:
Where R o represents the electrical resistance of the rock, R w represents the pore fluid resistance, L represents the length of the pore fluid-saturated rock, L w represents the length of the equivalent throat of the pore fluid-saturated rock, Representing the porosity of the rock.
Further, the method also comprises the following steps:
and predicting the longitudinal wave velocity of the rock in the reservoir according to the equivalent bulk modulus, the equivalent shear modulus and the equivalent density of the pore-saturated fluid rock.
The invention also provides a transverse wave prediction system based on the curvature of the pore canal, which comprises the following steps:
The model building module is used for building a rock equivalent volume model of the reservoir, wherein the rock is pore-fluid-saturated rock;
a first calculation module for characterizing a pore flatness of the rock in the reservoir by resistivity, porosity and pore tortuosity of the pore-saturated fluid rock based on the rock equivalent volume model;
the second calculation module is used for calculating and obtaining the elastic modulus of the rock skeleton and the rock matrix according to the first expression of the elastic modulus of the rock skeleton and the second expression of the related constant of the elastic modulus of the rock skeleton and the porosity;
the third calculation module is used for calculating the equivalent bulk modulus and the equivalent shear modulus of the pore-saturated fluid rock in the reservoir according to the elastic modulus of the rock skeleton and the rock matrix;
And the transverse wave prediction module is used for predicting the transverse wave speed of the rock in the reservoir according to the equivalent shear modulus and the equivalent density of the pore-saturated fluid rock.
Further, the model building module is specifically configured to:
The cross-sectional area of the saturated pore fluid rock is defined as A, the length of the saturated pore fluid rock is defined as L, the length of the equivalent throat of the saturated pore fluid rock is defined as L w, and the cross-sectional area of the equivalent throat of the saturated pore fluid rock is defined as A w.
Further, the first computing module is specifically configured to:
Based on the rock equivalent volume model, constructing a first relation of rock resistivity saturated with pore fluid according to a resistance parallel calculation formula;
determining a second relationship between the saturated water-free rock resistivity and the pore fluid resistivity from the first relationship;
and defining the pore canal curvature of the rock, and determining the porosity and the flat rate of the rock in the reservoir according to the porosity of the rock and the second relation.
The invention has the beneficial effects that: the rock equivalent volume model provided by the invention can be used for predicting the longitudinal wave speed and the transverse wave speed of the rock full of the pore fluid well, reduces the requirement on the input precision of porosity, mineral components and the like in the use process of the traditional Xu-White model, and has strong practicability.
Additional features and advantages of the invention will be set forth in the description which follows, and in part will be obvious from the description, or may be learned by practice of the invention. The objectives and other advantages of the invention may be realized and attained by the structure particularly pointed out in the written description and claims hereof as well as the appended drawings.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions of the prior art, the following description will briefly explain the drawings used in the embodiments or the description of the prior art, and it is obvious that the drawings in the following description are some embodiments of the present invention, and other drawings can be obtained according to these drawings without inventive effort for a person skilled in the art.
FIG. 1 shows a flow diagram of a method for shear wave prediction based on tunnel tortuosity according to an embodiment of the invention;
FIG. 2 shows a rock equivalent volume model of a reservoir according to an embodiment of the invention;
FIG. 3 is a schematic flow diagram of a shear wave prediction system based on tunnel curvature according to an embodiment of the present invention;
FIG. 4 is a graph showing a comparison of shear wave predictions with conventional method predictions in accordance with an embodiment of the present invention;
FIG. 5 is a graph showing a cross wave prediction result error analysis according to a conventional method;
FIG. 6 is a schematic diagram of a shear wave prediction error analysis according to an embodiment of the present invention.
Detailed Description
For the purpose of making the objects, technical solutions and advantages of the embodiments of the present invention more apparent, the technical solutions of 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, and it is apparent that the described embodiments are some embodiments of the present invention, but not all embodiments of the present invention. All other embodiments, which can be made by those skilled in the art based on the embodiments of the invention without making any inventive effort, are intended to be within the scope of the invention.
The embodiment of the invention provides a transverse wave prediction method and a system based on pore tortuosity, which are used for calculating the pore flat rate of a reservoir based on analysis of an equivalent volume model, calculating the equivalent volume modulus and shear modulus of a pore-saturated fluid rock by combining a K-T model and an Xu-White model, and further calculating the velocity of longitudinal waves and transverse waves.
As shown in fig. 1, a transverse wave prediction method based on the tortuosity of a duct comprises the following steps: constructing a rock equivalent volume model of the reservoir; characterizing the pore flatness of the rock in the reservoir by resistivity, porosity and pore tortuosity of the pore-saturated fluid rock based on the rock equivalent volume model; according to a first expression of the elastic modulus of the rock skeleton and a second expression of the correlation constant of the elastic modulus of the rock skeleton and the porosity, calculating to obtain the elastic modulus of the rock skeleton and the rock matrix; calculating the equivalent bulk modulus and the equivalent shear modulus of the pore-saturated fluid rock in the reservoir according to the rock skeleton and the rock matrix elastic modulus; predicting the shear wave velocity of the rock in the reservoir according to the equivalent shear modulus and the equivalent density of the pore-saturated fluid rock, wherein the shear wave velocity is as follows:
S1, constructing a rock equivalent volume model of a reservoir, wherein the rock is pore-fluid-saturated rock.
Specifically, as shown in fig. 2, the rock equivalent volume model of the reservoir defines that the sectional area of the pore-saturated fluid rock is a, the length of the pore-saturated fluid rock is L, the length of the equivalent throat is L w, and the sectional area of the equivalent throat is a w.
S2, representing the porosity and the flat rate of the rock in the reservoir by the resistivity, the porosity and the pore curvature of the rock saturated with the pore fluid based on the rock equivalent volume model.
For pure rock, the rock framework may be considered to be substantially non-conductive, with only the fluid in the rock pores being conductive. Whereas the rock openings are curved, the current also flows in a meandering manner in the rock. Therefore, the embodiment of the invention simplifies the rock equivalent volume model according to the current flowing condition.
Specifically, the resistivity and pore tortuosity of the supersaturated pore fluid rock based on the rock equivalent volume model are used for characterizing the pore flatness of the rock in the reservoir, and the method comprises the following steps:
S21, constructing a first relation equation of rock resistivity saturated with pore fluid according to a resistance parallel calculation formula based on a rock equivalent volume model, wherein the first relation equation is specifically as follows:
Where r 0、rma and r w are water-saturated pure rock resistance, rock matrix resistance, and pore fluid resistance, respectively.
S22, determining a second relation between the resistivity of the saturated pure rock and the resistivity of the pore fluid according to the first relation.
In general, it can be considered that the resistance of the rock skeleton tends to infinity, and the first relation is solved according to the relation between the resistivity and the resistance, so as to obtain a second relation, which is specifically as follows:
I.e.
Where R o represents the electrical resistance of the rock and R w represents the pore fluid resistance.
S23, defining the pore canal curvature of the rock, and determining the pore flatness of the rock in the reservoir according to the porosity of the rock and the second relation.
Wherein the curvature of the pore canal of the rock is thatI.e. the extent of bending of the aperture throat.
Wherein the porosity of the rockThe method comprises the following steps:
and determining the porosity flat rate of the rock in the reservoir according to the porosity of the rock and a second relation, wherein the porosity flat rate of the rock in the reservoir is specifically as follows:
S231, determining the relation between the pore canal curvature of the rock, the rock resistance, the pore fluid resistance and the porosity of the rock according to the porosity of the rock and the second relation, wherein the relation is specifically as follows:
s232, obtaining logarithm of pore curvature to obtain porosity alpha of the rock, wherein the pore flatness alpha is specifically as follows:
Through calculation, the porosity and the flat rate in sandstone are about 0.1 and 0.02-0.05 respectively.
S3, establishing a first expression of the elastic modulus of the rock framework by using a K-T model, wherein the elastic modulus comprises a bulk modulus and a shear modulus.
In the step, assuming that the Poisson's ratio of the rock does not change along with the porosity, a K-T equation based on a DEM theory can transform the problem of a linear ordinary differential equation set to obtain a first expression of the elastic modulus of the rock skeleton, which is specifically as follows:
In the method, in the process of the invention, And/>When the porosity is/>, respectivelyWhen the rock framework has bulk modulus and shear modulus; k m and G m are when the porosity is/>, respectivelyBulk and shear moduli of the rock matrix when; p and q are a set of coefficients related only to the porosity alpha and not to the porosity when the poisson ratio of the rock framework is constant.
And S4, establishing a second expression according to the relation between the correlation constant of the elastic modulus of the rock framework and the porosity and the flat rate.
In the step, introducing an Eshelby tensor, establishing a second expression according to the relation between the correlation constant of the elastic modulus of the rock skeleton and the porosity and the flat rate, wherein the second expression is specifically as follows:
wherein, T iijj, F are functions of pore flat rate respectively, and are defined in Eshelby tensor T ijkl; s, SH represents sandstone and shale, respectively; v l represents the volume percent of sandstone and shale to the rock matrix and α l represents the pore volume.
S5, according to the first expression and the second expression, the elastic modulus of the rock framework and the rock matrix is obtained, and the method specifically comprises the following steps:
And S51, calculating the relation between the bulk modulus and the shear modulus of the rock skeleton and the bulk modulus and the shear modulus of the rock matrix through a Krief model and a prism model.
S52, calculating to obtain the bulk modulus and the shear modulus of the rock skeleton and the bulk modulus and the shear modulus of the rock matrix according to the first expression and the second expression based on the calculation result of the relation between the bulk modulus and the shear modulus of the rock skeleton and the rock matrix.
S6, calculating the equivalent bulk modulus and the equivalent shear modulus of the pore-saturated fluid rock in the reservoir according to the rock skeleton and the rock matrix elastic modulus by using a Gassmann equation, wherein the equivalent bulk modulus and the equivalent shear modulus are specifically as follows:
G＝Gd
Wherein K is the equivalent bulk modulus of the pore-saturated fluid rock; g is the equivalent shear modulus of the pore-fluid-saturated rock; k d is the bulk modulus of the rock framework; k m is the bulk modulus of the rock matrix; k f is the bulk modulus of the pore fluid; g d is the shear modulus of the rock backbone, where K f＝LwAw.
S7, predicting the transverse wave speed of the rock in the reservoir according to the equivalent shear modulus and the equivalent density of the rock saturated with the pore fluid by using a Gassmann equation.
Where ρ is the equivalent density of the pore-saturated fluid rock; v s is the shear wave velocity of the pore-fluid-saturated rock.
Wherein ρ is the equivalent density of the pore-saturated fluid rock as follows:
Where ρ f is the density of the pore fluid, ρ d is the density of the rock framework, and ρ f、ρd can be calculated by willie on the basis of consulting reservoir fluid parameters and mineral density.
S8, predicting the longitudinal wave velocity of the rock in the reservoir according to the equivalent bulk modulus, the equivalent shear modulus and the equivalent density of the pore-saturated fluid rock by using a Gassmann equation, wherein the method comprises the following steps of:
Where v p is the longitudinal wave velocity of the pore-fluid-saturated rock.
As shown in fig. 3, based on the method for predicting the transverse wave based on the curvature of the duct, the embodiment of the invention also provides a transverse wave prediction system based on the curvature of the duct, which comprises a model building module, a first calculation module, a second calculation module, a third calculation module and a transverse wave prediction module.
The model building module is used for building a rock equivalent volume model of the reservoir, wherein the rock is pore-fluid-saturated rock; a first calculation module for characterizing a pore flatness of the rock in the reservoir by resistivity, porosity and pore tortuosity of the pore-saturated fluid rock based on the rock equivalent volume model; the second calculation module is used for calculating and obtaining the elastic modulus of the rock skeleton and the rock matrix according to the first expression of the elastic modulus of the rock skeleton and the second expression of the related constant of the elastic modulus of the rock skeleton and the porosity; the third calculation module is used for calculating the equivalent bulk modulus and the equivalent shear modulus of the pore-saturated fluid rock in the reservoir according to the elastic modulus of the rock skeleton and the rock matrix; and the transverse wave prediction module is used for predicting the transverse wave speed of the rock in the reservoir according to the equivalent shear modulus and the equivalent density of the pore-saturated fluid rock.
As shown in fig. 4, compared with the conventional method prediction result, the conventional method prediction result and the actually measured transverse wave curve can only be matched to 90%, while the transverse wave prediction result and the actually measured transverse wave curve of the embodiment of the invention can be matched to 96%, which indicates that the transverse wave prediction result of the embodiment of the invention is more accurate and achieves the expected effect.
As shown in fig. 5 and fig. 6, compared with the error of the transverse wave prediction result of the embodiment of the present invention and the transverse plate prediction result of the conventional method, the transverse wave speed prediction value and the actual measurement value of the embodiment of the present invention are closer and more accurate.
The rock equivalent volume model provided by the embodiment of the invention can better predict the longitudinal wave speed and the transverse wave speed of the rock full of the pore fluid, reduces the requirements on the input precision of porosity, mineral components and the like in the use process of the traditional Xu-White model, has strong practicability, and provides strong technical support for the mechanical evaluation and the seismic inversion of the stratum rock lacking transverse wave logging data.
Although the invention has been described in detail with reference to the foregoing embodiments, it will be understood by those of ordinary skill in the art that: the technical scheme described in the foregoing embodiments can be modified or some technical features thereof can be replaced by equivalents; such modifications and substitutions do not depart from the spirit and scope of the technical solutions of the embodiments of the present invention.
Claims (11)
1. The transverse wave prediction method based on the curvature of the pore canal is characterized by comprising the following steps of:
constructing a rock equivalent volume model of the reservoir, wherein the rock is pore-fluid-saturated rock;
Characterizing the pore flatness of the rock in the reservoir by resistivity, porosity and pore tortuosity of the pore-saturated fluid rock based on the rock equivalent volume model;
According to a first expression of the elastic modulus of the rock skeleton and a second expression of the correlation constant of the elastic modulus of the rock skeleton and the porosity, calculating to obtain the elastic modulus of the rock skeleton and the rock matrix;
calculating the equivalent bulk modulus and the equivalent shear modulus of the pore-saturated fluid rock in the reservoir according to the rock skeleton and the rock matrix elastic modulus;
And predicting the shear wave velocity of the rock in the reservoir according to the equivalent shear modulus and the equivalent density of the pore-saturated fluid rock.
2. The method of claim 1, wherein the constructing of the pore-saturated fluid rock equivalent volume model of the reservoir is performed by:
The cross-sectional area of the saturated pore fluid rock is defined as A, the length of the saturated pore fluid rock is defined as L, the length of the equivalent throat of the saturated pore fluid rock is defined as L w, and the cross-sectional area of the equivalent throat of the saturated pore fluid rock is defined as A w.
3. The method of channel tortuosity based shear wave prediction according to claim 2 wherein the pore flatness of the rock in the reservoir is characterized by the resistivity, porosity and channel tortuosity of the pore-saturated fluid rock based on the rock equivalent volume model is specified as follows:
Based on the rock equivalent volume model, constructing a first relation of rock resistivity saturated with pore fluid according to a resistance parallel calculation formula;
determining a second relationship between the saturated water-free rock resistivity and the pore fluid resistivity from the first relationship;
and defining the pore canal curvature of the rock, and determining the porosity and the flat rate of the rock in the reservoir according to the porosity of the rock and the second relation.
4. A method of channel tortuosity based shear wave prediction according to claim 3 wherein the first relationship is specifically:
Where r 0、rma and r w are water-saturated pure rock resistance, rock matrix resistance, and pore fluid resistance, respectively.
5. The method of claim 4, wherein the second relation between the resistivity of saturated pure rock and the resistivity of pore fluid is determined according to the first relation, and is specifically as follows:
defining the rock skeleton resistance to be infinite, solving a first relation according to the relation between the resistivity and the resistance to obtain a second relation, wherein the specific relation is as follows:
Where R o represents the electrical resistance of the rock, R w represents the pore fluid resistance, A represents the cross-sectional area of the saturated pore fluid rock, L represents the length of the saturated pore fluid rock, L w represents the length of the equivalent throat of the saturated pore fluid rock, and A w represents the cross-sectional area of the equivalent throat of the saturated pore fluid rock.
6. The method for predicting shear waves based on tunnel curvature according to claim 5, wherein the tunnel curvature of the rock isPorosity of rock/>The method comprises the following steps:
where A represents the cross-sectional area of the saturated pore fluid rock, L represents the length of the saturated pore fluid rock, L w represents the length of the equivalent throat of the saturated pore fluid rock, and A w represents the cross-sectional area of the equivalent throat of the saturated pore fluid rock.
7. The method of claim 6, wherein the porosity of the rock in the reservoir is determined based on the porosity of the rock and a second relationship, and wherein the method comprises:
according to the porosity of the rock and the second relation, the relation between the pore canal curvature of the rock, the rock resistance, the pore fluid resistance and the porosity of the rock is determined, and the relation is specifically as follows:
the logarithm of the curvature of the pore canal is obtained, and the porosity flat rate alpha of the rock is obtained, specifically as follows:
Where R o represents the electrical resistance of the rock, R w represents the pore fluid resistance, L represents the length of the pore fluid-saturated rock, L w represents the length of the equivalent throat of the pore fluid-saturated rock, Representing the porosity of the rock.
8. The bending-based shear wave prediction method according to any of claims 1-7, further comprising the steps of:
and predicting the longitudinal wave velocity of the rock in the reservoir according to the equivalent bulk modulus, the equivalent shear modulus and the equivalent density of the pore-saturated fluid rock.
9. A transverse wave prediction system based on tunnel curvature, comprising:
The model building module is used for building a rock equivalent volume model of the reservoir, wherein the rock is pore-fluid-saturated rock;
a first calculation module for characterizing a pore flatness of the rock in the reservoir by resistivity, porosity and pore tortuosity of the pore-saturated fluid rock based on the rock equivalent volume model;
the second calculation module is used for calculating and obtaining the elastic modulus of the rock skeleton and the rock matrix according to the first expression of the elastic modulus of the rock skeleton and the second expression of the related constant of the elastic modulus of the rock skeleton and the porosity;
the third calculation module is used for calculating the equivalent bulk modulus and the equivalent shear modulus of the pore-saturated fluid rock in the reservoir according to the elastic modulus of the rock skeleton and the rock matrix;
And the transverse wave prediction module is used for predicting the transverse wave speed of the rock in the reservoir according to the equivalent shear modulus and the equivalent density of the pore-saturated fluid rock.
10. The system for channel tortuosity-based shear wave prediction of claim 9 wherein the model building module is specifically configured to:
The cross-sectional area of the saturated pore fluid rock is defined as A, the length of the saturated pore fluid rock is defined as L, the length of the equivalent throat of the saturated pore fluid rock is defined as L w, and the cross-sectional area of the equivalent throat of the saturated pore fluid rock is defined as A w.
11. The system of claim 10, wherein the first computing module is configured to:
Based on the rock equivalent volume model, constructing a first relation of rock resistivity saturated with pore fluid according to a resistance parallel calculation formula;
determining a second relationship between the saturated water-free rock resistivity and the pore fluid resistivity from the first relationship;
and defining the pore canal curvature of the rock, and determining the porosity and the flat rate of the rock in the reservoir according to the porosity of the rock and the second relation.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310548829.4A CN116819607B (en) | 2023-05-16 | 2023-05-16 | Transverse wave prediction method and system based on tunnel curvature |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310548829.4A CN116819607B (en) | 2023-05-16 | 2023-05-16 | Transverse wave prediction method and system based on tunnel curvature |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN116819607A CN116819607A (en) | 2023-09-29 |
| CN116819607B true CN116819607B (en) | 2024-04-19 |
Family
ID=88117508
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202310548829.4A Active CN116819607B (en) | 2023-05-16 | 2023-05-16 | Transverse wave prediction method and system based on tunnel curvature |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN116819607B (en) |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2013149623A1 (en) * | 2012-04-01 | 2013-10-10 | Entreprise Nationale De Geophysique Enageo | Method for quantitatively evaluating the fluid tortuosity and the characteristics of the solid and of the fluids in a heterogeneous reservoir |
| CN103954999A (en) * | 2014-05-06 | 2014-07-30 | 成都理工大学 | Novel shear wave velocity prediction method suitable for low-porosity sand shale stratum |
| CN106202879A (en) * | 2016-06-29 | 2016-12-07 | 中国石油化工股份有限公司 | P-and s-wave velocity dynamic estimation method during carbon dioxide drive |
| CN115793048A (en) * | 2022-11-16 | 2023-03-14 | 河海大学 | Shale oil bed system compact reservoir transverse wave velocity prediction method |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| AU2009234101B2 (en) * | 2008-04-09 | 2014-01-09 | Exxonmobil Upstream Research Company | Method for generating anisotropic resistivity volumes from seismic and log data using a rock physics model |
-
2023
- 2023-05-16 CN CN202310548829.4A patent/CN116819607B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2013149623A1 (en) * | 2012-04-01 | 2013-10-10 | Entreprise Nationale De Geophysique Enageo | Method for quantitatively evaluating the fluid tortuosity and the characteristics of the solid and of the fluids in a heterogeneous reservoir |
| CN103954999A (en) * | 2014-05-06 | 2014-07-30 | 成都理工大学 | Novel shear wave velocity prediction method suitable for low-porosity sand shale stratum |
| CN106202879A (en) * | 2016-06-29 | 2016-12-07 | 中国石油化工股份有限公司 | P-and s-wave velocity dynamic estimation method during carbon dioxide drive |
| CN115793048A (en) * | 2022-11-16 | 2023-03-14 | 河海大学 | Shale oil bed system compact reservoir transverse wave velocity prediction method |
Non-Patent Citations (2)
| Title |
|---|
| 基于孔隙分类理论的自相容模型横波速度预测方法;熊晓军;李翔;刘阳;简世凯;;石油物探(第02期);全文 * |
| 测井资料Xu-White模型预测横波速度的一些新观点;石双虎;邓志文;白光宇;李红星;蔡敏贵;;地震工程学报(第04期);全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN116819607A (en) | 2023-09-29 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Maddux et al. | Turbulent flow over three‐dimensional dunes: 2. Fluid and bed stresses | |
| US20200211126A1 (en) | Prediction method for shale oil and gas sweet spot region, computer device and computer readable storage medium | |
| CN110954948B (en) | Physical property parameter inversion method and system for rock physical constraint reservoir | |
| CN104252007B (en) | A kind of compatibility rock physicses modeling method | |
| CN103954999B (en) | A kind of S-Wave Velocity Predicted Method being applicable to low-porosity sand-mudstone formation | |
| Chen et al. | Porosity prediction from well logs using back propagation neural network optimized by genetic algorithm in one heterogeneous oil reservoirs of Ordos Basin, China | |
| CN106680878A (en) | Transverse wave speed estimation method based on improved Biot coefficient analysis | |
| CN110133720A (en) | A kind of S-Wave Velocity Predicted Method based on statistics petrophysical model | |
| CN110426751A (en) | Method for predicting transverse wave time difference by using logging data | |
| Guo et al. | Tide-induced groundwater head fluctuation in coastal multi-layered aquifer systems with a submarine outlet-capping | |
| CN115793048B (en) | Shale oil reservoir system tight reservoir shear wave speed prediction method | |
| CN109581490A (en) | A kind of Wave velocity estimation method and system of compact reservoir | |
| Miah et al. | Model development for shear sonic velocity using geophysical log data: Sensitivity analysis and statistical assessment | |
| CN116819607B (en) | Transverse wave prediction method and system based on tunnel curvature | |
| Wang et al. | Lithology classification and porosity estimation of tight gas reservoirs with well logs based on an equivalent multi-component model | |
| Carling et al. | Turbulent flow across a natural compound channel | |
| Izumotani et al. | Estimating rock properties using seismic refraction survey data: a case study in an abandoned road tunnel | |
| CN115685335A (en) | Longitudinal and transverse wave velocity prediction method and device, computer equipment and storage medium | |
| CN114428372B (en) | Self-adaptive rock physical modeling method | |
| Bashmakov et al. | Some features of fluid filtration in a hydraulic fracture under transient well operation conditions | |
| CN109490988A (en) | Establish the method for being suitable for the rock physics new model of hard rock | |
| AlKawai et al. | Integrating statistical rock physics and pressure and thermal history modeling to map reservoir lithofacies in the deepwater Gulf of Mexico | |
| Weijun et al. | Prediction of shear wave velocity in shale reservoir based on logging data and machine learning | |
| Alhumoud | Analysis and evaluation of flood routing using Muskingum method | |
| Liungman et al. | Modeling and observations of deep water renewal and entrainment in a Swedish sill fjord |
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 | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |