CN109298443A - Two-fluid injection stream model elastic wave evaluation method and computer readable storage medium - Google Patents
Two-fluid injection stream model elastic wave evaluation method and computer readable storage medium Download PDFInfo
- Publication number
- CN109298443A CN109298443A CN201710612448.2A CN201710612448A CN109298443A CN 109298443 A CN109298443 A CN 109298443A CN 201710612448 A CN201710612448 A CN 201710612448A CN 109298443 A CN109298443 A CN 109298443A
- Authority
- CN
- China
- Prior art keywords
- fluid
- wave
- equation
- model
- frequency
- 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.)
- Granted
Links
- 239000012530 fluid Substances 0.000 title claims abstract description 178
- 238000002347 injection Methods 0.000 title claims abstract description 61
- 239000007924 injection Substances 0.000 title claims abstract description 61
- 238000011156 evaluation Methods 0.000 title claims abstract description 16
- 239000011148 porous material Substances 0.000 claims abstract description 27
- 239000006185 dispersion Substances 0.000 claims abstract description 14
- 238000006073 displacement reaction Methods 0.000 claims description 41
- 239000007787 solid Substances 0.000 claims description 10
- 238000010008 shearing Methods 0.000 claims description 8
- 230000035699 permeability Effects 0.000 claims description 7
- 238000004590 computer program Methods 0.000 claims description 3
- 230000007246 mechanism Effects 0.000 abstract description 28
- 239000011435 rock Substances 0.000 abstract description 15
- 238000000034 method Methods 0.000 abstract description 13
- 230000021715 photosynthesis, light harvesting Effects 0.000 abstract description 3
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 12
- 239000007921 spray Substances 0.000 description 11
- 238000002474 experimental method Methods 0.000 description 5
- 229920006395 saturated elastomer Polymers 0.000 description 5
- 230000008901 benefit Effects 0.000 description 4
- 230000000052 comparative effect Effects 0.000 description 4
- 230000006835 compression Effects 0.000 description 4
- 238000007906 compression Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 4
- 230000000694 effects Effects 0.000 description 4
- 230000008878 coupling Effects 0.000 description 3
- 238000010168 coupling process Methods 0.000 description 3
- 238000005859 coupling reaction Methods 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 239000007790 solid phase Substances 0.000 description 3
- 239000000243 solution Substances 0.000 description 3
- 230000008859 change Effects 0.000 description 2
- 238000005213 imbibition Methods 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 239000011159 matrix material Substances 0.000 description 2
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical group 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 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 150000001875 compounds Chemical class 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 239000013078 crystal Substances 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 239000010433 feldspar Substances 0.000 description 1
- 230000004927 fusion Effects 0.000 description 1
- 230000002070 germicidal effect Effects 0.000 description 1
- 230000003455 independent Effects 0.000 description 1
- 229910052500 inorganic mineral Inorganic materials 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 239000007788 liquid Substances 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 239000011707 mineral Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 239000004576 sand Substances 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
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/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
- G01V1/303—Analysis for determining velocity profiles or travel times
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)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
Disclose a kind of two-fluid injection stream model elastic wave evaluation method and computer readable storage medium.This method may include: to establish pore media kinetics equation based on single fluid, single hole porosity injection flow model, and combined hole clearance flow weight conservation equation obtains Fluid conservation equation, and then establishes two-fluid injection flow model；Flow model is sprayed based on two-fluid, derivative and average pressure by pore media kinetics equation, wave field pressure variations about radius obtain the equation of higher order of frequency and wave number；The equation of higher order of frequency and wave number is solved in frequency domain, obtains dispersion curve and attenuation curve, and then analyzes influence of the two-fluid to elastic wave.The present invention establishes while the petrophysical model comprising macroscopic view-Jie's sight-microcosmic three scales fractional saturation double porosity media and patch injection stream mechanism, realize integrated interpretation wavefield energy dissipation mechanism, and based on this spread speed of model prediction seismic wave in porous rocks medium.
Description
Technical field
The present invention relates to oil gas technical field of physical geography, more particularly, to a kind of two-fluid injection stream model elasticity
Wave evaluation method and computer readable storage medium.
Background technique
In unconventional oil and gas resource exploration, people often face subsurface reservoir rock pore structure complexity, pore-fluid
Fractional saturation, dynamic response rule difficult points, the applicability of common seismic theory of wave propagation such as non-linear face the challenge, thus
It brings the explanation of seismic observation data difficult, brings and seriously affect to Exploration of Oil And Gas.In order to consider pore-fluid to rock
The influence of middle seismic wave, Biot establish the theory of isotropism pore media Propagation of Elastic Wave for the first time, predict hole for the first time
There is three kinds of waves in gap medium i.e.: shear wave, fast longitudinal wave and Slow P-wave.Biot theory has attracted numerous scholars to put into after occurring
The research of the theory, using.The frequency dispersion and decaying of the wave caused in Biot theory predictive of macroscopic fluid flow, but use
When this macroscopic fluid flow mechanism estimates the wave attenuation in supersonic band pore media with dispersion phenomenon work, discovery is many
In the case of prediction result it is more obvious than experimental result relatively low.Many scholars improve Biot model, wherein microcosmic jet flow
(Squirt) mechanism obtains the approval of numerous scholars, and is used to explain the velocity dispersion of seismic wave and decaying in high Mid Frequency
Phenomenon.(Dvorkin J, the Nur A.Dynamic poroelasticity:A unified model with such as Dvorkin
The squirt and the Biot mechanisms.Geophysics, 1993,58 (4): 524-533.) propose one kind
New poroelasticity dynamic model is theoretical.The theory combines the Biot mechanism in macro-scale and the Squirt on micro-scale
Mechanism, therefore referred to as BISQ is theoretical.Relative to Biot theory, the maximum of the theory is the difference is that it not only allows for flowing
Body and wave propagate the fluid on the same direction and flow (Biot mechanism), it is also contemplated that fluid and wave propagate the flowing in vertical direction
(Squirt mechanism).BISQ model have in terms of the rock Experimental Ultrasonic data explanation of laboratory relative to Biot pore model compared with
It is big to improve, but it is still relatively low compared to actual measurement in the seismic velocity of the earthquake frequency range theoretical prediction, decaying.Work as compressional wave
When passing through multiphase porous media, on meso-scale (the heterogeneity characteristic dimension i.e. in pore media is greater than pore-size,
But it is much smaller than wavelength dimension) and local fluid flow can be generated, so as to cause elastic wave energy decaying and velocity dispersion.Jie's sight office
The phenomenon that portion's flow mechanism causes largely appears in sedimentary rock, can explain many important phenomenons, so as to cause researcher
Concern.White(White,J.E.,Computed Seismic Speeds and Attenuation in Rocks with
Partial Gas Saturation, 1975, Geophysics, 40 (2), 224-232.) to establish first meso-scale non-
The fluid flow model of uniform pores rock, the model are the part water saturation pore models that a kind of inside includes bubble, also by
Referred to as patchy saturation.Pride and Berryman (Pride, S.R., and J.G.Berryman, Linear dynamics
of double-porosity dual-permeability materials.I.Governing equations and
Acoustic attenuation, 2003, Phys Rev E, 68 (3), 036603) to have developed two kinds of different aperture rocks compound
The patchy saturation of body, Sun and Ba, Carcione (Sun, W., J.Ba, and J.M.Carcione, Theory of wave
propagation in partially saturated double-porosity rocks:a triple-layer
Patchy model, Geophysical Journal International, 2016,205 (1), 22-37.) it proposes simultaneously
Three layers of patchy saturation comprising solid skeletal and the heterogeneous influence of pore-fluid, make velocity of wave theoretical expectation values and multiple frequency ranges
Experimental data all realize well coincide.But this kind of theories often do not account for the non-of pore media skeleton substance itself
The influence of uniform property, especially double porosity and double elastic attribute, it is non-that obtained final result is also independent from reservoir
Even property feature, can not fundamentally reflect the complex wave field feature of pore media.In the research of the mechanics of hole containing fluid, macroscopic view
The Biot of scale is theoretical, Jie sees patchy saturation and the injection flow model of micro-scale is the important power of solid stream interaction
Mechanism.But people often use these three theoretical treatment problems respectively, up to the present, macroscopical Biot is theoretical, Jie sees
Across scale unified model of the fractional saturation Patchy model together with microcosmic injection stream Model Fusion does not occur also.Therefore, having must
Develop a kind of two-fluid injection stream model elastic wave evaluation method.
The information for being disclosed in background of invention part is merely intended to deepen the reason to general background technique of the invention
Solution, and it is known to those skilled in the art existing to be not construed as recognizing or imply that the information is constituted in any form
Technology.
Summary of the invention
The invention proposes a kind of two-fluid injection stream model elastic wave evaluation method and computer readable storage medium,
By establishing while stream mechanism can be sprayed comprising macroscopic view-Jie's sight-microcosmic three scales fractional saturation double porosity media and patch
Petrophysical model, realize integrated interpretation wavefield energy dissipation mechanism, and based on this model proposition predict seismic wave in hole
The method of spread speed in gap rock medium.
A kind of two-fluid injection stream model elastic wave evaluation method is proposed according to the present invention.The method may include:
Flow model is sprayed based on single fluid, single hole porosity, establishes pore media kinetics equation, combined hole clearance flow weight conservation side
Journey obtains Fluid conservation equation, and then establishes two-fluid injection flow model；Flow model is sprayed based on the two-fluid, passes through institute
The derivative and average pressure of pore media kinetics equation, wave field pressure variations about radius are stated, the height of frequency and wave number is obtained
Rank equation；The equation of higher order of the frequency and wave number is solved in frequency domain, obtains dispersion curve and attenuation curve, and then is analyzed double
Influence of the fluid to elastic wave.
Preferably, the Fluid conservation equation are as follows:
Wherein, p1,p2It is the pressure in two kinds of fluids, ZxIt is the local flow between Immiscible fluid patch, uxIt is solid
Displacement, UxIt is displacement of fluid, UrIt is radial fluid displacement, φ1=v1φ10,φ2=v2φ20, v1,v2It is body shared by two kinds of fluids
Product ratio, φ10,φ20It is two kinds of porositys present in medium,
Preferably, derivative of the wave field pressure variations about radius are as follows:
Wherein, p01,p02,It is the amplitude of wave field variable, p01,p02For wave pressure variable,For fluid
Density, ω are angular frequency, and k is wave number,μmIt is fluid viscosity coefficient, κmIt is skeleton permeability,
M=1,2.
Preferably, the average pressure are as follows:
Wherein,For average pressure,R is the radius of hole runner.
Preferably, the equation of higher order of the frequency and wave number are as follows:
Wherein, b23=b32=0, Wherein,Gm,νmIt is modulus of shearing and Poisson's ratio, m=1,2.
A kind of computer readable storage medium is proposed according to the present invention, is stored thereon with computer program, wherein described
It is performed the steps of when program is executed by processor based on single fluid, single hole porosity injection flow model, establishes pore media power
Equation is learned, combined hole clearance flow weight conservation equation obtains Fluid conservation equation, and then establishes two-fluid injection flow model；Base
Flow model is sprayed in the two-fluid, passes through the pore media kinetics equation, derivative of the wave field pressure variations about radius
With average pressure, the equation of higher order of frequency and wave number is obtained；The equation of higher order of the frequency and wave number is solved in frequency domain, is obtained
Dispersion curve and attenuation curve, and then analyze influence of the two-fluid to elastic wave.
Preferably, the Fluid conservation equation are as follows:
Wherein, p1,p2It is the pressure in two kinds of fluids, ZxIt is the local flow between Immiscible fluid patch, uxIt is solid
Displacement, UxIt is displacement of fluid, UrIt is radial fluid displacement, φ1=v1φ10,φ2=v2φ20, v1,v2It is body shared by two kinds of fluids
Product ratio, φ10,φ20It is two kinds of porositys present in medium,
Preferably, derivative of the wave field pressure variations about radius are as follows:
Wherein, p01,p02,It is the amplitude of wave field variable, p01,p02For wave pressure variable,For fluid
Density, ω are angular frequency, and k is wave number,μmIt is fluid viscosity coefficient, κmIt is skeleton permeability,
M=1,2.
Preferably, the average pressure are as follows:
Wherein,For average pressure,R is the radius of hole runner.
Preferably, the equation of higher order of the frequency and wave number are as follows:
Wherein, b23=b32=0, Wherein,Gm,νmIt is modulus of shearing and Poisson's ratio, m=1,2.
The beneficial effect is that: (1) microcosmic injection stream model mechanism fractional saturation patch in is introduced, while in patch
Still retain local flow mechanism on interface, breach the one fluid model of conventional spray stream, in fractional saturation double porosity media
Injection stream mechanism is introduced, two-fluid injection flow model is obtained；(2) by Patchy model wave equation (macroscopic view-Jie sees) and injection stream
Mechanism (microcosmic) combines, and realizes the wave equation that macroscopic view-Jie's sight-micromechanism combines, breaches conventional spray stream mould
The type limitations with vertical two displacements in parallel introduce two kinds of fluid boundary local flows displacements, different fluid compression may be implemented
The coupling of rate.
Methods and apparatus of the present invention has other characteristics and advantages, these characteristics and advantages are attached from what is be incorporated herein
It will be apparent in figure and subsequent specific embodiment, or will be in the attached drawing being incorporated herein and subsequent specific reality
It applies in mode and is stated in detail, the drawings and the detailed description together serve to explain specific principles of the invention.
Detailed description of the invention
Exemplary embodiment of the present is described in more detail in conjunction with the accompanying drawings, of the invention is above-mentioned and other
Purpose, feature and advantage will be apparent, wherein in exemplary embodiments of the present invention, identical reference label is usual
Represent same parts.
Fig. 1 shows the flow chart of the step of two-fluid injection stream model elastic wave evaluation method according to the present invention.
Fig. 2 a and Fig. 2 b respectively illustrate the signal of single fluid, single hole porosity injection flow model and two-fluid injection flow model
Figure.
Fig. 3 show the velocity of longitudinal wave prediction according to an embodiment of the invention to French Vosgian sandstone with
The schematic diagram of the Comparative result of experimental observation number.
Fig. 4 shows the velocity of longitudinal wave prediction and experiment according to an embodiment of the invention to Fort Union sandstone
Observe the schematic diagram of the Comparative result of number.
Specific embodiment
The present invention will be described in more detail below with reference to accompanying drawings.Although showing the preferred embodiment of the present invention in attached drawing,
However, it is to be appreciated that may be realized in various forms the present invention and should not be limited by the embodiments set forth herein.On the contrary, providing
These embodiments are of the invention more thorough and complete in order to make, and can will fully convey the scope of the invention to ability
The technical staff in domain.
Fig. 1 shows the flow chart of the step of two-fluid injection stream model elastic wave evaluation method according to the present invention.
In this embodiment, two-fluid injection stream model elastic wave evaluation method according to the present invention may include: step
101, flow model is sprayed based on single fluid, single hole porosity, establishes pore media kinetics equation, combined hole clearance flow weight conservation
Equation obtains Fluid conservation equation, and then establishes two-fluid injection flow model；Step 102, flow model is sprayed based on two-fluid,
Derivative and average pressure by pore media kinetics equation, wave field pressure variations about radius, acquisition frequency and wave number
The equation of higher order；Step 103, the equation of higher order of frequency and wave number is solved in frequency domain, obtains dispersion curve and attenuation curve, in turn
Analyze influence of the two-fluid to elastic wave.
The embodiment is by establishing while comprising macroscopic view-Jie's sight-microcosmic three scales fractional saturation double porosity media and spot
Block sprays the petrophysical model of stream mechanism, realizes integrated interpretation wavefield energy dissipation mechanism, and propose prediction based on this model
The method of spread speed of the seismic wave in porous rocks medium.
The following detailed description of the specific steps of the elastic wave evaluation method of two-fluid according to the present invention injection flow model.
In one example, flow model is sprayed based on single fluid, single hole porosity, establishes pore media kinetics equation, tied
Pore-fluid mass-conservation equation is closed, Fluid conservation equation is obtained, and then establishes two-fluid injection flow model.
Fig. 2 a and Fig. 2 b respectively illustrate the signal of single fluid, single hole porosity injection flow model and two-fluid injection flow model
Figure.
Specifically, as shown in Figure 2 a, it is to exist in hole runner that single fluid, single hole porosity, which spray the physical form of flow model,
A kind of fluid, there are two types of fluxus formaes for this fluid tool: the global flow along runner axial direction and the injection stream along runner radial direction.Such as
Shown in Fig. 2 b, it is in hole runner there are two kinds of fluids that two-fluid, which sprays the physical form of flow model, and this fluid is immiscible,
With interface, therefore there are three types of fluxus formaes for tool: along the global flow of runner axial direction, along the injection stream and fluid circle of runner radial direction
The local flow in face moves.Single fluid, single hole porosity injection flow model have 2 independents variable of axial displacement and radial displacement, these
Independent variable meets hydrodynamic equations group；Two-fluid injection flow model has more the local flow displacement an of fluid boundary, because
This introduces the kinetics equation that local flow displacement meets other than the kinetics equation group that single fluid injection flow model meets
Group matches equation group and the number of arguments, and local flow displacement meets linear displacement principle of stacking, root with other two displacements
According to the kinetics equation relational expression that flow field total displacement meets, flow model can be sprayed from single fluid and export two-fluid injection stream
Model.Flow model is sprayed based on single hole porosity, single fluid, injection stream mechanism is introduced inside fractional saturation patch, along fluctuation
Two kinds of fluids, respectively fluid 1 and fluid 2 are introduced in the flowing in direction, both fluids are parallel to wave field in addition to being respectively provided with
Direction and vertically and outside the flowing in wave field direction, considers the different compression ratios of two kinds of fluids, occurs on two kinds of fluid interface faces
Traditional single fluid injection flow equation is generalized to two-fluid situation by one third variable.
In one example, Fluid conservation equation are as follows:
Wherein, p1,p2It is the pressure in two kinds of fluids, ZxIt is the local flow between Immiscible fluid patch, uxIt is solid
Displacement, UxIt is displacement of fluid, UrIt is radial fluid displacement, φ1=v1φ10,φ2=v2φ20, v1,v2It is body shared by two kinds of fluids
Product ratio, φ10,φ20It is two kinds of porositys present in medium,
Specifically, flow model is sprayed based on two-fluid, according to pore media mechanical knowledge, axial displacement kinetics equation table
It is shown as:
Wherein,M is solid skeletal coefficient of elasticity, a=1-K/Ks, K is skeleton bulk modulus, KsIt is solid phase
Bulk modulus.β*=Q2/Q1.Here KfIt is fluid modulus, Rm,QmIt (m=1,2) is elasticity
Coefficient, φ10,φ20It is two kinds of porositys present in medium.p1,p2It is the pressure in two kinds of fluids, quality coefficient meets Here ρfIt is fluid density, ρ0
It is solid phase density.r10It is saturated with fluid patch
Radius.
Radial displacement kinetics equation is expressed as formula (2):
Wherein, UrIndicate radial fluid displacement,Indicate the density of fluid 1,Indicate the density of fluid 2,It indicates
Second derivative of 2 radial displacement of fluid about the time,It is first order derivative of 1 radial displacement of fluid about the time,It indicates
Second derivative of 1 radial displacement of fluid about the time,It is first order derivative of 2 radial displacement of fluid about the time.
Pore-fluid mass-conservation equation is expressed as formula (3):
Wherein,It is x directional flow change rate,It is radial flow change rate, t is the time.
Formula (3) linearisation is obtained into formula (4):
Wherein, ZxIt is the local flow between Immiscible fluid patch.According to It obtains
Porosity, which is obtained, about the derivative of time isThe differential form of fluid density isIt willWithIt substitutes into formula (4), it is formula (5) that abbreviation, which obtains Fluid conservation equation:
Wherein, c0It is fluid acoustic speed, KfIt is fluid modulus, KsIt is solid phase bulk modulus, K is solid skeletal volume mould
Amount.Rm,Qm, MmIt is coefficient of elasticity.
According to axial displacement kinetics equation (1), radial displacement kinetics equation (2) and Fluid conservation equation (5), simultaneous
The kinetics equation of saturation double porosity media is obtained, as two-fluid sprays flow model.
In one example, flow model is sprayed based on two-fluid, passes through pore media kinetics equation, wave field pressure variations
About the derivative and average pressure of radius, the equation of higher order of frequency and wave number is obtained.
In one example, derivative of the wave field pressure variations about radius are as follows:
Wherein, p01,p02,It is the amplitude of wave field variable, p01,p02For wave pressure variable,For fluid
Density, ω are angular frequency, and k is wave number,μmIt is fluid viscosity coefficient, κmIt is skeleton permeability,
M=1,2.
Specifically,ux,Zx,p1,p2Plane wave solution can be expressed as eι(ωt-k·x)Form, i.e.,
For formula (6):
Wherein, ω, k are frequency and wave number, C1,C3,p01,p02It is the amplitude of wave field variable,
The derivative that wave field pressure variations are then obtained about radius meets (7):
In one example, average pressure are as follows:
Wherein,For average pressure,R is the radius of hole runner.
In one example, the equation of higher order of frequency and wave number are as follows:
Wherein, b23=b32=0, Wherein,Gm,νmIt is modulus of shearing and Poisson's ratio, m=1,2.
In one example, the equation of higher order of frequency and wave number is solved in frequency domain, obtains dispersion curve and attenuation curve,
And then analyze influence of the two-fluid to elastic wave.
Specifically, obtaining Bessel equation by formula (7) is formula (8):
The solution of formula (8) is formula (9):
Wherein, J0Indicate Bessel function,
Then Fluid pressure is formula (10):
The average pressure of fluid channel is formula (11):
If wave field variable meets plane wave form, average pressure and wave field pressure variations are brought into public affairs about the derivative of radius
Formula (1), (2), (4) obtain formula (12):
According to plane wave analysis, by ei(ωt-kx)It substitutes into formula (12), the equation of higher order for obtaining frequency and wave number is formula
(13):
Wherein, b23=b32=0, Wherein,Gm,νmIt is modulus of shearing and Poisson's ratio, m=1,2.
Complex velocity is defined asVelocity of longitudinal wave frequency dispersion is vP=Re (v), attenuation of P-wave are expressed as
And becauseSo the dispersion curve and attenuation curve of the velocity of longitudinal wave in blowhole medium are acquired according to (13),
By observing the variation of wavefield velocity, influence of the two-fluid to elastic wave is analyzed.
This method introduces microcosmic injection stream model mechanism in fractional saturation patch, while still protecting on patch interface
Local flow mechanism is stayed, the one fluid model of conventional spray stream is breached, injection stream machine is introduced in fractional saturation double porosity media
System obtains two-fluid injection flow model；Patchy model wave equation (macroscopic view-Jie sees) and injection stream mechanism (microcosmic) are combined
Come, realize the wave equation that combines of macroscopic view-Jie's sight-micromechanism, breaches conventional spray flow model in parallel and vertical two
The limitation of displacement introduces two kinds of fluid boundary local flow displacements, the coupling of different fluid compression ratio may be implemented.
Two concrete application examples are given below in the scheme and its effect of the embodiment of the present invention for ease of understanding.This field
It should be understood to the one skilled in the art that the example is only for the purposes of understanding the present invention, any detail is not intended to be limited in any way
The system present invention.
Using example 1
This application example using Bacri and Salin (1986) obtain by the water saturated French Vosgian sand of oil
Rock ultrasonic velocity [Bacri and Salin, 1986], the sandstone porosity 21%, passes through imbibition and fluid-discharge method realizes liquid
Body saturation.In fluid-discharge method, sandstone is fully saturated by water, and then oil is injected into sample, and oil content is at 33% or so
Stop increasing, this is the water saturation influence that is bound.In imbibition method, then the fully saturated oil of sample injects water in sample,
Residual oil saturation until reaching 35%, the frequency of acoustic velocity measutement are 350 kilo hertzs.Bacri and Salin[Bacri and
Salin, 1986] rock parameter is as shown in table 1.
Table 1
Fig. 3 show the velocity of longitudinal wave prediction according to an embodiment of the invention to French Vosgian sandstone with
The schematic diagram of the Comparative result of experimental observation number, wherein Sw indicates water saturation, and Vp represents velocity of longitudinal wave, and unit is meter per second
(m/s).For French Vosgian sandstone frequency range ultrasonic experiment data, different injection stream features has been respectively adopted
Scale predicts velocity of longitudinal wave, as a result as shown in Figure 3.When characteristic dimension is smaller (0.01mm, 0.1mm), predetermined speed curve
Almost it is overlapped；(10mm, 50mm) predetermined speed curve co-insides when characteristic dimension is larger；Predetermined speed of larger characteristic dimension is bent
Line and experimental data are coincide fine.
When characteristic dimension is larger, the influence of jet flow gradually weakens, this shows jet flow in pore media at this time
It is less obvious, mainly based on macroscopic view flowing (Biot model) and local flow (Patchy model) mechanism.Meanwhile in current frequency
Under (350kHz), since frequency is higher, pore-fluid is in non-relaxed state, and jet flow can not occur on a large scale at this time, indulge
Wave velocity is gradually consistent with Biot model in high frequency, this also illustrates why the predicted value under larger jet flow size with
True value coincide well, and predetermined speed is almost unrelated with injection stream size at this time.
Using example 2
This application example observes data using the velocity of wave of the low hole unsaturation sandstone of low frequency (Fort Union sandstone)
(Murphy, JOURNAL OF GEOPHYSICAL RESEARCH, 1984), the result predicted with the present invention compare and analyze.
Rock parameter are as follows: matrix volume modulus 35GPa, skeleton bulk modulus 7.14Gpa, skeleton modulus of shearing 9.06Gpa, water body product module
Measure 2.25GPa, volume of air modulus 0.8MPa, water viscosity 0.001Pa*s, air viscosity 0.00001Pa*s, matrix averag density
2.65g/cm3, water density 0.997g/cm3, atmospheric density 0.1g/cm3, porosity 0.085, permeability 0.5mD.Fort
Union sandstone has the feature that (1) includes more loose material (i.e. feldspar and carg), and content is more than 25%；
(2) with the interstitial space or mineral rubble of more opening；(3) sandstone crystal grain diameter is between 0.125 to 0.15 millimeter.
Fig. 4 shows the velocity of longitudinal wave prediction and experiment according to an embodiment of the invention to Fort Union sandstone
Observe the schematic diagram of the Comparative result of number, wherein Sw indicates water saturation, and Vp represents velocity of longitudinal wave, and unit is meter per second.It adopts
With this method, prediction is carried out to Fort Union sandstone low-frequency range (5kHz) sound wave experiment data and has been calculated, discovery this item hair
Bright prediction of speed and Germicidal efficacy coincide preferably, and when between water saturation 20%-80%, experimental data is distributed in feature spray
The prediction of speed range of jet stream size 1mm-10mm, jet flow effect is obvious at this time；When water saturation is greater than 80%, experiment
Data variation is more violent, and when close to 100%, velocity of wave is steeply risen, and the prediction curve of injection stream size 10mm is more close
Experimental data, this shows that jet flow effect decreases when being transitioned into full saturation from fractional saturation；Injection stream size 0.01mm
It is overlapped with the prediction curve of 0.1mm, and predetermined speed is far below experimental data, this shows when injection stream is small-sized, injection
It is almost unrelated with injection stream size to flow effect.
Date comprision shows that this method has the characteristics that reflect that jet flow is strong and weak under different frequency and saturation degree,
Internal flow feature of the open rock under different fluid saturation can be more accurately portrayed, is had compared with other conventional methods bright
Aobvious advantage.
In conclusion introducing microcosmic injection stream model mechanism in fractional saturation patch, while on patch interface still
So retain local flow mechanism, breach the one fluid model of conventional spray stream, injection is introduced in fractional saturation double porosity media
Stream mechanism obtains two-fluid injection flow model；By Patchy model wave equation (macroscopic view-Jie sees) and (microcosmic) knot of injection stream mechanism
Altogether, it realizes the wave equation that macroscopic view-Jie's sight-micromechanism combines, it is in parallel and vertical to breach conventional spray flow model
The limitation of two displacements introduces two kinds of fluid boundary local flow displacements, the coupling of different fluid compression ratio may be implemented.
A kind of computer readable storage medium is proposed according to the present invention, is stored thereon with computer program, wherein described
It is performed the steps of when program is executed by processor based on single fluid, single hole porosity injection flow model, establishes pore media power
Equation is learned, combined hole clearance flow weight conservation equation obtains Fluid conservation equation, and then establishes two-fluid injection flow model；Base
Flow model, derivative and mean pressure by pore media kinetics equation, wave field pressure variations about radius are sprayed in two-fluid
Power obtains the equation of higher order of frequency and wave number；The equation of higher order of frequency and wave number is solved in frequency domain, dispersion curve is obtained and declines
Subtract curve, and then analyzes influence of the two-fluid to elastic wave.
In one example, Fluid conservation equation are as follows:
Wherein, p1,p2It is the pressure in two kinds of fluids, ZxIt is the local flow between Immiscible fluid patch, uxIt is solid
Displacement, UxIt is displacement of fluid, UrIt is radial fluid displacement, φ1=v1φ10,φ2=v2φ20, v1,v2It is body shared by two kinds of fluids
Product ratio, φ10,φ20It is two kinds of porositys present in medium,
In one example, derivative of the wave field pressure variations about radius are as follows:
Wherein, p01,p02,It is the amplitude of wave field variable, p01,p02For wave pressure variable,For fluid
Density, ω are angular frequency, and k is wave number,μmIt is fluid viscosity coefficient, κmIt is skeleton permeability,
M=1,2.
In one example, average pressure are as follows:
Wherein,For average pressure,R is the radius of hole runner.
In one example, the equation of higher order of frequency and wave number are as follows:
Wherein, b23=b32=0, Wherein,Gm,νmIt is modulus of shearing and Poisson's ratio, m=1,2.
Various embodiments of the present invention are described above, above description is exemplary, and non-exclusive, and
It is not limited to disclosed each embodiment.Without departing from the scope and spirit of illustrated each embodiment, for this skill
Many modifications and changes are obvious for the those of ordinary skill in art field.
Claims (10)
1. a kind of two-fluid injection stream model elastic wave evaluation method, comprising:
Flow model is sprayed based on single fluid, single hole porosity, establishes pore media kinetics equation, combined hole clearance flow weight conservation
Equation obtains Fluid conservation equation, and then establishes two-fluid injection flow model；
Flow model is sprayed based on the two-fluid, by the pore media kinetics equation, wave field pressure variations about radius
Derivative and average pressure, obtain frequency and wave number the equation of higher order；
The equation of higher order of the frequency and wave number is solved in frequency domain, obtains dispersion curve and attenuation curve, and then analyze double fluid
Influence of the body to elastic wave.
2. two-fluid injection stream model elastic wave evaluation method according to claim 1, wherein the Fluid conservation equation
Are as follows:
Wherein, p1,p2It is the pressure in two kinds of fluids, ZxIt is the local flow between Immiscible fluid patch, uxIt is solid displacement,
UxIt is displacement of fluid, UrIt is radial fluid displacement, φ1=v1φ10,φ2=v2φ20, v1,v2It is volume ratio shared by two kinds of fluids,
φ10,φ20It is two kinds of porositys present in medium,
3. two-fluid injection stream model elastic wave evaluation method according to claim 2, wherein the wave field pressure variations
Derivative about radius are as follows:
Wherein, p01,p02,It is the amplitude of wave field variable, p01,p02For wave pressure variable, ρf1,ρf2For fluid density,
ω is angular frequency, and k is wave number,μmIt is fluid viscosity coefficient, κmIt is skeleton permeability, m=1,
2。
4. two-fluid injection stream model elastic wave evaluation method according to claim 3, wherein the average pressure are as follows:
Wherein,For average pressure,R is the radius of hole runner.
5. two-fluid injection stream model elastic wave evaluation method according to claim 4, wherein the frequency and wave number
The equation of higher order are as follows:
Wherein, Wherein,Gm,νmIt is modulus of shearing and Poisson's ratio, m=1,2.
6. a kind of computer readable storage medium, is stored thereon with computer program, wherein when described program is executed by processor
It performs the steps of
Flow model is sprayed based on single fluid, single hole porosity, establishes pore media kinetics equation, combined hole clearance flow weight conservation
Equation obtains Fluid conservation equation, and then establishes two-fluid injection flow model；
Flow model is sprayed based on the two-fluid, by the pore media kinetics equation, wave field pressure variations about radius
Derivative and average pressure, obtain frequency and wave number the equation of higher order；
The equation of higher order of the frequency and wave number is solved in frequency domain, obtains dispersion curve and attenuation curve, and then analyze double fluid
Influence of the body to elastic wave.
7. computer readable storage medium according to claim 6, wherein the Fluid conservation equation are as follows:
Wherein, p1,p2It is the pressure in two kinds of fluids, ZxIt is the local flow between Immiscible fluid patch, uxIt is solid displacement,
UxIt is displacement of fluid, UrIt is radial fluid displacement, φ1=v1φ10,φ2=v2φ20, v1,v2It is volume ratio shared by two kinds of fluids,
φ10,φ20It is two kinds of porositys present in medium,
8. computer readable storage medium according to claim 7, wherein wave field pressure variations the leading about radius
Number are as follows:
Wherein, p01,p02,It is the amplitude of wave field variable, p01,p02For wave pressure variable, ρf1,ρf2For fluid density,
ω is angular frequency, and k is wave number,μmIt is fluid viscosity coefficient, κmIt is skeleton permeability, m=1,
2。
9. computer readable storage medium according to claim 8, wherein the average pressure are as follows:
Wherein,For average pressure,R is the radius of hole runner.
10. computer readable storage medium according to claim 9, wherein the equation of higher order of the frequency and wave number are as follows:
Wherein, Wherein,Gm,νmIt is modulus of shearing and Poisson's ratio, m=1,2.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710612448.2A CN109298443B (en) | 2017-07-25 | 2017-07-25 | Method for evaluating elastic wave of double-fluid jet flow model and computer readable storage medium |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710612448.2A CN109298443B (en) | 2017-07-25 | 2017-07-25 | Method for evaluating elastic wave of double-fluid jet flow model and computer readable storage medium |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109298443A true CN109298443A (en) | 2019-02-01 |
| CN109298443B CN109298443B (en) | 2020-12-11 |
Family
ID=65167315
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710612448.2A Active CN109298443B (en) | 2017-07-25 | 2017-07-25 | Method for evaluating elastic wave of double-fluid jet flow model and computer readable storage medium |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109298443B (en) |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111562613A (en) * | 2020-04-15 | 2020-08-21 | 东华理工大学 | Method for analyzing seismic wave frequency-dependent reflection coefficient of thin reservoir or mutual reservoir model |
| CN111597763A (en) * | 2020-04-09 | 2020-08-28 | 东华理工大学 | Cross multi-scale flow pore medium full-band elastic wave frequency dispersion attenuation analysis method |
| CN111812709A (en) * | 2020-07-08 | 2020-10-23 | 中国石油大学(北京) | Method, device and equipment for establishing multi-scale wave induced flow model |
| CN113283195A (en) * | 2021-06-04 | 2021-08-20 | 集美大学 | Shower outflow parameter prediction method, terminal device and storage medium |
| CN113391345A (en) * | 2021-06-02 | 2021-09-14 | 清华大学 | Method and device for predicting wave propagation characteristics, electronic device, and storage medium |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2005103766A2 (en) * | 2004-04-23 | 2005-11-03 | Schlumberger Canada Limited | Method and system for monitoring of fluid-filled domains in a medium based on interface waves propagating along their surfaces |
| CN103412336A (en) * | 2013-07-22 | 2013-11-27 | 中国石油天然气股份有限公司 | Method for predicting velocity of longitudinal wave of rock system in heterogeneous reservoir |
| CN103984027A (en) * | 2014-03-28 | 2014-08-13 | 清华大学 | Rock longitudinal wave speed prediction method based on ellipsoid double porosity model |
| CN104360383A (en) * | 2014-11-12 | 2015-02-18 | 中国石油大学(华东) | Method and system for predicting seismic wave attenuation |
| CN104570084A (en) * | 2013-10-29 | 2015-04-29 | 中国石油化工股份有限公司 | Cross-scale seismic rock physical attenuation model and method for predicating attenuation and dispersion |
-
2017
- 2017-07-25 CN CN201710612448.2A patent/CN109298443B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2005103766A2 (en) * | 2004-04-23 | 2005-11-03 | Schlumberger Canada Limited | Method and system for monitoring of fluid-filled domains in a medium based on interface waves propagating along their surfaces |
| CN103412336A (en) * | 2013-07-22 | 2013-11-27 | 中国石油天然气股份有限公司 | Method for predicting velocity of longitudinal wave of rock system in heterogeneous reservoir |
| CN104570084A (en) * | 2013-10-29 | 2015-04-29 | 中国石油化工股份有限公司 | Cross-scale seismic rock physical attenuation model and method for predicating attenuation and dispersion |
| CN103984027A (en) * | 2014-03-28 | 2014-08-13 | 清华大学 | Rock longitudinal wave speed prediction method based on ellipsoid double porosity model |
| CN104360383A (en) * | 2014-11-12 | 2015-02-18 | 中国石油大学(华东) | Method and system for predicting seismic wave attenuation |
Non-Patent Citations (1)
| Title |
|---|
| 孙卫涛 等: "孔隙介质弹性波频散-衰减理论模型", 《地球物理学进展》 * |
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111597763A (en) * | 2020-04-09 | 2020-08-28 | 东华理工大学 | Cross multi-scale flow pore medium full-band elastic wave frequency dispersion attenuation analysis method |
| CN111562613A (en) * | 2020-04-15 | 2020-08-21 | 东华理工大学 | Method for analyzing seismic wave frequency-dependent reflection coefficient of thin reservoir or mutual reservoir model |
| CN111562613B (en) * | 2020-04-15 | 2023-03-31 | 东华理工大学 | Method for analyzing seismic wave frequency-dependent reflection coefficient of thin reservoir or mutual reservoir model |
| CN111812709A (en) * | 2020-07-08 | 2020-10-23 | 中国石油大学(北京) | Method, device and equipment for establishing multi-scale wave induced flow model |
| CN111812709B (en) * | 2020-07-08 | 2022-03-04 | 中国石油大学(北京) | Method, device and equipment for establishing multi-scale wave induced flow model |
| CN113391345A (en) * | 2021-06-02 | 2021-09-14 | 清华大学 | Method and device for predicting wave propagation characteristics, electronic device, and storage medium |
| CN113283195A (en) * | 2021-06-04 | 2021-08-20 | 集美大学 | Shower outflow parameter prediction method, terminal device and storage medium |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109298443B (en) | 2020-12-11 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109298443A (en) | Two-fluid injection stream model elastic wave evaluation method and computer readable storage medium | |
| Ba et al. | Rock anelasticity due to patchy saturation and fabric heterogeneity: A double double‐porosity model of wave propagation | |
| Denneman et al. | Reflection and transmission of waves at a fluid/porous-medium interface | |
| Ba et al. | Biot‐Rayleigh theory of wave propagation in double‐porosity media | |
| Diallo et al. | Acoustic wave propagation in saturated porous media: reformulation of the Biot/Squirt flow theory | |
| Fabricius et al. | Elastic moduli of dry and water-saturated carbonates—Effect of depositional texture, porosity, and permeability | |
| Huang et al. | Brittleness index and seismic rock physics model for anisotropic tight-oil sandstone reservoirs | |
| Vogelaar et al. | Exact expression for the effective acoustics of patchy-saturated rocks | |
| Albers | Analysis of the propagation of sound waves in partially saturated soils by means of a macroscopic linear poroelastic model | |
| ZHANG et al. | Seismic wave propagation equations of conglomerate reservoirs: A triple-porosity structure model | |
| Qi et al. | Influence of interface condition on reflection of elastic waves in fluid-saturated porous media | |
| Brajanovski et al. | A model for strong attenuation and dispersion of seismic P-waves in a partially saturated fractured reservoir | |
| Mesgouez et al. | Transient solution for multilayered poroviscoelastic media obtained by an exact stiffness matrix formulation | |
| CN109505590A (en) | The determination method and computer readable storage medium of shale gas reservoir pressure | |
| DU et al. | Wavefield propagation characteristics in the fracture‐induced anisotropic double‐porosity medium | |
| Zhang et al. | Propagation and attenuation of P-waves in patchy saturated porous media | |
| Goyal et al. | Reflection and transmission of inhomogeneous waves at the plane interface between two dissimilar swelling porous half-spaces | |
| Dahl et al. | Local and global fluid effects on sonic wave modes | |
| Wang et al. | Effect of local fluid flow on the reflection and transmission of elastic waves at an interface between an elastic solid and a double-porosity medium | |
| Kumari et al. | Reflection of inhomogeneous waves at the surface of a dissipative poroelastic media | |
| CN107764697A (en) | Gas potential detection method based on the progressive equation non-linear inversion of pore media | |
| Wang et al. | Effects of the outermost boundary on acoustic waves in an artificial cased borehole | |
| CN109407148A (en) | One kind being based on the formation information prediction technique of " V " type AVO feature | |
| Kudarova | Effective models for seismic wave propagation in porous media | |
| Li et al. | Wettability-dependent wave velocities and attenuation in granular porous media |
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 |