CN106383053A - Engineering mechanical parameter related brittleness index prediction method - Google Patents

Engineering mechanical parameter related brittleness index prediction method Download PDF

Info

Publication number
CN106383053A
CN106383053A CN201610730193.5A CN201610730193A CN106383053A CN 106383053 A CN106383053 A CN 106383053A CN 201610730193 A CN201610730193 A CN 201610730193A CN 106383053 A CN106383053 A CN 106383053A
Authority
CN
China
Prior art keywords
brittleness index
rit
regression
core sample
core
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN201610730193.5A
Other languages
Chinese (zh)
Inventor
周立宏
蒲秀刚
陈长伟
肖敦清
杨飞
韩国猛
马建英
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Research Institute Of Petroleum Exploration & Development Dagang Oil Field Of Cnpc
Original Assignee
Research Institute Of Petroleum Exploration & Development Dagang Oil Field Of Cnpc
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Research Institute Of Petroleum Exploration & Development Dagang Oil Field Of Cnpc filed Critical Research Institute Of Petroleum Exploration & Development Dagang Oil Field Of Cnpc
Priority to CN201610730193.5A priority Critical patent/CN106383053A/en
Publication of CN106383053A publication Critical patent/CN106383053A/en
Pending legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N3/00Investigating strength properties of solid materials by application of mechanical stress
    • G01N3/08Investigating strength properties of solid materials by application of mechanical stress by applying steady tensile or compressive forces
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2203/00Investigating strength properties of solid materials by application of mechanical stress
    • G01N2203/0058Kind of property studied
    • G01N2203/0076Hardness, compressibility or resistance to crushing
    • G01N2203/0087Resistance to crushing

Landscapes

  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)

Abstract

Belonging to the technical field of petroleum and natural gas exploration, the invention relates to an engineering mechanical parameter related brittleness index prediction method. The method includes the steps of: 1) core sampling design and sampling; 2) core sample mineral component testing; 3) testing of core sample Young's modulus and Poisson's ratio; 4) establishment of a mathematical model between rock mineral components and engineering mechanical parameters, and determining the relationship between the mineral components and the rock mechanical parameters; 5) calculation of the engineering mechanical parameters of all samples; and 6) calculation of the rock brittleness index. The method acquires the relationship between the rock mechanical parameters and mineral components and the relationship between the mineral content and the rock brittleness index through mathematical models. The method has the advantages of short test period, simple sampling operation and cost saving, can guide the fracturability evaluation of the rock stratum more efficiently, and provides technical support for oil and gas exploration and development.

Description

Brittleness index prediction method related to engineering mechanical parameters
Technical Field
The invention relates to the technical field of petroleum and natural gas exploration, in particular to a brittleness index prediction method related to engineering mechanical parameters.
Background
Rock brittleness is an important parameter in the evaluation of fracability during the development of tight oil. Engineering mechanical parameters representing rock brittleness are mainly Young modulus and Poisson ratio, and as the Young modulus and the Poisson ratio of a tested stratum must depend on core data and are expensive, generally, less data are obtained. In early studies, the brittleness index of rock samples was often calculated using the level of the main brittle mineral (quartz) content, or the brittle mineral combination (quartz feldspar, dolomite, calcite, etc.). On one hand, as the content of the minerals and the mineral combination and the engineering mechanical parameters do not establish necessary correlation, no good corresponding relation exists between the calculated brittleness index and the rock fracturing performance; on the other hand, engineering mechanical parameters are different among different minerals, and simple accumulation results among brittle minerals cannot be effectively applied to fracture evaluation of rocks.
Disclosure of Invention
In order to solve the problem of determining the brittleness index in the fracturing evaluation process of the compact oil and gas reservoir, the invention aims to provide a brittleness index prediction method which is related to engineering mechanics and is convenient and fast to calculate.
In order to realize the purpose of the invention, the invention adopts the following technical scheme:
a brittleness index prediction method related to engineering mechanical parameters is characterized by comprising the following steps:
1) sampling the core to obtain a core sample;
2) testing the mineral composition of the core sample;
3) testing the Young modulus and Poisson ratio of the core sample;
4) calculating the brittleness index of the core sample according to the Young modulus and the Poisson ratio;
5) establishing a mathematical model between the mineral composition and the brittleness index, and determining the relationship between the mineral composition and the brittleness index;
6) the brittleness index was calculated using the mineral composition.
Wherein, in the step 1), core samples used for composition test are uniformly distributed and taken out in the core section, and the core samples can be used for testing rock composition and engineering mechanical parameters.
Wherein in step 2), the composition of the core sample and its relative content are determined using X-ray diffraction (XRD).
In the step 3), a rock mechanics triaxial experiment is performed on the core sample to obtain a stress-strain curve of the core sample, so that the Young modulus and the Poisson ratio of the core sample are obtained.
In the step 4), calculating the brittleness index of the core sample according to the Young modulus and the Poisson ratio, wherein the calculation method comprises the following steps: according to the normalization results of the Young modulus and the Poisson ratio, defining the root mean square of the Young modulus and the Poisson ratio as a brittleness index;
YBrit=(yi-ymin)/(ymax-ymin)×100
BBrit=(bmax-bi)/(bmax-bmin)×100
B r i t = Y B r i t 2 + B B r i t 2 2
in the formula: y isBritNormalized young's modulus; y isiAs a measure of Young's modulus, ymaxIs the maximum value of Young's modulus, yminMinimum value of Young's modulus, BBritA normalized Poisson's ratio; biIs a measure of poisson's ratio; bmaxIs the maximum poisson ratio; bminIs the minimum poisson ratio; b isritIs a brittleness index and is dimensionless.
Wherein, in step 5), a mathematical model between the mineral composition and the brittleness index is established, and the relationship between the mineral composition and the brittleness index is determined. And, step 5) comprises the following substeps:
5a) the model assumes that: assuming that m core samples exist in the model, and each core sample contains n mineral components; the composition of the mineral constituents of the core sample is regarded as RnSample space, the mineral composition of the m core samples then forming a point set, the mineral composition matrix A of the sample spacemn
A m n = a 11 a 12 ... a 1 n a 21 a 22 ... a 2 n . . . . . . . . . . . . a m 1 a m 2 ... a m n
At RpIn space, a one-dimensional subspace F passing through the origin can be found1Represents the one-dimensional subspace F1Is directed by a unit vector u1∈RnTo define; at F1Upper i-th brittleness index valueCan be observed from the value Mi∈RnMapping to u from corresponding points1To obtain F1On the upper partThe coordinates of (a) are given by:
B rit i = M i T u 1 | | u 1 | | = M i T u 1
defining an optimal line F using a "least squares method1: find u1∈RPTo minimize the following:
Σ i = 1 n | | M i - B rit i | | 2
according to the Pythagorean theorem:minimization of the problem Equivalent to maximization
The problem becomes to find in | | | u1Maximization under the constraint of 1 | | | >Obtaining:
B rit 1 B rit 2 . . . B rit n = M 1 T u 1 M 2 T u 1 M m T u 1 = Au 1
this problem is re-expressed as: find when | | u1U when 1 | |)1∈RnTo maximize the quadratic term (Au)1)TAu1) Or
m a x u 1 T u 1 u 1 T M T Mu 1
Mixing mineral component MiTo a brittleness index BritIs treated as a multivariate linear regression model
Brit=β1m12m2+…+βnmn
Formula (III) β1β2…βnIs a regression coefficient
In the formula b1b2…bnAre estimates of the regression coefficients.
5b) Parameter estimation
B r i t = B rit 1 B rit 2 . . . B rit n M = m 11 m 12 ... m 1 n m 21 m 22 ... m 2 n . . . . . . . . . . . . m m 1 m m 2 ... m m n B = b 1 b 2 . . . b m
B=(M′M)-1X′Y
5c) Hypothesis testing
Correlation coefficient-measure the degree of coincidence of regression equation and original data:
sum of the squares of total deviations SS:referred to as BritDispersion of (a). All areIs called BritSum of the squares of the total deviations SS:
SSresidthe sum of the squares of the residuals reflects the experimental valuesWith values calculated according to the regression equationThe smaller the total deviation, the better the regression effect; SSRegressionThe regression sum of squares reflects M and BritB is caused by the linear relationship ofritThe larger the magnitude of the change, the better the regression effect;
⑤r2correlation test
SSRegressionThe larger, M and BritThe more important the regression relationship of r2The closer to 1, the other SSresidThe smaller the linear relationship, the better;
sixthly, significance test of the regression equation:
to test whether there is a significant linear relationship between the dependent and independent variables in the model, a statistic is constructed:
for a given significance level α, determine the rejection region F > Fα(k, n-k-1). And calculating a statistic value, and judging whether to reject the original hypothesis.
Wherein in step 6), b is determined according to step 5)1,b2…bnAnd calculating the brittleness index of the core sample.
Compared with the prior art, the brittleness index calculation method related to engineering mechanical parameters has the following beneficial effects:
1) the method solves the problems encountered in the engineering fracturing process from the angle of engineering mechanical tests. The fracturing performance of the core is directly related to the engineering mechanical property, and the test result can be directly used for making a fracturing scheme and is an intuitive expression of the fracturing performance. The problem is solved from the perspective of engineering mechanics, and the practicability of the invention is increased.
2) And searching for determinants of engineering mechanical parameters according to the mineral component content of the rock sample. The engineering mechanical properties of rock are determined by the mineral composition and relative content thereof. The content of minerals with different engineering mechanical properties contributes differently to the brittleness of the rock. The mineral with high Young modulus and small Poisson ratio has higher brittleness index. Through the establishment of the mathematical model, the relation between the engineering mechanical property and the mineral content in the research area can be found out.
3) Engineering mechanical tests are more expensive than analysis of mineral composition in terms of sample requirements and time periods. The invention utilizes a small amount of engineering mechanical tests and combines with the conventional test project achievements to efficiently and quickly obtain the brittleness index of the research area, thereby saving the cost on the basis of ensuring the data validity.
Drawings
FIG. 1 is a model map of the mapping between mineral content and friability index.
FIG. 2 photo of rock mechanics triaxial pressure test core sample before fracturing.
FIG. 3 is a photo of a fractured rock mechanical triaxial pressure test core sample.
FIG. 4 rock brittleness index calculated by different methods for a well.
Detailed Description
The brittleness index calculation method of the present invention is further described below with reference to specific examples to help those skilled in the art to more completely, accurately and deeply understand the inventive concept and technical scheme of the present invention; it is to be understood that the description in specific embodiments is intended to be illustrative, and not restrictive, of the scope of the invention, which is defined by the appended claims.
Example 1
The brittleness index prediction method of the embodiment comprises the following steps:
1. designing and sampling a core sample;
the sample is ensured to be consistent in depth according to sampling requirements of XRD whole rock X-ray diffraction analysis and engineering mechanical tests.
2. Testing mineral components of the core sample;
XRD diffraction test analysis of the sample shows that the rock mineral composition in the compact rock stratum is various, and the main components are quartz, calcite, analcite, clay, dolomite, feldspar, pyrite and the like. The mineral composition test data are shown in table 1:
table 1 XRD analysis of mineral content
Sample numbering Depth (m) Horizon Quartz crystal Calcite Analcime (analcime) Clay Dolomite Feldspar Pyrite
1 2939.53 Ek2 1 13 2 48 18 2 15 2
2 2953.31 Ek2 1 16 3 27 13 18 20 3
3 2955.57 Ek2 1 16 5 29 14 8 21 7
4 2966.03 Ek2 1 14 18 9 7 39 12 1
5 2974.07 Ek2 1 10 5 14 11 52 7 1
6 3025.10 Ek2 2 17 13 24 11 6 28 1
7 3032.09 Ek2 2 16 14 14 12 28 15 1
8 3032.26 Ek2 2 14 5 16 10 43 11 1
9 3059.18 Ek2 2 14 8 28 13 23 13 1
10 3187.22 Ek2 3 15 6 6 20 26 26 1
11 3197.13 Ek2 3 21 8 6 20 18 18 9
12 3209.06 Ek2 3 18 7 0 25 11 38 1
13 3234.94 Ek2 3 19 0 9 15 41 16 0
14 3239.22 Ek2 3 14 0 26 33 13 13 1
15 3297.81 Ek2 3 14 10 8 10 48 9 1
16 3310.64 Ek2 4 28 14 21 12 5 19 1
17 3315.08 Ek2 4 14 13 17 15 27 13 1
18 3318.20 Ek2 4 36 15 6 9 9 25 0
19 3354.35 Ek2 4 40 8 0 3 0 49 0
20 3360.54 Ek2 4 39 14 0 2 0 45 0
21 3376.96 Ek2 4 37 12 0 3 0 48 0
22 3385.59 Ek2 4 43 4 0 10 0 0 0
3. Testing engineering mechanical parameters (Young modulus and Poisson ratio) of a core sample;
the poisson ratio and the rock compressive strength of the rock are obtained by measuring the longitudinal and transverse deformation of the rock test piece in a regular shape under the action of triaxial pressure.
Using a servo-controlled rock mechanics triaxial experimental setup, the core sample as shown in fig. 2 was placed in the pressure chamber of a triaxial autoclave, a certain lateral pressure (σ _3 ═ 25MPa) was applied, and then a vertical pressure (σ _1) was applied until the core sample was destroyed as shown in fig. 3. The instrument simultaneously measures axial and radial stress-strain data to respectively obtain axial and radial stress-strain relations. The Young's modulus and Poisson's ratio of each sample were obtained from the stress-strain relationship, as shown in Table 2.
TABLE 2 Young's modulus and Poisson's ratio for core samples
Sample numbering Depth before homing (m) Horizon Confining pressure (MPa) Poisson ratio Young's modulus (GPa)
1 2939.53 Ek2 1 25 0.379 15.62
2 2953.31 Ek2 1 25 0.417 13.96
3 2955.57 Ek2 1 25 0.309 19.33
4 2966.03 Ek2 1 25 0.313 23.3
5 2974.07 Ek2 1 25 0.265 23.62
6 3025.10 Ek2 2 25 0.326 21.85
7 3032.09 Ek2 2 25 0.193 27.97
8 3032.26 Ek2 2 25 0.279 30.82
9 3059.18 Ek2 2 25 0.274 26.19
10 3187.22 Ek2 3 25 0.341 14.59
11 3197.13 Ek2 3 25 0.226 22.15
12 3209.06 Ek2 3 25 0.402 9.64
13 3234.94 Ek2 3 25 0.271 28.1
14 3239.22 Ek2 3 25 0.277 13.72
15 3297.81 Ek2 3 25 0.184 21.61
16 3310.64 Ek2 4 25 0.259 31.57
17 3315.08 Ek2 4 25 0.213 25.78
18 3318.20 Ek2 4 25 0.198 43.72
19 3354.35 Ek2 4 25 0.206 22.14
20 3360.54 Ek2 4 25 0.197 26.78
21 3376.96 Ek2 4 25 0.238 25.92
22 3385.59 Ek2 4 25 0.311 26.89
3. Calculating the brittleness index of the sample according to the Young modulus and the Poisson ratio;
using the calculation formula of brittleness index:
YBrit=(yi-ymin)/(ymax-ynin)×100
BBrit=(bmax-bi)/(bmax-bmin)×100
B r i t = Y B r i t 2 + B B r i t 2 2
in the formula: y isBritIs the normalized young's modulus; y ismaxIs the maximum value of Young's modulus; y isminIs the Young's modulus minimum; b isBritIs the normalized Poisson's ratio; bmaxIs the maximum poisson ratio; bminIs the minimum poisson ratio; b isritIs a brittleness index and is dimensionless.
From table 2, it can be found that:
ymin=9.64;ymax=43.72;bmax=0.417;bmin=0.184
determining brittleness index Y of samples No. 1-22Brit1As shown in table 3.
Core sample brittleness index of table 31-22
Sample numbering Depth before homing (m) YBrit1 YBrit2 Sample numbering Depth before homing (m) YBrit1
1 2939.53 16.93 7.51 12 3209.06 3.22
2 2953.31 6.34 0.00 13 3234.94 58.41
3 2955.57 37.39 21.83 14 3239.22 36.03
4 2966.03 42.36 24.00 15 3297.81 67.56
5 2974.07 53.13 32.23 16 3310.64 66.08
6 3025.10 37.44 20.80 17 3315.08 67.46
7 3032.09 74.96 47.34 18 3318.20 97.00
8 3032.26 60.69 35.04 19 3354.35 63.62
9 3059.18 54.97 32.56 20 3360.54 72.36
10 3187.22 23.57 13.11 21 3376.96 62.30
11 3197.13 59.34 37.68 22 3385.59 48.05
4. Mathematical model between engineering mechanical parameters and brittleness index
The brittle mineral content is an important factor affecting the brittleness and fracturability of fine-grained sedimentary rocks. The higher the content of quartz and other brittle minerals is, the stronger the rock brittleness is, induced cracks are easy to form in the hydraulic fracturing process, and a complex crack network is formed, thereby being beneficial to the exploitation of compact oil. Brittleness index in the narrow sense Brit1Quartz mineral content is used to represent:
generalized brittleness index Brit2Comprising quartz, feldspar, dolomite, calcite, pyrite and analcite:
in the classical brittleness index and the generalized brittleness index, coefficients of minerals such as quartz are set to 1, and in order to make a model comparable, coefficients of quartz are set to 1 first, and then coefficients of other minerals are calculated respectively in the process of performing regression analysis.
The mathematical model established by the step 5 in the invention content part of the specification
In the formula b1b2…bnAre estimates of the regression coefficients.
Wherein,
B r i t = B rit 1 B rit 2 . . . B rit n M = m 11 m 12 ... m 1 n m 21 m 22 ... m 2 n . . . . . . . . . . . . m m 1 m m 2 ... m m n B = b 1 b 2 . . . b m
B=(M′M)-1X′Y
and establishing a model suitable for local brittleness index evaluation by using the mathematical model. The mineral components of samples nos. 1 to 22 and Brit3 were introduced into the SPSS database using SPSS multivariate statistical software, and subjected to multivariate linear regression to obtain the following models:
by comparing the mineral composition with the brittleness index Brit3Multiple linear regression analysis of (a) to (b) to obtain an evaluation model of friability index for the study area (table 4). The DW value of a regression model of the regression equation is 1.984, is close to 2, and shows that the rock mineral components are independent; r value of the regression model is0.903,R2A value of 0.815, indicating a brittleness index Brit3There is a good linear relationship with the mineral composition, and the model passes the significance hypothesis. The regression equation between brittleness index and mineral composition in the study area is:
Brit3quartz +0.629 dolomite +0.521 feldspar +0.25 calcite +0.204 pyrite +0.18 analcite +0.021 clay.
5. And calculating the rock brittleness index by using the mineral components.
The total 1219 total rock X-ray diffraction analysis data of a certain well core-taking section are calculated according to the classical brittleness index, the generalized brittleness index and the brittleness index calculation model established in the step 5, and the classical brittleness index, the generalized brittleness index and the regression equation brittleness index are respectively calculated, and the result is shown in FIG. 4.
The classical brittleness index is distributed between 0 and 0.5, and the average is 0.2; the generalized brittleness index is distributed between 0 and 1, and the average value is 0.85; the brittleness index of the time is distributed between 0 and 1, and the average value is 0.50. The classical brittleness index and the generalized brittleness index have certain defects in use, the classical brittleness index is small, the generalized brittleness value is large, and the classical brittleness index and the generalized brittleness index are often ignored or overestimated for the compressibility of the alloy when the engineering property of the alloy is judged. The brittleness index calculated using the regression equation can better describe the brittleness characteristics of fine grained sedimentary rock, and is of great significance for the preferred engineered dessert.
TABLE 4 model for evaluation of brittleness index in research area
Model summaryo,d
a. Prediction variables: clay, calcite, pyrite, dolomite, analcite, feldspar
b. For regression through the origin (no intercept model), the R-side can measure the proportion of variability in the dependent variable near the origin (explained by the regression). For models containing an intercept, this cannot be compared to the R-square.
c. Dependent variable: index of brittleness
d. Linear regression through the origin
Coefficient of performancea,b
a. Dependent variable: index of brittleness
b. Linear regression through the origin
It is obvious to those skilled in the art that the present invention is not limited to the above embodiments, and it is within the scope of the present invention to adopt various insubstantial modifications of the method concept and technical scheme of the present invention, or to directly apply the concept and technical scheme of the present invention to other occasions without modification.

Claims (7)

1. A brittleness index prediction method related to engineering mechanical parameters is characterized by comprising the following steps:
1) sampling the core to obtain a core sample;
2) testing the mineral composition of the core sample;
3) testing the Young modulus and Poisson ratio of the core sample;
4) calculating the brittleness index of the core sample according to the Young modulus and the Poisson ratio;
5) establishing a mathematical model between the mineral composition and the brittleness index, and determining the relationship between the mineral composition and the brittleness index;
6) the brittleness index was calculated using the mineral composition.
2. The brittleness index prediction method of claim 1, wherein: in the step 1), core samples used for composition testing are uniformly distributed and taken out in the core section, and the core samples can be used for testing rock compositions and engineering mechanical parameters.
3. The brittleness index prediction method of claim 1, wherein: in step 2), the composition of the core sample and its relative content were determined using X-ray diffraction.
4. The brittleness index prediction method of claim 1, wherein: in the step 3), a rock mechanics triaxial experiment is carried out on the core sample to obtain a stress-strain curve of the core sample, so that the Young modulus and the Poisson ratio of the core sample are obtained.
5. The brittleness index prediction method of claim 1, wherein: in the step 4), calculating the brittleness index of the core sample according to the Young modulus and the Poisson ratio, wherein the calculation method comprises the following steps: according to the normalization results of the Young modulus and the Poisson ratio, defining the root mean square of the Young modulus and the Poisson ratio as a brittleness index;
YBrit=(yi-ymin)/(ymax-ymin)×100
B B r i t = ( b max - b i ) / ( b max - b min ) × 100 B r i t = Y B r i t 2 + B B r i t 2 2
in the formula: y isBritNormalized young's modulus; y isiAs a measure of Young's modulus, ymaxIs the maximum value of Young's modulus, yminMinimum value of Young's modulus, BBritA normalized Poisson's ratio; biIs a measure of poisson's ratio; bmaxIs the maximum poisson ratio; bminIs the minimum poisson ratio; b isritIs a brittleness index and is dimensionless.
6. The brittleness index prediction method of claim 1, wherein: step 5) comprises the following substeps:
5a) the model assumes that: assuming that m core samples exist in the model, and each core sample contains n mineral components; the composition of the mineral constituents of the core sample is regarded as RnSample space, the mineral composition of the m core samples then forming a point set, the mineral composition matrix A of the sample spacemn
A m n = a 11 a 12 ... a 1 n a 21 a 22 ... a 2 n . . . . . . . . . . . . a m 1 a m 2 ... a m n
At RpIn space, a one-dimensional subspace F passing through the origin can be found1Represents the one-dimensional subspace F1Is directed by a unit vector u1∈RnTo define; at F1Upper i-th brittleness index valueCan be observed from the value Mi∈RnMapping to u from corresponding points1To obtain F1On the upper partThe coordinates of (a) are given by:
B rit i = M i T u 1 | | u 1 | | = M i T u 1
defining an optimal line F using a "least squares method1: find u1∈RPTo minimize the following:
Σ i = 1 n | | M i - B rit i | | 2
according to the Pythagorean theorem:minimization of the problem Equivalent to maximization
The problem becomes to find in | | | u1Maximization under the constraint of 1 | | | >U of (a)1∈RnObtaining:
B rit 1 B rit 2 . . . B rit m = M 1 T u 1 M 2 T u 1 M m T u 1 = Au 1
this problem is re-expressed as: find when | | u1U when 1 | |)1∈RnTo maximize the quadratic term (Au)1)TAu1) Or
max u 1 T u 1 u 1 T M T Mu 1
Mixing mineral component MiTo a brittleness index BritIs treated as a multivariate linear regression model
Brit=β1m12m2+…+βnmn
Formula (III) β1β2…βnIs a regression coefficient
In the formula b1b2…bnAre estimates of the regression coefficients.
5b) Parameter estimation
B r i t = B rit 1 B rit 2 . . . B rit n M = m 11 m 12 ... m 1 n m 21 m 22 ... m 2 n . . . . . . . . . . . . m m 1 m m 2 ... m m n B = b 1 b 2 . . . b m
B=(M′M)-1X′Y
5c) Hypothesis testing
Correlation coefficient-measure the degree of coincidence of the regression equation with the original data:
sum of squares of total deviationsReferred to as BritDispersion of (a). All areIs called BritSum of the squares of the total deviations SS:
SSresidthe sum of the squares of the residuals reflects the experimental valuesWith values calculated according to the regression equationThe smaller the total deviation, the better the regression effect; SSRegressionThe regression sum of squares reflects M and BritB is caused by the linear relationship ofritThe larger the magnitude of the change, the better the regression effect;
②r2correlation test
SSRegressionThe larger the size of the tube is,m and BritThe more important the regression relationship of r2The closer to 1, the other SSresidThe smaller the linear relationship, the better;
③ the significance test of the regression equation:
to test whether there is a significant linear relationship between the dependent and independent variables in the model, a statistic is constructed:
for a given significance level α, determine the rejection region F > Fα(k, n-k-1). And calculating a statistic value, and judging whether to reject the original hypothesis.
7. The brittleness index prediction method of claim 1, wherein: in step 6), according to b determined in step 5)1,b2…bnAnd calculating the brittleness index of the core sample.
CN201610730193.5A 2016-08-29 2016-08-29 Engineering mechanical parameter related brittleness index prediction method Pending CN106383053A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610730193.5A CN106383053A (en) 2016-08-29 2016-08-29 Engineering mechanical parameter related brittleness index prediction method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610730193.5A CN106383053A (en) 2016-08-29 2016-08-29 Engineering mechanical parameter related brittleness index prediction method

Publications (1)

Publication Number Publication Date
CN106383053A true CN106383053A (en) 2017-02-08

Family

ID=57917210

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610730193.5A Pending CN106383053A (en) 2016-08-29 2016-08-29 Engineering mechanical parameter related brittleness index prediction method

Country Status (1)

Country Link
CN (1) CN106383053A (en)

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108519281A (en) * 2018-02-27 2018-09-11 中国石油天然气股份有限公司 Method, device and system for determining brittleness index of rock
CN108535112A (en) * 2017-03-03 2018-09-14 中国石油化工股份有限公司 A kind of experiment analytical method for shale samples compressibility research
CN111238939A (en) * 2020-03-17 2020-06-05 中联煤层气有限责任公司 Method and device for determining rock brittleness index
CN113030440A (en) * 2021-03-22 2021-06-25 西南石油大学 Method for predicting fracture brittleness index of shale in hydraulic fracturing process of shale
CN115266785A (en) * 2022-06-27 2022-11-01 西南石油大学 Shale oil reservoir brittleness evaluation method based on rock mechanical properties
US12050297B2 (en) 2020-09-11 2024-07-30 Saudi Arabian Oil Company Method and system for determining energy-based brittleness

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104406849A (en) * 2014-11-21 2015-03-11 中国石油天然气股份有限公司 Method and device for predicting brittleness of reservoir rock
CN104677734A (en) * 2015-02-12 2015-06-03 中国石油大学(华东) Multi-scale information integration tight formation brittleness index measuring method
CN105221141A (en) * 2014-06-23 2016-01-06 中国石油化工股份有限公司 A kind of mud shale brittleness index Forecasting Methodology
CN105221140A (en) * 2014-06-20 2016-01-06 中国石油化工股份有限公司 A kind ofly determine that shale formation can the method for pressure break sex index
CN105277671A (en) * 2014-06-18 2016-01-27 中国石油化工股份有限公司 Method used for determining shale formation brittleness index

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105277671A (en) * 2014-06-18 2016-01-27 中国石油化工股份有限公司 Method used for determining shale formation brittleness index
CN105221140A (en) * 2014-06-20 2016-01-06 中国石油化工股份有限公司 A kind ofly determine that shale formation can the method for pressure break sex index
CN105221141A (en) * 2014-06-23 2016-01-06 中国石油化工股份有限公司 A kind of mud shale brittleness index Forecasting Methodology
CN104406849A (en) * 2014-11-21 2015-03-11 中国石油天然气股份有限公司 Method and device for predicting brittleness of reservoir rock
CN104677734A (en) * 2015-02-12 2015-06-03 中国石油大学(华东) Multi-scale information integration tight formation brittleness index measuring method

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
CHANGLIANG FANG ET AL.: "Influence Factors of Fracability in Nonmarine Shale", 《SPE INTERNATIONAL》 *
徐赣川等: "基于岩石物理实验的页岩脆性测井评价方法", 《天然气工业》 *
范卓颖等: "致密地层岩石脆性指数的测井优化建模", 《石油学报》 *

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108535112A (en) * 2017-03-03 2018-09-14 中国石油化工股份有限公司 A kind of experiment analytical method for shale samples compressibility research
CN108535112B (en) * 2017-03-03 2021-11-30 中国石油化工股份有限公司 Experimental analysis method for shale sample compressibility research
CN108519281A (en) * 2018-02-27 2018-09-11 中国石油天然气股份有限公司 Method, device and system for determining brittleness index of rock
CN108519281B (en) * 2018-02-27 2020-08-11 中国石油天然气股份有限公司 Method, device and system for determining brittleness index of rock
CN111238939A (en) * 2020-03-17 2020-06-05 中联煤层气有限责任公司 Method and device for determining rock brittleness index
US12050297B2 (en) 2020-09-11 2024-07-30 Saudi Arabian Oil Company Method and system for determining energy-based brittleness
CN113030440A (en) * 2021-03-22 2021-06-25 西南石油大学 Method for predicting fracture brittleness index of shale in hydraulic fracturing process of shale
CN115266785A (en) * 2022-06-27 2022-11-01 西南石油大学 Shale oil reservoir brittleness evaluation method based on rock mechanical properties
CN115266785B (en) * 2022-06-27 2024-07-05 西南石油大学 Shale oil reservoir brittleness evaluation method based on rock mechanical properties

Similar Documents

Publication Publication Date Title
CN106383053A (en) Engineering mechanical parameter related brittleness index prediction method
CN105445440B (en) A kind of shale compressibility evaluation method based on landwaste microscopic feature
US10571605B2 (en) Method and system for modeling and simulating a fractured reservoir
CN103278436B (en) Quantitative characterization method of low penetration double-medium sandstone oil reservoir microscopic aperture structure
Clarkson et al. Nanopores to megafractures: Current challenges and methods for shale gas reservoir and hydraulic fracture characterization
Yanjun et al. Division of coalbed methane desorption stages and its significance
Clarkson et al. A new low-permeability reservoir core analysis method based on rate-transient analysis theory
US20230220755A1 (en) Method and apparatus for predicting oil and gas yields in in-situ oil shale exploitation
CN103198363B (en) Reservoir gas production amount prediction method and device based on CT pore analysis
CN108804849B (en) Rock mechanical parameter evaluation method based on structural complexity
US20200240251A1 (en) Method for characterizing complexity of rock fracture based on fractal dimension and device thereof
CN104007482A (en) Shale rock physical model method based on anisotropic effective field
CN110108856B (en) Method for predicting maturity of organic matters of hydrocarbon source rock under overpressure background
CN108169179B (en) Method for determining effective conditions for evaluating biogenic source constitution of hydrocarbon source rock by using n-alkanes
CN109632594B (en) Method for characterizing pore throat characteristics of tight reservoir on basis of high-pressure mercury intrusion multi-scale
Su et al. Quantitative study on hydrocarbon expulsion mechanism based on micro-fracture
CN112129802B (en) Quantitative analysis method for pore volume increment of hydrated shale in different scales
CN110245686A (en) A kind of lithology method for quickly identifying calculating quartzy percentage contents
CN110568160A (en) Comprehensive evaluation method and device for oil and gas reservoir rock
CN114970235B (en) Geomechanical heterogeneity-anisotropy modeling method for fractured reservoir
Esatyana et al. Nanoindentation of shale cuttings and its application to core measurements
Lu et al. Gas flow characteristics in shale fractures after supercritical CO2 soaking
Song et al. Shale softening degree and rate induced by fracturing fluid under THMC coupling condition
Li et al. Combining macro and micro experiments to reveal the real-time evolution of permeability of shale
CN106501149A (en) A kind of porosity calculation method that is analyzed based on coal petrography industrial components

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20170208

WD01 Invention patent application deemed withdrawn after publication