CN113221361B - Shale reservoir fracture width change calculation method under closing pressure - Google Patents
Shale reservoir fracture width change calculation method under closing pressure Download PDFInfo
- Publication number
- CN113221361B CN113221361B CN202110539939.5A CN202110539939A CN113221361B CN 113221361 B CN113221361 B CN 113221361B CN 202110539939 A CN202110539939 A CN 202110539939A CN 113221361 B CN113221361 B CN 113221361B
- Authority
- CN
- China
- Prior art keywords
- layer
- propping agent
- crack
- layers
- pressure
- 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.)
- Active
Links
- 230000008859 change Effects 0.000 title claims abstract description 25
- 238000004364 calculation method Methods 0.000 title claims abstract description 8
- 238000000034 method Methods 0.000 claims abstract description 45
- 238000007906 compression Methods 0.000 claims abstract description 36
- 230000006835 compression Effects 0.000 claims abstract description 35
- 239000002245 particle Substances 0.000 claims abstract description 17
- 238000001125 extrusion Methods 0.000 claims abstract description 11
- 239000010410 layer Substances 0.000 claims description 99
- 238000004458 analytical method Methods 0.000 claims description 6
- 239000011148 porous material Substances 0.000 claims description 5
- 239000002356 single layer Substances 0.000 claims description 3
- 239000003795 chemical substances by application Substances 0.000 abstract description 48
- 230000008569 process Effects 0.000 abstract description 7
- 230000009471 action Effects 0.000 abstract description 4
- 238000002347 injection Methods 0.000 abstract description 3
- 239000007924 injection Substances 0.000 abstract description 3
- 238000004590 computer program Methods 0.000 description 7
- 238000010586 diagram Methods 0.000 description 6
- 238000005056 compaction Methods 0.000 description 4
- 230000000694 effects Effects 0.000 description 4
- 230000006870 function Effects 0.000 description 4
- 238000012545 processing Methods 0.000 description 4
- 238000011161 development Methods 0.000 description 3
- 239000012530 fluid Substances 0.000 description 3
- 238000004519 manufacturing process Methods 0.000 description 2
- 239000011435 rock Substances 0.000 description 2
- 239000000919 ceramic Substances 0.000 description 1
- 230000005489 elastic deformation Effects 0.000 description 1
- 238000002091 elastography Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 239000004576 sand Substances 0.000 description 1
- 238000007789 sealing Methods 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 238000010206 sensitivity analysis Methods 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
- E21B49/00—Testing the nature of borehole walls; Formation testing; Methods or apparatus for obtaining samples of soil or well fluids, specially adapted to earth drilling or wells
-
- 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
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A10/00—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE at coastal zones; at river basins
- Y02A10/40—Controlling or monitoring, e.g. of flood or hurricane; Forecasting, e.g. risk assessment or mapping
Landscapes
- 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 provides a calculation method for the width change of a crack under shale reservoir closing pressure, which comprises the following steps: radius R of proppant particle acquisition 1 Modulus of elasticity E 1 Poisson ratio v 1 The method comprises the steps of carrying out a first treatment on the surface of the And the initial width L of the crack to be filled, the modulus of elasticity E of the crack wall 2 Poisson ratio v 2 Acquiring acquisition data; respectively calculating extrusion forces of the propping agent A layer and the propping agent B layer under closed pressure according to acquired data and a Hertz contact theory; respectively calculating total pressure of the two layers A and the plurality of layers B according to a self-defined extrusion force formula; and calculating the obtained total pressure of the two layers A and the plurality of layers B to obtain the crack width variation. Based on Hertz contact theory, the method establishes a model of the propping action of the propping agent in the fracture under the fracture closing pressure, and compared with the existing model, the model considers the injection of the propping agent and the elastic compression between propping agents, and is closer to the actual physical process; the accuracy of the obtained crack width variation is higher.
Description
Technical Field
The invention relates to a fracture width change calculation method, in particular to a fracture width change calculation method under shale reservoir closing pressure.
Background
The description of the width of the crack is the key for evaluating the flow conductivity of the crack, and influences the development efficiency of unconventional oil and gas resources. Typically, when proppants are filled into a fracture, the fracture will undergo a change in width at the closing pressure, and the proppants resist the change in fracture width. Therefore, it is of great importance to build a theoretical model for predicting crack changes. The hydraulic fracturing technology is common in unconventional oil and gas reservoir development, a fracture network consisting of natural fractures and induced fractures is formed after the hydraulic fracturing process, a fluid seepage channel is provided, and propping agents are injected into the fracture network along with fracturing fluid; once the well is opened and put into production, the pore pressure can be obviously reduced, the width of the crack can be reduced, part of propping agent can be discharged back, and part of propping agent stays in the crack, and the propping effect is excited along with the change of the width of the crack, so that the opening degree of the crack is maintained; propping action of proppants under closing pressure can affect fracture width and conductivity, related to development efficiency. Therefore, it is of great importance to study the pressure resistance of proppants, describing the variation of fracture width under closure pressure.
In fact, during the compression process, proppants in the fracture may be compressed or embedded, resulting in a change in fracture width; the relation between the embedding degree of proppants and the closing pressure, the concentration of proppants and the rock mechanical properties is studied by students in 1998 as early as possible; in the aspect of theoretical research, a mathematical model of the proppant injection amount under the sealing pressure is established through analysis, and the sensitivity of the proppant particle size and the rock Young modulus is discussed; deducing a proppant embedding depth model by using an elasticity theory, and carrying out sensitivity analysis on basic parameters; the influence of the particle size and the layer number of the propping agent on the fracture conductivity is characterized by adopting a numerical method.
These studies have discussed the impact factors on the proppant's embedding process and fracture conductivity, but most are based on rigid sphere assumptions, ignoring contact deformation between proppants. The variation in slit width at closing pressure cannot be accurately described. In particular, a theoretical model that takes into account proppant elastography needs to be built.
Disclosure of Invention
The invention aims to provide a calculation method for the width change of a crack under the closing pressure of a shale reservoir, which can be used for calculating the width change of the crack under the closing pressure of the shale reservoir.
The invention aims at realizing the technical scheme, which comprises the following steps:
the method for calculating the width change of the crack of the shale reservoir under the closing pressure comprises the following specific steps of:
1) Radius R of proppant particle acquisition 1 Modulus of elasticity E 1 Poisson ratio v 1 The method comprises the steps of carrying out a first treatment on the surface of the And the initial width L of the crack to be filled, the modulus of elasticity E of the crack wall 2 Poisson ratio v 2 Acquiring acquisition data;
2) Respectively calculating extrusion forces of the propping agent A layer and the propping agent B layer under the closing pressure according to the acquired data of the step 1) and the Hertz contact theory; respectively calculating total pressure of the two layers A and the plurality of layers B according to a self-defined extrusion force formula;
3) And (3) calculating the variation of the width of the crack according to the total pressure of the two layers A and the plurality of layers B obtained in the step (2).
Further, the specific method for calculating the total compression of the propping agent and the fracture wall contact layer A layer in the step 2) is as follows:
the propping agent is directly contacted with the crack wall to form a layer A, and the normal force N of the propping agent of the single layer A is as follows according to the Hertz contact theory:
in the formula (1), R 1 For radius of propping agent, R 2 Radius of crack wall, E 1 Modulus of elasticity, E, of the proppant 2 Modulus of elasticity, v, of the fracture wall 1 Poisson's ratio, v as proppant 2 Poisson ratio, alpha, of fracture wall 1 Is the compression value of layer a;
radius R because the fracture wall tends to be planar 2 → infinity, formula (1) can be simplified as:
the proppants are closely arranged in multiple layers, each propping agent is extruded by surrounding 8 particles, the contact area formed by extrusion of each particle is approximately a regular quadrangle, wherein any side length is the particle diameter of 2r, and the area of a single triangle is as follows:
total contact area:
the proppants are elastic spheres, each of which is subjected to a pressure F equal to its normal force N for the a-layer proppants:
in formula (5), F is the pressure of the proppant, P is the closing pressure, P c For crack closure pressure, R 1 Is the proppant diameter. Based on force analysis, crack closure pressure P c The method comprises the following steps:
P c =P o -P p (6)
in formula (6), P o To overburden pressure, P p Is pore pressure;
compression value alpha of the a-layer proppant 1 The method comprises the following steps:
the upper and lower layers of propping agent are respectively contacted with the top layer and the bottom layer of the crack wall, so the total compression amount h caused by the propping agent of the layer A 1 The method comprises the following steps:
further, the specific method for calculating the total compression of the proppants and the proppants contact layer B in the step 2) is as follows:
the propping agent and propping agent contact layer is a B layer, and according to the Hertz contact theory, the normal force N' of the propping agent of the single B layer is as follows:
in the formula (9), R 1 Radius of propping agent, E 1 As the elastic modulus of the proppant, v 1 Poisson ratio, alpha, of proppant 2 A compression value for a single B layer;
the pressure F' to which the single proppants of the layer B are subjected by stress analysis is as follows:
compression value alpha of single B layer proppant 2 The method comprises the following steps:
according to the fracture width L and the proppant radius R 1 And calculating the total number of layers of the propping agent as n, wherein the total compression amount of the layer B is as follows:
further, the specific method for calculating the crack width change in the step 2) is as follows:
the fracture width change h is the sum of the total compression of the propping agent and the layer A and the total compression of the layer B:
due to the adoption of the technical scheme, the invention has the following advantages:
1. based on Hertz contact theory, the model of the propping action of the propping agent in the fracture under the fracture closing pressure is established, and compared with the existing model, the model considers the injection of the propping agent and the elastic compression between propping agents, and is closer to the actual physical process.
2. According to the method, the fracture width variation is obtained by calculating the compression amount of the propping agent and the fracture contact layer and the elastic compression between propping agents, and the accuracy of the obtained fracture width variation is higher.
Additional advantages, objects, and features of the invention will be set forth in part in the description which follows and in part will become apparent to those having ordinary skill in the art upon examination of the following or may be learned from practice of the invention. The objectives and other advantages of the invention will be realized and attained by the structure particularly pointed out in the written description and claims hereof.
Drawings
The drawings of the present invention are described below.
Fig. 1 is a schematic structural view of the present invention.
Fig. 2 is a compression model of proppant in a fracture at closure pressure.
Fig. 3 shows the predicted result of the crack width change.
Detailed Description
The invention is further described below with reference to the drawings and examples.
The method for calculating the width change of the crack under the closing pressure of the shale reservoir is shown in (a) - (c) in fig. 2, proppants are injected into the crack of the shale reservoir, the proppants in the crack are spherical, the proppants are mutually connected and distributed, the compaction effect of the closing pressure on the proppants is particle elastic deformation, the proppants are piled up in (a) in fig. 2, the proppants are deformed under the compaction effect, the proppants in (B) in fig. 2 are models after being compressed, the proppants in (c) in fig. 2 are models after being compressed, the proppants comprise two single contact layers of upper proppants and lower proppants and crack wall contact layers and a plurality of double contact layers of middle proppants, the single contact layer is defined as a layer A, the double contact layers are a layer B, the material performance of the proppants is kept stable in the deformation process, the compaction effect is contributed by the two parts of the layer A and the layer B, the stable shape of the proppants after being compressed is hexagonal, and the method for calculating the width change of the crack specifically comprises the following steps:
1) Radius R of proppant particle acquisition 1 Modulus of elasticity E 1 Poisson ratio v 1 The method comprises the steps of carrying out a first treatment on the surface of the And the initial width L of the crack to be filled, the modulus of elasticity E of the crack wall 2 Poisson ratio v 2 Acquiring acquisition data;
2) Respectively calculating extrusion forces of the propping agent A layer and the propping agent B layer under the closing pressure according to the acquired data of the step 1) and the Hertz contact theory; respectively calculating total pressure of the two layers A and the plurality of layers B according to a self-defined extrusion force formula;
3) And (3) calculating the variation of the width of the crack according to the total pressure of the two layers A and the plurality of layers B obtained in the step (2).
As an embodiment of the present invention, the specific method for calculating the total compression amount of the proppants and the fracture wall contact layer a layer in the step 2) is as follows:
the propping agent is directly contacted with the crack wall to form a layer A, and the normal force N of the propping agent of the single layer A is as follows according to the Hertz contact theory:
in the formula (14), R 1 For radius of propping agent, R 2 Radius of crack wall, E 1 Modulus of elasticity, E, of the proppant 2 Modulus of elasticity, v, of the fracture wall 1 Poisson's ratio, v as proppant 2 Poisson ratio, alpha, of fracture wall 1 Is the compression value of layer a;
radius R because the fracture wall tends to be planar 2 → infinity, formula (1) can be simplified as:
the proppants are closely arranged in multiple layers, each propping agent is extruded by surrounding 8 particles, the contact area formed by extrusion of each particle is approximately a regular quadrangle, wherein any side length is the particle diameter of 2r, and the area of a single triangle is as follows:
total contact area:
the proppants are elastic spheres, each of which is subjected to a pressure F equal to its normal force N for the a-layer proppants:
in the formula (18), F is the pressure of the propping agent, P is the closing pressure, P c For crack closure pressure, R 1 Is the proppant diameter;
under the action of fracture closure, the propping agent is extruded and deformed, the influence of pressure on the compressed propping agent is balanced, and the fracture wall is subjected to the closure pressure P c Influence of a value equal to the overburden pressure P o And pore pressure P p The difference between the compression force N and the closing pressure P of the propping agent c Maintaining balance, wherein gravity and buoyancy of proppants in the fracture are ignored in the compaction process; fracture closure pressure P c The method comprises the following steps:
P c =P o -P p (19)
in the formula (19), P o To overburden pressure, P p Is pore pressure;
compression value alpha of the a-layer proppant 1 The method comprises the following steps:
the upper and lower layers of propping agent are respectively contacted with the top layer and the bottom layer of the crack wall, so the total compression amount h caused by the propping agent of the layer A 1 The method comprises the following steps:
as an embodiment of the present invention, the specific method for calculating the total compression of the proppant and the proppant contacting layer B layer in step 2) is:
the propping agent and propping agent contact layer is a B layer, and according to the Hertz contact theory, the normal force N' of the propping agent of the single B layer is as follows:
in the formula (22), R 1 Radius of propping agent, E 1 As the elastic modulus of the proppant, v 1 Poisson ratio, alpha, of proppant 2 A compression value for a single B layer;
as shown in fig. 2 (d), the pressure F' to which the single proppant of layer B is subjected by stress analysis is:
compression value alpha of single B layer proppant 2 The method comprises the following steps:
according to the fracture width L and the proppant radius R 1 And calculating the total number of layers of the propping agent as n, wherein the total compression amount of the layer B is as follows:
as an embodiment of the present invention, the specific method for calculating the crack width change in step 2) is as follows:
the fracture width change h is the sum of the total compression of the propping agent and the layer A and the total compression of the layer B:
as shown in fig. 3, the ceramic proppant is selected, the elastic modulus is 11306MPa, and the poisson ratio is 0.2. The average diameters of the proppants were 0.55mm (20-40 mesh), 0.27mm (40-60 mesh) and 0.17mm (70-100 mesh), respectively. The modulus of elasticity and poisson's ratio of the fracture matrix were 0.25 and 8000MPa, respectively. The closing pressure is 0-70 MPa. By the above method, calculation of the crack width change is performed, and the predicted structure is shown in fig. 3 (Theoretical model).
In order to study the influence factors of the hydraulic support fracture conductivity of the shale reservoir, the factors of the type of propping agent, the grain size, the sand content, the circulating stress, the fracturing fluid backflow and the like are considered in the document Liu Xuewei, namely the shale reservoir hydraulic support fracture conductivity influence factor [ J ]. The broken block oil-gas field 2020,27 (3): 394-398', and an indoor experiment is carried out. The variation of the fracture width under closing pressure was studied using proppants of sizes 0.55mm (20-40 mesh), 0.27mm (40-60 mesh) and 0.17mm (70-100 mesh), as shown in fig. 3 (Liu, 2020).
As can be seen from fig. 3, the variation in fracture width predicted by the methods described herein at different closure pressures for different particle sizes of ceramsite proppants is consistent with experimental data (Liu, 2020). Therefore, the model provided herein is accurate, conforms to the physical mechanism of propping of the proppants in the fracture, and facilitates the prediction of fracture width under closed pressure conditions.
It will be appreciated by those skilled in the art that embodiments of the present application may be provided as a method, system, or computer program product. Accordingly, the present application may take the form of an entirely hardware embodiment, an entirely software embodiment, or an embodiment combining software and hardware aspects. Furthermore, the present application may take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, and the like) having computer-usable program code embodied therein.
The present application is described with reference to flowchart illustrations and/or block diagrams of methods, apparatus (systems) and computer program products according to embodiments of the application. It will be understood that each flow and/or block of the flowchart illustrations and/or block diagrams, and combinations of flows and/or blocks in the flowchart illustrations and/or block diagrams, can be implemented by computer program instructions. These computer program instructions may be provided to a processor of a general purpose computer, special purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be stored in a computer-readable memory that can direct a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer-readable memory produce an article of manufacture including instruction means which implement the function specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
Finally, it should be noted that: the above embodiments are only for illustrating the technical aspects of the present invention and not for limiting the same, and although the present invention has been described in detail with reference to the above embodiments, it should be understood by those of ordinary skill in the art that: modifications and equivalents may be made to the specific embodiments of the invention without departing from the spirit and scope of the invention, which is intended to be covered by the claims.
Claims (2)
1. The method for calculating the width change of the crack of the shale reservoir under the closing pressure of the shale reservoir is characterized in that the propping agent comprises two single contact layers of an upper propping agent and a lower propping agent and a crack wall contact layer and a plurality of double contact layers of an intermediate propping agent and the propping agent, wherein the single contact layer is defined as an A layer, the double contact layer is defined as a B layer, and the method for calculating the width change of the crack comprises the following specific steps:
1) Radius R of proppant particle acquisition 1 Modulus of elasticity E 1 Poisson ratio v 1 The method comprises the steps of carrying out a first treatment on the surface of the And the initial width L of the crack to be filled, the modulus of elasticity E of the crack wall 2 Poisson ratio v 2 Acquiring acquisition data;
2) Respectively calculating extrusion forces of the propping agent A layer and the propping agent B layer under the closing pressure according to the acquired data of the step 1) and the Hertz contact theory; respectively calculating total pressure of the two layers A and the plurality of layers B according to a self-defined extrusion force formula;
3) Obtaining the crack width variation according to the total pressure calculation of the two layers A and the plurality of layers B obtained in the step 2);
the specific method for calculating the total compression of the two A layers in the step 2) is as follows:
the propping agent is directly contacted with the crack wall to form a layer A, and the normal force N of the propping agent of the single layer A is as follows according to the Hertz contact theory:
in the formula (1), R 1 For radius of propping agent, R 2 Radius of crack wall, E 1 Modulus of elasticity, E, of the proppant 2 Modulus of elasticity, v, of the fracture wall 1 Poisson's ratio, v as proppant 2 Poisson ratio, alpha, of fracture wall 1 Is layer AIs a compressed value of (2);
radius R because the fracture wall tends to be planar 2 → infinity, formula (1) can be simplified as:
the proppants are closely arranged in multiple layers, each propping agent is extruded by surrounding 8 particles, the contact area formed by extrusion of each particle is approximately a regular quadrangle, wherein any side length is the particle diameter of 2r, and the area of a single triangle is as follows:
total contact area:
the proppants are elastic spheres, each of which is subjected to a pressure F equal to its normal force N for the a-layer proppants:
in formula (5), F is the pressure of the proppant, P is the closing pressure, P c For crack closure pressure, R 1 Is the proppant diameter; based on force analysis, crack closure pressure P c The method comprises the following steps:
P c =P o -P p (6)
in formula (6), P o To overburden pressure, P p Is pore pressure;
compression value alpha of the a-layer proppant 1 The method comprises the following steps:
the upper and lower layers of propping agent are respectively contacted with the top layer and the bottom layer of the crack wall, so the total compression amount h caused by the propping agent of the layer A 1 The method comprises the following steps:
the specific method for calculating the total compression of a plurality of B layers in the step 2) is as follows:
the propping agent and propping agent contact layer is a B layer, and according to the Hertz contact theory, the normal force N' of the propping agent of the single B layer is as follows:
in the formula (9), R 1 Radius of propping agent, E 1 As the elastic modulus of the proppant, v 1 Poisson ratio, alpha, of proppant 2 A compression value for a single B layer;
the pressure F' to which the single proppants of the layer B are subjected by stress analysis is as follows:
compression value alpha of single B layer proppant 2 The method comprises the following steps:
according to the fracture width L and the proppant radius R 1 And calculating the total number of layers of the propping agent as n, wherein the total compression amount of the layer B is as follows:
2. the method for calculating the width change of the crack under the closing pressure of the shale reservoir according to claim 1, wherein the specific method for calculating the width change of the crack in the step 3) is as follows:
the fracture width change h is the sum of the total compression of the propping agent and the layer A and the total compression of the layer B:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110539939.5A CN113221361B (en) | 2021-05-18 | 2021-05-18 | Shale reservoir fracture width change calculation method under closing pressure |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110539939.5A CN113221361B (en) | 2021-05-18 | 2021-05-18 | Shale reservoir fracture width change calculation method under closing pressure |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113221361A CN113221361A (en) | 2021-08-06 |
CN113221361B true CN113221361B (en) | 2023-06-23 |
Family
ID=77092638
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110539939.5A Active CN113221361B (en) | 2021-05-18 | 2021-05-18 | Shale reservoir fracture width change calculation method under closing pressure |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113221361B (en) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113591328B (en) * | 2021-08-24 | 2022-04-12 | 西南石油大学 | Deep shale hydraulic fracture long-term support opening change calculation method |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110443933A (en) * | 2019-09-06 | 2019-11-12 | 重庆三峡学院 | One kind carrying out detection device and method to number based on pressure source |
CN111720105A (en) * | 2020-08-04 | 2020-09-29 | 西南石油大学 | Diversion capacity prediction method considering particle size distribution of proppant |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10935684B2 (en) * | 2015-10-28 | 2021-03-02 | Halliburton Energy Services, Inc. | Near real-time return-on-fracturing-investment optimization for fracturing shale and tight reservoirs |
-
2021
- 2021-05-18 CN CN202110539939.5A patent/CN113221361B/en active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110443933A (en) * | 2019-09-06 | 2019-11-12 | 重庆三峡学院 | One kind carrying out detection device and method to number based on pressure source |
CN111720105A (en) * | 2020-08-04 | 2020-09-29 | 西南石油大学 | Diversion capacity prediction method considering particle size distribution of proppant |
Non-Patent Citations (3)
Title |
---|
考虑支撑剂及其嵌入程度对支撑裂缝导流能力影响的数值模拟;吴国涛;胥云;杨振周;杨立峰;张静;;天然气工业(第05期);72-75 * |
考虑支撑剂变形的压后支撑缝宽预测新模型;李勇明;程垒明;周文武;;科学技术与工程(第06期);112-118 * |
裂缝型储层渗透率应力敏感评价及对生产影响;曹耐;雷刚;董平川;李虹;沈茂丁;吴亚红;;断块油气田(第02期);63-67 * |
Also Published As
Publication number | Publication date |
---|---|
CN113221361A (en) | 2021-08-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110410054B (en) | Prediction method for time-space flow conductivity of fracturing fracture of coal-bed gas well | |
Lacy et al. | Fracture width and embedment testing in soft reservoir sandstone | |
CN109902918B (en) | Shale fracturing self-supporting fracture network permeability calculation method | |
CN114048695B (en) | Effective shale gas seam net volume inversion method based on flowback data | |
Fan et al. | Experimental and numerical investigations of the role of proppant embedment on fracture conductivity in narrow fractures (includes associated errata) | |
Wang et al. | Numerical investigation of the effects of proppant embedment on fracture permeability and well production in Queensland coal seam gas reservoirs | |
CN113221361B (en) | Shale reservoir fracture width change calculation method under closing pressure | |
Xu et al. | Effectiveness and time variation of induced fracture volume: Lessons from water flowback analysis | |
CN113530536B (en) | Method and system for evaluating efficiency of tight sandstone gas reservoir horizontal well fracture reservoir | |
Ju et al. | Adaptive finite element-discrete element method for numerical analysis of the multistage hydrofracturing of horizontal wells in tight reservoirs considering pre-existing fractures, hydromechanical coupling, and leak-off effects | |
Usui et al. | Effect of poroelasticity on hydraulic fracture interactions | |
CN111950112B (en) | Carbonate reservoir dynamic analysis method suitable for bottom closure | |
Xu et al. | Estimating compressibility of complex fracture networks in unconventional reservoirs | |
Yu et al. | Impact of asymmetric stimulated rock volume on casing deformation in multi-stage fracturing; a case study | |
CN106199712B (en) | A kind of method and device in definite pressure break casing deformation region | |
Geng et al. | Experimental and numerical investigation of Non-Darcy flow in propped hydraulic fractures: Identification and characterization | |
Zhu et al. | A fracture conductivity model for channel fracturing based on lattice-Boltzmann-method and computational-fluid-dynamics | |
Tian et al. | Quantifying micro-proppants crushing rate and evaluating propped micro-fractures | |
CN107193042B (en) | Man-made fracture modeling method and device, computer storage medium based on pressure break phase | |
CN112035993B (en) | Bottom constant-pressure carbonate reservoir testing and evaluating method | |
Ren et al. | Shale gas effective fracture network volume prediction and analysis based on flow back data: A case study of southern Sichuan Basin shale | |
CN106097133A (en) | A kind of coal seam water content and aquifer yield Forecasting Methodology | |
CN111950111B (en) | Dynamic analysis method suitable for open-bottom carbonate reservoir | |
Khan et al. | Hydraulic fracture conductivity in shale reservoirs | |
Hu et al. | Numerical simulation study on formation mechanism of void distributions of proppant particles |
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 |