CN105386753A - Method for constructing pseudo capillary pressure curves by using NMR (nuclear magnetic resonance) logging - Google Patents
Method for constructing pseudo capillary pressure curves by using NMR (nuclear magnetic resonance) logging Download PDFInfo
- Publication number
- CN105386753A CN105386753A CN201510711430.9A CN201510711430A CN105386753A CN 105386753 A CN105386753 A CN 105386753A CN 201510711430 A CN201510711430 A CN 201510711430A CN 105386753 A CN105386753 A CN 105386753A
- Authority
- CN
- China
- Prior art keywords
- mercury
- centerdot
- capillary pressure
- pseudo
- magnetic resonance
- 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.)
- Pending
Links
- 238000005481 NMR spectroscopy Methods 0.000 title claims abstract description 72
- 238000000034 method Methods 0.000 title claims abstract description 72
- 239000011148 porous material Substances 0.000 claims abstract description 104
- 238000001228 spectrum Methods 0.000 claims abstract description 26
- 230000009466 transformation Effects 0.000 claims abstract description 6
- QSHDDOUJBYECFT-UHFFFAOYSA-N mercury Chemical compound [Hg] QSHDDOUJBYECFT-UHFFFAOYSA-N 0.000 claims description 134
- 229910052753 mercury Inorganic materials 0.000 claims description 134
- 239000011435 rock Substances 0.000 claims description 51
- 238000002474 experimental method Methods 0.000 claims description 14
- 238000005259 measurement Methods 0.000 claims description 11
- 230000008569 process Effects 0.000 claims description 9
- 239000011159 matrix material Substances 0.000 claims description 7
- 238000000205 computational method Methods 0.000 claims description 4
- 238000005315 distribution function Methods 0.000 claims description 4
- 230000015572 biosynthetic process Effects 0.000 claims description 3
- 239000000463 material Substances 0.000 claims description 2
- 238000011156 evaluation Methods 0.000 abstract description 4
- 238000002347 injection Methods 0.000 abstract 7
- 239000007924 injection Substances 0.000 abstract 7
- 238000011158 quantitative evaluation Methods 0.000 abstract 1
- 238000010586 diagram Methods 0.000 description 10
- 230000035699 permeability Effects 0.000 description 9
- 238000012360 testing method Methods 0.000 description 8
- 238000006243 chemical reaction Methods 0.000 description 6
- 238000004458 analytical method Methods 0.000 description 3
- 238000013461 design Methods 0.000 description 3
- 230000011218 segmentation Effects 0.000 description 3
- 238000011161 development Methods 0.000 description 2
- 230000018109 developmental process Effects 0.000 description 2
- 230000005311 nuclear magnetism Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000000052 comparative effect Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000000280 densification Methods 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000005251 gamma ray Effects 0.000 description 1
- 230000000877 morphologic effect Effects 0.000 description 1
- 238000004445 quantitative analysis Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
- 238000009736 wetting Methods 0.000 description 1
Classifications
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B47/00—Survey of boreholes or wells
- E21B47/06—Measuring temperature or pressure
Landscapes
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Geology (AREA)
- Mining & Mineral Resources (AREA)
- Geophysics (AREA)
- Environmental & Geological Engineering (AREA)
- Fluid Mechanics (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geochemistry & Mineralogy (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention discloses a method for constructing pseudo capillary pressure curves by using NMR (nuclear magnetic resonance) logging. The method comprises the steps as follows: according to NMR logging data measured actually and limited capillary pressure data, seven different relaxation time is set; NMR T2 spectra are divided into eight parts; the percentage content of porosity components of the eight parts is calculated; a transformation model between the mercury-injection saturability and the percentage content of the NMR porosity components under different mercury-injection pressures is built; according to the transformation model, the mercury-injection saturability under different mercury-injection pressures is calculated; according to the calculated mercury-injection saturability and the mercury-injection pressure corresponding to the mercury-injection saturability, the pseudo capillary pressure curves are constructed; and according to the pseudo capillary pressure curves, pore throat radius distribution is obtained, evaluation parameters of a pore structure of a reservoir stratum are calculated, and the purpose of constructing the pseudo capillary pressure curves by using NMR logging data accurately and carrying out continuous quantitative evaluation on the pore structure of the reservoir stratum is achieved.
Description
Technical field
The invention belongs to reservoir evaluation field, particularly one utilizes nuclear magnetic resonance log to construct pseudo-capillary pressure curve method.
Background technology
Along with country is to the continuous increase of energy demand, wide variety of conventional reservoir is developed, and its output reduces day by day.Strengthen the exploration and development to low porosity and low permeability-fine and close reservoir, become country and solved oil gas supply shortage, guaranteed the important channel of energy security.Low porosity and low permeability-fine and close reservoir often has that reservoir densification, resistivity contrast are low, complex pore structure and the feature such as non-homogeneity is strong, cause conventional reservoir evaluation method bad to its effect, reservoir potential is in a large number mistaken for water layer or nonreservoir, reservoir capability forecasting difficulty.In order to improve the exploration efficiency for complicated reservoir, reduce development risk, carrying out quantitative assessment to Reservoir Pore Structure is effective method.
At present, capillary pressure curve evaluates Reservoir Pore Structure, the most effective data of division of reservoir type.But limit by experiment objective condition, the capillary pressure data obtained by core experiment is very limited, its continuous and quantitative cannot be utilized to evaluate Reservoir Pore Structure.The object of Reservoir Pore Structure is evaluated in order to realize continuous and quantitative, method the most frequently used at present utilizes NMR logging data continuous structure to go out pseudo-capillary pressure curve exactly, and utilizes it to replace rock core capillary pressure curve to carry out the pore structure that continuous and quantitative evaluates actual reservoir.
About the method utilizing NMR logging data to construct pseudo-capillary pressure curve has following several in the document delivered at present:
1, based on the pseudo-capillary pressure curve building method of Swanson parameter.See in June, 2008 " AppliedGeophysics " magazine, the article of " the ANewMethodtoConstructReservoirCapillaryPressureCurvesUsi ngNMRLoggingDataandItsApplication " of people's works such as XiaoLiang, describes the pseudo-capillary pressure curve building method based on Swanson parameter.The basic step of the method comprises: first, the capillary pressure curve data experimentally measured, and sets up the permeability of reservoir interpretation model based on Swanson parameter, utilizes it from NMR logging data, calculate the permeability of reservoir; Then, according to the morphological feature of capillary pressure curve, set up different enter the dependency relation entered between mercury saturation and nuclear-magnetism degree of porosity and permeability under mercury pressure, calculate into mercury saturation to utilize the porosity and permeability parameter of rock, finally, enter mercury saturation in conjunction with what calculate and enter mercury pressure, realizing the object of pseudo-capillary pressure curve reconstruct.
2, two-dimentional homalographic scale conversion coefficient method.See in February, 2009 " logging technique " magazine, the article of " application of nuclear magnetic resonance log in Reservoir Pore Structure is evaluated " of people's works such as Shao Weizhi, describes the method adopting the pseudo-capillary pressure curve of two-dimentional homalographic scale conversion coefficient method structure nuclear-magnetism.The basic step of the method comprises: first according to the relation between pore throat radius and capillary pressure, the capillary pressure curve of measurement is converted into pore throat radius distribution, utilizes the differential principle of similitude, determine the T of every block sample
2transverse conversion coefficient between spectrum distributes with corresponding pore throat radius; Then, utilize segmentation homalographic scale method, determine the T of every block sample
2longitudinal conversion coefficient between spectrum distributes with corresponding pore throat radius; Finally, the relation between transverse conversion coefficient and longitudinal conversion coefficient and rock porosity and permeability is set up, to realize utilizing nuclear magnetic resonance log to construct the object of pseudo-capillary pressure curve.
3, piecewise nonlinear power function scale method.See the patent of inventing rectify in January, 2010 etc. " a kind of method utilizing NMR logging data continuous and quantitative to evaluate Reservoir Pore Structure ", describe the method adopting piecewise nonlinear power function scale method to construct pseudo-capillary pressure curve from NMR logging data in the beginning of spring.Its general principle is: the parameter first reflecting reservoir difference according to rock porosity and computing permeability one
and utilize this parameter that reservoir is divided into four classes, for each type reservoir, adopt different segmentation power functions to construct pseudo-capillary pressure curve respectively.For the first kind, Equations of The Second Kind and the 3rd class reservoir, adopt the method for segmentation, in large pore throat and little pore throat section, adopt different power functions respectively, and for the 4th class reservoir, adopt single power function by nuclear magnetic resonance T
2spectrum is converted into pore throat radius distribution, then according to the relation between pore throat radius and capillary pressure, constructs pseudo-capillary pressure curve.
Existing method constructs in the process of pseudo-capillary pressure curve utilizing NMR logging data, has all used two intermediate parameters: the porosity and permeability of rock.But, in actual reservoir evaluation, these two intermediate parameters are all difficult to Obtaining Accurate, especially contain in gas reservoir in fine and close reservoir and complex lithology, porosity and permeability parameter is more difficult to effective acquisition, constructs pseudo-capillary pressure curve add difficulty to utilizing said method.Meanwhile, Permeability Parameters can not from nuclear magnetic resonance T
2directly obtain in spectrum, need to adopt special computational methods, the introducing of these computational processes, human error will certainly be brought, and these errors finally can be delivered in the process of establishing of pseudo-capillary pressure curve tectonic model, cause utilizing the pseudo-capillary pressure curve error of existing method construct comparatively large, do not reach the object utilizing NMR logging data accurate quantitative analysis to evaluate Reservoir Pore Structure.
Summary of the invention
In order to overcome now methodical deficiency, the invention provides a kind of method utilizing nuclear magnetic resonance log to construct pseudo-capillary pressure curve, realize the NMR logging data according to actual measurement and limited capillary pressure data, the object of the pseudo-capillary pressure curve of exact configuration.
For solving the problem, the present invention adopts following technical scheme:
According to NMR logging data and the limited capillary pressure data of actual measurement, set 7 different relaxation times, by nuclear magnetic resonance T
2spectrum is divided into 8 parts, calculate the degree of porosity component percentages of 8 parts, foundation difference enters the transformation model entered between mercury saturation and NMR porosity component percentages under mercury pressure, according to gained transformation model, calculate and different enter to enter mercury saturation under mercury pressure, according to calculate enter mercury saturation and with it correspondence enter mercury pressure, construct pseudo-capillary pressure curve, and according to pseudo-capillary pressure curve, obtain pore throat radius distribution, calculate Reservoir Pore Structure evaluating, realize utilizing the pseudo-capillary pressure curve of NMR logging data exact configuration and continuous and quantitative to evaluate the object of Reservoir Pore Structure.
The method of the invention comprises the following steps:
1) on the basis analyzing rock core pressure mercury characteristic of experiment and operating process, determine that M is entered mercury pressure;
2) utilize nuclear magnetic resonance logging instrument, target reservoir formation is carried out measuring also inverting and obtains the nuclear magnetic resonance log T of actual measurement
2spectrum;
3) given 7 different T
2relaxation time, by nuclear magnetic resonance T
2spectrum is divided into 8 parts, then calculates 8 degree of porosity component percentages;
4) set up M difference enter mercury pressure under the functional relation entered between mercury saturation and 8 degree of porosity component percentages, and utilize it to calculate from NMR logging data to enter the corresponding M of mercury pressure with M and enter mercury saturation value;
5) enter mercury saturation according to what calculate and given enter mercury pressure value, drawing the pseudo-capillary pressure curve of structure;
6) the pseudo-capillary pressure curve of structure, is utilized to obtain reservoir pore throat radius distribution, and calculate average pore throat radius, maximum pore throat radius, median radius, middle duty pressure, replacement pressure, equal value coefficient, the Reservoir Pore Structure evaluatings such as sorting coefficient.
Described step 1) M that relates to enters mercury pressure and layouts as follows:
P
c(i)＝0.005×2
i-1,i＝1,2,......,M
In formula: P
ci () is i-th and enters mercury pressure, unit is MPa.
Described step 3) given 7 T
2relaxation time is respectively T
2(1)=1.0ms, T
2(2)=3.0ms, T
2(3)=10.0ms, T
2(4)=33.0ms, T
2(5)=100.0ms, T
2(6)=300.0ms, T
2(7)=1000.0ms.For convenience of description, set: minimum T
2relaxation time is designated as T
2(0), maximum T
2relaxation time is designated as T
2(8).
Described step 3) the minimum T that defines
2relaxation time T
2(0)=0.3ms; Maximum T
2relaxation time T
2(8)=3000.0ms.
Described step 3) 8 degree of porosity component percentages relating to are respectively X1, X2, X3, X4, X5, X6, X7 and X8, and its computational methods are:
In formula, S (t) is and T
2the degree of porosity distribution function that relaxation time is relevant.
Described step 4) functional relation entered between mercury saturation and degree of porosity component percentages that relates to is as follows:
In formula
Represent that M difference is entered and enter mercury saturation under mercury pressure, form is percentage %;
Described step 4) in when utilizing NMR logging data to calculate into mercury saturation, only calculate into mercury pressure P
c(i) > 0.08MPa part enter mercury saturation, for entering mercury pressure P
ci the mercury saturation setting of entering of ()≤0.08MPa part equals 0.
Described step 6) in utilize pseudo-capillary pressure curve to obtain reservoir pore throat radius distribution and the method calculating Reservoir Pore Structure evaluating is carried out to wait method in gerneral institutes of higher education's Tenth Five-year plan teaching material " physics of oil layer " of works described in 209-233 page according to Yang Sheng.
The invention has the beneficial effects as follows: utilize NMR logging data to construct the pseudo-capillary pressure curve of continuous print, and utilize pseudo-capillary pressure curve obtain reservoir pore throat radius distribution and calculate Reservoir Pore Structure evaluating, and then realize the object evaluating Reservoir Pore Structure based on NMR logging data continuous and quantitative.
Accompanying drawing explanation
Below in conjunction with drawings and Examples, the invention will be further described.
Fig. 1 is that one provided by the invention utilizes NMR logging data to construct pseudo-capillary pressure curve method flow diagram.
Fig. 2 is the pressure mercury capillary pressure curve schematic diagram of the four kinds of different pore structures core samples in reservoir, domestic southwest that the embodiment of the present invention provides.
Fig. 3 is the nuclear magnetic resonance T of the four kinds of different pore structures core samples in reservoir, domestic southwest that the embodiment of the present invention provides
2spectrum schematic diagram.
Fig. 4 is that the pressure mercury capillary pressure curve of the better rock of pore structure that the embodiment of the present invention provides and the pseudo-capillary pressure curve form of structure contrast schematic diagram.
Fig. 5 is that the pressure mercury capillary pressure curve of the medium rock of pore structure that the embodiment of the present invention provides and the pseudo-capillary pressure curve form of structure contrast schematic diagram.
Fig. 6 is that the pressure mercury capillary pressure curve of the poor rock of pore structure that the embodiment of the present invention provides and the pseudo-capillary pressure curve form of structure contrast schematic diagram.
Fig. 7 be the embodiment of the present invention provide utilize the method for the invention Reservoir Pore Structure evaluating that the pseudo-capillary pressure curve that constructs and the pseudo-capillary pressure curve of utilization calculate from the NMR logging data of actual measurement and core experiment Comparative result design sketch.
Detailed description of the invention
Theory analysis
Below with the pressure mercury capillary pressure curve of certain 4 pieces of representative core sample in reservoir of Southwestern China area and nuclear magnetic resonance T
2spectrum is example, and general principle of the present invention and thinking are described.
See Fig. 2, be depicted as the pressure mercury capillary pressure experimental result of 4 pieces of samples, wherein, X-axis is for entering mercury saturation, and unit is percentage %; Y-axis is for entering mercury pressure, and unit is MPa.The form of the pressure mercury capillary pressure curve of 4 pieces of samples is completely different, and its pore structure representing rock also exists notable difference.The capillary pressure curve of No. 1 sample is arranged in the lower left corner of figure, and the pore structure representing rock is best; The capillary pressure curve of No. 2 samples is positioned on the capillary pressure curve of No. 1 sample, represents its pore structure comparatively No. 1 sample difference; The capillary pressure curve of No. 3 samples is positioned on the capillary pressure curve of No. 1 and No. 2 sample, under the capillary pressure curve being positioned at No. 4 samples, represents its pore structure comparatively No. 1 sample and No. 2 sample differences, and is better than No. 4 samples; The capillary pressure curve of No. 4 samples is arranged in the top of figure, represents its pore structure the poorest.
Pressure mercury capillary pressure experimental data corresponding to 4 pieces of samples is as shown in table 1.As shown in Table 1, carrying out pressing in the process of mercury experiment, to all samples apply to enter mercury pressure all identical, namely in Fig. 2, the Y-axis coordinate of all capillary pressure curves is all equal.The applying mode entering mercury pressure is shown below:
P
c(i)＝0.005×2
i-1,i＝1,2,…,13
In formula: P
ci () is that i-th of applying enters mercury pressure, unit is MPa.
Due to all core samples carry out applying when pressing mercury experiment to enter mercury pressure all identical, therefore, the form of the capillary pressure curve of core sample is mainly subject to control into the size of mercury saturation.When apply identical enter mercury pressure, for the good rock of pore structure, press-in blowhole space to enter mercury quantity more, it is corresponding that to enter mercury saturation higher; Otherwise, for the rock that pore structure is poor, press-in blowhole space to enter mercury quantity less, it is corresponding that to enter mercury saturation lower.
As can be seen here, only need to calculate difference and enter and enter mercury saturation under mercury pressure, in conjunction with fixing apply enter mercury pressure, can capillary pressure curve be reconstructed.
Table 1 four pieces of representative core sample pressure mercury capillary pressure experimental data tables
For core sample, when apply enter mercury pressure less (being less than 0.08MPa) when, almost the mercury of wetting phase is squeezed into the interstitial space of rock nothing but, namely when entering mercury pressure and being less than 0.08MPa, rock to enter mercury saturation less, substantially can ignore, and the form of this part capillary pressure curve is less for the contribution evaluating rock pore structure.Therefore, the present invention propose only estimate into mercury pressure be greater than 0.08MPa part enter mercury saturation to construct pseudo-capillary pressure curve, being less than 0.08MPa part for entering mercury pressure, setting and equaling 0 into mercury saturation.
See Fig. 3, be depicted as the nuclear magnetic resonance T of 4 pieces of samples
2the experimental result of spectrum, wherein X-axis is T
2in the relaxation time, unit is ms; Y-axis is relative amplitude, and unit is v/v.
Comparison diagram 2 and Fig. 3 can find, the good rock of pore structure, as in Fig. 2 No. 1 sample, and rock NMR T corresponding in figure 3
2keep right in the position of spectrum main peak, and T
2the wider distribution of spectrum, is mainly distributed between 2 ~ 1600ms, is distributed as master with macrovoid, and the proportion of macrovoid component shared by rock total pore space is larger; For No. 2 samples that pore structure is slightly poor, its nuclear magnetic resonance T
2the distribution of spectrum is narrow compared with No. 1 sample, and the distribution of main peak keeps left relatively, the proportion increase that fine pore component is shared in rock total pore space; For worse No. 3 samples of pore structure, T
2the distribution of spectrum compared with No. 1 and No. 2 samples narrow, the position of main peak is moved to the left further, and the proportion of fine pore component shared by rock total pore space increases further, and the proportion shared by macrovoid component reduces relatively; For No. 4 samples that pore structure is the poorest, its T
2the position of spectrum main peak is arranged in the Far Left of Fig. 3, and rock is mainly distributed as master with fine pore, and the proportion of macrovoid component shared by rock total pore space is minimum.
The above analysis can find, the proportion that in rock, each hole component is shared in rock total pore space has certain indicative function to rock pore structure.Therefore, its different from rock enter entering mercury saturation and should have stronger correlation under mercury pressure.In theory, identical enter under mercury pressure, for the good rock of pore structure, enter mercury saturation higher, the proportion shared in total pore space of macrovoid component is comparatively large, and component proportion shared in total pore space in fine pore is less; Otherwise, for the rock that pore structure is poor, enter mercury saturation and reduce, the proportion reduction that macrovoid component is shared in total pore space, and the proportion increase that fine pore component is shared in total pore space.
The present invention is on the basis of above-mentioned experimental analysis, the percentage composition that proposition rock each several part hole component accounts for rock total pore space component to enter mercury saturation under calculating difference and entering mercury pressure, then combine calculate enter mercury saturation and given enter mercury pressure value to construct the method for pseudo-capillary pressure curve.
In order to visual rationing characterizes different aperture component proportion shared in rock total pore space, adopt the T that 7 different
2relaxation time T
2(1), T
2(2), T
2(3), T
2(4), T
2(5), T
2and T (6)
2(7), by rock NMR T
2spectrum is divided into 8 intervals, then calculates 8 degree of porosity component percentages according to the radiometer of rock porosity component sum and total porosity in 8 intervals respectively.Concrete computational process is as follows:
In formula, S (t) is and T
2the degree of porosity distribution function that relaxation time is relevant.
After calculating above-mentioned 8 degree of porosity component percentages, set up and different enter to enter functional relation between mercury saturation and above-mentioned 8 degree of porosity component percentages under mercury pressure, realization utilizes NMR logging data to calculate into mercury saturation, then combines and fixing enter the object that mercury pressure value constructs pseudo-capillary pressure curve.
Embodiment 1
See Fig. 1, one utilizes nuclear magnetic resonance log to construct pseudo-capillary pressure curve method, and step is as follows:
1), according to feature and the operating process of pressing mercury experiment, determine that 13 are entered mercury pressure;
2), utilize nuclear magnetic resonance logging instrument, target reservoir formation is carried out measuring also inverting and obtains the nuclear magnetic resonance log T of actual measurement
2spectrum;
3), given 7 different T
2relaxation time T
2(1), T
2(2), T
2(3), T
2(4), T
2(5), T
2and T (6)
2(7), by nuclear magnetic resonance T
2spectrum is divided into 8 parts, then calculates 8 degree of porosity component percentages;
4), set up 13 differences enter mercury pressure under the functional relation entered between mercury saturation and 8 degree of porosity component percentages, and utilize it from NMR logging data, calculate 13 to enter mercury saturation value;
5), with calculate 13 enter mercury saturation value for abscissa, enter mercury pressure value for ordinate with given 13, set up plane right-angle coordinate, draw pseudo-capillary pressure curve;
6), according to the relation between pore throat radius and capillary pressure, obtain the distribution of reservoir rocks pore throat radius, and calculate average pore throat radius, maximum pore throat radius, median radius, middle duty pressure, replacement pressure, all Reservoir Pore Structure such as value coefficient and sorting coefficient evaluating, to realize the object that continuous and quantitative evaluates Reservoir Pore Structure.
Described step 1) relate to 13 enter mercury pressure value and adopt following formula to determine:
P
c(i)＝0.005×2
i-1,i＝1,2,…13
In formula: P
ci () is that i-th of applying enters mercury pressure, unit is MPa.
Described step 3) given 7 T
2relaxation time is respectively T
2(1)=1.0ms, T
2(2)=3.0ms, T
2(3)=10.0ms, T
2(4)=33.0ms, T
2(5)=100.0ms, T
2(6)=300.0ms and T
2(7)=1000.0ms.
Described step 3) the minimum T that defines
2relaxation time T
2(0)=0.3ms; Maximum T
2relaxation time T
2(8)=3000.0ms.
Described step 3) 8 degree of porosity component percentages relating to are respectively X1, X2, X3, X4, X5, X6, X7 and X8, and its computational methods are:
In formula, S (t) is and T
2the degree of porosity distribution function that relaxation time is relevant.
Described step 4) relate to 13 enter mercury saturation and adopt following formulae discovery:
In formula,
For coefficient matrix undetermined;
The defining method of described coefficient matrix and constant matrices is as follows:
Drill through a part of representative core sample, carry out nuclear magnetic resonance and the experiment of pressure mercury, first, according to 7 given T simultaneously
2relaxation time, by nuclear magnetic resonance T
2spectrum is divided into 8 parts, and calculates 8 rock porosity component percentages as independent variable; Then, according to pressure mercury capillary pressure experimental data, determine that 13 are entered mercury saturation value as dependent variable; Finally, adopt the method for multiple linear statistical regression, determine the numerical value of coefficient matrix and constant matrices respectively.
The numerical value of first to fourth row of described coefficient matrix and constant matrices is all set to 0, and the mercury saturation value of entering of calculating is also 0.
Embodiment 2
Nuclear magnetic resonance log is utilized to construct pseudo-capillary pressure curve method according to of the present invention, 20 pieces of core samples that reservoir, domestic southwest drills through are processed, demarcate the numerical value of coefficient matrix and constant matrices in mercury saturation design formulas, obtain the model of continuous structure capillary pressure curve.Utilize the NMR logging data of model to actual measurement of demarcating to process, obtain the pseudo-capillary pressure curve of continuous print.In order to evaluate the reliability of capillary pressure curve tectonic model, have read the pseudo-capillary pressure curve with 20 pieces of corresponding degree of depth of core sample, and itself and the pressure mercury capillary pressure curve of testing are contrasted.
See Fig. 4 to Fig. 6, list that pore structure is better, pore structure is medium respectively and the pseudo-capillary pressure curve comparison diagram of three class representative cores pressure mercury capillary pressure curves that pore structure is poor and structure.As we can see from the figure, for various dissimilar reservoir rocks, the pressure mercury capillary pressure curve that the pseudo-capillary pressure curve utilizing the method for the invention to construct and core experiment obtain all coincide better.Table 2 is the error statistics table between 20 pieces of core samples Reservoir Pore Structure evaluating of adopting the pseudo-capillary pressure curve of structure to calculate and the rock core Reservoir Pore Structure evaluating of pressing mercury to test to obtain.Can see from statistics, error between the Reservoir Pore Structure evaluating that the pseudo-capillary pressure curve utilizing the method for the invention to obtain calculates and the result of pressing mercury to test to obtain by rock core is less, its relative error is all less than 10.0%, meets the requirement of actual Reservoir Pore Structure quantitative assessment.
The Reservoir Pore Structure evaluating error statistics table that table 220 piece core sample distinct methods calculates
See Fig. 7, the pseudo-capillary pressure curve constructed for utilizing the method for the invention and the Reservoir Pore Structure evaluating of calculating and rock core press the comparison diagram of mercury experimental result.Design sketch shown in Fig. 7 is divided into Shi Erdao, and in figure, first comprises gamma ray curve (GR) and CAL (CAL), is mainly used in identifying effective sandstone reservoir; Second is depth track, unit m; 3rd road is deep lateral resistivity curve (RT) and shallow side direction resistivity curve (RXO); 4th road comprises density log (DEN) curve, neutron well logging (CNL) curve and acoustic travel time logging (AC) curve, is mainly used in the degree of porosity calculating reservoir; 5th road comprises the nuclear magnetic resonance log T of actual measurement
2spectrum T2_DIST; 6th road comprises the nuclear magnetic resonance log T that the method provided according to the embodiment of the present invention utilizes actual measurement
2the pseudo-capillary pressure curve of spectrum continuous structure and rock core press the contrast of mercury capillary pressure curve, in figure, black curve is the pseudo-capillary pressure curve utilizing NMR logging data continuous structure, in order to the convenience of curve display, in the present invention, pseudo-capillary pressure curve is added to maximum enter the state of mercury pressure 20.48MPa, discrete curve is then rock core pressure mercury capillary pressure curve; 7th road is depicted as the reservoir pore throat radius distribution utilizing the pseudo-capillary pressure curve of method construct shown in the embodiment of the present invention to convert to and the comparison diagram utilizing the pore throat radius of pressing mercury capillary pressure curve to obtain to distribute, and curve discrete in figure is the pore throat radius distribution obtained according to rock core pressure mercury capillary pressure curve.As can be seen from the result that the 6th road and the 7th road show, the distribution of the pore throat radius of the pseudo-capillary pressure curve of the method construct described in the embodiment of the present invention and acquisition is utilized to press mercury experimental result to have good uniformity with rock core.The average pore throat radius that the pseudo-capillary pressure curve that in 8th road, RM is the method continuous structure utilizing the embodiment of the present invention to provide calculates, CRM is that rock core pressure mercury tests the average pore throat radius obtained; The median radius that the pseudo-capillary pressure curve that in 9th road, RC50 is the method continuous structure utilizing the embodiment of the present invention to provide calculates, CR50 is that rock core pressure mercury tests the median radius obtained; The maximum pore throat radius that the pseudo-capillary pressure curve that in tenth road, RMAX is the method continuous structure utilizing the embodiment of the present invention to provide calculates, CRMAX is that rock core pressure mercury tests the maximum pore throat radius obtained.10th PD in one is the replacement pressure that the pseudo-capillary pressure curve of the method continuous structure utilizing the embodiment of the present invention to provide calculates, and CPD is that rock core pressure mercury tests the replacement pressure obtained; Duty pressure during the pseudo-capillary pressure curve that in 12 road, PC50 is the method continuous structure utilizing the embodiment of the present invention to provide calculates, CP50 is that rock core pressure mercury tests the middle duty pressure obtained.Upper as can be seen from figure, the Reservoir Pore Structure evaluating utilizing the pseudo-capillary pressure curve of the method for the invention continuous structure to calculate is close to rock core pressure mercury experimental result, and this illustrates, utilizes the method for the invention can by nuclear magnetic resonance log T
2spectrum converts pseudo-capillary pressure curve to, continuously to obtain reservoir pore throat radius distribution and Reservoir Pore Structure evaluating accurately.
Last it is noted that obviously, above-described embodiment is only for example of the present invention is clearly described, and the restriction not to embodiment.For those of ordinary skill in the field, can also make other changes in different forms on the basis of the above description.Here exhaustive without the need to also giving all embodiments.And thus the apparent change of extending out or variation be still among protection scope of the present invention.
Claims (9)
1. utilize nuclear magnetic resonance log to construct a pseudo-capillary pressure curve method, it is characterized in that: described method, according to the NMR logging data of actual measurement and limited pressure mercury capillary pressure data, sets 7 different relaxation times, by nuclear magnetic resonance T
2spectrum is divided into 8 parts, calculate the degree of porosity component percentages of 8 parts, foundation difference enters the transformation model entered between mercury saturation and NMR porosity component percentages under mercury pressure, according to gained transformation model, calculate and different enter to enter mercury saturation under mercury pressure, according to calculate enter mercury saturation and with it correspondence enter mercury pressure, construct pseudo-capillary pressure curve, and according to pseudo-capillary pressure curve, obtain pore throat radius distribution, calculate Reservoir Pore Structure evaluating, realize utilizing the pseudo-capillary pressure curve of NMR logging data exact configuration and continuous and quantitative to evaluate the object of Reservoir Pore Structure.
2. one as claimed in claim 1 utilizes nuclear magnetic resonance log to construct pseudo-capillary pressure curve method, it is characterized in that, said method comprising the steps of:
1) on the basis analyzing rock core pressure mercury characteristic of experiment and operating process, determine that M is entered mercury pressure;
2) utilize nuclear magnetic resonance logging instrument, target reservoir formation is carried out measuring also inverting and obtains the nuclear magnetic resonance log T of actual measurement
2spectrum;
3) given 7 different T
2relaxation time, by nuclear magnetic resonance T
2spectrum is divided into 8 parts, then calculates 8 degree of porosity component percentages;
4) set up M difference enter mercury pressure under the functional relation entered between mercury saturation and 8 degree of porosity component percentages, and utilize it to calculate from NMR logging data to enter the corresponding M of mercury pressure with M and enter mercury saturation value;
5) enter mercury saturation according to what calculate and given enter mercury pressure value, drawing the pseudo-capillary pressure curve of structure;
6) utilize the pseudo-capillary pressure curve of structure to obtain reservoir pore throat radius distribution, and calculate average pore throat radius, maximum pore throat radius, median radius, middle duty pressure, replacement pressure, equal value coefficient, the Reservoir Pore Structure evaluatings such as sorting coefficient.
3. one as claimed in claim 1 or 2 utilizes nuclear magnetic resonance log to construct pseudo-capillary pressure curve method, it is characterized in that: described step 1) relate to M enter mercury pressure and layout as follows:
P
c(i)＝0.005×2
i-1,i＝1,2,......,M
In formula: P
ci () is i-th and enters mercury pressure, unit is MPa.
4. one as claimed in claim 1 or 2 utilizes nuclear magnetic resonance log to construct pseudo-capillary pressure curve method, it is characterized in that: described step 3) given 7 T
2relaxation time is respectively T
2(1)=1.0ms, T
2(2)=3.0ms, T
2(3)=10.0ms, T
2(4)=33.0ms, T
2(5)=100.0ms, T
2(6)=300.0ms, T
2(7)=1000.0ms, for convenience of description, setting: minimum T
2relaxation time is designated as T
2(0), maximum T
2relaxation time is designated as T
2(8).
5. one as claimed in claim 1 or 2 utilizes nuclear magnetic resonance log to construct pseudo-capillary pressure curve method, it is characterized in that: described step 3) the minimum T that defines
2relaxation time T
2(0)=0.3ms; Maximum T
2relaxation time T
2(8)=3000.0ms.
6. one as claimed in claim 1 or 2 utilizes nuclear magnetic resonance log to construct pseudo-capillary pressure curve method, it is characterized in that: described step 3) 8 degree of porosity component percentages relating to are respectively X1, X2, X3, X4, X5, X6, X7 and X8, and its computational methods are:
In formula, S (t) is and T
2the degree of porosity distribution function that relaxation time is relevant.
7. one as claimed in claim 1 or 2 utilizes nuclear magnetic resonance log to construct pseudo-capillary pressure curve method, it is characterized in that: described step 4) functional relation entered between mercury saturation and degree of porosity component percentages that relates to is as follows:
In formula
Represent that M difference is entered and enter mercury saturation under mercury pressure, form is percentage %;
8. one as claimed in claim 1 or 2 utilizes nuclear magnetic resonance log to construct pseudo-capillary pressure curve method, it is characterized in that: described step 4) in when utilizing NMR logging data to calculate into mercury saturation, only calculate into mercury pressure P
c(i) > 0.08MPa part enter mercury saturation, for entering mercury pressure P
ci the mercury saturation setting of entering of ()≤0.08MPa part equals 0.
9. one as claimed in claim 1 or 2 utilizes nuclear magnetic resonance log to construct pseudo-capillary pressure curve method, it is characterized in that: described step 6) in utilize pseudo-capillary pressure curve to obtain reservoir pore throat radius distribution and the method calculating Reservoir Pore Structure evaluating is carried out to wait method in gerneral institutes of higher education's Tenth Five-year plan teaching material " physics of oil layer " of works described in 209-233 page according to Yang Sheng.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510711430.9A CN105386753A (en) | 2015-10-28 | 2015-10-28 | Method for constructing pseudo capillary pressure curves by using NMR (nuclear magnetic resonance) logging |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510711430.9A CN105386753A (en) | 2015-10-28 | 2015-10-28 | Method for constructing pseudo capillary pressure curves by using NMR (nuclear magnetic resonance) logging |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN105386753A true CN105386753A (en) | 2016-03-09 |
Family
ID=55419455
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510711430.9A Pending CN105386753A (en) | 2015-10-28 | 2015-10-28 | Method for constructing pseudo capillary pressure curves by using NMR (nuclear magnetic resonance) logging |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105386753A (en) |
Cited By (16)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106053315A (en) * | 2016-06-06 | 2016-10-26 | 中国石油大学(华东) | Method for classifying bioclastic limestone reservoir pore structures |
| CN106351652A (en) * | 2016-11-25 | 2017-01-25 | 中国地质大学（北京） | Shape correcting method for nuclear magnetic resonance logging T2 spectrum containing hydrocarbon reservoir layers |
| CN106777649A (en) * | 2016-12-09 | 2017-05-31 | 博明（北京）能源技术有限公司 | A kind of fractured reservoir pore structure method for quantitatively evaluating |
| CN106950606A (en) * | 2017-03-20 | 2017-07-14 | 中石化石油工程技术服务有限公司 | A kind of evaluation method of RESERVOIR PORE STRUCTURE |
| CN108763648A (en) * | 2018-04-26 | 2018-11-06 | 中国石油大学(北京) | Method and apparatus based on nuclear magnetic resonance T2 distributed acquisition capillary pressure curves |
| CN109100812A (en) * | 2018-07-25 | 2018-12-28 | 中国石油大学(北京) | Blowhole fractal dimension evaluation method and device based on nuclear magnetic resonance |
| CN109726430A (en) * | 2018-11-21 | 2019-05-07 | 中国石油天然气集团有限公司 | A kind of pseudo- capillary pressure curve method of nuclear-magnetism turn based on T2 spectrum corner position segmentation |
| CN111721685A (en) * | 2019-03-22 | 2020-09-29 | 中国石油化工股份有限公司 | Method and system for spectral conversion of capillary pressure curve of complex reservoir T2 |
| CN111999331A (en) * | 2019-05-27 | 2020-11-27 | 中国石油天然气股份有限公司 | Capillary pressure determination method and device based on nuclear magnetic test |
| CN112098449A (en) * | 2020-09-10 | 2020-12-18 | 中国石油大学(北京) | Method and device for acquiring pressure curve of reservoir capillary |
| CN112098291A (en) * | 2019-06-18 | 2020-12-18 | 中国石油化工股份有限公司 | Electric imaging porosity spectrum-based permeability calculation method and system |
| CN112881449A (en) * | 2021-01-22 | 2021-06-01 | 中国石油大学(北京) | Gas breakthrough pressure measuring method and correcting method |
| CN112983394A (en) * | 2021-02-07 | 2021-06-18 | 中国石油天然气股份有限公司 | Curve construction method and device based on logging data and storage medium |
| CN113818867A (en) * | 2021-09-17 | 2021-12-21 | 成都理工大学 | Method, system, medium, equipment and application for constructing pseudo capillary pressure curve |
| CN114233284A (en) * | 2021-09-28 | 2022-03-25 | 重庆科技学院 | Compact reservoir high-water-resistance layer identification method based on nuclear magnetic apparent free water porosity inversion |
| CN115506758A (en) * | 2021-06-04 | 2022-12-23 | 中国石油天然气股份有限公司 | Compact reservoir drainage pressure determination method, device, equipment and storage medium |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6008645A (en) * | 1997-03-11 | 1999-12-28 | Conoco Inc. | Prediction of permeability from capillary pressure curves derived from nuclear magnetic resonance pore size distributions |
| CN102141637A (en) * | 2010-01-28 | 2011-08-03 | 中国石油天然气股份有限公司 | Method for continuously and quantitatively evaluating pore structure of reservoir by utilizing nuclear magnetic resonance logging data |
| CN103267721A (en) * | 2013-05-03 | 2013-08-28 | 中国石油天然气集团公司 | Method for evaluating water containing characteristic and occurrence state of compact sandstone storage layer aperture |
| CN103279647A (en) * | 2013-05-03 | 2013-09-04 | 中国石油天然气集团公司 | Method for calculating tight sandstone reservoir permeability based on hole feature parameters |
-
2015
- 2015-10-28 CN CN201510711430.9A patent/CN105386753A/en active Pending
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6008645A (en) * | 1997-03-11 | 1999-12-28 | Conoco Inc. | Prediction of permeability from capillary pressure curves derived from nuclear magnetic resonance pore size distributions |
| CN102141637A (en) * | 2010-01-28 | 2011-08-03 | 中国石油天然气股份有限公司 | Method for continuously and quantitatively evaluating pore structure of reservoir by utilizing nuclear magnetic resonance logging data |
| CN103267721A (en) * | 2013-05-03 | 2013-08-28 | 中国石油天然气集团公司 | Method for evaluating water containing characteristic and occurrence state of compact sandstone storage layer aperture |
| CN103279647A (en) * | 2013-05-03 | 2013-09-04 | 中国石油天然气集团公司 | Method for calculating tight sandstone reservoir permeability based on hole feature parameters |
Non-Patent Citations (2)
| Title |
|---|
| CHONG ZHANG 等: "A Novel Method to Construct Capillary Pressure Curves by Using NMR Log Data and Its Application in Reservoir Evaluation", 《2009 KUWAIT INTERNATIONAL PETROLEUM CONFERENCE AND EXHIBITION》 * |
| 刘晓鹏 等: "储层毛管压力曲线构造方法及其应用", 《西南石油大学学报(自然科学版)》 * |
Cited By (27)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106053315A (en) * | 2016-06-06 | 2016-10-26 | 中国石油大学(华东) | Method for classifying bioclastic limestone reservoir pore structures |
| CN106053315B (en) * | 2016-06-06 | 2019-02-05 | 中国石油大学(华东) | A kind of biocalcarenite RESERVOIR PORE STRUCTURE classification method |
| CN106351652A (en) * | 2016-11-25 | 2017-01-25 | 中国地质大学（北京） | Shape correcting method for nuclear magnetic resonance logging T2 spectrum containing hydrocarbon reservoir layers |
| CN106351652B (en) * | 2016-11-25 | 2019-03-29 | 中国地质大学（北京） | One kind T2 of nuclear magnetic resonance log containing hydrocarbon reservoir composes morphological correction method |
| CN106777649A (en) * | 2016-12-09 | 2017-05-31 | 博明（北京）能源技术有限公司 | A kind of fractured reservoir pore structure method for quantitatively evaluating |
| CN106777649B (en) * | 2016-12-09 | 2020-06-12 | 博明（北京）能源技术有限公司 | Quantitative evaluation method for pore structure of fractured reservoir |
| CN106950606A (en) * | 2017-03-20 | 2017-07-14 | 中石化石油工程技术服务有限公司 | A kind of evaluation method of RESERVOIR PORE STRUCTURE |
| CN108763648A (en) * | 2018-04-26 | 2018-11-06 | 中国石油大学(北京) | Method and apparatus based on nuclear magnetic resonance T2 distributed acquisition capillary pressure curves |
| CN108763648B (en) * | 2018-04-26 | 2020-05-12 | 中国石油大学(北京) | Method and device for acquiring capillary pressure curve based on nuclear magnetic resonance T2 distribution |
| CN109100812A (en) * | 2018-07-25 | 2018-12-28 | 中国石油大学(北京) | Blowhole fractal dimension evaluation method and device based on nuclear magnetic resonance |
| CN109100812B (en) * | 2018-07-25 | 2020-06-05 | 中国石油大学(北京) | Rock pore fractal dimension evaluation method and device based on nuclear magnetic resonance |
| CN109726430A (en) * | 2018-11-21 | 2019-05-07 | 中国石油天然气集团有限公司 | A kind of pseudo- capillary pressure curve method of nuclear-magnetism turn based on T2 spectrum corner position segmentation |
| CN111721685A (en) * | 2019-03-22 | 2020-09-29 | 中国石油化工股份有限公司 | Method and system for spectral conversion of capillary pressure curve of complex reservoir T2 |
| CN111721685B (en) * | 2019-03-22 | 2023-11-28 | 中国石油化工股份有限公司 | Method and system for converting T2 spectrum of complex reservoir into capillary pressure curve |
| CN111999331B (en) * | 2019-05-27 | 2023-10-31 | 中国石油天然气股份有限公司 | Capillary pressure determining method and device based on nuclear magnetic test |
| CN111999331A (en) * | 2019-05-27 | 2020-11-27 | 中国石油天然气股份有限公司 | Capillary pressure determination method and device based on nuclear magnetic test |
| CN112098291A (en) * | 2019-06-18 | 2020-12-18 | 中国石油化工股份有限公司 | Electric imaging porosity spectrum-based permeability calculation method and system |
| CN112098449A (en) * | 2020-09-10 | 2020-12-18 | 中国石油大学(北京) | Method and device for acquiring pressure curve of reservoir capillary |
| CN112098449B (en) * | 2020-09-10 | 2021-09-28 | 中国石油大学(北京) | Method and device for acquiring pressure curve of reservoir capillary |
| CN112881449A (en) * | 2021-01-22 | 2021-06-01 | 中国石油大学(北京) | Gas breakthrough pressure measuring method and correcting method |
| CN112881449B (en) * | 2021-01-22 | 2021-12-28 | 中国石油大学(北京) | Gas breakthrough pressure measuring method and correcting method |
| CN112983394A (en) * | 2021-02-07 | 2021-06-18 | 中国石油天然气股份有限公司 | Curve construction method and device based on logging data and storage medium |
| CN115506758B (en) * | 2021-06-04 | 2024-05-28 | 中国石油天然气股份有限公司 | Method, device, equipment and storage medium for determining drainage and driving pressure of tight reservoir |
| CN115506758A (en) * | 2021-06-04 | 2022-12-23 | 中国石油天然气股份有限公司 | Compact reservoir drainage pressure determination method, device, equipment and storage medium |
| CN113818867A (en) * | 2021-09-17 | 2021-12-21 | 成都理工大学 | Method, system, medium, equipment and application for constructing pseudo capillary pressure curve |
| CN114233284B (en) * | 2021-09-28 | 2022-07-15 | 重庆科技学院 | Compact reservoir high-water-resistance layer identification method based on nuclear magnetic apparent free water porosity inversion |
| CN114233284A (en) * | 2021-09-28 | 2022-03-25 | 重庆科技学院 | Compact reservoir high-water-resistance layer identification method based on nuclear magnetic apparent free water porosity inversion |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105386753A (en) | Method for constructing pseudo capillary pressure curves by using NMR (nuclear magnetic resonance) logging | |
| CN104278989B (en) | Method for obtaining saturation index of low-porosity and low-permeability reservoir | |
| CN104750896B (en) | A kind of fractured-cavernous carbonate reservoir method for numerical simulation | |
| CN103279647B (en) | A kind of method calculating Sandstone Gas Reservoir permeability based on pore character parameter | |
| CN105626009B (en) | Quantitative evaluation method for single-well water injection oil replacement effect of fracture-cave carbonate reservoir | |
| CN104634804B (en) | Method for determining relative permeability of reservoir by utilizing nuclear magnetic resonance T2 spectrum | |
| CN104819923B (en) | Low permeability sandstone reservoir pore structure quantitative inversion method based on nuclear magnetic resonance | |
| CN105353419A (en) | Nuclear magnetism capillary pressure curve construction method based on reservoir stratum classification | |
| CN104806232B (en) | A kind of method for determining porosity lower limit of fracture | |
| CN104453874B (en) | Glutenite reservoir oil saturation calculation method based on nuclear magnetic resonance | |
| CN105781539A (en) | Saturability well logging calculation method of tight oil and gas reservoir | |
| CN103590827A (en) | Dense clastic rock natural gas well productivity prediction method based on reservoir classification | |
| CN104134101A (en) | Low-permeability reservoir natural gas productivity prediction method | |
| CN104865614A (en) | Complicated reservoir fluid identification method based on variable skeleton parameter | |
| CN107780923B (en) | Method for establishing and simulating water saturation model based on argillaceous correction | |
| CN104712330A (en) | Well logging permeability interpretation method | |
| CN107767063A (en) | Quantitative evaluation parameter and determination method for geological ecological environment quality | |
| CN105447762A (en) | Calculation method for low permeability reservoir flooding information of fluid replacement | |
| CN104863574A (en) | Fluid identification method applicable to compact sandstone reservoir | |
| CN110533237A (en) | A kind of sandstone reservoir oily PRODUCTION FORECASTING METHODS | |
| Qi et al. | An analysis of upscaling techniques for reservoir simulation | |
| CN105804737A (en) | Method for solving formation porosity on basis of iterative algorithm | |
| Liu et al. | A simplified and efficient method for water flooding production index calculations in low permeable fractured reservoir | |
| CN105928827B (en) | Method for calculating optimal water content of fine-grained soil compaction test | |
| CN104514553B (en) | A kind of recognition methods and its application of rock-fragment sandstone layer Mesosphere |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20160309 |
|
| RJ01 | Rejection of invention patent application after publication |