CN113887150B - Method for estimating length of characteristic jet flow of compact sandstone - Google Patents
Method for estimating length of characteristic jet flow of compact sandstone Download PDFInfo
- Publication number
- CN113887150B CN113887150B CN202111101795.1A CN202111101795A CN113887150B CN 113887150 B CN113887150 B CN 113887150B CN 202111101795 A CN202111101795 A CN 202111101795A CN 113887150 B CN113887150 B CN 113887150B
- Authority
- CN
- China
- Prior art keywords
- rock
- fracture
- pressure
- characteristic
- length
- 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 37
- 239000011435 rock Substances 0.000 claims abstract description 115
- 230000035699 permeability Effects 0.000 claims abstract description 20
- 230000014509 gene expression Effects 0.000 claims abstract description 19
- 238000005259 measurement Methods 0.000 claims abstract description 11
- 239000011148 porous material Substances 0.000 claims description 80
- 239000012530 fluid Substances 0.000 claims description 54
- 239000011159 matrix material Substances 0.000 claims description 15
- 238000009792 diffusion process Methods 0.000 claims description 13
- 230000008859 change Effects 0.000 claims description 12
- 239000007787 solid Substances 0.000 claims description 12
- 238000012360 testing method Methods 0.000 claims description 11
- 238000006073 displacement reaction Methods 0.000 claims description 10
- 238000002474 experimental method Methods 0.000 claims description 10
- 239000007788 liquid Substances 0.000 claims description 10
- 239000002245 particle Substances 0.000 claims description 9
- 229920006395 saturated elastomer Polymers 0.000 claims description 9
- 230000001186 cumulative effect Effects 0.000 claims description 5
- 230000008878 coupling Effects 0.000 claims description 4
- 238000010168 coupling process Methods 0.000 claims description 4
- 238000005859 coupling reaction Methods 0.000 claims description 4
- 230000006355 external stress Effects 0.000 claims description 4
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 claims description 3
- 230000003993 interaction Effects 0.000 claims description 3
- 230000000704 physical effect Effects 0.000 claims description 3
- 230000035882 stress Effects 0.000 claims description 3
- 238000000691 measurement method Methods 0.000 claims description 2
- 239000004576 sand Substances 0.000 claims 1
- 238000011160 research Methods 0.000 abstract description 6
- 206010017076 Fracture Diseases 0.000 description 126
- 208000010392 Bone Fractures Diseases 0.000 description 106
- 239000006185 dispersion Substances 0.000 description 5
- 230000000694 effects Effects 0.000 description 5
- 230000006872 improvement Effects 0.000 description 4
- 238000012821 model calculation Methods 0.000 description 3
- 239000010453 quartz Substances 0.000 description 3
- VYPSYNLAJGMNEJ-UHFFFAOYSA-N silicon dioxide Inorganic materials O=[Si]=O VYPSYNLAJGMNEJ-UHFFFAOYSA-N 0.000 description 3
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 2
- 230000007423 decrease Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000001035 drying Methods 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 238000006677 Appel reaction Methods 0.000 description 1
- 102100021385 H/ACA ribonucleoprotein complex subunit 1 Human genes 0.000 description 1
- 101000771075 Homo sapiens Cyclic nucleotide-gated cation channel beta-1 Proteins 0.000 description 1
- 101000819109 Homo sapiens H/ACA ribonucleoprotein complex subunit 1 Proteins 0.000 description 1
- 240000007594 Oryza sativa Species 0.000 description 1
- 235000007164 Oryza sativa Nutrition 0.000 description 1
- 208000013201 Stress fracture Diseases 0.000 description 1
- 238000009825 accumulation Methods 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 239000004568 cement Substances 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 229920001971 elastomer Polymers 0.000 description 1
- 238000005485 electric heating Methods 0.000 description 1
- 239000010433 feldspar Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000005498 polishing Methods 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 235000009566 rice Nutrition 0.000 description 1
- 238000009738 saturating Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/02—Generating seismic energy
- G01V1/133—Generating seismic energy using fluidic driving means, e.g. highly pressurised fluids; using implosion
-
- 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. for interpretation or for event detection
- G01V1/30—Analysis
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Remote Sensing (AREA)
- Theoretical Computer Science (AREA)
- Geophysics (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geology (AREA)
- Environmental & Geological Engineering (AREA)
- Acoustics & Sound (AREA)
- Mathematical Analysis (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Physics (AREA)
- Mathematical Optimization (AREA)
- Fluid Mechanics (AREA)
- Computing Systems (AREA)
- Algebra (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
Abstract
The invention discloses a method for estimating the length of a characteristic jet flow of tight sandstone. The research on the distribution characteristics of the cracks in the rock and the characteristic jet length characteristics is carried out, the influence of parameters such as permeability and crack density on the characteristic jet length is analyzed, and finally, a quantitative expression of the characteristic jet length is obtained, so that an effective method is provided for determining the characteristic jet length of the compact sandstone. The invention fully considers the relation between the characteristic jet flow length of the compact sandstone and the fracture distribution characteristic in the rock and forms a method for determining the characteristic jet flow length. The method characterizes the characteristic jet flow length through measurable parameters, so that the characteristic jet flow length is converted to be acquired through an indirect measurement mode.
Description
Technical Field
The invention relates to the field of seismic rock physics, in particular to a method for estimating the characteristic jet flow length of tight sandstone.
Background
The tight sandstone pore permeability is low, the micro-fracture in the rock develops, strong heterogeneity is presented, and the difficulty of reservoir prediction and fluid identification is greatly increased. The micro-pore structure of the reservoir is not only closely related to the rock physical characteristics, but also is a main factor influencing the oil and gas distribution of the reservoir. Therefore, the research on the characteristics of the micro-pore structure can not only define the basic characteristics and the reservoir performance of the reservoir, but also explain the basis of the reservoir oil-gas enrichment rule (white bin, etc., 2013). The presence of fissures inside the rock has an effect not only on the elasticity of the rock, but also on the flow of fluids in the rock (Carcione, 2007; Muller et al, 2010). In order to study the influence of fractures on the elastic properties of rocks, researchers studied equivalent medium theory and petrophysical experiments (Song et al, 2016; Yin et al, 2017). Based on the method, David and Zimmmerman (2012) predict the characteristics of pore distribution (DZ model for short) corresponding to fractures with different aspect ratios from the ultrasonic velocity of dry rock by using an MT model (Mori and Tanaka,1973) and a DEM model (Berrayman,1992), and analyze the rule that the characteristic of the longitudinal and transverse waves of saturated rock changes along with the effective pressure.
In order to study the wave response characteristics of seismic waves as they pass through complex pore media, researchers have proposed a series of theoretical models (Ba et al, 2011; 2017; Muller et al, 2010; Zhang et al, 2019; Zhang et al, 2020). With the progress of theoretical research and the accumulation of data measured in the field, it is gradually recognized that local flow caused by heterogeneity is the main cause of elastic wave attenuation frequency dispersion (yankui et al, 2000). The BISQ model (Dvorkin,1993) is based on the Biot model (Biot,1956a,1956b,1962), a new petrophysical parameter 'characteristic jet length' is introduced to describe a jet flow (Squirt flow), and the combination of macro scale and micro scale is realized. The fast longitudinal wave velocity frequency dispersion predicted by the model is obviously larger than the calculation result of the Biot model. However, the fast longitudinal wave velocity estimated by the model is smaller than the predicted value of the Gassmann model (Gassmann,1951) at lower frequencies, and is consistent with the prediction result of the Biot model at higher frequencies, which are not consistent with the physical fact. Mavko and Jizba (1991) propose a jet relation (M-J model) for calculating the elastic modulus of a high frequency unrelaxed wet rock skeleton based on a modified rock skeleton. The M-J model predicts a higher frequency compressional velocity value than that calculated by the Biot model, but the model predicts a compressional velocity value greater than that predicted by the Gassmann model at lower frequencies, the jet characteristic frequency is consistent with that of the Biot flow, and the formula is not applicable to the case where the pores contain gas or dry rock. Aiming at the limitation of the M-J model, Gurevich et al (2010) applies a pore structure model proposed by Murphy et al (1986), considers the pores between the contact surfaces of adjacent particles as soft pores, generalizes the M-J model to be suitable for the gas saturation condition and establishes a jet flow model with intermediate frequency dispersion, so that the model is expanded to the whole frequency range. Dvorkin et al (1995) rebuilds a complex modulus model of the fluid saturated rock jet based on the BISQ model principle in combination with a modified rock skeleton jet theoretical model, and finds an expression of the complex elastic modulus by considering one-dimensional radial flow. The predicted longitudinal wave velocity of the deduced complex modulus at the low-frequency limit is consistent with the result calculated by the Gassmann theory. However, this model does not consider the effect of Biot flow on the propagation characteristics of elastic waves, and the longitudinal wave velocity calculated by this model is larger than the theoretical maximum value at high frequencies (the longitudinal wave velocity calculated by the Biot model when all rock fractures are closed) (Wu et al, 2020).
Dvorkin et al (1993) first introduced a characteristic jet length to study the propagation characteristics of elastic waves in a reservoir, a parameter that cannot be obtained by experimental measurements and is on the order of magnitude of the rock particle size. Dvorkin et al (1994) then gives the physical meaning of the characteristic jet length, geometrically, this parameter represents the radius of the cylinder axially aligned with the wave propagation axis. Marketons and Best (2010) rely on ultrasonic experimental data of Best (1992) to study that a characteristic jet flow length has a certain power series relation with the viscosity of liquid. After the BISQ model was proposed, the researchers performed a lot of research work on the model, but focused on improving the BISQ model (Diallo and Appel., 2000; Diallo et al, 2003; Dvorkin, 1995; Wu et al, 2019; Wu et al,2020), jet anisotropy (Yang and Zhang,2000,2002) and wave field simulation (Yanwidend et al, 2002; Shenqing and Yankee, 2004), etc., without proposing a method for determining the characteristic jet length.
Disclosure of Invention
Aiming at the fact that the characteristic jet flow length is an important parameter of a relevant jet flow theoretical model, an effective establishing method is not formed at present, the invention provides a method for estimating the characteristic jet flow length of the compact sandstone, aiming at the compact sandstone, a fracture distribution parameter and a characteristic jet flow length value are calculated, the relation between the fracture distribution parameter and the characteristic jet flow length value is analyzed, and a quantitative relation formula of the characteristic jet flow length is established.
In order to solve the problems of the prior art, the invention adopts the technical scheme that:
a method of estimating a tight sandstone characteristic jet length, comprising the steps of:
and 7, selecting the main fracture aspect ratio, permeability, matrix modulus, fluid volume modulus and fracture density to characterize the characteristic jet flow length to form a characteristic jet flow length quantitative relation, wherein the characteristic jet flow length quantitative relation of the compact sandstone is as follows:
and fitting values of the parameters a and b according to the measurement result of the ultrasonic experiment, wherein the values have certain difference for different rock samples.
The improvement is that step 1 obtains ultrasonic velocity experimental data of pressure change of the tight sandstone, and the steps are as follows: and testing physical property parameters of the rock of the compact sandstone sample under the pressure change in a dry state and a water-saturated state by using an ultrasonic experimental measurement method, and calculating longitudinal wave speed and transverse wave speed in the dry and water-saturated states.
The improvement is that the step 2 of estimating the fracture density of the tight sandstone comprises the following steps:
the model establishes a quantitative relationship between the elastic modulus and the hard pores in the rock, and the equivalent elastic modulus can be expressed as:
wherein,andthe equivalent bulk modulus and the equivalent shear modulus of the rock, K, are respectively when only the hard pores are contained0And mu0Respectively the bulk modulus and shear modulus of the rock particles, phisIs the hard porosity of the rock, and M and N are the shape factors of the hard porosity, defined as follows:
wherein the shape of the hard pores is default to elliptical pores, γ is the aspect ratio of the elliptical pores, upsilon is the poisson ratio of the rock particles, and upsilon ═ 3K0-2μ0)/(6K0+2μ0) In the formula, the expression of g is as follows:
if the fractures are added to a pore medium containing hard pores, the content of the fractures is small, no interaction exists between the fractures and the pores, and the equivalent elastic modulus of the rock can be expressed as:
wherein,andrespectively the equivalent bulk modulus and the equivalent shear modulus when the rock contains cracks and hard pores,is the rock poisson's ratio with only hard pores, Γ is the cumulative fracture density; because the fracture is approximately and completely closed under higher effective pressure, and the interior of the rock only contains hard pores, the optimal aspect ratio of the hard pores can be obtained by a formula (1) and a least square method;
the pressure dependence of the effective bulk modulus of the rock is related to the fracture density, so when the fracture density is known, the elastic modulus of the rock can be obtained through the formulas (5) and (6), and conversely, the fracture density at each pressure can be calculated through the elastic modulus at each effective pressure by using a least square method.
As an improvement, the step 3 of estimating the aspect ratio of the main fractures of the tight sandstone comprises the following steps:
quantitative relationship of cumulative fracture density with effective pressure change
In the formula, gammaiIs the initial fracture density at zero pressure, p is the pressure,is a pressure constant of the same order of magnitude as the pressure;
dividing the pressure interval delta p of the dry rock test into n parts, and solving the fracture density of each pressure by using a formula (7):
wherein k is more than or equal to 0 and less than or equal to n;
dividing the pressure interval of the dry rock test at equal intervals delta p into n parts, wherein each part is delta p, and the fracture aspect ratio at each pressure point is as follows:
wherein,is the effective poisson's ratio at high pressure,is the effective young's modulus at high pressure;
the relationship between fracture porosity and fracture density is:
and (3) solving the fracture aspect ratio corresponding to the maximum fracture density by using the formula (10), wherein the fracture aspect ratio is the main fracture aspect ratio.
Preferably, the step of deriving the MFS model at step 4 is as follows:
assuming an ideal cylinder of radius R, the relationship between the deformation of the framework (e ═ du/dx) and the fracture and pore pressures is as follows:
wherein u is the solid displacement, c0Is the acoustic velocity, p, of the liquidfIs the fluid density, Qc=K0/(1-Kmsd/K0-φc),Kmsd=(1/K0-1/Khp+1/Kd)-1,KdIs the bulk modulus, phi, of the dry rockcIs the fracture porosity, KhpIs the bulk modulus of the dry rock when the fracture is fully closed;
assuming that the liquid inside the fracture can only flow in the direction perpendicular to the wave propagation direction, the fluid conservation of mass equation is:
wherein t is time and w is displacement of the liquid;
substituting the equations (11) and (12) into the equation (13) can obtain:
the relationship between fluid velocity and pore pressure gradient can be derived:
where η is the fluid viscosity and κ is the permeability.
The simple harmonic expression is:
exp (-i ω t) (16) substituting equation (16) into equations (14) and (15) can be obtained:
wherein, KfIs the bulk modulus of the fluid;
from equations (17) and (18), the differential equation of the pressure is derived as follows:
at the fracture boundary, the pore pressure is always constant, so at R, p, dp, 0, the result of solving equation (19) is:
wherein R is the characteristic jet length, λ2=iωηφc/κ(1/Kf+1/(φcQc)),J0Is a zero order bessel function;
the average fluid pressure is:
wherein, J1Is a first order Bessel function;
in the modified rock, the external stress and displacement and pore mean pressure are as follows:
dσ=Kmsdde-(1-Kmsd/K0)dpav (22)
modified bulk modulus of solids Kms:
Substituting equations (21) and (22) into equation (23) yields:
the following relationship exists between the dry pore volume, matrix modulus and dry rock modulus according to the definition of bulk modulus, and the dry bulk modulus of the modified rock can be derived.
The expression for modifying the dry shear modulus of the rock is as follows:
wherein, mudIs the shear modulus of dry rock according toand Johnston,1981) wave number versus wave velocity and attenuation, and deriving expressions for longitudinal wave velocity and attenuation factor.
ρaFor fluid-solid coupling density, ρsIs the density of the matrix and is,for a characteristic frequency, MdIs the uniaxial stress modulus of the rock skeleton,
as a refinement, the step 5 of calculating the characteristic jet length using MFS comprises the steps of: and (4) calculating the characteristic jet flow lengths of different pressures by using the MFS model obtained in the step 4, and preferably selecting the most appropriate characteristic jet flow length of each sample by using a least square method.
As an improvement, the source steps of the compact sandstone characteristic jet flow length quantitative relation formula in the step 7 are as follows:
the critical diffusion length of the fluid is given according to the characteristic frequency of the jet flow as follows:
from the equations (31) and (32), the critical diffusion length of the fluid is related to the fracture aspect ratio, permeability and other parameters, and is expressed as:
when the fluid diffusion length is less than the critical diffusion length, the pore fluid pressure has sufficient time to equilibrate with the fluid flow, whereas the pore fluid pressure has insufficient time to equilibrate.
The method is characterized in that the main fracture aspect ratio, permeability, matrix modulus, fluid volume modulus and fracture density are selected to represent the characteristic jet flow length, a characteristic jet flow length quantitative relation is formed, and the dense sandstone characteristic jet flow length quantitative relation is as follows:
and fitting values of the parameters a and b according to the measurement result of the ultrasonic experiment, wherein the values have certain difference for different rock samples.
Has the advantages that:
compared with the existing theory, the method for estimating the characteristic jet flow length of the tight sandstone combines the ultrasonic experimental data of the tight sandstone, adopts a DZ model to analyze the fracture distribution characteristics in the tight sandstone rock, adopts an MFS model to calculate the characteristic jet flow length, determines the influence factors of the characteristic jet flow length, researches the relationship between the characteristic jet flow length and the fracture distribution, mainly analyzes the relationship between the characteristic jet flow length and parameters such as fracture density, fracture aspect ratio, permeability and the like, and finally forms a quantitative relational expression of the characteristic jet flow length of the tight sandstone. The invention provides a certain basis for establishing the characteristic jet flow length of the compact sandstone, and the specific effects are as follows:
a) estimating fracture distribution parameters in the rock based on the ultrasonic experimental data of the compact sandstone with pressure change;
b) calculating the characteristic jet flow length by combining an MFS model;
c) and analyzing the relation between the characteristic jet flow length and the fracture parameters through compact sandstone experimental data, and providing an expression of the characteristic jet flow length.
Drawings
FIG. 1 is a schematic flow diagram of a method of estimating a characteristic jet length of tight sandstone in accordance with the present invention;
FIG. 2 is the longitudinal and transverse wave velocity ratio of dry and water-saturated state of 12 pieces of tight sandstone;
FIG. 3 is a graph of the relationship between the density and pressure of 12 tight sandstone fractures;
FIG. 4 is a plot of fracture density versus most dominant fracture aspect ratio for sample GAR 7;
FIG. 5 is a diagram of a modified rock skeleton;
FIG. 6 is the predicted compressional velocity of the MFS model;
FIG. 7 is a graph of model calculated values versus experimental values for sample GAR7 at various pressures;
FIG. 8 is a graph of the results of comparing the predicted compressional velocity of different characteristic jet lengths of GAR7 with experimental values;
FIG. 9 is a graph of characteristic jet length versus various parameters;
FIG. 10 is a plot of characteristic jet length versus fracture density, major fracture aspect ratio, and permeability.
Detailed Description
The invention is further described below with reference to the accompanying drawings. The following examples are only for illustrating the technical solutions of the present invention more clearly, and the protection scope of the present invention is not limited thereby.
The method aims at the fact that the characteristic jet flow length is an important factor of a jet flow theoretical model, firstly, fracture parameters (fracture density, fracture porosity and main fracture aspect ratio) are calculated by using a DZ model, secondly, a derived MFS model is adopted, the characteristic jet flow length is calculated by using the model, and finally, influence factors of the characteristic jet flow length are analyzed, and an expression for representing the characteristic jet flow length through the fracture parameters is formed by combining the physical significance of the characteristic jet flow length.
A method for estimating the characteristic jet flow length of tight sandstone comprises the following specific steps as shown in figure 1:
The invention aims at a work area of a Sanhua junjia river group on a Guangan gas field in the middle of a Sichuan basin, and 12 compact sandstone samples are selected in the work area. The compact sandstone sample is mainly rock debris quartz sandstone and a small amount of siliceous quartz sandstone, and consists of quartz, feldspar, rock debris, cement and the like, wherein the porosity is mainly distributed in a range of 3-14%, and the modulus of the compact sandstone matrix is 39 GPa.
And (3) processing the compact rock sample into cylinders with the diameters of 25mm and the heights of 25-50mm, and polishing and grinding two ends of each cylinder. The ultrasonic testing device disclosed by the invention is used for carrying out ultrasonic experimental measurement on 12 compact sandstone samples by using an ultrasonic pulse testing device, adopts an experimental device disclosed in the literature (Guomeishiu et al, 2009; Yan et al, 2011), and consists of a reaction container, a temperature control unit, a confining pressure control unit, a pore pressure control unit, a sound wave testing unit and the like.
For experimental observation of dried samples, samples were dried by the method of the literature (wei yi jun et al,2020), the rock samples were sealed with rubber sleeves and placed in an experimental apparatus, effective pressures (5MPa, 10MPa, 15MPa, 20MPa, 25MPa, 30MPa and 35MPa) were applied to each sample (here, 12 samples were applied in sequence, and 7 pressures were applied to each sample), the fluid for applying confining pressure was heated by an electric heating wire in the apparatus to 80 ℃ and maintained for half an hour, ultrasonic experimental observation was performed, and the ultrasonic waveform passing through the rock samples was recorded. And (4) observing a water saturation experiment, namely saturating the compact sandstone sample with water by using a vacuumizing and pressurizing saturation method, performing ultrasonic experiment measurement according to a drying experiment process, and recording the ultrasonic waveform passing through the rock sample. And calculating the longitudinal and transverse wave velocities of the compact sandstone under different effective pressures based on the extracted initial values of the longitudinal and transverse wave waveforms. Specific petrophysical parameters of the samples are shown in table 1.
TABLE 1 physical Properties of tight sandstone rock
Sample number | Porosity (%) | Permeability (mD) | Particle Density (kg/m)3) |
GAR11 | 3.03 | 0.001 | 2687 |
GA3 | 3.46 | 0.005 | 2694 |
GAR1 | 4.19 | 0.0005 | 2676 |
GA10 | 4.89 | 0.004 | 2691 |
GA6 | 6.26 | 0.046 | 2672 |
GAR6 | 6.33 | 0.047 | 2665 |
GA8 | 8.55 | 0.082 | 2670 |
GAR7 | 8.65 | 0.028 | 2668 |
GAR12 | 8.97 | 0.14 | 2662 |
GA1 | 13.26 | 1.21 | 2659 |
GAR8 | 13.35 | 1.320 | 2653 |
GA2 | 13.91 | 1.370 | 2660 |
The relationship between the velocity of the longitudinal and transverse waves and the pressure of the dry and water-saturated states of 12 groups of tight sandstone rock samples is reflected in fig. 2. The longitudinal and transverse wave velocities in the saturated water state are greater than those in the dry state, and in the two different states, similar changes are observed between the longitudinal and transverse wave velocities and the ambient pressure, both the longitudinal and transverse wave velocities increase with increasing pressure, the rate of change of velocity with increasing pressure is greater when the pressure is lower, and the velocity slows with increasing pressure (Carcione and Quiroga-Goode, 1995; Dengdeng et al 2015; Song Negteng et al 2015). The main reason for this is that since the fractures close successively during the pressure increase, the effect of pressure on the hard pores is very small and negligible (Chen 39065and Huang Fang, 2001; David and Zimmerman, 2012; Shapiro, 2003). Based on the method, the fracture distribution characteristics can be obtained by utilizing the speed-pressure change relation.
Establishing a quantitative relation between the elastic modulus and hard pores in the rock based on a Mori-Tanaka model, wherein the equivalent elastic modulus can be expressed as:
wherein,andthe equivalent bulk modulus and the equivalent shear modulus of the rock, K, are respectively when only the hard pores are contained0And mu0Respectively the bulk modulus and shear modulus of the rock particles, phisIs the hard porosity of the rock, and M and N are the shape factors of the hard porosity, defined as follows:
wherein the shape of the hard pores is default to elliptical pores, gamma is the aspect ratio of the elliptical pores, upsilon is the Poisson ratio of rock particles, and upsilon is (3K)0-2μ0)/(6K0+2μ0) In the formula, the expression of g is as follows:
if the fractures are added to a pore medium containing hard pores, the content of the fractures is small, no interaction exists between the fractures and the pores, and the equivalent elastic modulus of the rock can be expressed as:
wherein,is the poisson ratio of the rock when only hard pores are present,and withThe fracture density is the equivalent volume modulus and the equivalent shear modulus when the rock contains fractures and hard pores, wherein gamma is the accumulated fracture density, and the fractures are approximately and completely closed under higher effective pressure, and the rock only contains the hard pores, so that the optimal aspect ratio of the hard pores can be obtained by a formula (1) and a least square method.
The pressure dependence of the effective bulk modulus of the rock is related to fracture density, so when the fracture density is known, the elastic modulus of the rock can be obtained through formulas (5) and (6), conversely, the fracture density can be obtained through the elastic modulus under each effective pressure, and the accumulated fracture density under each pressure can be calculated by using a least square method.
quantitative relationship of cumulative fracture density as a function of effective pressure:
in the formula, gammaiIs the initial fracture density at zero pressure, p is the pressure,is a pressure constant of the same order of magnitude as the pressure.
Dividing pressure intervals of the dry rock test into n parts at equal intervals, and solving fracture density of each pressure by using a formula (7):
wherein k is more than or equal to 0 and less than or equal to n.
Dividing the pressure interval of the dry rock test at equal intervals delta p into n parts, wherein each part is delta p, and the fracture aspect ratio at each pressure point is as follows:
wherein,is the effective poisson's ratio at high pressure,is the effective young's modulus at high pressure.
Fracture porosity is related to fracture density by (David and Zimmerman, 2012):
and (3) solving the fracture aspect ratio corresponding to the maximum fracture density by using the formula (10), wherein the fracture aspect ratio is the main fracture aspect ratio.
The aspect ratio of fractures in actual reservoir rock is not a fixed value, but is continuously distributed within a certain range (dunne, 2015). However, at a certain fixed pressure, the major fracture aspect ratio is a fixed value. By major fracture aspect ratio is meant the fracture aspect ratio corresponding to the maximum fracture porosity (Sun et al,2020), and thus the invention selects the major fracture aspect ratio, and figure 4 shows the fracture aspect ratio of the major fracture at different pressures for the sample GAR 7. In the pressure range of 5-35MPa, when the pressure is different, the fracture aspect ratio is similar to the distribution form of the fracture porosity, but the fracture porosity corresponding to each fracture aspect ratio is reduced when the pressure is higher, and the reduction of the fracture porosity is more obvious when the pressure is lower. Mainly due to the increased pressure, which closes part of the fracture.
And 4, deriving an MFS model:
walsh (1965), Mavko and Jizaba (1991) and Shapiro (2003) believe that the pores of a rock can be divided into hard and soft pores, with hard pores constituting the primary porosity of the rock. Soft pores (fissures) have a very small aspect ratio and a very narrow space compared to pores inside rocks, and easily expand and contract under external force, thereby forming a jet flow.
The MFS model considers the heterogeneity of pore media, divides pores into fractures and hard pores, meanwhile, incorporates fractures (jet flows) into a solid framework, has a structure shown in figure 5, corrects the rock framework (similar to a new rock framework), calculates the bulk modulus and shear modulus of the corrected rock, and analyzes the reservoir elastic wave propagation characteristics.
Using the Dvorkin and Nur (1993) concept, assuming an ideal cylinder of radius R, the relationship between the deformation of the framework (e ═ du/dx) and the fracture and pore pressures is as follows:
where u is the solid displacement and x is parallel to the wave propagation direction. c. C0Is the acoustic velocity, p, of the liquidfIs the fluid density, Qc=K0/(1-Kmsd/K0-φc),Kmsd=(1/K0-1/Khp+1/Kd)-1,KdIs the bulk modulus, phi, of the dry rockcIs the fracture porosity, KhpIs the bulk modulus of the dry rock when the fracture is fully closed.
Assuming that the liquid inside the fracture can only flow perpendicular to the direction of wave propagation, the fluid conservation of mass equation is therefore:
where t is the time and w is the displacement of the liquid.
Substituting the equations (11) and (12) into the equation (13) can obtain:
the relationship between fluid velocity and pore pressure gradient can be derived by Darcy's law:
where η is the fluid viscosity and κ is the permeability.
The simple harmonic expression is:
exp (-i ω t) (16) substituting equation (16) into equations (14) and (15) yields:
wherein, KfIs the bulk modulus of the fluid.
From equations (17) and (18), the differential equation of the pressure is obtained as follows:
at the fracture boundary, the pore pressure is always constant, so at R, p, dp, 0, the result of solving equation (19) is:
wherein R is the characteristic jet length, λ2=iωηφc/κ(1/Kf+1/(φcQc)),J0Is a zero order bessel function.
The average fluid pressure is:
wherein, J1Is a first order bessel function.
(Biot, 1941; Rice and clean, 1976; Dvorkin and Nur,1993) points out the relationship between external stress and displacement and pore pressure in rock, and thus, in modified rock, external stress and displacement and pore mean pressure are as follows:
dσ=Kmsdde-(1-Kmsd/K0)dpav (22)
modified bulk modulus of solids Kms:
Substituting equations (21) and (22) into equation (23) yields:
the dry bulk modulus of the modified rock can be derived from the definition of bulk modulus, the following relationship (Walsh,1965) between dry pore volume, matrix modulus and dry rock modulus.
Mavko and Jizba (1991) gives the expression for modifying the dry shear modulus of rock as follows:
wherein, mudIs the shear modulus of dry rock according toand Johnston,1981) wave number versus wave velocity and attenuation, and deriving expressions for longitudinal wave velocity and attenuation factor.
ρaFor fluid-solid coupling density, ρsIs the density of the matrix and is,is the characteristic frequency. MdIs the uniaxial stress modulus of the rock skeleton,
the numerical simulation calculation parameters of the present invention are according to Carcione and Gurenvich (2011). The values of all parameters are respectively porosity of 0.2, solid framework bulk modulus of 18GPa, Poisson ratio of 0.15 and matrix density of 2650kg/m3Matrix bulk modulus 50GPa, permeability 2X 10-15m2Fluid density 1040kg/m3The viscosity of the fluid is 0.1cp, the bulk modulus of the fluid is 2.25GPa, and the fluid-solid coupling density is 420kg/m3The characteristic jet length is 2mm, the fracture porosity is 0.0002, the bulk modulus of the dry rock is 20GPa when no fracture exists, and the aspect ratio of the fracture is 0.0008. And calculating the longitudinal wave velocity and drawing a curve graph.
Analysis of FIG. 6 reveals that the two models (MFS model and Dvorkin (1995)) in the graph are in agreement with the Biot model calculations of longitudinal wave velocity in the low frequency band. However, the calculated compressional velocity of the Dvorkin (1995) model exceeds the maximum value (Vpmax), thus limiting the application of the model. The main reason for this is that the boundary conditions are not set properly. The maximum value is obtained by setting the dry bulk modulus of the corrected rock to the bulk modulus (K) at which the fracture is completely closedmd=Khp) Then, the velocity of the longitudinal wave is calculated by substituting the value into the Biot equation, and the result is the maximum value of the velocity of the longitudinal wave.
The MFS model calculations proposed by the present invention show that the longitudinal wave velocity approaches the maximum value and does not exceed the maximum value as the frequency increases, which is consistent with the theoretical research results in the industry (Dvorkin and Nur, 1993; Joge, 1997).
In conclusion, in the MFS model provided by the invention, the calculated value of the longitudinal wave velocity is consistent with the calculated value of the Gassmann-Biot model in the low frequency band, and the longitudinal wave velocity calculated in the high frequency band does not exceed the theoretical maximum value. Therefore, the model provided by the invention can reflect the reservoir elastic wave propagation characteristics more accurately and comprehensively.
And 5, calculating the characteristic jet flow length by using an MFS model:
based on the MFS model derived in step 4, the characteristic jet flow lengths for different pressures are calculated by the model, and the most suitable characteristic jet flow length for each sample is optimized by using the least square method.
FIG. 7 is a graph showing the dispersion characteristics of GAR7 samples at different pressures, wherein the pressure is in the range of 5-35MPa and the characteristic jet length values are: 0.055mm, 0.051mm, 0.049mm, 0.044mm, 0.041mm, 0.038mm and 0.030 mm. The relevant parameters of the fluid were calculated by the Batzle-Wang equation (Batzle and Wang, 1992). Based on the drying longitudinal and transverse wave speeds and densities measured by experiments, the shear modulus and the bulk modulus of the dry rock can be calculated.
As can be seen from fig. 7, a characteristic jet length can be found at each pressure, so that the predicted value of the model is consistent with the experimentally measured value, and therefore, the model can effectively describe the characteristics of the velocity and the pressure change of the tight sandstone compressional wave. The primary reason for the decrease in the dispersion of compressional velocity in tight sandstone with increasing pressure is the change in the internal micro-pore structure of the rock caused by the increase in pressure.
FIG. 8 shows the comparison of compressional wave velocity with experimental values for the jet length prediction samples GAR7 with different characteristics. The characteristic jet flow length is increased, the longitudinal wave speed predicted by the model is increased, the relation between the pressure and the model predicted value and the experimentally measured value is consistent, namely the pressure is increased, the speed is increased rapidly when the pressure is lower, and the speed is increased slowly when the pressure is higher. With increasing pressure, the appropriate characteristic jet length decreases. The main reason for this is that when the pressure increases, the microcracks, which have small aspect ratios, close. As can be seen from the comparison of the model calculated values and the experimental measurements for tight sandstone samples in fig. 8, the rock samples can be characterized by a fixed characteristic jet length under different pressure conditions. It follows that the characteristic jet length is a fundamental property of rock (Dvorkin and Nur, 1993).
and (5) analyzing the influence of the fracture parameters on the length of the characteristic jet flow through the computed fracture distribution parameters and the characteristic jet flow length in the step 2-5.
Figure 9 is a graph of the characteristic jet length of 12 pieces of tight sandstone plotted against various parameters. The characteristic jet length increases with increasing permeability, increasing fracture density, and smaller major fracture aspect ratios.
fracture density, which is a parameter used to describe the number of fractures, and fracture aspect ratio, which is a parameter indicative of the fracture shape, are the main parameters that characterize the microstructure of the fracture, both parameters affecting the jet effect (dawn, 2010). The fracture density and fracture aspect ratio can be obtained through model calculation according to experimental data, and the characteristic jet flow length is a variable introduced for researching the local movement of fluid in the fracture. Analysis has shown that for rock samples, there is a characteristic jet flow length that is most appropriate over a range of pressures, while at a fixed pressure, the major fracture aspect ratio and fracture density are fixed values. Permeability is an indication of how easily fluid in the pore medium flows in its pores and can be measured directly by an instrument to obtain this parameter, and analysis has shown that permeability also affects the characteristic jet length. Thus, the characteristic jet length is related to parameters such as fracture density, major fracture aspect ratio, and permeability.
In order to comprehensively clarify the relation between the characteristic jet flow length and each parameter, the invention introduces the characteristic frequency of the jet flow and the critical fluid diffusion length to more clearly characterize the characteristic of the characteristic jet flow length.
The authors indicate that jet eigenfrequency is closely related to bulk modulus, fracture aspect ratio and fluid viscosity (O' Connell and Budiansky, 1978; Jones et al, 1980; Li et al,2018), and that the expression:
the critical diffusion length of the fluid is given according to the characteristic frequency of the jet flow as follows:
from equations (31) and (32), the critical diffusion length of the fluid is related to the fracture aspect ratio and permeability parameters, and is expressed as:
when the fluid diffusion length is less than the critical diffusion length, the pore fluid pressure has sufficient time to equilibrate with the fluid flow, whereas the pore fluid pressure has insufficient time to equilibrate.
Based on previous studies (Pride, 2004; Dvorkin and Nur,1993) and the present invention analysis, it was suggested that parameters such as major fracture aspect ratio, permeability, matrix modulus, bulk modulus of fluid and fracture density, etc. characterize the characteristic jet length, the formula proposed in the present invention is:
and fitting values of the parameters a and b according to the measurement result of the ultrasonic experiment, wherein the values have certain difference for different rock samples. The invention fits the values of a and b to 0.03865 and 63.42, respectively, based on experimental data for 12 tight sandstones. Accordingly, equation (34) can be expressed as:
figure 10 visually depicts the quantitative relationship between the characteristic jet length and various parameters.
In conclusion, the invention fully considers the relationship between the characteristic jet flow length of the tight sandstone and the fracture distribution characteristics in the rock, and forms a method for determining the characteristic jet flow length. The method characterizes the characteristic jet flow length through measurable parameters, so that the characteristic jet flow length is converted to be acquired through an indirect measurement mode.
The above description is only a preferred embodiment of the present invention, and the scope of the present invention is not limited thereto, and any simple modifications or equivalent substitutions of the technical solutions that can be obviously obtained by those skilled in the art within the technical scope of the present invention are within the scope of the present invention.
Claims (5)
1. A method of estimating the characteristic jet length of tight sand comprising the steps of:
step 1, acquiring ultrasonic velocity experimental data of pressure change of tight sandstone;
step 2, estimating the fracture density of the compact sandstone;
step 3, estimating the aspect ratio of main fractures of the compact sandstone;
step 4, constructing an MFS model;
and 5, calculating the characteristic jet flow length by using an MFS model, specifically: calculating the characteristic jet flow length of different pressures by using the MFS model constructed in the step 4, and preferably selecting the characteristic jet flow length of each sample by using a least square method;
step 6, analyzing the relation between the characteristic jet flow length and the fracture distribution parameters;
step 7, selecting the main fracture aspect ratio gamma, the permeability kappa and the matrix bulk modulus K0Fluid bulk modulus KfAnd characterizing the characteristic jet flow length R with the accumulated fracture density gamma to form a characteristic jet flow length quantitative relation, wherein the characteristic jet flow length quantitative relation of the compact sandstone is as follows:
fitting values of the parameters a and b according to the measurement result of the ultrasonic experiment, wherein the values have respective differences for different rock samples;
wherein, the step 4 of constructing the MFS model comprises the following steps:
setting an ideal cylinder with radius R, the relationship between the deformation of the framework and the fracture and pore pressure is as follows:
wherein p is pressure, u is solid displacement, c0Is the acoustic velocity, p, of the liquidfIs the fluid density, Qc=K0/(1-Kmsd/K0-φc),Kmsd=(1/K0-1/Khp+1/Kd)-1,KdIs the bulk modulus, phi, of the dry rockcIs the fracture porosity, KhpIs the bulk modulus of the dry rock when the fracture is fully closed;
assuming that the liquid inside the fracture can only flow in the direction perpendicular to the wave propagation direction, the fluid conservation of mass equation is:
wherein t is time and w is displacement of the liquid;
substituting the equations (11) and (12) into the equation (13) can obtain:
the relationship between fluid velocity and pore pressure gradient can be derived:
wherein η is fluid viscosity and κ is permeability;
the simple harmonic expression is:
exp(-iωt) (16)
substituting equation (16) into equations (14) and (15) yields:
wherein, KfIs the bulk modulus of the fluid;
from equations (17) and (18), the differential equation of the pressure is derived as follows:
at the fracture boundary, the pore pressure is always constant, so at R, p, dp, 0, the result of solving equation (19) is:
wherein R is the characteristic jet length, λ2=iωηφc/κ(1/Kf+1/(φcQc)),J0Is a zero order Bessel function;
the average fluid pressure is:
wherein, J1Is a first order Bessel function;
in the modified rock, the external stress and displacement and pore mean pressure are as follows:
dσ=Kmsdde-(1-Kmsd/K0)dpav (22)
modified bulk modulus of solids Kms:
Substituting equations (21) and (22) into equation (23) yields:
the following relationships exist between the bulk modulus definition, the dry pore volume, the matrix bulk modulus and the dry rock modulus, and the dry bulk modulus of the modified rock can be derived
The expression for modifying the dry shear modulus of the rock is as follows:
wherein, mudThe shear modulus of the dry rock is obtained by obtaining expressions of longitudinal wave velocity and attenuation factor according to the relationship between wave number and wave velocity and attenuation
2. the method for estimating the characteristic jet flow length of tight sandstone according to claim 1, wherein the step 1 is to obtain ultrasonic velocity experimental data of the pressure change of tight sandstone, and comprises the following steps: and testing physical property parameters of the rock of the compact sandstone sample under the pressure change in a dry state and a water-saturated state by using an ultrasonic experimental measurement method, and calculating longitudinal wave speed and transverse wave speed in the dry and water-saturated states.
3. The method for estimating the tight sandstone characteristic jet flow length as claimed in claim 1, wherein the step 2 of estimating the tight sandstone fracture density comprises the following steps:
establishing a quantitative relationship between the elastic modulus and the hard pores in the rock, the equivalent elastic modulus can be expressed as:
wherein,andthe equivalent bulk modulus and the equivalent shear modulus of the rock, K, are respectively when only the hard pores are contained0And mu0Respectively the bulk modulus and shear modulus of the rock particles, phisIs the hard porosity of the rock, and M and N are the shape factors of the hard porosity, defined as follows:
wherein the shape of the hard pores is default to elliptical pores, gamma is the aspect ratio of the elliptical pores, upsilon is the Poisson ratio of rock particles, and upsilon is (3K)0-2μ0)/(6K0+2μ0) In the formula, the expression of g is as follows:
if the fractures are added to a pore medium containing hard pores, the content of the fractures is small, no interaction exists between the fractures and the hard pores, and the equivalent elastic modulus of the rock can be expressed as:
wherein,andrespectively the equivalent bulk modulus and the equivalent shear modulus when the rock contains cracks and hard pores,is the rock poisson's ratio with only hard pores, Γ is the cumulative fracture density; because the fracture is approximately and completely closed under the effective pressure, the rock only contains hard pores, and the optimal aspect ratio of the hard pores is obtained by adopting a least square method through a formula (1);
the pressure correlation of the effective bulk modulus of the rock is related to fracture density, when the fracture density is known, the equivalent bulk modulus and the equivalent shear modulus of the rock can be obtained through formulas (5) and (6), and conversely, the fracture density under each pressure can be calculated through the equivalent bulk modulus and the equivalent shear modulus under each effective pressure by using a least square method.
4. The method for estimating the characteristic jet length of tight sandstone according to claim 1, wherein the step 3 of estimating the aspect ratio of the main fracture of tight sandstone comprises the following steps:
quantitative relationship of cumulative fracture density as a function of effective pressure:
in the formula, gammaiIs the initial fracture density at zero pressure, p is the pressure,is a pressure constant of the same order of magnitude as the pressure;
dividing the pressure interval delta p of the dry rock test into n parts, and solving the fracture density of each pressure by using a formula (7):
wherein k is more than or equal to 0 and less than or equal to n;
dividing the pressure interval of the dry rock test at equal intervals delta p into n parts, wherein each part is delta p, and the fracture aspect ratio at each pressure point is as follows:
wherein,is the effective poisson's ratio at high pressure,is the effective young's modulus at high pressure;
the relationship between fracture porosity and fracture density is:
and (3) solving the fracture aspect ratio corresponding to the maximum fracture density by using a formula (10), namely the main fracture aspect ratio.
5. The method for estimating the length of the tight sandstone characteristic jet flow according to claim 1, wherein the quantitative length relation of the tight sandstone characteristic jet flow length in the step 7 is derived by the following steps: the jet characteristic frequency is closely related to the matrix bulk modulus, the main fracture aspect ratio and the fluid viscosity, and the expression is as follows:
the critical diffusion length of the fluid is given according to the characteristic frequency of the jet flow as follows:
as can be seen from equations (31) and (32), the critical diffusion length of the fluid is expressed as:
when the fluid diffusion length is less than the critical diffusion length, the pore fluid pressure has enough time to reach equilibrium through the fluid flow, otherwise, the pore fluid pressure has insufficient time to reach equilibrium;
selecting the main fracture aspect ratio, permeability, matrix bulk modulus, fluid bulk modulus and accumulated fracture density to characterize the characteristic jet flow length to form a characteristic jet flow length quantitative relation, wherein the characteristic jet flow length quantitative relation of the compact sandstone is as follows:
and fitting the values of the parameters a and b according to the measurement result of the ultrasonic experiment, wherein the values have respective differences for different rock samples.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111101795.1A CN113887150B (en) | 2021-09-18 | 2021-09-18 | Method for estimating length of characteristic jet flow of compact sandstone |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111101795.1A CN113887150B (en) | 2021-09-18 | 2021-09-18 | Method for estimating length of characteristic jet flow of compact sandstone |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113887150A CN113887150A (en) | 2022-01-04 |
CN113887150B true CN113887150B (en) | 2022-05-06 |
Family
ID=79009950
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111101795.1A Active CN113887150B (en) | 2021-09-18 | 2021-09-18 | Method for estimating length of characteristic jet flow of compact sandstone |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113887150B (en) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116840912B (en) * | 2022-09-13 | 2024-07-19 | 上海电子信息职业技术学院 | Method for predicting longitudinal wave velocity in partially saturated pore medium with different fracture aspect ratios |
CN116227139B (en) * | 2022-12-23 | 2023-08-29 | 河海大学 | Method for measuring and calculating transmission properties of hole-crack rock based on two-part Hooke's law |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109116420A (en) * | 2018-10-16 | 2019-01-01 | 河海大学 | A kind of pore media velocity of longitudinal wave containing crack and decaying prediction technique |
CN109374497A (en) * | 2018-10-29 | 2019-02-22 | 河海大学 | A kind of rock micropore structure test method |
CN109471168A (en) * | 2018-11-06 | 2019-03-15 | 河海大学 | The prediction technique of velocity of longitudinal wave and decaying in a kind of hole fissuted medium |
CN113075728A (en) * | 2021-02-26 | 2021-07-06 | 河海大学 | Method for establishing compact sandstone multi-scale three-dimensional rock physical drawing board |
-
2021
- 2021-09-18 CN CN202111101795.1A patent/CN113887150B/en active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109116420A (en) * | 2018-10-16 | 2019-01-01 | 河海大学 | A kind of pore media velocity of longitudinal wave containing crack and decaying prediction technique |
CN109374497A (en) * | 2018-10-29 | 2019-02-22 | 河海大学 | A kind of rock micropore structure test method |
CN109471168A (en) * | 2018-11-06 | 2019-03-15 | 河海大学 | The prediction technique of velocity of longitudinal wave and decaying in a kind of hole fissuted medium |
CN113075728A (en) * | 2021-02-26 | 2021-07-06 | 河海大学 | Method for establishing compact sandstone multi-scale three-dimensional rock physical drawing board |
Non-Patent Citations (2)
Title |
---|
致密岩石多频段岩石物理实验研究;李越;《中国优秀硕士学位论文全文数据库 基础科学辑》;20200215;第3-4章 * |
致密碳酸盐岩跨频段岩石物理实验及频散分析;李闯等;《地球物理学报》;20200228;第628-634页 * |
Also Published As
Publication number | Publication date |
---|---|
CN113887150A (en) | 2022-01-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN113887150B (en) | Method for estimating length of characteristic jet flow of compact sandstone | |
CN109471168B (en) | Method for predicting velocity and attenuation of longitudinal wave in pore fracture medium | |
CN109116420B (en) | Method for predicting longitudinal wave velocity and attenuation of pore medium containing fractures | |
CN114236609B (en) | Prediction method for longitudinal wave velocity and attenuation of partially saturated hole fracture medium | |
CN109374497B (en) | Rock micro-pore structure testing method | |
Mockovčiaková et al. | Study of the relation between the static and dynamic moduli of rocks | |
Sergeev et al. | Experimental research of the colloidal systems with nanoparticles influence on filtration characteristics of hydraulic fractures | |
CN113075728B (en) | Method for establishing compact sandstone multi-scale three-dimensional rock physical drawing board | |
CN104020276A (en) | Determination method for mechanical parameters of transverse isotropy shale reservoir rocks | |
CN112505772A (en) | Method for inverting rock pore distribution characteristics by utilizing pore and fracture medium elastic wave theory | |
Sun et al. | Thermal modulation of nonlinear ultrasonic wave for concrete damage evaluation | |
Knabe et al. | Permeability characterization on tight gas samples using pore pressure oscillation method | |
Khlybov | Studying the effect of microscopic medium inhomogeneity on the propagation of surface waves | |
Rajaoalison et al. | Mechanical properties of sandstone using non-destructive method | |
Ma et al. | Velocity dispersion and fluid substitution in sandstone under partially saturated conditions | |
Wang et al. | Fracture mechanical properties and identification of principal stress direction from the rock frequency domain | |
Murphy III | Seismic to ultrasonic velocity drift: intrinsic absorption and dispersion in crystalline rock | |
Poplygin et al. | Influence of frequency of wave action on oil production | |
Spalvier et al. | The stress-induced surface wave velocity variations in concrete | |
Wang et al. | Study on the Effect of Bedding Plane and Loading Style on Fracture Process Zone in Coal | |
Chowdhury et al. | Seismic behaviour of CO2 saturated Fontainebleau sandstone under in situ conditions | |
CN116840912B (en) | Method for predicting longitudinal wave velocity in partially saturated pore medium with different fracture aspect ratios | |
Nikolenko et al. | The Effect of Limestone Porosity on the Velocity of P-and S-Waves under Mechanical and Thermal Loading | |
Li et al. | Fracture extension behavior and micro-mechanical damage mechanism under different CO2 fracturing methods in continental shale reservoirs | |
He et al. | Laboratory study of acoustic velocity in different types of rocks at seismic frequency band |
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 |