CN110275206A - A kind of crack-pore type rock physics Elastic forming board - Google Patents
A kind of crack-pore type rock physics Elastic forming board Download PDFInfo
- Publication number
- CN110275206A CN110275206A CN201910358736.9A CN201910358736A CN110275206A CN 110275206 A CN110275206 A CN 110275206A CN 201910358736 A CN201910358736 A CN 201910358736A CN 110275206 A CN110275206 A CN 110275206A
- Authority
- CN
- China
- Prior art keywords
- rock
- crack
- porosity
- wave
- reservoir
- 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
- 239000011435 rock Substances 0.000 title claims abstract description 101
- 239000011148 porous material Substances 0.000 title claims abstract description 30
- 239000012530 fluid Substances 0.000 claims abstract description 30
- 229910052500 inorganic mineral Inorganic materials 0.000 claims abstract description 7
- 239000011707 mineral Substances 0.000 claims abstract description 7
- 239000011159 matrix material Substances 0.000 claims description 24
- 239000004576 sand Substances 0.000 claims description 18
- 238000010008 shearing Methods 0.000 claims description 10
- 230000000875 corresponding effect Effects 0.000 claims description 9
- 238000006073 displacement reaction Methods 0.000 claims description 9
- 239000010433 feldspar Substances 0.000 claims description 9
- 239000010453 quartz Substances 0.000 claims description 9
- VYPSYNLAJGMNEJ-UHFFFAOYSA-N silicon dioxide Inorganic materials O=[Si]=O VYPSYNLAJGMNEJ-UHFFFAOYSA-N 0.000 claims description 9
- 238000012937 correction Methods 0.000 claims description 8
- 230000005540 biological transmission Effects 0.000 claims description 6
- 239000004615 ingredient Substances 0.000 claims description 6
- 238000001615 p wave Methods 0.000 claims description 6
- 230000002596 correlated effect Effects 0.000 claims description 5
- 238000011161 development Methods 0.000 claims description 4
- 238000002474 experimental method Methods 0.000 claims description 4
- 241000208340 Araliaceae Species 0.000 claims description 3
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 claims description 3
- 235000003140 Panax quinquefolius Nutrition 0.000 claims description 3
- 239000004927 clay Substances 0.000 claims description 3
- 239000000470 constituent Substances 0.000 claims description 3
- 230000008878 coupling Effects 0.000 claims description 3
- 238000010168 coupling process Methods 0.000 claims description 3
- 238000005859 coupling reaction Methods 0.000 claims description 3
- 230000000694 effects Effects 0.000 claims description 3
- 230000005284 excitation Effects 0.000 claims description 3
- 238000002637 fluid replacement therapy Methods 0.000 claims description 3
- 235000008434 ginseng Nutrition 0.000 claims description 3
- 230000003993 interaction Effects 0.000 claims description 3
- 239000000463 material Substances 0.000 claims description 3
- 230000035699 permeability Effects 0.000 claims description 3
- 238000012545 processing Methods 0.000 claims description 3
- 238000010200 validation analysis Methods 0.000 claims description 3
- 230000015572 biosynthetic process Effects 0.000 claims 1
- 238000005381 potential energy Methods 0.000 claims 1
- 229920006395 saturated elastomer Polymers 0.000 abstract description 3
- 230000000704 physical effect Effects 0.000 abstract description 2
- 238000000034 method Methods 0.000 description 5
- 239000004575 stone Substances 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
- 230000009977 dual effect Effects 0.000 description 2
- 239000000284 extract Substances 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 2
- 230000008901 benefit Effects 0.000 description 1
- 150000001875 compounds Chemical class 0.000 description 1
- 239000006185 dispersion Substances 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 210000002569 neuron Anatomy 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 239000003209 petroleum derivative Substances 0.000 description 1
- 230000004044 response Effects 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
-
- 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
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/40—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging
- G01V1/44—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging using generators and receivers in the same well
- G01V1/48—Processing data
- G01V1/50—Analysing data
Landscapes
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Remote Sensing (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)
Abstract
The invention discloses a kind of crack-pore type rock physics Elastic forming boards, belong to exploration rock geophysics field.For the tight sandstone reservoir for studying high saturated air, microfissure is developed, heterogeneity is strong, which is the elasticity modulus that mixed mineral is calculated using Voigt-Reuss-Hill model;Skeleton elasticity modulus containing crack, porous rocks is described using differential EFFECTIVE MEDIUM (DEM) model, by Biot-Rayleigh theoretical description hole partial fluid feature, and discloses the relationship of tight sandstone reservoir elastic parameter and physical property.
Description
Technical field
The invention belongs to explore rock geophysics field, and in particular to a kind of crack-pore type rock physics springform
Plate, by establishing the tight sandstone reservoir that rock physics Elastic forming board studies high saturated air, microfissure is developed, heterogeneity is strong.
Background technique
Tight sand generally has the geologic feature in low hole, the development of hypotonic and microcrack, and shows very strong non-equal
Even property.Compared with conventional sandstone reservoir, tight sandstone reservoir has apparent petrophysical property, permeation fluid mechanics properties
Difference.Smith etc. (2010) points out that crack in tight sand is the principal element for influencing seimic wave velocity, and microfissure is deposited
During the heterogeneity of the quantum of output that will affect tight sandstone reservoir oil gas, while rock interior will lead to elastic wave propagation
The loss (Carcione and Picotti, 2006) of energy, and the complexity of practical pore structure increases reservoir gas-bearing property prediction
Difficulty.Therefore, how more accurate, reasonable prediction and description Heterogeneous reservoir are always petroleum exploration technical research
Hot and difficult issue.
Petrophysical model is the bridge that earthquake information is converted into reservoir properties information.Biot and Gassmann is proposed double
Phase medium theory of wave propagation, the theoretical hypothesis fluid are uniformly distributed in porous rock, have been widely used in ground since self-forming
Seismic exploration and rock engineering field, but there are biggish deviations for the description of actual reservoir medium.Later White, Nur and
Mavko and Dvorkin considers the heterogeneity of rock interior, has carried out to Biot-Gassmann theory perfect.Berryman
Gassmann equation is generalized in the compound porous material containing two kinds of pore structures with Milton, but not by local flow
Dynamic and double-pore structure connects well.Pride etc. is based on volume approximation on the average method, derives isotropism double porosity media
In fluid neuron network equation.Achievement based on Pride, the double-porosity system description that Ba Jing etc. discusses saturation single fluid are real
The reasonability of border rock medium wave propagation phenomenon, and realized based on Biot-Rayleigh wave equation to the more of unsaturation rock
Scale issue modeling and engineer application.Xu and White (1995) is based on the experiment analysis results of (1986) water bearing sands such as Han,
Other models, Wyllie equation and Gassmann equation are comprehensively utilized, sand, the mud stone mixed model of proposition are widely used in reservoir
Shear wave speed calculates and Rock physical analysis, but the computational accuracy of the model is affected by input parameter, therefore in industry
The practical application effect on boundary is not fine.Avseth etc. is based on Biot-Gassmann theory and constructs rock physics template, and answers
For petroleum-gas prediction.Ruiz etc. (2008) considers the influence in crack based on the soft pore model that tight sand proposes, can be fine
Matching tight sand petrophysical model.Yan etc. (2011) is based on effective media theory and experimental analysis gives Chuan Zhongdi
Qu Zhong, low hole sandstone skeleton model, and point out that matrix modulus, pore shape and pore components can all influence skeleton modulus.
End of blade south etc. (2015) is based on differential EFFECTIVE MEDIUM THEORY and the theoretical building petrophysical model of patch shape saturation is applied to distinguish
Sand shale.The reservoir rock physical model that Wang great Xing (2016) is established based on the analysis of experimental data of Soviet Union's Sulige gas field tight sand
It is successfully applied to gas distribution prediction.Guo Mengqiu etc. (2018) is based on dual dual pore structure model analysis sandstone containing fluid tight
Longitudinal wave frequency dispersion and decay characteristics.Non- Xian etc. (2018) establishes more rulers of tight sand based on improved random patchy saturation
Rock physics template is spent to detect applied to gas-bearing layer.Yang Peijie (2018) is based on rock theoretical basis and proposes a kind of Simultaneous Retrieving sand
The method of mudstone porosity and water saturation, this method reduce the multi-solutions of reservoir prediction.
For conventional gas and oil, tight sand oil-gas reservoir must use unconventional thinking and technology, be caused by research
The heterogeneity of close sandstone reservoir is influenced caused by seismic response, and then from according to earthquake Data Inversion tight sand oil gas
Reservoir heterogeneity and its distribution.
Summary of the invention
The present invention is to overcome conventional petrophysical model only to consider single factors, and provide a kind of crack -- pore type rock
Elastic physical property template.For the strong tight sandstone reservoir of high saturated air, microfissure development, heterogeneity, by establishing rock object
It manages Elastic forming board and studies tight sandstone reservoir.
The present invention solves its technical problem and is achieved through the following technical solutions:
A kind of crack-pore type rock physics Elastic forming board, establishment step are as follows:
(1) quartz, the ingredient of feldspar and its content in Rock Matrix are obtained according to geologic information, obtains the bullet of Rock Matrix
Property parameter and density;
(2) double porosity media wave transmission controe equation is introduced;
(3) Voigt-Reuss-Hill model, DEM model and Biot-Rayleigh establishing equation crack-pore type are combined
Petrophysical model;
(4) elastic parameter impact analysis and the building of rock physics chart board are carried out;
(5) porosity and crack content prediction are carried out by example, and make inversion result figure, validation template it is effective
Property.
Further, the step (2) specifically:
The interregional local fluid flow interaction of different aperture is introduced into strain energy, kinetic energy, establishes corresponding gesture
Energy function, energetic function and dissipative function, and then derive double porosity media wave transmission controe equation are as follows:
U in formula, U(1),U(2)Respectively the average grain displacement of the dry skeleton of rock, the fluid phase 1 of fluid in main body framework,
The displacement of the fluid phase 2 of fluid, ε, ζ in microfissure(1),ζ(2)It is corresponding 3 displacements Divergence Field;Indicate earthquake wave excitation mistake
The local fluid deformation increment generated in journey, rock interior is developed due to the heterogeneity of pore structure two class different apertures,
φ10And φ20It is the local porosity of main body framework and crack skeleton, R12For the radius of microfissure;φ1,φ2It is two class holes
Absolute porosity;ρf,η,κ1For the density of fluid, viscosity and permeability, A, N, Q1、R1、Q2With R2For elastic parameter, ρ11、ρ12、
ρ13、ρ22With ρ33For density parameter, b1With b2For Dissipation Parameters.
Further, the step (3) specifically:
1) quartz, the ingredient of feldspar and its content in Rock Matrix are obtained according to geologic information, uses Voigt-Reuss-
Hill model seeks the elasticity modulus M of Rock MatrixVRH:
In formula: fi、MiIndicate volume fraction, the elasticity modulus of i-th kind of mineral constituent, MVExpression is sought using Voigt model
Rock Matrix elasticity modulus, MRIndicate the Rock Matrix elasticity modulus sought using Reuss model;
2) dry hole, crack are added in Rock Matrix using DEM model, obtain the elasticity ginseng of the dry skeleton of rock
Several and density, the differential equation group of the coupling of equivalent volume and modulus of shearing are as follows:
Wherein primary condition is K*(0)=K1, μ*(0)=μ1, K1, μ1The bulk modulus and shearing mould of=initial major phase material
Amount, K2, μ2The bulk modulus and modulus of shearing of=the inclusion being gradually added into, the content of y=inclusion, for fluid inclusion
With empty inclusion, y is equal to porosity, P*iAnd Q*iIt is for self-compatibility equivalent elastic modulus μSC *And KSC *Background media
In i-th kind of component the mineral form factor;
3) fluid replacement is carried out using Biot-Rayleigh equation, gas is added to the dry skeleton of rock, obtain full gas rock
The velocity of longitudinal wave and shear wave velocity of stone finally obtain description crack -- the heterogeneous diplopore petrophysical model of porosity reservoir.
Further, elastic Analysis of Parameter Effect in the step (4) specifically:
Based on established petrophysical model, sunykatuib analysis porosity and influence of the crack content to elastic parameter, root
According to the log data curve matching figure of gas-bearing reservoir, P-S wave velocity ratio and being positively correlated property of porosity are obtained, with velocity of longitudinal wave
Negatively correlated property.
Further, rock physics chart board is formed in the step (4) specifically:
Rock physics modeling is carried out to work area tight sand target zone, using quartz, feldspar and clay as matrix, is used
DEM model obtains the dry skeleton elasticity modulus of the rock containing porosity, microfissure rock, is obtained using Biot-Rayleigh equation
The velocity of longitudinal wave and shear wave velocity of different frequency range form crack -- pore type petrophysics chart board based on the above parameter, then
Chart board correction is carried out using measured data, that is, experimental data, log data.
Further, the step (5) specifically:
Reservoir porosity, the crack content of well survey line are crossed by quantitative interpretation, and prestack inversion is carried out to interval of interest first
The data volume of P-wave impedance and P-S wave velocity ratio is obtained, the value of p-wave impedance and P-S wave velocity ratio that then inverting obtains
It extracts, obtains wave impedance and P-S wave velocity ratio two dimensional cross-section;
On rock physics Elastic forming board after obtained P-S wave velocity ratio and p-wave impedance value to be projected to Data correction,
Within the scope of the reservoir parameter of template, the nearest template lattice point of range data point is judged, and its porosity, crack are contained into numerical quantity
As reservoir parameter corresponding to the data point;
In target zone, non-reservoir is done away from larger situation to the boundary difference of data point and template and is handled, except template
It is not exploration targets including the too low compact reservoir of porosity, can directly does non-reservoir processing;
Establish tight sand crack -- porous rocks physical template, using Rock experiment observation and well-log information to rock object
Reason template is verified.
The invention has the benefit that
The present invention
1. comparing other conventional rock physical models, which can reasonably be applied to tight sandstone reservoir porosity
And in the explanation of crack content.
2. the model can more accurate, reasonable prediction and description Heterogeneous reservoir.
3. the tight sand rock physics template that new parameter crack content constructs is based on earthquake data before superposition and carries out reservoir
Porosity and crack content prediction advance the development of non-homogeneous saturation reservoir porosity and crack content prediction.
Detailed description of the invention
Fig. 1 is rock physics modeling procedure schematic diagram of the present invention;
Fig. 2 is work area gas-bearing layer log data curve matching schematic diagram of the present invention;
Fig. 3 is the relationship of skeleton bulk modulus (a), modulus of shearing (b) and porosity and crack content;
Fig. 4 is the relationship of velocity of longitudinal wave (a), shear wave velocity (b) and porosity and crack content;
Fig. 5 is the relationship of P-S wave velocity ratio Yu porosity and crack content;
Fig. 6 is the relationship of Poisson's ratio Yu porosity and crack content;
Fig. 7 is log data correction;
Fig. 8 is experimental data correction;
Fig. 9 is P-S wave velocity ratio (a) and P-wave impedance (b) two-dimension earthquake section;
Figure 10 is crack content (a) and porosity (b) inversion result.
Specific embodiment
Below by specific embodiment, the invention will be further described, and it is not limit that following embodiment, which is descriptive,
Qualitatively, this does not limit the scope of protection of the present invention.
A kind of crack-pore type rock physics Elastic forming board, establishment step are as follows:
(1) quartz, the ingredient of feldspar and its content in Rock Matrix are obtained according to geologic information, obtains the bullet of Rock Matrix
Property parameter and density;
(2) double porosity media wave transmission controe equation is introduced:
The interregional local fluid flow interaction of different aperture is introduced into strain energy, kinetic energy, establishes corresponding gesture
Energy function, energetic function and dissipative function, and then derive double porosity media wave transmission controe equation are as follows:
U in formula, U(1),U(2)Respectively the average grain displacement of the dry skeleton of rock, the fluid phase 1 of fluid in main body framework,
The displacement of the fluid phase 2 of fluid, ε, ζ in microfissure(1),ζ(2)It is corresponding 3 displacements Divergence Field;Indicate earthquake wave excitation mistake
The local fluid deformation increment generated in journey, rock interior is developed due to the heterogeneity of pore structure two class different apertures,
It is hole and crack, φ respectively10And φ20It is the local porosity of main body framework and crack skeleton, R12For the radius of microfissure;
φ1,φ2It is the absolute porosity of two class holes;ρf,η,κ1For the density of fluid, viscosity and permeability, A, N, Q1、R1、Q2With R2
For elastic parameter, ρ11、ρ12、ρ13、ρ22With ρ33For density parameter, b1With b2For Dissipation Parameters.
(3) Voigt-Reuss-Hill model, DEM model and Biot-Rayleigh establishing equation crack-pore type are combined
Petrophysical model:
1) quartz, the ingredient of feldspar and its content in Rock Matrix are obtained according to geologic information, uses Voigt-Reuss-
Hill model seeks the elasticity modulus M of Rock MatrixVRH:
In formula: fi、MiIndicate volume fraction, the elasticity modulus of i-th kind of mineral constituent, MVExpression is sought using Voigt model
Rock Matrix elasticity modulus, MRIndicate the Rock Matrix elasticity modulus sought using Reuss model;
2) dry hole, crack are added in Rock Matrix using DEM model, obtain the elasticity ginseng of the dry skeleton of rock
Several and density, the differential equation group of the coupling of equivalent volume and modulus of shearing are as follows:
Wherein primary condition is K*(0)=K1, μ*(0)=μ1, K1, μ1The bulk modulus and shearing mould of=initial major phase material
Amount, K2, μ2The bulk modulus and modulus of shearing of=the inclusion being gradually added into, the content of y=inclusion, for fluid inclusion
With empty inclusion, y is equal to porosity, P*iAnd Q*iIt is for self-compatibility equivalent elastic modulus μSC *And KSC *Background media
In i-th kind of component the mineral form factor;
3) fluid replacement is carried out using Biot-Rayleigh equation, gas is added to the dry skeleton of rock, obtain full gas rock
The velocity of longitudinal wave and shear wave velocity of stone finally obtain description crack -- the heterogeneous diplopore petrophysical model of porosity reservoir.
(4) elastic parameter impact analysis and the building of rock physics chart board are carried out:
The resilient nature and porosity and crack content of rock are closely related.In general, porosity and crack content
Variation will have a direct impact on the variation of p-and s-wave velocity.Rock total porosity and crack content are chosen to analyze it to elastic parameter
Influence help to be better described with analysis tight sandstone reservoir Characteristics of Seismic Wave Propagation.
Based on established petrophysical model, sunykatuib analysis porosity and influence of the crack content to elastic parameter, root
According to the log data curve matching figure of gas-bearing reservoir, P-S wave velocity ratio and being positively correlated property of porosity are obtained, with velocity of longitudinal wave
Negatively correlated property.
Rock physics modeling is carried out to work area tight sand target zone, using quartz, feldspar and clay as matrix, is used
DEM model obtains the dry skeleton elasticity modulus of the rock containing porosity, microfissure rock, is obtained using Biot-Rayleigh equation
The velocity of longitudinal wave and shear wave velocity of different frequency range form crack -- pore type petrophysics chart board based on the above parameter, then
Chart board correction is carried out using measured data, that is, experimental data, log data.
Need to adjust each input parameter on the basis of theory before exporting final rock physics chart board, make rock physics chart board with
The regularity of distribution of data point is consistent.
(5) porosity and crack content prediction are carried out by example, and make inversion result figure, validation template it is effective
Property:
Reservoir porosity, the crack content of well survey line are crossed by quantitative interpretation, and prestack inversion is carried out to interval of interest first
The data volume of P-wave impedance and P-S wave velocity ratio is obtained, the value of p-wave impedance and P-S wave velocity ratio that then inverting obtains
It extracts, obtains wave impedance and P-S wave velocity ratio two dimensional cross-section;
On rock physics Elastic forming board after obtained P-S wave velocity ratio and p-wave impedance value to be projected to Data correction,
Within the scope of the reservoir parameter of template, the nearest template lattice point of range data point is judged, and its porosity, crack are contained into numerical quantity
As reservoir parameter corresponding to the data point;
In target zone, non-reservoir is done away from larger situation to the boundary difference of data point and template and is handled, except template
It is not exploration targets including the too low compact reservoir of porosity, can directly does non-reservoir processing;
Based on practical geologic information, tight sand crack is established -- porous rocks physical template.It is observed using Rock experiment
And well-log information verifies rock physics template, needs to adjust respectively on the basis of theory before exporting final rock physics chart board
Parameter is inputted, keeps rock physics chart board consistent with the regularity of distribution of data point, sees Fig. 8.The template can preferably be applied to cause
During close elements of sandstone porosity and crack content are explained.
The above is only a preferred embodiment of the present invention, it is noted that for the ordinary skill people of the art
For member, various improvements and modifications may be made without departing from the principle of the present invention, these improvements and modifications are also answered
It is considered as protection scope of the present invention.
Claims (6)
1. a kind of crack-pore type rock physics Elastic forming board, it is characterised in that: establishment step is as follows:
(1) quartz, the ingredient of feldspar and its content in Rock Matrix are obtained according to geologic information, obtains the elasticity ginseng of Rock Matrix
Several and density;
(2) double porosity media wave transmission controe equation is introduced;
(3) Voigt-Reuss-Hill model, DEM model and Biot-Rayleigh establishing equation crack-pore type rock are combined
Physical model;
(4) elastic parameter impact analysis and the building of rock physics chart board are carried out;
(5) porosity and crack content prediction are carried out by example, and makes inversion result figure, the validity of validation template.
2. a kind of crack-pore type rock physics Elastic forming board as described in claim 1, it is characterised in that: the step (2)
Specifically:
The interregional local fluid flow interaction of different aperture is introduced into strain energy, kinetic energy, establishes corresponding potential energy letter
Number, energetic function and dissipative function, and then derive double porosity media wave transmission controe equation are as follows:
U in formula, U(1),U(2)Respectively the average grain displacement of the dry skeleton of rock, in main body framework fluid fluid phase 1, fine fisssure
The displacement of the fluid phase 2 of fluid, ε, ζ in gap(1),ζ(2)It is corresponding 3 displacements Divergence Field;During expression earthquake wave excitation
The local fluid deformation increment of generation, rock interior have two class different apertures, φ due to the heterogeneity development of pore structure10With
φ20It is the local porosity of main body framework and crack skeleton, R12For the radius of microfissure;φ1,φ2It is the absolute of two class holes
Porosity;ρf,η,κ1For the density of fluid, viscosity and permeability, A, N, Q1、R1、Q2With R2For elastic parameter, ρ11、ρ12、ρ13、ρ22
With ρ33For density parameter, b1With b2For Dissipation Parameters.
3. a kind of crack-pore type rock physics Elastic forming board as described in claim 1, it is characterised in that: the step (3)
Specifically:
1) quartz, the ingredient of feldspar and its content in Rock Matrix are obtained according to geologic information, uses Voigt-Reuss-Hill
Model seeks the elasticity modulus M of Rock MatrixVRH:
In formula: fi、MiIndicate volume fraction, the elasticity modulus of i-th kind of mineral constituent, MVIndicate the rock sought using Voigt model
Ground mass matter elasticity modulus, MRIndicate the Rock Matrix elasticity modulus sought using Reuss model;
2) dry hole, crack are added in Rock Matrix using DEM model, obtain the dry skeleton of rock elastic parameter and
The differential equation group of the coupling of density, equivalent volume and modulus of shearing is as follows:
Wherein primary condition is K*(0)=K1, μ*(0)=μ1, K1, μ1The bulk modulus and modulus of shearing of=initial major phase material,
K2, μ2The bulk modulus and modulus of shearing of=the inclusion being gradually added into, the content of y=inclusion, for fluid inclusion and sky
Inclusion, y are equal to porosity, P*iAnd Q*iIt is for self-compatibility equivalent elastic modulus μSC *And KSC *Background media in i-th
The mineral form factor of kind component;
3) fluid replacement is carried out using Biot-Rayleigh equation, gas is added to the dry skeleton of rock, obtain full gas rock
Velocity of longitudinal wave and shear wave velocity finally obtain description crack -- the heterogeneous diplopore petrophysical model of porosity reservoir.
4. a kind of crack-pore type rock physics Elastic forming board as described in claim 1, it is characterised in that: the step (4)
Middle elasticity Analysis of Parameter Effect specifically:
Based on established petrophysical model, sunykatuib analysis porosity and influence of the crack content to elastic parameter, according to containing
The log data curve matching figure of gas reservoir, obtains P-S wave velocity ratio and being positively correlated property of porosity, with velocity of longitudinal wave in negative
Correlation.
5. a kind of crack-pore type rock physics Elastic forming board as claimed in claim 4, it is characterised in that: the step (4)
Middle formation rock physics chart board specifically:
Rock physics modeling is carried out to work area tight sand target zone, using quartz, feldspar and clay as matrix, using DEM mould
Type obtains the dry skeleton elasticity modulus of the rock containing porosity, microfissure rock, obtains different frequencies using Biot-Rayleigh equation
The velocity of longitudinal wave and shear wave velocity of section form crack -- pore type petrophysics chart board based on the above parameter, then using real
Measured data, that is, experimental data, log data carry out chart board correction.
6. a kind of crack-pore type rock physics Elastic forming board as described in claim 1, it is characterised in that: the step (5)
Specifically:
Reservoir porosity, the crack content of well survey line are crossed by quantitative interpretation, and prestack inversion is carried out to interval of interest first and is obtained
The value of the data volume of P-wave impedance and P-S wave velocity ratio, the p-wave impedance and P-S wave velocity ratio that then inverting obtains is extracted
Out, wave impedance and P-S wave velocity ratio two dimensional cross-section are obtained;
On rock physics Elastic forming board after obtained P-S wave velocity ratio and p-wave impedance value to be projected to Data correction, in mould
Within the scope of the reservoir parameter of plate, judge the nearest template lattice point of range data point, and using its porosity, crack containing numerical quantity as
Reservoir parameter corresponding to the data point;
In target zone, non-reservoir is done away from larger situation to the boundary difference of data point and template and is handled, except template includes
The too low compact reservoir of porosity is not exploration targets, can directly do non-reservoir processing;
Establish tight sand crack -- porous rocks physical template, using Rock experiment observation and well-log information to rock physics mould
Plate is verified.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910358736.9A CN110275206B (en) | 2019-08-12 | 2019-08-12 | Fracture-pore rock physical elastic template |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910358736.9A CN110275206B (en) | 2019-08-12 | 2019-08-12 | Fracture-pore rock physical elastic template |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110275206A true CN110275206A (en) | 2019-09-24 |
CN110275206B CN110275206B (en) | 2021-09-28 |
Family
ID=67959708
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910358736.9A Active CN110275206B (en) | 2019-08-12 | 2019-08-12 | Fracture-pore rock physical elastic template |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110275206B (en) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110703322A (en) * | 2019-10-10 | 2020-01-17 | 清华大学 | Wave propagation processing method, device and equipment |
CN112379416A (en) * | 2020-10-13 | 2021-02-19 | 北京恒泰兴科信息技术有限公司 | Method and device for predicting transverse wave through coal rock physical modeling and electronic equipment |
CN112505772A (en) * | 2020-12-10 | 2021-03-16 | 中国石油大学(华东) | Method for inverting rock pore distribution characteristics by utilizing pore and fracture medium elastic wave theory |
CN114674934A (en) * | 2022-02-18 | 2022-06-28 | 河海大学 | Method for establishing theoretical model of change of wave velocity of saturated heavy oil rock along with temperature |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102854531A (en) * | 2012-09-11 | 2013-01-02 | 中国石油天然气股份有限公司 | Multi-scale rock physical charting method and device for detecting reservoir hydrocarbon |
CN103424772A (en) * | 2012-05-24 | 2013-12-04 | 中国石油化工股份有限公司 | Reservoir shear wave velocity prediction method based on rock physics |
CN103576195A (en) * | 2013-10-28 | 2014-02-12 | 西北大学 | Method for forecasting fissured medium transverse wave velocity varying with pressure |
CN103984027A (en) * | 2014-03-28 | 2014-08-13 | 清华大学 | Rock longitudinal wave speed prediction method based on ellipsoid double porosity model |
CN105095631A (en) * | 2014-05-21 | 2015-11-25 | 中国石油化工股份有限公司 | Shale anisotropic rock physical modeling method |
US20160109593A1 (en) * | 2014-10-17 | 2016-04-21 | Vimal SAXENA | Methods and systems for generating percolated rock physics models for predicting permeability and petrophysical quantities |
CN106054248A (en) * | 2016-07-15 | 2016-10-26 | 河海大学 | Earthquake rock physical inversion method based on large area tight reservoir |
CN109655940A (en) * | 2017-10-12 | 2019-04-19 | 中国石油化工股份有限公司 | Shale anisotropic rock physical model modeling method |
-
2019
- 2019-08-12 CN CN201910358736.9A patent/CN110275206B/en active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103424772A (en) * | 2012-05-24 | 2013-12-04 | 中国石油化工股份有限公司 | Reservoir shear wave velocity prediction method based on rock physics |
CN102854531A (en) * | 2012-09-11 | 2013-01-02 | 中国石油天然气股份有限公司 | Multi-scale rock physical charting method and device for detecting reservoir hydrocarbon |
CN103576195A (en) * | 2013-10-28 | 2014-02-12 | 西北大学 | Method for forecasting fissured medium transverse wave velocity varying with pressure |
CN103984027A (en) * | 2014-03-28 | 2014-08-13 | 清华大学 | Rock longitudinal wave speed prediction method based on ellipsoid double porosity model |
CN105095631A (en) * | 2014-05-21 | 2015-11-25 | 中国石油化工股份有限公司 | Shale anisotropic rock physical modeling method |
US20160109593A1 (en) * | 2014-10-17 | 2016-04-21 | Vimal SAXENA | Methods and systems for generating percolated rock physics models for predicting permeability and petrophysical quantities |
CN106054248A (en) * | 2016-07-15 | 2016-10-26 | 河海大学 | Earthquake rock physical inversion method based on large area tight reservoir |
CN109655940A (en) * | 2017-10-12 | 2019-04-19 | 中国石油化工股份有限公司 | Shale anisotropic rock physical model modeling method |
Non-Patent Citations (3)
Title |
---|
LIANG ZUO,等: "Elastic Properties of Polycrystals in the Voigt-Reuss-Hill Approximation", 《J. APPL. CRYST.》 * |
印兴耀,等: "储层地震岩石物理建模研究现状与进展", 《石油物探》 * |
许孝凯,等: "基于岩石物理的致密气层识别方法研究", 《应用声学》 * |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110703322A (en) * | 2019-10-10 | 2020-01-17 | 清华大学 | Wave propagation processing method, device and equipment |
CN112379416A (en) * | 2020-10-13 | 2021-02-19 | 北京恒泰兴科信息技术有限公司 | Method and device for predicting transverse wave through coal rock physical modeling and electronic equipment |
CN112379416B (en) * | 2020-10-13 | 2024-02-06 | 北京恒泰兴科信息技术有限公司 | Method and device for predicting transverse waves through coal rock physical modeling and electronic equipment |
CN112505772A (en) * | 2020-12-10 | 2021-03-16 | 中国石油大学(华东) | Method for inverting rock pore distribution characteristics by utilizing pore and fracture medium elastic wave theory |
CN114674934A (en) * | 2022-02-18 | 2022-06-28 | 河海大学 | Method for establishing theoretical model of change of wave velocity of saturated heavy oil rock along with temperature |
Also Published As
Publication number | Publication date |
---|---|
CN110275206B (en) | 2021-09-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110275206A (en) | A kind of crack-pore type rock physics Elastic forming board | |
EP1820137B1 (en) | Integrated anisotropic rock physics model | |
CN103713320B (en) | Organic-matter-rich mud shale rock physical model establishing method | |
CN102426390B (en) | Method for determining reserve volume of nonhomogeneous sandstone reservoir | |
Vernik et al. | Ultrasonic velocity and anisotropy of hydrocarbon source rocks | |
CN102870087B (en) | The method and system of fluid limited bulk emulation | |
CN102445708B (en) | Rich argillaceous sand stone three-dimension equivalent speed prediction model | |
CN102508296B (en) | Method and device for analyzing dispersion and attenuation of unsaturated double-porosity medium earthquake waves | |
CN106054248A (en) | Earthquake rock physical inversion method based on large area tight reservoir | |
CN104007485B (en) | One sets up complicated hole mud shale petrophysical model method | |
CN105487122B (en) | For the method and system for the sensitiveness for determining the fluid identification of reservoir factor | |
CN110333549A (en) | Well point fracture aperture prediction technique based on well-log information Dispersion Characteristics | |
CN104570065B (en) | Method for quantitatively inverting porosity by using seismic wave impedance | |
CN105093351B (en) | The method for recognizing reservoir micro-fractures | |
CN115793048A (en) | Shale oil bed system compact reservoir transverse wave velocity prediction method | |
CN104330822A (en) | Method and device for determining remaining oil and gas distribution through coupling four-dimensional seismic inversion | |
Li et al. | Nonlinear simultaneous inversion of pore structure and physical parameters based on elastic impedance | |
CN110954942A (en) | Method and system for making chart with clay content identification function | |
CN105717535A (en) | Variable-parameter fracture model material and preparing method thereof | |
CN104459771A (en) | Reservoir gas-bearing semi-quantitative recognition method based on frequency division AVO inversion | |
Ren et al. | Anisotropic rock-physics modeling for the Haynesville Shale | |
Shukla et al. | Effect of capillary pressure on seismic velocities and attenuation | |
CN107764697A (en) | Gas potential detection method based on the progressive equation non-linear inversion of pore media | |
RU2740597C1 (en) | Method for forecasting change in production rate of production wells during propagation of elastic vibrations in bottomhole formation zone | |
Lin et al. | AVO forwarding modeling in two-phase media: Multiconstrained matrix mineral modulus 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 |