CN109917102A - A kind of brittleness evaluation method considering rock rupture complexity - Google Patents
A kind of brittleness evaluation method considering rock rupture complexity Download PDFInfo
- Publication number
- CN109917102A CN109917102A CN201910081147.0A CN201910081147A CN109917102A CN 109917102 A CN109917102 A CN 109917102A CN 201910081147 A CN201910081147 A CN 201910081147A CN 109917102 A CN109917102 A CN 109917102A
- Authority
- CN
- China
- Prior art keywords
- young
- modulus
- rock
- brittleness
- normalize
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Abstract
The present invention relates to a kind of brittleness evaluation methods for considering rock rupture complexity, belong to the brittleness evaluation field of rock.In order to solve the problems in the existing technology, the present invention provides a kind of brittleness evaluation methods for considering rock rupture complexity, comprising the following steps: step A, acquires rock rupture sample, and collects the basic parameter of rock；Step B, dilative angle, Young's modulus and peak strain are normalized by following formula respectively: step C, the brittleness of rock core is evaluated according to data obtained above and following formula.The present invention can be than more fully characterizing the Brittleness of rock.
Description
Technical field
The present invention relates to a kind of brittleness evaluation methods for considering rock rupture complexity, belong to the brittleness evaluation neck of rock
Domain.
Background technique
The exploitation of shale (oil) gas reservoir, key technology first is that using hydraulic fracturing technology seam, improve producing well
Involve range and reservoir its permeability.Mud shale Brittleness is related to its compressibility, and therefore, rock brittleness becomes its feature
Evaluate importance.
The quantitative evaluation method of existing rock brittleness mainly uses the evaluation method based on brittle mineral --- Jarvie etc.
Deng there are defects below for this method: (1) cannot be used for the rupture process of evaluation deep layer shale；(2) it is non-to ignore intrinsic fracture etc.
The influence of brittleness factors；(3) mineral constituent is not the leading factor of shale rupture process；(4) rock rupture Morphological Characterization is excessively
It is subjective.
Summary of the invention
In order to solve the problems in the existing technology, the crisp of rock rupture complexity is considered the present invention provides a kind of
Property evaluation method.
The present invention solves technical solution provided by above-mentioned technical problem: a kind of to consider the crisp of rock rupture complexity
Property evaluation method, comprising the following steps:
Step A, rock rupture sample is acquired, and collects the basic parameter of rock；
Step B, dilative angle, Young's modulus and peak strain are normalized by following formula respectively:
In formula: E is Young's modulus；EminFor minimum Young's modulus；EmaxFor maximum Young's modulus；EnTo normalize Young mould
Amount；ψ is dilative angle；ψmaxFor maximum dilative angle；ψminFor minimum dilative angle；ψnTo normalize dilative angle；εpFor peak strain；
εpmaxFor peak-peak strain；εpminFor minimum peak strain；εpnTo normalize overstrain；
Step C, the brittleness of rock core is evaluated according to data obtained above and following formula, judgement schematics are as follows:
BI=W1En+W2ψn+W3εpn
W1+W2+W3=1
In formula: EnTo normalize Young's modulus；ψnTo normalize dilative angle；εpnTo normalize overstrain；BIRefer to for brittleness
Number；W1For the weight coefficient for normalizing Young's modulus；W2For the weight coefficient for normalizing dilative angle；W3To normalize overstrain
Weight coefficient.
Further technical solution is that basic parameter includes dilative angle, Young's modulus and the peak value of rock in the step A
Strain.
Further technical solution is that the weight coefficient of normalization Young's modulus is 0.262 in the step C；Normalization
The weight coefficient of dilative angle is 0.353；The weight coefficient for normalizing overstrain is 0.385.
Beneficial effects of the present invention: the present invention can be than more fully characterizing the Brittleness of rock.
Detailed description of the invention
Fig. 1 is different phase schematic diagram of each parameter on stress-strain diagram in the present invention.
Specific embodiment
Invention is further described in detail with reference to the accompanying drawings and examples.
The brittleness evaluation method of a kind of consideration rock rupture complexity of the invention as shown in Figure 1:, including following step
It is rapid:
Step A, rock rupture sample is acquired, and collects the basic parameter of rock, basic parameter includes that cutting for rock is swollen
Angle, Young's modulus and peak strain；
Step B, dilative angle, Young's modulus and peak strain are normalized by following formula respectively:
In formula: E is Young's modulus；EminFor minimum Young's modulus；EmaxFor maximum Young's modulus；EnTo normalize Young mould
Amount；ψ is dilative angle；ψmaxFor maximum dilative angle；ψminFor minimum dilative angle；ψnTo normalize dilative angle；εpFor peak strain；
εpmaxFor peak-peak strain；εpminFor minimum peak strain；εpnTo normalize overstrain；
Step C, the brittleness of rock core is evaluated according to data obtained above and following formula, judgement schematics are as follows:
BI=W1En+W2ψn+W3εpn
W1+W2+W3=1
In formula: EnTo normalize Young's modulus；ψnTo normalize dilative angle；εpnTo normalize overstrain；BIRefer to for brittleness
Number；W1For the weight coefficient for normalizing Young's modulus；W2For the weight coefficient for normalizing dilative angle；W3To normalize overstrain
Weight coefficient.
Weight coefficient is calculated using grey correlation theory, gray system theory, which is proposed, carries out grey correlation to each subsystem
The concept for spending analysis, seeks the numerical relation in system between each subsystem by certain method.For between two systems
Factor, the degree of association is defined as at any time or different object and the measurement of relevance size that changes.What if two factors changed
Trend is with uniformity, that is, is that the two degree of association is higher；Conversely, then lower.Currently, correlation analysis has been applied to science
The every field of research, application are very extensive.
It mainly comprises the steps that
1. determining subnumber column and reference sequence
It will affect brittle each factor to arrange as subnumber, subnumber column expression formula are as follows:
ri=(ri1,ri2,......,rim)
r0=(r1,r2,......,rm)
2. calculating utility function
For the interference for eliminating physical quantity unit, data are handled using normalized method.According to preliminary analysis knot
Fruit is energy to failure, dilative angle and overstrain respectively for three parameters relevant to brittleness.Not according to statistical result early period
Together, it is handled respectively according to following two method:
The more bigger more excellent type index of a, the calculating of utility function are as follows:
The smaller more excellent type index of b, the calculating of utility function are as follows:
Wherein, (rij)min(rij)max--- the minimum value and maximum value of sample.
It can thus be concluded that Jacobian matrix:
R=[bij]m×n
3. seeking incidence coefficient
Wherein:
Δi=| x0(j)-xi(j)|
ρ ∈ (0 ,+∞) is resolution ratio, and ρ is smaller, and resolving power is bigger, and the value interval of ρ is [0,1].
4. solving the degree of association
The incidence coefficient and the degree of association of subsequence are all larger than 0, and the degree of association of subsequence should be the flat of each incidence coefficient
Mean value, therefore available:
5. solving weight coefficient
All parameters have opposite impacts on degree to brittleness, for the influence size for comparing parameters, using it is each because
The specific gravity that the plain degree of association accounts for total correlation degree calculates its weight coefficient:
For energy to failure and overstrain, using smaller more excellent type utility function, for rupturing complexity coefficient and dilative angle,
Using more bigger more excellent type utility function, processing result is that the weight coefficient of normalization Young's modulus is 0.262；Normalize dilative angle
Weight coefficient be 0.353；The weight coefficient for normalizing overstrain is 0.385.
Therefore, establishing is suitble to the brittleness evaluation method of shale as follows:
BI=0.262En+0.353ψn+0.385εpn
The above is not intended to limit the present invention in any form, although the present invention takes off through the foregoing embodiment
Show, however, it is not intended to limit the invention, any person skilled in the art, is not departing from technical solution of the present invention range
It is interior, made when the technology contents using the disclosure above and change or be modified to the equivalent embodiments of equivalent variations a bit, but it is all not
Be detached from technical solution of the present invention content, according to the technical essence of the invention it is to the above embodiments it is any it is simple modification,
Equivalent variations and modification, all of which are still within the scope of the technical scheme of the invention.
Claims (3)
1. a kind of brittleness evaluation method for considering rock rupture complexity, which comprises the following steps:
Step A, rock rupture sample is acquired, and collects the basic parameter of rock；
Step B, dilative angle, Young's modulus and peak strain are normalized by following formula respectively:
In formula: E is Young's modulus；EminFor minimum Young's modulus；EmaxFor maximum Young's modulus；EnTo normalize Young's modulus；ψ
For dilative angle；ψmaxFor maximum dilative angle；ψminFor minimum dilative angle；ψnTo normalize dilative angle；εpFor peak strain；εpmaxFor
Peak-peak strain；εpminFor minimum peak strain；εpnTo normalize overstrain；
Step C, the brittleness of rock core is evaluated according to data obtained above and following formula, judgement schematics are as follows:
BI=W1En+W2ψn+W3εpn
W1+W2+W3=1
In formula: EnTo normalize Young's modulus；ψnTo normalize dilative angle；εpnTo normalize overstrain；BIFor brittleness index；
W1For the weight coefficient for normalizing Young's modulus；W2For the weight coefficient for normalizing dilative angle；W3For the power for normalizing overstrain
Weight coefficient.
2. a kind of brittleness evaluation method for considering rock rupture complexity according to claim 1, which is characterized in that institute
State dilative angle, Young's modulus and the peak strain that basic parameter in step A includes rock.
3. a kind of brittleness evaluation method for considering rock rupture complexity according to claim 1, which is characterized in that institute
Stating and normalizing the weight coefficient of Young's modulus in step C is 0.262；The weight coefficient for normalizing dilative angle is 0.353；Normalization
The weight coefficient of overstrain is 0.385.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910081147.0A CN109917102A (en) | 2019-01-28 | 2019-01-28 | A kind of brittleness evaluation method considering rock rupture complexity |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910081147.0A CN109917102A (en) | 2019-01-28 | 2019-01-28 | A kind of brittleness evaluation method considering rock rupture complexity |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN109917102A true CN109917102A (en) | 2019-06-21 |
Family
ID=66960998
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910081147.0A Pending CN109917102A (en) | 2019-01-28 | 2019-01-28 | A kind of brittleness evaluation method considering rock rupture complexity |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109917102A (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110715859A (en) * | 2019-10-23 | 2020-01-21 | 成都理工大学 | Brittleness index evaluation method based on elastic-plastic deformation and fracture strength |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105156103A (en) * | 2015-09-29 | 2015-12-16 | 西南石油大学 | Debris-core-borehole-reservoir multiscale shale reservoir three-dimensional fracturing evaluation method |
| CN106547034A (en) * | 2016-11-09 | 2017-03-29 | 西南石油大学 | A kind of method for calculating compact reservoir rock brittleness index |
| CN108615102A (en) * | 2016-12-12 | 2018-10-02 | 中国石油天然气股份有限公司 | A method of the fine and close oil gas pressure break of evaluation forms network fracture ability |
-
2019
- 2019-01-28 CN CN201910081147.0A patent/CN109917102A/en active Pending
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105156103A (en) * | 2015-09-29 | 2015-12-16 | 西南石油大学 | Debris-core-borehole-reservoir multiscale shale reservoir three-dimensional fracturing evaluation method |
| CN106547034A (en) * | 2016-11-09 | 2017-03-29 | 西南石油大学 | A kind of method for calculating compact reservoir rock brittleness index |
| CN108615102A (en) * | 2016-12-12 | 2018-10-02 | 中国石油天然气股份有限公司 | A method of the fine and close oil gas pressure break of evaluation forms network fracture ability |
Non-Patent Citations (1)
| Title |
|---|
| 楼裕胜著: "《商业银行竞争力评价方法研究》", 31 October 2007, 浙江大学出版社 * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110715859A (en) * | 2019-10-23 | 2020-01-21 | 成都理工大学 | Brittleness index evaluation method based on elastic-plastic deformation and fracture strength |
| CN110715859B (en) * | 2019-10-23 | 2020-09-08 | 成都理工大学 | Brittleness index evaluation method based on elastic-plastic deformation and fracture strength |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108009716B (en) | Multi-level evaluation method for influence factors of horizontal well volume fracturing effect | |
| WO2020220745A1 (en) | Quantitative grading method for lost-circulation material selection optimization in well completion of well drilling | |
| CN105626025B (en) | A kind of compressibility evaluation method of shale reservoir pressure break | |
| CN108229811B (en) | Method for evaluating water injection effect of fractured-vuggy carbonate reservoir | |
| CN111738620B (en) | Well leakage risk prediction and leakage stoppage decision system and method based on association rules | |
| CN110714753A (en) | Shale reservoir compressibility evaluation method | |
| CN104832166A (en) | Initial productivity prediction method of shale gas horizontal well | |
| CN108960651A (en) | A kind of integrated evaluating method of densification oil-gas reservoir multistage fracturing horizontal well completion efficiency | |
| CN104747185A (en) | Heterogeneous reservoir stratum synthetic classifying evaluation method | |
| WO2014166244A1 (en) | Method for evaluating water inrush weakness of coal seam floor | |
| CN109184660A (en) | A kind of shale gas evaluating production capacity method based on Production log data | |
| Zifei et al. | A study on remaining oil distribution in a carbonate oil reservoir based on reservoir flow units | |
| CN106250984A (en) | The determination methods of the oil water relation pattern of oil well and device | |
| CN109917102A (en) | A kind of brittleness evaluation method considering rock rupture complexity | |
| He et al. | Quantitative evaluation and influencing factors analysis of the brittleness of deep shale reservoir based on multiply rock mechanics experiments | |
| CN108661628A (en) | One kind being based on the preferred engineering dessert quantitative evaluation method of parameter | |
| CN111242490A (en) | Low-permeability reservoir characteristic master control factor extraction method and device and readable storage medium | |
| CN109670729A (en) | A kind of top plate aquifer water well evaluation method | |
| CN105719065A (en) | Comprehensive evaluation method of complex oil reservoir reserves quality classification | |
| Chi et al. | Prediction of top-coal caving and drawing characteristics by the analytic hierarchy process-fuzzy discrimination method in extra-thick coal seams | |
| CN111611704B (en) | Method and device for establishing water-drive series curve plate and readable storage medium | |
| Deyong et al. | Coal and gas outburst early-warning technology and application based on AHP and extension theory | |
| Liu et al. | Research and application of sidewall stability predic-tion method based on analytic hierarchy process and fuzzy integrative evaluation method | |
| Gong et al. | Productivity Prediction Methods of Multi-Fractured Horizontal Wells in Low-Permeability Reservoirs: Progress and Challenges | |
| CN105808936A (en) | Method for quickly judging and identifying danger degree of landslide |
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 | ||
| RJ01 | Rejection of invention patent application after publication | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20190621 |