CN112304770A - Method and system for quantitatively analyzing fracture complexity after fracturing - Google Patents
Method and system for quantitatively analyzing fracture complexity after fracturing Download PDFInfo
- Publication number
- CN112304770A CN112304770A CN201910634085.1A CN201910634085A CN112304770A CN 112304770 A CN112304770 A CN 112304770A CN 201910634085 A CN201910634085 A CN 201910634085A CN 112304770 A CN112304770 A CN 112304770A
- Authority
- CN
- China
- Prior art keywords
- complexity
- judgment conditions
- equal
- class
- judging
- 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
- 238000000034 method Methods 0.000 title claims abstract description 34
- 238000010276 construction Methods 0.000 claims abstract description 12
- 238000004364 calculation method Methods 0.000 claims description 7
- 238000003860 storage Methods 0.000 claims description 4
- 238000012986 modification Methods 0.000 abstract description 8
- 230000004048 modification Effects 0.000 abstract description 8
- 230000009286 beneficial effect Effects 0.000 abstract description 5
- 230000000694 effects Effects 0.000 abstract description 5
- 238000013461 design Methods 0.000 abstract description 4
- 238000005457 optimization Methods 0.000 abstract description 4
- 206010017076 Fracture Diseases 0.000 description 47
- 208000010392 Bone Fractures Diseases 0.000 description 41
- 238000011156 evaluation Methods 0.000 description 8
- 239000007789 gas Substances 0.000 description 7
- 238000011158 quantitative evaluation Methods 0.000 description 5
- 238000004458 analytical method Methods 0.000 description 3
- 238000012512 characterization method Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
- 238000001914 filtration Methods 0.000 description 2
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 2
- 230000010355 oscillation Effects 0.000 description 2
- 238000005086 pumping Methods 0.000 description 2
- 238000004445 quantitative analysis Methods 0.000 description 2
- BVKZGUZCCUSVTD-UHFFFAOYSA-L Carbonate Chemical compound [O-]C([O-])=O BVKZGUZCCUSVTD-UHFFFAOYSA-L 0.000 description 1
- 206010010214 Compression fracture Diseases 0.000 description 1
- 238000013211 curve analysis Methods 0.000 description 1
- 238000007405 data analysis Methods 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 239000003345 natural gas Substances 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 239000003209 petroleum derivative Substances 0.000 description 1
- 238000003825 pressing Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 239000011435 rock Substances 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N3/00—Investigating strength properties of solid materials by application of mechanical stress
- G01N3/08—Investigating strength properties of solid materials by application of mechanical stress by applying steady tensile or compressive forces
- G01N3/10—Investigating strength properties of solid materials by application of mechanical stress by applying steady tensile or compressive forces generated by pneumatic or hydraulic pressure
- G01N3/12—Pressure testing
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2203/00—Investigating strength properties of solid materials by application of mechanical stress
- G01N2203/003—Generation of the force
- G01N2203/0042—Pneumatic or hydraulic means
- G01N2203/0044—Pneumatic means
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2203/00—Investigating strength properties of solid materials by application of mechanical stress
- G01N2203/0058—Kind of property studied
- G01N2203/006—Crack, flaws, fracture or rupture
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)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
Abstract
The invention provides a method and a system for quantitatively analyzing fracture complexity after fracturing, and belongs to the field of shale gas fracturing. The method quantitatively evaluates the fracture complexity after fracturing through the wave crest number of GdP/dG curve, the maximum slope of GdP/dG curve first part and standard deviation of second part; the first part refers to a curve part which is included from an origin to a first peak and then to a first trough in an GdP/dG curve; the second part refers to the part of the GdP/dG curve from the first trough to the end point. The method can be used for quantitatively analyzing the complexity of the fractured cracks, can clearly know the fracturing modification degree of each section, and is beneficial to guiding fracturing design optimization, fracturing construction optimization and improving the modification effect.
Description
Technical Field
The invention belongs to the field of shale gas fracturing, and particularly relates to a method and a system for quantitatively analyzing fracture complexity after fracturing.
Background
At present, the technology for diagnosing and analyzing the fracture after hydraulic fracturing is widely applied to various oil and gas reservoirs such as shale, tight sandstone, carbonate rock and the like. The post-fracturing fracture diagnosis technology is mainly used for predicting fracture size, complexity, fracture extension conditions and the like. The shale gas well fracturing belongs to volume fracturing, and the core of the volume fracturing is to greatly improve the fracture complexity to increase the modification volume and improve the yield after fracturing, so that the fracture complexity analysis is very important for evaluating the shale gas well fracturing effect.
In the prior post-fracturing analysis, the fracture complexity is qualitatively analyzed by analyzing the fluctuation degree of a fracturing construction pressure curve and the fluctuation degree of a G function curve (the fluctuation degree of the G function is closely related to the post-fracturing fracture complexity), the method is helpful for judging the post-fracturing fracture complexity to a certain extent, but the misjudgment condition is easy to occur, because the apparent fluctuation degree of the curve has certain deceptiveness, the fluctuation degree difference shown in different data ranges is large; this also results in the inability to form a relatively constant fracture complexity evaluation system. Therefore, it is necessary to provide a method for quantitatively analyzing the complexity of the post-compression fracture to solve the above problems.
Chinese patent publication CN106769463A discloses a quantitative characterization method for fracture complexity after core pressing, which is mainly used for quantitatively evaluating the fracture complexity according to the inclination angle and the area of the fracture of the pressed core and belongs to a laboratory evaluation method for core size; the Chinese document 'quantitative characterization of shale volume fracturing complex fractures' (petroleum and natural gas geology, 2017, stage 01) discloses a quantitative characterization method of shale volume fracturing complex fractures, which quantitatively characterizes complex fractures after shale fracturing through 4 established typical fracture distribution forms and fracturing physical simulation experiments; the Chinese document ' evaluation research of complexity of fracture after fracturing based on G function curve analysis ' (' scientific technology and engineering ' 2016 (02 th year) ') discloses evaluation of complexity of fracture after fracturing, which analyzes the filtration behavior of fracturing fluid when a natural fracture is opened on the basis of the classic G function theory, corrects the filtration coefficient of the closed stage of the natural fracture by introducing a free variable omega and establishes a pressure drop model considering the opening of the natural fracture; the curve shapes of the first derivative dp/dG, the superimposed derivative ISIP-Gdp/dG and Gdp/dG for a simple main fracture, a natural fracture communicating with a plurality of natural fractures are compared and analyzed. The document corrects a G function theory, considers more influences of fluid loss characteristics under the opening of a natural crack on the G function, judges the complexity of the crack according to a corrected G function curve, and judges whether to communicate the natural crack and form a complex crack or only form a crack system mainly based on a main crack by simply using the fluctuation of the G function in vision (namely, qualitative rather than quantitative) when judging the complexity of the crack, but does not provide a method and an idea for quantitatively evaluating the fluctuation of the G function, and does not form a specific index.
Disclosure of Invention
The invention aims to solve the problems in the prior art and provides a method and a system for quantitatively analyzing the complexity of fractured fractures, so that the complexity of the fractured fractures is systematically known, the fracturing design and construction are guided to be optimized, and the efficient yield increase of shale gas wells is realized.
The invention is realized by the following technical scheme:
a method for quantitatively analyzing fracture complexity after fracturing is used for quantitatively evaluating the fracture complexity after fracturing through the wave crest number of GdP/dG curve, the maximum slope of a first part of GdP/dG curve and the standard deviation of a second part;
the first part refers to a curve part which is included from an origin to a first peak and then to a first trough in an GdP/dG curve;
the second part refers to the part of the GdP/dG curve from the first trough to the end point.
The method comprises the following steps:
(1) collecting second point data of fracturing construction, and drawing an GdP/dG curve graph according to the second point data of fracturing construction;
(2) adjusting GdP/dG curve chart to make its abscissa range between [0,0.1 ];
(3) obtaining GdP/dG peak number；
(4) Calculating GdP/dG curve maximum slope and second standard deviation;
(5) and obtaining fracture complexity according to the wave peak number, the maximum slope of the first part and the standard deviation of the second part.
The operation of calculating the maximum slope of the first part in the step (4) includes:
the ordinate H corresponding to the highest point of the first part (i.e. the highest point of the first peak) is read from the data of the GdP/dG curvemaxAnd the abscissa Tmax；
The maximum slope k of the first portion is calculated using the following equation:
the operation of calculating the standard deviation of the second part in the step (4) comprises:
(a) calculating an average of the ordinate of the second portionValue of
Wherein; x is the number ofiIs the ordinate of each data point in the second part, i ═ 1, 2, 3 · · n; n is the number of data points of the second portion;
(b) calculate the variance of the second part:
(c) calculate the standard deviation of the second part:
the operation of the step (5) comprises the following steps:
a1, if three judgment conditions of equal complexity are simultaneously satisfied, judging the system to be equal complexity;
a2, if any one of three judgment conditions of four equal complexity is satisfied, judging that the three judgment conditions are four equal complexity;
a3, if three judgment conditions of secondary complexity are simultaneously met, judging the secondary complexity; if one or two of the three judgment conditions with first-class complexity are met, and the other judgment conditions with second-class complexity are met, judging that the second-class complexity is achieved;
a4, if three judgment conditions of three equal complexities are simultaneously satisfied, judging that the three equal complexities are three equal complexities; if one or two of the three judgment conditions with equal complexity are met, and the other judgment conditions with equal complexity are met, judging the three judgment conditions to be equal complexity; if one of three judgment conditions with first-class complexity, one of three judgment conditions with second-class complexity and one of three judgment conditions with third-class complexity are met, judging that the three-class complexity is third-class complexity; if one or two of the judgment conditions with the second-class complexity are met, and the other judgment conditions with the third-class complexity are met, judging that the third-class complexity is achieved;
the three first-class complex judgment conditions are as follows: the number of peaks is more than 12, the maximum slope of the first part is more than 100000, and the standard deviation of the second part is more than 40;
the two complex three judgment conditions are: the number of peaks is greater than 9 and equal to or less than 12, the maximum slope of the first portion is greater than 60000 and equal to or less than 100000, and the standard deviation of the second portion is greater than 25 and equal to or less than 40;
three equally complicated three judgment conditions are: the number of peaks is greater than 5 and equal to or less than 9, the maximum slope of the first part is greater than 12000 and equal to or less than 60000, and the standard deviation of the second part is greater than or equal to 15 and equal to or less than 25;
four equally complicated three judgment conditions are: the number of peaks is not more than 5, the maximum slope of the first part is not more than 12000, and the standard deviation of the second part is not more than 15.
The invention also provides a system for quantitatively analyzing the complexity of the fractured cracks, which comprises the following components:
the curve drawing unit is used for drawing an GdP/dG curve according to the second point data of the fracturing construction;
the abscissa adjusting unit is connected with the graph drawing unit and is used for adjusting the range of the abscissa of the GdP/dG graph output by the graph drawing unit to be between [0,0.1 ];
the peak number counting unit is connected with the abscissa adjusting unit and is used for obtaining the peak number on the GdP/dG curve;
the maximum slope calculation unit is connected with the abscissa adjustment unit and is used for obtaining the maximum slope of the first part;
the standard deviation calculation unit is connected with the abscissa adjusting unit and used for obtaining the standard deviation of the second part;
and the crack complexity judging unit is respectively connected with the peak number counting unit, the maximum slope calculating unit and the standard deviation calculating unit and is used for judging the crack complexity according to the peak number, the maximum slope of the first part and the standard deviation of the second part.
The maximum slope calculating unit reads the ordinate H corresponding to the highest point of the first part from the data of the GdP/dG curve output by the abscissa adjusting unitmaxAnd the abscissa TmaxThe maximum slope k of the first portion is then calculated using the following equation:
the standard deviation calculating unit reads the ordinate of each data point of the second portion from the data of the GdP/dG curve output from the abscissa adjusting unit, and then calculates the standard deviation of the second portion using the following formula:
wherein; sigma is the standard deviation, S is the variance,
is the average of the ordinate of the second part, xiIs the ordinate of each data point in the second part, i ═ 1, 2, 3 · · n; n is the number of data points of the second portion.
The crack complexity determination unit is configured to:
if three judgment conditions of equal complexity are simultaneously met, judging the system to be equal complexity;
if any one of three judgment conditions with four equal complexities is met, judging that the three judgment conditions are four equal complexities;
if three judgment conditions of second-class complexity are simultaneously met, judging that the second-class complexity is second-class complexity; if one or two of the three judgment conditions with first-class complexity are met, and the other judgment conditions with second-class complexity are met, judging that the second-class complexity is achieved;
if three judging conditions with three equal complexities are simultaneously met, judging that the three equal complexities are achieved; if one or two of the three judgment conditions with equal complexity are met, and the other judgment conditions with equal complexity are met, judging the three judgment conditions to be equal complexity; if one of three judgment conditions with first-class complexity, one of three judgment conditions with second-class complexity and one of three judgment conditions with third-class complexity are met, judging that the three-class complexity is third-class complexity; and if one or two of the judgment conditions with the second-class complexity are met, and the other judgment conditions with the third-class complexity are met, judging that the system is three-class complex.
The present invention also provides a computer readable storage medium storing at least one program executable by a computer, the at least one program, when executed by the computer, causing the computer to perform the steps in the method of quantitatively analyzing post-fracture complexity of the present invention.
Compared with the prior art, the invention has the beneficial effects that: the method can be used for quantitatively analyzing the complexity of the fractured cracks, can clearly know the fracturing modification degree of each section, and is beneficial to guiding fracturing design optimization, fracturing construction optimization and improving the modification effect.
Drawings
FIG. 1 is a block diagram of the steps of the method of the present invention;
FIG. 2 is a graph of a G function according to an embodiment of the present invention;
FIG. 3 is a schematic diagram of the system of the present invention.
Detailed Description
The invention is described in further detail below with reference to the accompanying drawings:
1) determining GdP/dG function curve fluctuation times
Intuitively, the fluctuation degree of the GdP/dG curve (a G function curve) is closely related to the number of peaks (troughs) and the data oscillation amplitude in the GdP/dG curve (as shown in FIG. 2), and the two are intuitively closely related to the range of coordinate axes, and the data oscillation amplitude can be judged through quantitative data analysis. In order to judge the fluctuation times (namely, determine the peak number) of the GdP/dG curve on the same scale, the invention limits the range of the Nolte G time on the abscissa corresponding to the GdP/dG curve to [0,0.1] (as shown in FIG. 2), and then determines the fluctuation times (namely, the peak number) in a counting way.
2) GdP/dG function curve fluctuation degree evaluation
The GdP/dG function curve is divided into two parts, the first part is the curve part from the origin to the first peak to the first trough (i.e. including the first peak), and the remaining part (i.e. the curve part from the first trough back to the end) is the second part (also called the first peak back part). The first part is positioned in a short time after the pump is stopped, the speed of data change often reflects the connectivity between cracks and the extensibility of the cracks after the pump is stopped, the faster the data change, the better the connectivity and the extensibility of the cracks, and for quantitative evaluation, the slope size corresponding to the maximum value of the first peak is taken as the judgment standard.
The fluctuation degree of the second part of the curve often reflects the connectivity of the crack and the number of open and closed parts, and the stronger fluctuation degree indicates the higher complexity of the crack. The present invention evaluates the fluctuation degree by the standard deviation of the curve data, and the larger the standard deviation is, the stronger the fluctuation is.
3) Shale air pressure post-fracture complexity quantitative evaluation index
Based on the first and second points, the invention establishes a quantitative evaluation index of fracture complexity, as shown in table 1:
TABLE 1
In table 1, if only one of the three conditions, i.e., the number of peaks, the maximum slope of the first portion, and the standard deviation of the second portion, satisfies the current level, the current level is determined even if the other conditions satisfy the previous level, that is, if any one of the three conditions does not satisfy the current level, the next level is determined. The method comprises the following specific steps:
three equally complicated determination conditions are: the number of peaks is more than 12, the maximum slope of the first part is more than 100000, and the standard deviation of the second part is more than 40;
the two complex three judgment conditions are: the number of peaks is greater than 9 and equal to or less than 12, the maximum slope of the first portion is greater than 60000 and equal to or less than 100000, and the standard deviation of the second portion is greater than 25 and equal to or less than 40;
three equally complicated three judgment conditions are: the number of peaks is greater than 5 and equal to or less than 9, the maximum slope of the first part is greater than 12000 and equal to or less than 60000, and the standard deviation of the second part is greater than or equal to 15 and equal to or less than 25;
four equally complicated three judgment conditions are: the number of peaks is less than or equal to 5, the maximum slope of the first part is less than or equal to 12000, and the standard deviation of the second part is less than 15;
the judging steps are as follows:
a1, if three judgment conditions of equal complexity are simultaneously satisfied, judging the system to be equal complexity;
a2, if any one of three judgment conditions of four equal complexity is satisfied, judging that the three judgment conditions are four equal complexity;
a3, if three judgment conditions of secondary complexity are simultaneously met, judging the secondary complexity; if one or two of the three judgment conditions with first-class complexity are met, and the other judgment conditions with second-class complexity are met, judging that the second-class complexity is achieved;
a4, if three judgment conditions of three equal complexities are simultaneously satisfied, judging that the three equal complexities are three equal complexities; if one or two of the three judgment conditions with equal complexity are met, and the other judgment conditions with equal complexity are met, judging the three judgment conditions to be equal complexity; if one of three judgment conditions with first-class complexity, one of three judgment conditions with second-class complexity and one of three judgment conditions with third-class complexity are met, judging that the three-class complexity is third-class complexity; and if one or two of the judgment conditions with the second-class complexity are met, and the other judgment conditions with the third-class complexity are met, judging that the system is three-class complex.
The examples of the invention are as follows:
the method is applied to evaluation of the shale gas well pressure in the south of the Sichuan, and as shown in figure 1, the method comprises the following steps:
1. sorting fracturing construction second point data (data comprises wellhead pressure, pumping displacement, pumping time, proppant concentration and particle size, which are generally derived by an instrument truck during fracturing operation.) and drawing an GdP/dG curve graph by using Meyer fracturing design software (or other software and modes), as shown in FIG. 2;
2. the GdP/dG curve was adjusted to have an abscissa range of [0,0.1 ]: the curve is placed over the interval [0,0.1], with the left-most value (minimum) on the abscissa of the curve placed at 0, and the right-most value (maximum) on the abscissa placed at 0.1. (essentially, the curve is laterally compressed or laterally expanded).
3. The number of peaks on the GdP/dG curve is obtained, for example, the number of peaks can be directly counted, or the number of peaks on the curve can be obtained by using the existing method.
4. GdP/dG curve data are derived, and the maximum slope of the first part and the standard deviation of the second part are calculated:
the ordinate H corresponding to the highest point of the first part (i.e. the highest point of the first peak) is read from the data of the GdP/dG curvemaxAnd the abscissa TmaxThe maximum slope k of the first portion is calculated using the following equation:
the standard deviation of the second part was calculated using the following formula:
wherein; sigma is the standard deviation, S is the variance,
is the average of the ordinate of the second part, xiIs the ordinate of each data point in the second part, i ═ 1, 2, 3 · · n; n is the number of data points of the second portion.
5. The fracture complexity of the fractured section is evaluated according to the quantitative evaluation index of the fracture complexity shown in the table 1, and the result is shown in the table 2:
| number of fracturing stages | Number of peaks | Maximum slope of the first part | Standard deviation of the second part | |
| 1 | 3 | 8704.78 | 7.82 | Complexity of four equal parts |
| 2 | 2 | 2013.43 | 7.16 | Complexity of four equal parts |
| 5 | 1 | 1900.48 | 20.58 | Complexity of four equal parts |
| 7 | 13 | 63215.31 | 33.29 | Second degree of complexity |
| 8 | 13 | 13303.98 | 16.14 | Complexity of three equal parts |
| 9 | 3 | 5817.79 | 16.55 | Complexity of four equal parts |
| 10 | 13 | 31417.16 | 29.86 | Complexity of three equal parts |
| 11 | 12 | 44010.54 | 31.54 | Complexity of three equal parts |
| 12 | 3 | 5919.22 | 28.64 | Complexity of four equal parts |
| 13 | 4 | 9699.33 | 28.64 | Complexity of four equal parts |
| 15 | 4 | 15714.29 | 12.38 | Complexity of four equal parts |
TABLE 2
In the experiment, 15 fracturing stages are used in total, and each fracturing stage is evaluated independently, namely 15 GdP/dG curves are made, and then the fracture complexity of the 15 fracturing stages is evaluated respectively.
As shown in fig. 3, the present invention also provides a system for quantitatively analyzing fracture complexity after fracturing, comprising:
the graph drawing unit 100 is used for drawing an GdP/dG graph according to the second point data of fracturing construction;
an abscissa adjusting unit 200, connected to the graph plotting unit 100, for adjusting the range of the abscissa of the GdP/dG graph output by the graph plotting unit to [0,0.1 ];
a peak count unit 300, connected to the abscissa adjusting unit 200, for obtaining the number of peaks on the GdP/dG curve;
a maximum slope calculation unit 400 connected to the abscissa adjustment unit 200, for obtaining a maximum slope of the first portion;
a standard deviation calculation unit 500 connected to the abscissa adjustment unit 200 for obtaining a standard deviation of the second part;
the crack complexity determining unit 600 is connected to the peak count counting unit 300, the maximum slope calculating unit 400, and the standard deviation calculating unit 500, respectively, and is configured to determine the crack complexity according to the peak count, the maximum slope of the first portion, and the standard deviation of the second portion.
The maximum slope calculating unit 400 reads the ordinate H corresponding to the highest point of the first part from the data of the GdP/dG curve outputted from the abscissa adjusting unit 200maxAnd the abscissa TmaxThe maximum slope k of the first portion is then calculated using the following equation:
the standard deviation calculation unit 500 reads the ordinate of each data point of the second portion from the data of the GdP/dG curve output from the abscissa adjustment unit 200, and then calculates the standard deviation of the second portion using the following formula:
wherein; sigma is the standard deviation, S is the variance,
is the average of the ordinate of the second part, xiIs the ordinate of each data point in the second part, i ═ 1, 2, 3 · · n; n is the number of data points of the second portion.
The fracture complexity determination unit 600 is configured to:
if three judgment conditions of equal complexity are simultaneously met, judging the system to be equal complexity;
if any one of three judgment conditions with four equal complexities is met, judging that the three judgment conditions are four equal complexities;
if three judgment conditions of second-class complexity are simultaneously met, judging that the second-class complexity is second-class complexity; if one or two of the three judgment conditions with first-class complexity are met, and the other judgment conditions with second-class complexity are met, judging that the second-class complexity is achieved;
if three judging conditions with three equal complexities are simultaneously met, judging that the three equal complexities are achieved; if one or two of the three judgment conditions with equal complexity are met, and the other judgment conditions with equal complexity are met, judging the three judgment conditions to be equal complexity; if one of three judgment conditions with first-class complexity, one of three judgment conditions with second-class complexity and one of three judgment conditions with third-class complexity are met, judging that the three-class complexity is third-class complexity; and if one or two of the judgment conditions with the second-class complexity are met, and the other judgment conditions with the third-class complexity are met, judging that the system is three-class complex.
The present invention also provides a computer readable storage medium storing at least one program executable by a computer, the at least one program, when executed by the computer, causing the computer to perform the steps in the method of quantitatively analyzing post-fracture complexity of the present invention.
The invention analyzes the limitation of the evaluation aspect of the fracture complexity after shale gas pressure at present, and provides a quantitative analysis method for the fracture complexity. The method is based on a G function analysis curve, and three elements which can reflect the complexity of the crack most in the G function curve are preferably researched and selected: and the wave peak number, the maximum slope of the first wave peak and the standard deviation of the data after the first wave peak are calculated and analyzed to form a fracture complexity quantitative evaluation system. The method belongs to a method for diagnosing and quantitatively analyzing cracks after hydraulic fracturing. The method is mainly applied to the field of shale post-pressure evaluation, can clearly perform quantitative analysis on the complexity of the post-pressure fracture, is beneficial to knowing the fracture forming effect of fracturing modification, and is beneficial to guiding fracturing design and improving the modification effect.
The above-described embodiment is only one embodiment of the present invention, and it will be apparent to those skilled in the art that various modifications and variations can be easily made based on the application and principle of the present invention disclosed in the present application, and the present invention is not limited to the method described in the above-described embodiment of the present invention, so that the above-described embodiment is only preferred, and not restrictive.
Claims (10)
1. A method for quantitatively analyzing fracture complexity after fracturing is characterized by comprising the following steps: the method quantitatively evaluates the fracture complexity after fracturing through the number of wave peaks of GdP/dG curve, the maximum slope of the first part of GdP/dG curve and the standard deviation of the second part;
the first part refers to a curve part which is included from an origin to a first peak and then to a first trough in an GdP/dG curve;
the second part refers to the part of the GdP/dG curve from the first trough to the end point.
2. The method of quantitatively analyzing fracture complexity after fracturing of claim 1, wherein: the method comprises the following steps:
(1) collecting second point data of fracturing construction, and drawing an GdP/dG curve graph according to the second point data of fracturing construction;
(2) adjusting GdP/dG curve chart to make its abscissa range between [0,0.1 ];
(3) obtaining GdP/dG peak number；
(4) Calculating GdP/dG curve maximum slope and second standard deviation;
(5) and obtaining fracture complexity according to the wave peak number, the maximum slope of the first part and the standard deviation of the second part.
3. The method of quantitatively analyzing fracture complexity after fracturing of claim 2, wherein: the operation of calculating the maximum slope of the first part in the step (4) includes:
the ordinate H corresponding to the highest point in the first part is read from the data of the GdP/dG curvemaxAnd the abscissa Tmax；
The maximum slope k of the first portion is calculated using the following equation:
4. the method of quantitatively analyzing fracture complexity after fracturing of claim 2, wherein: the operation of calculating the standard deviation of the second part in the step (4) comprises:
(a) calculating an average of the ordinate of the second portion
Wherein; x is the number ofiIs the ordinate of each data point in the second part, i ═ 1, 2, 3 · · n; n is the number of data points of the second portion;
(b) calculate the variance of the second part:
(c) calculate the standard deviation of the second part:
5. the method of quantitatively analyzing fracture complexity after fracturing of claim 2, wherein: the operation of the step (5) comprises the following steps:
a1, if three judgment conditions of equal complexity are simultaneously satisfied, judging the system to be equal complexity;
a2, if any one of three judgment conditions of four equal complexity is satisfied, judging that the three judgment conditions are four equal complexity;
a3, if three judgment conditions of secondary complexity are simultaneously met, judging the secondary complexity; if one or two of the three judgment conditions with first-class complexity are met, and the other judgment conditions with second-class complexity are met, judging that the second-class complexity is achieved;
a4, if three judgment conditions of three equal complexities are simultaneously satisfied, judging that the three equal complexities are three equal complexities; if one or two of the three judgment conditions with equal complexity are met, and the other judgment conditions with equal complexity are met, judging the three judgment conditions to be equal complexity; if one of three judgment conditions with first-class complexity, one of three judgment conditions with second-class complexity and one of three judgment conditions with third-class complexity are met, judging that the three-class complexity is third-class complexity; if one or two of the judgment conditions with the second-class complexity are met, and the other judgment conditions with the third-class complexity are met, judging that the third-class complexity is achieved;
the three first-class complex judgment conditions are as follows: the number of peaks is more than 12, the maximum slope of the first part is more than 100000, and the standard deviation of the second part is more than 40;
the two complex three judgment conditions are: the number of peaks is greater than 9 and equal to or less than 12, the maximum slope of the first portion is greater than 60000 and equal to or less than 100000, and the standard deviation of the second portion is greater than 25 and equal to or less than 40;
three equally complicated three judgment conditions are: the number of peaks is greater than 5 and equal to or less than 9, the maximum slope of the first part is greater than 12000 and equal to or less than 60000, and the standard deviation of the second part is greater than or equal to 15 and equal to or less than 25;
four equally complicated three judgment conditions are: the number of peaks is not more than 5, the maximum slope of the first part is not more than 12000, and the standard deviation of the second part is not more than 15.
6. A system for quantitatively analyzing fracture complexity after fracturing, comprising: the system comprises:
the curve drawing unit is used for drawing an GdP/dG curve according to the second point data of the fracturing construction;
the abscissa adjusting unit is connected with the graph drawing unit and is used for adjusting the range of the abscissa of the GdP/dG graph output by the graph drawing unit to be between [0,0.1 ];
the peak number counting unit is connected with the abscissa adjusting unit and is used for obtaining the peak number on the GdP/dG curve;
the maximum slope calculation unit is connected with the abscissa adjustment unit and is used for obtaining the maximum slope of the first part;
the standard deviation calculation unit is connected with the abscissa adjusting unit and used for obtaining the standard deviation of the second part;
and the crack complexity judging unit is respectively connected with the peak number counting unit, the maximum slope calculating unit and the standard deviation calculating unit and is used for judging the crack complexity according to the peak number, the maximum slope of the first part and the standard deviation of the second part.
7. The system for quantitatively analyzing fracture complexity after fracturing of claim 6, wherein: the maximum slope calculating unit reads out a first part from the GdP/dG curve data outputted from the abscissa adjusting unitCorresponding to the highest point ofmaxAnd the abscissa TmaxThe maximum slope k of the first portion is then calculated using the following equation:
8. the system for quantitatively analyzing fracture complexity after fracturing of claim 6, wherein: the standard deviation calculating unit reads the ordinate of each data point of the second portion from the data of the GdP/dG curve output from the abscissa adjusting unit, and then calculates the standard deviation of the second portion using the following formula:
wherein; sigma is the standard deviation, S is the variance,
is the average of the ordinate of the second part, xiIs the ordinate of each data point in the second part, i ═ 1, 2, 3 · · n; n is the number of data points of the second portion.
9. The system for quantitatively analyzing fracture complexity after fracturing of claim 6, wherein: the crack complexity determination unit is configured to:
if three judgment conditions of equal complexity are simultaneously met, judging the system to be equal complexity;
if any one of three judgment conditions with four equal complexities is met, judging that the three judgment conditions are four equal complexities;
if three judgment conditions of second-class complexity are simultaneously met, judging that the second-class complexity is second-class complexity; if one or two of the three judgment conditions with first-class complexity are met, and the other judgment conditions with second-class complexity are met, judging that the second-class complexity is achieved;
if three judging conditions with three equal complexities are simultaneously met, judging that the three equal complexities are achieved; if one or two of the three judgment conditions with equal complexity are met, and the other judgment conditions with equal complexity are met, judging the three judgment conditions to be equal complexity; if one of three judgment conditions with first-class complexity, one of three judgment conditions with second-class complexity and one of three judgment conditions with third-class complexity are met, judging that the three-class complexity is third-class complexity; and if one or two of the judgment conditions with the second-class complexity are met, and the other judgment conditions with the third-class complexity are met, judging that the system is three-class complex.
10. A computer-readable storage medium characterized by: the computer readable storage medium stores at least one program executable by a computer, the at least one program when executed by the computer causing the computer to perform the steps in the method of quantitatively analyzing post-fracture complexity of any one of claims 1-5.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910634085.1A CN112304770B (en) | 2019-07-15 | 2019-07-15 | Method and system for quantitatively analyzing complexity of fracture after fracturing |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910634085.1A CN112304770B (en) | 2019-07-15 | 2019-07-15 | Method and system for quantitatively analyzing complexity of fracture after fracturing |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112304770A true CN112304770A (en) | 2021-02-02 |
| CN112304770B CN112304770B (en) | 2024-04-05 |
Family
ID=74328729
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910634085.1A Active CN112304770B (en) | 2019-07-15 | 2019-07-15 | Method and system for quantitatively analyzing complexity of fracture after fracturing |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112304770B (en) |
Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH0777498A (en) * | 1993-08-05 | 1995-03-20 | Hideki Matsuoka | Crack detecting method for object |
| CN2530038Y (en) * | 2002-01-25 | 2003-01-08 | 杨学政 | strong pulse gas energy fracturing device |
| JP2005338056A (en) * | 2004-05-26 | 2005-12-08 | Ching Yun Univ | Method of measuring crack in inside of reinforced concrete beam by p-wave penetration, and structure safety evaluation method |
| CN101726750A (en) * | 2009-11-30 | 2010-06-09 | 中国矿业大学(北京) | 'Three-map method' for evaluating ground fissures |
| CN106204497A (en) * | 2016-07-20 | 2016-12-07 | 长安大学 | A kind of pavement crack extraction algorithm based on smooth smoothed curve and matched curve |
| CN106769463A (en) * | 2016-12-20 | 2017-05-31 | 中国石油天然气集团公司 | Crack complexity quantitatively characterizing method after a kind of rock core pressure |
| US20170247995A1 (en) * | 2015-05-07 | 2017-08-31 | Baker Hughes Incorporated | Evaluating far field fracture complexity and optimizing fracture design in multi-well pad development |
| CN107702999A (en) * | 2017-11-13 | 2018-02-16 | 中国科学院地质与地球物理研究所 | Miniature hydraulic pressure-break net quantitative testing device and method |
| CN109426689A (en) * | 2017-08-31 | 2019-03-05 | 中国石油化工股份有限公司 | The evaluation method and system in fractured horizontal well crack |
-
2019
- 2019-07-15 CN CN201910634085.1A patent/CN112304770B/en active Active
Patent Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH0777498A (en) * | 1993-08-05 | 1995-03-20 | Hideki Matsuoka | Crack detecting method for object |
| CN2530038Y (en) * | 2002-01-25 | 2003-01-08 | 杨学政 | strong pulse gas energy fracturing device |
| JP2005338056A (en) * | 2004-05-26 | 2005-12-08 | Ching Yun Univ | Method of measuring crack in inside of reinforced concrete beam by p-wave penetration, and structure safety evaluation method |
| CN101726750A (en) * | 2009-11-30 | 2010-06-09 | 中国矿业大学(北京) | 'Three-map method' for evaluating ground fissures |
| US20170247995A1 (en) * | 2015-05-07 | 2017-08-31 | Baker Hughes Incorporated | Evaluating far field fracture complexity and optimizing fracture design in multi-well pad development |
| CN106204497A (en) * | 2016-07-20 | 2016-12-07 | 长安大学 | A kind of pavement crack extraction algorithm based on smooth smoothed curve and matched curve |
| CN106769463A (en) * | 2016-12-20 | 2017-05-31 | 中国石油天然气集团公司 | Crack complexity quantitatively characterizing method after a kind of rock core pressure |
| CN109426689A (en) * | 2017-08-31 | 2019-03-05 | 中国石油化工股份有限公司 | The evaluation method and system in fractured horizontal well crack |
| CN107702999A (en) * | 2017-11-13 | 2018-02-16 | 中国科学院地质与地球物理研究所 | Miniature hydraulic pressure-break net quantitative testing device and method |
Non-Patent Citations (6)
| Title |
|---|
| 何平;王建;陈柯伶;彭立才;张博;郭双生;: "巴喀油气田八道湾组气藏储集层裂缝定量分析", 新疆石油地质 * |
| 唐颖;邢云;李乐忠;张滨海;蒋时馨;: "页岩储层可压裂性影响因素及评价方法", 地学前缘, no. 05 * |
| 杨发荣;吴国铭;李熙;呤魃;刘华勋;吴国铨;: "孔洞型碳酸盐岩T_2谱分维值与微观结构参数关系探究", 科学技术与工程, no. 19 * |
| 赵文;张遂安;孙志宇;赵亚东;杨杨;: "基于G函数曲线分析的压后裂缝复杂性评估研究", 科学技术与工程, no. 33 * |
| 邢杨义;赵立强;刘平礼;龚云蕾;邢杨梅;: "裂缝性储层水平井压裂过程中起裂条件研究", 石油天然气学报, no. 05 * |
| 陈添;汪志明;杨刚;: "煤岩T型缝压裂实验及压力曲线分析", 特种油气藏, no. 03 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112304770B (en) | 2024-04-05 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107506948B (en) | Shale oil gas comprehensive yield analysis method based on dynamic drainage volume | |
| WO2020192675A1 (en) | Productivity prediction method for fractured horizontal well in tight oil reservoir | |
| CN104775810A (en) | Method for evaluating compressibility of shale gas reservoir | |
| CN111456709B (en) | Horizontal well multistage fracturing segmentation clustering method based on logging curve | |
| CN105631184A (en) | Oil and gas well production decline analysis method and system | |
| WO2021004442A1 (en) | Method for correcting crack fracturing parameter result for microseismic interpretation | |
| CN106337680A (en) | Method of determining recoverable reserves in oil reservoir production wells | |
| CN109977586B (en) | Sectional clustering method and device for volume fractured horizontal well | |
| CN109003248A (en) | A kind of characterizing method of pulveryte laminated texture | |
| CN110895629A (en) | Evaluation method for oil and gas reservoir exploitation degree | |
| CN108661628A (en) | One kind being based on the preferred engineering dessert quantitative evaluation method of parameter | |
| CN112304770A (en) | Method and system for quantitatively analyzing fracture complexity after fracturing | |
| CN113762614A (en) | Shale gas well productivity prediction method based on static oil and gas field data | |
| CN108664679A (en) | A kind of oil/gas well creation data analysis method | |
| CN109426689A (en) | The evaluation method and system in fractured horizontal well crack | |
| CN109736792B (en) | Natural gas reservoir and automatic identification method of fluid-containing property thereof | |
| CN105095552B (en) | A kind of common probability estimation of distribution parameters flow and analysis method | |
| CN113627640A (en) | Productivity well testing prediction method and system for fracture-cavity type oil reservoir oil and gas well | |
| CN115964961A (en) | Horizontal well water plugging and well selection decision-making method | |
| CN112836860B (en) | Method and system for determining dynamic yield of fractured well in whole period of yield decreasing stage | |
| CN116413831A (en) | Multi-well synthetic seismic record automatic calibration and construction interpretation and reservoir prediction method | |
| CN111894545B (en) | Method for determining proppant pumping scheme | |
| CN113627069A (en) | Well testing dynamic yield evaluation method and system for fracture-cavity type oil reservoir oil and gas well | |
| CN112145166A (en) | Underground condition identification and pre-judgment method in fracturing process of shale gas horizontal well | |
| CN113221232A (en) | Method for optimizing design of fracture sandstone thick-layer hydraulic fracturing construction parameters |
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 |