CN107290782A - Reservoir porosity, water saturation and shale content parameter Simultaneous Inversion new method - Google Patents
Reservoir porosity, water saturation and shale content parameter Simultaneous Inversion new method Download PDFInfo
- Publication number
- CN107290782A CN107290782A CN201610191787.3A CN201610191787A CN107290782A CN 107290782 A CN107290782 A CN 107290782A CN 201610191787 A CN201610191787 A CN 201610191787A CN 107290782 A CN107290782 A CN 107290782A
- Authority
- CN
- China
- Prior art keywords
- reservoir
- parameter
- water saturation
- shale content
- porosity
- 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
- 238000000034 method Methods 0.000 title claims abstract description 58
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 title claims abstract description 49
- 239000011435 rock Substances 0.000 claims abstract description 53
- 238000004458 analytical method Methods 0.000 claims abstract description 8
- 230000006870 function Effects 0.000 claims description 43
- 238000010008 shearing Methods 0.000 claims description 19
- 238000004422 calculation algorithm Methods 0.000 claims description 17
- 230000002068 genetic effect Effects 0.000 claims description 12
- 238000009826 distribution Methods 0.000 claims description 6
- 239000011159 matrix material Substances 0.000 claims description 5
- 238000003860 storage Methods 0.000 claims description 5
- 241001274660 Modulus Species 0.000 claims description 4
- 238000004364 calculation method Methods 0.000 claims description 4
- 238000000605 extraction Methods 0.000 claims description 4
- 239000004576 sand Substances 0.000 claims description 4
- 239000011148 porous material Substances 0.000 claims description 3
- 239000002023 wood Substances 0.000 claims description 3
- 235000008331 Pinus X rigitaeda Nutrition 0.000 claims description 2
- 235000011613 Pinus brutia Nutrition 0.000 claims description 2
- 241000018646 Pinus brutia Species 0.000 claims description 2
- 239000004575 stone Substances 0.000 claims description 2
- 235000013399 edible fruits Nutrition 0.000 claims 1
- 239000004744 fabric Substances 0.000 claims 1
- 238000011160 research Methods 0.000 abstract description 8
- 230000000192 social effect Effects 0.000 abstract description 2
- 239000012530 fluid Substances 0.000 description 24
- 238000005516 engineering process Methods 0.000 description 9
- 230000000694 effects Effects 0.000 description 7
- 238000011161 development Methods 0.000 description 5
- 230000000704 physical effect Effects 0.000 description 4
- 241000208340 Araliaceae Species 0.000 description 2
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 2
- 235000003140 Panax quinquefolius Nutrition 0.000 description 2
- 238000013459 approach Methods 0.000 description 2
- 238000006243 chemical reaction Methods 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 235000008434 ginseng Nutrition 0.000 description 2
- 238000005457 optimization Methods 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 229920006395 saturated elastomer Polymers 0.000 description 2
- 239000004215 Carbon black (E152) Substances 0.000 description 1
- 238000000342 Monte Carlo simulation Methods 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000000205 computational method Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 239000003989 dielectric material Substances 0.000 description 1
- 230000002349 favourable effect Effects 0.000 description 1
- 229930195733 hydrocarbon Natural products 0.000 description 1
- 150000002430 hydrocarbons Chemical class 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000012827 research and development Methods 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 241000894007 species Species 0.000 description 1
- 238000000528 statistical test Methods 0.000 description 1
- 238000002948 stochastic simulation Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 230000010415 tropism 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
- G01V1/306—Analysis for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/62—Physical property of subsurface
- G01V2210/624—Reservoir parameters
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/62—Physical property of subsurface
- G01V2210/624—Reservoir parameters
- G01V2210/6244—Porosity
-
- 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/6246—Permeability
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)
Abstract
The present invention provides a kind of reservoir porosity, water saturation and shale content parameter Simultaneous Inversion new method, and the reservoir porosity, water saturation and shale content parameter Simultaneous Inversion new method include:Step 1, collection elastic parameter, log data and core data information;Step 2, Rock Elastic Parameters are combined with reservoir physical parameter, set up petrophysical model;Step 3, the object function of reservoir porosity, water saturation and the parametric inversion of shale content parameter three is set up, and solves object function;Step 4, by solving object function, output reservoir porosity, water saturation and shale content parametric inversion result.Rock physicses and geostatistical analysis are combined by this method, carry out the research that reservoir properties recognize new method, the physical parameters such as porosity, the saturation degree of reservoir are obtained, the objectivity and accuracy of reservoir properties estimated result are improved, with important economy and social effect.
Description
Technical field
The present invention relates to the reservoir geophysics field of geophysical exploration, a kind of reservoir is especially related to
Porosity, water saturation and shale content parameter Simultaneous Inversion new method.
Background technology
Poststack wave impedance inversion technique starts to occur in 1970s, at that time to the research of seismic inversion
Simply based on the one-dimensional wave impedance inversion of poststack based on convolution model.Extensively should it be obtained in the nineties
With treatment technology is gradually improved and moved to maturity in substantial amounts of practice, around one-dimensional wave impedance inversion
All kinds of algorithms and application achievements emerge in an endless stream.But with the development of wave impedance inversion technique, people are gradually
Recognize some defects and deficiency present in wave impedance inversion:Using the loss of full angle multiple stacking and mould
The information of some reflection lithology and oil-gas possibility has been pasted, seismic data reflection reservoir characteristic change is weakened
Sensitiveness, causes poststack seismic inversion data and the application of seismic properties data to there is multiresolution issue;Separately
Outside, wave impedance, which is drilled, can only provide the parameters such as the seldom P-wave impedance of species, in Study In Reservoir physical property and stream
Be restricted in terms of body etc..
With the progress of computer technology and going deep into for exploration and development process, geophysical work person wish from
Subsurface information as much as possible is extracted in abundant seismic data, this is accomplished by Prestack seismic data inverting
Studied.Since 1980s, people constantly explore and develop Prestack seismic data inverting
Theory and method, AVO inversion techniques are one of which.AVO technologies are with skew according to reflected amplitude
Away from the subsurface lithologic that is reflected of changing rule and its property of pore-fluid carry out Direct Prediction of Oil
With a technology of estimation formation lithology parameter.The theoretical foundation of AVO invertings is Zoeppritz equations, is led to
A variety of reduced forms can be obtained by crossing the simplified equation, can be obtained using linearly or nonlinearly inversion method
The parameter that many post-stack inversions can not be obtained, especially to reservoir properties and the more sensitive ginseng of change of fluid
Number.Therefore, AVO invertings are gradually applied in production with its distinctive advantage, and achieve good answer
Use effect.
In recent years, people are had made intensive studies based on P-wave And S data to fluid identification of reservoir.
Smith and Gidlow (1987) propose to be folded by different weights function using pre stack data first
Plus obtain fluid factor and pseudo- Poisson's-ratio section to predict lithology and fluid.Goodway et al. (1997) is carried
A kind of technology recognized for Fluid Anomalies, i.e. lambda-mu-rho technologies are gone out.Hedlin(2000)
The P-wave And S impedance information that works based on Murphy et al. (1993) has been incorporated in earthquake simultaneously, is proposed
The method of hole modulus.Hilterman (2001) describes the concept of fluid factor and summarized
Goodway and Hedlin et al. achievement.If it is considered that during the rock of multibore fluid saturation, situation becomes
Considerably complicated, Biot (1941) and Gassmann (1951) consider this problem respectively.Krief(1990)
Et al. point out, both approaches are all derived identical conclusion.Before Russell et al. (2003) is summarized
The viewpoint of people, is entered using Biot-Gassmann equations to the velocity of longitudinal wave equation under saturation fluid condition
Rewriting is gone, has obtained fluid factor and carry out fluid identification.Gidlow and Smith (2003) are according to prestack AVO
Analysis, it is proposed that the concept of fluid factor angle and cross plot angle, is calculated using both angles
To obtain fluid factor.Ning Zhonghua etc. (2006) is on the basis of analysis and summary forefathers' method, it is proposed that high
The fluid factor method of sensitivity.Mark etc. (2006) proposes the concept of Poisson impedance, this concept connection
Poisson's ratio and density attributes have been closed, fluid can be more effectively distinguished than single Poisson's ratio or density parameter.
Zhang Guangzhi etc. (2011) makes full use of the difference between them, fluid to believe from angular-trace gather data
Breath and difference of the framework information on yardstick and direction, and the multiple dimensioned property of Curvelet conversion and multi-party
Tropism, it is proposed that the method that application angle fluid factor attribute carries out fluid identification.Malong etc. (2011)
Gas saturation method is calculated based on the Russell fluid factors proposed, and then evaluates loose sand reservoir
Gas-bearing property.Stream based on pre-stack seismic ripple modulus direct inversion in length and breadth in (2012) such as the emerging credits of print are proposed
Body detecting method.In terms of reservoir properties research, at present more than be angle from statistical learning, lead to
Cross what the method for supervised learning and stochastic simulation was realized, such as many attribution inversions and geostatistical inversion,
The shortcoming of these methods is that geology and geophysical theoretical foundation are not enough, more a kind of simply mathematics
Conversion, when leaving the control point for exporting statistical relationship, statistical relationship can become unreliable, thing
Property identification result can also become it is unreliable.
In a word, with the continuous progress of seismic technology, the field of seismic study is extended to from the exploration phase
Development phase, lithology classification and fluid identification are deep into by structure interpretation, and prestack inversion and physical properties of rock
The extraction and Fluid Identification Method research of parameter play the effect of key wherein.For this, we have invented one
Kind new reservoir porosity, water saturation and shale content parameter Simultaneous Inversion new method, solve with
Upper technical problem.
The content of the invention
Rock physicses and geostatistical analysis are combined it is an object of the invention to provide one kind, carry out reservoir
Physical property recognizes the research of new method, obtains the reservoir pore space of the physical parameters such as porosity, the saturation degree of reservoir
Degree, water saturation and shale content parameter Simultaneous Inversion new method.
The purpose of the present invention can be achieved by the following technical measures:Reservoir porosity, water saturation and
Shale content parameter Simultaneous Inversion new method, the reservoir porosity, water saturation and shale content parameter
Simultaneous Inversion new method includes:Step 1, collection elastic parameter, log data and core data information;
Step 2, Rock Elastic Parameters are combined with reservoir physical parameter, set up petrophysical model;Step 3,
The object function of reservoir porosity, water saturation and the parametric inversion of shale content parameter three is set up, and is asked
Solve object function;Step 4, by solving object function, output reservoir porosity, water saturation and
Shale content parametric inversion result.
The purpose of the present invention can be also achieved by the following technical measures:
In step 1, elastic parameter is obtained by prestack inversion, prestack inversion is by analyzing prestack road
Collection carrys out the earthquake analysis method of Study of The Underground lithology, by means of Zoeppritz equations or its approximate expression, most
Elastic parameter is finally inversed by eventually, the bulk modulus of the elastic parameter being finally inversed by including rock, modulus of shearing and close
Spend information.
In step 1, log data includes these letters of reservoir porosity, water saturation and shale content
Breath;Core data provides the auxiliary information of log data, including rock type, porosity.
In step 2, first, Rock Elastic Parameters are combined with reservoir physical parameter, set up and determine
The petrophysical model of property;Secondly, for the reservoir of different Storage categories, with certainty rock physicses mould
Based on type, by statistical information, the petrophysical model of statistics is set up, statistics petrophysical model is used
To describe all possible reservoir conditions realization not embodied in well logging.
In step 2, shown in the petrophysical model set up such as formula (1):
[K, μ, ρ]=fRPM(φ,sw,C)+ε (1)
Wherein, fRPMCertainty petrophysical model is represented, K, μ, ρ represents Rock Elastic Parameters, including body
Product module amount, modulus of shearing, density, φ, sw, C represents reservoir physical parameter, including porosity, aqueous full
With degree, shale content, ε represents the error of actual observation data.
In step 2, the statistics petrophysical model of foundation is to set up reservoir porosity, water saturation
Relational expression between bulk modulus, modulus of shearing and the density of shale content and rock;Tight sand mould
Inter-granular contact of the type based on Hertz-Mindlin is theoretical, and the model is estimated the bulk modulus K of dry rock and cut
Shear modulu μ;The bulk modulus and modulus of shearing of Rock Matrix are averagely obtained by Voigt-Reuss-Hill;
The bulk modulus and modulus of shearing of dry rock are obtained by improved Hashin-Shtrikman models;Mixed flow
The bulk modulus of body is given by Wood formula, the theoretical p-and s-wave velocity such as public affairs obtained by petrophysical model
Shown in formula (2):
Wherein, Vp represents velocity of longitudinal wave, and Vs represents shear wave velocity, and K represents bulk modulus, and μ represents to cut
Shear modulu, ρ represents density.
In step 3, the object function of physical parameter extraction, object function are set up according to bayesian algorithm
The maximum a posteriori probability for being reservoir physical parameter under known elasticity Parameter Conditions distribution, is expressed as:
[φ, sw, C]=arg Max { P ([φ, sw, C]j|[K,μ,ρ])} (3)
Wherein, φ, sw, C is respectively porosity, water saturation, shale content, and P () represents that probability is close
Degree, arg Max represent maximizing, and K represents bulk modulus, and μ represents modulus of shearing, and ρ represents density,
P([φ,sw,C]j| [K, μ, ρ]) be reservoir physical parameter Posterior distrbutionp, the meaning of the formula is, it is known that
This elastic parameter, calculates porosity, water saturation and shale content and belongs to each reservoir thing respectively
The Posterior probability distribution of property parameter class;When posterior probability takes maximum, the reservoir physical parameter is most
Whole inversion result;According to actual conditions, inverting target reservoir physical parameter can be porosity, aqueous full
With any combination of degree and shale content, Rock Elastic Parameters can also be other elastic parameters, such as moor
Pine ratio, Lame parameter.
In step 3, solution object function is combined by genetic algorithm and DSMC, first
Generation initialization colony, and the fitness function of each individual is calculated, if result is satisfied with, heredity is calculated
Method terminates, if result is dissatisfied, updates the colony of initialization, is emulated by Monte Carlo, and
The fitness function of each individual is calculated again, colony, so circulation is updated, until result of calculation satisfaction
Untill.
In step 4, by solving object function, at the same obtain reservoir porosity, water saturation and
Three parameters of shale content, three parameters are exported, and carry out further layer description.
Reservoir porosity, water saturation and shale content parameter Simultaneous Inversion new method in the present invention,
The angle analyzed from pre-stack seismic, rock physicses and geostatistical analysis are combined, pass through rock
Physics sets up the relation between porosity, water saturation and elastic parameter, and theory is estimated using Bayes
The data such as prior information, likelihood function, elastic parameter are effectively combined, and derives reservoir porosity and contains
The object function of inverting while water saturation, is combined finally by genetic algorithm and Monte Carlo simulation
Method solve the object function, it is final to obtain reservoir porosity, water saturation and shale simultaneously and contain
Three parameters are measured, more information are provided for the Comprehensive Study of Reservoir, the objective of reservoir properties estimated result is improved
Property and accuracy, reduce investment risk, for geological personnel estimation oil and gas reserves, determine that well location is provided more
Reliable foundation, with important economy and social effect.
Brief description of the drawings
Fig. 1 is newly square for reservoir porosity, water saturation and the shale content parameter Simultaneous Inversion of the present invention
The flow chart of one specific embodiment of method;
Fig. 2 is the structure flow chart of petrophysical model in a specific embodiment of the invention;
Fig. 3 is Bayes's parameter Estimation flow chart in the specific embodiment of the present invention;
Fig. 4 is that genetic algorithm in the specific embodiment of the present invention and DSMC are combined solution mesh
Scalar functions flow chart.
Embodiment
For enable the present invention above and other objects, features and advantages become apparent, it is cited below particularly go out
Preferred embodiment, and coordinate shown in accompanying drawing, it is described in detail below.
As shown in figure 1, Fig. 1 is same for the reservoir porosity of the present invention, water saturation and shale content parameter
When inverting new method flow chart.
In step 101, collection elastic parameter, log data and core data information.
Wherein elastic parameter is obtained by prestack inversion, prestack inversion be it is a kind of by analyze prestack trace gather come
The earthquake analysis method of Study of The Underground lithology, by means of Zoeppritz equations or its approximate expression, finally may be used
To be finally inversed by elastic parameter, the bulk modulus, modulus of shearing and density of rock is mainly used to believe here
Breath.Log data mainly includes the information, porosity such as reservoir porosity, water saturation and shale content
Data will be used when calculating Rock Matrix modulus, skeleton modulus of elasticity, fluid-mixing modulus, aqueous
Saturation degree is that, for calculating fluid-mixing modulus, shale content data are for calculating Rock Matrix modulus.
The effect of core data is to provide the auxiliary information of log data, such as rock type, porosity etc..
In step 102, petrophysical model is set up.
First, Rock Elastic Parameters are combined with reservoir physical parameter, set up deterministic rock physicses
Model.Secondly, for the reservoir of different Storage categories, based on certainty petrophysical model, lead to
Statistical information is crossed, the petrophysical model of statistics is set up.Count one important effect of petrophysical model
The all possible reservoir conditions for being exactly used for describing not embody in well logging are realized.
One of main purpose of seismic interpretation is to determine that can aqueous saturated rock or hydrocarbonaceous saturated rock
Produce favourable reflection.In order to complete this task, it is necessary to estimate water saturation state and hydrocarbon saturation state
Between rock behavio(u)r difference.Thus need to set up some basic rock physics relations.These relations
Including the empirical equation and theory relation between rock behavio(u)r and elastic parameter, it is also required to use ripple biography in addition
Broadcast model.The achievement in research of petrophysical property mainly be from seismic data extract subsurface rock and
The property of its saturation fluid establishes physical basis.On the other hand, seismic wave characteristics and rock, fluid are understood
The relation of property, helps to simulate propagation law of the seismic wave in complex dielectrics.Rock physicses research
Main purpose be understand lithology, porosity, pore pressure, fluid type, saturation degree, anisotropy,
Temperature and frequency etc. are to compressional wave, shear wave velocity and the influence of decay in rock.
First, the Rock Elastic Parameters that prestack inversion is obtained are combined with reservoir physical parameter, set up true
Qualitatively petrophysical model.Secondly as subsurface reservoir condition is complicated and changeable, deterministic rock
Stone physical model is difficult to the situation of accurate simulation underground different reservoir, for the reservoir of different Storage categories,
Based on certainty petrophysical model, by statistical information, statistics petrophysical model is set up.System
Meter petrophysical model one it is important effect be exactly be used for describe well logging on do not embody be possible to
Reservoir conditions realize that constructed petrophysical model such as formula (1) is shown.
[K, μ, ρ]=fRPM(φ,sw,C)+ε (1)
Wherein, fRPMCertainty petrophysical model is represented, K, μ, ρ represents Rock Elastic Parameters (such as volume
Modulus, modulus of shearing, density), φ, sw, C represent reservoir physical parameter (porosity, water saturation,
Shale content), ε represents the error of actual observation data.
For the reservoir of different Storage categories, based on certainty petrophysical model, believed by counting
Breath, sets up statistics petrophysical model.Reservoir porosity, water saturation and shale is mainly set up to contain
Amount and the relational expression between bulk modulus, modulus of shearing and the density of rock.Tight sand model is based on
Hertz-Mindlin inter-granular contact is theoretical, and the model can estimate bulk modulus K and the shearing of dry rock
Modulus μ.The bulk modulus and modulus of shearing of Rock Matrix are averagely obtained by Voigt-Reuss-Hill.It is dry
The bulk modulus and modulus of shearing of rock are obtained by improved Hashin-Shtrikman models.Fluid-mixing
Bulk modulus given by Wood formula, the theoretical p-and s-wave velocity such as formula obtained by petrophysical model
(2) shown in, it is as shown in Figure 2 that petrophysical model builds flow.
Wherein, Vp represents velocity of longitudinal wave, and Vs represents shear wave velocity, and K represents bulk modulus, and μ represents to cut
Shear modulu, ρ represents density.
In step 103, the object function of three parametric inversions is set up.
The two of the core of this patent are establish reservoir porosity, water saturation and shale content same
When inverting object function, so that the elastic parameter of three parameters and rock is set up into a kind of contact, the mesh
Scalar functions are the key and core of three parameter Simultaneous Inversions.Different reservoir, the prior information of its reservoir properties
It is different, meanwhile, likelihood function is also different, accordingly, it would be desirable to be set up according to different reservoirs
Different inversion for physical properties object functions.One of key point is the acquisition of prior distribution information, generally assumes that elder generation
Test information and obey gauss hybrid models, the two of key point is the acquisition of likelihood function, generally assumes that likelihood letter
Number obeys single Gaussian Profile.
In mathematical statistics field, because viewpoint is different, various schools are formd.Its Main Schools have early in
The classical school that there is in 19th century and Bayesian schools.Generally, familiar classical school theory and
Method, such as least square method, point estimation, maximal possibility estimation.The application of these theoretical and methods
Quite extensively, but as statistics is widely used in the fields such as natural science, economic research, people
Gradually it is found that the application value of bayesian theory.By the research and development of decades, Bayesian schools
Formed and develop into one it is being made a strong impact in statistics, may be with Frequency school side by side
Group, it is seen that, bayes method is strictly current statistical study hotspot.Bayes method can be with
Say that prior information and likelihood function are effectively combined, maximized by posteriority parameter, carry out the estimation of parameter,
As shown in Figure 3.
The object function of physical parameter extraction is set up according to bayesian algorithm, object function is joined for reservoir properties
Maximum a posteriori probability distribution of the number under known elasticity Parameter Conditions, is expressed as:
[φ, sw, C]=arg Max { P ([φ, sw, C]j|[K,μ,λ,ρ])} (3)
Wherein, P ([φ, sw, C]j| [K, μ, λ, ρ]) be reservoir physical parameter Posterior distrbutionp.The meaning of the formula
Justice is, it is known that this elastic parameter, and porosity, water saturation and shale content are calculated respectively and is belonged to
The Posterior probability distribution of each reservoir physical parameter class.When posterior probability takes maximum, the reservoir properties
Parameter is final inversion result.According to actual conditions, inverting target reservoir physical parameter can be hole
Any combination of degree, water saturation and shale content, Rock Elastic Parameters can also be other elasticity ginsengs
Number, such as Poisson's ratio, Lame parameter etc..
In step 104, object function is solved.
The three of the core of this patent are, due to object function have it is very strong non-linear, by traditional excellent
Change method hardly results in the solution of problem, and this patent has used the new approaches that a kind of object function is solved, will covered
The advantage of special Caro method and genetic algorithm is combined solution object function, can further improve three parameters
The efficiency and stability of inverting.
Because object function has very strong non-linear, problem is hardly resulted in by traditional optimization method
Solution, accordingly, it would be desirable to study new thinking to carry out the solution of object function, is calculated present invention employs heredity
Method that method and DSMC are combined solves object function.
Genetic algorithm begins at the eighties of last century sixties.It is initially that some biologists utilize computer pair
Genetic system is simulated.University of Michigan of U.S. professor Holland is opened by biologist's analog result
Hair, first Application simulates operator to study adaptability problem.Subsequent Holland develops a kind of programming skill
Art, its basic thought is the mode using similar natural selection come the program of designing a calculating machine.In Software for Design
In, adaptability is the different choice based on monitor program, the program not good by constantly rejecting effect, is allowed
The good program of those Solve problems is more and more dominant, so that system finally adapts to arbitrary environment.
Later this algorithm is applied in wider problem by more scholars, such as function optimization and categorizing system
Deng, and achieve good effect.
DSMC is also known as random sampling skill or statistical test method.Since over half a century, by
In the development and the invention of electronic computer of science and technology, this method is carried as a kind of independent method
Out, and first applied in the experiment and development of nuclear weapon.DSMC is a kind of meter
Calculation method, but have very big difference with prevailing value computational methods.It is based on Probability Statistics Theory
A kind of method.The characteristics of things can more realistically be described due to DSMC and Physical Experiment
Journey, solves the insoluble problem of some numerical methods, thus the application field of this method is increasingly extensive.
DSMC solves the approximate solution of mathematics or physical problem by way of stochastical sampling, it
Advantage be algorithmic stability, have the disadvantage that efficiency is low, precision is poor.Genetic algorithm is a kind of global optimizing
Algorithm, advantage is can to obtain globally optimal solution, have the disadvantage for it is very complicated the problem of, the result of solution
It is unstable.Therefore, genetic algorithm and DSMC are combined, solve such nonlinear problem,
Further improve the efficiency and stability of solution.Genetic algorithm and DSMC are combined solution target letter
Number flow is as shown in Figure 4.Initialization colony, and the fitness function of each individual of calculating are firstly generated,
If result is satisfied with, genetic algorithm terminates, if result is dissatisfied, updates the colony of initialization, leads to
Cross Monte Carlo to be emulated, and calculate the fitness function of each individual again, update colony, so
Circulation, untill result of calculation is satisfied with.
In step 105, three parametric inversion results are exported.
By solving object function, it is possible to contain while obtaining reservoir porosity, water saturation and shale
Three parameters are measured, are now exported three parameters, it is possible to further layer description is carried out, so as to carry
The objectivity and accuracy of high reservoir properties estimated result, reduce investment risk, are geological personnel estimation oil
Gas reserves, determine that well location provides more structurally sound foundation.
Claims (9)
1. reservoir porosity, water saturation and shale content parameter Simultaneous Inversion new method, it is characterised in that
The reservoir porosity, water saturation and shale content parameter Simultaneous Inversion new method include:
Step 1, collection elastic parameter, log data and core data information;
Step 2, Rock Elastic Parameters are combined with reservoir physical parameter, set up petrophysical model;
Step 3, the mesh of reservoir porosity, water saturation and the parametric inversion of shale content parameter three is set up
Scalar functions, and solve object function;
Step 4, by solving object function, output reservoir porosity, water saturation and shale content
Parametric inversion result.
2. reservoir porosity according to claim 1, water saturation and shale content parameter Simultaneous Inversion
New method, it is characterised in that in step 1, elastic parameter, prestack inversion are obtained by prestack inversion
It is by analyzing prestack trace gather come the earthquake analysis method of Study of The Underground lithology, by means of Zoeppritz side
Journey or its approximate expression, are finally finally inversed by elastic parameter, and the elastic parameter being finally inversed by includes the volume mould of rock
Amount, modulus of shearing and density information.
3. reservoir porosity according to claim 1, water saturation and shale content parameter Simultaneous Inversion
New method, it is characterised in that in step 1, log data includes reservoir porosity, water saturation
With these information of shale content;Core data provides the auxiliary information of log data, including rock type,
Porosity.
4. reservoir porosity according to claim 1, water saturation and shale content parameter Simultaneous Inversion
New method, it is characterised in that in step 2, first, by Rock Elastic Parameters and reservoir physical parameter
It is combined, sets up deterministic petrophysical model;Secondly, for the reservoir of different Storage categories, with
Based on certainty petrophysical model, by statistical information, the petrophysical model of statistics, system are set up
The all possible reservoir conditions that meter petrophysical model is used for describing not embody in well logging are realized.
5. reservoir porosity according to claim 4, water saturation and shale content parameter Simultaneous Inversion
New method, it is characterised in that in step 2, shown in the petrophysical model set up such as formula (1):
[K, μ, ρ]=fRPM(φ,sw,C)+ε (1)
Wherein, fRPMCertainty petrophysical model is represented, K, μ, ρ represents Rock Elastic Parameters, including body
Product module amount, modulus of shearing, density, φ, sw, C represents reservoir physical parameter, including porosity, aqueous full
With degree, shale content, ε represents the error of actual observation data.
6. reservoir porosity according to claim 4, water saturation and shale content parameter Simultaneous Inversion
New method, it is characterised in that in step 2, the statistics petrophysical model of foundation is to set up reservoir hole
Pass between bulk modulus, modulus of shearing and the density of porosity, water saturation and shale content and rock
It is formula;Inter-granular contact of the tight sand model based on Hertz-Mindlin is theoretical, and the model estimates dry rock
The bulk modulus K and modulus of shearing μ of stone;The bulk modulus and modulus of shearing of Rock Matrix by
Voigt-Reuss-Hill is averagely obtained;The bulk modulus and modulus of shearing of dry rock are by improved
Hashin-Shtrikman models are obtained;The bulk modulus of fluid-mixing is given by Wood formula, by rock
Shown in the theoretical p-and s-wave velocity such as formula (2) that physical model is obtained:
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>V</mi>
<mi>P</mi>
</msub>
<mo>=</mo>
<msqrt>
<mfrac>
<mrow>
<mi>K</mi>
<mo>+</mo>
<mfrac>
<mn>3</mn>
<mn>4</mn>
</mfrac>
<mi>&mu;</mi>
</mrow>
<mi>&rho;</mi>
</mfrac>
</msqrt>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>V</mi>
<mi>S</mi>
</msub>
<mo>=</mo>
<msqrt>
<mfrac>
<mi>&mu;</mi>
<mi>&rho;</mi>
</mfrac>
</msqrt>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
Wherein, Vp represents velocity of longitudinal wave, and Vs represents shear wave velocity, and K represents bulk modulus, and μ represents to cut
Shear modulu, ρ represents density.
7. reservoir porosity according to claim 1, water saturation and shale content parameter Simultaneous Inversion
New method, it is characterised in that in step 3, the mesh of physical parameter extraction is set up according to bayesian algorithm
Scalar functions, object function is the reservoir physical parameter maximum a posteriori probability under known elasticity Parameter Conditions point
Cloth, is expressed as:
[φ, sw, C]=argMax { P ([φ, sw, C]j|[K,μ,ρ])} (3)
Wherein, φ, sw, C is respectively porosity, water saturation, shale content, and P () represents that probability is close
Degree, argMax represents maximizing, and K represents bulk modulus, and μ represents modulus of shearing, and ρ represents density,
P([φ,sw,C]j| [K, μ, ρ]) be reservoir physical parameter Posterior distrbutionp, the meaning of the formula is, it is known that
This elastic parameter, calculates porosity, water saturation and shale content and belongs to each reservoir thing respectively
The Posterior probability distribution of property parameter class;When posterior probability takes maximum, the reservoir physical parameter is most
Whole inversion result;According to actual conditions, inverting target reservoir physical parameter can be porosity, aqueous full
With any combination of degree and shale content, Rock Elastic Parameters can also be other elastic parameters, such as moor
Pine ratio, Lame parameter.
8. reservoir porosity according to claim 1, water saturation and shale content parameter Simultaneous Inversion
New method, it is characterised in that in step 3, is combined by genetic algorithm and DSMC and asked
Object function is solved, initialization colony is firstly generated, and calculates the fitness function of each individual, if knot
Fruit is satisfied, then genetic algorithm terminates, if result is dissatisfied, updates the colony of initialization, special by covering
Caro is emulated, and calculates the fitness function of each individual again, updates colony, so circulation,
Untill result of calculation is satisfied with.
9. reservoir porosity according to claim 1, water saturation and shale content parameter Simultaneous Inversion
New method, it is characterised in that in step 4, by solving object function, while obtaining reservoir pore space
Degree, three parameters of water saturation and shale content, three parameters are exported, and carry out further reservoir
Description.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610191787.3A CN107290782B (en) | 2016-03-30 | 2016-03-30 | Reservoir porosity, water saturation and shale content parameter Simultaneous Inversion new method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610191787.3A CN107290782B (en) | 2016-03-30 | 2016-03-30 | Reservoir porosity, water saturation and shale content parameter Simultaneous Inversion new method |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107290782A true CN107290782A (en) | 2017-10-24 |
CN107290782B CN107290782B (en) | 2019-07-12 |
Family
ID=60087066
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610191787.3A Active CN107290782B (en) | 2016-03-30 | 2016-03-30 | Reservoir porosity, water saturation and shale content parameter Simultaneous Inversion new method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107290782B (en) |
Cited By (21)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108548765A (en) * | 2018-03-27 | 2018-09-18 | 中国石油天然气集团有限公司 | A kind of porosity calculation method for the soil skeleton that becomes sticky |
CN110333534A (en) * | 2019-06-19 | 2019-10-15 | 中国地质大学(武汉) | A kind of Bayes's time shift AVO inversion method and system based on Biot theory |
CN110554064A (en) * | 2019-07-26 | 2019-12-10 | 中国石油大学(华东) | Method for accurately estimating hydrate saturation in marine sediment based on dielectric properties |
CN110580390A (en) * | 2019-09-04 | 2019-12-17 | 电子科技大学 | geostatistical stochastic inversion method based on improved genetic algorithm and information entropy |
CN111323823A (en) * | 2019-12-27 | 2020-06-23 | 中国石油天然气股份有限公司 | Method and system for determining logging porosity curve |
CN111399042A (en) * | 2020-03-20 | 2020-07-10 | 清华大学 | Reservoir physical property parameter prediction method and electronic equipment |
CN111856583A (en) * | 2020-06-12 | 2020-10-30 | 中国石油天然气股份有限公司 | Reservoir water saturation quantitative prediction method based on porosity constraint condition |
CN111897011A (en) * | 2020-07-27 | 2020-11-06 | 中国石油大学(华东) | Reservoir pore characteristic determination method, device and equipment |
CN112149282A (en) * | 2020-08-28 | 2020-12-29 | 中国石油天然气集团有限公司 | Physical calculation method and system for natural gas hydrate saturation rock in well |
CN112415616A (en) * | 2019-08-22 | 2021-02-26 | 中国石油天然气集团有限公司 | Deep-buried reservoir porosity inversion method and device |
CN113031068A (en) * | 2021-02-24 | 2021-06-25 | 浙江大学 | Reflection coefficient accurate base tracking prestack seismic inversion method |
CN113253347A (en) * | 2021-05-14 | 2021-08-13 | 中南大学 | Shale reservoir AVO inversion characterization method and system based on VTI medium |
CN113433587A (en) * | 2020-03-23 | 2021-09-24 | 中国石油天然气股份有限公司 | Rock physics inversion method and system |
CN113779897A (en) * | 2021-08-05 | 2021-12-10 | 中国石油大学(华东) | Method and device for predicting physical property parameters of hydrate reservoir and storage medium |
CN113970787A (en) * | 2020-07-22 | 2022-01-25 | 中国石油化工股份有限公司 | Physical property parameter inversion method, physical property parameter inversion device, computer equipment and storage medium |
CN114563821A (en) * | 2020-11-27 | 2022-05-31 | 中国石油天然气股份有限公司 | Intelligent seismic wave impedance inversion label data manufacturing method and device |
CN115079261A (en) * | 2022-06-06 | 2022-09-20 | 吉林大学 | Compact sandstone gas reservoir evaluation method based on multi-parameter quantitative interpretation template |
CN115586573A (en) * | 2022-09-15 | 2023-01-10 | 河海大学 | Dynamic constraint physical property parameter seismic inversion method for tight sandstone reservoir |
US11578596B2 (en) | 2021-07-08 | 2023-02-14 | Saudi Arabian Oil Company | Constrained natural fracture parameter hydrocarbon reservoir development |
CN116165701A (en) * | 2022-12-08 | 2023-05-26 | 中国矿业大学 | Pore structure parameter pre-stack seismic direct inversion method for carbonate reservoir |
CN116299696A (en) * | 2023-02-08 | 2023-06-23 | 中海石油(中国)有限公司深圳分公司 | Method for simultaneously and quantitatively predicting hydrocarbon source rock shale content, TOC, porosity and water saturation |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104516017A (en) * | 2013-09-29 | 2015-04-15 | 中国石油化工股份有限公司 | Carbonate rock physical parameter seismic inversion method |
CN104570127A (en) * | 2013-10-09 | 2015-04-29 | 中国石油化工股份有限公司 | Method for simultaneously inverting porosity and shaliness by using seismic wave velocity |
-
2016
- 2016-03-30 CN CN201610191787.3A patent/CN107290782B/en active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104516017A (en) * | 2013-09-29 | 2015-04-15 | 中国石油化工股份有限公司 | Carbonate rock physical parameter seismic inversion method |
CN104570127A (en) * | 2013-10-09 | 2015-04-29 | 中国石油化工股份有限公司 | Method for simultaneously inverting porosity and shaliness by using seismic wave velocity |
Non-Patent Citations (6)
Title |
---|
DARIO GRANA ET AL.: "Probabilistic petrophysical-properties estimation integrating statistical rock physics with seismic inversion", 《GEOPHYSICS》 * |
XINGYAO YIN ET AL.: "A Novel Prestack AVO Inversion and its Application", 《SEG LAS VEGAS 2008 ANNUAL MEETING》 * |
印兴耀等: "基于弹性阻抗的储层物性参数预测方法", 《地球物理学报》 * |
王小丹等: "统计岩石物理在储层物性反演中的应用", 《中国地球科学联合学术年会2014》 * |
胡华锋等: "基于贝叶斯分类的储层物性参数联合反演方法", 《石油物探》 * |
谭茂金: "《有机页岩测井岩石物理》", 30 November 2015, 北京:石油工业出版社 * |
Cited By (34)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108548765B (en) * | 2018-03-27 | 2020-11-06 | 中国石油天然气集团有限公司 | Porosity calculation method for clay-changing framework |
CN108548765A (en) * | 2018-03-27 | 2018-09-18 | 中国石油天然气集团有限公司 | A kind of porosity calculation method for the soil skeleton that becomes sticky |
CN110333534B (en) * | 2019-06-19 | 2020-10-30 | 中国地质大学(武汉) | Bayes time-shifting AVO inversion method and system based on Biot theory |
CN110333534A (en) * | 2019-06-19 | 2019-10-15 | 中国地质大学(武汉) | A kind of Bayes's time shift AVO inversion method and system based on Biot theory |
CN110554064A (en) * | 2019-07-26 | 2019-12-10 | 中国石油大学(华东) | Method for accurately estimating hydrate saturation in marine sediment based on dielectric properties |
CN110554064B (en) * | 2019-07-26 | 2021-12-21 | 中国石油大学(华东) | Method for accurately estimating hydrate saturation in marine sediment based on dielectric properties |
CN112415616A (en) * | 2019-08-22 | 2021-02-26 | 中国石油天然气集团有限公司 | Deep-buried reservoir porosity inversion method and device |
CN110580390A (en) * | 2019-09-04 | 2019-12-17 | 电子科技大学 | geostatistical stochastic inversion method based on improved genetic algorithm and information entropy |
CN110580390B (en) * | 2019-09-04 | 2021-05-11 | 电子科技大学 | Geostatistical stochastic inversion method based on improved genetic algorithm and information entropy |
CN111323823A (en) * | 2019-12-27 | 2020-06-23 | 中国石油天然气股份有限公司 | Method and system for determining logging porosity curve |
CN111323823B (en) * | 2019-12-27 | 2022-05-10 | 中国石油天然气股份有限公司 | Method and system for determining logging porosity curve |
CN111399042A (en) * | 2020-03-20 | 2020-07-10 | 清华大学 | Reservoir physical property parameter prediction method and electronic equipment |
CN111399042B (en) * | 2020-03-20 | 2021-02-02 | 清华大学 | Reservoir physical property parameter prediction method and electronic equipment |
CN113433587B (en) * | 2020-03-23 | 2023-09-26 | 中国石油天然气股份有限公司 | Rock physical inversion method and system |
CN113433587A (en) * | 2020-03-23 | 2021-09-24 | 中国石油天然气股份有限公司 | Rock physics inversion method and system |
CN111856583A (en) * | 2020-06-12 | 2020-10-30 | 中国石油天然气股份有限公司 | Reservoir water saturation quantitative prediction method based on porosity constraint condition |
CN113970787A (en) * | 2020-07-22 | 2022-01-25 | 中国石油化工股份有限公司 | Physical property parameter inversion method, physical property parameter inversion device, computer equipment and storage medium |
CN111897011A (en) * | 2020-07-27 | 2020-11-06 | 中国石油大学(华东) | Reservoir pore characteristic determination method, device and equipment |
CN111897011B (en) * | 2020-07-27 | 2024-03-22 | 中国石油大学(华东) | Reservoir pore characteristic determination method, device and equipment |
CN112149282B (en) * | 2020-08-28 | 2024-05-28 | 中国石油天然气集团有限公司 | Rock physical calculation method and system for saturation of natural gas hydrate in well |
CN112149282A (en) * | 2020-08-28 | 2020-12-29 | 中国石油天然气集团有限公司 | Physical calculation method and system for natural gas hydrate saturation rock in well |
CN114563821B (en) * | 2020-11-27 | 2024-05-28 | 中国石油天然气股份有限公司 | Intelligent seismic wave impedance inversion tag data manufacturing method and device |
CN114563821A (en) * | 2020-11-27 | 2022-05-31 | 中国石油天然气股份有限公司 | Intelligent seismic wave impedance inversion label data manufacturing method and device |
CN113031068A (en) * | 2021-02-24 | 2021-06-25 | 浙江大学 | Reflection coefficient accurate base tracking prestack seismic inversion method |
CN113253347A (en) * | 2021-05-14 | 2021-08-13 | 中南大学 | Shale reservoir AVO inversion characterization method and system based on VTI medium |
US11578596B2 (en) | 2021-07-08 | 2023-02-14 | Saudi Arabian Oil Company | Constrained natural fracture parameter hydrocarbon reservoir development |
CN113779897B (en) * | 2021-08-05 | 2022-10-18 | 中国石油大学(华东) | Method and device for predicting physical property parameters of hydrate reservoir and storage medium |
CN113779897A (en) * | 2021-08-05 | 2021-12-10 | 中国石油大学(华东) | Method and device for predicting physical property parameters of hydrate reservoir and storage medium |
CN115079261A (en) * | 2022-06-06 | 2022-09-20 | 吉林大学 | Compact sandstone gas reservoir evaluation method based on multi-parameter quantitative interpretation template |
CN115079261B (en) * | 2022-06-06 | 2023-07-11 | 吉林大学 | Compact sandstone gas reservoir evaluation method based on multiparameter quantitative interpretation template |
CN115586573A (en) * | 2022-09-15 | 2023-01-10 | 河海大学 | Dynamic constraint physical property parameter seismic inversion method for tight sandstone reservoir |
CN116165701A (en) * | 2022-12-08 | 2023-05-26 | 中国矿业大学 | Pore structure parameter pre-stack seismic direct inversion method for carbonate reservoir |
CN116165701B (en) * | 2022-12-08 | 2024-02-20 | 中国矿业大学 | Pore structure parameter pre-stack seismic direct inversion method for carbonate reservoir |
CN116299696A (en) * | 2023-02-08 | 2023-06-23 | 中海石油(中国)有限公司深圳分公司 | Method for simultaneously and quantitatively predicting hydrocarbon source rock shale content, TOC, porosity and water saturation |
Also Published As
Publication number | Publication date |
---|---|
CN107290782B (en) | 2019-07-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107290782B (en) | Reservoir porosity, water saturation and shale content parameter Simultaneous Inversion new method | |
CN104516017B (en) | A kind of carbonate rock petrophysical parameter seismic inversion method | |
Iturrarán-Viveros et al. | Artificial neural networks applied to estimate permeability, porosity and intrinsic attenuation using seismic attributes and well-log data | |
CN105629303B (en) | Based on petrophysical pre-stack fracture quantitative forecasting technique and system | |
CN105182424B (en) | A kind of method and apparatus based on patchy saturation quantitative forecast reservoir porosity | |
CN104950331B (en) | A kind of porosity of sand mud reservoir and the earthquake prediction method of shale content | |
CN104502997B (en) | A kind of method of utilization fracture spacing curve prediction fracture spacing body | |
Anifowose et al. | Integrating seismic and log data for improved petroleum reservoir properties estimation using non-linear feature-selection based hybrid computational intelligence models | |
CN108037528B (en) | Porosity prediction method and system of few wellblock based on statistics rock physics modeling | |
Zhang et al. | Improvement of petrophysical workflow for shear wave velocity prediction based on machine learning methods for complex carbonate reservoirs | |
CN111596978A (en) | Web page display method, module and system for lithofacies classification by artificial intelligence | |
CN110333551A (en) | Dolostone reservoirs prediction technique, system and the storage medium combined based on well shake | |
Shi et al. | Brittleness index prediction in shale gas reservoirs based on efficient network models | |
CN107203005A (en) | A kind of method that quantification calculates crack characterising parameter | |
CN106501854A (en) | Oil gas drilling target geological quantifying risk evaluation methodology | |
CN104570101A (en) | AVO (amplitude versus offset) three-parameter inversion method based on particle swarm optimization | |
CN110058323A (en) | A kind of tight sand formation brittleness index calculation method | |
US20210381362A1 (en) | Method and apparatus for estimating lithofacies by learning well logs | |
CN106597543A (en) | Stratigraphic sedimentary facies division method | |
Gavidia et al. | Utilizing integrated artificial intelligence for characterizing mineralogy and facies in a pre-salt carbonate reservoir, Santos Basin, Brazil, using cores, wireline logs, and multi-mineral petrophysical evaluation | |
Hussain et al. | Machine learning-a novel approach to predict the porosity curve using geophysical logs data: An example from the Lower Goru sand reservoir in the Southern Indus Basin, Pakistan | |
Lee et al. | Interpreting the effects of shale rock properties on seismic anisotropy by statistical and machine learning methods | |
Hosseini et al. | A Meta attribute for reservoir permeability classification using well logs and 3D seismic data with probabilistic neural network. | |
Qian et al. | A rock physics driven Bayesian inversion for TOC in the Fuling shale gas reservoir | |
Malvić et al. | Neural networks in petroleum geology as interpretation tools |
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 |