CN104122182A - Method for obtaining effective thickness lower limit of mine reservoir - Google Patents

Method for obtaining effective thickness lower limit of mine reservoir Download PDF

Info

Publication number
CN104122182A
CN104122182A CN201410307584.7A CN201410307584A CN104122182A CN 104122182 A CN104122182 A CN 104122182A CN 201410307584 A CN201410307584 A CN 201410307584A CN 104122182 A CN104122182 A CN 104122182A
Authority
CN
China
Prior art keywords
power function
relationship formula
function relationship
curvature
core
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
Application number
CN201410307584.7A
Other languages
Chinese (zh)
Other versions
CN104122182B (en
Inventor
谢玉洪
周家雄
何胜林
胡向阳
张恒荣
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
WUHAN TIMES GEOSMART TECHNOLOGY Co Ltd
China National Offshore Oil Corp CNOOC
CNOOC China Ltd Zhanjiang Branch
Original Assignee
WUHAN TIMES GEOSMART TECHNOLOGY Co Ltd
China National Offshore Oil Corp CNOOC
CNOOC China Ltd Zhanjiang Branch
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by WUHAN TIMES GEOSMART TECHNOLOGY Co Ltd, China National Offshore Oil Corp CNOOC, CNOOC China Ltd Zhanjiang Branch filed Critical WUHAN TIMES GEOSMART TECHNOLOGY Co Ltd
Priority to CN201410307584.7A priority Critical patent/CN104122182B/en
Publication of CN104122182A publication Critical patent/CN104122182A/en
Application granted granted Critical
Publication of CN104122182B publication Critical patent/CN104122182B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Electroluminescent Light Sources (AREA)

Abstract

The invention discloses a method for obtaining effective thickness lower limit of a mine reservoir. According to the basic principle, the method comprises the following steps: taking a core from a reservoir to be tested, and respectively performing laboratory test and analysis on the core, so as to obtain porosity data and permeability data of the core; injecting mercury into the core to obtain capillary pressure corresponding to the condition that the injected mercury saturation is 50%, namely mercury saturation medium-value pressure; fitting the mercury saturation medium-value pressure with the core porosity and the core permeability, respectively, so as to obtain a regression relation; and solving a maximum curvature value of a power function curve by use of a first-order derivative, a second-order derivative and the like, so that the corresponding porosity and permeability are the lower limit values. The method is used for calculating the lower limit values of the reservoir physical property by use of mercury saturation medium-value pressure parameters which are capable of well reflecting a pore throat structure of rock, thereby covering the shortage that a conventional method is capable of obtaining the lower limits of the reservoir physical property only by use of an empirical statistical method, and therefore, the accuracy of the obtained lower limits of the reservoir physical property is improved, and meanwhile, the efficiency of obtaining the lower limits of the reservoir physical property is improved.

Description

The acquisition methods of mine reservoir net thickness lower limit
Technical field
The present invention relates to Geophysical Logging field, refer to particularly a kind of acquisition methods of mine reservoir net thickness lower limit.
Background technology
In oilfield explorating developing process, need logging technology personnel carry out timely well logging interpretation and mark off Effective Reservoirs (oil reservoir, gas-bearing formation) the underground logging data obtaining, the key of work is to determine the standard (also determining porosity and permeability lower limit) of dividing Effective Reservoirs.The acquisition methods of existing mine reservoir net thickness lower limit is oil-bearing occurrence method, in the detailed record of the method and the company standard Q/HXJ3026-93 of Nanhai Western Co., Chinese Marine Petroleum, but technician finds that the accuracy of the mine reservoir net thickness lower limit that said method obtains is lower in long production practices, tracing it to its cause is mainly that the subjectivity of the method is larger, has used the process of a lot of estimations.In addition, in the method, an important measurement parameter is oiliness, for this method in gas field inapplicable.
Summary of the invention
Object of the present invention is exactly the acquisition methods that a kind of mine reservoir net thickness lower limit will be provided, the accuracy of the mine reservoir net thickness lower limit that utilizes the method to significantly improve to obtain, and the method is applicable equally to gas field, has stronger versatility.
For realizing this object, the acquisition methods of the designed mine reservoir net thickness lower limit of the present invention, is characterized in that, it comprises the steps:
Step 1: coring in stratum to be measured, this rock core is carried out respectively laboratory assay analysis and is obtained core porosity data and the core permeability data of this rock core, then inject mercury to described rock core, corresponding capillary pressure when obtaining injecting mercury saturation and reaching 50%, i.e. duty pressure in mercury saturation;
Step 2: duty pressure in core porosity data and mercury saturation is set up to the first power function relationship formula y 1=a 1x 1 b1, wherein, x 1for core porosity data, y 1for duty pressure in mercury saturation, a 1and b 1be power function constant; Duty pressure in core permeability data and mercury saturation is set up to the second power function relationship formula y 2=a 2x 2 b2, wherein, x 2for core permeability data, y 2for duty pressure in mercury saturation, a 2and b 2be power function constant;
Step 3: the method by nonlinear fitting is carried out power function fitting to the first power function relationship formula and the second power function relationship formula and asked for power function constant a 1, b 1, a 2and b 2;
Step 4: in the first power function relationship formula, curvature maximum represents mercury saturation intermediate value pressure change rate maximum, is power function curve catastrophe point herein, and X-axis numerical value corresponding to this power function curve catastrophe point is core porosity lower limit;
In the second power function relationship formula, curvature maximum represents mercury saturation intermediate value pressure change rate maximum, is power function curve catastrophe point herein, and X-axis numerical value corresponding to this power function curve catastrophe point is core permeability lower limit;
Step 5: by the first power function relationship formula y 1=a 1x 1 b1the curvature equation of substitution the first power function relationship formula;
k 1 = | y 1 ′ ′ | ( 1 + y 1 ′ 2 ) 3 2
Wherein, y 1' be that the single order of the first power function relationship formula is led, y 1" second order that is the first power function relationship formula is led;
Obtain the curvature of the first power function relationship formula
k 1 = | a 1 b 1 ( b 1 - 1 ) x 1 ( b 1 - 2 ) | ( 1 + a 1 2 b 1 2 x 1 ( 2 b 1 - 2 ) ) 3 2
By the second power function relationship formula y 2=a 2x 2 b2the curvature equation of substitution the second power function relationship formula;
k 2 = | y 2 ′ ′ | ( 1 + y 2 ′ 2 ) 3 2
Wherein, y 2' be that the single order of the second power function relationship formula is led, y 2" second order that is the second power function relationship formula is led;
Obtain the curvature of the second power function relationship formula:
k 2 = | a 2 b 2 ( b 2 - 1 ) x 2 ( b 2 - 2 ) | ( 1 + a 2 2 b 2 2 x 2 ( 2 b 2 - 2 ) ) 3 2
Step 6: calculate the maximal value of the curvature of the first power function relationship formula, the curvature of the first power function relationship formula is asked to first order derivative according to following formula:
k 1 ′ = a 1 b 1 ( b 1 - 1 ) ( b 1 - 2 ) x 1 ( b 1 - 3 ) ( 1 + a 1 2 b 1 2 x 1 ( 2 b 1 - 2 ) ) 3 2 - a 1 b 1 ( b 1 - 1 ) x 1 ( b 1 - 2 ) 3 2 ( 1 + a 1 2 b 1 2 x 1 ( 2 b 1 - 2 ) ) 1 2 a 1 2 b 1 2 2 ( b 1 - 1 ) x 1 ( 2 b 1 - 3 ) ( 1 + a 1 2 b 1 2 x 1 ( 2 b 1 - 2 ) ) 3
First order derivative in the curvature of the first power function relationship formula is 1 o'clock, the x obtaining 1value be the maximal value of the curvature of the first power function relationship formula;
Calculate the maximal value of the curvature of the second power function relationship formula, the curvature of the second power function relationship formula asked to first order derivative according to following formula:
k 2 ′ = a 2 b 2 ( b 2 - 1 ) ( b 2 - 2 ) x 2 ( b 2 - 3 ) ( 1 + a 2 2 b 2 2 x 2 ( 2 b 2 - 2 ) ) 3 2 - a 2 b 2 ( b 2 - 1 ) x 2 ( b 2 - 2 ) 3 2 ( 1 + a 2 2 b 2 2 x 2 ( 2 b 2 - 2 ) ) 1 2 a 2 2 b 2 2 2 ( b 2 - 1 ) x 2 ( 2 b 2 - 3 ) ( 1 + a 2 2 b 2 2 x 2 ( 2 b 2 - 2 ) ) 3
First order derivative in the curvature of the second power function relationship formula is 1 o'clock, the x obtaining 2value be the maximal value of the curvature of the second power function relationship formula;
Step S7: the relational expression of the curvature first order derivative of above-mentioned the first power function relationship formula is reduced to:
k 1 ′ = a 1 b 1 ( b 1 - 1 ) x 1 ( b 1 - 3 ) [ ( b 1 - 2 ) ( 1 + a 1 2 b 1 2 x 1 ( 2 b 1 - 2 ) ) 3 2 - 3 a 1 2 b 1 2 ( b 1 - 1 ) x 1 ( 2 b 1 - 2 ) ( 1 + a 1 2 b 1 2 x 1 ( 2 b 1 - 2 ) ) 1 2 ] ( 1 + a 1 2 b 1 2 x 1 ( 2 b 1 - 2 ) ) 3
The relational expression of the curvature first order derivative of above-mentioned the second power function relationship formula is reduced to:
k 2 ′ = a 2 b 2 ( b 2 - 1 ) x 2 ( b 2 - 3 ) [ ( b 2 - 2 ) ( 1 + a 2 2 b 2 2 x 2 ( 2 b 2 - 2 ) ) 3 2 - 3 a 2 2 b 2 2 ( b 2 - 1 ) x 2 ( 2 b 2 - 2 ) ( 1 + a 2 2 b 2 2 x 2 ( 2 b 2 - 2 ) ) 1 2 ] ( 1 + a 2 2 b 2 2 x 2 ( 2 b 2 - 2 ) ) 3
Step S8: for the relational expression of the curvature first order derivative of the first power function relationship formula, due to a in the first power function relationship formula 1and b 1all non-vanishing, if a 1and b 1have one to be zero, the first power function relationship formula becomes y 1=0, now the curvature of the first power function relationship formula does not have maximum point; And b 1be not 1, if b 1be that 1 the first power function relationship formula becomes y 1=a 1x 1, there is no point of maximum curvature; Meanwhile, x 1be greater than zero, due to x 1represent core porosity data, so x 1must be greater than zero, therefore the first order derivative of the curvature of the first power function relationship formula is that the condition of 1 o'clock is:
( b 1 - 2 ) + ( b 1 - 2 ) a 1 2 b 1 2 x 1 ( 2 b 1 - 2 ) - 3 ( b 1 - 1 ) a 1 2 b 1 2 x 1 ( 2 b 1 - 2 ) = 0
Have:
x 1 ( 2 b 1 - 2 ) = b 1 - 2 ( 2 b 1 - 1 ) a 1 2 b 1 2
By a 1and b 1substitution above formula is also obtained x 1can obtain core porosity data lower limit;
For the relational expression of the curvature first order derivative of the second power function relationship formula, due to a in the second power function relationship formula 2and b 2all non-vanishing, if a 2and b 2have one to be zero, the second power function relationship formula becomes y 2=0, now the curvature of the second power function relationship formula does not have maximum point; And b 2be not 1, if b 2be that 1 the second power function relationship formula becomes y 2=a 2x 2, there is no point of maximum curvature; Meanwhile, x 2be greater than zero, due to x 2represent core permeability data, so x 2must be greater than zero, therefore the first order derivative of the curvature of the second power function relationship formula is that the condition of 1 o'clock is:
( b 2 - 2 ) + ( b 2 - 2 ) a 2 2 b 2 2 x 2 ( 2 b 2 - 2 ) - 3 ( b 2 - 1 ) a 2 2 b 2 2 x 2 ( 2 b 2 - 2 ) = 0
Have:
x 2 ( 2 b 2 - 2 ) = b 2 - 2 ( 2 b 2 - 1 ) a 2 2 b 2 2
By a 2and b 2substitution above formula is also obtained x 2can obtain core permeability data lower limit.
Beneficial effect of the present invention:
The present invention can better react pore throat structure (rock interior pore constriction structure) mercury saturation intermediate value pressure parameter by adopting calculates reservoir properties (core porosity data and core permeability data) lower limit, make up the deficiency that conventional method can only obtain by empirical statistics method reservoir properties lower limit, improve the accuracy of the reservoir properties lower limit obtaining, improved the efficiency of reservoir properties lower limit simultaneously.The present invention assesses better more stable the lower limit of reservoirs acquisition methods is provided for reservoir properties lower limit in oilfield prospecting developing.
Brief description of the drawings
Fig. 1 is duty pressure and core porosity data fitting relational graph in mercury saturation in the present invention;
Embodiment
Below in conjunction with the drawings and specific embodiments, the present invention is described in further detail:
An acquisition methods for mine reservoir net thickness lower limit, it comprises the steps:
Step 1: coring in stratum to be measured, this rock core is carried out respectively laboratory assay analysis and is obtained core porosity data and the core permeability data of this rock core, then inject mercury to described rock core, corresponding capillary pressure when obtaining injecting mercury saturation and reaching 50%, i.e. duty pressure in mercury saturation;
Step 2: duty pressure in core porosity data and mercury saturation is set up to the first power function relationship formula y 1=a 1x 1 b1, wherein, x 1for core porosity data, y 1for duty pressure in mercury saturation, a 1and b 1be power function constant; Duty pressure in core permeability data and mercury saturation is set up to the second power function relationship formula y 2=a 2x 2 b2, wherein, x 2for core permeability data, y 2for duty pressure in mercury saturation, a 2and b 2be power function constant;
Step 3: the method by nonlinear fitting is carried out power function fitting to the first power function relationship formula and the second power function relationship formula and asked for power function constant a 1, b 1, a 2and b 2;
Step 4: in the first power function relationship formula, curvature maximum represents mercury saturation intermediate value pressure change rate maximum, is power function curve catastrophe point herein, and X-axis numerical value corresponding to this power function curve catastrophe point is core porosity lower limit;
In the second power function relationship formula, curvature maximum represents mercury saturation intermediate value pressure change rate maximum, is power function curve catastrophe point herein, and X-axis numerical value corresponding to this power function curve catastrophe point is core permeability lower limit;
Step 5: by the first power function relationship formula y 1=a 1x 1 b1the curvature equation of substitution the first power function relationship formula;
k 1 = | y 1 ′ ′ | ( 1 + y 1 ′ 2 ) 3 2
Wherein, y 1' be that the single order of the first power function relationship formula is led, y 1" second order that is the first power function relationship formula is led;
Obtain the curvature of the first power function relationship formula
k 1 = | a 1 b 1 ( b 1 - 1 ) x 1 ( b 1 - 2 ) | ( 1 + a 1 2 b 1 2 x 1 ( 2 b 1 - 2 ) ) 3 2
By the second power function relationship formula y 2=a 2x 2 b2the curvature equation of substitution the second power function relationship formula;
k 2 = | y 2 ′ ′ | ( 1 + y 2 ′ 2 ) 3 2
Wherein, y 2' be that the single order of the second power function relationship formula is led, y 2" second order that is the second power function relationship formula is led;
Obtain the curvature of the second power function relationship formula:
k 2 = | a 2 b 2 ( b 2 - 1 ) x 2 ( b 2 - 2 ) | ( 1 + a 2 2 b 2 2 x 2 ( 2 b 2 - 2 ) ) 3 2
Step 6: calculate the maximal value of the curvature of the first power function relationship formula, the curvature of the first power function relationship formula is asked to first order derivative according to following formula:
k 1 ′ = a 1 b 1 ( b 1 - 1 ) ( b 1 - 2 ) x 1 ( b 1 - 3 ) ( 1 + a 1 2 b 1 2 x 1 ( 2 b 1 - 2 ) ) 3 2 - a 1 b 1 ( b 1 - 1 ) x 1 ( b 1 - 2 ) 3 2 ( 1 + a 1 2 b 1 2 x 1 ( 2 b 1 - 2 ) ) 1 2 a 1 2 b 1 2 2 ( b 1 - 1 ) x 1 ( 2 b 1 - 3 ) ( 1 + a 1 2 b 1 2 x 1 ( 2 b 1 - 2 ) ) 3
First order derivative in the curvature of the first power function relationship formula is 1 o'clock, the x obtaining 1value be the maximal value of the curvature of the first power function relationship formula;
Calculate the maximal value of the curvature of the second power function relationship formula, the curvature of the second power function relationship formula asked to first order derivative according to following formula:
k 2 ′ = a 2 b 2 ( b 2 - 1 ) ( b 2 - 2 ) x 2 ( b 2 - 3 ) ( 1 + a 2 2 b 2 2 x 2 ( 2 b 2 - 2 ) ) 3 2 - a 2 b 2 ( b 2 - 1 ) x 2 ( b 2 - 2 ) 3 2 ( 1 + a 2 2 b 2 2 x 2 ( 2 b 2 - 2 ) ) 1 2 a 2 2 b 2 2 2 ( b 2 - 1 ) x 2 ( 2 b 2 - 3 ) ( 1 + a 2 2 b 2 2 x 2 ( 2 b 2 - 2 ) ) 3
First order derivative in the curvature of the second power function relationship formula is 1 o'clock, the x obtaining 2value be the maximal value of the curvature of the second power function relationship formula;
Step S7: the relational expression of the curvature first order derivative of above-mentioned the first power function relationship formula is reduced to:
k 1 ′ = a 1 b 1 ( b 1 - 1 ) x 1 ( b 1 - 3 ) [ ( b 1 - 2 ) ( 1 + a 1 2 b 1 2 x 1 ( 2 b 1 - 2 ) ) 3 2 - 3 a 1 2 b 1 2 ( b 1 - 1 ) x 1 ( 2 b 1 - 2 ) ( 1 + a 1 2 b 1 2 x 1 ( 2 b 1 - 2 ) ) 1 2 ] ( 1 + a 1 2 b 1 2 x 1 ( 2 b 1 - 2 ) ) 3
The relational expression of the curvature first order derivative of above-mentioned the second power function relationship formula is reduced to:
Step S8: for the relational expression of the curvature first order derivative of the first power function relationship formula, due to a in the first power function relationship formula 1and b 1all non-vanishing, if a 1and b 1there is one to be zero,
k 2 ′ = a 2 b 2 ( b 2 - 1 ) x 2 ( b 2 - 3 ) [ ( b 2 - 2 ) ( 1 + a 2 2 b 2 2 x 2 ( 2 b 2 - 2 ) ) 3 2 - 3 a 2 2 b 2 2 ( b 2 - 1 ) x 2 ( 2 b 2 - 2 ) ( 1 + a 2 2 b 2 2 x 2 ( 2 b 2 - 2 ) ) 1 2 ] ( 1 + a 2 2 b 2 2 x 2 ( 2 b 2 - 2 ) ) 3
The first power function relationship formula becomes y 1=0, now the curvature of the first power function relationship formula does not have maximum point; And b 1be not 1, if b 1be that 1 the first power function relationship formula becomes y 1=a 1x 1, there is no point of maximum curvature; Meanwhile, x 1be greater than zero, due to x 1represent core porosity data, so x 1must be greater than zero, therefore the first order derivative of the curvature of the first power function relationship formula is 1 o'clock
( b 1 - 2 ) + ( b 1 - 2 ) a 1 2 b 1 2 x 1 ( 2 b 1 - 2 ) - 3 ( b 1 - 1 ) a 1 2 b 1 2 x 1 ( 2 b 1 - 2 ) = 0 Condition be:
Have:
x 1 ( 2 b 1 - 2 ) = b 1 - 2 ( 2 b 1 - 1 ) a 1 2 b 1 2
By a 1and b 1substitution above formula is also obtained x 1can obtain core porosity data lower limit;
For the relational expression of the curvature first order derivative of the second power function relationship formula, due to the second power
( b 2 - 2 ) + ( b 2 - 2 ) a 2 2 b 2 2 x 2 ( 2 b 2 - 2 ) - 3 ( b 2 - 1 ) a 2 2 b 2 2 x 2 ( 2 b 2 - 2 ) = 0
A in functional relation 2and b 2all non-vanishing, if a 2and b 2have one to be zero, the second power function relationship formula becomes y 2=0, now the curvature of the second power function relationship formula does not have maximum point; And b 2be not 1, if b 2be that 1 the second power function relationship formula becomes y 2=a 2x 2, there is no point of maximum curvature; Meanwhile, x 2be greater than zero, due to x 2represent core permeability data, so x 2must be greater than zero, therefore the first order derivative of the curvature of the second power function relationship formula is that the condition of 1 o'clock is:
Have:
x 2 ( 2 b 2 - 2 ) = b 2 - 2 ( 2 b 2 - 1 ) a 2 2 b 2 2
By a 2and b 2substitution above formula is also obtained x 2can obtain core permeability data lower limit.
Ultimate principle of the present invention is coring in stratum to be measured, this rock core is carried out respectively laboratory assay analysis and is obtained core porosity data and the core permeability data of this rock core, then inject mercury to described rock core, corresponding capillary pressure when obtaining injecting mercury saturation and reaching 50%, it is duty pressure in mercury saturation, duty pressure in mercury saturation and core porosity, core permeability are distinguished to matching, obtain the regression relation (power function equation formula) of the two.Ask for the curvature maximal value of this power function curve by the mode such as first order derivative, second derivative, corresponding factor of porosity, permeability is lower limit (curvature maximum shows that curve changes maximum herein, utilizes this feature as boundary value) herein.
In the first well and the second well, sample respectively the measurement of duty pressure in carrying out core porosity and mercury saturation, the core porosity of the first well and the second well and mercury saturation intermediate value force value are drawn, obtain two graphs of a relation, then data point is carried out to matching, obtain duty pressure and core porosity data fitting relational graph in mercury saturation as described in Figure 1, in Fig. 1, at the first power function relationship formula y 1=1E+09x 1 -7.507the curvature maximum of (expression way that E+09 is scientific notation, 9 powers that its value equals 10) has corresponding core porosity lower limit.In figure, R 2represent duty pressure and core porosity data in mercury saturation related coefficient square, in the larger expression mercury saturation of this value, the correlativity of duty pressure and core porosity data is better, be that the two correlativity is better, return the formula obtaining more reliable, the lower limit of calculating is more accurate.
The content that this instructions is not described in detail belongs to the known prior art of professional and technical personnel in the field.

Claims (1)

1. an acquisition methods for mine reservoir net thickness lower limit, is characterized in that, it comprises the steps:
Step 1: coring in stratum to be measured, this rock core is carried out respectively laboratory assay analysis and is obtained core porosity data and the core permeability data of this rock core, then inject mercury to described rock core, corresponding capillary pressure when obtaining injecting mercury saturation and reaching 50%, i.e. duty pressure in mercury saturation;
Step 2: duty pressure in core porosity data and mercury saturation is set up to the first power function relationship formula y 1=a 1x 1 b1, wherein, x 1for core porosity data, y 1for duty pressure in mercury saturation, a 1and b 1be power function constant; Duty pressure in core permeability data and mercury saturation is set up to the second power function relationship formula y 2=a 2x 2 b2, wherein, x 2for core permeability data, y 2for duty pressure in mercury saturation, a 2and b 2be power function constant;
Step 3: the method by nonlinear fitting is carried out power function fitting to the first power function relationship formula and the second power function relationship formula and asked for power function constant a 1, b 1, a 2and b 2;
Step 4: in the first power function relationship formula, curvature maximum represents mercury saturation intermediate value pressure change rate maximum, is power function curve catastrophe point herein, and X-axis numerical value corresponding to this power function curve catastrophe point is core porosity lower limit;
In the second power function relationship formula, curvature maximum represents mercury saturation intermediate value pressure change rate maximum, is power function curve catastrophe point herein, and X-axis numerical value corresponding to this power function curve catastrophe point is core permeability lower limit;
Step 5: by the first power function relationship formula y 1=a 1x 1 b1the curvature equation of substitution the first power function relationship formula;
k 1 = | y 1 ′ ′ | ( 1 + y 1 ′ 2 ) 3 2
Wherein, y 1' be that the single order of the first power function relationship formula is led, y 1" second order that is the first power function relationship formula is led;
Obtain the curvature of the first power function relationship formula
k 1 = | a 1 b 1 ( b 1 - 1 ) x 1 ( b 1 - 2 ) | ( 1 + a 1 2 b 1 2 x 1 ( 2 b 1 - 2 ) ) 3 2
By the second power function relationship formula y 2=a 2x 2 b2the curvature equation of substitution the second power function relationship formula;
k 2 = | y 2 ′ ′ | ( 1 + y 2 ′ 2 ) 3 2
Wherein, y 2' be that the single order of the second power function relationship formula is led, y 2" second order that is the second power function relationship formula is led;
Obtain the curvature of the second power function relationship formula:
k 2 = | a 2 b 2 ( b 2 - 1 ) x 2 ( b 2 - 2 ) | ( 1 + a 2 2 b 2 2 x 2 ( 2 b 2 - 2 ) ) 3 2
Step 6: calculate the maximal value of the curvature of the first power function relationship formula, the curvature of the first power function relationship formula is asked to first order derivative according to following formula:
k 1 ′ = a 1 b 1 ( b 1 - 1 ) ( b 1 - 2 ) x 1 ( b 1 - 3 ) ( 1 + a 1 2 b 1 2 x 1 ( 2 b 1 - 2 ) ) 3 2 - a 1 b 1 ( b 1 - 1 ) x 1 ( b 1 - 2 ) 3 2 ( 1 + a 1 2 b 1 2 x 1 ( 2 b 1 - 2 ) ) 1 2 a 1 2 b 1 2 2 ( b 1 - 1 ) x 1 ( 2 b 1 - 3 ) ( 1 + a 1 2 b 1 2 x 1 ( 2 b 1 - 2 ) ) 3
First order derivative in the curvature of the first power function relationship formula is 1 o'clock, the x obtaining 1value be the maximal value of the curvature of the first power function relationship formula;
Calculate the maximal value of the curvature of the second power function relationship formula, the curvature of the second power function relationship formula asked to first order derivative according to following formula:
k 2 ′ = a 2 b 2 ( b 2 - 1 ) ( b 2 - 2 ) x 2 ( b 2 - 3 ) ( 1 + a 2 2 b 2 2 x 2 ( 2 b 2 - 2 ) ) 3 2 - a 2 b 2 ( b 2 - 1 ) x 2 ( b 2 - 2 ) 3 2 ( 1 + a 2 2 b 2 2 x 2 ( 2 b 2 - 2 ) ) 1 2 a 2 2 b 2 2 2 ( b 2 - 1 ) x 2 ( 2 b 2 - 3 ) ( 1 + a 2 2 b 2 2 x 2 ( 2 b 2 - 2 ) ) 3
First order derivative in the curvature of the second power function relationship formula is 1 o'clock, the x obtaining 2value be the maximal value of the curvature of the second power function relationship formula;
Step S7: the relational expression of the curvature first order derivative of above-mentioned the first power function relationship formula is reduced to:
k 1 ′ = a 1 b 1 ( b 1 - 1 ) x 1 ( b 1 - 3 ) [ ( b 1 - 2 ) ( 1 + a 1 2 b 1 2 x 1 ( 2 b 1 - 2 ) ) 3 2 - 3 a 1 2 b 1 2 ( b 1 - 1 ) x 1 ( 2 b 1 - 2 ) ( 1 + a 1 2 b 1 2 x 1 ( 2 b 1 - 2 ) ) 1 2 ] ( 1 + a 1 2 b 1 2 x 1 ( 2 b 1 - 2 ) ) 3
The relational expression of the curvature first order derivative of above-mentioned the second power function relationship formula is reduced to:
k 2 ′ = a 2 b 2 ( b 2 - 1 ) x 2 ( b 2 - 3 ) [ ( b 2 - 2 ) ( 1 + a 2 2 b 2 2 x 2 ( 2 b 2 - 2 ) ) 3 2 - 3 a 2 2 b 2 2 ( b 2 - 1 ) x 2 ( 2 b 2 - 2 ) ( 1 + a 2 2 b 2 2 x 2 ( 2 b 2 - 2 ) ) 1 2 ] ( 1 + a 2 2 b 2 2 x 2 ( 2 b 2 - 2 ) ) 3
Step S8: for the relational expression of the curvature first order derivative of the first power function relationship formula, due to a in the first power function relationship formula 1and b 1all non-vanishing, if a 1and b 1have one to be zero, the first power function relationship formula becomes y 1=0, now the curvature of the first power function relationship formula does not have maximum point; And b 1be not 1, if b 1be that 1 the first power function relationship formula becomes y 1=a 1x 1, there is no point of maximum curvature; Meanwhile, x 1be greater than zero, due to x 1represent core porosity data, so x 1must be greater than zero, therefore the first order derivative of the curvature of the first power function relationship formula is that the condition of 1 o'clock is:
( b 1 - 2 ) + ( b 1 - 2 ) a 1 2 b 1 2 x 1 ( 2 b 1 - 2 ) - 3 ( b 1 - 1 ) a 1 2 b 1 2 x 1 ( 2 b 1 - 2 ) = 0
Have:
x 1 ( 2 b 1 - 2 ) = b 1 - 2 ( 2 b 1 - 1 ) a 1 2 b 1 2
By a 1and b 1substitution above formula is also obtained x 1can obtain core porosity data lower limit;
For the relational expression of the curvature first order derivative of the second power function relationship formula, due to a in the second power function relationship formula 2and b 2all non-vanishing, if a 2and b 2have one to be zero, the second power function relationship formula becomes y 2=0, now the curvature of the second power function relationship formula does not have maximum point; And b 2be not 1, if b 2be that 1 the second power function relationship formula becomes y 2=a 2x 2, there is no point of maximum curvature; Meanwhile, x 2be greater than zero, due to x 2represent core permeability data, so x 2must be greater than zero, therefore the first order derivative of the curvature of the second power function relationship formula is that the condition of 1 o'clock is:
( b 2 - 2 ) + ( b 2 - 2 ) a 2 2 b 2 2 x 2 ( 2 b 2 - 2 ) - 3 ( b 2 - 1 ) a 2 2 b 2 2 x 2 ( 2 b 2 - 2 ) = 0
Have:
x 2 ( 2 b 2 - 2 ) = b 2 - 2 ( 2 b 2 - 1 ) a 2 2 b 2 2
By a 2and b 2substitution above formula is also obtained x 2can obtain core permeability data lower limit.
CN201410307584.7A 2014-06-30 2014-06-30 The acquisition methods of mine reservoir net thickness lower limit Expired - Fee Related CN104122182B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410307584.7A CN104122182B (en) 2014-06-30 2014-06-30 The acquisition methods of mine reservoir net thickness lower limit

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410307584.7A CN104122182B (en) 2014-06-30 2014-06-30 The acquisition methods of mine reservoir net thickness lower limit

Publications (2)

Publication Number Publication Date
CN104122182A true CN104122182A (en) 2014-10-29
CN104122182B CN104122182B (en) 2016-04-13

Family

ID=51767681

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410307584.7A Expired - Fee Related CN104122182B (en) 2014-06-30 2014-06-30 The acquisition methods of mine reservoir net thickness lower limit

Country Status (1)

Country Link
CN (1) CN104122182B (en)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104912550A (en) * 2015-05-05 2015-09-16 中国海洋石油总公司 Method for quantitatively calculating reservoir fluid producing profile by nuclear magnetic resonance well logging information
CN109358377A (en) * 2018-11-30 2019-02-19 中国石油大港油田勘探开发研究院 A kind of determining clastic rock Heterogeneous reservoir at hiding lower limit method
CN114184529A (en) * 2020-09-15 2022-03-15 中国石油化工股份有限公司 Method, device, electronic equipment and medium for acquiring irreducible water saturation
CN114961711A (en) * 2021-02-24 2022-08-30 中国石油化工股份有限公司 Method for determining effective thickness limit in low-permeability reservoir reserve calculation

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4785882A (en) * 1987-06-24 1988-11-22 Mobil Oil Corporation Enhanced hydrocarbon recovery
EP1519186A1 (en) * 2003-09-27 2005-03-30 Bayerische Motoren Werke Aktiengesellschaft Test method and apparatus for air-permeable materials
CN102645678A (en) * 2012-05-08 2012-08-22 中国石油大学(华东) Lower limit computation method for reservoir forming physical property of effective reservoir layer under restraint of reservoir forming power and pore structure
CN103278436A (en) * 2013-02-01 2013-09-04 西安石油大学 Quantitative characterization method of low penetration double-medium sandstone oil reservoir microscopic aperture structure
CN103645516A (en) * 2013-11-20 2014-03-19 中国石油大学(北京) Method of determining oil and gas productivity based on oil and gas control effects of petrophysical facies
CN103778328A (en) * 2014-01-10 2014-05-07 中国石油大学(华东) Reservoir sensitive pore throat extracting method based on standard-deviation analysis

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4785882A (en) * 1987-06-24 1988-11-22 Mobil Oil Corporation Enhanced hydrocarbon recovery
EP1519186A1 (en) * 2003-09-27 2005-03-30 Bayerische Motoren Werke Aktiengesellschaft Test method and apparatus for air-permeable materials
CN102645678A (en) * 2012-05-08 2012-08-22 中国石油大学(华东) Lower limit computation method for reservoir forming physical property of effective reservoir layer under restraint of reservoir forming power and pore structure
CN103278436A (en) * 2013-02-01 2013-09-04 西安石油大学 Quantitative characterization method of low penetration double-medium sandstone oil reservoir microscopic aperture structure
CN103645516A (en) * 2013-11-20 2014-03-19 中国石油大学(北京) Method of determining oil and gas productivity based on oil and gas control effects of petrophysical facies
CN103778328A (en) * 2014-01-10 2014-05-07 中国石油大学(华东) Reservoir sensitive pore throat extracting method based on standard-deviation analysis

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
崔永斌: "有效储层物性下限值的确定方法", 《国外测井技术》 *
张春: "有效储层基质物性下限确定方法", 《油气地球物理》 *
曾伟等: "川东嘉二段孔隙层下限及分类与评价", 《矿物岩石》 *

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104912550A (en) * 2015-05-05 2015-09-16 中国海洋石油总公司 Method for quantitatively calculating reservoir fluid producing profile by nuclear magnetic resonance well logging information
CN104912550B (en) * 2015-05-05 2018-08-03 中国海洋石油集团有限公司 A kind of NMR logging data quantitatively calculates reservoir production profile method
CN109358377A (en) * 2018-11-30 2019-02-19 中国石油大港油田勘探开发研究院 A kind of determining clastic rock Heterogeneous reservoir at hiding lower limit method
CN114184529A (en) * 2020-09-15 2022-03-15 中国石油化工股份有限公司 Method, device, electronic equipment and medium for acquiring irreducible water saturation
CN114184529B (en) * 2020-09-15 2024-05-07 中国石油化工股份有限公司 Method, device, electronic equipment and medium for acquiring irreducible water saturation
CN114961711A (en) * 2021-02-24 2022-08-30 中国石油化工股份有限公司 Method for determining effective thickness limit in low-permeability reservoir reserve calculation

Also Published As

Publication number Publication date
CN104122182B (en) 2016-04-13

Similar Documents

Publication Publication Date Title
CN101892837B (en) Formation factor determining method and oil saturation determining method
CN103279647B (en) A kind of method calculating Sandstone Gas Reservoir permeability based on pore character parameter
CN104453874B (en) Glutenite reservoir oil saturation calculation method based on nuclear magnetic resonance
CN104268377B (en) A kind of friability index defining method based on coal petrography industry component
CN104634804B (en) A kind of method utilizing nuclear magnetic resonance T 2 spectrum to determine reservoir relative permeability
CN104806232B (en) A kind of method for determining porosity lower limit of fracture
CN101413388A (en) Method for obtaining oil-water common-layer original oil-containing saturation degree and method for estimating non-test oil-water common-layer original oil-containing saturation degree
MX2012001961A (en) Reservoir architecture and connectivity analysis.
CN104564042A (en) Method for evaluating brittleness of shale reservoir
CN105510203A (en) Method for determination of sandstone oil reservoir oil-water relative permeability under different temperature gradients
CN104899411B (en) A kind of reservoir productivity prediction model method for building up and system
CN108107182A (en) Free gas and adsorbed gas content determination method in a kind of shale gas
CN104122182B (en) The acquisition methods of mine reservoir net thickness lower limit
WO2006120366A8 (en) Methods for analysis of pressure response in underground formations
CN104948176A (en) Method for identifying carbonate reservoir fractures based on permeability increasing rate
CN105930932A (en) Gas index-based shale-gas-layer standardized open-flow capacity obtaining method
CN106355571A (en) Determination method and device for dolostone reservoir quality
CN105372298A (en) Analysis method for formation rock mineral constituents
CN106897531A (en) A kind of method for quantitatively evaluating of hyposmosis limestone reservoir permeability
CN106771071B (en) A kind of sealing core drilling saturation correction method mutually seeped based on grease
CN105804737A (en) Method for solving formation porosity on basis of iterative algorithm
CN105952446B (en) A kind of petroleum gas gas-bearing formation component content measuring method
CN104374679A (en) Method for determining permeability contribution value of reservoir
CN105003258B (en) A kind of acquisition methods of HTHP gas-bearing formation methane stream density matrix parameter
GB2460210A (en) Fluid conductivity measurement tool and methods

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20160413

Termination date: 20210630