Low primary and secondary recoveries of original oil in place from modern unconventional reservoirs beg for utilization of tertiary recovery techniques. Enhanced oil recovery (EOR) via cyclic gas injection (“huff ‘n’ puff”) has indeed enhanced the oil recovery in many fields, and many of those projects have also been documented in industry technical papers/case studies. However, the need remains to document new techniques in new reservoirs. This paper documents a small-scale EOR pilot project in the eastern Eagle Ford and shows promising well results.

In preparation for the pilot, full characterization of the oil and injection gas was done along with laboratory testing to identify the miscibility properties of the two fluids. Once the injection well facility design was completed, a series of progressively larger gas volumes were injected followed by correspondingly longer production times. Fluids in the returning liquid and gas streams were monitored for compositional changes, and the learnings from each cycle led to adjustments and facility changes to improve the next cycle.

After completing five injection/withdrawal cycles in the pilot, a few key observations can be made. The implementation of cyclic gas injection can be both a technical and a commercial success early in its life if reasonable cost controls are implemented and the scope is kept manageable. The process has proved to be both repeatable and predictable, allowing for future economic modeling to be used to help determine timing of subsequent injection cycles. A key component of the success of this pilot has been the availability of small compressors capable of the high pressures required for these projects and learning how to implement cost saving facility designs that still meet high safety standards.

The Lower Eagle Ford in the northern Eagle Ford extension, north of the San Marcos Arch, is similar in age and time of deposition to the Eagle Ford in south Texas. Across this basin, we see higher siliceous content and lower carbonate contents when compared to the Eagle Ford deposition south of the Arch (Vallabhaneni et al. 2016). The subject well is on the northern end of this extension, located in Madison County, Texas, USA. Well logs, core, and production results show that the reservoir quality across the basin is quite comparable. The target reservoir in Madison County is a siliceous mudstone that has a total organic carbon content in the range of 3 to 8% as observed in cores. Average reservoir pressure gradients place it in the overpressured category at 0.71 psi/ft with a bottomhole static temperature of ~235°F. Oil produced ranges from 39° to 42° API with solution gas/oil ratios from 500 to 700 scf/bbl. Long-term water production is low, settling at ~20% of total fluid production within a few months of the start of production.

Development in the Lower Eagle Ford has accelerated in the past few years as large public exploration and production operations made their entrance into the subbasin. Before the advent of horizontal drilling, there were many penetrations in vertical wells, but there were very few completions due to the characteristically low permeability, and those that were attempted were uneconomic. As operators came to this area attracted by the potential of horizontal development of the tight conventional Woodbine sands, which lie 400 to 600 ft true vertical depth above the Eagle Ford, they also began to experiment with the Eagle Ford shale. The progression of technique for completion has been the primary driver in the improving results from the reservoir over time, a very similar pattern to the one we have seen in every major horizontal shale basin in the US. The standard progression of longer laterals, large proppant to small, gelled fluids to slick, increasing sand volumes, and decreasing stage spacing has generally led to higher initial production and estimated ultimate recoveries in the search for higher returns.

Miscible gasflooding is a proven viable process that enhances oil recovery in conventional and low-permeability reservoirs (Rao 2001; Gozalpour et al. 2005; Sohrabi et al. 2008; Kovscek et al. 2008; Arshad et al. 2009; Wang and Gu 2011; Anand et al. 2019; Sie and Nguyen 2020a, 2020b; Liu et al. 2022). Several laboratory investigations and field tests of the efficiency of different gas injection strategies for EOR in shale formations have been reported in the literature over the past few years (Hawthorne et al. 2014; Sheng 2015; Hoffman and John 2016; Alfarge et al. 2017, 2018; Tovar et al. 2018a; Hamdi et al. 2018, 2021; Hoffman 2018; Jia et al. 2019; Carlsen et al. 2020a, 2020b; Huang et al. 2022). The hydrocarbon gas huff ‘n’ puff process was found to be the most promising strategy based on both reservoir simulation and field-scale pilot testing (Alfarge et al. 2018; Hoffman and John 2016; Hoffman 2018). However, the use of field gas, as in shale huff ‘n’ puff, has not been well explored because, unlike CO2 which is another commonly used huff ‘n’ puff, field gas is a mixture of hydrocarbon components that have different interactions with organic matter–rich shales, and thus a laboratory evaluation of cyclic injection of field gas in shales would be difficult. It is commonly observed that the composition of natural gas has been oversimplified in most previous laboratory huff ‘n’ puff studies (Li et al. 2017a; Hawthorne et al. 2017).

An experimental study of huff ‘n’ puff processes by Hawthorne et al. (2017) showed that both CO2 and ethane could recover nearly all of the oil initially in Bakken shale cores. The authors of that study also concluded that the use of methane above its minimum miscibility pressure with the crude oil reduced oil recovery by only 10%. It has also been well established that enriching methane with natural gas liquids in huff ‘n’ puff could significantly reduce the minimum miscibility pressure and enhance oil extraction (Tran et al. 2021). Song et al. (2021) conducted huff ‘n’ puff experiments on Upper Devonian Duvernay core samples and observed that adding propane to the injected mixture of methane and ethane resulted in an increase in oil recovery. However, the oil recovery factor does not monotonically increase with the injection pressure. The experimental work by Thomas et al. (2020a, 2020b) showed that the increase of recovery factor for pure methane cyclic injection became marginal as the huff pressure was raised above an optimal value. The compositional complexity of crude oil could also strongly impact the performance of a gas huff ‘n’ puff process. For example, the oil recovery in shale samples with methane cyclic injection could be as high as 80% of gas condensate but reduce by more than 30% for a volatile oil (Thomas et al. 2020a). The impact of the initial solution gas before gas injection has also been found important (Mahzari et al. 2021)

The industry has long understood that shale reservoirs have characteristically low primary recoveries (US Energy Information Administration 2016). This was confirmed by the petrophysical work done during primary development. As the companies in the area searched for ways to improve the overall economics of shale development, this area was identified as a good candidate to test the effectiveness of huff ‘n’ puff recovery. Due to the low development density in the subject area, many of the well-to-well interference problems experienced in overdeveloped reservoirs could be avoided. Another significant factor in the decision to implement the project was the availability of relatively low-cost rich injection gas that was controlled by the operator of the well. To determine in which well to first implement the process, the available wells were screened, and wells that had low-intensity initial completions or that were producing more than 100 barrels of oil per day (BOPD) were eliminated. Injecting into a well producing more than 100 BOPD was an economic risk the company was not willing to take, as the outcome of the experiment was unknown and it was feared that the process may damage the wellbore or reservoir. In screening for a suitable example well, it was recognized that the size of the primary estimated ultimate recovery would probably be indicative of the potential recoveries that could be expected from an EOR process. In the final selection, the Montana 1 H (42-313-31220) was the oldest well that had a full size fracturing job applied to it, had not demonstrated any well interference with nearby wells during stimulation, and was not near any known faults. At the time the first injection began, the subject well had been producing for 4.3 years and had a cumulative production of 152,000 bbl of crude oil, 47 million ft3, and 48,000 bbl of water. The minimum breakdown pressure observed while stimulating this well was 8,087 psi as observed through the instantaneous shut-in pressures after each stimulation. Measurements taken after the well primary production began indicate a reservoir pressure of approximately 6,174 psi.

Once a decision was made to implement the pilot, experiments were designed to understand the fluid-fluid interactions that would guide many of the design and sizing decisions for the pilot. Once the physical interactions were understood, engineering design sized and scoped appropriate equipment for the pilot. It was important that the equipment selected would use the known parameters and allow field experiments to be designed that would allow for quick learning cycles and modification. The remainder of this paper will serve to describe the process of experimental design both at the laboratory and field level and to highlight some significant results and observations.

Materials

The injection gas used for all laboratory tests in this work consists of 66.20 mol% methane, 18.65 mol% ethane, 10.48 mol% propane, 3.55% n-butane, and less than 1 mol% n-pentane. A core plug sample (2 in. long and 1 in. in radius) taken from a shale formation was used in the gas huff ‘n’ puff experiments. The permeability of this core plug is expected to be very close to 0.05 md. This permeability was obtained from a permeability test (i.e., the rock-eval pyrolysis) using a twin core plug taken from the same well at nearly the sample depth. It is uncertain that this measured permeability includes fracture permeability. The core porosity of 4.98% was determined based on the bulk volume of the core plug and the amount of dead crude oil saturated in the core plug before the field gas extraction experiment. The specific gravity of the dead crude oil used in this work is 0.8167 at 60°F. This oil sample was taken from a shale oil formation producing at a very low gas/oil ratio. It is noted that the use of dead crude oil in the huff ‘n’ puff experiments in this work would not precisely capture the interaction of actual crude oil at very low gas/oil ratio in the reservoir with the injection gas. The composition of the dead crude oil sample was measured by standard condensate analysis (GPA-2103) and is presented in the  Appendix.

Minimum Capillary Pressure Test

Fig. 1a shows a schematic of the experimental setup used to determine the minimum capillary pressure (MCP). The improvements on this design were based on the work by Hawthorne et al. (2016). Three capillary tubes (World Precision Instruments) with inner diameters of 1.12, 0.84, and 0.68 mm, respectively, were used in the experiment. The capillary tubes and a piston were placed into a custom-designed high-pressure sapphire visualization cell. The crude oil with the desired volume was transferred into the cell while the piston was at the bottom of the cell and the cell was heated up to the reservoir temperature (235°F). The synthetic field gas was then injected into the cell until the system reached a target pressure. The hydrocarbon liquid column heights in the capillary tubes (e.g., Figs. 1b and 1c) were monitored and recorded by a camera. The MCP is an estimate of the minimum pressure above which the liquid imbibition into the capillary tubes is completely inhibited because the interfacial tension between the gas and the oil at pressures above the MCP becomes so small that the capillary pressure could be negligible (De Gennes et al. 2004). Thus, the effect of system pressure on the oil extraction during a huff ‘n’ puff process could be better evaluated without a significant effect of capillary pressure (Sie and Nguyen 2020c).

Fig. 1

(a) Capillary rise experiment setup, (b) a photo of the capillary tube at low pressure, and (c) a photo of the capillary tubes when full miscibility was almost established.

Fig. 1

(a) Capillary rise experiment setup, (b) a photo of the capillary tube at low pressure, and (c) a photo of the capillary tubes when full miscibility was almost established.

Close modal

Oil Extraction Experiment

The residual oil in the reservoir core was removed using the gas mixture and the following procedure. The core was placed in a high-pressure cell in the oven, and then the system was vacuumed for an hour. After the system reached 235°F, the gas was injected into the cell until the pressure of the system reached 4,500 psia. This pressure was selected based on the MCP as discussed in the Results and Discussion section. The core was soaked for 3 days at this pressure to remove the residual oil and then placed into a glass vial at atmospheric pressure for another 12 hours to remove the residual gas. Next, the clean core was saturated with the dead crude oil and soaked in the synthetic field gas in the high-pressure cell at 235°F and 4,500 psia. The same gas mixture was used to flush the cell at 4,500 psia after a 24-hour soaking period to remove the produced oil from the cell. This soak ‘n’ flush process was repeated until no further oil production was observed.

The beginning points of the facility design included consideration of both the pressure required to stimulate the reservoir as well as an acceptable flow rate of gas into the reservoir. Due to the cyclic nature of this process, a mobile equipment design was specified. Above average gas pressures required for the process demanded an in-depth review of equipment and procedures to ensure safety. The details of each of these points and the impact on the facility design are discussed.

Pressure Considerations

The pressure required to stimulate the reservoir needed to be above the MCP of approximately 4,500 psia (refer to Results and Discussion), but significantly below the breakdown pressure of 8,087 psig. These pressures are substantially above conventional field gas pressures (1,200 psig) and rapidly demanded a thorough search for unconventional compression equipment as traditional field gas compressors are limited to approximately 1,200 psig. Compressors achieving higher pressures have historically been almost exclusively purpose-designed and -fabricated large stationary compressors (Jacobs 2019). However, a relatively new compressor design capable of developing pressures of up to 5,500 psig for a relatively new artificial lift method, high-pressure gas lift , entered the market and offered promise. The compressor is a critical piece of equipment for the project, influencing a number of important factors in the success of this project.

Gas Flow Rate Considerations

Gas flow rates to be injected into the reservoir were bound by two factors: On the low end is gas compressor throughput, which was approximately 1.5 MMscf/D, and on the high end is gas availability from neighboring wells, which was limited to approximately 4 MMscf/D. To minimize the time to pressurize the reservoir, two compressors were installed in parallel. Acknowledging this project as a pilot and the potential for learning and needed changes, the idea of scalability (both up and down) was highly desirable.

Equipment Portability

The desire existed in this pilot to try the gas EOR process across a number of various wells in the field. In lieu of building an expensive and long gas gathering and/or distribution system, utilizing mobile equipment was deemed necessary. Skidded equipment (compressors, production separators, flow meters, etc.) must be movable by flatbed trucks with rolling tailboards and winches as opposed to requiring cranes. Connection of the compressor process gas system to the wellhead was done with 1502 flow iron, which provided the practical benefit of eliminating the need for welding and hydrotesting operations during installation. Standard seals were replaced with Viton seals to raise the temperature limit of the flow iron to over 300°F. In keeping with the rapid mobilization/demobilization theme, utility connections, such as scrubber dump lines were made via flexible high-pressure hydraulic hoses. In the end, a modular and easily transportable equipment approach was taken and became valuable to the economic viability of the project.

Facility Safety

A process hazard analysis or hazard and operability analysis is insightful in any new project; however, given the above average pressures here, these analyses were deemed imperative and revealed some necessary safeguards. Because lower pressure equipment previously existed, in a move of extreme caution, blinds were installed on the valving in the production tree wing to isolate the production facility from high pressures during injection phases.

Given the high pressures combined with the high Btu value of the gas, hydrate formation was also a concern. Referencing the adjacent hydrate formation chart (Fig. 2), it can readily be seen that as gas pressures and gas specific gravities increase, so does the temperature at which hydrates can form. In process pressures of 4,000 psia and gas specific gravities of 0.8, hydrates can form in excess of 75°F. Thus, hydrate formation is conceivable at even springtime and falltime ambient temperatures. Elevating gas temperatures and insulating pipes to maintain elevated temperatures in gas lines becomes critical for reliable and safe operations. In cold weather months, use of methanol as an antifreeze agent is a further aid and added insurance.

Fig. 2

Phase behavior diagram for gas hydrate formation.

Fig. 2

Phase behavior diagram for gas hydrate formation.

Close modal

Compressor Design

The mobility of and, specifically, the rapid deployment/redeployment capabilities of the compressor skid are key. A self-contained compressor package that houses oil tanks and fuel supply systems on-skid reduces set time and costs. Requirements for off-skid utilities/services (i.e., instrument air supply, electricity, etc.) should also be eliminated. A truly mobile package should only require three off-skid connections: (i) inlet gas, (ii) gas discharge, and (iii) scrubber drains.

As mentioned above in the Facility Safety section and further discussed elsewhere (Elmer and Elmer 2017), elevated gas temperatures are needed to help manage hydrates. The compression process inherently generates heat, which can be useful. Most compression packages contain oversized gas cooling, negating the positive impact the compression process can contribute to managing hydrates. In EOR applications, gas injection temperatures tend only to be limited by component limitations, which in many cases is the wellhead, for which a 250°F limit is quite typical. As such, little to no process gas cooling in the compression process is oftentimes possible and is the case with the compression equipment selected for this project.

A related matter is dealing with hydrates on the compression skid itself, specifically in the compressor utility systems. Scrubber dump systems are notorious for creating hydrates during the dump process and benefit from heat tracing and insulation. Further, the fuel gas system, now that it is on-skid and with its large pressure drops (and consequently cooling effects brought on from the Joule-Thompson effect), also warrants added heat from heat tracing. This heat tracing is easily provided by heat exchangers and use of coolant from the engine jacket water system.

Gas-Oil Interaction and Efficiency of Oil Extraction in Reservoir Core

Fig. 3 shows the results of the capillary rise experiment. The MCP was 4,429.76 ± 138.89 psia, which was estimated by extrapolating the linear regression lines to the pressure at which the oil column height is zero (Hawthorne et al. 2016). When the field gas becomes more miscible with the crude oil, the interfacial tension between them approaches zero (Rao and Lee 2003; Hawthorne et al. 2017). Thus, the MCP is expected to be close to the minimum miscibility pressure. It is noted that the minimum miscibility pressure is often determined by either the slimtube experiment or using an equation of state, which is very different from the capillary raise method for the MCP measurement described above. It is commonly used to evaluate the efficiency of multiple-contact miscibility processes in conventional gasflooding. Based on this result, an oil extraction experiment was conducted at 235°F and 4,500 psia. This experiment confirmed that more than 80% of the initial dead oil in place could be recovered at this reservoir condition. It was also observed that the preponderance of oil recovery occurred in the first two days of extraction.

Fig. 3

Oil column height in three capillary tubes [1.5, 1.2, and 1 mm in outer diameter (OD)] as a function of pressure.

Fig. 3

Oil column height in three capillary tubes [1.5, 1.2, and 1 mm in outer diameter (OD)] as a function of pressure.

Close modal

During the huff stage, the injection gas in the fracture invades the matrix due to a large pressure difference between the fracture and the matrix and thus creates gas-oil mixing zones in the matrix. Thermodynamic equilibration redistributes the phase and composition in the mixing zone (Anand et al. 2019; Tang et al. 2020). During the soaking stage, thermodynamic equilibration redistributes the components at the fluid interface in the mixing zone according to the local composition. The injection gas components dissolve into the crude oil-rich phase and migrate deeper into the matrix by concentration-driven molecular diffusion. Meanwhile, the crude oil components are brought to the vapor phase and diffuse from the matrix into the fracture (Guo et al. 2009; Tharanivasan et al. 2006). During the pressure depletion (puff) stage, as the pressure in the fracture reduces to the puff pressure, the resulting pressure gradient induces the convective flow. Also, increasing gas expansion as pressure continues to drop could carry more crude oil from the matrix to the fracture (Mahzari et al. 2019, 2020). When the huff ‘n’ puff process is operated at high interfacial tension (e.g., below the MCP), the lower huff pressure provides less driving forces for the mass transfer between the oil and the injection gas (Tang et al. 2020; Sie et al. 2018, 2019; Anand et al. 2019). As the interfacial tension is reduced at elevated huff pressure (e.g., above the MCP), this mass transfer is greatly enhanced.

Field Cyclic Gas Injection Results

The pilot was designed with a staged approach. Keeping the early injection cycles small minimized the impact of reservoir containment loss and the impact to nearby producers should injection gas break through unexpectedly. The early cycles were kept small as the following series of questions was answered experimentally: What will the pressure response to injection look like? Is there an appreciable fluid composition change between injection and production? How much of the injected gas will be recovered? How fast will oil recovery happen?Table 1 summarizes the design parameters for five gas injection cycles and the associated production responses. Note that the incremental oil recovery was estimated based on comparing the oil rate at the end of each cycle with the oil rate projected by a decline curve analysis had the primary production continued.

Table 1

Injection parameters and production responses for five gas injection cycles.

CyclePeak Surface Pressure
(psi)
Peak Bottomhole Pressure
(psi)
Injected Volume
(MMscf)
Gas Recovery
(%)
Net Gas Efficiency
(Mscf/bbl)
Incremental Oil Recovery
(bbl)
Projected Oil Rate*
(BOPD)
Actual Oil Rate*
(BOPD)
Days on InjectionDays on Production
3,880 4,739 79 2.03 836 36 58 10 25 
4,325 5,169 16 75 5.33 818 35 46 16 30 
4,423 5,333 20 81 5.20 804 33 40 23 111 
5,420 5,417 97 62 5.64 6,510 29 36 63 291 
5,450 5,449 121 66.8+ 4.27 9,383 26 29 58 384 
CyclePeak Surface Pressure
(psi)
Peak Bottomhole Pressure
(psi)
Injected Volume
(MMscf)
Gas Recovery
(%)
Net Gas Efficiency
(Mscf/bbl)
Incremental Oil Recovery
(bbl)
Projected Oil Rate*
(BOPD)
Actual Oil Rate*
(BOPD)
Days on InjectionDays on Production
3,880 4,739 79 2.03 836 36 58 10 25 
4,325 5,169 16 75 5.33 818 35 46 16 30 
4,423 5,333 20 81 5.20 804 33 40 23 111 
5,420 5,417 97 62 5.64 6,510 29 36 63 291 
5,450 5,449 121 66.8+ 4.27 9,383 26 29 58 384 
*

At the end of cycle.

Injection Pressure Responses

The first injection cycles did demonstrate that injection gas was contained and recoverable in large enough quantities to allow the volumes on subsequent cycles to be increased. It was apparent that the injection pressure would have to be raised if large volumes were to be used. The containment was indicated by the consistent rise in injection pressure as injection continued during each cycle. Additionally, the injection well is directly underneath a parallel horizontal well in a sand formation ~500 ft true vertical depth above the injection well and no increase in gas production was observed from this or any surrounding wells.

The facility was operated at its maximum injection capacity throughout the pilot. In the beginning of each injection, the available gas rate from the field limited the injection rate; in the middle of each period, the experiment was limited by the horsepower available for compression; and toward the end, the experiment was limited by the maximum allowable operating pressure of the surface facilities. The injection behavior across all five cycles has been remarkably stable and repeatable (Fig. 4). With the initial 20 MMscf of gas pressurizing up the wellbore and near-wellbore region before establishing a relatively stable injection pressure, much larger volumes of gas have been forced into the far-field regions. Some spread between Cycles 4 and 5 can be observed in the injection pressures, though these are largely due to differences in the daily injection rate at various times. While a high injection rate was planned to be held constant until maximum pressures were achieved, this was not always possible, though operational improvements were made with each cycle that improved injection throughput and equipment run times. Basic design improvements such as adding a slug catcher to catch any liquids coming into the system off the pipeline and strategically placing methanol injection throughout the dump piping were very helpful.

Fig. 4

Calculated bottomhole pressure (BHP) vs. injected gas volume for five gas injection cycles. The initial reservoir pressure (P*) is also shown.

Fig. 4

Calculated bottomhole pressure (BHP) vs. injected gas volume for five gas injection cycles. The initial reservoir pressure (P*) is also shown.

Close modal

Produced Oil Gravity

The first three cycles were pressure limited by the maximum allowable operating pressure of the system at the time, with the last two cycles being volume limited after the improvement in the system maximum allowable operating pressure. Once an injection cycle was ended, the wellhead and flow lines had their high-pressure sides blinded off and the production side blinds removed. The well was then flowed at a constant gas rate until produced liquids hit. Due to the constraints of the existing production facility, the maximum sustainable flowback rate was ~1 MMscf/D, during which time a line heater was used along with methanol and an extra medium pressure separator to minimize hydrate formation. Once liquids began to hit, a fixed choke was used to allow pressure to slowly continue to bleed off. The transition from rate limited flow to fixed choke flow and finally to starting the gas lift system when tubing pressures got below 200 psig allowed a gradual bottomhole pressure drawdown to be achieved. Initial flowback during each withdrawal (puff) cycle was 100% injection gas, with liquids gradually increasing in volume as the pressures decreased. Oil samples showed that while the early oils made were quite light, these oils were being dropped out of the gas stream due to the extreme pressure drop from the wellhead pressure to the facility pressures early in the flowback. The API gravity and color change in the produced fluids can be observed in Fig. 5. Once the oil gravity stabilized during each withdrawal period, the gravity and quality of the oil produced were indistinguishable from oil produced before the test.

Fig. 5

Produced oil gravity vs. time during Flowback Cycle 4. The top and bottom images (left to right) show the color change in produced fluids over time.

Fig. 5

Produced oil gravity vs. time during Flowback Cycle 4. The top and bottom images (left to right) show the color change in produced fluids over time.

Close modal

Produced Gas Composition

Gas composition during the injection and production cycle was monitored and can be seen in Fig. 6. The injection gas for the project is a stream of wet produced gas from nearby wells. The Btu content is slightly lower than that of the native gas from the subject well. It has been established that in some cyclic EOR processes, the light ends of the reservoir oil may vaporize into the high-pressure gas stream, leaving behind a heavier and heavier fraction of the oil with every cycle (Sie and Nguyen 2020c, 2022). Based on the sampling of the produced gas stream throughout this pilot test, there is no indication that the gas composition materially changed by being cycled into and out of the well as a working fluid. There is also no indication that any significant fraction of native gas was mixing into the injection gas. This would indicate that the oil being extracted was dead oil without a significant solution gas component.

Fig. 6

Gas compositions for the injection gas, two native gas samples taken before the beginning of cyclic gas injection, and four produced gas samples from the first (Cycles 1-1 and 1-2), third, and fourth withdrawal cycles.

Fig. 6

Gas compositions for the injection gas, two native gas samples taken before the beginning of cyclic gas injection, and four produced gas samples from the first (Cycles 1-1 and 1-2), third, and fourth withdrawal cycles.

Close modal

Oil Production Rate

Both the injection and production responses were stable (Fig. 7), and when compared on a time scale, the production profiles exhibited the characteristics of oil production from an injection gas-crude oil mixing zone. During the huff phase, an injection gas-crude oil mixing zone is generated in the shale matrix, and an induced fracture network referred to as a stimulated reservoir volume is formed. The injection gas components partition from the injection gas into the crude oil driven by thermodynamic phase behavior and continue to penetrate deeper into the stimulated reservoir volume. The peak oil production rate and production rate decline are strongly influenced by the oil content and the extent and volume of the mixing zone, which has been shown to increase with system pressure (Vallabhaneni et al. 2016; Li et al. 2017b; Sie and Nguyen 2022). This is consistent with the liquid production profiles for the five injection cycles (Fig. 7). The volume, rate, and decline of the oil production response appear to be proportional to the injection volumes, with the higher sustained peaks correlating with the larger later injection volumes. The days to peak oil production did continue to get longer for each subsequent cycle as the mixing zone continued to propagate deeper into the stimulated reservoir volume. Based on the core experiments that showed high first-cycle recovery, the delayed peak through cycles could indicate that each cycle is bringing oil incrementally farther from the wellbore. Further investigation could yield useful insights into the reservoir and fluid characteristics and allow quantitative conclusions to be drawn from future experiments. Some attempts have been made to identify oil production mechanisms from the fracture network and the matrix by inverse modeling of these data.

Fig. 7

Oil production rate vs. time for five gas injection cycles.

Fig. 7

Oil production rate vs. time for five gas injection cycles.

Close modal

Gas Production Rate

Initial gas production rate was limited by the ability of the surface facility to heat and process the gas. The solvent extraction process that characterizes huff ‘n’ puff in this implementation is dependent on pressure to maintain the desired mass transfer between the injection gas and the oil in the mixing zone. Thus, slower flowback rates were selected to continuously promote the far-field mass transfer and to reduce the effect of competing gas flow on the relative permeability of oil due to excessive pressure depletion. An example of the implementation of this strategy can be seen in Fig. 8, where the initial choke at the surface was adjusted to try and maintain a constant gas rate during the initial flowback period while the flowback was primarily gas. As the fluid fraction of the flowback increased, the gas rate began to fall precipitously due to the rapidly increasing hydrostatic pressure from liquids entering the tubing. Once this rapid fall was observed and confirmed by increasing fluid rates, the choke was held constant until tubing pressure declined to the point that gas lift was needed. This conservative approach seeks to minimize that rate of drawdown in the well while allowing liquids to be produced at reasonable rates.

Fig. 8

Oil and net gas production rates vs. time for Cycle 5. Produced oil gravity and precyclic gas injection oil production rate are also shown.

Fig. 8

Oil and net gas production rates vs. time for Cycle 5. Produced oil gravity and precyclic gas injection oil production rate are also shown.

Close modal

A primary economic driver is the recovery of the injected gas as it has an economic value that will be lost if the gas is not recovered. It is hard to be certain what the final recovery of any of the cycles would be because no production phase has yet been run until all the gas was recovered or until the well stopped producing gas. From Fig. 9, we can see that when the production phase ended for the next cycle to begin, there was still meaningful gas recovery taking place. If the recovery profiles were projected out to approximate a final gas recovery, it would be in the range of 70 to 80% of the injected volume.

Fig. 9

Gas production in terms of (a) the percentage of injection gas recovered and (b) net gas production rate in Mscf/D from each of the five gas injection cycles.

Fig. 9

Gas production in terms of (a) the percentage of injection gas recovered and (b) net gas production rate in Mscf/D from each of the five gas injection cycles.

Close modal

It should also be noted that during Cycle 5 of this project, a gas measurement error caused the underreporting of produced gas from approximately Day 150 until Day 275. The appropriate gas volumes would be a curve projected through the gap caused by the bad data and would positively affect both Figs. 8 and 9. For accuracy, no massaging of the data is done to account for this error in the data presented in this paper.

When comparing the gas flowback from cycle to cycle (Fig. 9b), the slope is a mirror image of the injection plots (Fig. 4), where a slope change can clearly be observed as flow transitions from near-wellbore to far-field or matrix-dominated flow regimes. The amount of gas recovered before liquid production started was also larger with each cycle. Whether this is a function of solely the larger volumes of injection gas each cycle placed or it is also a function of fluids having to migrate from deeper in the formation in each subsequent cycle is difficult to ascertain but is an important area of interest that deserves further study.

Cumulative Production

The overall effects of the implementation of the huff ‘n’ puff process on this well can be seen in Fig. 10. The cumulative baseline projection, obtained from a decline curve analysis using a reserve estimation software, is shown and makes clear that this trial produced incremental oil over the baseline curve. When trying to decide how long to flow a cycle, these data iarehelpful because they highlight the lost production during injection and how that is made up for during the production phase. The short initial cycles prove not to be overly effective at changing the trajectory of the cumulative production plot. It is because the production time is not very long compared to the injection time during which no oil is produced. The larger longer cycles prove to be the most effective at recovering incremental volumes over the baseline projection. Future cycles will use the profiles created here to optimize the size and timing of the cycles to maximize the recovery and net present value.

Fig. 10

Cumulative oil production from the huff ‘n’ puff pilot compared to the projected baseline cumulative oil production. Oil production rate from the pilot used to calculate the cumulative production is also shown.

Fig. 10

Cumulative oil production from the huff ‘n’ puff pilot compared to the projected baseline cumulative oil production. Oil production rate from the pilot used to calculate the cumulative production is also shown.

Close modal

These data represent the first well in the northern Eagle Ford extension to have cyclic gas injection implemented; this is one of the few single well examples of huff ‘n’ puff implementation. Preparations for the pilot were focused on several areas, including (1) gaining practical insights that could be used in real time to optimize the process, (2) creating high-frequency learning loops that incorporated the real-time results from the laboratory, and (3) conducting field experiments and driving design and implementation changes. This process of continuous incorporation of learnings allowed several step changes that yielded a result that displayed incrementally better results with each cycle of implementation.

Laboratory experiments were designed to provide understanding and insight into the nature and thresholds of the fluid-fluid interaction as well as the fluid-rock interactions. The MCP for the crude oilfield gas system was about 4,500 psia at the reservoir temperature. An oil recovery test on a shale sample indicates that more than 80% of the initial dead oil in place could be recovered by field gas injection at MCP.

The results of the field experiments answered some basic questions about the containment of injected gas within the zone and allowed for increasing volumes to be utilized. Both on injection and withdrawal, the system displayed remarkable consistency and repeatability in its behavior. The injection pressure response was shown to be predictable based on the volume injected and indicative of the flow regime that the injectant front was propagating through. Likewise, the withdrawal phases demonstrated that under similar flowback procedures, similar flow profiles can be achieved even with incrementally larger volumes of injected gas. The recovery of very high API oil rapidly abated and reverted to normal reservoir quality fluids as the pressures declined during each withdrawal cycle. A critical observation was that good quantities of gas were recovered each cycle and that the composition of that gas was not substantially changed during the injection and withdrawal processes.

The implementation of the huff ‘n’ puff process in this well was not only a technical success in that it produced incrementally more oil than was predicted based on its decline curve, but it was also an economic success from its early stages. Lessons learned in this project have already allowed for several other low-cost tests and allowed the reduction of capital implementation cost to be reduced from approximately USD 250,000 to less than approximately USD 100,000. One of the lessons that enabled this success was designing a system that is flexible, modular, and can be easily adapted as a project progresses.

The authors acknowledge MD America Energy LLC for the support of this study. Special thanks are expressed to Chao-yu Sie, a former PhD student at The University of Texas at Austin, for conducting the laboratory experiments.

Appendix A

Table A-1

Dead oil composition.

Methane 0.06% 
Ethane 0.36% 
Propane 2.34% 
Iso-butane 0.97% 
n-butane 4.04% 
Iso-pentane 3.04% 
n-pentane 4.28% 
i-hexanes 3.48% 
n-hexane 2.74% 
2,2,4-Trimethylpentane 0.07% 
Benzene 0.35% 
Heptanes 10.13% 
Toluene 1.09% 
Octanes 10.30% 
Ethylbenzene 0.31% 
Xylenes 1.57% 
Nonanes 7.05% 
Decanes Plus 46.25% 
Total 100.00% 
Methane 0.06% 
Ethane 0.36% 
Propane 2.34% 
Iso-butane 0.97% 
n-butane 4.04% 
Iso-pentane 3.04% 
n-pentane 4.28% 
i-hexanes 3.48% 
n-hexane 2.74% 
2,2,4-Trimethylpentane 0.07% 
Benzene 0.35% 
Heptanes 10.13% 
Toluene 1.09% 
Octanes 10.30% 
Ethylbenzene 0.31% 
Xylenes 1.57% 
Nonanes 7.05% 
Decanes Plus 46.25% 
Total 100.00% 

This paper (SPE 209429) was accepted for presentation at the SPE Improved Oil Recovery Conference, Virtual, 25–29 April 2022, and revised for publication. Original manuscript received for review 2 May 2022. Revised manuscript received for review 21 January 2023. Paper peer approved 23 January 2023.

Alfarge
,
D.
,
Wei
,
M.
, and
Bai
,
B
.
2017
.
IOR Methods in Unconventional Reservoirs of North America: Comprehensive Review
.
Paper presented at the
SPE Western Regional Meeting
,
Bakersfield, California, USA
, 23–27 April. SPE-185640-MS. 10.2118/185640-MS.
Alfarge
,
D.
,
Alsaba
,
M.
,
Wei
,
M
. et al. 
.
2018
.
Miscible Gases Based EOR in Unconventional Liquids Rich Reservoirs: What We Can Learn
.
Paper presented at the
SPE International Heavy Oil Conference and Exhibition
,
Kuwait City, Kuwait
, 10–12 December. SPE-193748-MS. 10.2118/193748-MS.
Anand
,
N.
,
Tang
,
B.
,
Nguyen
,
B
. et al. 
.
2019
.
Effects of Rate, Temperature, and Solvent Type on Vapor/Oil Gravity Drainage (VOGD) in Fractured Reservoirs
.
SPE J.
24
(
3
):
973
987
. SPE-190350-PA. 10.2118/190350-PA.
Arshad
,
A.
,
Al-Majed
,
A. A.
,
Menouar
,
H
. et al. 
.
2009
.
Carbon Dioxide (CO2) Miscible Flooding in Tight Oil Reservoirs: A Case Study
.
Paper presented at the
Kuwait International Petroleum Conference and Exhibition
,
Kuwait City, Kuwait
, 14–16 December. SPE-127616-MS. 10.2118/127616-MS.
Carlsen
,
M. L.
,
Dahouk
,
M. M.
,
Mydland
,
S
. et al. 
.
2020
a.
Compositional Tracking: Predicting Wellstream Compositions in Tight Unconventionals
.
Paper presented at the
International Petroleum Technology Conference
,
Dhahran, Kingdom of Saudi Arabia
, 13–15 January. IPTC-19596-MS. 10.2523/IPTC-19596-MS.
Carlsen
,
M. L.
,
Mydland
,
S.
,
Dahouk
,
M. M
. et al. 
.
2020
b.
Compositional Tracking of A Huff-n-Puff Project in the Eagle Ford: A Second Take
.
Paper presented at the
SPE/AAPG/SEG Unconventional Resources Technology Conference
,
Virtual
, 20–22 July. URTEC-2020-2869-MS. 10.15530/urtec-2020-2869.
Guo
,
P.
,
Wang
,
Z.
,
Shen
,
P
. et al. 
.
2009
.
Molecular Diffusion Coefficients of the Multicomponent Gas-Crude Oil Systems under High Temperature and Pressure
.
Ind Eng Chem Res
48
(
19
):
9023
9027
. 10.1021/ie801671u.
De Gennes
,
P.-G.
,
Brochard-Wyart
,
F.
, and
Quéré
,
D
.
2004
.
Capillarity and Wetting Phenomena
,
35
38
.
New York, USA
:
Springer
. 10.1007/978-0-387-21656-0.
Gozalpour
,
F.
,
Ren
,
S. R.
, and
Tohidi
,
B
.
2005
.
CO2 EOR and Storage in Oil Reservoir
.
Oil Gas Sci Technol
60
(
3
):
537
546
. 10.2516/ogst:2005036.
Hamdi
,
H.
,
Clarkson
,
C. R.
,
Ghanizadeh
,
A
. et al. 
.
2018
.
Huff-n-Puff Gas Injection Performance in Shale Reservoirs: A Case Study From Duvernay Shale in Alberta, Canada
.
Paper presented at the
SPE/AAPG/SEG Unconventional Resources Technology Conference
,
Houston, Texas, USA
, 23–25 July. URTEC-2018-2902835. 10.15530/urtec-2018-2902835.
Hamdi
,
H.
,
Clarkson
,
C. R.
,
Esmail
,
A
. et al. 
.
2021
.
Optimizing the Huff ‘n’ Puff Gas Injection Performance in Shale Reservoirs Considering the Uncertainty: A Duvernay Shale Example
.
SPE Res Eval & Eng
24
(
1
):
219
237
. SPE-195438-PA. 10.2118/195438-PA.
Hawthorne
,
S. B.
,
Gorecki
,
C. D.
,
Sorensen
,
J. A.
et al. 
.
2014
.
Hydrocarbon Mobilization Mechanisms Using CO2 in an Unconventional Oil Play
.
Energy Procedia
63
:
7717
7723
. 10.1016/j.egypro.2014.11.805.
Hawthorne
,
S. B.
,
Miller
,
D. J.
,
Grabanski
,
C. B
. et al. 
.
2017
.
Measured Crude Oil MMPs with Pure and Mixed CO2, Methane, and Ethane, and Their Relevance to Enhanced Oil Recovery from Middle Bakken and Bakken Shales
.
Paper presented at the
SPE Unconventional Resources Conference
,
Calgary, Alberta, Canada
, 15–16 February. SPE-185072-MS. 10.2118/185072-MS.
Elmer
,
W. G.
and
Elmer
,
D. C
.
2017
.
Two-Stage Compression With Elevated Cooler-Discharge Temperatures Improves Wellsite Gas Lift Operations
.
SPE Prod & Oper
32
(
2
):
159
167
. SPE-181773-PA. 10.2118/181773-PA.
Hawthorne
,
S. B.
,
Miller
,
D. J.
,
Jin
,
L
. et al. 
.
2016
.
Rapid and Simple Capillary-Rise/Vanishing Interfacial Tension Method To Determine Crude Oil Minimum Miscibility Pressure: Pure and Mixed CO2, Methane, and Ethane
.
Energy Fuels
30
(
8
):
6365
6372
. 10.1021/acs.energyfuels.6b01151.
Hoffman
,
B. T
.
2018
.
Huff-N-Puff Gas Injection Pilot Projects in the Eagle Ford
.
Paper presented at the
SPE Canada Unconventional Resources Conference
,
Calgary, Alberta, Canada
, 13–14 March. SPE-189816-MS. 10.2118/189816-MS.
Hoffman
,
B. T.
and
John
,
J. E
.
2016
.
Improved Oil Recovery IOR Pilot Projects in the Bakken Formation
.
Paper presented at the
SPE Low Perm Symposium
,
Denver, Colorado, USA
, 5–6 May. SPE-180270-MS. 10.2118/180270-MS.
Huang
,
X.
,
Gu
,
L.
,
Li
,
S.
et al. 
.
2022
.
Absolute Adsorption of Light Hydrocarbons on Organic-Rich Shale: An Efficient Determination Method
.
Fuel
308
:
121998
. 10.1016/j.fuel.2021.121998.
Jacobs
,
T
.
2019
.
Shale EOR Delivers, So Why Won’t the Sector Go Big?
J Pet Technol
71
(
5
):
37
41
. 10.2118/0519-0037-JPT.
Jia
,
B.
,
Tsau
,
J.-S.
, and
Barati
,
R
.
2019
.
A Review of the Current Progress of CO2 Injection EOR and Carbon Storage in Shale Oil Reservoirs
.
Fuel
236
:
404
427
. 10.1016/j.fuel.2018.08.103.
Kovscek
,
A. R.
,
Tang
,
G. Q.
, and
Vega
,
B
.
2008
.
Experimental Investigation of Oil Recovery From Siliceous Shale by CO2 Injection
.
Paper presented at the
SPE Annual Technical Conference and Exhibition
,
Denver, Colorado, USA
, 21–24 September. SPE-115679-MS. 10.2118/115679-MS.
Li
,
L.
,
Zhang
,
Y.
, and
Sheng
,
J. J
.
2017
b.
Effect of the Injection Pressure on Enhancing Oil Recovery in Shale Cores during the CO2 Huff-n-Puff Process When It Is above and below the Minimum Miscibility Pressure
.
Energy Fuels
31
(
4
):
3856
3867
. 10.1021/acs.energyfuels.7b00031.
Li
,
L.
,
Sheng
,
J. J.
, and
Xu
,
J
.
2017
a.
Gas Selection for Huff-n-Puff EOR in Shale Oil Reservoirs Based upon Experimental and Numerical Study
.
Paper presented at the
SPE Unconventional Resources Conference
,
Calgary, Alberta, Canada
, 15–16 February. SPE-185066-MS. 10.2118/185066-MS.
Liu
,
Y.
,
Rui
,
Z.
,
Yang
,
T
. et al. 
.
2022
.
Using Propanol as an Additive to CO2 for Improving CO2 Utilization and Storage in Oil Reservoirs
.
Appl Energy
311
:
118640
. 10.1016/j.apenergy.2022.118640.
Mahzari
,
P.
,
Mitchell
,
T. M.
,
Jones
,
A. P
. et al. 
.
2021
.
Novel Laboratory Investigation of Huff-n-Puff Gas Injection for Shale Oils under Realistic Reservoir Conditions
.
Fuel
284
(
15
):
118950
. 10.1016/j.fuel.2020.118950.
Mahzari
,
P.
,
Jones
,
A. P.
, and
Oelkers
,
E. H
.
2020
.
Impact of In-Situ Gas Liberation for Enhanced Oil Recovery and CO2 Storage in Liquid-Rich Shale Reservoirs
.
Energy Sources A: Recovery Util Environ Eff
1
(
1
):
1
21
. 10.1080/15567036.2020.1815907.
Mahzari
,
P.
,
Mitchell
,
T.
,
Jones
,
A
. et al. 
.
2019
.
A New Mechanism for Enhanced Oil Recovery by CO2 in Shale Oil Reservoirs
.
Paper presented at the
IOR 2019–20th European Symposium on Improved Oil Recovery
,
Pau, France
. 10.3997/2214-4609.201900100.
Rao
,
D
.
2001
.
Gas Injection EOR- A New Meaning in the New Millennium
.
J Can Pet Technol
40
(
2
):
11
18
. PETSOC-01-02-DAS. 10.2118/01-02-DAS.
Rao
,
D. N.
and
Lee
,
J. I
.
2003
.
Determination of Gas-Oil Miscibility Conditions by Interfacial Tension Measurements
.
J Colloid Interface Sci
262
(
2
):
474
482
. 10.1016/S0021-9797(03)00175-9.
Sheng
,
J. J
.
2015
.
Enhanced Oil Recovery in Shale Reservoirs by Gas Injection
.
J Nat Gas Sci Eng
22
:
252
259
. 10.1016/j.jngse.2014.12.002.
Sie
,
C.
and
Nguyen
,
Q. P
.
2020
c.
Field Gas Huff-n-Puff for Improving Oil Recovery from Eagle Ford Shale Reservoirs
.
Paper presented at the
SPE Improved Oil Recovery Conference
,
Virtual
, 31 August– 4 September. SPE-200471-MS. 10.2118/200471-MS.
Vallabhaneni
,
S.
,
Olszewski
,
T. D.
,
Pope
,
M. C
. et al. 
.
2016
.
Facies and Stratigraphic Interpretation of the Eaglebine Play in Central Texas
.
GCAGS J
5
:
25
46
.
Sie
,
C.
,
Nguyen
,
B.
,
Verlaan
,
M
. et al. 
.
2018
.
Viscous Oil Recovery and In Situ Deasphalting in Fractured Reservoirs: Part 1. The Effect of Solvent Injection Rate
.
Energy Fuels
32
(
1
):
360
372
. 10.1021/acs.energyfuels.7b03375.
Sie
,
C.
,
Nguyen
,
B.
,
Castellanos Diaz
,
O
. et al. 
.
2019
.
Viscous Oil Recovery and In-Situ Deasphalting in Fractured Reservoirs – Part 2: Effect of Solvent Type and Temperature
.
Fuel
247
:
294
301
. 10.1016/j.fuel.2019.03.044.
Sie
,
C.
and
Nguyen
,
Q. P
.
2020
b.
A Pore-Scale Study of Non-Aqueous Foam for Improving Hydrocarbon Miscible Flooding in Water-Sensitive Tight Oil Formations
.
J Pet Sci Eng
195
. 10.15530/urtec-2019-955.
Sie
,
C.
and
Nguyen
,
Q. P
.
2020
a.
A Non-Aqueous Foam Concept for Improving Hydrocarbon Miscible Flooding in Low Permeability Oil Formations
.
Fuel
288
:
119732
. 10.1016/j.fuel.2020.119732.
Song
,
C.
,
Clarkson
,
C. R.
,
Hamdi
,
H
. et al. 
.
2021
.
Comparison of CO2 and Lean Gas Cyclic Injection (“Huff-n-Puff”) in Artificially-Fractured Shale Core Samples
.
Paper presented at the
SPE/AAPG/SEG Unconventional Resources Technology Conference
,
Houston, Texas, USA
, 26–28 July. URTEC-2021-5309-MS. 10.15530/urtec-2021-5309.
Sohrabi
,
M.
,
Danesh
,
A.
,
Tehrani
,
D. H.
et al. 
.
2008
.
Microscopic Mechanisms of Oil Recovery By Near-Miscible Gas Injection
.
Transp Porous Med
72
(
3
):
351
367
. 10.1007/s11242-007-9154-z.
Tang
,
B.
,
Anand
,
N.
,
Nguyen
,
B
. et al. 
.
2020
.
Effect of Pressure on Vapor/Oil Gravity Drainage in Fractured Reservoirs
.
SPE J.
25
(
1
):
197
211
. SPE-190251-PA. 10.2118/190251-PA.
Tharanivasan
,
A. K.
,
Yang
,
C.
, and
Gu
,
Y
.
2006
.
Measurements of Molecular Diffusion Coefficients of Carbon Dioxide, Methane, and Propane in Heavy Oil under Reservoir Conditions
.
Energy Fuels
20
(
6
):
2509
2517
. 10.1021/ef060080d.
Thomas
,
F. B.
,
Piwowar
,
M.
,
Noroozi
,
M
. et al. 
.
2020
b.
Experimental Measurements of Montney and Duvernay Gas-Cycling Enhanced Oil Recovery GCEOR
.
Paper presented at the
SPE Canada Unconventional Resources Conference
,
Virtual
, 28 September–2 October. SPE-199994-MS. 10.2118/199994-MS.
Thomas
,
F. B.
,
Noroozi
,
M.
,
Piwowar
,
M
. et al. 
.
2020
a.
A Comparison of Primary Depletion and GCEOR in the Montney Formation: Volatile Oil and Rich Gas Condensate
.
Paper presented at the
SPE/AAPG/SEG Unconventional Resources Technology Conference
,
Virtual
, 20–22 July. URTEC-2020-3266. 10.15530/urtec-2020-3266.
Tovar
,
F. D.
,
Barrufet
,
M. A.
, and
Schechter
,
D. S
.
2018
a.
Gas Injection for EOR in Organic Rich Shale. Part I: Operational Philosophy
.
Paper presented at the
SPE Improved Oil Recovery Conference
,
Tulsa, Oklahoma, USA
, 14–18 April. SPE-190323-MS. 10.2118/190323-MS.
Tovar
,
F. D.
,
Barrufet
,
M. A.
, and
Schechter
,
D. S
.
2018
b.
Gas Injection for EOR in Organic Rich Shales. Part II: Mechanisms of Recovery
.
Paper presented at the
SPE/AAPG/SEG Unconventional Resources Technology Conference
,
Houston, Texas, USA
, 23–25 July. URTEC-2018-2903026. 10.15530/urtec-2018-2903026.
Wang
,
X.
and
Gu
,
Y
.
2011
.
Oil Recovery and Permeability Reduction of a Tight Sandstone Reservoir in Immiscible and Miscible CO2 Flooding Processes
.
Ind Eng Chem Res
50
(
4
):
2388
2399
. 10.1021/ie1016046.
Tran
,
S.
,
Yassin
,
M. R.
,
Eghbali
,
S
. et al. 
.
2021
.
Quantifying Oil-Recovery Mechanisms during Natural-Gas Huff ‘n’ Puff Experiments on Ultratight Core Plugs
.
SPE J
26
(
1
):
498
514
. 10.2118/200341-PA.
US Energy Information Administration
.
2016
.
Annual Energy Outlook 2016 with Projections to 2040
.
DOE report DOE/EIA-0383
.
Sie
,
C.
and
Nguyen
,
Q. P
.
2022
.
Laboratory Investigations on Field Gas Huff-n-Puff for Improving Oil Recovery in Eagle Ford Shale-Effect of Operating Conditions
.
Energy Fuels
36
(
1
):
195
209
. 10.1021/acs.energyfuels.1c03003.