CN109116420B - Method for predicting longitudinal wave velocity and attenuation of pore medium containing fractures - Google Patents
Method for predicting longitudinal wave velocity and attenuation of pore medium containing fractures Download PDFInfo
- Publication number
- CN109116420B CN109116420B CN201811200996.5A CN201811200996A CN109116420B CN 109116420 B CN109116420 B CN 109116420B CN 201811200996 A CN201811200996 A CN 201811200996A CN 109116420 B CN109116420 B CN 109116420B
- Authority
- CN
- China
- Prior art keywords
- fluid
- attenuation
- longitudinal wave
- bisq
- fracture
- 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
- 239000011148 porous material Substances 0.000 title claims abstract description 50
- 238000000034 method Methods 0.000 title claims abstract description 28
- 239000012530 fluid Substances 0.000 claims abstract description 88
- 239000011435 rock Substances 0.000 claims abstract description 52
- 230000010355 oscillation Effects 0.000 claims abstract description 7
- 238000006073 displacement reaction Methods 0.000 claims description 26
- 230000014509 gene expression Effects 0.000 claims description 21
- 239000007787 solid Substances 0.000 claims description 17
- 239000000126 substance Substances 0.000 claims description 12
- 230000035699 permeability Effects 0.000 claims description 7
- 239000011159 matrix material Substances 0.000 claims description 5
- 230000008859 change Effects 0.000 claims description 4
- 238000004458 analytical method Methods 0.000 claims description 3
- 238000009795 derivation Methods 0.000 claims description 3
- 239000007788 liquid Substances 0.000 claims description 3
- 238000005507 spraying Methods 0.000 claims description 2
- 230000000694 effects Effects 0.000 abstract description 3
- 230000008569 process Effects 0.000 abstract description 2
- 238000004088 simulation Methods 0.000 abstract description 2
- 239000006185 dispersion Substances 0.000 description 6
- 238000005259 measurement Methods 0.000 description 4
- 230000002051 biphasic effect Effects 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 3
- 230000007246 mechanism Effects 0.000 description 3
- 238000002474 experimental method Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 229920006395 saturated elastomer Polymers 0.000 description 2
- 240000007594 Oryza sativa Species 0.000 description 1
- 235000007164 Oryza sativa Nutrition 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000003776 cleavage reaction Methods 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000005284 excitation Effects 0.000 description 1
- 230000004807 localization Effects 0.000 description 1
- 238000000386 microscopy Methods 0.000 description 1
- 235000009566 rice Nutrition 0.000 description 1
- 230000007017 scission Effects 0.000 description 1
- 239000002344 surface layer Substances 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/30—Analysis
- G01V1/306—Analysis for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/62—Physical property of subsurface
- G01V2210/624—Reservoir parameters
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/62—Physical property of subsurface
- G01V2210/624—Reservoir parameters
- G01V2210/6244—Porosity
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
The invention provides a method for predicting the longitudinal wave velocity and attenuation of a pore medium containing fractures, and belongs to the field of seismic rock physics. According to the invention, according to the BISQ model principle, the oscillation of the fluid in the fracture in the direction vertical to the wave propagation direction is considered. Starting from the description of two important characteristic parameters of fracture aspect ratio and fracture density, the BISQ model is improved based on the law of conservation of fluid mass, an improved theoretical model (Re-BISQ model for short) is deduced, the characteristics of the improved theoretical model are analyzed by numerical simulation, and experimental data of sandstone in Nanampton, England and compact sandstone in Su Li gas field in Ordos basin are used for verifying the effect of the invention. The result shows that the influence of the internal fracture of the underground rock in the wave propagation process is fully considered in the prediction of the longitudinal wave velocity and the attenuation, and an important basis is provided for reservoir prediction.
Description
Technical Field
The invention relates to the field of seismic rock physics, in particular to a method for predicting the longitudinal wave velocity and attenuation of a pore medium containing fractures.
Background
Reservoir media elastic wave propagation theory has undergone a move from single pore media to more complex pore media theory. Research on the wave propagation theory of single-pore media began at the first 50 s of the 20 th century, and Gassmann (1951) discussed the effect of the presence of pore fluid in single-pore media on the elastic properties of the media, without considering the relative motion between the pore fluid and the scaffold. Biot (1956 a; 1956 b; 1962) proposed the elastic wave kinetic equation of a fluid-solid coupled biphasic system, which lays the foundation of theoretical studies on biphasic media and predicts the existence of slow longitudinal waves in biphasic media. However, the Biot theoretical model also has the disadvantage (White, 1975; Mavko and Nur,1975) that it is difficult to explain the phenomena of high dispersion and strong attenuation of the waves (Dvorkin and Nur, 1993). The main reason for this is that the theoretical model assumes that all pores inside the rock are uniform (i.e. single pore medium model), and that the fluid under the excitation of the wave only undergoes pipe-in-pipe laminar flow in the wave propagation direction without other reactions, and these basic assumptions do not correspond to the complex situation inside the actual heterogeneous rock.
To address the above problems, some scholars consider the effect of the micro-flow of the pore fluid on wave propagation. Mavkoand Nur (1975) proposed a jet mechanism based on pore micro-geometry, successfully explaining the phenomena of high frequency dispersion and strong attenuation of waves, but this theory relies heavily on pore micro-geometry and is thus difficult to popularize and apply, thus also causing the cleavage of Biot flow mechanism and jet mechanism, which is not in accordance with physical reality. Dvorkin and Nur (1993) have established a BISQ (Biot/Squick) model based on the assumption that the pores contain saturated fluid, taking into account the Biot flow and the squrit flow in a mechanical model, so that the phenomena of high dispersion and strong attenuation of the waves are well explained. The BISQ model is improved by avoiding the characteristic jet flow length, but the predicted longitudinal wave velocity dispersion and attenuation of the BISQ model are shifted to a high frequency band along with the shift of the rock permeability, which is contrary to the result predicted by the BISQ model and is not consistent with the general knowledge of local fluid flow of scholars.
Pride (2004) indicates that the BISQ model defines pores in rock as being spatially continuous uniformly distributed pores, an assumption that is not consistent with physical reality. Dawn (2011) states that the BISQ theory does not relate to two important parameters in fractured media: fracture density and aspect ratio, and thus dawn derives the wave equation from two important parameters of the fracture. However, the model only increases the fracture parameters and cannot change the dispersion and attenuation of the fast longitudinal wave, so that the model parameters become more complex.
Disclosure of Invention
The invention aims to overcome the defects in the prior art and provides a method for predicting the longitudinal wave velocity and attenuation of a pore medium containing fractures.
The method is based on the BISQ model, combines the density and the aspect ratio characteristics of the fractures in the rock, mainly considers the influence of the coin-shaped fractures on wave propagation, and improves the BISQ model from the perspective of fluid mass conservation, so that the BISQ model is not complicated, and the wave propagation characteristics of the reservoir can be well reflected.
The invention provides a method for predicting the longitudinal wave velocity and attenuation of a pore medium containing fractures, which is characterized by comprising the following steps of:
step S1, characterizing the internal pore structure of the rock;
step S2, characterizing the wave propagation characteristics in the rock with cracks;
step S3, estimating the proportion of the porosity of the rock fractures to the total porosity;
step S4, deducing a Re-BISQ model;
step S5, predicting the longitudinal wave velocity and the inverse quality factor Q by using the Re-BISQ model-1;
Step S6, predicting the longitudinal wave velocity and the inverse quality factor Q by using the BISQ model-1;
Step S7, result analysis.
In step S1, the internal pore structure of the rock is characterized by using the morphology of the rock slice under a microscope.
In step S2, the propagation characteristics of the seismic waves inside the rock are represented by the propagation characteristic image of the elastic waves in the cylinder sponge.
In step S3, an expression of the ratio of fracture porosity to total porosity:
phi is the total porosity, phicIs the fracture porosity, D is the ratio of the fracture porosity to the total porosity; the fracture porosity expression isZeta is the aspect ratio of the fissures and ε is the fissure density.
Through steps S4-S6, longitudinal wave velocity and inverse quality factor of the two models BISQ and Re-BISQ predictions are calculated.
The differences between the two models BISQ and Re-BISQ in predicting the longitudinal wave velocity and the inverse quality factor are comparatively analyzed by step S7.
In step S4, in the case of one dimension, the equation for the oscillation of the fluid in the two-phase medium can be obtained as follows:
where t denotes time, κ and η denote porosity, permeability and fluid viscosity, respectively, u and w denote solid displacement and fluid displacement in the x-direction, respectively, and ρs、ρfAnd ρaRepresenting solids, fluids, and add-on densities; m is the Biot elasticity parameter; p is the fluid pressure;is the pore elasticity parameter, KdryIs a solid skeleton bulk modulus, K0Is the matrix bulk modulus;
when the rock is squeezed, the fluid inside the rock fracture is sprayed to the hard pores, so that a spraying flow is generated, and the following equation can be obtained by combining the conservation of the mass of the fluid:
wherein the content of the first and second substances,v is the ratio of the fracture to the total porosity, v is the displacement of the fluid in the direction r; according to David and Zimmerman (2012), the expression of the fracture porosity isZeta is the aspect ratio of the crack, epsilon is the crack density;
the relationship between porosity and framework and fluid is as follows:
substituting (4) and (5) into (3),
wherein the content of the first and second substances,γ=α-φ,and c is0Is the speed of sound in the liquid;
through Lagrange's equation and combining Biot (1956) dissipation function, the fluid dissipation function of the oscillation perpendicular to the propagation direction of the longitudinal wave is obtained, namely the fluid dissipation fracture dissipation function relationship is as follows:
the r-direction fluid displacement expression and the change form of the fracture fluid pressure are respectively set as follows:
v(x,r,t)=v0(r)ei(lx-ωt), (8a)
P(x,r,t)=P0(r)ei(lx-ωt). (8b)
substituting the formula (8) into the formula (7), the fluid pressure gradient equation in the vertical longitudinal wave propagation direction is derived as follows:
the solid displacement and the displacement of the fluid in the x direction are set as follows:
u(x,t)=C1ei(lx-ωt), (10a)
w(x,t)=C2ei(lx-ωt), (10b)
wherein, C1And C2Is a constant;
substituting equations (7) - (9) into equation (5) yields the following relationship between fluid and pressure:
solving the equation using the constant voltage boundary results in:
wherein, J0Is a function of the zero order bessel function,
the average fluid pressure can be derived from equations (8a), (8b), (10a) and (10 b):
wherein, J1Is a first order Bessel function;
the partial derivative is calculated from the above equation (11), so that the following relation can be obtained:
assuming the relationship:
in step S5, the longitudinal wave velocity V is obtainedpAnd attenuation factor α expression;
a1,2=ωIm(X1,2) (16)
where ω is the angular velocity, ωc=ηφ/κρf,ρ1=(1-φ)ρs,ρ2=φρf,
Expression (16) is expressed in inverse quality factor form as:
compared with the prior art, the invention has the following beneficial effects:
a) considering the influence of the cracks in the rock on the wave propagation characteristics;
b) the BISQ theory is combined, and a wave propagation theoretical equation containing the fracture is deduced;
c) through analyzing the experimental data of the rock samples in the two regions, the method provided by the invention is further proved to be capable of reflecting the propagation characteristics of the internal waves of the rock.
Drawings
FIG. 1 is a schematic flow chart of a method for predicting longitudinal wave velocity and attenuation of a porous medium containing fractures according to the present invention;
FIG. 2 is a graph of the wave propagation characteristics characterizing the presence of fractures within the rock, where black represents the jet, white represents the solid skeleton, the cylinder represents the cylindrical sponge assumed by Dvorkin and Nur (1993), the transverse arrows in the cylinder represent the jet direction and are perpendicular to the wave propagation direction, the vertical single arrow represents the Biot flow direction and coincides with the wave propagation direction, and R represents the characteristic jet length;
FIG. 3 is graphs showing the changes of the longitudinal wave velocity (graphs a and c) and the attenuation (graphs b and d) predicted by the three models (Biot, BISQ and Re-BISQ) with frequency, wherein the broken line represents the longitudinal wave velocity or the attenuation value predicted by the Biot model, the double-dashed line represents the longitudinal wave velocity or the attenuation value predicted by the BISQ model, and the chain-dashed line represents the longitudinal wave velocity or the attenuation value predicted by the Re-BISQ model;
FIG. 4 is a graph showing the relationship between permeability and attenuation of conventional sandstone compressional waves in south Amphon, England predicted by two models (BISQ and Re-BISQ), wherein a dotted line represents compressional wave attenuation values experimentally measured by Kliomentos and McCann (1990), a solid line represents compressional wave attenuation values predicted by the Re-BISQ model, and a dotted line represents compressional wave attenuation values predicted by the BISQ model;
fig. 5 is a relation between the velocity and the porosity of the tight sandstone compressional predicted by two models, wherein a dotted line represents the compressional wave attenuation value measured by the experiment (wanda xing, 2016), a solid line represents the compressional wave attenuation value predicted by the Re-BISQ model, and a dotted line represents the compressional wave attenuation value predicted by the BISQ model.
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 of the invention deduces a wave propagation theoretical equation aiming at the influence of the fracture on the wave propagation, compares the difference of the two models in predicting the velocity and attenuation of the longitudinal wave, and further elaborates the effectiveness of the invention by combining the conventional sandstone and the compact sandstone.
The invention provides a method for predicting the longitudinal wave velocity and attenuation of a pore medium containing fractures, which specifically comprises the following steps as shown in figure 1:
step S1, characterizing the internal pore structure of the rock:
the rock medium observed from the earth's surface layer contains a large number of pores and fissures inside, and the observation of rock slices by microscopy shows that the aspect ratio of the fissures is very small compared to the pores inside the rock, increasing the complexity of the elastic properties of the rock.
Step S2, characterizing the propagation characteristics of the internal wave of the rock containing the cracks:
before introducing the Re-BISQ model, we first recognized a schematic of the model, as shown in figure 2. Dvorkin and Nur (1993) first studied Biot flow and local fluid flow in the same mechanical model, creating an instructive BISQ (Biot/Squrit) model. As shown in fig. 2(a), the actual model considered by Dvorkin and Nur is a cylindrical sponge, and the fluid pressure outside the sponge is always a constant, which we set to zero directly for the convenience of the study. Under the action of plane waves, a characteristic cylinder is selected in the rock, the axial direction of which is parallel to the propagation direction of longitudinal waves, the radius R of which is a characteristic parameter which is an order of magnitude dependent on the rock pore size and is called the characteristic jet length, Dvorkin and Nur (1993) assume that this parameter is a fundamental parameter of the rock, independent of the frequency and the properties of the fluid. Macroscopic Biot flow occurs along the axial direction. When seismic waves press against the rock, the fluid in the pores also oscillates perpendicular to the direction of wave propagation, flowing from the fracture to the hard pores (Mavko and Nur,1975), generating jets. From the established conservation of fluid mass relationship, it can be seen that Dvorkin and Nur (1993) do not distinguish between fractures and pores in rock when calculating the jet flow. In fact, the part of the oscillating fluid is mainly the fluid in the fracture, and the specific flowing condition is shown in fig. 2 (b).
Step S3, estimating the ratio of the porosity of the rock fractures to the total porosity:
the BISQ model does not separate the pores and fractures inside the rock, and considers that the jet flow is generated in the total pores, but actually the jet flow only acts in the fractures, and in order to more accurately characterize the wave propagation characteristic, the ratio of the rock fracture porosity to the total porosity needs to be calculated, and the expression is as follows:
according to David and Zimmerman (2012), the expression of the fracture porosity isZeta is the aspect ratio of the fissures and ε is the fissure density.
Step S4, derivation of Re-BISQ model:
biot (1956a), in the one-dimensional case, the equation for the oscillation of a fluid-containing medium in a two-phase medium can be obtained as follows:
where t denotes time, κ and η denote porosity, permeability and fluid viscosity, respectively, u and w denote solid displacement and fluid displacement in the x-direction, respectively, and ρs、ρfAnd ρaRepresenting solids, fluids, and add-on densities; m is the Biot elasticity parameter; p is the fluid pressure;is the pore elasticity parameter. Wherein KdryIs a solid skeleton bulk modulus, K0Is the matrix bulk modulus.
When the rock is squeezed, the fluid inside the rock fracture is ejected towards the hard pores, creating a jet (Mavkoand Nur,1975), which, combined with the conservation of fluid mass, yields the following equation:
wherein the content of the first and second substances,v is the fraction of the total porosity of the fracture and v is the displacement of the fluid in the direction r. According to David and Zimmerman (2012), the expression of the fracture porosity isZeta is the aspect ratio of the fissures and ε is the fissure density.
Biot (1941), Rice and clean (1976) indicate that the relationship between porosity and scaffold and fluid is as follows:
substituting (4) and (5) into (3),
wherein the content of the first and second substances,γ=α-φ,and c is0Is the speed of sound in the liquid. Q and F are intermediate variables only and have no practical meaning.
Equation of dynamics
Through Lagrange's equation and combining Biot (1956) dissipation function, the fluid dissipation function of the oscillation perpendicular to the propagation direction of the longitudinal wave is obtained, namely the fluid dissipation fracture dissipation function relationship is as follows:
the r-direction fluid displacement expression and the change form of the fracture fluid pressure are respectively set as follows:
v(x,r,t)=v0(r)ei(lx-ωt), (8a)
P(x,r,t)=P0(r)ei(lx-ωt). (8b)
wherein the content of the first and second substances,(both of equations 8a and 8b) is the wave number, i represents the imaginary unit, vttRepresenting one derivative of the displacement of the fluid in the r direction, vtRepresenting the quadratic derivative of the fluid displacement in the r direction. V is0Representing the mode of fluid displacement in the r direction.
Substituting the formula (8) into the formula (7), the fluid pressure gradient equation in the vertical longitudinal wave propagation direction is derived as follows:
wherein the content of the first and second substances,is a characteristic angular frequency, P0Is a die of pressure.
Fluid pressure
We believe that the displacement of the solid and the displacement of the fluid in the x-direction are influenced by the average of the pressure and the displacement of the fluid in the r-direction, which is influenced only by the local flow in the spherical Biot flow. Thus, the solid displacement and the displacement of the fluid in the x direction are set to be:
u(x,t)=C1ei(lx-ωt), (10a)
w(x,t)=C2ei(lx-ωt), (10b)
wherein, C1And C2Is a constant.
Substituting equations (7) - (9) into equation (5) yields the following relationship between fluid and pressure:
solving the equation using the constant voltage boundary results in:
wherein, J0Is a function of the zero order bessel function,
(λReis an intermediate variable)
The average fluid pressure can be derived from equations (8a), (8b), (10a) and (10 b):
wherein, J1Is a first order bessel function.
The partial derivative is calculated from the above equation (11), so that the following relation can be obtained:
assuming the relationship:
(FReis an intermediate variable)
Step S5, predicting phase velocity and inverse quality factor Q by using Re-BISQ model-1:
Average local fluid pressure PavThe actual fluid pressure kinetic equation (1b) can be substituted for (Dvorkin and Chur, 1993). From the relationship between (Toksoz and Johnston,1981) wave number and wave velocity and attenuation, the velocity (V) of the longitudinal wave is obtainedp) And the expression of attenuation factor (α). The specific derivation process is described in (Dvorkin and Nur,1993) appendix。
a1,2=ωIm(X1,2) (16)
Where ω is the angular velocity, ρ1=(1-φ)ρs,ρ2=φρf,ρ1,ρ2Are all density constants.
Expression (16) is expressed in inverse quality factor form as:
note: except for parameters with a definite meaning, other parameters are intermediate variables.
Step S6, predicting phase velocity and inverse quality factor Q by using BISQ model-1:
The expressions of the longitudinal wave velocity and the attenuation (i.e., the inverse quality factor) are respectively expressed by expressions (15) and (17), and F in expressions (18) and (19)ReAnd F is replaced.
Step S7, result analysis:
the numerical simulation calculation parameters of the present invention are derived from Berryman (1980). The values of all parameters are respectively phi is 0.15, and the volume modulus K of the solid framework isdry=16×109Pa, poisson ratio υ 0.15, matrix density ρs=2650kg/m3,
Matrix bulk modulus K0=38×109Pa, permeability kappa of 1 × 10-15m2Fluid density ρf=1000kg/m3Fluid viscosity η ═ 0.001Pa · s, bulk modulus K of the fluidf=2.25×109Pa, fluid-solid coupling density rhoa=420kg/m3The characteristic jet length R is 1mm, the fracture aspect ratio ζ is 0.02(Mavko et al, 1998; Dongning et al 2014), and the fracture density ε is 0.15. The velocity and attenuation of longitudinal waves are calculated through steps S4-S5, and a graph is drawn.
FIG. 3 shows the fast longitudinal wave velocity (graph a) and attenuation (graph b), and the slow longitudinal wave velocity (graph c) and attenuation (graph d) as a function of frequency predicted by the Biot model, the BISQ model and the Re-BISQ model. The attenuation peak frequencies of the Re-BISQ model and the BISQ model are equal, and under the frequency, the fast longitudinal wave speed and the attenuation calculated by the Re-BISQ model are both larger than the value calculated by the BISQ model; as shown in fig. 3(b), when the calculation result is satisfied with the angular frequency ω → ∞ as a result of the calculation, the fast longitudinal wave velocity and the attenuation value are the same as those calculated by the Biot model. It can be seen that the fast longitudinal wave velocity values calculated by the two models are not greatly different, the calculated value of the Re-BISQ model is slightly larger than the BISQ model value, and the difference value is not more than 1%. In the low frequency band (10)3-105Hz) the calculated attenuation value of the Re-BISQ model is larger than the calculated value of the BISQ model, namely the attenuation value of the low-frequency band is increased. In the high frequency range (greater than 10)5Hz), the attenuation value of the Re-BISQ model is slightly smaller than the value calculated by the BISQ model, namely the attenuation of a high frequency band is reduced, the value of the attenuation is closer to the attenuation value calculated by the Biot model along with the increase of the frequency, and the Re-BISQ model can distinguish squit attenuation from Biot attenuation. This is mainly because, in the low frequency band, due to the elongated fissures contained in the rock, when the seismic waves press against the rock, the fluid in the pores flows from the fissures towards the hard pores (Mavko and Nur,1975), generating jets. The jet stream causes the rock to produce a greater attenuation in the mesoscale range, i.e. the coin-type attenuation value is greater than the sphere-type attenuation value (Pride, 2004). At high frequencies, the fluid flow within the rock is dominated by macroscopic Biot flow, i.e. the fluid pressure equilibrium within the rock occurs in the wavelength scale range (Pride, 2004). Thus, at high frequencies within the rock, macroscopic Biot flow dominates, and at low frequencies, mesoscopic localization dominatesA stream is an squrit stream. Therefore, the Re-BISQ model is superior to the BISQ model and the Biot model in predicting the internal attenuation of the rock.
When the frequency band of the Re-BISQ model is lower, the calculated slow longitudinal wave speed and attenuation are larger than those of the BISQ model and the Biot model, and the difference between the calculated slow longitudinal wave speed and attenuation and the calculated slow longitudinal wave speed and attenuation is very small. And the Re-BSQ is very close to the calculated value of the BISQ model, and the calculated results of the three models are equal when the frequency is higher.
In order to further verify the effectiveness and the correctness of the method provided by the invention, two different experimental data of conventional sandstone and compact sandstone are respectively selected, and the relationship between the prediction results and the experimental values of the two models is analyzed. The following detailed description of specific embodiments of the present invention is made with reference to the accompanying drawings and the embodiments:
example one
In this example, 10 samples of saturated conventional sandstone samples from southampton, uk were compared with theoretical predicted values:
10 pieces of conventional sandstone with the experimental measurement frequency of 1MHz, the confining pressure of 40MPa and the bulk modulus K of pore fluidf=2.25×109Pa, fluid density ρf=1000Kg/m3. The porosity phi is 15 +/-2%, and the permeability is less than 100mD.
Fig. 4 is a graph showing the relationship between the actually measured values (Klimentos and McCann,1990) and the model-predicted longitudinal wave attenuation values, in which the solid line represents the longitudinal wave attenuation values predicted by the Re-BISQ model, the broken line represents the longitudinal wave attenuation values predicted by the BISQ model according to the steps, and the chain line represents the experimentally measured longitudinal wave attenuation values. The prediction result of the Re-BISQ model is closer to the experimental measurement value. Compared with measured values, the prediction result of the method is well matched with experimental data, and the method shows that the method more accurately reflects the attenuation characteristic of the conventional sandstone.
Example two
In this example, the experimental data and the theoretical predicted value of 10 tight sandstone samples in the tondors basin surrige gas field are compared:
experimental setup, experimental conditions, experimental protocol and literature for 10 tight sandstones (wushiyang et al, 2000; wang daxing et al,2006) substantially identical was used. The stratum conditions adopted in the experiment are 105 ℃ and 29 MPa. The experimental measurement frequency was 1 MHz. Bulk modulus K of pore fluidf=2.25×109Pa, fluid density ρf=1000Kg/m3。
Fig. 5 shows a schematic diagram of the relationship between the experimentally measured value (wang da xing, 2016) and the model-predicted longitudinal wave attenuation value, in which the solid line represents the Re-BISQ model-predicted longitudinal wave attenuation value, the broken line represents the BISQ model-predicted longitudinal wave attenuation value according to the procedure, and the dot-dash line represents the experimentally measured longitudinal wave attenuation value. The prediction result of the Re-BISQ model is closer to the experimental measurement value. Compared with measured values, the prediction result of the method is well matched with experimental data, and the method shows that the method more accurately reflects the attenuation characteristic of the compact sandstone.
The invention is researched on the basis of the research results of predecessors that the fracture has an important influence (pride, 2004; dawn, 2011) on the propagation of the elastic wave, a wave propagation theoretical equation is deduced, and the influence of the fracture on the wave propagation is analyzed on sandstone and tight sandstone samples in two different areas, so that the influence of rock heterogeneity on wave dispersion and attenuation is further explained.
The above description is only a preferred embodiment of the present invention, and it should be noted that, for those skilled in the art, several modifications and variations can be made without departing from the technical principle of the present invention, and these modifications and variations should also be regarded as the protection scope of the present invention.
Claims (7)
1. A method for predicting the longitudinal wave velocity and attenuation of a pore medium containing fractures is characterized by comprising the following steps:
step S1, characterizing the internal pore structure of the rock;
step S2, characterizing the wave propagation characteristics in the rock with cracks;
step S3, estimating the proportion of the porosity of the rock fractures to the total porosity;
in step S4, in the case of one dimension, the equation for the oscillation of the fluid in the two-phase medium can be obtained as follows:
where t denotes time, κ and η denote permeability and fluid viscosity, respectively, u and w denote displacement of solid and fluid in the x-direction, respectively, and ρs、ρfAnd ρaRepresenting solids, fluids, and add-on densities; m is the Biot elasticity parameter; p is the fluid pressure;is the pore elasticity parameter; wherein KdryIs a solid skeleton bulk modulus, K0Is the matrix bulk modulus;
when the rock is squeezed, the fluid inside the rock fracture is sprayed to the hard pores, so that a spraying flow is generated, and the following equation can be obtained by combining the conservation of the mass of the fluid:
wherein the content of the first and second substances,the ratio of the total porosity of the fracture, v is the displacement of the fluid in the direction r, and according to David and Zimmerman, the expression for the fracture porosity isZeta is the aspect ratio of the crack, epsilon is the crack density;
the relationship between porosity and framework and fluid is as follows:
substituting (4) and (5) into (3),
wherein the content of the first and second substances,and c is0Being the speed of sound in the liquid, Q and F are intermediate variables; through Lagrange's equation and combining with Biot dissipation function, obtain the fluid dissipation function of the oscillation perpendicular to the propagation direction of the longitudinal wave, namely the fluid dissipation function relation is as follows:
the r-direction fluid displacement expression and the change form of the fracture fluid pressure are respectively set as follows:
v(x,r,t)=v0(r)ei(lx-ωt)(8a)
P(x,r,t)=P0(r)ei(lx-ωt)(8b)
wherein the content of the first and second substances,is the wavenumber, i denotes the imaginary unit, vttRepresenting one derivative of the displacement of the fluid in the r direction, vtRepresenting the secondary derivation of the r-direction fluid displacement; v is0A mode representing a fluid displacement in the r direction;
substituting the formula (8) into the formula (7), the fluid pressure gradient equation in the vertical longitudinal wave propagation direction is derived as follows:
wherein the content of the first and second substances,is a characteristic angular frequency, P0A die that is a press;
the solid displacement and the displacement of the fluid in the x direction are set as follows:
u(x,t)=C1ei(lx-ωt), (10a)
w(x,t)=C2ei(lx-ωt),(10b)
wherein, C1And C2Is a constant;
substituting equations (7) - (9) into equation (5) yields the following relationship between fluid and pressure:
solving the equation using the constant voltage boundary results in:
wherein, J0Is a function of the zero order bessel function,
λReis an intermediate variable;
the average fluid pressure can be derived from equations (8a), (8b), (10a) and (10 b):
wherein, J1Is a first order Bessel function;
the partial derivative is calculated from the above equation (11), so that the following relation can be obtained:
assuming the relationship:
FReis an intermediate variable;
step S5, predicting the longitudinal wave velocity and the inverse quality factor Q by using the Re-BISQ model-1;
Step S6, predicting the longitudinal wave velocity and the inverse quality factor Q by using the BISQ model-1;
Step S7, result analysis.
2. The method for predicting the longitudinal wave velocity and the attenuation of the pore medium containing the fractures according to claim 1, wherein in the step S1, the internal pore structure characteristics of the rock are characterized by utilizing the shape of the rock slice under a microscope.
3. The method for predicting the velocity and attenuation of longitudinal waves of a pore medium containing fractures according to claim 1, wherein in step S2, the propagation characteristics of seismic waves inside rocks are characterized by using the propagation characteristic image of elastic waves in a cylinder sponge.
4. The method for predicting the longitudinal wave velocity and attenuation of a pore medium containing fractures as claimed in claim 1, wherein in the step S3, the expression of the ratio of fracture porosity to total porosity is as follows:
5. The method for predicting the longitudinal wave velocity and the attenuation of the pore medium containing the fissure as claimed in claim 1, wherein the longitudinal wave velocity and the inverse quality factor predicted by the two models BISQ and Re-BISQ are calculated through steps S4-S6.
6. The method as claimed in claim 1, wherein the difference between the two models of BISQ and Re-BISQ in predicting the longitudinal wave velocity and the inverse quality factor is comparatively analyzed by step S7.
7. The method for predicting the longitudinal wave velocity and attenuation of a pore medium containing fractures according to claim 1, wherein in step S5, the longitudinal wave velocity V is obtainedpAnd attenuation factor α expression;
a1,2=ωIm(X1,2) (16)
where ω is the angular velocity, ωc=ηφ/κρf,ρ1=(1-φ)ρs,ρ2=φρf,
Expression (16) is expressed in inverse quality factor form as:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811200996.5A CN109116420B (en) | 2018-10-16 | 2018-10-16 | Method for predicting longitudinal wave velocity and attenuation of pore medium containing fractures |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811200996.5A CN109116420B (en) | 2018-10-16 | 2018-10-16 | Method for predicting longitudinal wave velocity and attenuation of pore medium containing fractures |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109116420A CN109116420A (en) | 2019-01-01 |
CN109116420B true CN109116420B (en) | 2020-02-21 |
Family
ID=64854650
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811200996.5A Active CN109116420B (en) | 2018-10-16 | 2018-10-16 | Method for predicting longitudinal wave velocity and attenuation of pore medium containing fractures |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109116420B (en) |
Families Citing this family (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110276091B (en) * | 2019-04-26 | 2022-10-14 | 中国石油化工股份有限公司 | Elastic wave response model modeling method based on rock multi-pore structure |
CN110456412B (en) * | 2019-07-09 | 2021-06-08 | 河海大学 | Method for identifying fluid saturation of carbonate reservoir based on post-stack seismic data |
CN112630829A (en) * | 2019-10-08 | 2021-04-09 | 中国石油化工股份有限公司 | Method and system for analyzing tight sandstone elastic wave attenuation property |
CN110703322B (en) * | 2019-10-10 | 2020-11-13 | 清华大学 | Wave propagation processing method, device and equipment |
CN111695308B (en) * | 2020-05-22 | 2021-04-27 | 河海大学 | Method for predicting change characteristics of longitudinal wave speed along with temperature change |
CN112505772B (en) * | 2020-12-10 | 2022-05-31 | 中国石油大学(华东) | Method for inverting rock pore distribution characteristics by utilizing pore and fracture medium elastic wave theory |
CN113391345B (en) * | 2021-06-02 | 2022-09-30 | 清华大学 | Method and device for predicting wave propagation characteristics, electronic device, and storage medium |
CN113887150B (en) * | 2021-09-18 | 2022-05-06 | 河海大学 | Method for estimating length of characteristic jet flow of compact sandstone |
CN114021498B (en) * | 2021-11-05 | 2022-10-11 | 中国矿业大学 | High-efficiency numerical simulation method for predicting elastic modulus of multiphase pore medium |
CN114236609B (en) * | 2021-12-17 | 2022-07-19 | 河海大学 | Prediction method for longitudinal wave velocity and attenuation of partially saturated hole fracture medium |
CN114674934B (en) * | 2022-02-18 | 2023-06-13 | 河海大学 | Method for establishing theoretical model of saturated heavy oil rock wave velocity changing along with temperature |
CN116840912A (en) * | 2022-09-13 | 2023-10-03 | 上海电子信息职业技术学院 | Method for predicting longitudinal wave velocity in partially saturated pore medium with different fracture aspect ratios |
CN115774286B (en) * | 2022-12-06 | 2023-09-05 | 上海电子信息职业技术学院 | Longitudinal wave speed prediction method for pore medium containing viscoelastic fluid and application of longitudinal wave speed prediction method |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP3568304B2 (en) * | 1995-12-06 | 2004-09-22 | 株式会社超高温材料研究所 | Analysis method of thermal constant using laser flash method |
CN102508296B (en) * | 2011-11-14 | 2013-07-31 | 中国石油天然气股份有限公司 | Method and device for analyzing dispersion and attenuation of unsaturated double-porosity medium earthquake waves |
CN103412336B (en) * | 2013-07-22 | 2015-12-09 | 中国石油天然气股份有限公司 | The velocity of longitudinal wave Forecasting Methodology of rock system in a kind of heterogeneous reservoir |
CN104360383B (en) * | 2014-11-12 | 2017-03-22 | 中国石油大学(华东) | Method and system for predicting seismic wave attenuation |
US9599743B2 (en) * | 2015-04-29 | 2017-03-21 | Baker Hughes Incorporated | Density measurements using detectors on a pulsed neutron measurement platform |
CN106959469A (en) * | 2017-04-14 | 2017-07-18 | 中国石油天然气股份有限公司 | The speed and decay analog analysing method and device of seismic wave |
-
2018
- 2018-10-16 CN CN201811200996.5A patent/CN109116420B/en active Active
Also Published As
Publication number | Publication date |
---|---|
CN109116420A (en) | 2019-01-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109116420B (en) | Method for predicting longitudinal wave velocity and attenuation of pore medium containing fractures | |
CN109471168B (en) | Method for predicting velocity and attenuation of longitudinal wave in pore fracture medium | |
CN110276091B (en) | Elastic wave response model modeling method based on rock multi-pore structure | |
CN111695308B (en) | Method for predicting change characteristics of longitudinal wave speed along with temperature change | |
CN112505772B (en) | Method for inverting rock pore distribution characteristics by utilizing pore and fracture medium elastic wave theory | |
CN114236609B (en) | Prediction method for longitudinal wave velocity and attenuation of partially saturated hole fracture medium | |
CN110275206A (en) | A kind of crack-pore type rock physics Elastic forming board | |
CN113887150B (en) | Method for estimating length of characteristic jet flow of compact sandstone | |
Wang et al. | Modelling of low-frequency acoustic wave propagation in dilute gas-bubbly liquids | |
Santos et al. | Ultrasonic spectral and complexity measurements on brine and oil saturated rocks | |
Sadouki | RETRACTED: Experimental characterization of human cancellous bone via the first ultrasonic reflected wave–Application of Biot’s theory | |
Li et al. | Pulsating pressurization of two-phase fluid in a pipe filled with water and a little gas | |
Shagapov et al. | On the Theory of Acoustic Sounding of Hydraulic-Fracturing Cracks Perpendicular to the Well | |
CN115774286B (en) | Longitudinal wave speed prediction method for pore medium containing viscoelastic fluid and application of longitudinal wave speed prediction method | |
Du et al. | Wave propagation characteristics in frozen saturated soil | |
CN112649854A (en) | Seismic wave frequency dispersion and attenuation prediction method and system based on dual-scale model | |
CN116522810B (en) | Reservoir frequency dispersion attenuation analysis method based on Norris-KG equivalent medium modeling | |
Wu et al. | Anelasticity of porous rocks containing microcracks: a reformulation of the BISQ model. | |
Poesio et al. | Removal of particle bridges from a porous material by ultrasonic irradiation | |
CN116840912A (en) | Method for predicting longitudinal wave velocity in partially saturated pore medium with different fracture aspect ratios | |
Wu et al. | A patchy-saturated rock physics model for tight sandstone based on microscopic pore structures | |
CN112630829A (en) | Method and system for analyzing tight sandstone elastic wave attenuation property | |
He et al. | Laboratory experiments and theoretical study of pressure and fluid influences on acoustic response in tight rocks with pore microstructure | |
CN111597763B (en) | Cross multi-scale flow pore medium full-band elastic wave dispersion attenuation analysis method | |
Tafkirte et al. | A numerical investigation of longitudinal wave for thermo-acoustoelastic effects on fluid-embedded in two viscoelastic layer plates at higher temperatures |
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 |