CN113622905A - Shale reservoir brittleness evaluation method based on multi-factor comprehensive analysis - Google Patents
Shale reservoir brittleness evaluation method based on multi-factor comprehensive analysis Download PDFInfo
- Publication number
- CN113622905A CN113622905A CN202110819772.8A CN202110819772A CN113622905A CN 113622905 A CN113622905 A CN 113622905A CN 202110819772 A CN202110819772 A CN 202110819772A CN 113622905 A CN113622905 A CN 113622905A
- Authority
- CN
- China
- Prior art keywords
- shale
- brittleness
- value
- hardness
- factor
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000011156 evaluation Methods 0.000 title claims abstract description 28
- 238000004458 analytical method Methods 0.000 title claims abstract description 23
- 229910052500 inorganic mineral Inorganic materials 0.000 claims abstract description 53
- 239000011707 mineral Substances 0.000 claims abstract description 53
- 238000004364 calculation method Methods 0.000 claims abstract description 26
- 238000012360 testing method Methods 0.000 claims abstract description 19
- 238000000034 method Methods 0.000 claims abstract description 15
- 238000002474 experimental method Methods 0.000 claims abstract description 10
- 230000008569 process Effects 0.000 claims abstract description 7
- 238000002441 X-ray diffraction Methods 0.000 claims abstract description 6
- 239000011435 rock Substances 0.000 claims description 17
- 238000010606 normalization Methods 0.000 claims description 13
- 239000002734 clay mineral Substances 0.000 claims description 6
- 238000005259 measurement Methods 0.000 claims description 6
- 239000010453 quartz Substances 0.000 claims description 6
- VYPSYNLAJGMNEJ-UHFFFAOYSA-N silicon dioxide Inorganic materials O=[Si]=O VYPSYNLAJGMNEJ-UHFFFAOYSA-N 0.000 claims description 6
- BVKZGUZCCUSVTD-UHFFFAOYSA-L Carbonate Chemical compound [O-]C([O-])=O BVKZGUZCCUSVTD-UHFFFAOYSA-L 0.000 claims description 3
- BVKZGUZCCUSVTD-UHFFFAOYSA-N carbonic acid Chemical compound OC(O)=O BVKZGUZCCUSVTD-UHFFFAOYSA-N 0.000 claims description 3
- 238000002360 preparation method Methods 0.000 abstract description 3
- 239000011159 matrix material Substances 0.000 description 13
- 238000012545 processing Methods 0.000 description 4
- 239000003079 shale oil Substances 0.000 description 3
- 239000002131 composite material Substances 0.000 description 2
- 238000012937 correction Methods 0.000 description 2
- 238000005553 drilling Methods 0.000 description 2
- 239000000203 mixture Substances 0.000 description 2
- 238000012935 Averaging Methods 0.000 description 1
- 108010076282 Factor IX Proteins 0.000 description 1
- 108010023321 Factor VII Proteins 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000002349 favourable effect Effects 0.000 description 1
- 238000007542 hardness measurement Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000035699 permeability Effects 0.000 description 1
- 208000010787 postorgasmic illness syndrome Diseases 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Images
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
- E21B49/00—Testing the nature of borehole walls; Formation testing; Methods or apparatus for obtaining samples of soil or well fluids, specially adapted to earth drilling or wells
-
- 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
- E21B49/00—Testing the nature of borehole walls; Formation testing; Methods or apparatus for obtaining samples of soil or well fluids, specially adapted to earth drilling or wells
- E21B49/02—Testing the nature of borehole walls; Formation testing; Methods or apparatus for obtaining samples of soil or well fluids, specially adapted to earth drilling or wells by mechanically taking samples of the soil
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A10/00—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE at coastal zones; at river basins
- Y02A10/40—Controlling or monitoring, e.g. of flood or hurricane; Forecasting, e.g. risk assessment or mapping
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Geology (AREA)
- Mining & Mineral Resources (AREA)
- Physics & Mathematics (AREA)
- Environmental & Geological Engineering (AREA)
- Fluid Mechanics (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geochemistry & Mineralogy (AREA)
- Soil Sciences (AREA)
- Analysing Materials By The Use Of Radiation (AREA)
Abstract
The invention provides a shale reservoir brittleness evaluation method based on multi-factor comprehensive analysis, which comprises the steps of obtaining mineral brittleness index, hardness value and TOC content of a shale sample by utilizing an X-ray diffraction experiment, a hardness tester test and TOC determination, weighting the three factors by utilizing an analytic hierarchy process to obtain a comprehensive brittleness index calculation formula, and finally calculating to obtain the comprehensive brittleness index, so that the shale reservoir brittleness can be accurately evaluated, the sample preparation difficulty is low, the cost is low, and the shale core sample after the test is completely stored.
Description
Technical Field
The invention relates to the field of shale reservoir brittleness evaluation, in particular to a shale reservoir brittleness evaluation method based on multi-factor comprehensive analysis.
Background
In recent years, shale oil and gas have become hot spots for exploration and development of unconventional oil and gas resources on a global scale. The physical characteristics of low porosity and low permeability enable the successful exploitation of shale oil and gas to depend on multi-stage volume fracturing transformation to a great extent, and the brittleness of the rock determines the difficulty degree of fracturing and the form of fracturing to a great extent, and is a main factor influencing the fracturing quality and the fracturing effect. Therefore, the brittleness of shale becomes an important indicator of selecting a favorable fracture interval.
Currently, in shale oil and gas reservoir evaluation, the most common brittleness evaluation method is to calculate a brittleness index by using elasticity parameters of shale. The method is characterized in that the Young modulus and the Poisson ratio of the shale are subjected to normalization processing, and the brittleness index of the shale is calculated by taking the average value of the Young modulus and the Poisson ratio. Compared with a mineral composition calculation method, the method is more scientific and reliable, but the method has the defects that the mechanical property test is carried out through a uniaxial fracturing experiment or a triaxial fracturing experiment to obtain parameters such as Young modulus, Poisson ratio and the like, the fracturing experiment can cause the damage of a shale sample, and the subsequent test analysis is difficult to carry out.
Disclosure of Invention
In view of the above, the invention provides a shale reservoir brittleness evaluation method based on multi-factor comprehensive analysis. The mineral brittleness index, the hardness value and the TOC content of the shale sample are obtained by utilizing an X-ray diffraction experiment, a hardness tester for testing and TOC determination, the three factors are weighted by utilizing an analytic hierarchy process, and the comprehensive brittleness index is finally calculated, so that the shale reservoir brittleness can be accurately evaluated, the sample preparation difficulty is low, the cost is low, and the shale core sample after testing is completely stored.
The technical scheme of the invention is realized as follows: the invention provides a shale reservoir brittleness evaluation method based on multi-factor comprehensive analysis, which is characterized by comprising the following steps of:
s1, selecting a proper shale sample;
s2, measuring the relative content of the whole rock minerals of the shale sample, and calculating the brittleness index of the minerals;
s3, measuring the hardness value of the shale sample;
s4, measuring the TOC content of the shale sample;
s5, weighting the mineral brittleness index, the hardness value and the TOC content by using an analytic hierarchy process to obtain a comprehensive brittleness index calculation formula;
and S6, calculating the comprehensive brittleness index of shale samples at different depths according to a comprehensive brittleness index calculation formula, and evaluating the brittleness of the shale reservoir.
On the basis of the above technical solution, preferably, in step S2,
in the formula, BI1Is mineral brittleness index, dimensionless, wQuartzIs the mass fraction of quartz, wClay mineralIs the mass fraction of clay mineral, wCarbonate rockIs the mass fraction of carbonate rock.
On the basis of the above technical solution, preferably, in step S2, the relative content measurement mode of the shale sample whole rock minerals is as follows: and obtaining the relative content of the whole rock minerals of the shale sample through an X-ray diffraction experiment.
On the basis of the above technical solution, preferably, in step S3, the instrument for measuring the hardness value of the shale sample is a hardness tester, the hardness tester is used to test a plurality of groups of hardness values, and the average hardness value is taken as the reference value of the hardness of the sample.
Further preferably, in step S3, before the hardness value of the shale sample is measured, a standard hardness block is selected, the hardness is measured by using a hardness tester, and the test value is subjected to zeroing correction.
Based on the above technical solution, preferably, in step S4, the instrument for measuring the TOC content of the shale sample is a TOC determinator.
Based on the above technical solution, preferably, in step S5, the comprehensive brittleness index calculation formula is:
in the formula, BI2Is a comprehensive brittleness index;
Iiis an index of each itemWeighting coefficient of (I)1Is a weighting coefficient of the mineral brittleness index, I2As a weighting factor for hardness values, I3A weighting factor for the TOC content;
cinormalized values for the various factors after normalization, c1Is a normalized value of mineral brittleness index, c2As a normalized value of the hardness value, c3Normalized values for TOC content.
More preferably, in the calculation formula of the comprehensive brittleness index in step S5, the mineral brittleness index and hardness value are normalized positively, and the TOC content is normalized negatively.
More preferably, the step of normalizing the mineral brittleness index parameter value, the hardness parameter value and the TOC content parameter value comprises:
measuring parameter values of three factors, namely mineral brittleness index, hardness value and TOC content of shale samples at different depths, and calculating a standardized value of each factor parameter value through a standardized formula;
the forward normalization formula is:
the negative normalization formula is:
in the formula, PiAnd PjThe standardized value of the parameter value of each factor is dimensionless; xiAnd XjA parameter value that is a factor; xmaxAnd XminThe maximum value and the minimum value of the parameter values of the factors.
More preferably, in step S5, the comprehensive brittleness index calculation formula is:
BI2=0.0574×c1+0.5783×c2+0.3643×c3。
compared with the prior art, the shale reservoir brittleness evaluation method based on multi-factor comprehensive analysis has the following beneficial effects:
(1) according to the shale reservoir brittleness evaluation method based on multi-factor comprehensive analysis, three factors including the mineral brittleness index, the hardness value and the TOC content are weighted by using an analytic hierarchy process to obtain a comprehensive brittleness index calculation formula, the comprehensive brittleness index is finally calculated, accurate evaluation of shale reservoir brittleness is achieved, sample preparation difficulty is low, cost is low, and a shale core sample after test is completely stored;
(2) according to the shale reservoir brittleness evaluation method based on multi-factor comprehensive analysis, the mineral brittleness index parameter value, the hardness parameter value and the TOC content parameter value are standardized, so that the comprehensive brittleness index is calculated more accurately.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings used in the description of the embodiments or the prior art will be briefly described below, it is obvious that the drawings in the following description are only some embodiments of the present invention, and for those skilled in the art, other drawings can be obtained according to the drawings without creative efforts.
FIG. 1 is a schematic flow chart of a shale reservoir brittleness evaluation method based on multi-factor comprehensive analysis according to the present invention;
FIG. 2 is a cross-plot of the hardness versus the brittleness index of the elastic parameter for shale samples of the present invention;
FIG. 3 is a cross-plot of mineral friability index and an elastic parameter friability index for shale samples of the present invention;
FIG. 4 is a cross-plot of TOC content versus friability index, an elastic parameter, for shale samples of the present invention;
FIG. 5 is a cross-plot of the composite friability index and the elasticity parameter friability index for shale samples of the present invention.
Detailed Description
The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be obtained by a person skilled in the art without any inventive step based on the embodiments of the present invention, are within the scope of the present invention.
As shown in fig. 1, the embodiment of the invention discloses a shale reservoir brittleness evaluation method based on multi-factor comprehensive analysis, which comprises the following steps:
s1, selecting a proper shale sample;
s2, measuring the relative content of the whole rock minerals of the shale sample, and calculating the brittleness index of the minerals;
s3, measuring the hardness value of the shale sample;
s4, measuring the TOC content of the shale sample;
s5, weighting the mineral brittleness index, the hardness value and the TOC content by using an analytic hierarchy process to obtain a comprehensive brittleness index calculation formula;
and S6, calculating the comprehensive brittleness index of shale samples at different depths according to a comprehensive brittleness index calculation formula, and evaluating the brittleness of the shale reservoir.
Specifically, in step S2,
in the formula, BI1Is mineral brittleness index and has no dimension; w is aQuartzIs the mass fraction of quartz, wClay mineralIs the mass fraction of clay mineral, wCarbonate rockIs the mass fraction of carbonate rock.
The relative content measurement mode of the shale sample whole rock minerals is as follows: and obtaining the relative content of the whole rock minerals of the shale sample through an X-ray diffraction experiment.
In step S3, the instrument for measuring the hardness value of the shale sample is a hardness tester.
Because the shale has strong heterogeneity and the sensitivity of the hardness tester is strong, the test points need to be moved to complete multiple groups of tests when the hardness test is carried out. And testing multiple groups of hardness values by using a hardness tester, and taking the average hardness value as a reference value of the hardness of the sample. Typically 10 points are tested per sample and the average hardness value is taken.
Specifically, the surface of the shale sample used for testing needs to be flat, smooth and clean, and a large rock sample is optimal. The main principle of the hardness tester in hardness testing is to trigger the impact body to impact the surface of a sample, record the speed before impact and the speed after rebound, and express the hardness value by the ratio of the rebound speed to the impact speed. Therefore, in the impact process, firstly, the surface of the sample is ensured to be flat and clean, concave-convex surfaces and rough surfaces are avoided, and secondly, the sample is strictly fixed to prevent shaking. The hardness tester can be an Equotip Bambino hardness tester.
In order to prevent errors in measurement of the hardness value of the shale sample, a standard hardness block is selected before the hardness value of the shale sample is measured, a hardness tester is used for measurement, and the test value is subjected to zero resetting correction.
In step S4, the instrument for measuring the TOC content of the shale sample is a TOC meter.
In step S5, the comprehensive brittleness index calculation formula is:
in the formula, BI2Is a comprehensive brittleness index;
Iias weighting factors for the indices, I1Is a weighting coefficient of the mineral brittleness index, I2As a weighting factor for hardness values, I3A weighting factor for the TOC content;
cinormalized values for the various factors after normalization, c1Is a normalized value of mineral brittleness index, c2As a normalized value of the hardness value, c3Normalized values for TOC content.
Further, in the calculation formula of the comprehensive brittleness index of step S5, the mineral brittleness index and hardness value are normalized positively, and the TOC content is normalized negatively.
Further, the step of standardizing the mineral brittleness index parameter value, the hardness parameter value and the TOC content parameter value comprises the following steps:
measuring parameter values of three factors, namely a mineral brittleness index, a hardness value and TOC content of shale samples at different depths, calculating a standardized value of each factor parameter value through a standardized formula, and enabling the value range of each parameter value to be 0-100;
the forward normalization formula is:
the negative normalization formula is:
in the formula, PiAnd PjThe standardized value of the parameter value of each factor is dimensionless; xiAnd XjA parameter value that is a factor; xmaxAnd XminThe maximum value and the minimum value of the parameter values of the factors.
Further, in step S5, the comprehensive brittleness index calculation formula is:
BI2=0.0574×c1+0.5783×c2+0.3643×c3。
example (b):
take the Dingshan region DY5 well in Sichuan basin as an example.
(1) And observing the drilling shale core system, and sampling the drilling shale cores at different depths.
(2) Dynamic Young modulus and Poisson's ratio are obtained by calculation by utilizing dipole sound wave and density logging data, and the calculation formula is as follows:
wherein POIS represents Poisson's ratio, YMOD represents Young's modulus, DEN represents density, and V representspRepresenting the velocity, V, of longitudinal wavessRepresenting the shear wave velocity.
Carrying out normalization processing on the Young modulus and the Poisson ratio, and calculating the brittleness index BI of the elastic parameter3The correlation formula is as follows:
in the formula, YMODmaxAt the maximum value of Young's modulus, YMODminIs the minimum of Young's modulus, POISmaxAt the maximum value of the Poisson's ratio, POISminIs the minimum value of the Poisson's ratio, EBITo normalize the treated Young's modulus, vBITo normalize the Poisson's ratio after processing, BI3Is the brittleness index of the elastic parameter.
(3) And analyzing the mineral composition of each sample by using an X-ray diffraction experiment, and calculating by using a brittleness index calculation formula to obtain a mineral brittleness index.
(4) And (3) carrying out hardness test on the core sample by using an Equotip Bambino hardness tester, and testing a plurality of points on each hand specimen sample by considering that the sample has strong heterogeneity, and averaging hardness values.
(5) And determining the TOC content of each core sample by TOC.
(6) The shale samples were analyzed for their hardness, mineral friability index and TOC content in relation to the friability index, a shale reservoir elastic parameter, and the results are shown in figures 2-4.
(7) The results obtained from the above experiments and after standardization are shown in table 1, wherein the mineral brittleness index and hardness are positive indicators, the TOC content is a negative indicator, and the formula used for standardization is as follows:
TABLE 1 mineral friability index, hardness number, TOC content and standardized processing results for well core samples
(8) Combining the correlation coefficient and the experience recognition of each evaluation index and according to the judgment matrix aijAnd comparing every two of the judgment matrixes by a calibration method and assigning values to obtain a judgment matrix. Wherein, the judgment matrix aijThe scaling method is shown in Table 2 and the decision matrix is shown in Table 3.
TABLE 2 decision matrix aijScale method
Scale | Means of |
1 | Showing the same importance of the two factors compared |
3 | Indicating that one factor is slightly more important than the other factor when compared to the |
5 | Indicating that one factor is significantly more important than the other factor when compared to the other factor |
7 | Indicating that one factor is more important than the other factor |
9 | Indicating that one factor is extremely important compared to the |
2,4,6,8 | Median value of the above two adjacent judgments |
TABLE 3 decision matrix
Evaluation parameters | Index of mineral brittleness | Hardness of | TOC content |
Index of |
1 | 1/8 | 1/8 |
Hardness of | 8 | 1 | 2 |
TOC content | 8 | 1/2 | 1 |
(9) And (3) checking the consistency of the judgment matrix, wherein the checking steps are as follows:
calculating a consistency index CI, wherein the calculation formula is as follows:
in the formula, λmaxJudging the maximum eigenvalue of the matrix; n is the order of the decision matrix.
From the table lookup, when the matrix order n is 3, the random consistency index RI is 0.52.
And thirdly, calculating a consistency ratio CR, wherein the calculation formula is as follows:
CR is calculated to be 0.0516<0.1, so the consistency of the decision matrix is acceptable.
(10) Calculating a weight vector W using a square root methodiThe calculation formula is as follows:
in the formula, aijTo determine the elements in the matrix, where i, j ═ 1, 2, …, n; n isThe order of the matrix. The calculation results are shown in table 4.
TABLE 4 weight coefficient calculation results
Evaluation parameters | Index of mineral brittleness | Hardness of | TOC content | Weight coefficient |
Index of |
1 | 1/8 | 1/8 | 0.0574 |
Hardness of | 8 | 1 | 2 | 0.5783 |
TOC content | 8 | 1/2 | 1 | 0.3643 |
(11) Construction of the comprehensive brittleness index BI using the above equation 42CalculatingThe formula is as follows:
in the formula, BI2The index is a comprehensive brittleness index and has no dimension; i isiA weighting coefficient for each index, i being 1, 2, 3; c. CiThe standard indexes are dimensionless; c. C1Is a normalized value of mineral brittleness index, c2As a normalized value of the hardness value, c3Normalized values for TOC content.
(12) Analysis of shale samples for comprehensive brittleness index BI2And the brittleness index BI of the elastic parameter3The results are shown in FIG. 5.
FIGS. 2-4 are hardness, mineral friability index BI of DY5 well shale samples1And TOC content and brittleness index BI of elastic parameter3The cross-section of (1). Passing hardness parameter value and mineral brittleness index BI of each shale sample1And the TOC content parameter value and the brittleness index BI of the elastic parameter of each shale sample3In the distribution in the figure, it can be seen that the hardness and the mineral brittleness index BI of the well shale sample1And TOC content and brittleness index BI of elastic parameter3The correlation coefficients between the two are respectively 0.66, 0.42 and 0.35, and have certain correlation, but the correlation is poor.
FIG. 5 is the composite friability index BI of DY5 well shale samples2And the brittleness index BI of the elastic parameter3The cross-section of (1). The result shows that the comprehensive brittleness index BI of the well shale sample2And the brittleness index BI of the elastic parameter3The correlation coefficient between them was 0.81. Thus, the mineral brittleness index BI of the shale sample by using the analytic hierarchy process1Weighting the hardness value and the TOC content, and finally calculating the obtained comprehensive brittleness index BI2And the brittleness index BI of the elastic parameter3Highest correlation, comprehensive brittleness index BI2The brittleness of the shale can be well characterized.
The above description is only for the purpose of illustrating the preferred embodiments of the present invention and is not to be construed as limiting the invention, and any modifications, equivalents, improvements and the like that fall within the spirit and principle of the present invention are intended to be included therein.
Claims (10)
1. A shale reservoir brittleness evaluation method based on multi-factor comprehensive analysis is characterized by comprising the following steps:
s1, selecting a proper shale sample;
s2, measuring the relative content of the whole rock minerals of the shale sample, and calculating the brittleness index of the minerals;
s3, measuring the hardness value of the shale sample;
s4, measuring the TOC content of the shale sample;
s5, weighting the mineral brittleness index, the hardness value and the TOC content by using an analytic hierarchy process to obtain a comprehensive brittleness index calculation formula;
and S6, calculating the comprehensive brittleness index of shale samples at different depths according to a comprehensive brittleness index calculation formula, and evaluating the brittleness of the shale reservoir.
2. The shale reservoir brittleness evaluation method based on the multi-factor comprehensive analysis of claim 1, wherein: in the step S2, in the step S,
in the formula, BI1Is mineral brittleness index, dimensionless, wQuartzIs the mass fraction of quartz, wClay mineralIs the mass fraction of clay mineral, wCarbonate rockIs the mass fraction of carbonate rock.
3. The shale reservoir brittleness evaluation method based on the multi-factor comprehensive analysis as claimed in claim 1, wherein in step S2, the relative content measurement mode of the shale sample whole rock minerals is: and obtaining the relative content of the whole rock minerals of the shale sample through an X-ray diffraction experiment.
4. The shale reservoir brittleness evaluation method based on the multi-factor comprehensive analysis of claim 1, wherein: in step S3, the instrument for measuring the hardness value of the shale sample is a hardness tester, the hardness tester is used to test a plurality of hardness values, and the average hardness value is taken as the reference value of the hardness of the sample.
5. The shale reservoir brittleness evaluation method based on multi-factor comprehensive analysis according to claim 4, wherein: in step S3, before the hardness value of the shale sample is measured, a standard hardness block is selected, a hardness tester is used for measurement, and the zero calibration is performed on the test value.
6. The shale reservoir brittleness evaluation method based on the multi-factor comprehensive analysis of claim 1, wherein: in step S4, the instrument for measuring the TOC content of the shale sample is a TOC meter.
7. The shale reservoir brittleness evaluation method based on multi-factor comprehensive analysis according to claim 1, wherein in step S5, the comprehensive brittleness index calculation formula is:
in the formula, BI2Is a comprehensive brittleness index;
Iias weighting factors for the indices, I1Is a weighting coefficient of the mineral brittleness index, I2As a weighting factor for hardness values, I3A weighting factor for the TOC content;
cinormalized values for the various factors after normalization, c1Is a normalized value of mineral brittleness index, c2As a normalized value of the hardness value, c3Normalized values for TOC content.
8. The shale reservoir brittleness evaluation method based on the multi-factor comprehensive analysis of claim 7, wherein: in the calculation formula of the comprehensive brittleness index of the step S5, the mineral brittleness index and the hardness value are subjected to positive standardization, and the TOC content is subjected to negative standardization.
9. The shale reservoir brittleness evaluation method based on the multifactor comprehensive analysis according to claim 8, wherein the step of standardizing the mineral brittleness index parameter value, the hardness parameter value and the TOC content parameter value comprises the following steps:
measuring parameter values of three factors, namely mineral brittleness index, hardness value and TOC content of shale samples at different depths, and calculating a standardized value of each factor parameter value through a standardized formula;
the forward normalization formula is:
the negative normalization formula is:
in the formula, PiAnd PjThe standardized value of the parameter value of each factor is dimensionless; xiAnd XjA parameter value that is a factor; xmaxAnd XminThe maximum value and the minimum value of the parameter values of the factors.
10. The shale reservoir brittleness evaluation method based on multi-factor comprehensive analysis according to claim 9, wherein in step S5, the comprehensive brittleness index calculation formula is:
BI2=0.0574×c1+0.5783×c2+0.3643×c3。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110819772.8A CN113622905B (en) | 2021-07-20 | 2021-07-20 | Shale reservoir brittleness evaluation method based on multi-factor comprehensive analysis |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110819772.8A CN113622905B (en) | 2021-07-20 | 2021-07-20 | Shale reservoir brittleness evaluation method based on multi-factor comprehensive analysis |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113622905A true CN113622905A (en) | 2021-11-09 |
CN113622905B CN113622905B (en) | 2023-05-02 |
Family
ID=78380264
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110819772.8A Active CN113622905B (en) | 2021-07-20 | 2021-07-20 | Shale reservoir brittleness evaluation method based on multi-factor comprehensive analysis |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113622905B (en) |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2016041189A1 (en) * | 2014-09-19 | 2016-03-24 | 杨顺伟 | Method for evaluating shale gas reservoir and seeking desert area |
CN106248494A (en) * | 2016-08-29 | 2016-12-21 | 中国石油化工股份有限公司江汉油田分公司石油工程技术研究院 | A kind of method for shale gas well reservoir fragility overall merit |
CN108009705A (en) * | 2017-11-07 | 2018-05-08 | 中国石油大学(华东) | A kind of shale reservoir compressibility evaluation method based on support vector machines technology |
CN108073786A (en) * | 2017-12-27 | 2018-05-25 | 中国石油化工股份有限公司江汉油田分公司勘探开发研究院 | Shale compressibility evaluation method and device |
CN109407150A (en) * | 2018-12-06 | 2019-03-01 | 中国石油大学(北京) | Based on the petrophysical shale reservoir compressibility means of interpretation of statistics and system |
CN110714753A (en) * | 2019-08-09 | 2020-01-21 | 中国华能集团清洁能源技术研究院有限公司 | Shale reservoir compressibility evaluation method |
CN113090258A (en) * | 2021-05-25 | 2021-07-09 | 中国石油天然气股份有限公司 | Deep shale gas horizontal well fracturing differentiation design method based on logging data |
-
2021
- 2021-07-20 CN CN202110819772.8A patent/CN113622905B/en active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2016041189A1 (en) * | 2014-09-19 | 2016-03-24 | 杨顺伟 | Method for evaluating shale gas reservoir and seeking desert area |
CN106248494A (en) * | 2016-08-29 | 2016-12-21 | 中国石油化工股份有限公司江汉油田分公司石油工程技术研究院 | A kind of method for shale gas well reservoir fragility overall merit |
CN108009705A (en) * | 2017-11-07 | 2018-05-08 | 中国石油大学(华东) | A kind of shale reservoir compressibility evaluation method based on support vector machines technology |
CN108073786A (en) * | 2017-12-27 | 2018-05-25 | 中国石油化工股份有限公司江汉油田分公司勘探开发研究院 | Shale compressibility evaluation method and device |
CN109407150A (en) * | 2018-12-06 | 2019-03-01 | 中国石油大学(北京) | Based on the petrophysical shale reservoir compressibility means of interpretation of statistics and system |
CN110714753A (en) * | 2019-08-09 | 2020-01-21 | 中国华能集团清洁能源技术研究院有限公司 | Shale reservoir compressibility evaluation method |
CN113090258A (en) * | 2021-05-25 | 2021-07-09 | 中国石油天然气股份有限公司 | Deep shale gas horizontal well fracturing differentiation design method based on logging data |
Non-Patent Citations (1)
Title |
---|
王乔等: "页岩脆性的表征方法及主控因素", 《断块油气田》 * |
Also Published As
Publication number | Publication date |
---|---|
CN113622905B (en) | 2023-05-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107917865B (en) | compact sandstone reservoir multi-parameter permeability prediction method | |
CN111381292B (en) | Logging interpretation method and device for predicting sandstone hydrocarbon-bearing reservoir | |
RU2485553C1 (en) | Method of estimating fracture porosity based on borehole seismic data | |
Ghorbani et al. | Application of Leeb hardness test in prediction of dynamic elastic constants of sedimentary and igneous rocks | |
CN113216945B (en) | Quantitative evaluation method for permeability of tight sandstone reservoir | |
CN109901238A (en) | A kind of High stress zone resistivity correction method based on the experiment of stress difference resistivity | |
Sengupta et al. | On sieving and settling techniques for sand analysis | |
Carreira et al. | Nondestructive evaluation of Corymbia citriodora logs by means of the free transverse vibration test | |
CN113622905B (en) | Shale reservoir brittleness evaluation method based on multi-factor comprehensive analysis | |
CN112945772A (en) | Engineering rock body mechanical property analysis method under water rock circulation | |
CN115356772B (en) | Method for evaluating brittleness of continental facies shale gas reservoir by considering interlayer type | |
Teder et al. | Investigation of the physical-mechanical properties of timber using ultrasound examination | |
CN106680375A (en) | Air-coupled ultrasonic testing method for determining elastic modulus of material | |
CN115391739A (en) | Quantitative calculation method and system for crack permeability | |
CN112082918B (en) | Method, device and equipment for determining porosity | |
Fondjo et al. | Models for predicting the suction of heaving compacted soils using geotechnical properties | |
CN109520461B (en) | Method for determining statistical sample number of series-size rock structural surface roughness samples | |
CN108507923A (en) | Rock core magnetic nuclear resonance analyzer porosity measurement precision is examined and bearing calibration | |
CN116859455A (en) | Complex reservoir anisotropy parameter prediction method and device, electronic equipment and medium | |
CN112526622A (en) | Pseudo-nuclear magnetic echo data calculation method based on imaging logging porosity spectrum | |
CN112798652A (en) | Method for predicting collapsibility of loess by using resistivity | |
Zeigler et al. | A practical method for assigning uncertainty and improving the accuracy of alpha-ejection corrections and eU concentrations in apatite (U-Th)/He chronology | |
Fernández-Serrano et al. | Longitudinal, transverse and ultrasound vibration for the prediction of stiffness using models incorporating features in Pinus sylvestris timber | |
Röshoff et al. | Strategy for a Rock Mechanics Site Descriptive Model. Development and testing of the empirical approach | |
Fernández-Serrano et al. | Modulus of Rupture Prediction in Pinus sylvestris with Multivariate Models Constructed with Resonance, Ultrasound, and Wood Heterogeneity Variables. |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
TA01 | Transfer of patent application right |
Effective date of registration: 20230410 Address after: 430000 Lu Mill Road, Hongshan District, Wuhan, Hubei Province, No. 388 Applicant after: CHINA University OF GEOSCIENCES (WUHAN CITY) Applicant after: CHINA PETROLEUM & CHEMICAL Corp. Address before: 430000 Lu Mill Road, Hongshan District, Wuhan, Hubei Province, No. 388 Applicant before: CHINA University OF GEOSCIENCES (WUHAN CITY) |
|
TA01 | Transfer of patent application right | ||
GR01 | Patent grant | ||
GR01 | Patent grant |