CN109100801B - A kind of high-pressure medium seimic wave propagation analogy method - Google Patents
A kind of high-pressure medium seimic wave propagation analogy method Download PDFInfo
- Publication number
- CN109100801B CN109100801B CN201810721895.6A CN201810721895A CN109100801B CN 109100801 B CN109100801 B CN 109100801B CN 201810721895 A CN201810721895 A CN 201810721895A CN 109100801 B CN109100801 B CN 109100801B
- Authority
- CN
- China
- Prior art keywords
- rock
- modulus
- wave
- pressure
- value
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 34
- 239000011435 rock Substances 0.000 claims abstract description 92
- 239000011148 porous material Substances 0.000 claims abstract description 17
- 238000004088 simulation Methods 0.000 claims abstract description 17
- 238000005311 autocorrelation function Methods 0.000 claims abstract description 16
- 239000002245 particle Substances 0.000 claims abstract description 6
- 238000000605 extraction Methods 0.000 claims abstract description 4
- 238000002474 experimental method Methods 0.000 claims description 22
- 238000010008 shearing Methods 0.000 claims description 13
- 230000008859 change Effects 0.000 claims description 10
- 229910052500 inorganic mineral Inorganic materials 0.000 claims description 8
- 239000011707 mineral Substances 0.000 claims description 8
- 238000002591 computed tomography Methods 0.000 claims description 7
- 238000012360 testing method Methods 0.000 claims description 5
- 230000008569 process Effects 0.000 claims description 4
- 230000007613 environmental effect Effects 0.000 claims description 3
- 230000003068 static effect Effects 0.000 claims description 3
- 241000208340 Araliaceae Species 0.000 claims 1
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 claims 1
- 235000003140 Panax quinquefolius Nutrition 0.000 claims 1
- 235000008434 ginseng Nutrition 0.000 claims 1
- 235000013399 edible fruits Nutrition 0.000 description 5
- 238000010586 diagram Methods 0.000 description 3
- 230000001788 irregular Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 239000007787 solid Substances 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 1
- 238000012512 characterization method Methods 0.000 description 1
- 230000000052 comparative effect Effects 0.000 description 1
- 239000013078 crystal Substances 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000005284 excitation Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000002347 injection Methods 0.000 description 1
- 239000007924 injection Substances 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 230000000704 physical effect Effects 0.000 description 1
- 238000005381 potential energy Methods 0.000 description 1
- 230000011514 reflex Effects 0.000 description 1
- 239000000243 solution Substances 0.000 description 1
- 239000004575 stone Substances 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 230000017260 vegetative to reproductive phase transition of meristem Effects 0.000 description 1
- 210000001260 vocal cord Anatomy 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/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
- G01V2210/6242—Elastic parameters, e.g. Young, Lamé or Poisson
-
- 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)
- Geophysics And Detection Of Objects (AREA)
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
Abstract
The present invention relates to a kind of high-pressure medium seimic wave propagation analogy methods, comprising the following steps: 1) establishes the novel acoustic elasticity formula combined with parameter of pore structure；2) equivalent elastic modulus and elastic wave velocity under different pressures state are calculated using novel acoustic elasticity formula；3) structured data between extraction rock interior particle: auto-correlation length；4) it is directed to each pressure state, the equivalent elastic modulus that step 2) is obtained respectively is as the background value of rock medium elasticity modulus, and the auto-correlation length extracted according to step 3), in conjunction with auto-correlation function, elasticity modulus disturbance is added in background value, establishes the rock numerical model under each pressure state；5) it carries out finite element Simulation of Seismic Wave respectively in each rock numerical model, and is directed to deep reservoir characteristic, extract reservoir scatter attenuation Qc ‑1Value, wherein Qc ‑1Indicate Coda Attenuation coefficient.
Description
Technical field
The present invention relates to a kind of high-pressure medium seimic wave propagation analogy methods, belong to seismic exploration technique field.
Background technique
With seismic exploration technique research gradually deeply, seismic prospecting is gradually from superficial part reservoir exploration to Deep Oil-gas
Hiding and complex oil and gas reservoir exploration transformation, deep and complicated structure reservoir by ambient pressure, stress influence compared with
Greatly, deep reservoir dynamic mechanics parameter and shallow state are entirely different, the seismic prospecting for deep reservoir, and reservoir properties are multiple
Miscellaneous, reflection signal is weak, therefore, it is necessary to first establish reservoir properties model under true complex pressure environment, and according to numerical simulation,
Select the information that can identify reservoir characteristic.Deep reservoir rock is influenced seriously by pressure and complex stress.Conventional earthquake
Wave theory does not consider as the influence caused by physical properties of rock of environmental pressure and tectonic stress, so that conventional survey seismology
Method can not obtain true deep reservoir structure, can not also obtain the seismic profile of true reflection deep underground structure.
Acoustic elasticity method is a kind of method of rock elastic property under calculating complex pressure environment.This method passes through research
Change of the pressure to Rock Matrix elastic potential energy, to calculate the change of rock elasticity constant.Digital core technology, by rock
Stone carries out CT scan, can show rock interior structure.By digital core technology, rock interior structure can be studied to rock bullet
Property property influence, also, combined by digital core technology with numerical simulation, rock interior structure can be reacted to seismic wave
The influence of field scattering signatures.Under conditions of deep reflex signal is not strong, seismic wave field scattered information is studied, is with studying deep
The important method of flowering structure.
But the variation of the main study of rocks crystal inside lattice potential of Vocal cord injection at present, rock interior particle phase
Position, pore structure variation do not consider the influence that rock elastic property generates.Further, since deep reservoir reflects signal
It is weak, it is therefore desirable to study deep reservoir structure using wave field scattering signatures.
Summary of the invention
In view of the above-mentioned problems, the object of the present invention is to provide a kind of high-pressure medium seimic wave propagation analogy method, this method
It fully takes into account rock interior particle relative position, the influence that pore structure variation generates rock elastic property, can give
Rock elasticity constant, rock interior Scattering of Elastic Wave decay characteristics are calculated under conditions of level pressure force field parameter, are surveyed for deep reservoir
Spy lays the foundation.
To achieve the above object, the invention adopts the following technical scheme: a kind of high-pressure medium seimic wave propagation analogy method,
The following steps are included: 1) establish the novel acoustic elasticity formula combined with parameter of pore structure；2) novel acoustic elasticity formula is utilized
Calculate the equivalent elastic modulus and elastic wave velocity under different pressures state；3) structured data between extraction rock interior particle: from
Correlation length；4) it is directed to each pressure state, the equivalent elastic modulus for respectively obtaining step 2) is as rock medium elasticity modulus
Background value, and the auto-correlation length extracted according to step 3) is added elasticity modulus in background value and is disturbed in conjunction with auto-correlation function
It is dynamic, establish the rock numerical model under each pressure state；5) finite element seismic wave field is carried out respectively in each rock numerical model
Simulation, and it is directed to deep reservoir characteristic, extract reservoir scatter attenuation Qc -1Value, wherein Qc -1Indicate Coda Attenuation coefficient.
2, a kind of high-pressure medium seimic wave propagation analogy method as described in claim 1, it is characterised in that: complete to walk
It is rapid 4) and implementation steps 5) before, first choose rock numerical model under a part of pressure state and carry out finite element seismic wave field
Simulation, and extract reservoir scatter attenuation Qc -1Value, meanwhile, choose the Q that at least six petrophysics experiment obtainsc -1Data point, by wave
The Q that field simulation obtainsc -1The Q that value is obtained with petrophysics experiment respectivelyc -1It compares, if relative error illustrates less than 5%
Model parameter is correct.
The Q that at least six petrophysics experiment obtainsc -1Data point is different the experimental data under pressure state.
Novel acoustic elasticity formula in the step 1) are as follows:
In formula, ρ is the density of rock, VpAnd VsIt is equivalent rock longitudinal and shear wave velocity of wave.LdrsAnd μdrsIt is dry rock in soft hole
Longitudinal and shear wave modulus when closure, KdrsIt is bulk modulus of the dry rock in the soft closing of pores.WithIt is rock in soft hole
Three rank elastic constants of gap closure, P is the effective pressure of rock local environment.φc0It is the porosity of the soft hole of rock interior,
θc、θcL、θcμ、And θcKIt is all amount relevant with soft hole geometric shape.
Each determination method for parameter is as follows in the novel acoustic elasticity formula:
1. measuring rock density ρ, high pressure petrophysics experiment is carried out, pressure limit is 5MPa~60MPa；
2. choosing the data point in 50MPa~60MPa, velocity of wave is measured, calculates longitudinal and shear wave modulusAnd in pressure
When for 60MPa, static modulus test experiments are done, measure bulk modulus, bulk modulus can be used as K at this timedrs, environmental pressure at this time
Height, the soft closing of pores, soft hole influence are not present, and linear change is presented with effective pressure P in longitudinal and shear wave modulus, due to KdrsIt is normal
Number, longitudinal and shear wave modulus is with P/KdrsAlso it changes linearly, at this time the soft closing of pores, soft hole influence is not present, φc0、θc、θcL、
θcμ、And θcKIt is zero, chooses two data points of P=50MPa and P=60MPa and form equation group:
It solves above-mentioned equation group and obtains Ldrs、μdrs、Four amounts；
3. being fitted with pressure history according to the whole speed of experiment, seeking parameter phic0、θc、θcL、θcμ、And θcK, the related coefficient of the theoretical curve and actual curve that are obtained by novel hole acoustic elasticity formula reaches
When 98%, the value that the above parametric inversion obtains can be used as the soft hole geometric parameter of this rocks, for describing soft hole to rock
Equivalent elastic constant and velocity of wave influence.
In the step 3), extracting from correlation length is realized by digital cores technology.
The process for establishing rock numerical model is as follows:
1. determining the auto-correlation function of control disturbance distribution, scattering in view of middle sight disturbance on seismic exploration influences maximum, choosing
Exponential type auto-correlation function Φ is selected, prominent medium measure feature:
In formula, a is the auto-correlation length of the secondary mineral distribution of rock obtained according to CT scan, and value is to work as auto-correlation
When functional value falls to 1/e (e ≈ 2.71828), the value of r；R is data point position and core centre of slice in numerical model
Relative distance；The covariance C (x, z) and auto-correlation function relationship of disturbance are as follows:
Φ (x, z)=C (x, z)/σ2,
In formula, σ2It is the variance of disturbance, for moderate sorted sandstone, disturbs variance and choose 35%, pass through variance and association
Variance, it may be determined that disturbance quantity ε；
2. the equivalent elastic constant determined using novel hole acoustic elasticity formula is as Model Background value, model is by following at this time
Formula calculates；
Wherein:
In formula, K is the bulk modulus at any point (x, z) in space in model；μ is that in space any point in model
The modulus of shearing of (x, z)；The abscissa at any point in x representation space；The ordinate at any point in z representation space；μ0For
Modulus of shearing background value K0For bulk modulus background value, μ0And K0It is calculated and is obtained by novel hole acoustic elasticity formula；ερ、εK、εμ
It is disturbance quantity, specifically, ερThe density perturbation generated for rock heterogeneity；εKThe bulk modulus generated for rock heterogeneity
Disturbance；εμThe modulus of shearing disturbance generated for rock heterogeneity.
The invention adopts the above technical scheme, which has the following advantages: 1, the present invention is by constructing novel hole sound
Elastic formula can describe under pressure environment, the variation of the equivalent elastic modulus of pore media, compared to acoustic elasticity formula before,
The elastic modulus change that can preferably predict deep reservoir has very important significance for oil-gas exploration and development tool.2, this hair
It is bright to obtain the disturbance distribution of rock interior equivalent elastic modulus by extracting rock forming mineral distributed architecture using digital core technology
Situation, and in the form of disturbance quantity, it is superimposed upon in the variation of equivalent elastic modulus caused by acoustoelastic effect, in this way, can
To establish more true deep reservoir numerical model.3, the present invention uses finite element numerical simulation, can adapt to irregular shape
Shape boundary, meets free boundary condition automatically.4, the present invention can be obtained by establishing more true deep reservoir numerical model
More true deep reservoir earthquake record, to use scatter attenuation Qc -1The physical quantitys such as value detection deep reservoir structure is established
Basis.5, due to scatter attenuation Qc -1Value and the in-fighting decaying during seimic wave propagation are mutually indepedent, not by in-fighting decaying shadow
It rings, therefore, for reflection signal by the serious seismic signal of in-fighting decaying, using scatter attenuation Qc -1Value can more preferable identifying purpose
Reservoir.
Detailed description of the invention
Fig. 1 is certain petrophysical data measured data (circle), and seeks LdrsμdrsWithUsed in four amounts
Data matched curve (dotted line), wherein figure (a) is longitudinal wave modulus with pressure history, and scheming (b) is that modulus of shearing becomes with pressure
Change curve；
Fig. 2 is that certain petrophysical data measured data (circle) linear fit result (dotted line) and hole acoustic elasticity formula are pre-
It surveys result (solid line), figure (a) is longitudinal wave modulus with pressure history, and scheming (b) is modulus of shearing with pressure history；
Fig. 3 is certain rock sample slice and its corresponding CT scan number core picture；Wherein, figure (a) is rock-like
Product slice, figure (b) is CT scan number core picture,；
Fig. 4 is the auto-correlation function situation of change and auto-correlation length of rock sample slice；
Fig. 5 is rock sample numerical model, and white five-pointed star is excitaton source, and white triangles are receiving point；
Fig. 6 is finite element node valuation method schematic diagram；
Fig. 7 is triangular element assignment method schematic diagram；
Fig. 8 is the comparison of numerical simulation and experiment test result: where figure (a) is the longitudinal wave comparison knot under 60MPa pressure
Fruit, figure (b) is the longitudinal wave comparing result under 50MPa pressure, and figure (c) is the longitudinal wave comparing result under 40MPa pressure, figure
It (d) is the longitudinal wave comparing result under 30MPa pressure, figure (e) is the longitudinal wave comparing result under 20MPa pressure, and figure (f) is
Longitudinal wave comparing result under 15MPa pressure, figure (g) is the longitudinal wave comparing result under 10MPa pressure, and figure (h) is in 5MPa
Longitudinal wave comparing result under pressure；Scheming (i) is the shear wave comparing result under 60MPa pressure, and figure (j) is under 50MPa pressure
Shear wave comparing result, figure (k) be the shear wave comparing result under 40MPa pressure, figure (l) be the shear wave under 30MPa pressure
Comparing result, figure (m) is the shear wave comparing result under 20MPa pressure, and figure (n) is the shear wave comparison knot under 15MPa pressure
Fruit, figure (o) is the shear wave comparing result under 10MPa pressure, and figure (p) is the shear wave comparing result under 5MPa pressure；
Fig. 9 is longitudinal wave coda wave amplitude linearity section, wherein figure (a) is that effective pressure 30MPa-60MPa linearity range chooses knot
Fruit, figure (b) are that effective pressure 05MPa-20MPa linearity range chooses result；
Figure 10 is shear wave coda wave amplitude linearity section, wherein figure (a) is that effective pressure 30MPa-60MPa linearity range chooses knot
Fruit, figure (b) are that effective pressure 05MPa-20MPa linearity range chooses result；
Figure 11 is that numerical simulation and experiment tests coda wave scattering QcIt is worth Comparative result: where figure (a) is the scattering of longitudinal wave coda wave
QcValue comparison；Scheming (b) is that shear wave coda wave scatters QcValue comparison.
Specific embodiment
The present invention is described in detail below with reference to the accompanying drawings and embodiments.
The invention proposes a kind of high-pressure medium seimic wave propagation analogy methods, comprising the following steps:
1) the novel acoustic elasticity formula combined with parameter of pore structure is established, as follows:
In formula, ρ is the density of rock, VpAnd VsIt is equivalent rock longitudinal and shear wave velocity of wave.LdrsAnd μdrsIt is dry rock in soft hole
Longitudinal and shear wave modulus when closure, KdrsIt is bulk modulus of the dry rock in the soft closing of pores.WithIt is rock in soft hole
Three rank elastic constants of gap closure, P is the effective pressure of rock local environment.φc0It is the porosity of the soft hole of rock interior,
θc、θcL、θcμ、And θcKIt is all amount relevant with soft hole geometric shape, θcSoft hole geometric shape is represented to body
The influence of product variation, θcLRepresent influence of the soft hole geometric shape to longitudinal wave modulus, θcμSoft hole geometric shape is represented to shearing
The influence of modulus,Represent influence and θ of the soft hole geometric shape to three rank elastic constantscKRepresent soft hole geometry
Influence of the form to bulk modulus is usually obtained by the method for being fitted experimental data to make these parameters convenient for calculating.
The determination process of each parameter is as follows in above-mentioned novel acoustic elasticity formula:
1. measuring rock density ρ, high pressure petrophysics experiment is carried out, pressure limit is 5MPa~60MPa；
2. choosing the data point in 50MPa~60MPa, velocity of wave is measured, calculates longitudinal and shear wave modulusAnd it is pressing
When power is 60MPa, static modulus test experiments are done, measure bulk modulus, bulk modulus can be used as K at this timedrs, environment pressure at this time
Power is high, the soft closing of pores, and soft hole influence is not present, and linear change (such as Fig. 1 dotted line is presented with effective pressure P in longitudinal and shear wave modulus
It is shown), due to KdrsIt is constant, longitudinal and shear wave modulus is with P/KdrsAlso it changes linearly, at this time the soft closing of pores, soft hole influences not
In the presence of φc0、θc、θcL、θcμ、And θcKIt is zero, chooses the two data point composition sides P=50MPa and P=60MPa
Journey group:
It solves above-mentioned equation group and obtains Ldrs、μdrs、Four amounts；
3. it is fitted according to whole (effective pressure is changed to 60MPa from the 5MPa) speed of experiment with pressure history,
Seek parameter phic0、θc、θcL、θcμ、And θcK, the theoretical curve and reality that are obtained by novel hole acoustic elasticity formula
When the related coefficient of curve reaches 98% (as shown in Figure 2), the value that the above parametric inversion obtains can be used as the soft hole of this rocks
Geometric parameter, for describing influence of the soft hole to the equivalent elastic constant and velocity of wave of rock.
2) equivalent elastic modulus and elastic wave velocity (such as figure under different pressures state are calculated using novel acoustic elasticity formula
Shown in 2 solid black lines).
3) structured data between extraction rock interior particle: auto-correlation length.
It can specifically be realized by digital cores technology: core being sliced and carries out CT scan, determines secondary mineral space
Distribution calculates the auto-correlation length for controlling secondary mineral distribution；Fig. 3 (a) is core slice, is identified by RGB, core can be obtained
The spatial distribution (shown in such as Fig. 3 (b)) of internal mineral, by the spatial distribution of core internal mineral, the auto-correlation that can must be sliced
Function, curve is the dotted line in Fig. 4, and when auto-correlation function decays to 1/e (e=2.71828), the value of r is that auto-correlation is long
Degree.
4) it is directed to each pressure state, the equivalent elastic modulus for respectively obtaining step 2) is as rock medium elasticity modulus
Background value, and the auto-correlation length extracted according to step 3) add elasticity modulus in background value and disturb in conjunction with auto-correlation function
It is dynamic, establish the rock numerical model under each pressure state.
The process for establishing rock numerical model is as follows:
1. determining the auto-correlation function of control disturbance distribution, scattering in view of middle sight disturbance on seismic exploration influences maximum, choosing
Exponential type auto-correlation function Φ is selected, prominent medium measure feature:
In formula, a is the auto-correlation length of the secondary mineral distribution of rock obtained according to CT scan, and value is to work as auto-correlation
When functional value falls to 1/e (e ≈ 2.71828), the value of r；R is data point position and core centre of slice in numerical model
Relative distance；The covariance C (x, z) and auto-correlation function relationship of disturbance are as follows:
Φ (x, z)=C (x, z)/σ2,
In formula, σ2It is the variance of disturbance, for moderate sorted sandstone, disturbs variance and choose 35%, pass through variance and association
Variance, it may be determined that disturbance quantity ε；
2. the equivalent elastic constant determined using novel hole acoustic elasticity formula is as Model Background value, model is by following at this time
Formula calculates；
Wherein:
In formula, K is the bulk modulus at any point (x, z) in space in model；μ is that in space any point in model
The modulus of shearing of (x, z)；The abscissa at any point in x representation space；The ordinate at any point in z representation space；μ0For
Modulus of shearing background value K0For bulk modulus background value, μ0And K0It is calculated and is obtained by novel hole acoustic elasticity formula；ερ、εK、εμ
It is disturbance quantity, specifically, ερThe density perturbation generated for rock heterogeneity；εKThe bulk modulus generated for rock heterogeneity
Disturbance；εμThe modulus of shearing disturbance generated for rock heterogeneity；The numerical model of foundation is Fig. 6, size and rock physics
Experimental rock size is consistent.
5) it carries out finite element Simulation of Seismic Wave respectively in each rock numerical model, and is directed to deep reservoir characteristic, mention
Take reservoir scatter attenuation Qc -1Value, wherein Qc -1Indicate Coda Attenuation coefficient.
Preferably, finite element uses triangular element, and to meet irregular obstacle body, after mesh generation, node parameter is assigned
Value method is as shown in fig. 6, point O is triangular element node, a1、a2、a3、a4For the node unit node in core numerical simulation
Parameter, C0It is determined by following equation:
Wherein, C1、C2、C3、C4For the elastic parameter of core numerical model node.Fig. 7 is triangular element schematic diagram,
Elastic parameter is the average value of three node parameters.
Rock scatter attenuation passes through rock scatter attenuation Qc -1Value characterization, calculation method are as follows:
1. rock scatter attenuation Qc -1Value definition is；
Wherein, A is the coda wave amplitude in rock wave field record, and f is the centre frequency of incident wavelet, when t continues for coda wave
Between.H (f) is the geometrical factor of scattering, Cong Shizhong, it can be seen that logarithm is taken to the amplitude of wave recording, by seeking amplitude logarithm
The slope of the value part linear with the time, can be obtained scatter attenuation Qc -1Value.
2. the waveform and experimental waveform simulated are handled as follows: and the incidence wave waveform comparison of excitation, it receives
The waveform arrived, initial part is similar with incident waveform, is direct wave, casts out, and belongs to coda wave part after direct wave, takes to amplitude
Absolute value, and logarithm is sought, linear segment is found, slope is sought, according to slope, seeks scatter attenuation Qc -1Value.
Respectively wave-amplitude takes the amplitude after logarithm to change over time curve to Fig. 9 and Figure 10 in length and breadth, takes in solid box
Near line linearity curve does linear fit, seeks scattering Qc -1Value, can obtain Figure 11, and Figure 11 compared the scattering Q that experiment measurement obtainscValue with
And the scattering Q of simulationcValue seeks experiment QcValue and simulation QcRelated coefficient between value, related coefficient are greater than 98%, illustrate mould
Type is true.
In above-described embodiment, step 4) and implementation steps 5 are completed) before, it can first choose a part of pressure state
Under rock numerical model carry out finite element Simulation of Seismic Wave, and extract reservoir scatter attenuation Qc -1Value, meanwhile, choose at least 6
The Q that a petrophysics experiment obtainsc -1Data point, the Q that wave-field simulation is obtainedc -1The Q that value is obtained with petrophysics experiment respectivelyc -1It compares, if relative error less than 5%, illustrates that model parameter is correct.
Fig. 8 represents certain numerical simulation result and true measured data compares, and solid line is analog result, and dotted line is experiment knot
Fruit, direct wave first arrival time and through wave morphology are consistent.
The various embodiments described above are merely to illustrate the present invention, and wherein the implementation steps etc. of method may be changed,
All equivalents and improvement carried out based on the technical solution of the present invention, should not exclude in protection scope of the present invention
Except.
Claims (6)
1. a kind of high-pressure medium seimic wave propagation analogy method, comprising the following steps:
1) the novel acoustic elasticity formula combined with parameter of pore structure, the novel acoustic elasticity formula are established are as follows:
In formula, ρ is the density of rock, VpAnd VsIt is equivalent rock longitudinal and shear wave velocity of wave；LdrsAnd μdrsIt is dry rock in the soft closing of pores
When longitudinal and shear wave modulus, KdrsIt is bulk modulus of the dry rock in the soft closing of pores；WithIt is that rock is closed in soft hole
The three rank elastic constants closed, P is the effective pressure of rock local environment；φc0It is the porosity of the soft hole of rock interior, θc、θcL、
θcμ、And θcKIt is all amount relevant with soft hole geometric shape, θcSoft hole geometric shape is represented to volume change
Influence, θcLRepresent influence of the soft hole geometric shape to longitudinal wave modulus, θcμSoft hole geometric shape is represented to modulus of shearing
It influences,Represent influence of the soft hole geometric shape to three rank elastic constants, θcKRepresent soft hole geometric shape pair
The influence of bulk modulus；
2) equivalent elastic modulus and elastic wave velocity under different pressures state are calculated using novel acoustic elasticity formula；
3) structured data between extraction rock interior particle: auto-correlation length；
4) it is directed to each pressure state, the equivalent elastic modulus for respectively obtaining step 2) is as the background of rock medium elasticity modulus
Value, and the auto-correlation length extracted according to step 3) are added elasticity modulus disturbance in background value, are built in conjunction with auto-correlation function
Found the rock numerical model under each pressure state；
5) it carries out finite element Simulation of Seismic Wave respectively in each rock numerical model, and is directed to deep reservoir characteristic, extract storage
Layer scatter attenuation Qc -1Value, wherein Qc -1Indicate Coda Attenuation coefficient.
2. a kind of high-pressure medium seimic wave propagation analogy method as described in claim 1, it is characterised in that: complete step 4)
And implementation steps 5) before, the rock numerical model first chosen under a part of pressure state carries out finite element seismic wave field mould
It is quasi-, and extract reservoir scatter attenuation Qc -1Value, meanwhile, choose the Q that at least six petrophysics experiment obtainsc -1Data point, by wave field
Simulate obtained Qc -1The Q that value is obtained with petrophysics experiment respectivelyc -1It compares, if relative error illustrates mould less than 5%
Shape parameter is correct.
3. a kind of high-pressure medium seimic wave propagation analogy method as claimed in claim 2, it is characterised in that: described at least 6
The Q that a petrophysics experiment obtainsc -1Data point is different the experimental data under pressure state.
4. a kind of high-pressure medium seimic wave propagation analogy method as described in claim 1, it is characterised in that: the novel sound
Each determination method for parameter is as follows in elastic formula:
1. measuring rock density ρ, high pressure petrophysics experiment is carried out, pressure limit is 5MPa~60MPa；
2. choosing the data point in 50MPa~60MPa, velocity of wave is measured, calculates longitudinal and shear wave modulusAnd it is in pressure
When 60MPa, static modulus test experiments are done, measure bulk modulus, bulk modulus is as K at this timedrs, environmental pressure is high at this time, soft
The closing of pores, soft hole influence are not present, and linear change is presented with effective pressure P in longitudinal and shear wave modulus, due to KdrsIt is constant, indulges
Shear wave modulus is with P/KdrsAlso it changes linearly, at this time the soft closing of pores, soft hole influence is not present, φc0、θc、θcL、θcμ、And θcKIt is zero, chooses two data points of P=50MPa and P=60MPa and form equation group:
It solves above-mentioned equation group and obtains Ldrs、μdrs、Four amounts；
3. being fitted with pressure history according to the whole speed of experiment, seeking parameter phic0、θc、θcL、θcμ、With
θcK, when the related coefficient of the theoretical curve and actual curve that are obtained by novel hole acoustic elasticity formula reaches 98%, the above ginseng
The value that number inverting obtains, as the soft hole geometric parameter of this rocks, for describing soft hole to the equivalent elastic constant of rock
With the influence of velocity of wave.
5. a kind of high-pressure medium seimic wave propagation analogy method as described in claim 1, it is characterised in that: in the step 3)
In, extracting from correlation length is realized by digital cores technology.
6. a kind of high-pressure medium seimic wave propagation analogy method as described in claim 1, it is characterised in that: the step 4)
In, the process for establishing rock numerical model is as follows:
1. determining the auto-correlation function of control disturbance distribution, maximum on seismic exploration scattering influence in view of middle sight disturbance, selection refers to
Number type auto-correlation function Φ, prominent medium measure feature:
In formula, a is the auto-correlation length of the secondary mineral distribution of rock obtained according to CT scan, and value is to work as auto-correlation function
When value falls to 1/e, the value of r, wherein e ≈ 2.71828；R is data point position and core centre of slice in numerical model
Relative distance；The covariance C (x, z) and auto-correlation function relationship of disturbance are as follows:
Φ (x, z)=C (x, z)/σ2
In formula, σ2It is the variance of disturbance, for moderate sorted sandstone, disturbs variance and choose 35%, by variance and covariance,
Determine disturbance quantity ε；
2. the equivalent elastic constant determined using novel hole acoustic elasticity formula is as Model Background value, model is by following formula at this time
It calculates；
Wherein:
In formula, K is the bulk modulus at any point (x, z) in space in model；μ is that in space any point (x, z) in model
Modulus of shearing；The abscissa at any point in x representation space；The ordinate at any point in z representation space；μ0For shearing mould
Measure background value, K0For bulk modulus background value, μ0And K0It is calculated and is obtained by novel hole acoustic elasticity formula；ερ、εK、εμIt is disturbance
Amount, specifically, ερThe density perturbation generated for rock heterogeneity；εKThe bulk modulus disturbance generated for rock heterogeneity；εμ
The modulus of shearing disturbance generated for rock heterogeneity.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810721895.6A CN109100801B (en) | 2018-07-04 | 2018-07-04 | A kind of high-pressure medium seimic wave propagation analogy method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810721895.6A CN109100801B (en) | 2018-07-04 | 2018-07-04 | A kind of high-pressure medium seimic wave propagation analogy method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109100801A CN109100801A (en) | 2018-12-28 |
| CN109100801B true CN109100801B (en) | 2019-08-23 |
Family
ID=64845611
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810721895.6A Active CN109100801B (en) | 2018-07-04 | 2018-07-04 | A kind of high-pressure medium seimic wave propagation analogy method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109100801B (en) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110703322B (en) * | 2019-10-10 | 2020-11-13 | 清华大学 | Wave propagation processing method, device and equipment |
Family Cites Families (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| SU811168A1 (en) * | 1979-04-10 | 1981-03-07 | Научно-Производственное Объединение"Геофизика" | Method of determining elastic wave propagation speed in basin bottom ground |
| DE69932932D1 (en) * | 1999-10-22 | 2006-10-05 | Jason Geosystems B V | Method of determining the elastic parameters and rock composition of subterranean formations using seismic data |
| US6681185B1 (en) * | 2002-07-26 | 2004-01-20 | Eseis | Method of seismic signal processing |
| WO2012087176A2 (en) * | 2010-12-21 | 2012-06-28 | Schlumberger Holdings Limited | Method for computing wavefields |
| CN102288986B (en) * | 2011-05-16 | 2013-01-30 | 中国石油大学(北京) | Method for obtaining elastic modulus of carbonate rock reservoir under earthquake scale |
| CN103235338B (en) * | 2013-03-25 | 2015-06-10 | 中国石油大学(华东) | Rock fracture parameter inversion method |
| CN103462643B (en) * | 2013-09-29 | 2015-02-11 | 深圳市开立科技有限公司 | Shear wave speed measurement method, device and system |
| CN105092595B (en) * | 2015-08-31 | 2018-03-02 | 哈尔滨工业大学（威海） | Photoacoustic elastic imaging method and device applied to rail examination |
-
2018
- 2018-07-04 CN CN201810721895.6A patent/CN109100801B/en active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN109100801A (en) | 2018-12-28 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US11098565B2 (en) | Method for estimating permeability of fractured rock formations from induced slow fluid pressure waves | |
| Gholami et al. | Brittleness of gas shale reservoirs: A case study from the north Perth basin, Australia | |
| US9482771B2 (en) | Method of indicating the presence of gas hydrate and shallow gas in deepwater environment | |
| CN106597544B (en) | Method and device for predicting brittleness of compact oil and gas reservoir | |
| Hartzell et al. | Effects of 3D random correlated velocity perturbations on predicted ground motions | |
| CN103163553B (en) | Earthquake hydrocarbon detection method and device based on multiple pore medium model | |
| CN101231346A (en) | Method for estimating coal, rock mass physical mechanics parameter through seismic wave velocity | |
| CN105182424B (en) | Method and device for quantitatively predicting reservoir porosity based on patch saturation model | |
| Chehrazi et al. | Pore-facies as a tool for incorporation of small-scale dynamic information in integrated reservoir studies | |
| CN111399044A (en) | Reservoir permeability prediction method and device and storage medium | |
| Wang et al. | Analysis and estimation of an inclusion-based effective fluid modulus for tight gas-bearing sandstone reservoirs | |
| CN109100801B (en) | A kind of high-pressure medium seimic wave propagation analogy method | |
| CN110456412A (en) | A method of carbonate reservoir fluid saturation is identified based on post-stack seismic data | |
| Guo et al. | A fluid discrimination method based on Gassmann-Brie-Patchy Equation full waveform simulations and time-frequency analysis | |
| Cuixia et al. | Sensitivity parameters of tight sand gas: a case study of Lower Cretaceous Yingcheng Formation of Yingtai gas field in Songliao Basin, NE China | |
| US11346968B2 (en) | Estimation of reservoir flow properties from seismic data | |
| Li et al. | Study on the scales of heterogeneous geologic bodies in random media | |
| Jiang et al. | The constructing of pore structure factor in carbonate rocks and the inversion of reservoir parameters | |
| CN102809761A (en) | Control method for acquiring speed characteristics of coal rock | |
| Bianchi et al. | A Study of Progressive Failure in Porous Rocks Using Numerical and Experimental Modeling | |
| Yang et al. | Coherence estimation algorithm using Kendall’s concordance measurement on seismic data | |
| CN109901221A (en) | A kind of seismic data anisotropy modeling method based on NMO velocity parameter | |
| CN107764697A (en) | Gas potential detection method based on the progressive equation non-linear inversion of pore media | |
| Wang et al. | Fluids discrimination by ray-path elastic impedance inversion: A successful case from Sulige tight gas field | |
| Wang et al. | Direct prediction method of fracturing ability in shale formations based on pre-stack seismic inversion |
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 |