CN112307601A - Complex reservoir fracturing property evaluation method - Google Patents
Complex reservoir fracturing property evaluation method Download PDFInfo
- Publication number
- CN112307601A CN112307601A CN202011087617.3A CN202011087617A CN112307601A CN 112307601 A CN112307601 A CN 112307601A CN 202011087617 A CN202011087617 A CN 202011087617A CN 112307601 A CN112307601 A CN 112307601A
- Authority
- CN
- China
- Prior art keywords
- core sample
- rock
- fracture
- reservoir
- core
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000011156 evaluation Methods 0.000 title claims abstract description 24
- 239000011435 rock Substances 0.000 claims abstract description 57
- 238000002474 experimental method Methods 0.000 claims abstract description 33
- 238000000034 method Methods 0.000 claims abstract description 21
- 238000012360 testing method Methods 0.000 claims abstract description 14
- 238000012937 correction Methods 0.000 claims abstract description 12
- 238000013210 evaluation model Methods 0.000 claims abstract description 8
- 230000006378 damage Effects 0.000 claims abstract description 6
- 238000004364 calculation method Methods 0.000 claims description 28
- 238000001228 spectrum Methods 0.000 claims description 16
- 150000001875 compounds Chemical class 0.000 claims description 6
- 238000011160 research Methods 0.000 claims description 5
- 238000010219 correlation analysis Methods 0.000 claims description 3
- 238000011161 development Methods 0.000 abstract description 10
- 238000000611 regression analysis Methods 0.000 abstract description 7
- 230000009286 beneficial effect Effects 0.000 abstract description 4
- 206010017076 Fracture Diseases 0.000 description 25
- BVKZGUZCCUSVTD-UHFFFAOYSA-L Carbonate Chemical compound [O-]C([O-])=O BVKZGUZCCUSVTD-UHFFFAOYSA-L 0.000 description 5
- 239000011148 porous material Substances 0.000 description 2
- 238000007792 addition Methods 0.000 description 1
- 230000004075 alteration Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- BVKZGUZCCUSVTD-UHFFFAOYSA-N carbonic acid Chemical compound OC(O)=O BVKZGUZCCUSVTD-UHFFFAOYSA-N 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000012417 linear regression Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000000704 physical effect Effects 0.000 description 1
- 239000004575 stone Substances 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 239000013589 supplement Substances 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B43/00—Methods or apparatus for obtaining oil, gas, water, soluble or meltable materials or a slurry of minerals from wells
- E21B43/25—Methods for stimulating production
- E21B43/26—Methods for stimulating production by forming crevices or fractures
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Mining & Mineral Resources (AREA)
- Life Sciences & Earth Sciences (AREA)
- Theoretical Computer Science (AREA)
- Geology (AREA)
- Geochemistry & Mineralogy (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Fluid Mechanics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Environmental & Geological Engineering (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
The invention discloses a method for evaluating the fracturing property of a complex reservoir, and particularly relates to the field of oil and gas exploration and development of the complex reservoir. The method comprises the steps of selecting a plurality of core samples in a complex reservoir, sequentially performing a ground stress test experiment, a rock triaxial mechanics experiment and a rock destruction experiment on one part of the core samples to determine the horizontal ground stress difference, brittleness and residual stress of the core samples, performing a Brazilian fracture experiment on the other part of the core samples to determine the fracture toughness and tensile strength of the core samples, determining the fitting degree between each parameter and the fracture pressure of the rock samples based on regression analysis, calculating the weighting coefficient of each parameter, calculating a fracture correction factor by using logging information, integrating each measured parameter and the fracture correction factor of the rock samples, establishing a complex reservoir fracturability evaluation model, and evaluating the fracturability of the complex reservoir. The method improves the accuracy of fracturing evaluation of the complex reservoir, and is beneficial to efficient development of complex oil and gas resources.
Description
Technical Field
The invention relates to the field of complex reservoir oil and gas exploration and development, in particular to a complex reservoir fracturing property evaluation method.
Background
With the continuous increase of the demand of oil and gas resources and the development of technologies, unconventional oil and gas resources are taken as important supplements of conventional oil and gas resources and are widely concerned by countries in the world. Because the unconventional development mode has difficulty in effectively developing a dense reservoir which is deeply buried and has poor physical properties, the formation of a complex fracture network by hydraulic fracturing to improve the flow conductivity has become a main means for the commercial development of the unconventional reservoir nowadays. The reservoir compressibility evaluation optimizes the fracturing interval by assessing the difficulty of reservoir reformation, and provides a reference for establishing an efficient and economical multistage pressure scheme. The method for evaluating the fracturing performance of the reservoir, which is commonly used at the present stage, is developed based on rock brittleness, namely the higher the rock brittleness of the reservoir is, the stronger the fracturing performance of the reservoir is.
On the other hand, complex reservoirs are widely distributed globally, with the exploratory reserves of carbonate reservoirs occupying half of the total reserves, and most are large-scale high-production oil and gas fields. However, the carbonate pore cavern is developed, strong in heterogeneity and hard in texture, and the difficulty of reservoir fracturing in a specific construction process is difficult to accurately evaluate by adopting a single brittleness index of the existing reservoir compressibility evaluation, so that certain obstacles are brought to engineering development. Meanwhile, more and more researchers propose that the stratum with high brittleness coefficient is not necessarily an excellent fracturable interval, and factors such as natural fractures, diagenesis, crustal stress and the like are also important evaluation parameters of the fracturability of the reservoir.
Therefore, accurate evaluation of the fracturing performance of the complex reservoir is difficult by adopting the existing method, and a fracturing performance evaluation method suitable for the complex reservoir needs to be provided.
Disclosure of Invention
The invention provides a method for evaluating the fracturing property of a complex reservoir, which aims at solving the problem that the fracturing property of the complex reservoir is difficult to accurately evaluate at the present stage, and is particularly suitable for carbonate reservoirs. The method improves the rationality of fracturing evaluation of the complex reservoir and is beneficial to efficient development of complex oil and gas resources.
The invention adopts the following technical scheme:
a method for evaluating the fracturing property of a complex reservoir specifically comprises the following steps:
step 1, selecting a plurality of core samples in a complex reservoir based on the complex reservoir in a research block;
step 2, performing an earth stress test experiment on the core sample, measuring the maximum earth stress and the minimum earth stress of the core sample, and calculating the horizontal earth stress difference and the horizontal earth stress difference influence factor of the core sample;
step 3, performing a rock triaxial mechanical experiment on the core sample by using a triaxial mechanical test system, measuring the Young modulus and Poisson ratio of the core sample, and calculating the brittleness of the core sample;
step 4, performing a rock destruction experiment on the core sample by using a triaxial mechanical testing system, measuring the residual strength and the fracture pressure of the core sample, and calculating a normalized residual strength influence factor of the core sample;
step 5, performing Brazilian splitting experiments on the core sample, measuring the fracture toughness and the tensile strength of the core sample, and calculating fracture toughness influence factors and tensile strength influence factors of the core sample;
step 6, performing single correlation analysis on the horizontal ground stress difference, brittleness, fracture toughness, tensile strength and residual strength of the rock core sample and the fracture pressure of the rock sample respectively, and determining the weighting coefficient of each parameter according to the fitting degree;
step 7, acquiring the normalized porosity and the normalized porosity spectrum mean square error of the target well according to the logging information in the research block, and calculating a fracture correction factor;
and 8, establishing a target reservoir fracturing evaluation model according to the weighting coefficients and the fracture correction factors of the parameters of the core sample, calculating the target reservoir fracturing evaluation parameters, and determining the fracturing of the target reservoir.
Preferably, in the step 1, a part of the core samples are used for performing a ground stress test experiment, a rock triaxial mechanical experiment and a rock failure experiment in sequence, and the other part of the core samples are used for performing a brazilian splitting experiment.
Preferably, in the step 2, the horizontal ground stress difference of the core sample is the difference between the maximum ground stress and the minimum ground stress of the rock sample;
horizontal ground stress difference influence factor K of core samplehThe calculation formula is as follows:
in the formula, σHRepresents the maximum ground stress of the core sample, in MPa; sigmahThe minimum stress of the core sample is expressed in MPa.
Preferably, in step 3, the brittleness BI of the core sample is calculated according to the following formula:
BI=(μ+E)/2 (2)
where μ represents the normalized poisson's ratio and E represents the normalized young's modulus, the calculation formula is:
in the formula, mu*Representing the measured poisson ratio of the core sample; mu.smaxRepresenting the maximum value of the Poisson's ratio of the core sample; mu.sminRepresenting the minimum value of the Poisson's ratio of the core sample; e*The measured Young modulus of the core sample is expressed in MPa; emaxThe maximum value of the Young modulus of the core sample is expressed in MPa; eminRepresents the minimum value of the Young's modulus of the core sample, and has a unit of MPa.
Preferably, in the step 4, the core sample is normalizedResidual intensity influencing factor KREThe calculation formula is as follows:
in the formula, σREThe residual strength of the core sample is expressed in MPa; sigmaBCThe fracture pressure of the core sample is expressed in MPa.
Preferably, in the step 5, the normalized fracture toughness influence factor K of the core sampleICThe calculation formula is as follows:
in the formula (I), the compound is shown in the specification,the measured fracture toughness of the core sample is expressed in kg/cm2;KIC_maxThe maximum value of the fracture toughness of the core sample is expressed in kg/cm2;KIC_minRepresents the minimum value of the fracture toughness of the core sample and has the unit of kg/cm2;
Normalized tensile strength impact factor K of core sampleTThe calculation formula is as follows:
in the formula (I), the compound is shown in the specification,the tensile strength of the measured core sample is expressed in MPa; sigmaT_maxThe maximum value of the tensile strength of the core sample is expressed in MPa; sigmaT_minRepresents the minimum value of tensile strength of the core sample in MPa.
Preferably, in step 6, the weighting coefficient calculation formula is:
wherein i represents a parameter number, wherein i ═ 1 represents the horizontal ground stress difference of the core sample, i ═ 2 represents the brittleness of the core sample, i ═ 3 represents the fracture toughness of the core sample, i ═ 4 represents the tensile strength of the core sample, and i ═ 5 represents the residual strength of the core sample; ziA weighting coefficient representing a parameter; ri 2The degree of fit of the parameter is indicated.
Preferably, in the step 7, the crack correction factor KCOThe calculation formula is as follows:
KCO=C1×C2 (9)
wherein the porosity C is normalized1The calculation formula is as follows:
where por represents the porosity of the target well in%; pormaxRepresents the maximum value of porosity in the logging data, in%; porminRepresents the minimum value of porosity in the logging data in%;
normalized porosity spectrum mean square error C2The calculation formula is as follows:
wherein porV represents the mean square error of the porosity spectrum of the target well; porVmaxRepresenting the maximum value of the mean square error of the porosity spectrum of the logging data; porVminThe minimum value of the mean square error of the porosity spectrum of the logging data is represented.
Preferably, in step 8, the target reservoir fracability evaluation model is:
FI=Z1Kh+Z2BI+Z3KIC+Z4σT+Z5KRE+KCO (12)
in the formula, FI represents a fracability evaluation parameter.
The invention has the following beneficial effects:
the method is based on continuous elastomechanics and linear elasticity theory, comprehensively considers the complex reservoir fracture characteristics, introduces four evaluation parameters of fracture toughness, horizontal ground stress difference, residual strength and tensile strength and fracture correction factors, solves the problem of poor evaluation accuracy of the deep fracturing property of the reservoir when a complex reservoir is evaluated by using a single brittle parameter, and makes up the defects of the existing single brittle parameter evaluation;
the method disclosed by the invention is more in line with the characteristics of complex reservoir pore cavern development, strong anisotropy and large influence of primary fractures on the fracturing property, and utilizes results of a ground stress test experiment, a rock triaxial mechanics experiment, a rock destruction experiment and a Brazilian fracture experiment to carry out linear regression analysis to calibrate the weighting coefficient of each parameter, so that the fracturing property of the complex reservoir is quantized, and the fracturing property evaluation of the complex reservoir is more accurate;
the method provided by the invention constructs the fracturing property evaluation model based on the characteristics of the target reservoir, realizes accurate evaluation of the fracturing property of the reservoir, has strong applicability, and is beneficial to efficient development and safe exploitation of oil and gas resources of the complex reservoir.
Drawings
Fig. 1 is a flow chart of a complex reservoir fracturing evaluation method.
FIG. 2 is a graph of regression analysis between the differential stress at rock sample level and fracture pressure.
FIG. 3 is a graph of regression analysis between rock sample brittleness and fracture pressure.
FIG. 4 is a graph of regression analysis between fracture toughness and fracture pressure for rock samples.
FIG. 5 is a graph of regression analysis between tensile strength and burst pressure for rock samples.
FIG. 6 is a graph of regression analysis between residual strength and fracture pressure of rock samples.
Detailed Description
The following will further explain the embodiments of the present invention by taking the example of a well in a carbonate block and the accompanying drawings as an example:
taking a certain example of a well in the ancient city area of Tandong as an example, the buried depth of a carbonate reservoir is 6088-6780 m, and the fracturing evaluation method for the complex reservoir provided by the invention is used for carrying out fracturing evaluation on the reservoir, and as shown in figure 1, the fracturing evaluation method specifically comprises the following steps:
step 1, selecting 26 standard-size rock core samples with the block size of 250mm from rock core samples collected at 9 depth points of a sample well of a carbonic acid rock block to perform rock physics experiments, wherein 18 rock core samples are used for sequentially performing a ground stress test experiment, a rock triaxial mechanics experiment and a rock destruction experiment, and 8 rock core samples are used for performing a Brazilian splitting experiment.
Step 2, performing a ground stress test experiment on the 18 rock core samples selected in the step 1, measuring the maximum ground stress and the minimum ground stress of the rock core samples, calculating the horizontal ground stress difference and the horizontal ground stress difference influence factor of the rock core samples, and calculating the horizontal ground stress difference influence factor KhThe calculation formula is as follows:
in the formula, σHRepresents the maximum ground stress of the core sample, in MPa; sigmahThe minimum stress of the core sample is expressed in MPa.
Step 3, performing rock triaxial mechanical experiment on the core sample by using a triaxial mechanical test system, measuring the Young modulus and the Poisson ratio of 18 core samples, and calculating the brittleness of the core sample according to the Young modulus and the Poisson ratio of the core sample, wherein the formula (2) is as follows:
BI=(μ+E)/2 (2)
where μ represents the normalized poisson's ratio and E represents the normalized young's modulus, the calculation formula is:
in the formula, mu*Representing the measured poisson ratio of the core sample; mu.smaxRepresenting the maximum value of the Poisson's ratio of the core sample; mu.sminRepresenting the minimum value of the Poisson's ratio of the core sample; e*The measured Young modulus of the core sample is expressed in MPa; emaxThe maximum value of the Young modulus of the core sample is expressed in MPa; eminRepresents the minimum value of the Young's modulus of the core sample, and has a unit of MPa.
Step 4, performing a rock destruction experiment on 18 rock core samples by using a triaxial mechanical testing system, measuring the residual strength and the fracture pressure of the rock core samples, and calculating normalized residual strength influence factors of the rock core samples, wherein the normalized residual strength influence factors are shown in formula (5):
normalized residual intensity influencing factor KREThe calculation formula is as follows:
in the formula, σREThe residual strength of the core sample is expressed in MPa; sigmaBCThe fracture pressure of the core sample is expressed in MPa.
Step 5, carrying out Brazilian splitting experiments on the 8 rock core samples selected in the step 1, measuring the fracture toughness and the tensile strength of the rock core samples, and calculating fracture toughness influence factors and tensile strength influence factors of the rock core samples as follows:
normalized fracture toughness influencing factor K of core sampleICThe calculation formula is as follows:
in the formula (I), the compound is shown in the specification,the measured fracture toughness of the core sample is expressed in kg/cm2;KIC_maxThe maximum value of the fracture toughness of the core sample is expressed in kg/cm2;KIC_minRepresents the minimum value of the fracture toughness of the core sample and has the unit of kg/cm2;
Normalized tensile strength impact factor K of core sampleTThe calculation formula is as follows:
in the formula (I), the compound is shown in the specification,the tensile strength of the measured core sample is expressed in MPa; sigmaT_maxThe maximum value of the tensile strength of the core sample is expressed in MPa; sigmaT_minRepresents the minimum value of tensile strength of the core sample in MPa.
Step 6, performing single correlation analysis on horizontal ground stress difference, brittleness, fracture toughness, tensile strength and residual strength of the rock core sample and fracture pressure of the rock sample respectively, wherein regression analysis results are shown in fig. 2 to 6;
determining the weighting coefficient of each parameter according to the fitting degree, wherein the calculation formula of the weighting coefficient is as follows:
wherein i represents a parameter number, wherein i ═ 1 represents the horizontal ground stress difference of the core sample, i ═ 2 represents the brittleness of the core sample, i ═ 3 represents the fracture toughness of the core sample, i ═ 4 represents the tensile strength of the core sample, and i ═ 5 represents the residual strength of the core sample; ziA weighting coefficient representing a parameter; ri 2Representing the degree of fit of the parameter;
in the embodiment, the fitting degree R of the horizontal ground stress difference and the fracture pressure of the rock sample1 20.3975 rockDegree of matching R of brittleness and fracture pressure of stone sample2 20.0416, the degree of match R of fracture toughness to fracture pressure of rock sample3 20.6713, rock sample tensile Strength to fracture pressure Fit R4 2The rock sample residual strength to fracture pressure fit is R0.72655 20.8323; and (3) substituting the weight coefficient into the formula (8) for calculation to obtain that the weight coefficient corresponding to the horizontal ground stress difference of the rock sample is 0.14, the weight coefficient corresponding to the brittleness of the rock sample is 0.09, the weight coefficient corresponding to the fracture toughness of the rock sample is 0.23, the weight coefficient corresponding to the tensile strength of the rock sample is 0.25, and the weight coefficient corresponding to the residual strength of the rock sample is 0.29.
And 7, acquiring the normalized porosity and the normalized porosity spectrum mean square error of the target well according to the logging information in the research block, and calculating a fracture correction factor by using the normalized porosity and the normalized porosity spectrum mean square error, wherein the formula is shown as (9):
crack correction factor KCOThe calculation formula is as follows:
KCO=C1×C2 (9)
wherein the porosity C is normalized1The calculation formula is as follows:
where por represents the porosity of the target well in%; pormaxRepresents the maximum value of porosity in the logging data, in%; porminRepresents the minimum value of porosity in the logging data in%;
normalized porosity spectrum mean square error C2The calculation formula is as follows:
wherein porV represents the mean square error of the porosity spectrum of the target well; porVmaxRepresentation loggingMaximum value of mean square error of porosity spectrum of the data; porVminThe minimum value of the mean square error of the porosity spectrum of the logging data is represented.
Step 8, establishing a fracturing property evaluation model of the target reservoir according to the weighting coefficients and the fracture correction factors of the parameters of the core sample, as shown in the formula (12):
FI=0.09Kh+0.14BI+0.23KIC+0.25σT+0.29KRE+KCO (12)
and calculating a fracturing evaluation parameter FI by using the fracturing evaluation model of the target reservoir, and evaluating the fracturing of the target reservoir.
It is to be understood that the above description is not intended to limit the present invention, and the present invention is not limited to the above examples, and those skilled in the art may make modifications, alterations, additions or substitutions within the spirit and scope of the present invention.
Claims (9)
1. The method for evaluating the fracturing property of the complex reservoir is characterized by comprising the following steps:
step 1, selecting a plurality of core samples in a complex reservoir based on the complex reservoir in a research block;
step 2, performing an earth stress test experiment on the core sample, measuring the maximum earth stress and the minimum earth stress of the core sample, and calculating the horizontal earth stress difference and the horizontal earth stress difference influence factor of the core sample;
step 3, performing a rock triaxial mechanical experiment on the core sample by using a triaxial mechanical test system, measuring the Young modulus and Poisson ratio of the core sample, and calculating the brittleness of the core sample;
step 4, performing a rock destruction experiment on the core sample by using a triaxial mechanical testing system, measuring the residual strength and the fracture pressure of the core sample, and calculating a normalized residual strength influence factor of the core sample;
step 5, performing Brazilian splitting experiments on the core sample, measuring the fracture toughness and the tensile strength of the core sample, and calculating fracture toughness influence factors and tensile strength influence factors of the core sample;
step 6, performing single correlation analysis on the horizontal ground stress difference, brittleness, fracture toughness, tensile strength and residual strength of the rock core sample and the fracture pressure of the rock sample respectively, and determining the weighting coefficient of each parameter according to the fitting degree;
step 7, acquiring the normalized porosity and the normalized porosity spectrum mean square error of the target well according to the logging information in the research block, and calculating a fracture correction factor;
and 8, establishing a target reservoir fracturing evaluation model according to the weighting coefficients and the fracture correction factors of the parameters of the core sample, calculating the target reservoir fracturing evaluation parameters, and determining the fracturing of the target reservoir.
2. The method for evaluating the fracturability of a complex reservoir according to claim 1, wherein in the step 1, a part of core samples are used for performing an earth stress test experiment, a rock triaxial mechanical experiment and a rock failure experiment in sequence, and another part of core samples are used for performing a brazilian fracture experiment.
3. The method for evaluating the fracability of a complex reservoir according to claim 1, wherein in the step 2, the horizontal ground stress difference of the core sample is the difference between the maximum ground stress and the minimum ground stress of the rock sample;
horizontal ground stress difference influence factor K of core samplehThe calculation formula is as follows:
in the formula, σHRepresents the maximum ground stress of the core sample, in MPa; sigmahThe minimum stress of the core sample is expressed in MPa.
4. The method for evaluating the fracability of a complex reservoir as claimed in claim 1, wherein in the step 3, the brittle BI of the core sample is calculated by the following formula:
BI=(μ+E)/2 (2)
where μ represents the normalized poisson's ratio and E represents the normalized young's modulus, the calculation formula is:
in the formula, mu*Representing the measured poisson ratio of the core sample; mu.smaxRepresenting the maximum value of the Poisson's ratio of the core sample; mu.sminRepresenting the minimum value of the Poisson's ratio of the core sample; e*The measured Young modulus of the core sample is expressed in MPa; emaxThe maximum value of the Young modulus of the core sample is expressed in MPa; eminRepresents the minimum value of the Young's modulus of the core sample, and has a unit of MPa.
5. The method for evaluating the fracability of a complex reservoir as claimed in claim 1, wherein in the step 4, the normalized residual strength influence factor K of the core sampleREThe calculation formula is as follows:
in the formula, σREThe residual strength of the core sample is expressed in MPa; sigmaBCThe fracture pressure of the core sample is expressed in MPa.
6. The method for evaluating the fracability of a complex reservoir as claimed in claim 1, wherein in the step 5, the normalized fracture toughness influence factor K of the core sampleICThe calculation formula is as follows:
in the formula (I), the compound is shown in the specification,the measured fracture toughness of the core sample is expressed in kg/cm2;KIC_maxThe maximum value of the fracture toughness of the core sample is expressed in kg/cm2;KIC_minRepresents the minimum value of the fracture toughness of the core sample and has the unit of kg/cm2;
Normalized tensile strength impact factor K of core sampleTThe calculation formula is as follows:
in the formula (I), the compound is shown in the specification,the tensile strength of the measured core sample is expressed in MPa; sigmaT_maxThe maximum value of the tensile strength of the core sample is expressed in MPa; sigmaT_minRepresents the minimum value of tensile strength of the core sample in MPa.
7. The method for evaluating the fracability of a complex reservoir according to claim 1, wherein in the step 6, the weighting coefficient calculation formula is as follows:
wherein i represents a parameter number, wherein i ═ 1 represents the horizontal ground stress difference of the core sample, i ═ 2 represents the brittleness of the core sample, i ═ 3 represents the fracture toughness of the core sample, i ═ 4 represents the tensile strength of the core sample, and i ═ 5 represents the residual strength of the core sample; ziRepresenting parametersThe weighting coefficient of (2); ri 2The degree of fit of the parameter is indicated.
8. The method for evaluating the fracability of a complex reservoir as claimed in claim 1, wherein in step 7, the fracture correction factor KCOThe calculation formula is as follows:
KCO=C1×C2 (9)
wherein the porosity C is normalized1The calculation formula is as follows:
where por represents the porosity of the target well in%; pormaxRepresents the maximum value of porosity in the logging data, in%; porminRepresents the minimum value of porosity in the logging data in%;
normalized porosity spectrum mean square error C2The calculation formula is as follows:
wherein porV represents the mean square error of the porosity spectrum of the target well; porVmaxRepresenting the maximum value of the mean square error of the porosity spectrum of the logging data; porVminThe minimum value of the mean square error of the porosity spectrum of the logging data is represented.
9. The method for evaluating the fracability of the complex reservoir according to claim 1, wherein in the step 8, the target reservoir fracability evaluation model is as follows:
FI=Z1Kh+Z2BI+Z3KIC+Z4σT+Z5KRE+KCO (12)
in the formula, FI represents a fracability evaluation parameter.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011087617.3A CN112307601B (en) | 2020-10-13 | 2020-10-13 | Complex reservoir fracturing property evaluation method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011087617.3A CN112307601B (en) | 2020-10-13 | 2020-10-13 | Complex reservoir fracturing property evaluation method |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112307601A true CN112307601A (en) | 2021-02-02 |
CN112307601B CN112307601B (en) | 2022-06-24 |
Family
ID=74489811
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011087617.3A Active CN112307601B (en) | 2020-10-13 | 2020-10-13 | Complex reservoir fracturing property evaluation method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112307601B (en) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112784437A (en) * | 2021-02-04 | 2021-05-11 | 中国人民解放军63863部队 | System for evaluating damage efficiency of air-defense missile to helicopter target |
CN113236242A (en) * | 2021-06-22 | 2021-08-10 | 中国石油天然气股份有限公司 | Method for determining maximum operating pressure of expansion gas storage |
CN113283108A (en) * | 2021-06-10 | 2021-08-20 | 西安石油大学 | Method and system for quantitatively evaluating fracturing property of shale oil reservoir |
CN113776963A (en) * | 2021-08-04 | 2021-12-10 | 中国石油大学(华东) | Method for calculating fracture toughness by using spherical indentation method |
CN114544367A (en) * | 2022-02-21 | 2022-05-27 | 西北大学 | Reservoir fracturing evaluation and fracturing scheme design method based on core experiment |
Citations (2)
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 |
CN111353218A (en) * | 2020-02-20 | 2020-06-30 | 西安石油大学 | Logging quantitative evaluation method for coal bed gas-dense gas reservoir compaction property |
-
2020
- 2020-10-13 CN CN202011087617.3A patent/CN112307601B/en active Active
Patent Citations (2)
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 |
CN111353218A (en) * | 2020-02-20 | 2020-06-30 | 西安石油大学 | Logging quantitative evaluation method for coal bed gas-dense gas reservoir compaction property |
Non-Patent Citations (2)
Title |
---|
张永浩: "页岩气储层可压裂性评价方法实验研究", 《中国优秀硕士论文库》 * |
郭冀宁: "L地区致密砂岩可压裂性评价", 《中国优秀硕士论文库》 * |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112784437A (en) * | 2021-02-04 | 2021-05-11 | 中国人民解放军63863部队 | System for evaluating damage efficiency of air-defense missile to helicopter target |
CN112784437B (en) * | 2021-02-04 | 2022-11-29 | 中国人民解放军63863部队 | System for evaluating damage efficiency of air-defense missile to helicopter target |
CN113283108A (en) * | 2021-06-10 | 2021-08-20 | 西安石油大学 | Method and system for quantitatively evaluating fracturing property of shale oil reservoir |
CN113283108B (en) * | 2021-06-10 | 2023-09-22 | 西安石油大学 | Method and system for quantitatively evaluating fracturing property of shale oil reservoir |
CN113236242A (en) * | 2021-06-22 | 2021-08-10 | 中国石油天然气股份有限公司 | Method for determining maximum operating pressure of expansion gas storage |
CN113236242B (en) * | 2021-06-22 | 2023-08-08 | 中国石油天然气股份有限公司 | Method for determining maximum operating pressure of capacity-expansion gas storage |
CN113776963A (en) * | 2021-08-04 | 2021-12-10 | 中国石油大学(华东) | Method for calculating fracture toughness by using spherical indentation method |
CN113776963B (en) * | 2021-08-04 | 2023-07-28 | 中国石油大学(华东) | Method for calculating fracture toughness by spherical indentation method |
CN114544367A (en) * | 2022-02-21 | 2022-05-27 | 西北大学 | Reservoir fracturing evaluation and fracturing scheme design method based on core experiment |
CN114544367B (en) * | 2022-02-21 | 2024-02-09 | 西北大学 | Reservoir fracturing property evaluation and fracturing scheme design method based on core experiment |
Also Published As
Publication number | Publication date |
---|---|
CN112307601B (en) | 2022-06-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112307601B (en) | Complex reservoir fracturing property evaluation method | |
CN102175832B (en) | Method for determining optimal saturation computing model for typical reservoir | |
CN108009705A (en) | A kind of shale reservoir compressibility evaluation method based on support vector machines technology | |
CN107038313B (en) | Layered crustal stress fine description method based on numerical value core | |
CN113820750B (en) | Method for quantitatively predicting mudstone structural cracks based on elastoplastic mechanics | |
CN105651966A (en) | Shale oil and gas high-quality reservoir stratum evaluation method and parameter determination method | |
US20070112518A1 (en) | Method and Apparatus for Measuring the Wettability of Geological Formations | |
CN105134189A (en) | Logging GeoMechanics Identify Reservoir (LogGMIR) method | |
CN106569288B (en) | Fractured reservoir quality evaluation method based on reservoir fracture effectiveness cluster analysis | |
Boutaleb et al. | Sedimentological, petrophysical, and geochemical controls on deep marine unconventional tight limestone and dolostone reservoir: Insights from the Cenomanian/Turonian oceanic anoxic event 2 organic-rich sediments, Southeast Constantine Basin, Algeria | |
CN113419284B (en) | Method for identifying physical facies double desserts of well logging rock based on cluster analysis | |
CN104749622A (en) | Method for quantitatively characterizing shale compressibility based on rock physics | |
AU2007275040B2 (en) | Method and system for interpreting borehole geological data | |
CN103867198A (en) | Method for distinguishing formation density of carbonatite natural gas reservoir | |
CN105114067A (en) | Lithology electrofacies method | |
CN108374657A (en) | Well breakpoint automatic identifying method | |
Wimmers et al. | Integration of sedimentology, petrophysics and rock typing as key to understanding a tight gas reservoir | |
CN115030707A (en) | Rapid evaluation method of oil shale dessert | |
CN109211745A (en) | A kind of restoration methods of the evolutionary process rich in organic matter mud shale porosity | |
CN110322363A (en) | Shale gas reservoir reconstruction volume calculation method and system | |
CN112832738A (en) | Clastic rock cumulative hydrocarbon generation strength determination method and dessert layer identification and evaluation method | |
CN109738955B (en) | Metamorphic rock lithology comprehensive judgment method based on component-structure classification | |
Estrada et al. | Pilot phase of the Aguada Pichana Este block, gas window | |
CN115099014A (en) | Natural gas well geological exploration reserves calculation method based on logging while drilling | |
CN110454155B (en) | Method for determining age of quaternary stratum by using susceptibility logging method |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |