You must log in to edit PetroWiki. Help with editing
Content of PetroWiki is intended for personal use only and to supplement, not replace, engineering judgment. SPE disclaims any and all liability for your use of such content. More information
Interpreting data from well to well tracer tests: Difference between revisions
No edit summary |
No edit summary |
||
Line 1: | Line 1: | ||
To obtain high-quality tracer-response curves that are the basis for the further interpretation, a well-designed sampling program is needed. In general, more samples will give the potential for extraction of more information from field tests. Too often, interpretation is difficult because of limited tracer data. | To obtain high-quality tracer-response curves that are the basis for the further interpretation, a well-designed sampling program is needed. In general, more samples will give the potential for extraction of more information from field tests. Too often, interpretation is difficult because of limited tracer data. | ||
The final objective of a [[ | The final objective of a [[Well_to_well_tracer_tests|well-to-well study]] is the interpretation of the response curves. A good analysis of the information given by the tracers, in combination with other available data, gives a better understanding of the flow in the reservoir, not just verification of communication between injector and producer. | ||
==Different levels of interpretation== | == Different levels of interpretation == | ||
The | The response curve can be analyzed from three different points of view or complexity levels. The simplest interpretation is the qualitative one. By observing the curves, the following pattern characteristics can be derived: injection-water arriving time (breakthrough); existence of high-permeability channels, barriers, and fractures between wells; communication between different layers; stratification in the same layer; and preferential flow directions in the reservoir. Furthermore, the cumulative response can be obtained by integrating the concentration vs. time curve if the production flow rate is known. From this new curve, the fraction of injection water reaching each producer is easily calculated. A common spreadsheet is the best way to make these calculations. This type of interpretation can be carried out without any advanced simulation. It is important to integrate the data obtained from the geological model, primary production data, 4D seismic data (if available), and tracer data in a systematic way. | ||
Finally, complex mathematical models like numerical finite-element simulators or streamtube modeling can be used to achieve a deeper analysis. Most of the commercially available simulators have tracer options with varying degrees of complexity. Not all simulators include important physical effects like partitioning, dispersion, and adsorption. A tracer simulator that can be coupled to full-field reservoir simulators has been developed recently. Simulation is a tool to improve the existing reservoir model; therefore, it is crucial to have access to the best available model to enable an efficient optimization based on the tracer-production data. | The second level of analysis uses basic mathematical models to fit simple response curves by means of theoretical expressions and to decompose complex responses in several simpler functions. In this way, partial residence times, as well as other parameters, can be determined for each function. Mathematical models also allow for the evaluation of important parameters, such as permeability, and make it possible to predict the behavior of unknown patterns. | ||
Finally, complex mathematical models like numerical finite-element simulators or streamtube modeling can be used to achieve a deeper analysis. Most of the commercially available simulators have tracer options with varying degrees of complexity. Not all simulators include important physical effects like partitioning, dispersion, and adsorption. A tracer simulator that can be coupled to full-field reservoir simulators has been developed recently. Simulation is a tool to improve the existing reservoir model; therefore, it is crucial to have access to the best available model to enable an efficient optimization based on the tracer-production data. | |||
== Evaluation based on hand calculation and adequate data presentation == | == Evaluation based on hand calculation and adequate data presentation == | ||
Qualitative interpretation of field data is illustrated by examples from the Snorre field in the North Sea. '''Fig. 1''' shows the segment of the field that was tested and the location of injectors (triangles) and producers (circles). <ref name="r1" /> A water-alternating-gas (WAG) injection program is being used in the field. | Qualitative interpretation of field data is illustrated by examples from the Snorre field in the North Sea. '''Fig. 1''' shows the segment of the field that was tested and the location of injectors (triangles) and producers (circles). <ref name="r1">_</ref> A water-alternating-gas (WAG) injection program is being used in the field. | ||
<gallery widths="300px" heights="200px"> | <gallery widths="300px" heights="200px"> | ||
Line 18: | Line 19: | ||
</gallery> | </gallery> | ||
'''Fig. 2''' shows a way to present field data for a gas-tracer program. The data are plotted as a time response. Alternatively, the data can be plotted as a volumetric response. To obtain a better understanding, it is important to plot gas-tracer-production data with GOR. The tracer data are plotted relative to the amount of tracer injected, which means that the concentrations presented for each individual tracer can be compared directly. | '''Fig. 2''' shows a way to present field data for a gas-tracer program. The data are plotted as a time response. Alternatively, the data can be plotted as a volumetric response. To obtain a better understanding, it is important to plot gas-tracer-production data with GOR. The tracer data are plotted relative to the amount of tracer injected, which means that the concentrations presented for each individual tracer can be compared directly. | ||
<gallery widths="300px" heights="200px"> | <gallery widths="300px" heights="200px"> | ||
Line 24: | Line 25: | ||
</gallery> | </gallery> | ||
From '''Fig. 2''', it is possible not only to find communication but also to extract information about the prevailing flow regime. The gas tracer injected into Well P25 breaks through without any significant increase in GOR, which means that the gas has not followed a free gas layer but has been dissolved in the oil. Three different gas tracers were injected in Well P28 at three different times. These tracers were produced only after a significant GOR was recorded. All the tracers were produced at almost the same time. This shows that gas has been trapped and that the gas in the individual gas slugs has subsequently filled this trap. At a certain time, the gas trap has been filled and the gas moves on to the producers. When reaching the producer, the gas contains all three different tracers injected over a 3-year period. | From '''Fig. 2''', it is possible not only to find communication but also to extract information about the prevailing flow regime. The gas tracer injected into Well P25 breaks through without any significant increase in GOR, which means that the gas has not followed a free gas layer but has been dissolved in the oil. Three different gas tracers were injected in Well P28 at three different times. These tracers were produced only after a significant GOR was recorded. All the tracers were produced at almost the same time. This shows that gas has been trapped and that the gas in the individual gas slugs has subsequently filled this trap. At a certain time, the gas trap has been filled and the gas moves on to the producers. When reaching the producer, the gas contains all three different tracers injected over a 3-year period. | ||
From the data in '''Fig. 2''', it is possible to carry out a spreadsheet calculation to estimate the recovery of the injected tracer in each individual production well. This shows that a large percentage of the injected gas has moved in different directions. Breakthrough time also may be estimated, either as the first sample that contains tracer above the detection limit or as an extrapolation backward of the initial phase of the production curve. However, the accuracy of this calculation is limited by the sampling frequency. | From the data in '''Fig. 2''', it is possible to carry out a spreadsheet calculation to estimate the recovery of the injected tracer in each individual production well. This shows that a large percentage of the injected gas has moved in different directions. Breakthrough time also may be estimated, either as the first sample that contains tracer above the detection limit or as an extrapolation backward of the initial phase of the production curve. However, the accuracy of this calculation is limited by the sampling frequency. | ||
To fully understand the flow behavior, it is important to add all information available. The picture may be modified by the knowledge of injection rates and production rates in the neighboring wells. | To fully understand the flow behavior, it is important to add all information available. The picture may be modified by the knowledge of injection rates and production rates in the neighboring wells. | ||
Because of maintenance and operational problems, production wells may be shut in for some periods. Because of fluid drift in the reservoir during shut-in, samples collected immediately after such periods may give unexpected results. This additional information can be exploited to better understand the pressure distribution in the field. Somaruga ''et al.''<ref name="r2" /> used well shut-in in a systematic way to obtain additional information. If the changes that may occur during shut-in are not considered, it may cause misinterpretation of the tracer-response curves. | Because of maintenance and operational problems, production wells may be shut in for some periods. Because of fluid drift in the reservoir during shut-in, samples collected immediately after such periods may give unexpected results. This additional information can be exploited to better understand the pressure distribution in the field. Somaruga ''et al.''<ref name="r2">_</ref> used well shut-in in a systematic way to obtain additional information. If the changes that may occur during shut-in are not considered, it may cause misinterpretation of the tracer-response curves. | ||
The water-tracer-response curves should be presented together with the water cut. '''Fig. 3''' shows water-tracer-response curves in the same well as discussed previously. Initially, a small water-cut increase is recorded. The tracer injected in Well P34 follows this increase. When the water cut starts to increase more rapidly, one of the other tracers starts production while the concentration of the first tracer vanishes. This method of presentation increases the understanding of the flow behavior and the water contribution from the different injectors. | The water-tracer-response curves should be presented together with the water cut. '''Fig. 3''' shows water-tracer-response curves in the same well as discussed previously. Initially, a small water-cut increase is recorded. The tracer injected in Well P34 follows this increase. When the water cut starts to increase more rapidly, one of the other tracers starts production while the concentration of the first tracer vanishes. This method of presentation increases the understanding of the flow behavior and the water contribution from the different injectors. | ||
<gallery widths="300px" heights="200px"> | <gallery widths="300px" heights="200px"> | ||
Line 38: | Line 39: | ||
</gallery> | </gallery> | ||
===Response curves=== | === Response curves === | ||
The time response is the graphic representation of the concentration of activity (after background subtraction and decay correction) as a function of time. Preprocessing the experimental data can smooth the response. From this curve, the cumulative response (recovered activity vs. time) is derived by a simple numeric integration. | |||
The time response is the graphic representation of the concentration of activity (after background subtraction and decay correction) as a function of time. Preprocessing the experimental data can smooth the response. From this curve, the cumulative response (recovered activity vs. time) is derived by a simple numeric integration. | |||
Concerning the cumulative response, the '''Eq. 1''' gives the cumulative tracer recovered up to an instant, ''t''<sub>''i''</sub>. | Concerning the cumulative response, the '''Eq. 1''' gives the cumulative tracer recovered up to an instant, ''t''<sub>''i''</sub>. | ||
[[File:Vol5 page 0672 eq 001.png]]....................(1) | [[File:Vol5 page 0672 eq 001.png|RTENOTITLE]]....................(1) | ||
where [[File:Vol5 page 0672 inline 001.png|RTENOTITLE]] = the total tracer recuperation up to ''t''<sub>''i''</sub>, [[File:Vol5 page 0672 inline 002-3.png|RTENOTITLE]] = the production water flow rate as a function of time, ''C''<sub>''t''</sub> = tracer concentration as a function of time, and ''t''<sub>''i''</sub> = the elapsed time (days after the injection). [[File:Vol5 page 0672 inline 002-3.png|RTENOTITLE]] often will be available as a discrete value representing each day. ''C''<sub>''t''</sub> will be available only as an individual measurement according to the sampling program. The accuracy of the total recovery will depend on the sampling frequency. | |||
=== Sweep volume === | |||
Instead of the time-response curves, the data are often presented as volumetric-response curves. The presentation is especially convenient for estimating the swept volume between the particular injector and producer pair. Swept volume, ''V''<sub>''S''</sub>, can be estimated from the product of the mean produced volume, [[File:Vol5 page 0673 inline 001-2.png|RTENOTITLE]], and the ratio of the mean volumetric water rate between the injector and the producer, ''Q''<sub>''m''</sub>, to the mean volumetric water production rate, ''Q''<sub>''p''</sub>.<ref name="r3">_</ref> | |||
Instead of the time-response curves, the data are often presented as volumetric-response curves. The presentation is especially convenient for estimating the swept volume between the particular injector and producer pair. Swept volume, ''V''<sub>''S''</sub>, can be estimated from the product of the mean produced volume, [[File:Vol5 page 0673 inline 001-2.png]], and the ratio of the mean volumetric water rate between the injector and the producer, ''Q''<sub>''m''</sub>, to the mean volumetric water production rate, ''Q''<sub>''p''</sub>.<ref name="r3" /> | |||
[[File:Vol5 page 0673 eq 001.png]]....................(2) | [[File:Vol5 page 0673 eq 001.png|RTENOTITLE]]....................(2) | ||
[[File:Vol5 page 0673 inline 001-2.png]] is calculated from the first moment of the produced-tracer concentration, ''C''. | [[File:Vol5 page 0673 inline 001-2.png|RTENOTITLE]] is calculated from the first moment of the produced-tracer concentration, ''C''. | ||
[[File:Vol5 page 0673 eq 002.png]]....................(3) | [[File:Vol5 page 0673 eq 002.png|RTENOTITLE]]....................(3) | ||
''Q''<sub>''m''</sub>, the flow rate between injector and producer, is estimated from the fraction of injected tracer produced at the well at the mean injection rate, ''Q''<sub>''i''</sub>. | ''Q''<sub>''m''</sub>, the flow rate between injector and producer, is estimated from the fraction of injected tracer produced at the well at the mean injection rate, ''Q''<sub>''i''</sub>. | ||
[[File:Vol5 page 0673 eq 003.png]]....................(4) | [[File:Vol5 page 0673 eq 003.png|RTENOTITLE]]....................(4) | ||
where ''m'' = the amount of tracer produced at a given well and ''M'' = the amount of tracer originally injected. The swept volume can be expressed in terms of the injected and produced flow rates and the fraction, ''m''/''M'', of water going from the injector to the producer. | where ''m'' = the amount of tracer produced at a given well and ''M'' = the amount of tracer originally injected. The swept volume can be expressed in terms of the injected and produced flow rates and the fraction, ''m''/''M'', of water going from the injector to the producer. | ||
[[File:Vol5 page 0674 eq 001.png]]....................(5) | [[File:Vol5 page 0674 eq 001.png|RTENOTITLE]]....................(5) | ||
Because of limited data, the tracer-production curve is composed of a discontinuous set of points and, in practice, the integral in '''Eq. 3''' is approximated by | Because of limited data, the tracer-production curve is composed of a discontinuous set of points and, in practice, the integral in '''Eq. 3''' is approximated by | ||
[[File:Vol5 page 0674 eq 002.png]]....................(6) | [[File:Vol5 page 0674 eq 002.png|RTENOTITLE]]....................(6) | ||
The calculation of swept volume often is complicated by a lack of data. The mean produced volume may need to be estimated on the basis of extrapolated tracer-production curves. In many cases, the sampling program is either not finished or it has been stopped without following the tail of the production curve. To fit the missing data in the tail of the production curve, an exponential decline approximation can be applied. | The calculation of swept volume often is complicated by a lack of data. The mean produced volume may need to be estimated on the basis of extrapolated tracer-production curves. In many cases, the sampling program is either not finished or it has been stopped without following the tail of the production curve. To fit the missing data in the tail of the production curve, an exponential decline approximation can be applied. | ||
[[File:Vol5 page 0674 eq 003.png]]....................(7) | [[File:Vol5 page 0674 eq 003.png|RTENOTITLE]]....................(7) | ||
where ''C''<sub>''e''</sub> = the measured tracer concentration at the value of V e from which the exponential fit starts. The slope of the line is represented by 1/''a''. If '''Eq. 7''' is substituted for ''C'' in '''Eq. 3''' and the integral replaced by a finite sum, the following expression is valid: | where ''C''<sub>''e''</sub> = the measured tracer concentration at the value of V e from which the exponential fit starts. The slope of the line is represented by 1/''a''. If '''Eq. 7''' is substituted for ''C'' in '''Eq. 3''' and the integral replaced by a finite sum, the following expression is valid: | ||
[[File:Vol5 page 0674 eq 004.png]]....................(8) | [[File:Vol5 page 0674 eq 004.png|RTENOTITLE]]....................(8) | ||
To find the swept volume in an incomplete data set, this equation is used for [[File:Vol5 page 0673 inline 001-2.png]] in '''Eq. 5'''. | To find the swept volume in an incomplete data set, this equation is used for [[File:Vol5 page 0673 inline 001-2.png|RTENOTITLE]] in '''Eq. 5'''. | ||
== Evaluation based on analytical solutions == | == Evaluation based on analytical solutions == | ||
'''Eq. 9''' gives the fundamental equation for tracer breakthrough in a streamtube. A correlation into a single breakthrough curve for various repeated homogeneous-flooding patterns has been found. <ref name="r4" /> In a layered system, the overall tracer-response curve is a combination of responses from individual layers. The individual-layer responses are predictable from the correlated pattern-breakthrough curve; however, the tracer arrival time at the production well and the concentration contribution from each layer are functions of the porosity, permeability, and thickness of that layer. Conversely, the decomposition of an actual tracer-production curve from a multilayered system into the individual-layer responses can yield individual-layer parameters. Computer algorithms exist that deconvolve the overall tracer profile into the individual-layer responses and compute ''ϕ'' h and ''kh''/Σ''kh'' of the individual layers. | '''Eq. 9''' gives the fundamental equation for tracer breakthrough in a streamtube. A correlation into a single breakthrough curve for various repeated homogeneous-flooding patterns has been found. <ref name="r4">_</ref> In a layered system, the overall tracer-response curve is a combination of responses from individual layers. The individual-layer responses are predictable from the correlated pattern-breakthrough curve; however, the tracer arrival time at the production well and the concentration contribution from each layer are functions of the porosity, permeability, and thickness of that layer. Conversely, the decomposition of an actual tracer-production curve from a multilayered system into the individual-layer responses can yield individual-layer parameters. Computer algorithms exist that deconvolve the overall tracer profile into the individual-layer responses and compute ''ϕ'' h and ''kh''/Σ''kh'' of the individual layers. | ||
[[File:Vol5 page 0668 eq 002.png]]....................(9) | [[File:Vol5 page 0668 eq 002.png|RTENOTITLE]]....................(9) | ||
== Evaluation of tracer-response curves with numerical simulation == | == Evaluation of tracer-response curves with numerical simulation == | ||
The most thorough analyses of tracer data are carried out in combination with a reservoir simulator applied for that particular field. Many of the applied reservoir simulators, based on numerical solution of flow equations applying finite-difference methods, have options for handling tracers. <ref name="r5" /> In addition, streamline models<ref name="r6" /><ref name="r7" /> exist that handle tracer flow. However, normal simulators have limitations. One problem is the limited possibilities for including physical properties such as partitioning between phases and adsorption to grain surfaces. This problem is of special importance for gas tracers. There is also a problem connected to handling of dispersion because the tracer pulse is, in its initial phase, much smaller than a typical gridblock size. | The most thorough analyses of tracer data are carried out in combination with a reservoir simulator applied for that particular field. Many of the applied reservoir simulators, based on numerical solution of flow equations applying finite-difference methods, have options for handling tracers. <ref name="r5">_</ref> In addition, streamline models<ref name="r6">_</ref><ref name="r7">_</ref> exist that handle tracer flow. However, normal simulators have limitations. One problem is the limited possibilities for including physical properties such as partitioning between phases and adsorption to grain surfaces. This problem is of special importance for gas tracers. There is also a problem connected to handling of dispersion because the tracer pulse is, in its initial phase, much smaller than a typical gridblock size. | ||
To overcome some of these problems, Sagen ''et al.''<ref name="r8">_</ref> developed a simulation module to improve the accuracy of the calculations. To resolve the influence of reservoir heterogeneities on the measured tracer response, an accurate numerical treatment of the tracer equation is needed. This is especially important when narrow tracer slugs are injected in a reservoir. This tracer module calculates tracer flow with an explicit method for integration of the convection part of the tracer equation. To reduce numerical smearing of the tracer pulses, the timestep for the tracer calculation is selected as large as possible but may still be much smaller than the timesteps of the host reservoir simulator, which most often uses an implicit formulation. In the tracer module, the main tool for reducing numerical dispersion is the use of a second-order numerical scheme for integrating the tracer equation. A separate grid-refinement option for tracer calculation is available. In combination, these methods lead to a good resolution of narrow slugs propagating through the reservoir. The method of separate grid refinement is far less time consuming than performing the whole reservoir simulation on a refined grid. This tracer simulation module has been coupled to some of the standard reservoir simulation tools. | |||
== Nomenclature == | |||
{| | {| | ||
|- | |- | ||
|'' | | ''a'' | ||
|= | | = | ||
| | | cross section, L<sup>2</sup> | ||
|- | |- | ||
|''C'' | | ''C'' | ||
|= | | = | ||
| | | concentration | ||
|- | |- | ||
|''C''<sub>'' | | ''C''<sub>''e''</sub> | ||
|= | | = | ||
|concentration | | measured tracer concentration | ||
|- | |- | ||
|''C''<sub>'' | | ''C''<sub>''O''</sub> | ||
|= | | = | ||
| | | concentration of tracer in oil phase | ||
|- | |- | ||
|'' | | ''C''<sub>''t''</sub> | ||
|= | | = | ||
| | | tracer concentration as a function of time | ||
|- | |- | ||
| | | ''m'' | ||
|= | | = | ||
| | | amount of tracer produced in a given well | ||
|- | |- | ||
| | | [[File:Vol5 page 0672 inline 001.png|RTENOTITLE]] | ||
|= | | = | ||
| | | total tracer recuperation up to ''t''<sub>''i''</sub> | ||
|- | |- | ||
|'' | | ''L'' | ||
|= | | = | ||
| | | distance, L | ||
|- | |- | ||
|'' | | ''M'' | ||
|= | | = | ||
| | | amount of tracer originally injected | ||
|- | |- | ||
|''Q''<sub>'' | | ''Q''<sub>''m''</sub> | ||
|= | | = | ||
|mean volumetric water | | mean volumetric water rate between injector and producer, L<sup>3</sup>/t | ||
|- | |- | ||
| | | ''Q''<sub>''p''</sub> | ||
|= | | = | ||
|production | | mean volumetric water production rate, L<sup>3</sup>/t | ||
|- | |- | ||
| | | Qpt | ||
|= | | = | ||
| | | production-water flow rate as a function of time, L<sup>3</sup>/t | ||
|- | |- | ||
|''s''<sub> | | ''s''<sub>1</sub> | ||
|= | | = | ||
| | | front location of tracer pulse in the streamline, L | ||
|- | |- | ||
|'' | | ''s''<sub>2</sub> | ||
|= | | = | ||
| | | end location of tracer pulse in the streamline, L | ||
|- | |- | ||
|'' | | ''t''<sub>''i''</sub> | ||
|= | | = | ||
| | | elapsed time for tracer ''i'', t | ||
|- | |- | ||
| | | ''V'' | ||
|= | | = | ||
| | | parameter that decides the shape of the curve for the actual system at a specific temperature | ||
|- | |- | ||
| | | [[File:Vol5 page 0673 inline 001-2.png|RTENOTITLE]] | ||
|= | | = | ||
|produced volume | | mean produced volume, L<sup>3</sup> | ||
|- | |- | ||
|''V''<sub>'' | | ''V''<sub>''e''</sub> | ||
|= | | = | ||
| | | produced volume from which the exponential fit starts, L<sup>3</sup> | ||
|- | |- | ||
|'' | | ''V''<sub>''S''</sub> | ||
|= | | = | ||
|standard | | retention volume of standard reference tracer, L<sup>3</sup> | ||
|- | |- | ||
|''ψ'' | | ''σ'' | ||
|= | | = | ||
|streamtube | | standard deviation | ||
|- | |||
| ''ψ'' | |||
| = | |||
| streamtube | |||
|} | |} | ||
==References== | == References == | ||
< | <references /> | ||
== Noteworthy papers in OnePetro == | |||
Use this section to list papers in OnePetro that a reader who wants to learn more should definitely read | |||
== External links == | |||
Use this section to provide links to relevant material on websites other than PetroWiki and OnePetro | |||
== | == See also == | ||
[[Well_to_well_tracer_tests|Well to well tracer tests]] | |||
[[Tracer_flow_in_porous_reservoir_rock|Tracer flow in porous reservoir rock]] | |||
[[ | |||
[[ | [[Planning_and_design_of_well_to_well_tracer_tests|Planning and design of well to well tracer tests]] | ||
[[ | [[Field_experience_with_well_to_well_tracer_tests|Field experience with well to well tracer tests]] | ||
[[ | [[PEH:Well-To-Well_Tracer_Tests]] | ||
[[ | [[Category:5.6.5 Tracer test analysis]] |
Revision as of 11:26, 4 June 2015
To obtain high-quality tracer-response curves that are the basis for the further interpretation, a well-designed sampling program is needed. In general, more samples will give the potential for extraction of more information from field tests. Too often, interpretation is difficult because of limited tracer data.
The final objective of a well-to-well study is the interpretation of the response curves. A good analysis of the information given by the tracers, in combination with other available data, gives a better understanding of the flow in the reservoir, not just verification of communication between injector and producer.
Different levels of interpretation
The response curve can be analyzed from three different points of view or complexity levels. The simplest interpretation is the qualitative one. By observing the curves, the following pattern characteristics can be derived: injection-water arriving time (breakthrough); existence of high-permeability channels, barriers, and fractures between wells; communication between different layers; stratification in the same layer; and preferential flow directions in the reservoir. Furthermore, the cumulative response can be obtained by integrating the concentration vs. time curve if the production flow rate is known. From this new curve, the fraction of injection water reaching each producer is easily calculated. A common spreadsheet is the best way to make these calculations. This type of interpretation can be carried out without any advanced simulation. It is important to integrate the data obtained from the geological model, primary production data, 4D seismic data (if available), and tracer data in a systematic way.
The second level of analysis uses basic mathematical models to fit simple response curves by means of theoretical expressions and to decompose complex responses in several simpler functions. In this way, partial residence times, as well as other parameters, can be determined for each function. Mathematical models also allow for the evaluation of important parameters, such as permeability, and make it possible to predict the behavior of unknown patterns.
Finally, complex mathematical models like numerical finite-element simulators or streamtube modeling can be used to achieve a deeper analysis. Most of the commercially available simulators have tracer options with varying degrees of complexity. Not all simulators include important physical effects like partitioning, dispersion, and adsorption. A tracer simulator that can be coupled to full-field reservoir simulators has been developed recently. Simulation is a tool to improve the existing reservoir model; therefore, it is crucial to have access to the best available model to enable an efficient optimization based on the tracer-production data.
Evaluation based on hand calculation and adequate data presentation
Qualitative interpretation of field data is illustrated by examples from the Snorre field in the North Sea. Fig. 1 shows the segment of the field that was tested and the location of injectors (triangles) and producers (circles). [1] A water-alternating-gas (WAG) injection program is being used in the field.
Fig. 1 – Well location in the central fault block (CFB) on the Snorre field. Arrows indicate the main producers of injected tracers.[1]
Fig. 2 shows a way to present field data for a gas-tracer program. The data are plotted as a time response. Alternatively, the data can be plotted as a volumetric response. To obtain a better understanding, it is important to plot gas-tracer-production data with GOR. The tracer data are plotted relative to the amount of tracer injected, which means that the concentrations presented for each individual tracer can be compared directly.
Fig. 2 – Gas tracer data from the Snorre field in the North Sea.[1]
From Fig. 2, it is possible not only to find communication but also to extract information about the prevailing flow regime. The gas tracer injected into Well P25 breaks through without any significant increase in GOR, which means that the gas has not followed a free gas layer but has been dissolved in the oil. Three different gas tracers were injected in Well P28 at three different times. These tracers were produced only after a significant GOR was recorded. All the tracers were produced at almost the same time. This shows that gas has been trapped and that the gas in the individual gas slugs has subsequently filled this trap. At a certain time, the gas trap has been filled and the gas moves on to the producers. When reaching the producer, the gas contains all three different tracers injected over a 3-year period.
From the data in Fig. 2, it is possible to carry out a spreadsheet calculation to estimate the recovery of the injected tracer in each individual production well. This shows that a large percentage of the injected gas has moved in different directions. Breakthrough time also may be estimated, either as the first sample that contains tracer above the detection limit or as an extrapolation backward of the initial phase of the production curve. However, the accuracy of this calculation is limited by the sampling frequency.
To fully understand the flow behavior, it is important to add all information available. The picture may be modified by the knowledge of injection rates and production rates in the neighboring wells.
Because of maintenance and operational problems, production wells may be shut in for some periods. Because of fluid drift in the reservoir during shut-in, samples collected immediately after such periods may give unexpected results. This additional information can be exploited to better understand the pressure distribution in the field. Somaruga et al.[2] used well shut-in in a systematic way to obtain additional information. If the changes that may occur during shut-in are not considered, it may cause misinterpretation of the tracer-response curves.
The water-tracer-response curves should be presented together with the water cut. Fig. 3 shows water-tracer-response curves in the same well as discussed previously. Initially, a small water-cut increase is recorded. The tracer injected in Well P34 follows this increase. When the water cut starts to increase more rapidly, one of the other tracers starts production while the concentration of the first tracer vanishes. This method of presentation increases the understanding of the flow behavior and the water contribution from the different injectors.
Fig. 3 – Water-tracer-response curve from the Snorre field in the North Sea.[1]
Response curves
The time response is the graphic representation of the concentration of activity (after background subtraction and decay correction) as a function of time. Preprocessing the experimental data can smooth the response. From this curve, the cumulative response (recovered activity vs. time) is derived by a simple numeric integration.
Concerning the cumulative response, the Eq. 1 gives the cumulative tracer recovered up to an instant, ti.
where = the total tracer recuperation up to ti, = the production water flow rate as a function of time, Ct = tracer concentration as a function of time, and ti = the elapsed time (days after the injection). often will be available as a discrete value representing each day. Ct will be available only as an individual measurement according to the sampling program. The accuracy of the total recovery will depend on the sampling frequency.
Sweep volume
Instead of the time-response curves, the data are often presented as volumetric-response curves. The presentation is especially convenient for estimating the swept volume between the particular injector and producer pair. Swept volume, VS, can be estimated from the product of the mean produced volume, , and the ratio of the mean volumetric water rate between the injector and the producer, Qm, to the mean volumetric water production rate, Qp.[3]
is calculated from the first moment of the produced-tracer concentration, C.
Qm, the flow rate between injector and producer, is estimated from the fraction of injected tracer produced at the well at the mean injection rate, Qi.
where m = the amount of tracer produced at a given well and M = the amount of tracer originally injected. The swept volume can be expressed in terms of the injected and produced flow rates and the fraction, m/M, of water going from the injector to the producer.
Because of limited data, the tracer-production curve is composed of a discontinuous set of points and, in practice, the integral in Eq. 3 is approximated by
The calculation of swept volume often is complicated by a lack of data. The mean produced volume may need to be estimated on the basis of extrapolated tracer-production curves. In many cases, the sampling program is either not finished or it has been stopped without following the tail of the production curve. To fit the missing data in the tail of the production curve, an exponential decline approximation can be applied.
where Ce = the measured tracer concentration at the value of V e from which the exponential fit starts. The slope of the line is represented by 1/a. If Eq. 7 is substituted for C in Eq. 3 and the integral replaced by a finite sum, the following expression is valid:
To find the swept volume in an incomplete data set, this equation is used for in Eq. 5.
Evaluation based on analytical solutions
Eq. 9 gives the fundamental equation for tracer breakthrough in a streamtube. A correlation into a single breakthrough curve for various repeated homogeneous-flooding patterns has been found. [4] In a layered system, the overall tracer-response curve is a combination of responses from individual layers. The individual-layer responses are predictable from the correlated pattern-breakthrough curve; however, the tracer arrival time at the production well and the concentration contribution from each layer are functions of the porosity, permeability, and thickness of that layer. Conversely, the decomposition of an actual tracer-production curve from a multilayered system into the individual-layer responses can yield individual-layer parameters. Computer algorithms exist that deconvolve the overall tracer profile into the individual-layer responses and compute ϕ h and kh/Σkh of the individual layers.
Evaluation of tracer-response curves with numerical simulation
The most thorough analyses of tracer data are carried out in combination with a reservoir simulator applied for that particular field. Many of the applied reservoir simulators, based on numerical solution of flow equations applying finite-difference methods, have options for handling tracers. [5] In addition, streamline models[6][7] exist that handle tracer flow. However, normal simulators have limitations. One problem is the limited possibilities for including physical properties such as partitioning between phases and adsorption to grain surfaces. This problem is of special importance for gas tracers. There is also a problem connected to handling of dispersion because the tracer pulse is, in its initial phase, much smaller than a typical gridblock size.
To overcome some of these problems, Sagen et al.[8] developed a simulation module to improve the accuracy of the calculations. To resolve the influence of reservoir heterogeneities on the measured tracer response, an accurate numerical treatment of the tracer equation is needed. This is especially important when narrow tracer slugs are injected in a reservoir. This tracer module calculates tracer flow with an explicit method for integration of the convection part of the tracer equation. To reduce numerical smearing of the tracer pulses, the timestep for the tracer calculation is selected as large as possible but may still be much smaller than the timesteps of the host reservoir simulator, which most often uses an implicit formulation. In the tracer module, the main tool for reducing numerical dispersion is the use of a second-order numerical scheme for integrating the tracer equation. A separate grid-refinement option for tracer calculation is available. In combination, these methods lead to a good resolution of narrow slugs propagating through the reservoir. The method of separate grid refinement is far less time consuming than performing the whole reservoir simulation on a refined grid. This tracer simulation module has been coupled to some of the standard reservoir simulation tools.
Nomenclature
References
Noteworthy papers in OnePetro
Use this section to list papers in OnePetro that a reader who wants to learn more should definitely read
External links
Use this section to provide links to relevant material on websites other than PetroWiki and OnePetro
See also
Tracer flow in porous reservoir rock
Planning and design of well to well tracer tests