CA2808858C - Wellbore real-time monitoring and analysis of fracture contribution - Google Patents
Wellbore real-time monitoring and analysis of fracture contribution Download PDFInfo
- Publication number
- CA2808858C CA2808858C CA2808858A CA2808858A CA2808858C CA 2808858 C CA2808858 C CA 2808858C CA 2808858 A CA2808858 A CA 2808858A CA 2808858 A CA2808858 A CA 2808858A CA 2808858 C CA2808858 C CA 2808858C
- Authority
- CA
- Canada
- Prior art keywords
- fractures
- time
- fractured intervals
- determining
- ratio
- 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
- 238000012544 monitoring process Methods 0.000 title abstract description 7
- 238000004458 analytical method Methods 0.000 title abstract description 5
- 238000004519 manufacturing process Methods 0.000 claims abstract description 77
- 238000000034 method Methods 0.000 claims abstract description 27
- 238000009826 distribution Methods 0.000 claims abstract description 20
- 229930195733 hydrocarbon Natural products 0.000 claims description 22
- 150000002430 hydrocarbons Chemical class 0.000 claims description 22
- 238000009530 blood pressure measurement Methods 0.000 claims description 12
- 230000001052 transient effect Effects 0.000 claims description 11
- 238000012545 processing Methods 0.000 claims description 6
- 238000005259 measurement Methods 0.000 abstract description 23
- 206010017076 Fracture Diseases 0.000 description 109
- 208000010392 Bone Fractures Diseases 0.000 description 39
- 239000004215 Carbon black (E152) Substances 0.000 description 13
- 238000005070 sampling Methods 0.000 description 13
- 239000012530 fluid Substances 0.000 description 9
- 230000035699 permeability Effects 0.000 description 5
- 230000008859 change Effects 0.000 description 4
- 238000012360 testing method Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 238000005553 drilling Methods 0.000 description 3
- 230000003287 optical effect Effects 0.000 description 3
- 208000006670 Multiple fractures Diseases 0.000 description 2
- 238000013459 approach Methods 0.000 description 2
- 238000009529 body temperature measurement Methods 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 238000004891 communication Methods 0.000 description 2
- 239000000203 mixture Substances 0.000 description 2
- 230000000638 stimulation Effects 0.000 description 2
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 2
- 238000013500 data storage Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 239000000835 fiber Substances 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000012552 review Methods 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 238000010200 validation analysis Methods 0.000 description 1
- 230000003442 weekly effect Effects 0.000 description 1
Classifications
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH DRILLING; MINING
- E21B—EARTH DRILLING, e.g. DEEP 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
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH DRILLING; MINING
- E21B—EARTH DRILLING, e.g. DEEP DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B47/00—Survey of boreholes or wells
- E21B47/10—Locating fluid leaks, intrusions or movements
- E21B47/103—Locating fluid leaks, intrusions or movements using thermal measurements
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH DRILLING; MINING
- E21B—EARTH DRILLING, e.g. DEEP DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B43/00—Methods or apparatus for obtaining oil, gas, water, soluble or meltable materials or a slurry of minerals from wells
- E21B43/25—Methods for stimulating production
- E21B43/26—Methods for stimulating production by forming crevices or fractures
Abstract
Methods and apparatus are provided for calculating production of each of a plurality of fractured intervals (or fractures) and monitoring changes in the fracture contribution with time. Such real-time monitoring and analysis may be performed by combining temperature distribution (and pressure) measurements, a real-time surface multiphase flow measurement, and an inflow model for each fractured interval (or fracture). In this manner, the industry may be able to understand the behavior of fractures and, in turn, optimize the number of stages (i.e., fractured intervals), the number of fractures, and the spacing between fractures and stages.
Description
. .
=
WELLBORE REAL-TIME MONITORING AND
ANALYSIS OF FRACTURE CONTRIBUTION
BACKGROUND OF THE INVENTION
Field of the Invention Embodiments of the present invention generally relate to hydrocarbon production and, more particularly, to determining the individual contribution of fractured intervals (or fractures) in time.
Description of the Related Art Various tools may be used in order to measure the contribution of the fractures within wellbores. Different services companies may run production logging tools, and chemical tracers may also be used to determine the fracture contribution.
However, these measurements may only provide a snapshot of what is happening at the moment the measurements are performed, and may change with time because conditions within the wellbore are transient.
SUMMARY OF THE INVENTION
Embodiments of the invention generally relate to allocating production of each of a plurality of fractured intervals (or fractures). This allocation may be performed by combining temperature distribution (and pressure) measurements, a real-time surface multiphase flow measurement, and an inflow model for each fractured interval (or fracture).
One embodiment of the invention is a method for determining production of hydrocarbons. The method generally includes determining a temperature distribution associated with a plurality of fractured intervals or fractures disposed along a well;
measuring a total flow rate for the well; modeling an inflow rate for each of the plurality of fractured intervals or fractures; and allocating production of each of the plurality of fractured intervals or fractures based on the temperature distribution, the total flow rate, and the inflow rates.
Another embodiment of the invention provides a system for determining production of hydrocarbons. The system generally includes a temperature sensing device configured to determine a temperature distribution associated with a plurality of fractured intervals or fractures disposed along a well, a flowmeter configured to measure a total flow rate for the well, and a processing unit. The processing unit is typically configured to model an inflow rate for each of the plurality of fractured intervals or fractures and to allocate production of each of the plurality of fractured intervals or fractures based on the temperature distribution, the total flow rate, and the inflow rates.
Yet another embodiment of the invention provides a system for determining production hydrocarbons. The system generally includes means for determining a temperature distribution associated with a plurality of fractured intervals or fractures disposed along a well; means for measuring a total flow rate for the well;
means for modeling an inflow rate for each of the plurality of fractured intervals or fractures; and means for allocating production of each of the plurality of fractured intervals or fractures based on the temperature distribution, the total flow rate, and the inflow rates.
BRIEF DESCRIPTION OF THE DRAWINGS
So that the manner in which the above recited features of the present invention can be understood in detail, a more particular description of the invention, briefly summarized above, may be had by reference to embodiments, some of which are illustrated in the appended drawings. It is to be noted, however, that the appended drawings illustrate only typical embodiments of this invention and are therefore not to be considered limiting of its scope, for the invention may admit to other equally effective embodiments.
FIG. 1 is a conceptual diagram of a system for producing hydrocarbons, the system having a pipe inside a casing and downhole tools positioned at various locations along the pipe, in accordance with an embodiment of the invention.
=
WELLBORE REAL-TIME MONITORING AND
ANALYSIS OF FRACTURE CONTRIBUTION
BACKGROUND OF THE INVENTION
Field of the Invention Embodiments of the present invention generally relate to hydrocarbon production and, more particularly, to determining the individual contribution of fractured intervals (or fractures) in time.
Description of the Related Art Various tools may be used in order to measure the contribution of the fractures within wellbores. Different services companies may run production logging tools, and chemical tracers may also be used to determine the fracture contribution.
However, these measurements may only provide a snapshot of what is happening at the moment the measurements are performed, and may change with time because conditions within the wellbore are transient.
SUMMARY OF THE INVENTION
Embodiments of the invention generally relate to allocating production of each of a plurality of fractured intervals (or fractures). This allocation may be performed by combining temperature distribution (and pressure) measurements, a real-time surface multiphase flow measurement, and an inflow model for each fractured interval (or fracture).
One embodiment of the invention is a method for determining production of hydrocarbons. The method generally includes determining a temperature distribution associated with a plurality of fractured intervals or fractures disposed along a well;
measuring a total flow rate for the well; modeling an inflow rate for each of the plurality of fractured intervals or fractures; and allocating production of each of the plurality of fractured intervals or fractures based on the temperature distribution, the total flow rate, and the inflow rates.
Another embodiment of the invention provides a system for determining production of hydrocarbons. The system generally includes a temperature sensing device configured to determine a temperature distribution associated with a plurality of fractured intervals or fractures disposed along a well, a flowmeter configured to measure a total flow rate for the well, and a processing unit. The processing unit is typically configured to model an inflow rate for each of the plurality of fractured intervals or fractures and to allocate production of each of the plurality of fractured intervals or fractures based on the temperature distribution, the total flow rate, and the inflow rates.
Yet another embodiment of the invention provides a system for determining production hydrocarbons. The system generally includes means for determining a temperature distribution associated with a plurality of fractured intervals or fractures disposed along a well; means for measuring a total flow rate for the well;
means for modeling an inflow rate for each of the plurality of fractured intervals or fractures; and means for allocating production of each of the plurality of fractured intervals or fractures based on the temperature distribution, the total flow rate, and the inflow rates.
BRIEF DESCRIPTION OF THE DRAWINGS
So that the manner in which the above recited features of the present invention can be understood in detail, a more particular description of the invention, briefly summarized above, may be had by reference to embodiments, some of which are illustrated in the appended drawings. It is to be noted, however, that the appended drawings illustrate only typical embodiments of this invention and are therefore not to be considered limiting of its scope, for the invention may admit to other equally effective embodiments.
FIG. 1 is a conceptual diagram of a system for producing hydrocarbons, the system having a pipe inside a casing and downhole tools positioned at various locations along the pipe, in accordance with an embodiment of the invention.
2 . .
FIG. 2 illustrates an ideal reservoir model with multiple fractures, in accordance with an embodiment of the invention.
FIG. 3 illustrates hydrocarbon production allocation from multiple wells, in accordance with an embodiment of the invention.
FIG. 4 illustrates hydrocarbon production allocation from a horizontal well with multiple fractured intervals, in accordance with an embodiment of the invention.
FIG. 5 is a flow diagram of example operations for allocating hydrocarbon production to multiple fractured intervals (or fractures), in accordance with an embodiment of the invention.
FIG. 6 illustrates a workflow for identifying and calculating the contribution of each fractured interval (or fracture), in accordance with an embodiment of the invention.
FIG. 7 illustrates an example plot of gas production versus number of contributing fractures, in accordance with an embodiment of the invention.
DETAILED DESCRIPTION
Embodiments of the invention provide techniques and apparatus for calculating production of each of a plurality of fractured intervals (or fractures) and monitoring changes in the fracture contribution with time. Such real-time monitoring and analysis may be based on a combination of different measurements in the wellbore, on the surface, and from a mathematical model, as described below.
In this manner, the industry may be able to understand the behavior of fractures and, in turn, optimize the number of stages (i.e., fractured intervals), the number of fractures, and the spacing between fractures and stages.
Referring to FIG. 1, there is shown a hydrocarbon production system 100 containing one or more production pipes 102 (also known as production tubing) that may extend downward through a casing 104 to one or more hydrocarbon sources 106 (e.g., reservoirs). An annulus 108 may exist between the pipe 102 and the
FIG. 2 illustrates an ideal reservoir model with multiple fractures, in accordance with an embodiment of the invention.
FIG. 3 illustrates hydrocarbon production allocation from multiple wells, in accordance with an embodiment of the invention.
FIG. 4 illustrates hydrocarbon production allocation from a horizontal well with multiple fractured intervals, in accordance with an embodiment of the invention.
FIG. 5 is a flow diagram of example operations for allocating hydrocarbon production to multiple fractured intervals (or fractures), in accordance with an embodiment of the invention.
FIG. 6 illustrates a workflow for identifying and calculating the contribution of each fractured interval (or fracture), in accordance with an embodiment of the invention.
FIG. 7 illustrates an example plot of gas production versus number of contributing fractures, in accordance with an embodiment of the invention.
DETAILED DESCRIPTION
Embodiments of the invention provide techniques and apparatus for calculating production of each of a plurality of fractured intervals (or fractures) and monitoring changes in the fracture contribution with time. Such real-time monitoring and analysis may be based on a combination of different measurements in the wellbore, on the surface, and from a mathematical model, as described below.
In this manner, the industry may be able to understand the behavior of fractures and, in turn, optimize the number of stages (i.e., fractured intervals), the number of fractures, and the spacing between fractures and stages.
Referring to FIG. 1, there is shown a hydrocarbon production system 100 containing one or more production pipes 102 (also known as production tubing) that may extend downward through a casing 104 to one or more hydrocarbon sources 106 (e.g., reservoirs). An annulus 108 may exist between the pipe 102 and the
3 . , casing 104. Each production pipe 102 may include one or more lateral sections (e.g., created by horizontal drilling) that branch off to access different hydrocarbon sources 106 or different areas of the same hydrocarbon source 106. The fluid mixture may flow from sources 106 to the well completion through the production pipes 102, as indicated by fluid flow 130. The production pipe 102 may include one or more tools 122 for performing various tasks (e.g., sensing parameters such as pressure or temperature) in, on, or adjacent a pipe or other conduit as the fluid mixtures flow through the production pipes 102. The tools 122 may be any type of downhole device, such as a flow control device (e.g., a valve), a sensor (e.g., a pressure, temperature or fluid flow sensor) or other instrument, an actuator (e.g., a solenoid), a data storage device (e.g., a programmable memory), a communication device (e.g., a transmitter or a receiver), etc.
Each tool 122 may be incorporated into an existing section of production pipe 102 or may be incorporated into a specific pipe section that is inserted in line with the production pipe 102. The distributed scheme of tools 122 shown in FIG. 1 may permit an operator of the system 100 to determine, for example, the level of depletion of the hydrocarbon reservoir. This information may permit the operator to monitor and intelligently control production of the hydrocarbon reservoir.
Advances in directional drilling (e.g., horizontal drilling as shown in FIG.
1) and reservoir stimulation techniques have dramatically increased gas production from wells drilled in shale reservoirs that were considered uneconomical not too long ago.
In spite of many advances in understanding the behavior of the production of this type of reservoir, many unknowns remain, such as determining the optimal length of horizontal sections, how many stages, and determining how many fractures are optimal. Particularly, it is difficult to predict productivity from cores, logs, drillstem tests (DSTs), or early well-production performance. Drainage volumes are uncertain, and well spacing is based on trial and error methods.
The use of microseismic and production logs has helped in the fracture evaluation to determine the drainage volume and fracture inflow. Microseismic can
Each tool 122 may be incorporated into an existing section of production pipe 102 or may be incorporated into a specific pipe section that is inserted in line with the production pipe 102. The distributed scheme of tools 122 shown in FIG. 1 may permit an operator of the system 100 to determine, for example, the level of depletion of the hydrocarbon reservoir. This information may permit the operator to monitor and intelligently control production of the hydrocarbon reservoir.
Advances in directional drilling (e.g., horizontal drilling as shown in FIG.
1) and reservoir stimulation techniques have dramatically increased gas production from wells drilled in shale reservoirs that were considered uneconomical not too long ago.
In spite of many advances in understanding the behavior of the production of this type of reservoir, many unknowns remain, such as determining the optimal length of horizontal sections, how many stages, and determining how many fractures are optimal. Particularly, it is difficult to predict productivity from cores, logs, drillstem tests (DSTs), or early well-production performance. Drainage volumes are uncertain, and well spacing is based on trial and error methods.
The use of microseismic and production logs has helped in the fracture evaluation to determine the drainage volume and fracture inflow. Microseismic can
4 . , provide useful information on the development of fracture symmetry, half-length, azimuth, width and height, and their dependence on the treatment parameters and reservoir characteristics. Additionally, these fracture geometries in conjunction with other measured or calculated parameters (e.g., rates, inflow models, etc.) can be used to better understand fracture modeling and production characteristics.
Review of production logs have indicated that only a percentage of the fractures are contributing to the production, and until now, only snapshots of the fracture contributions have been made. However, considering that this is a transient system (where fracture contributions typically change with time, typically for the first 15 to 20 months of production), a snapshot measurement is not sufficient to understand the behavior of the fractures and their contribution over time.
Accordingly, what is needed are techniques and apparatus for establishing which fractures (or at least which fractured intervals) are contributing and how much.
Due to the transient behavior, an ideal system would offer continuous, permanent, and real-time monitoring on key variables like production rates, pressure and temperature in an effort to determine the fracture contributions.
Procedures that integrate different types of measurements and calculations in "real time" may help to find and understand the behavior of the fractures and to optimize the number of stages, fractures, and spacing.
Embodiments of the invention provide methods and apparatus to optimize, or at least increase, the production of horizontal fractured wells in shale reservoirs, for example. By integrating different types of real-time measurements, methods described herein enable the optimization of the number of fractures, the spacing of fractures, and the length of the horizontal section by determining the contribution of the fracture stages (or the fractures) over time.
One way to solve this problem might be the installation of downhole flowmeters in each fracture stage. However, this can be a challenge operationally and may also be very costly and risky.
Review of production logs have indicated that only a percentage of the fractures are contributing to the production, and until now, only snapshots of the fracture contributions have been made. However, considering that this is a transient system (where fracture contributions typically change with time, typically for the first 15 to 20 months of production), a snapshot measurement is not sufficient to understand the behavior of the fractures and their contribution over time.
Accordingly, what is needed are techniques and apparatus for establishing which fractures (or at least which fractured intervals) are contributing and how much.
Due to the transient behavior, an ideal system would offer continuous, permanent, and real-time monitoring on key variables like production rates, pressure and temperature in an effort to determine the fracture contributions.
Procedures that integrate different types of measurements and calculations in "real time" may help to find and understand the behavior of the fractures and to optimize the number of stages, fractures, and spacing.
Embodiments of the invention provide methods and apparatus to optimize, or at least increase, the production of horizontal fractured wells in shale reservoirs, for example. By integrating different types of real-time measurements, methods described herein enable the optimization of the number of fractures, the spacing of fractures, and the length of the horizontal section by determining the contribution of the fracture stages (or the fractures) over time.
One way to solve this problem might be the installation of downhole flowmeters in each fracture stage. However, this can be a challenge operationally and may also be very costly and risky.
5 . =
Instead, considering the very low permeability of shale reservoirs (on the order of nanodarcys), it can be established that a reservoir is created only after fracturing.
If the spacing between fractures is correct (such that the fractures do not interfere with one another), the production allocation of each fracture stage (or fracture) may be calculated in an analogous way to that performed in a traditional field, where the total production rates are allocated to each production well using well testing measurements, done periodically with daily measurement information like wellhead pressure. In this particular case, by combining permanent downhole measurement of temperature (and one or more pressure measurements at the heel and the toe of the wellbore, for example), permanent wellhead flow measurement of the different phases, and a mathematical transient model of the production rates of each fracture, an acceptable production allocation can be made as a function of time. Because the system is transient, such allocation may be performed on a real-time basis.
In scenarios where the number of fractures is large, the idealized system 200 shown in FIG. 2 may be used to model the reservoir. In FIG. 2, multiple fractures 204, 206 are represented as spaced along and transverse to the horizontal well trajectory 202. Assuming fracturing conditions were the same, the length and width of each fracture in the fracture stage may be considered equal. These parallel fractures are formed in an area (e.g., a shale reservoir) with essentially zero permeability (as illustrated in the region 212 unshaded in FIG. 2), thereby forming a region 214 of modified permeability (shaded in FIG. 2), essentially creating a reservoir where none existed before. Although any number of fractures (Nfrac) may be formed with any spacing therebetween, five fractures are illustrated in the fracture stage of FIG. 2 (two external fractures 204 and three internal fractures 206) with equal fracture spacing. The fracture stage is defined by confining external boundaries 210. FIG. 2 shows that external fractures 204 are confined by virtual no-flow boundaries 208, which force the external fractures to have the same behavior as the internal fractures 206, and pure linear flow initially occurs. In shale gas reservoirs of nanodarcy permeability, pure linear flow opposite the fracture faces occurs for very long times.
Instead, considering the very low permeability of shale reservoirs (on the order of nanodarcys), it can be established that a reservoir is created only after fracturing.
If the spacing between fractures is correct (such that the fractures do not interfere with one another), the production allocation of each fracture stage (or fracture) may be calculated in an analogous way to that performed in a traditional field, where the total production rates are allocated to each production well using well testing measurements, done periodically with daily measurement information like wellhead pressure. In this particular case, by combining permanent downhole measurement of temperature (and one or more pressure measurements at the heel and the toe of the wellbore, for example), permanent wellhead flow measurement of the different phases, and a mathematical transient model of the production rates of each fracture, an acceptable production allocation can be made as a function of time. Because the system is transient, such allocation may be performed on a real-time basis.
In scenarios where the number of fractures is large, the idealized system 200 shown in FIG. 2 may be used to model the reservoir. In FIG. 2, multiple fractures 204, 206 are represented as spaced along and transverse to the horizontal well trajectory 202. Assuming fracturing conditions were the same, the length and width of each fracture in the fracture stage may be considered equal. These parallel fractures are formed in an area (e.g., a shale reservoir) with essentially zero permeability (as illustrated in the region 212 unshaded in FIG. 2), thereby forming a region 214 of modified permeability (shaded in FIG. 2), essentially creating a reservoir where none existed before. Although any number of fractures (Nfrac) may be formed with any spacing therebetween, five fractures are illustrated in the fracture stage of FIG. 2 (two external fractures 204 and three internal fractures 206) with equal fracture spacing. The fracture stage is defined by confining external boundaries 210. FIG. 2 shows that external fractures 204 are confined by virtual no-flow boundaries 208, which force the external fractures to have the same behavior as the internal fractures 206, and pure linear flow initially occurs. In shale gas reservoirs of nanodarcy permeability, pure linear flow opposite the fracture faces occurs for very long times.
6 =
The concept of Stimulated Reservoir Volume (SRV) is based on the premise that negligible flow occurs from beyond the fracture tips. The reservoir is created by the fracturing, and the reservoir size is limited by the length of the main fracture.
Production performance from the fractured reservoir may be based on the SRV, the fracture spacing, and the fracture conductivity.
The near-wellbore temperature distribution yielded by distributed temperature sensing (DTS) or multi-point or array temperature sensing (ATS) may be used to determine the relative amount of fluid that each perforation interval contributes. If this information is combined with one or more pressure measurements and a real-time surface multiphase flow measurement in conjunction with an inflow model for each fractured interval, a production allocation may be calculated for each fracture. This approach is analogous to a traditional well allocation where a daily aggregated measurement at the production plant is back-allocated to each well based on wellhead measurements like pressure, temperature, and well performance. The description below provides details on the use of these technologies to analyze the fracture behavior in horizontal wells in shale reservoirs, for example.
FIG. 3 illustrates a multi-well system 300 in an oil/gas production field, in which hydrocarbon production may be allocated to each of the wells. In this allocation process, periodical (e.g., 15 days to weeks or months) production well tests are performed on each individual well, and daily (or in some cases, every few hours) pressure (P) and/or temperature (T) measurements at or near the wellhead 302 of each well are registered. The produced fluids from each well may be collected at a manifold and then separated by a separator 310 into oil, gas, and water. Daily (or in some cases, every few hours or minutes) total flow rates of oil (Qo), gas (Qg), and water (Qw) may be measured. With the production well tests, using nodal analysis techniques, the well performance (P vs. Q relation) for each well at the wellhead 302 is calculated. The use of this wellhead performance with frequent wellhead pressure measurements allows the flow rates of each individual well to be determined.
The concept of Stimulated Reservoir Volume (SRV) is based on the premise that negligible flow occurs from beyond the fracture tips. The reservoir is created by the fracturing, and the reservoir size is limited by the length of the main fracture.
Production performance from the fractured reservoir may be based on the SRV, the fracture spacing, and the fracture conductivity.
The near-wellbore temperature distribution yielded by distributed temperature sensing (DTS) or multi-point or array temperature sensing (ATS) may be used to determine the relative amount of fluid that each perforation interval contributes. If this information is combined with one or more pressure measurements and a real-time surface multiphase flow measurement in conjunction with an inflow model for each fractured interval, a production allocation may be calculated for each fracture. This approach is analogous to a traditional well allocation where a daily aggregated measurement at the production plant is back-allocated to each well based on wellhead measurements like pressure, temperature, and well performance. The description below provides details on the use of these technologies to analyze the fracture behavior in horizontal wells in shale reservoirs, for example.
FIG. 3 illustrates a multi-well system 300 in an oil/gas production field, in which hydrocarbon production may be allocated to each of the wells. In this allocation process, periodical (e.g., 15 days to weeks or months) production well tests are performed on each individual well, and daily (or in some cases, every few hours) pressure (P) and/or temperature (T) measurements at or near the wellhead 302 of each well are registered. The produced fluids from each well may be collected at a manifold and then separated by a separator 310 into oil, gas, and water. Daily (or in some cases, every few hours or minutes) total flow rates of oil (Qo), gas (Qg), and water (Qw) may be measured. With the production well tests, using nodal analysis techniques, the well performance (P vs. Q relation) for each well at the wellhead 302 is calculated. The use of this wellhead performance with frequent wellhead pressure measurements allows the flow rates of each individual well to be determined.
7 . .
Ideally, the addition of all these individual well flow rates is the total production of the field, but for various reasons (e.g., well performance of each well can change over time), there is a difference between these values. To eliminate this difference, an allocation factor (K) is found using the relationship between the total flow rate (Qt) measured and the sum of the individual well flow rates (IQi) and may be subsequently used.
FIG. 4 illustrates a system 400 for allocating hydrocarbon produced from a horizontal well with multiple fractured intervals 402 along a horizontal well, in accordance with an embodiment of the invention. Although seven fractured intervals 402, each with five fractures 404, are shown in FIG. 4, any number of fractured intervals and any number of fractures per interval may be used. The system 400 also includes a multiphase real-time flowmeter 406 and a DTS cable 408 disposed downhole. The system may also include one or more sensors 410 for measuring pressure (P) and/or temperature (T), which may be disposed anywhere in the wellbore, such as in the vertical section as shown. The multiphase flowmeter may be installed at or adjacent the wellhead or within the wellbore and, for some embodiments, may be an optical flowmeter (e.g., an optical downhole flowmeter).
The DTS cable 408 may be installed adjacent the casing 104, as shown in FIG.
4.
Drawing an analogy to the multi-well system 300 of FIG. 3, each stage (i.e., fractured interval 402) in FIG. 4 is akin to a producing well. With the help of the variation of temperature and a transient inflow model, it is possible to calculate the production of each stage at any time. In fact, if the temperature variation is high enough to distinguish between fractures 404, it may also be possible to allocate the production of each particular fracture.
The analogy between production allocation for individual wells and stages (or fractures) is possible (i.e., each stage or fracture may be considered as an individual contributor to production) because, due to the low permeability of this type of reservoir (as described above with respect to FIG. 2), the communication between stages, and even between fractures, is negligible. The main characteristics of the
Ideally, the addition of all these individual well flow rates is the total production of the field, but for various reasons (e.g., well performance of each well can change over time), there is a difference between these values. To eliminate this difference, an allocation factor (K) is found using the relationship between the total flow rate (Qt) measured and the sum of the individual well flow rates (IQi) and may be subsequently used.
FIG. 4 illustrates a system 400 for allocating hydrocarbon produced from a horizontal well with multiple fractured intervals 402 along a horizontal well, in accordance with an embodiment of the invention. Although seven fractured intervals 402, each with five fractures 404, are shown in FIG. 4, any number of fractured intervals and any number of fractures per interval may be used. The system 400 also includes a multiphase real-time flowmeter 406 and a DTS cable 408 disposed downhole. The system may also include one or more sensors 410 for measuring pressure (P) and/or temperature (T), which may be disposed anywhere in the wellbore, such as in the vertical section as shown. The multiphase flowmeter may be installed at or adjacent the wellhead or within the wellbore and, for some embodiments, may be an optical flowmeter (e.g., an optical downhole flowmeter).
The DTS cable 408 may be installed adjacent the casing 104, as shown in FIG.
4.
Drawing an analogy to the multi-well system 300 of FIG. 3, each stage (i.e., fractured interval 402) in FIG. 4 is akin to a producing well. With the help of the variation of temperature and a transient inflow model, it is possible to calculate the production of each stage at any time. In fact, if the temperature variation is high enough to distinguish between fractures 404, it may also be possible to allocate the production of each particular fracture.
The analogy between production allocation for individual wells and stages (or fractures) is possible (i.e., each stage or fracture may be considered as an individual contributor to production) because, due to the low permeability of this type of reservoir (as described above with respect to FIG. 2), the communication between stages, and even between fractures, is negligible. The main characteristics of the
8 . .
fractures (e.g., length and width) may be considered equal in each stage, assuming fracturing conditions were the same. The inflow rate of each fracture will be computed by an analytical transient model and combined with the change in temperature (as determined by the DTS cable 408, for example) at each stage referenced to an initial condition prior to fracturing. In conjunction with the total flow rate (Qt) measured by the multiphase flowmeter 406, a production allocation for each stage (Qsi) (or each fracture) will be performed.
FIG. 5 is a flow diagram of example operations 500 for determining the contribution to hydrocarbon production of each fractured interval (or each fracture).
The operations 500 may begin, at 502, by determining a temperature distribution associated with a plurality of fractured intervals or fractures disposed along a well.
The temperature distribution may be determined by performing at least one of distributed temperature sensing (DTS) or array temperature sensing (ATS). The plurality of fractured intervals or fractures may be located in a shale reservoir, for example.
At 504, a total flow rate of a fluid (or any combination of fluids) produced by the well (i.e., the produced hydrocarbons) is measured. The total flow rate may be a total gas flow rate or a total oil flow rate, for example. For some embodiments, the total flow rate may be measured using a flowmeter disposed at the surface. For example, the flowmeter may be disposed at or adjacent a wellhead of the well.
An inflow rate is modeled at 506 for each of the plurality of fractured intervals or fractures. The inflow rate may be an inflow gas rate or an inflow oil rate, for example.
At 508, production of each of the plurality of fractured intervals or fractures is allocated based on the temperature distribution, the total flow rate, and the inflow rates. For some embodiments, allocating the production at 508 may include: (1) determining a first temperature value To at a first time to (e.g., before production starts) for each of the plurality of fractured intervals or fractures; (2) determining a second temperature value Tn at a second time tn (e.g., subsequent to the first time to)
fractures (e.g., length and width) may be considered equal in each stage, assuming fracturing conditions were the same. The inflow rate of each fracture will be computed by an analytical transient model and combined with the change in temperature (as determined by the DTS cable 408, for example) at each stage referenced to an initial condition prior to fracturing. In conjunction with the total flow rate (Qt) measured by the multiphase flowmeter 406, a production allocation for each stage (Qsi) (or each fracture) will be performed.
FIG. 5 is a flow diagram of example operations 500 for determining the contribution to hydrocarbon production of each fractured interval (or each fracture).
The operations 500 may begin, at 502, by determining a temperature distribution associated with a plurality of fractured intervals or fractures disposed along a well.
The temperature distribution may be determined by performing at least one of distributed temperature sensing (DTS) or array temperature sensing (ATS). The plurality of fractured intervals or fractures may be located in a shale reservoir, for example.
At 504, a total flow rate of a fluid (or any combination of fluids) produced by the well (i.e., the produced hydrocarbons) is measured. The total flow rate may be a total gas flow rate or a total oil flow rate, for example. For some embodiments, the total flow rate may be measured using a flowmeter disposed at the surface. For example, the flowmeter may be disposed at or adjacent a wellhead of the well.
An inflow rate is modeled at 506 for each of the plurality of fractured intervals or fractures. The inflow rate may be an inflow gas rate or an inflow oil rate, for example.
At 508, production of each of the plurality of fractured intervals or fractures is allocated based on the temperature distribution, the total flow rate, and the inflow rates. For some embodiments, allocating the production at 508 may include: (1) determining a first temperature value To at a first time to (e.g., before production starts) for each of the plurality of fractured intervals or fractures; (2) determining a second temperature value Tn at a second time tn (e.g., subsequent to the first time to)
9 . .
for each of the plurality of fractured intervals or fractures; (3) calculating a delta temperature value (ATn = Tn - To) for the second time tn for each of the plurality of fractured intervals or fractures by determining a difference between the first and second temperature values for each of the plurality of fractured intervals or fractures;
(4) calculating a first ratio (AT/Tg)n of the delta temperature value ATn for the second time tn for each of the plurality of fractured intervals or fractures to a geothermal temperature (Tg) at the second time tn (5) comparing the first ratio (AT/Tg),, for the second time tn to a maximum value of the first ratio over all previous times for each of the plurality of fractured intervals or fractures (6) for each of the plurality of fractured intervals or fractures, designating the first ratio for the second time tn as the maximum value of the first ratio over all previous times if the first ratio for the second time tn is greater than the previously designated maximum value (7) calculating a second ratio (AT/Tg)/(AT/Tg)max of the first ratio for the second time tõ for each of the plurality of fractured intervals or fractures to the currently designated maximum value of the first ratio over all previous times for each of the plurality of fractured intervals or fractures; (8) multiplying the second ratio for the second time tn with the modeled inflow rate corresponding to the second time tn for each of the plurality of fractured intervals or fractures; (9) summing results of the multiplication for each of the plurality of fractured intervals or fractures; (10) determining an allocation factor (K) by dividing the measured total flow rate corresponding to the second time tn by the sum; (11) applying the allocation factor (K) to the modeled inflow rate for each of the plurality of fractured intervals or fractures.
For some embodiments, the operations 500 may also include repeating the determining at 502, the measuring at 504, and the modeling at 506 within a period short enough to observe transient behavior of the plurality of fractured intervals or fractures. The determining, measuring, and/or modeling described above may be performed and repeated with any desired frequency (at any desired rate or periodicity). For example, the determining, measuring, and/or modeling may be performed continuously, hourly, daily, weekly, or with other frequencies.
For some embodiments, the operations 500 may also include determining one or more pressure measurements for the well. In this case, allocation of the production at 508 may also be based on the pressure measurements. The pressure measurements may be made by one or more pressure sensors located downhole, along the horizontal or vertical portion of the wellbore. The pressure sensors may be optical-based pressure sensors having one or more fiber Bragg gratings (FBGs) located therein.
FIG. 6 illustrates a workflow 600 for identifying and calculating the contribution of each fractured interval (or fracture), in accordance with an embodiment of the invention. For simplicity, the description below will focus on production allocation for each fractured interval. The workflow 600 can be easily expanded to production allocation for each fracture, as long as the temperature variation is high enough to distinguish between fractures.
In the workflow 600, the DTS (or ATS) data 602 is related to the geothermal gradient value for each stage 402. The cable 408 may be sampled with some periodicity to generate the data 602, leading to temperature measurements at certain sampling times (tn). For each sampling time (tn), the delta temperature (AT) between the temperature at the sampling time and at time to is calculated for each stage 402.
At 604, the AT values for each stage are divided by Tg to normalize the data.
For some embodiments, pressure measurements (e.g., taken by the sensors 410) may be used to ensure accuracy of the AT values for each stage (e.g., by correlation with the temperature measurements). At 606, a ratio ((ATTTg)/(AT/Tg)max) for the sampling time (tn) is calculated for each stage 402. The ratio for each stage is calculated by dividing the Tg-normalized AT value for this particular stage by the maximum Tg-normalized AT value over all previous times for this stage.
The AT value at time to is initially assumed to be the maximum Tg-normalized AT value, so the ratio in this case will be 1. The maximum AT value is stored for later validation of this assumption.
< .
At 608, inflow transient models are run to generate inflow rates for each stage 402 (indexed by "i"). The workflow 600 of FIG. 6 generates inflow gas rates for each stage (Qgfi), but inflow oil rates or both may also be used. The inflow transient models either produce the inflow rates at the sampling time (tn) as shown at 610, or interpolation or other techniques are used to determine inflow rates at the sampling time based on inflow rates produced for other times. At 612, the ratio at the sampling time (tn) for each stage calculated at 606 is multiplied with the modeled inflow rate for each stage from 610 corresponding to the sampling time.
As described above, surface multiphase measurements may be made at 614, for example, by the flowmeter 406, to generate one or more total flow rates (Qg, Qo, and/or Qw) for the well. The total flow rates may either be generated at the sampling time (tn) as shown at 616, or interpolation or other techniques may be used to determine the total flow rates at sampling time based on measurements taken at other times.
The results of the multiplications at 612 for each of the stages 402 at the sampling time (tr,) may be summed (ZQ'gfi). At 618, this sum may be compared to the total gas flow rate (Qg) corresponding to the sampling time (tn).
At the first sampling time (to), the ratio for each stage 402 calculated at 606 is multiplied by the Qgfi at to for each stage at 612, and the sum of all Qgfi values is compared to the Qg corresponding to to at 618. For this time to, it is being assumed that all fractures are contributing at their 100% capacities, unless the AT
value is zero, in the case of no contribution. For the next time tl, the value of AT, will be compared to the value of AT . If ATi is bigger, then a new maximum value is obtained. This new maximum value replaces the previous value, and in this case the contribution of this particular stage will be 100% during this period of time, and the assumption on the previous time step was wrong. A new calculation for to will be performed to correct the first assumption and similarly at any time that a new maximum value is found.
The workflow 600, operating on a "real-time" basis, will increase well productivity, helping to determine what is the optimal choke size to flow back the well and to have all fractures contributing (or to find out which fractures do not contribute at all). After this procedure is performed on different wells with a different number of stages and/or fractures, a normalized graph of production versus a number of contributing stages and/or fractures can be obtained and, based on these results, an optimal number of stages and/or fractures may be determined. A good relationship is expected of production versus number of contributing fractures, more consistent than the plot 700 of gas production versus number of contributing fractures shown in FIG.
7 (from Modeland N. etal., "Stimulation's Influence on Production in the Haynesville Shale: A Playwide Examination of Fracture-Treatment Variables that Show Effect on Production," SPE 148940 presented at Canadian Unconventional Resources Conference, 15-17 November 2011, Alberta, Canada).
As described above, the near-wellbore temperature distribution yielded by distributed temperature sensing (DTS) or multi-point or array temperature sensing (ATS) may be used to determine the relative amount of fluid that each perforation interval contributes. If this information is combined with a real-time surface multiphase flow measurement in conjunction with an inflow model for each fractured interval (and one or more pressure measurements), a production allocation may be calculated for each fractured interval or fracture. This approach is analogous to a traditional well allocation where a daily aggregated measurement at the production plant is back-allocated to each well based on wellhead measurements like pressure, temperature, and well performance.
While the foregoing is directed to embodiments of the present invention, other and further embodiments of the invention may be devised without departing from the basic scope thereof, and the scope thereof is determined by the claims that follow.
for each of the plurality of fractured intervals or fractures; (3) calculating a delta temperature value (ATn = Tn - To) for the second time tn for each of the plurality of fractured intervals or fractures by determining a difference between the first and second temperature values for each of the plurality of fractured intervals or fractures;
(4) calculating a first ratio (AT/Tg)n of the delta temperature value ATn for the second time tn for each of the plurality of fractured intervals or fractures to a geothermal temperature (Tg) at the second time tn (5) comparing the first ratio (AT/Tg),, for the second time tn to a maximum value of the first ratio over all previous times for each of the plurality of fractured intervals or fractures (6) for each of the plurality of fractured intervals or fractures, designating the first ratio for the second time tn as the maximum value of the first ratio over all previous times if the first ratio for the second time tn is greater than the previously designated maximum value (7) calculating a second ratio (AT/Tg)/(AT/Tg)max of the first ratio for the second time tõ for each of the plurality of fractured intervals or fractures to the currently designated maximum value of the first ratio over all previous times for each of the plurality of fractured intervals or fractures; (8) multiplying the second ratio for the second time tn with the modeled inflow rate corresponding to the second time tn for each of the plurality of fractured intervals or fractures; (9) summing results of the multiplication for each of the plurality of fractured intervals or fractures; (10) determining an allocation factor (K) by dividing the measured total flow rate corresponding to the second time tn by the sum; (11) applying the allocation factor (K) to the modeled inflow rate for each of the plurality of fractured intervals or fractures.
For some embodiments, the operations 500 may also include repeating the determining at 502, the measuring at 504, and the modeling at 506 within a period short enough to observe transient behavior of the plurality of fractured intervals or fractures. The determining, measuring, and/or modeling described above may be performed and repeated with any desired frequency (at any desired rate or periodicity). For example, the determining, measuring, and/or modeling may be performed continuously, hourly, daily, weekly, or with other frequencies.
For some embodiments, the operations 500 may also include determining one or more pressure measurements for the well. In this case, allocation of the production at 508 may also be based on the pressure measurements. The pressure measurements may be made by one or more pressure sensors located downhole, along the horizontal or vertical portion of the wellbore. The pressure sensors may be optical-based pressure sensors having one or more fiber Bragg gratings (FBGs) located therein.
FIG. 6 illustrates a workflow 600 for identifying and calculating the contribution of each fractured interval (or fracture), in accordance with an embodiment of the invention. For simplicity, the description below will focus on production allocation for each fractured interval. The workflow 600 can be easily expanded to production allocation for each fracture, as long as the temperature variation is high enough to distinguish between fractures.
In the workflow 600, the DTS (or ATS) data 602 is related to the geothermal gradient value for each stage 402. The cable 408 may be sampled with some periodicity to generate the data 602, leading to temperature measurements at certain sampling times (tn). For each sampling time (tn), the delta temperature (AT) between the temperature at the sampling time and at time to is calculated for each stage 402.
At 604, the AT values for each stage are divided by Tg to normalize the data.
For some embodiments, pressure measurements (e.g., taken by the sensors 410) may be used to ensure accuracy of the AT values for each stage (e.g., by correlation with the temperature measurements). At 606, a ratio ((ATTTg)/(AT/Tg)max) for the sampling time (tn) is calculated for each stage 402. The ratio for each stage is calculated by dividing the Tg-normalized AT value for this particular stage by the maximum Tg-normalized AT value over all previous times for this stage.
The AT value at time to is initially assumed to be the maximum Tg-normalized AT value, so the ratio in this case will be 1. The maximum AT value is stored for later validation of this assumption.
< .
At 608, inflow transient models are run to generate inflow rates for each stage 402 (indexed by "i"). The workflow 600 of FIG. 6 generates inflow gas rates for each stage (Qgfi), but inflow oil rates or both may also be used. The inflow transient models either produce the inflow rates at the sampling time (tn) as shown at 610, or interpolation or other techniques are used to determine inflow rates at the sampling time based on inflow rates produced for other times. At 612, the ratio at the sampling time (tn) for each stage calculated at 606 is multiplied with the modeled inflow rate for each stage from 610 corresponding to the sampling time.
As described above, surface multiphase measurements may be made at 614, for example, by the flowmeter 406, to generate one or more total flow rates (Qg, Qo, and/or Qw) for the well. The total flow rates may either be generated at the sampling time (tn) as shown at 616, or interpolation or other techniques may be used to determine the total flow rates at sampling time based on measurements taken at other times.
The results of the multiplications at 612 for each of the stages 402 at the sampling time (tr,) may be summed (ZQ'gfi). At 618, this sum may be compared to the total gas flow rate (Qg) corresponding to the sampling time (tn).
At the first sampling time (to), the ratio for each stage 402 calculated at 606 is multiplied by the Qgfi at to for each stage at 612, and the sum of all Qgfi values is compared to the Qg corresponding to to at 618. For this time to, it is being assumed that all fractures are contributing at their 100% capacities, unless the AT
value is zero, in the case of no contribution. For the next time tl, the value of AT, will be compared to the value of AT . If ATi is bigger, then a new maximum value is obtained. This new maximum value replaces the previous value, and in this case the contribution of this particular stage will be 100% during this period of time, and the assumption on the previous time step was wrong. A new calculation for to will be performed to correct the first assumption and similarly at any time that a new maximum value is found.
The workflow 600, operating on a "real-time" basis, will increase well productivity, helping to determine what is the optimal choke size to flow back the well and to have all fractures contributing (or to find out which fractures do not contribute at all). After this procedure is performed on different wells with a different number of stages and/or fractures, a normalized graph of production versus a number of contributing stages and/or fractures can be obtained and, based on these results, an optimal number of stages and/or fractures may be determined. A good relationship is expected of production versus number of contributing fractures, more consistent than the plot 700 of gas production versus number of contributing fractures shown in FIG.
7 (from Modeland N. etal., "Stimulation's Influence on Production in the Haynesville Shale: A Playwide Examination of Fracture-Treatment Variables that Show Effect on Production," SPE 148940 presented at Canadian Unconventional Resources Conference, 15-17 November 2011, Alberta, Canada).
As described above, the near-wellbore temperature distribution yielded by distributed temperature sensing (DTS) or multi-point or array temperature sensing (ATS) may be used to determine the relative amount of fluid that each perforation interval contributes. If this information is combined with a real-time surface multiphase flow measurement in conjunction with an inflow model for each fractured interval (and one or more pressure measurements), a production allocation may be calculated for each fractured interval or fracture. This approach is analogous to a traditional well allocation where a daily aggregated measurement at the production plant is back-allocated to each well based on wellhead measurements like pressure, temperature, and well performance.
While the foregoing is directed to embodiments of the present invention, other and further embodiments of the invention may be devised without departing from the basic scope thereof, and the scope thereof is determined by the claims that follow.
Claims (20)
1. A method for determining production of hydrocarbons, comprising:
determining a temperature distribution associated with a plurality of fractured intervals or fractures disposed along a well;
measuring a total flow rate for the well;
modeling an inflow rate for each of the plurality of fractured intervals or fractures; and allocating production of each of the plurality of fractured intervals or fractures based on the temperature distribution, the total flow rate, and the inflow rates.
determining a temperature distribution associated with a plurality of fractured intervals or fractures disposed along a well;
measuring a total flow rate for the well;
modeling an inflow rate for each of the plurality of fractured intervals or fractures; and allocating production of each of the plurality of fractured intervals or fractures based on the temperature distribution, the total flow rate, and the inflow rates.
2. The method of claim 1, further comprising repeating the determining, the measuring, and the modeling within a period short enough to observe transient behavior of the plurality of fractured intervals or fractures.
3. The method of claim 1, further comprising determining one or more pressure measurements for the well, wherein allocating the production is further based on the pressure measurements.
4. The method of claim 1, wherein determining the temperature distribution comprises performing at least one of distributed temperature sensing (DTS) or array temperature sensing (ATS).
5. The method of claim 1, wherein the measuring comprises measuring the total flow rate using a multiphase flowmeter.
6. The method of claim 1, wherein at least one of the determining, the measuring, or the modeling is performed daily.
7. The method of claim 1, wherein at least one of the determining, the measuring, or the modeling is performed continuously.
8. The method of claim 1, wherein allocating the production comprises:
determining a first temperature value at a first time for each of the plurality of fractured intervals or fractures;
determining a second temperature value at a second time for each of the plurality of fractured intervals or fractures;
calculating a delta temperature value for the second time for each of the plurality of fractured intervals or fractures by determining a difference between the first and second temperature values for each of the plurality of fractured intervals or fractures;
calculating a first ratio of the delta temperature value for the second time for each of the plurality of fractured intervals or fractures to a geothermal temperature;
comparing the first ratio for the second time to a maximum value of the first ratio over all previous times for each of the plurality of fractured intervals or fractures;
for each of the plurality of fractured intervals or fractures, designating the first ratio for the second time as the maximum value of the first ratio over all previous times if the first ratio for the second time is greater than a previously designated maximum value;
for each of the plurality of fractured intervals or fractures, calculating a second ratio of the first ratio for the second time to a currently designated maximum value of the first ratio over all previous times;
multiplying the second ratio for the second time with the modeled inflow rate corresponding to the second time for each of the plurality of fractured intervals or fractures;
summing results of the multiplication for each of the plurality of fractured intervals or fractures; and determining an allocation factor by dividing the measured total flow rate corresponding to the second time by the sum.
determining a first temperature value at a first time for each of the plurality of fractured intervals or fractures;
determining a second temperature value at a second time for each of the plurality of fractured intervals or fractures;
calculating a delta temperature value for the second time for each of the plurality of fractured intervals or fractures by determining a difference between the first and second temperature values for each of the plurality of fractured intervals or fractures;
calculating a first ratio of the delta temperature value for the second time for each of the plurality of fractured intervals or fractures to a geothermal temperature;
comparing the first ratio for the second time to a maximum value of the first ratio over all previous times for each of the plurality of fractured intervals or fractures;
for each of the plurality of fractured intervals or fractures, designating the first ratio for the second time as the maximum value of the first ratio over all previous times if the first ratio for the second time is greater than a previously designated maximum value;
for each of the plurality of fractured intervals or fractures, calculating a second ratio of the first ratio for the second time to a currently designated maximum value of the first ratio over all previous times;
multiplying the second ratio for the second time with the modeled inflow rate corresponding to the second time for each of the plurality of fractured intervals or fractures;
summing results of the multiplication for each of the plurality of fractured intervals or fractures; and determining an allocation factor by dividing the measured total flow rate corresponding to the second time by the sum.
9. The method of claim 8, wherein the first time occurs before the hydrocarbons are produced.
10. The method of claim 8, further comprising applying the allocation factor to the modeled inflow rate for each of the plurality of fractured intervals or fractures.
11. The method of claim 1, wherein the total flow rate comprises a total gas flow rate and wherein the inflow rates comprise inflow gas rates.
12. The method of claim 1, wherein the plurality of fractured intervals or fractures is located in a shale reservoir.
13. A system for determining production of hydrocarbons, comprising:
a temperature sensing device configured to determine a temperature distribution associated with a plurality of fractured intervals or fractures disposed along a well;
a flowmeter configured to measure a total flow rate for the well; and a processing unit configured to:
model an inflow rate for each of the plurality of fractured intervals or fractures; and allocate production of each of the plurality of fractured intervals or fractures based on the temperature distribution, the total flow rate, and the inflow rates.
a temperature sensing device configured to determine a temperature distribution associated with a plurality of fractured intervals or fractures disposed along a well;
a flowmeter configured to measure a total flow rate for the well; and a processing unit configured to:
model an inflow rate for each of the plurality of fractured intervals or fractures; and allocate production of each of the plurality of fractured intervals or fractures based on the temperature distribution, the total flow rate, and the inflow rates.
14. The system of claim 13, wherein the plurality of fractured intervals or fractures is located in a shale reservoir.
15. The system of claim 13, further comprising a pressure sensor configured to determine one or more pressure measurements for the well, wherein the processing unit is configured to allocate the production further based on the pressure measurements.
16. The system of claim 13, wherein the processing unit is configured to allocate the production by:
determining a first temperature value at a first time for each of the plurality of fractured intervals or fractures;
determining a second temperature value at a second time for each of the plurality of fractured intervals or fractures;
calculating a delta temperature value for the second time for each of the plurality of fractured intervals or fractures by determining a difference between the first and second temperature values for each of the plurality of fractured intervals or fractures;
calculating a first ratio of the delta temperature value for the second time for each of the plurality of fractured intervals or fractures to a geothermal temperature;
comparing the first ratio for the second time to a maximum value of the first ratio over all previous times for each of the plurality of fractured intervals or fractures;
for each of the plurality of fractured intervals or fractures, designating the first ratio for the second time as the maximum value of the first ratio over all previous times if the first ratio for the second time is greater than a previously designated maximum value;
for each of the plurality of fractured intervals or fractures, calculating a second ratio of the first ratio for the second time to a currently designated maximum value of the first ratio over all previous times;
multiplying the second ratio for the second time with the modeled inflow rate corresponding to the second time for each of the plurality of fractured intervals or fractures;
summing results of the multiplication for each of the plurality of fractured intervals or fractures; and determining an allocation factor by dividing the measured total flow rate corresponding to the second time by the sum.
determining a first temperature value at a first time for each of the plurality of fractured intervals or fractures;
determining a second temperature value at a second time for each of the plurality of fractured intervals or fractures;
calculating a delta temperature value for the second time for each of the plurality of fractured intervals or fractures by determining a difference between the first and second temperature values for each of the plurality of fractured intervals or fractures;
calculating a first ratio of the delta temperature value for the second time for each of the plurality of fractured intervals or fractures to a geothermal temperature;
comparing the first ratio for the second time to a maximum value of the first ratio over all previous times for each of the plurality of fractured intervals or fractures;
for each of the plurality of fractured intervals or fractures, designating the first ratio for the second time as the maximum value of the first ratio over all previous times if the first ratio for the second time is greater than a previously designated maximum value;
for each of the plurality of fractured intervals or fractures, calculating a second ratio of the first ratio for the second time to a currently designated maximum value of the first ratio over all previous times;
multiplying the second ratio for the second time with the modeled inflow rate corresponding to the second time for each of the plurality of fractured intervals or fractures;
summing results of the multiplication for each of the plurality of fractured intervals or fractures; and determining an allocation factor by dividing the measured total flow rate corresponding to the second time by the sum.
17. The system of claim 16, wherein the processing unit is further configured to apply the allocation factor to the modeled inflow rate for each of the plurality of fractured intervals or fractures.
18. The system of claim 13, wherein the temperature sensing device comprises a distributed temperature sensing (DTS) device or an array temperature sensing (ATS) device.
19. The system of claim 13, wherein the total flow rate comprises a total gas flow rate and wherein the inflow rates comprise inflow gas rates.
20. A system for determining production of hydrocarbons, comprising:
means for determining a temperature distribution associated with a plurality of fractured intervals or fractures disposed along a well;
means for measuring a total flow rate for the well;
means for modeling an inflow rate for each of the plurality of fractured intervals or fractures; and means for allocating production of each of the plurality of fractured intervals or fractures based on the temperature distribution, the total flow rate, and the inflow rates.
means for determining a temperature distribution associated with a plurality of fractured intervals or fractures disposed along a well;
means for measuring a total flow rate for the well;
means for modeling an inflow rate for each of the plurality of fractured intervals or fractures; and means for allocating production of each of the plurality of fractured intervals or fractures based on the temperature distribution, the total flow rate, and the inflow rates.
Applications Claiming Priority (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| US201261611924P | 2012-03-16 | 2012-03-16 | |
| US61/611,924 | 2012-03-16 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CA2808858A1 CA2808858A1 (en) | 2013-09-16 |
| CA2808858C true CA2808858C (en) | 2016-01-26 |
Family
ID=48049763
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CA2808858A Active CA2808858C (en) | 2012-03-16 | 2013-03-11 | Wellbore real-time monitoring and analysis of fracture contribution |
Country Status (7)
| Country | Link |
|---|---|
| US (1) | US20130245953A1 (en) |
| EP (1) | EP2639401B1 (en) |
| CN (1) | CN103306664A (en) |
| AR (1) | AR090353A1 (en) |
| AU (1) | AU2013201757B2 (en) |
| BR (1) | BR102013006266B1 (en) |
| CA (1) | CA2808858C (en) |
Families Citing this family (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US10400580B2 (en) | 2015-07-07 | 2019-09-03 | Schlumberger Technology Corporation | Temperature sensor technique for determining a well fluid characteristic |
| US10415382B2 (en) * | 2016-05-03 | 2019-09-17 | Schlumberger Technology Corporation | Method and system for establishing well performance during plug mill-out or cleanout/workover operations |
| US10303819B2 (en) | 2016-08-25 | 2019-05-28 | Drilling Info, Inc. | Systems and methods for allocating hydrocarbon production values |
| US11263370B2 (en) | 2016-08-25 | 2022-03-01 | Enverus, Inc. | Systems and methods for allocating hydrocarbon production values |
| US11892579B2 (en) * | 2016-09-30 | 2024-02-06 | Schlumberger Technology Corporation | Crosswell microseismic system |
| US10584577B2 (en) | 2018-03-13 | 2020-03-10 | Saudi Arabian Oil Company | In-situ reservoir depletion management based on surface characteristics of production |
| CN110318742B (en) * | 2018-03-30 | 2022-07-15 | 中国石油化工股份有限公司 | Method and system for determining fracture closure length based on fractured well production data |
| EP3837429A4 (en) * | 2018-08-16 | 2022-08-24 | Fervo Energy Company | Methods and systems to control flow and heat transfer between subsurface wellbores connected hydraulically by fractures |
| US11293280B2 (en) * | 2018-12-19 | 2022-04-05 | Exxonmobil Upstream Research Company | Method and system for monitoring post-stimulation operations through acoustic wireless sensor network |
| US11326440B2 (en) | 2019-09-18 | 2022-05-10 | Exxonmobil Upstream Research Company | Instrumented couplings |
| CN111255442B (en) * | 2020-01-14 | 2023-04-07 | 大庆油田有限责任公司 | Method for evaluating fracturing fracture by using interference well testing theory |
| CN112878982B (en) * | 2020-12-31 | 2022-03-01 | 西南石油大学 | Deep shale gas productivity prediction method considering long-term fracture conductivity |
| CN113187472B (en) * | 2021-05-11 | 2023-09-26 | 中国石油天然气股份有限公司 | Identification method of water-flooding development seepage dominant channel of layered sandstone reservoir |
Family Cites Families (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US3913398A (en) * | 1973-10-09 | 1975-10-21 | Schlumberger Technology Corp | Apparatus and method for determining fluid flow rates from temperature log data |
| FR2538849A1 (en) * | 1982-12-30 | 1984-07-06 | Schlumberger Prospection | METHOD AND DEVICE FOR DETERMINING THE FLOW PROPERTIES OF A FLUID IN A WELL FROM TEMPERATURE MEASUREMENTS |
| US6142229A (en) * | 1998-09-16 | 2000-11-07 | Atlantic Richfield Company | Method and system for producing fluids from low permeability formations |
| DZ3413A1 (en) * | 2000-09-12 | 2002-03-21 | Sofitech Nv | EVALUATION OF MULTI-LAYERED AMALGAMATED TANK AND HYDRAULIC FRACTURE PROPERTIES USING AMALGAMATED TANK PRODUCTION DATA AND PRODUCTION LOGGING INFORMATION |
| US6789937B2 (en) * | 2001-11-30 | 2004-09-14 | Schlumberger Technology Corporation | Method of predicting formation temperature |
| US7703525B2 (en) * | 2004-12-03 | 2010-04-27 | Halliburton Energy Services, Inc. | Well perforating and fracturing |
| US20110011595A1 (en) * | 2008-05-13 | 2011-01-20 | Hao Huang | Modeling of Hydrocarbon Reservoirs Using Design of Experiments Methods |
| CN201334901Y (en) * | 2008-07-01 | 2009-10-28 | 电子科大科园股份有限公司 | Safety real time monitoring system for gas drilling |
| EP2334904A1 (en) * | 2008-08-08 | 2011-06-22 | Altarock Energy, Inc. | Method for testing an engineered geothermal system using one stimulated well |
| US8463585B2 (en) * | 2009-05-22 | 2013-06-11 | Baker Hughes Incorporated | Apparatus and method for modeling well designs and well performance |
| US8788251B2 (en) * | 2010-05-21 | 2014-07-22 | Schlumberger Technology Corporation | Method for interpretation of distributed temperature sensors during wellbore treatment |
-
2013
- 2013-03-11 CA CA2808858A patent/CA2808858C/en active Active
- 2013-03-14 US US13/828,055 patent/US20130245953A1/en not_active Abandoned
- 2013-03-15 BR BR102013006266-9A patent/BR102013006266B1/en active IP Right Grant
- 2013-03-15 EP EP13159586.0A patent/EP2639401B1/en active Active
- 2013-03-15 AR ARP130100848A patent/AR090353A1/en active IP Right Grant
- 2013-03-15 AU AU2013201757A patent/AU2013201757B2/en active Active
- 2013-03-18 CN CN2013100850393A patent/CN103306664A/en active Pending
Also Published As
| Publication number | Publication date |
|---|---|
| AU2013201757B2 (en) | 2015-10-22 |
| BR102013006266A8 (en) | 2017-07-11 |
| CA2808858A1 (en) | 2013-09-16 |
| AR090353A1 (en) | 2014-11-05 |
| BR102013006266A2 (en) | 2015-07-07 |
| EP2639401A1 (en) | 2013-09-18 |
| US20130245953A1 (en) | 2013-09-19 |
| AU2013201757A1 (en) | 2013-10-03 |
| EP2639401B1 (en) | 2016-05-04 |
| CN103306664A (en) | 2013-09-18 |
| BR102013006266B1 (en) | 2021-02-17 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CA2808858C (en) | Wellbore real-time monitoring and analysis of fracture contribution | |
| US9341060B2 (en) | Method and system for permeability calculation using production logs for horizontal wells | |
| CA2922573C (en) | Controlling an injection treatment of a subterranean region based on stride test data | |
| EP2038809B1 (en) | Method for comparing and back allocating production | |
| US9574443B2 (en) | Designing an injection treatment for a subterranean region based on stride test data | |
| US9500076B2 (en) | Injection testing a subterranean region | |
| US9341557B2 (en) | Method and system for permeability calculation using production logs for horizontal wells, using a downhole tool | |
| US20120158310A1 (en) | Method of determining reservoir pressure | |
| EA015435B1 (en) | A method of modeling well technological indices | |
| US11556612B2 (en) | Predicting material distribution in a hydraulic fracturing treatment stage | |
| Abbasi | A comparative study of flowback rate and pressure transient behaviour in multifractured horizontal wells | |
| CN110945209A (en) | Improvements in or relating to injection wells | |
| Williams-Kovacs et al. | Analysis of multi-well and stage-by-stage flowback from multi-fractured horizontal wells | |
| Camilleri et al. | Delivering pressure transient analysis during drawdown on ESP wells: case studies and lessons learned | |
| Shurunov et al. | Application of the HW with MSHF investigations to manage the development of low-permeability reservoirs | |
| Gonzalez et al. | Wellbore real-time monitoring and analysis for shale reservoirs | |
| Zhan et al. | Using an innovative tool system to estimate in-situ permeability and pressure at multiple targets in a monitoring well in Permian basin | |
| US11885220B2 (en) | System to determine existing fluids remaining saturation in homogenous and/or naturally fractured reservoirs | |
| Wang | Processing and analysis of transient pressure from permanent down-hole gauges | |
| Tandon | Identification of productive zones in unconventional reservoirs | |
| Sun | Implementation and application of fracture diagnostic tools: fiber optic sensing and diagnostic fracture injection test (DFIT) | |
| Zeinabadybejestani | Advancing Design and Analysis of the Diagnostic Fracture Injection Test-Flowback Analysis ('DFIT-FBA') Method and Post-Fracture Pressure Decay (PFPD) Technique | |
| Brown | Investigating The Impact Of Offset Fracture Hits Using Rate Transient Analysis In The Bakken And Three Forks Formation, Divide County, North Dakota | |
| Denney | Effective Completion and Reservoir Management With Well Logs | |
| RODRIGUEZ | RESERVOIR CHARACTERIZATION USING RATE-TRANSIENT ANALYSIS IN THE EAGLE FORD SHALE |