CN107194104B - The method for calculating inanimate matter porosity in shale reservoir - Google Patents
The method for calculating inanimate matter porosity in shale reservoir Download PDFInfo
- Publication number
- CN107194104B CN107194104B CN201710421951.XA CN201710421951A CN107194104B CN 107194104 B CN107194104 B CN 107194104B CN 201710421951 A CN201710421951 A CN 201710421951A CN 107194104 B CN107194104 B CN 107194104B
- Authority
- CN
- China
- Prior art keywords
- formula
- porosity
- volume
- shale reservoir
- matter
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH DRILLING; MINING
- E21B—EARTH DRILLING, e.g. DEEP DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B49/00—Testing the nature of borehole walls; Formation testing; Methods or apparatus for obtaining samples of soil or well fluids, specially adapted to earth drilling or wells
Abstract
A method of inanimate matter porosity, total porosity and component in shale reservoir are calculated using log, are related to geological technique field.The method that inanimate matter porosity in shale reservoir is calculated using log, comprising: introduce assumed condition.According to the composition of shale reservoir, the volume-based model of shale reservoir is determined.The computation model of inanimate matter porosity is determined according to the volume-based model of relationship and shale reservoir between acoustic logging, density log curve and the composition of shale reservoir.A method of component in shale reservoir is calculated using log, calculated inanimate matter porosity and each component accounting in shale reservoir can be calculated using the relationship between acoustic logging, density log curve and the composition of shale reservoir.A method of total porosity in shale reservoir being calculated using log, each component accounting and inanimate matter porosity in calculated shale reservoir can calculate the total porosity in shale reservoir.
Description
Technical field
The present invention relates to geological technique fields, and in particular to a kind of side for calculating inanimate matter porosity in shale reservoir
Method.
Background technique
The rock texture composition of shale reservoir generally comprises: organic composition, inanimate matter ingredient and each interstitial space,
Wherein inanimate matter ingredient includes clay mineral content and non-clay mineral ingredient.It is more in shale reservoir compared to conventional reservoir
Organic matter (kerogen).There is free state natural gas in the hole in organic matter (kerogen), is also adsorbed in cheese root surface
There are a large amount of natural gases.Currently, first constructing volume-based model, then use hole during the conventional logging of conventional sandstone gas-bearing formation is explained at home
Porosity log is explained with Resistivity log.But shale gas Reservoir Minerals composition, porosity type and fluid are assigned
Form complexity is deposited, has very big difference with conventional gas-bearing formation.If indiscriminately imitating the log interpretation model of conventional sandstone gas-bearing formation, page can not be portrayed
The component difference of rock and conventional sandstone reservoir hinders the popularization and application of the log interpretation method based on volume-based model at the scene.
Summary of the invention
The purpose of the present invention is to provide it is a kind of using log calculate shale reservoir in inanimate matter porosity method,
It is by volume-based model and utilizes the relationship between acoustic logging, density log curve and the composition of shale reservoir can
To calculate inanimate matter porosity in shale reservoir.
The second object of the present invention is to provide a kind of method for calculating component in shale reservoir using log, utilize
Between the inanimate matter porosity and utilization acoustic logging, density log curve and the composition of shale reservoir being calculated
Relationship can calculate each component accounting in shale reservoir.
The third object of the present invention is to provide a kind of method for calculating total porosity in shale reservoir using log,
Each component accounting and inanimate matter porosity in calculated shale reservoir can calculate total in shale reservoir
Porosity.
The present invention solves its technical problem and adopts the following technical solutions to realize.
The present invention proposes a kind of method for calculating inanimate matter porosity in shale reservoir using log, comprising:
Introduce assumed condition: shale reservoir is made of inanimate matter and organic matter;It is organic containing only gas in the hole of organic matter
The porosity of matter is definite value;Containing only water in the hole of inanimate matter;The grain density of shale reservoir is definite value;
According to the composition of shale reservoir, the volume-based model of shale reservoir is determined;
According to the relationship and shale reservoir between acoustic logging, density log curve and the composition of shale reservoir
Volume-based model determines the computation model of inanimate matter porosity.
The present invention proposes a kind of method for calculating component in shale reservoir using log, using using log meter
The porosity φ for the inanimate matter that the method for calculating inanimate matter porosity in shale reservoir is calculatedpmAnd formula (1)-formula (9)
Arrangement can obtain:
By parameters and φpmIt brings formula (11) into, percentage of the organic matter volume in shale total volume is calculated
Vpk。
The present invention also proposes a kind of method for calculating total porosity in shale reservoir using log, the meter of total porosity
Calculate model are as follows: φt=Vpm×φpm+Vpk×φpk, the above-mentioned side that component in shale reservoir is calculated using log will be utilized
The V that method is calculatedpkAnd other parameters are brought into the computation model of total porosity, and the total porosity of shale reservoir is calculated
φt。
The beneficial effect of the embodiment of the present invention is: a kind of to calculate inanimate matter porosity in shale reservoir using log
Method by volume-based model and utilizes the pass between acoustic logging, density log curve and the composition of shale reservoir
System can calculate inanimate matter porosity in shale reservoir.
A method of using log calculate shale reservoir in component, calculated inanimate matter porosity with
And it can be calculated in shale reservoir using the relationship between acoustic logging, density log curve and the composition of shale reservoir
Each component accounting.
A method of total porosity in shale reservoir being calculated using log, in calculated shale reservoir
Each component accounting and inanimate matter porosity can calculate the total porosity in shale reservoir.
Specific embodiment
It in order to make the object, technical scheme and advantages of the embodiment of the invention clearer, below will be in the embodiment of the present invention
Technical solution be clearly and completely described.The person that is not specified actual conditions in embodiment, according to normal conditions or manufacturer builds
The condition of view carries out.Reagents or instruments used without specified manufacturer is the conventional production that can be obtained by commercially available purchase
Product.
Inanimate matter porosity, total pore space in shale reservoir are calculated using log to a kind of of the embodiment of the present invention below
The method of degree and component is specifically described.
A method of inanimate matter porosity in shale reservoir is calculated using log, comprising:
(1) introduce assumed condition: shale reservoir is made of inanimate matter and organic matter;Containing only gas in the hole of organic matter,
The porosity of organic matter is definite value;Containing only water in the hole of inanimate matter;The grain density of shale reservoir is definite value.
Table 1 is the composition and fluid bearing features schematic diagram of shale reservoir, more in shale reservoir compared to conventional reservoir
Organic matter (kerogen).There is free state natural gas in the hole in organic matter (kerogen), is also adsorbed in cheese root surface
A large amount of natural gases.
The component and fluid occurrence status of 1 shale of table
According to the shale composition in table 1, the component that can be obtained in shale logging module simplifies, as shown in table 2.
Component in 2 shale logging module of table simplifies
(2) according to the composition of shale reservoir, the volume-based model of shale reservoir is determined.
I.e. according to the component of the rock formation in table 2, the volume-based model of shale reservoir can be obtained are as follows:
Vpm+Vpk=1 formula (1);
Vpm=Vpmφpm+Vsm=Vpmφpm+Vpm(1-φpm) formula (2);
Vpk=Vpkφpk+Vsk=Vpkφpk+Vpk(1-φpk) formula (3);
In formula: VpkFor percentage of the organic matter volume in shale total volume;VpmIt is inanimate matter volume in shale total volume
In percentage;φpmFor the porosity of inanimate matter;φpkFor the porosity of organic matter, VsmFor inanimate matter solid volume in nothing
Percentage in machine plastid product;VskFor percentage of the volume in organic matter volume of organic matter solid.
(3) according to the relationship and shale storage between acoustic logging, density log curve and the composition of shale reservoir
The volume-based model of layer determines the computation model of inanimate matter porosity.
Relationship between acoustic logging, density log curve and shale reservoir are as follows:
ρpk=ρsk×(1-φpk)+ρpkfl×φpkFormula (4);
DTpk=DTsk×(1-φpk)+DTpkfl×φpkFormula (5);
ρpm=ρsm×(1-φpm)+ρpmfl×φpmFormula (6);
DTpm=DTsm×(1-φpm)+DTpmfl×φpmFormula (7);
ρb=ρpmflVpmφpm+ρsmVsm+ρskVsk+ρpkflVpkφpkFormula (8);
DT=DTpmflVpmφpm+DTsmVsm+DTskVsk+DTpkflVpkφpkFormula (9).
In formula, ρpkFor the density of organic matter;ρskFor the density of organic matter solid;ρpkflFor the fluid in the hole of organic matter
Density;DTpkFor the interval transit time of organic matter;DTskFor the interval transit time of organic matter solid;DTpkflTo be flowed in the hole of organic matter
Bulk acoustic wave value;ρpmFor the density of inanimate matter;ρsmFor the density of inanimate matter solid;ρpmflFor the Fluid in Pore density of inanimate matter;
DTpmFor the interval transit time of inanimate matter;DTsmFor the interval transit time of inanimate matter solid;DTpmflIt is close for the Fluid in Pore of organic matter
Degree;ρbFor density log curve;VsmFor percentage of the volume in inanimate matter volume of inanimate matter solid;VskFor organic matter solid
Percentage of the volume in organic matter volume;DT is acoustic travel time logging curve.
Formula (1)-formula (9) arrangement can obtain formula (10):
It, can be directly true according to door test data when the rock core information of shale reservoir and layer description data abundance
The input parameter of cover half type: ρsm、ρpmfl、ρsk、ρpkfl、DTsm、DTpmfl、DTsk、DTpkfl、φpk, can be obtained according to density log curve
To ρb, DT can be obtained according to acoustic travel time logging curve.ρ can be calculated in above parameters substitution formula (4)pk, substitute into
DT can be calculated in formula (5)pk.Again by DT, DTsm、DTpmfl、ρpk、ρsm、ρb、DTpkAnd ρpmflSubstituting into formula (10) can calculate
To the porosity φ of inanimate matterpm, according to equation φpm+φpk=1 can be calculated the porosity φ of organic matterpk。
A method of component in shale reservoir being calculated using log, formula (1)-formula (9) arrangement can obtain:It will utilize above-mentioned
Using the porosity φ for the inanimate matter that the method that log calculates inanimate matter porosity in shale reservoir is calculatedpm, and
DT、DTsm、DTpkAnd DTpmflParameter, which substitutes into formula (11), can be calculated percentage of the organic matter volume in shale total volume
Vpk.Percentage V of the inanimate matter volume in shale total volume can be calculated further according to formula (1)pm。
Formula (1)-formula (3) arrangement can obtain formula (12): Vsm=Vpm×(1-φpm);Formula (13): Vsk=Vpk×
(1-φpk).By VpmAnd φpmPercentage of the volume of inanimate matter solid in inanimate matter volume can be calculated by substituting into formula (12)
Compare Vsm;By φpkAnd VpkPercentage of the volume of organic matter solid in organic matter volume can be calculated by substituting into formula (13)
Vsk。
A method of total porosity in shale reservoir, the computation model of total porosity are as follows: φ are calculated using logt
=Vpm×φpm+Vpk×φpk, by what is be calculated using the above-mentioned method for calculating component in shale reservoir using log
VpkAnd other parameters are brought into the computation model of total porosity, and the total porosity φ of shale reservoir is calculatedt。
In above-mentioned calculating, V has been obtainedpm、φpm、φpkAnd Vpk, by Vpm、φpm、φpkAnd VpkThis four parameters
It substitutes into the computation model of total porosity, the total porosity φ of shale reservoir can be calculatedt。
In the calculating process of inanimate matter porosity, organic matter porosity, shale reservoir component and total porosity, need
Determining parameter shares 9, respectively ρsm、ρpmfl、ρsk、ρpkfl、DTsm、DTpmfl、DTsk、DTpkfl、φpk, but item indoors
Under part, this nine parameters are difficult to directly determine, it is therefore desirable to take optimization algorithm fitted model parameters.It completes indoors first
The porosity of shale reservoir rock core is tested, and carries out core Location.The porosity value that is obtained again with indoor test is as a comparison
Benchmark takes optimization algorithm to adjust this nine parameters, and the value of this nine parameters is determined after fitting.Using constraint in fit procedure
Optimization algorithm, objective function are as follows:Wherein n is shale samples number;
φtiFor the calculated value of i-th of shale total porosity;φtciFor the test value of i-th of shale total porosity.Work as objective function
Calculated result less than 0.001 when, meet computational accuracy, the adjusted value of parameters is final model parameter, by each ginseng
Several adjusted values brings the computation model of formula (1)-(13) and total porosity, the ρ obtained further according to density log curve intobWith
More accurate inanimate matter porosity, organic matter porosity, shale can be finally calculated in the DT that acoustic travel time logging curve obtains
Reservoir component and total porosity.
Feature and performance of the invention are described in further detail with reference to embodiments.
Embodiment 1
A method of inanimate matter porosity in shale reservoir is calculated using log, comprising:
(1) introduce assumed condition: shale reservoir is made of inanimate matter and organic matter;Containing only gas in the hole of organic matter,
The porosity of organic matter is definite value;Containing only water in the hole of inanimate matter;The grain density of shale reservoir is definite value.
(2) according to the composition of shale reservoir, the volume-based model of shale reservoir is determined.
The volume-based model of shale reservoir are as follows:
Vpm+Vpk=1 formula (1);
Vpm=Vpmφpm+Vsm=Vpmφpm+Vpm(1-φpm) formula (2);
Vpk=Vpkφpk+Vsk=Vpkφpk+Vpk(1-φpk) formula (3);
In formula: VpkFor percentage of the organic matter volume in shale total volume;VpmIt is inanimate matter volume in shale total volume
In percentage;φpmFor the porosity of inanimate matter;φpkFor the porosity of organic matter, VsmFor inanimate matter solid volume in nothing
Percentage in machine plastid product;VskFor percentage of the volume in organic matter volume of organic matter solid.
(3) according to the relationship and shale storage between acoustic logging, density log curve and the composition of shale reservoir
The volume-based model of layer determines the computation model of inanimate matter porosity.
Relationship between acoustic logging, density log curve and shale reservoir are as follows:
ρpk=ρsk×(1-φpk)+ρpkfl×φpkFormula (4);
DTpk=DTsk×(1-φpk)+DTpkfl×φpkFormula (5);
ρpm=ρsm×(1-φpm)+ρpmfl×φpmFormula (6);
DTpm=DTsm×(1-φpm)+DTpmfl×φpmFormula (7);
ρb=ρpmflVpmφpm+ρsmVsm+ρskVsk+ρpkflVpkφpkFormula (8);
DT=DTpmflVpmφpm+DTsmVsm+DTskVsk+DTpkflVpkφpkFormula (9).
In formula, ρpkFor the density of organic matter;ρskFor the density of organic matter solid;ρpkflFor the fluid in the hole of organic matter
Density;DTpkFor the interval transit time of organic matter;DTskFor the interval transit time of organic matter solid;DTpkflTo be flowed in the hole of organic matter
Bulk acoustic wave value;ρpmFor the density of inanimate matter;ρsmFor the density of inanimate matter solid;ρpmflFor the Fluid in Pore density of inanimate matter;
DTpmFor the interval transit time of inanimate matter;DTsmFor the interval transit time of inanimate matter solid;DTpmflIt is close for the Fluid in Pore of organic matter
Degree;ρbFor density log curve;VsmFor percentage of the volume in inanimate matter volume of inanimate matter solid;VskFor organic matter solid
Percentage of the volume in organic matter volume;DT is acoustic travel time logging curve.
Formula (1)-formula (9) arrangement can obtain formula (10):According to room
Interior testing data directly determines the input parameter of model: ρsm、ρpmfl、ρsk、ρpkfl、DTsm、DTpmfl、DTsk、DTpkfl、φpk, according to
ρ can be obtained in density log curveb, DT can be obtained according to acoustic travel time logging curve.It can by above parameters substitution formula (4)
ρ is calculatedpk, DT can be calculated by substituting into formula (5)pk.Again by DT, DTsm、DTpmfl、ρpk、ρsm、ρb、DTpkAnd ρpmflIt substitutes into
The porosity φ of inanimate matter can be calculated in formula (10)pm, according to equation φpm+φpk=1 can be calculated the hole of organic matter
Spend φpk。
A method of component in shale reservoir being calculated using log, formula (1)-formula (9) arrangement can obtain:It will utilize above-mentioned
Using the porosity φ for the inanimate matter that the method that log calculates inanimate matter porosity in shale reservoir is calculatedpm, and
DT、DTsm、DTpkAnd DTpmflParameter, which substitutes into formula (11), can be calculated percentage of the organic matter volume in shale total volume
Vpk.Percentage V of the inanimate matter volume in shale total volume can be calculated further according to formula (1)pm。
Formula (1)-formula (3) arrangement can obtain formula (12): Vsm=Vpm×(1-φpm);Formula (13): Vsk=Vpk×
(1-φpk).By VpmAnd φpmPercentage of the volume of inanimate matter solid in inanimate matter volume can be calculated by substituting into formula (12)
Compare Vsm;By φpkAnd VpkPercentage of the volume of organic matter solid in organic matter volume can be calculated by substituting into formula (13)
Vsk。
A method of total porosity in shale reservoir, the computation model of total porosity are as follows: φ are calculated using logt
=Vpm×φpm+Vpk×φpk, V will have been obtained beforepm、φpm、φpkAnd Vpk, by Vpm、φpm、φpkAnd VpkThis four parameter generations
Enter in the computation model of total porosity, the total porosity φ of shale reservoir is calculatedt。
Embodiment 2
The a kind of of the present embodiment calculates inanimate matter porosity, total porosity and component in shale reservoir using log
Method and embodiment 1 it is a kind of using log calculate shale reservoir in inanimate matter porosity, total porosity and component
Method it is essentially identical, the difference is that, the present embodiment is directed in room conditions, ρsm、ρpmfl、ρsk、ρpkfl、DTsm、
DTpmfl、DTsk、DTpkfl、φpkThis nine parameters are difficult the case where directly determining.
On the basis of embodiment 1, embodiment 2 also adds following steps: as the uncertain ρ of utilizationsm、ρpmfl、ρsk、
ρpkfl、DTsm、DTpmfl、DTsk、DTpkfl、φpkWhen porosity is calculated in parameter, the hole of shale reservoir rock core is first completed indoors
Porosity test, and carry out core Location.The benchmark of the porosity value obtained again with indoor test as a comparison, takes optimization to calculate
Method adjusts ρsm、ρpmfl、ρsk、ρpkfl、DTsm、DTpmfl、DTsk、DTpkfl、φpkThis nine parameters determine this nine parameters after fitting
Value.Constrained optimization algorithm, objective function are used in fit procedure are as follows:Wherein n is shale samples number;φtiFor i-th shale total porosity
Calculated value;φtciFor the test value of i-th of shale total porosity.It is full when the calculated result of objective function is less than 0.001
The adjusted value of sufficient computational accuracy, parameters is final model parameter, brings the adjusted value of parameters into formula (1)-
(13) and the computation model of total porosity, the ρ obtained further according to density log curvebIt is obtained with interval transit time log
More accurate inanimate matter porosity, organic matter porosity, shale reservoir component and total porosity can be finally calculated in DT.
In conclusion a kind of of the embodiment of the present invention calculates inanimate matter porosity in shale reservoir, total using log
The method of porosity and component, the method for calculating inanimate matter porosity in shale reservoir using log, passes through volume
Model and page can be calculated using the relationship between acoustic logging, density log curve and the composition of shale reservoir
Inanimate matter porosity in rock reservoir.Calculated inanimate matter porosity and utilization acoustic logging, density log
Relationship between curve and the composition of shale reservoir can calculate each component accounting in shale reservoir.It is calculated
Each component accounting and inanimate matter porosity in shale reservoir can calculate the total porosity in shale reservoir.
Embodiments described above is a part of the embodiment of the present invention, instead of all the embodiments.Reality of the invention
The detailed description for applying example is not intended to limit the range of claimed invention, but is merely representative of selected implementation of the invention
Example.Based on the embodiments of the present invention, obtained by those of ordinary skill in the art without making creative efforts
Every other embodiment, shall fall within the protection scope of the present invention.
Claims (7)
1. a kind of method for calculating inanimate matter porosity in shale reservoir using log characterized by comprising
Introduce assumed condition: the shale reservoir is made of inanimate matter and organic matter;Containing only gas in the hole of the organic matter,
The porosity of the organic matter is definite value;Containing only water in the hole of the inanimate matter;The grain density of the shale reservoir is fixed
Value;
According to the composition of the shale reservoir, the volume-based model of the shale reservoir is determined;And
According to the relationship and shale storage between acoustic logging, density log curve and the composition of the shale reservoir
The volume-based model of layer determines the computation model of the inanimate matter porosity;
The volume-based model of the shale reservoir are as follows:
Vpm+Vpk=1 formula (1);
Vpm=Vpmφpm+Vsm=Vpmφpm+Vpm(1-φpm) formula (2);
Vpk=Vpkφpk+Vsk=Vpkφpk+Vpk(1-φpk) formula (3);
In formula: VpkFor percentage of the organic matter volume in shale total volume;VpmIt is inanimate matter volume in shale total volume
Percentage;φpmFor the porosity of inanimate matter;φpkFor the porosity of organic matter, VsmFor inanimate matter solid volume in inanimate matter
Percentage in volume;VskFor percentage of the volume in organic matter volume of organic matter solid;
Relationship between the acoustic logging, the density log curve and the composition of the shale reservoir are as follows:
ρpk=ρsk×(1-φpk)+ρpkfl×φpkFormula (4);
DTpk=DTsk×(1-φpk)+DTpkfl×φpkFormula (5);
ρpm=ρsm×(1-φpm)+ρpmfl×φpmFormula (6);
DTpm=DTsm×(1-φpm)+DTpmfl×φpmFormula (7);
ρb=ρpmflVpmφpm+ρsmVsm+ρskVsk+ρpkflVpkφpkFormula (8);
DT=DTpmflVpmφpm+DTsmVsm+DTskVsk+DTpkflVpkφpkFormula (9);
In formula, ρpkFor the density of organic matter;ρskFor the density of organic matter solid;ρpkflIt is close for the fluid in the hole of organic matter
Degree;DTpkFor the interval transit time of organic matter;DTskFor the interval transit time of organic matter solid;DTpkflFor the Fluid in Pore of organic matter
Sound wave value;ρpmFor the density of inanimate matter;ρsmFor the density of inanimate matter solid;ρpmflFor the Fluid in Pore density of inanimate matter;
DTpmFor the interval transit time of inanimate matter;DTsmFor the interval transit time of inanimate matter solid;DTpmflIt is close for the Fluid in Pore of organic matter
Degree;ρbFor density log curve;VsmFor percentage of the volume in inanimate matter volume of inanimate matter solid;VskFor organic matter solid
Percentage of the volume in organic matter volume;DT is acoustic travel time logging curve;
The formula (the 1)-formula (9) arrangement can obtain:
It brings parameters into the formula (10), the porosity φ of inanimate matter is calculatedpm。
2. the method according to claim 1 for calculating inanimate matter porosity in shale reservoir using log, feature
It is,
It include: the porosity φ for the inanimate matter that will be calculatedpmAnd the formula (the 1)-formula (9) arrangement can
:
By the porosity φ of the inanimate matterpmThe formula (11) is substituted into other parameters, the organic matter volume is calculated
Percentage V in shale total volumepk。
3. the method according to claim 2 for calculating inanimate matter porosity in shale reservoir using log, feature
It is, the formula (the 1)-formula (3) arrangement can obtain;
Vsm=Vpm×(1-φpm) formula (12);
Vsk=Vpk×(1-φpk) formula (13);
In formula: parameters being brought into the formula (12) and the formula (13) respectively, it is solid to can be calculated the inanimate matter
Percentage V of the volume of body in inanimate matter volumesm, percentage of the volume of the organic matter solid in organic matter volume
Vsk。
4. the method according to claim 3 for calculating inanimate matter porosity in shale reservoir using log, feature
It is, the computation model of total porosity are as follows: φt=Vpm×φpm+Vpk×φpk, the V that will be calculatedpkAnd other parameters are brought into
In the computation model of the total porosity, the total porosity φ of the shale reservoir is calculatedt。
5. the method according to claim 4 for calculating inanimate matter porosity in shale reservoir using log, feature
Be, further include the steps that optimization algorithm fitted model parameters: by under laboratory conditions to the rock core of shale reservoir into
The benchmark of the total pore space angle value that row test obtains as a comparison, adjusts parameters using optimization algorithm, determines after fitting each
The value of a parameter finally determines final total porosity again.
6. the method according to claim 5 for calculating inanimate matter porosity in shale reservoir using log, feature
It is, the objective function of the optimization algorithm are as follows:
Wherein n is shale samples number;φtiFor i-th of shale total porosity
Calculated value;φtciFor the test value of i-th of shale total porosity.
7. the method according to claim 6 for calculating inanimate matter porosity in shale reservoir using log, feature
It is, when the calculated result of the objective function is less than 0.001, meets computational accuracy, the adjusted value of the parameters is
Final model parameter brings the adjusted value of the parameters into the formula (1)-formula (13) and total hole
The computation model of porosity can calculate the total porosity of shale reservoir using acoustic logging and density log curve.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710421951.XA CN107194104B (en) | 2017-06-07 | 2017-06-07 | The method for calculating inanimate matter porosity in shale reservoir |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710421951.XA CN107194104B (en) | 2017-06-07 | 2017-06-07 | The method for calculating inanimate matter porosity in shale reservoir |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107194104A CN107194104A (en) | 2017-09-22 |
CN107194104B true CN107194104B (en) | 2019-08-09 |
Family
ID=59876486
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710421951.XA Active CN107194104B (en) | 2017-06-07 | 2017-06-07 | The method for calculating inanimate matter porosity in shale reservoir |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107194104B (en) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108412488B (en) * | 2018-03-13 | 2021-09-10 | 中石化石油工程技术服务有限公司 | Logging method for rapidly determining organic porosity of shale gas reservoir |
CN110223282B (en) * | 2019-06-10 | 2021-05-14 | 中国石油大学(华东) | Automatic identification method and system for organic pores and inorganic pores of shale |
CN111122408A (en) * | 2019-12-23 | 2020-05-08 | 中国地质大学(武汉) | Method for estimating organic porosity of shale with pore diameter of 0.3-10nm |
CN111079313B (en) * | 2019-12-31 | 2021-12-17 | 成都理工大学 | Carbonate reservoir information classification processing method and information data processing terminal |
CN112304837B (en) * | 2020-06-24 | 2021-06-25 | 成都理工大学 | Method for judging shale oil reservoir organic matter microscopic occurrence structure |
CN113484220B (en) * | 2021-07-26 | 2022-07-22 | 中国石油大学(北京) | Method and device for determining organic matter and inorganic matter distribution of shale and electronic equipment |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104573344A (en) * | 2014-12-25 | 2015-04-29 | 中国海洋石油总公司 | Method for acquiring gas content of shale reservoir through well logging data |
CN105043957A (en) * | 2015-07-06 | 2015-11-11 | 成都理工大学 | Method for determining pore structure through classifying isothermal adsorption curve of mud shale |
CN105696998A (en) * | 2016-04-05 | 2016-06-22 | 成都理工大学 | Shale gas reservoir fracturing modification method |
-
2017
- 2017-06-07 CN CN201710421951.XA patent/CN107194104B/en active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104573344A (en) * | 2014-12-25 | 2015-04-29 | 中国海洋石油总公司 | Method for acquiring gas content of shale reservoir through well logging data |
CN105043957A (en) * | 2015-07-06 | 2015-11-11 | 成都理工大学 | Method for determining pore structure through classifying isothermal adsorption curve of mud shale |
CN105696998A (en) * | 2016-04-05 | 2016-06-22 | 成都理工大学 | Shale gas reservoir fracturing modification method |
Non-Patent Citations (2)
Title |
---|
砂岩-页岩互层气藏物质平衡方程构建与应用;杨宇 等;《煤炭学报》;20160131;第41卷(第1期);第174-180页 * |
页岩油气储层测井评价技术及应用;李霞 等;《天然气地球科学》;20150531;第26卷(第5期);第904-914页 * |
Also Published As
Publication number | Publication date |
---|---|
CN107194104A (en) | 2017-09-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107194104B (en) | The method for calculating inanimate matter porosity in shale reservoir | |
CN105986815B (en) | Method for identifying geological dessert in shale stratum | |
Winton | Do climate models underestimate the sensitivity of Northern Hemisphere sea ice cover? | |
CN105095631B (en) | A kind of shale anisotropic rock Method of Physical Modeling | |
CN105913494B (en) | The modeling of multi-scale facture fine geology and method for numerical simulation and device | |
Niwa et al. | A three-dimensional icosahedral grid advection scheme preserving monotonicity and consistency with continuity for atmospheric tracer transport | |
CN105089615B (en) | A kind of log data historical regression processing method based on reservoir model | |
CN112133377B (en) | Method and system for judging occurrence state of natural gas hydrate in sea area | |
CN111123354B (en) | Method and equipment for predicting dense gas layer based on frequency-dependent reflection amplitude attenuation | |
CN105891089B (en) | Reservoir permeability determines method and device | |
CN107329171A (en) | Depth Domain reservoir seismic inversion method and device | |
CN105425280A (en) | Prediction method for mineral modulus and pore structure | |
CN107609253B (en) | Carbonate rock deposition numerical simulation method | |
CN105863628B (en) | A kind of phase of development subtle hydrocarbon reservoir method | |
CN112946783B (en) | Hydrate saturation determination method, device and equipment | |
CN109829217A (en) | Pressure break Fractured Reservoir productivity simulation method and device | |
CN110703330B (en) | Carbonate rock logging correction and transverse wave velocity prediction method based on rock physics | |
CN112796738A (en) | Stratum permeability calculation method combining array acoustic logging and conventional logging | |
CN111208565B (en) | KT model-based hole seam parameter inversion method and device and storage medium | |
CN105317429A (en) | Method for determining content of shale stratum kerogen | |
CN108957554B (en) | Seismic inversion method in geophysical exploration | |
CN110320575A (en) | Method and device is determined based on the shale content of organic matter of petrophysical model | |
CN105301657B (en) | A kind of curve correcting method based on rock physicses meaning | |
CN117055114B (en) | Quantitative analysis method for free gas saturation of reservoir sediment | |
CN113189674A (en) | Method and system for estimating saturation of natural gas hydrate |
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 |