Smart Profile Control by Salt-Reversible Flocculation of Cationic

Oct 22, 2013 - In the profile modification test, this flocculation system diverts nearly 100% chase water into the low permeable tube compared to 20% ...
3 downloads 0 Views 2MB Size
Article pubs.acs.org/EF

Smart Profile Control by Salt-Reversible Flocculation of Cationic Microgels and Polyacrylamide Guanghui Li,†,* Guicai Zhang,† and Lei Wang‡ †

School of Petroleum Engineering, China University of Petroleum, Qingdao, Shandong 266555, People’s Republic of China Department of Petroleum Engineering, Colorado School of Mines, Golden, Colorado 80401, United States



ABSTRACT: A flocculation system is detailed to deal with the size matching problem between pore throats and microgels for profile modification in high permeable formation. This system consists of partially hydrolyzed polyacrylamide (HPAM) and cationic microgels synthesized by inverse suspension polymerization with acrylamide and methylacryloxylethyl trimethyl ammonium chloride as monomers. By adding HPAM into cationic microgel suspensions, flocculation occurs because of electrostatic attraction. Then it reversibly restructures from bulk to an insulated microgel particle at a critical salt concentration (CSC). In the profile modification test, this flocculation system diverts nearly 100% chase water into the low permeable tube compared to 20% by using microgels only. The performance improvement testifies the mechanism that HPAM flocculates around the cationic microgels to make up the interstitials between microgels and irregular pore throats. This flocculation plug can be removed by water with salinity above CSC due to the desorption of HPAM on cationic microgels.

1. INTRODUCTION Due to the washing out of excessive injection, most matured water-flooding reserviors are sufferring large permeability variance, yielding water channels of high permeability, and leaving significant oil unswept. Therefore, facile profile modification treatment for high permeable formation is more urgently needed than ever before. Considerable attention has been paid to gel technology with varied formulations for different reservoir conditions, usually including bulk gel,1 colloidal dispersion gel,2,3 and preperformed gel.4,5 However, most of these methods have been proved only marginally effective when dealing with water channeling in high permeable zone.6−8 In 1996, the Industry Consortium (BP, Chevron, Texaco, and Nalco Company) developed a novel “Bright Water” system, which was a microgel emulsion with particle diameters around 0.1−3 μm.9 And this system has been practiced for more than 85 treatments since the first application in Minas Field, Indonesia. But all of the successful applications published are restricted to the formation with permeability lower than 1500 × 10−3 μm2 (mD).10−13 For more permeable formation with large permeability range, size mismatching between microgels and pore throats probably discounts its performance. To deal with this mismatching problem, Yao et al.14 defined a matching factor, which is the ratio of particle diameter to average pore-throat diameter. They found that great profile-control performance could be achieved at the narrowly optimal matching factor 1.35−1.55; while beyond this range, the mismatching occurred and resulted into poor effect in profile control. This mismatching problem can be schematized by the large permeability of interstitial between irregular pore throats and trapped microgels, which is illustrated in Figure 1(a). To solve the above-mentioned problem, partially hydrolyzed polyacrylamide (HPAM) is carefully considered due to its flocculation effect with cationic particle. The flocculation phenomenon between microgels and opposite polyelectrolyte © 2013 American Chemical Society

Figure 1. The schematic diagram of plugging effect by microgels in water channel.

has been extensively studied in water treatment, papermaking, and biomedical engineering.15−18 And the flocs can restructure and change phase behaviors with component properties17,19 and environmental conditions.18,20 This restructuring feature gives an inspiration for solving the problem in an environmentally-sensitive way. In this work, microgels with a series of cationic degrees (CD, defined by mole fraction of cationic monomer in all monomers) were synthesized by inverse suspension polymerization with acrylamide (AM) and methylacryloxylethyl trimethyl ammonium chloride (DMC) as monomers. These particles can first narrow pore throats by retention on the pore walls. And then HPAM injected subsequently flocculates around these particles, making up the interstitials to solve the mismatching problem, as shown in Figure 1(b). The properties of these particles and the plugging mechanisms are validated in the following sections. Received: September 18, 2013 Revised: October 22, 2013 Published: October 22, 2013 6632

dx.doi.org/10.1021/ef401881q | Energy Fuels 2013, 27, 6632−6636

Energy & Fuels

Article

silicon sand. The diameter of glass bead lies within 170−430 μm with an average value of 340 μm. The corresponding mean diameter of the pore throat is calculated to be 52.63 μm on the premise of hexagonal close packing, in which the pore-throat diameter is 1/6.46 times of the average bead diameter.23 The tube packing method was described by Ding, et al.24 The test procedure is as follows: (1) inject water; (2) inject 0.5 pore volume (PV) microgel suspension with mass concentration of 0.4%; (3) inject 0.5 PV HPAM with concentration of 2000 mg/L; and (4) switch the flow direction of the low permeable tube to avoid blocking of the sand filter net at the end of tube after the previous two slugs, then continue waterflood until flow pressure is steady. A control experiment was done with an additional procedure after the four slugs. This procedure was that the highly permeable tube was flooded with 2%NaCl solution for 5PV. The tube permeabilities in the control experiment were 5323 and 382 mD, respectively. In all processes, pressure was recorded by pressure transducer and outflow was weighted to determine the flow fraction of these two tubes.

2. MATERIALS AND METHODS 2.1. Microgel Synthesis. Acrylamide (analytically pure) was recrystallized by chloroform and 78 wt % water solution of methylacryloxylethyl trimethyl ammonium chloride (DMC) was extracted by benzene. Methylene bisacrylamide (MBA) was used as the cross-linker. Span80 (Sorbitan monooleate, chemically pure) and 2,2′-Azobis [2-(2-imidazolin-2-yl)propane] dihydrochloride(Va-044) were used as dispersant and initiator, respectively; Industrial white oil, with a boiling range between 150 and 300 °C, was chosen as continuous phase. Partially hydrolyzed polyacrylamide (HPAM) was determined with the intrinsic viscosity of 2.087 L/g in 1 mol/L NaCl solution at 30 °C and hydrolysis degree about 25.2%. The inverse suspension polymerization was carried out as the following procedures. The oil phase was prepared by 0.60 g Span80 and 59.40 g white oil. The aqueous phase consisted of 50% comonomer (AM/DMC), 1% MBA, 0.05% Va-044, and 48.95% water (based on the total weight of water phase). Then, oil phase was put in four neck round-bottomed flask and stirred at 450 rpm in 50 °C water bath. Nitrogen purging was carried out for 15 min before 30 g aqueous phase was added by dropwise using a syringe pump at a rate of 1.0 mL/min. The reaction mixture was maintained for another 2 h after the addition process. Samples were withdrawn at the end of the reaction using a hypodermic syringe to measure the particle size distribution. In order to avoid irreproducibility, the sample point was labeled on the reactor and fixed at 3 cm away from the impeller shaft. At the end of the reaction, stirring was stopped and the two phases, white oil and copolymer particles, were obtained from top and bottom, respectively. Then these particles were rinsed by acetone and dried under a primary vacuum at 50 °C for 24 h. The final products were white powder and used for other experiments. 2.2. Morphology and Particle Size. Morphology and particle size were analyzed by observing under microscope. The microscope accuracy is 0.4 μm, and amplification factors are 40, 100, 400, and 1000. At least 400 particles were numerated to calculate the mean diameter and particle size distribution. 2.3. Zeta Potential. The products were diluted by freshly deionized water to 0.1% in mass concentration. The particle mobility was recorded by visual microscope of electrophoresis meter at 20 °C. The microsphere products with cationic degree of 5%, 10%, 15%, 20%, 25%, 30%, and 40% were all determined. For each product 10 runs were taken. Then the zeta potential was further calculated by Smoluchowski formula which was built in the attached software of the meter. 2.4. Flocculation and Adsorption. The HPAM solution with concentration of 10 000 mg/L was added by dropwise into 50 mL 0.1 wt % microgels suspension, and flocculation formed immediately due to the electrostatic interaction. This flocculation can disappear due to the addition of sufficient salt and violent oscillation, so the critical salt concentration (CSC) for this phase change was evaluated for microgels with varied CDs. Additionally, the adsorption amount of HPAM on the particles with 25% CD was determined as a function of HPAM concentration at 20 °C by a semiquantitative method. The principle is that the viscosity of HPAM is proportional to its concentration in dilute concentration regime.21,22 This method was performed as the following procedures. The above flocculation system was oscillated for 2 h in thermostat, which allowed the sufficient reaction. Then it was centrifugated at 4000 rpm for 30 min to get supernatant. The viscosity of supernatant was measured by Brookfield Viscometer LVDVII plus pro with ULA spindle at 20 °C. According to the viscosity-concentration standard curve of HPAM, the adsorption capacity was obtained regardless of possible preferential adsorption for HPAM with different molecular weight. 2.5. Profile-Control Test. The objective of the combined HPAM and cationic microgels is to selectively block highly permeable channels and divert subsequent water to sweep low permeability zone. So this test was implemented by paralleled sand-pack tubes of 30 cm in length and 2.5 cm in diameter at 20 °C. These two tubes were packed to get different permeabilities, the higher one of 5092 mD with 40−60 mesh glass beads and the lower one of 301.3mD with 80−120 mesh

3. RESULTS AND DISCUSSION 3.1. Morphology and Particle Size. A series of microgels with cationic degrees of 5%, 10%, 15%, 20%, 25%, 30%, and 40% was synthesized. The particle morphology of 25% CD microgels is shown as a typical example in Figure 2. It can be

Figure 2. The morphology photo of microgels with 25% cationic degree. The scale bar is 200 μm.

seen that most microgel particles are closely spherical and have a wide distribution. The 25% CD microgels were used for the following profile control test. As shown in Figure 3, the swollen size distribution ranges from several to more than 100 μm, and conforms to logarithm normal distribution. It can be observed

Figure 3. Size distribution of swollen microgels with 25% cationic degree. 6633

dx.doi.org/10.1021/ef401881q | Energy Fuels 2013, 27, 6632−6636

Energy & Fuels

Article

restructuring of flocs, it opens a way using this flocculating system to control profile with formation salinity lower than CSC and remove the flocculation by injecting higher salinity water. Further, this structure evolution also leads to capacity change of HPAM adsorption on cationic microgels, as shown in Figure 6. The presented curve indicates that the adsorption process of

as well that the mean diameter is 63 μm with standard deviation about 0.518. 3.2. Flocculation and Adsorption. Due to the electrostatic attraction between cationic microgels and HPAM, microgel suspension flocculated rapidly when it was dripped by HPAM solution. These flocculation structures varied at different salt concentrations. Taking 25% CD microgels as an example, the floc structures transform from bulk [in Figure 4(inset b)] to insulated particles [in Figure 4(inset a)] at the

Figure 6. The relationship of HPAM adsorption capacity to its concentration at 20 °C.

Figure 4. The evolution of flocculation morphology for microgels with varied cationic degrees as a function of salt concentration. Insets (a) and (b) are photos of the flocs with salt concentration higher and lower than restructuring salt concentration, respectively. The scale bar is 200 μm in these photos.

HPAM on microgels can be represented by the adsorption curve of Langmuir type. In this process, the HPAM fast adsorbs on microgel surface, initially with a flat configuration. As the polymer concentration increases, the subsequently adsorbed polymer molecules may “crowd” the flat adsorbed polymer molecules into the spacial configuration. Considerable rearrangement of the polymer segments on the surface is likely to take place during this crowding stage and maximum adsorption may be reached slowly.26 According to the adsorption curve, the maximum adsorbing capacity is about 409 mg/g in 0 wt % NaCl solution and lowers to 237 mg/g in 2 wt % NaCl solution, which indicates the desorption feature of HPAM under high salinity. It is reasonable that the desorbed HPAM can be easily carried away within removing flocculation by higher salinity water. According to this adsorption curve, the HPAM usage for sufficiently flocculating 0.4 wt % microgels is designed to be 2000 mg/L in next profile control test. 3.3. The Profile-Control Test. Figure 7 shows the pressure response (bottom) and flow fractions (top) in the profile control test. From Figure 7, the initial injection water goes along the high permeability tube with fairly low pressure

critical salt concentration (CSC). When diluted by distilled water, the isolated floc reverses back to bulk state. A similar phenomenon is found for microgels with different CDs, and the CSC increases with the CD, as indicated by the line in Figure 4. This structure evolution is consistent with rearrangement process of adsorbed polyelectrolyte on the polystyrene particles investigated by Adachi, et al.15 And it can be explained that when salt concentration is lower than CSC, the polyelectrolyte coil has a spatially adsorbing conformation with tails and loops extending outside the electrostatic double layer of microgels, so a bridging effect occurs; and beyond this concentration, polyelectrolyte coils tend to be a flat adsorbing conformation, resulting in the insulated floc structure.17,19 Besides, in Figure 5, the zeta potential increases to the maximum at 30% CD, corresponding to the tendency of CSC.

Figure 5. The zeta potential changes as a function of cationic degree.

The reason is that the higher Zeta potential, the denser binding sites between microgel surface and polyelectrolyte, rendering that HPAM prefers being adsorbed in the patch-like conformation,25 so higher salinity is required to compress polymer coils into a flat geometry. In view of reversible

Figure 7. The pressure response (bottom) and flow fraction (top) change with injection volume. 6634

dx.doi.org/10.1021/ef401881q | Energy Fuels 2013, 27, 6632−6636

Energy & Fuels

Article

permeable tube smoothly during injection. While injecting 2% NaCl solution, a wave-type pressure is observed before the injection volume reaches 1 PV, and afterward it gradually decreases to 0.074 MPa, which is very close to the stable injection pressure of 0.065 MPa of microgels. This wave-type pressure indicates the advance of flocs deep into the tube. Similar pressure curves are found on preperformed gel studied by Bai et al27 and microgels with optimal matching factor studied by Yao et al.14 The continuous decrease in pressure indicates that the flocculation becomes movable because of HPAM desorption in water of high salinity. As is stated above, the flocculation produced by HPAM and cationic microgels has potential to modify profile for the high permeability formation, and avoid the mismatching problem between particles and pore throats in using microgels only.

response of 0.004 MPa; then, microgels are injected and get ten times of initial pressure response, but the flow fractions are stable at 80% for high permeability tube and 20% for low permeability tube. This poor performance indicates an insufficient plugging due to the mismatching as shown in Figure 1(a), even though the mean diameter of microgels is 1.2 times larger than the mean pore throat of glass beads. In the HPAM slug, the pressure gradually rises up to 0.45 MPa, possibly attributing to the flocculation reaction, but the flow fraction of low permeable tube reversely goes down due to the entrance block of the low permeable tube by this reaction. In the final slug, water breaks through these flocculation block formed in the previous two slugs when the pressure goes up to 0.3 MPa, and then the flow is nearly completely diverted to the low permeable tube, demonstrating great performance of profile modification. After the test, three samples were taken from the location 1, 2, 3 of the highly permeable tube in Figure 8(a). The

4. CONCLUSIONS To deal with the mismatching problem between microgels and pore throats in highly permeable formation, a smart flocculation system sensitive to salinity was successfully designed. This system consisted of HPAM and cationic microgels synthesized by inverse suspension polymerization. The floc structures formed by electrostatic attraction between these microgels and HPAM reversibly change from bulk to an insulated microgel at critical salt concentration (CSC). The adsorption curve of HPAM on cationic microgels conforms to the Langmuir-type adsorption law and its adsorption capacity on 25% CD microgels decreases from 409 mg/g at 0% NaCl to 237 mg/g at 2% NaCl. This flocculation system shows far better performance in a paralleled tube test than that of microgels only, which proves the assumed mechanism that HPAM flocculates around cationic microgels to plug the interstitials between microgels and pore throats. This flocculation plug can be removed by water with salinity above the CSC.

Figure 8. The schematic setup of paralleled sand pack and adsorption morphology corresponding to the sample points. Part (a) is the schematic of paralleled sand pack; parts (b), (c), and (d) are photos of samples taken in high permeable tube from points 1, 2, and 3, respectively. The scale bar is 200 μm in the photos.



morphologies of these samples are presented in Figure 8(b−d). It can be seen in the Figure 8(b),(c) that few cationic microgels are detained, which rules out possible performance contribution from the redistribution of microgels in depth; while a mount of flocculation is observed, as shown in the red circle of Figure 8(d), consolidating that the injected HPAM flocculates around the cationic microgels and the interstitials are plugged as the assumed mechanism in Figure 1(b). Figure 9 shows the pressure contrast of two steps: microgel injection and 2% NaCl injection in the control experiment. As shown in Figure 9, the microgels transport in the highly

AUTHOR INFORMATION

Corresponding Author

*Tel: +86 532 86981178; e-mail: [email protected]. Notes

The authors declare no competing financial interest.



ACKNOWLEDGMENTS The authors appreciate the great help from the Enhanced Oil Recovery Center in Oilfield Chemistry Department, China University of Petroleum.



REFERENCES

(1) Dai, C. l.; You, Q.; Yuhong, X.; Long, H.; Ya, C.; Fulin, Z. Case Study on Polymer Gel to Control Water Coning for Horizontal Well in Offshore Oilfield. Presented at Offshore Technology Conference, Houston, Texas, U.S.A., 2−5 May, 2011; Paper SPE 21125. (2) Coste, J.-P.; Liu, Y.; Bai, B.; Li, Y.; Shen, P.; Wang, Z.; Zhu, G. InDepth Fluid Diversion by Pre-Gelled Particles. Laboratory Study and Pilot Testing, Presented at SPE/DOE Improved Oil Recovery Symposium, Tulsa, Oklahoma, USA, 3−5 April, 2000. Paper SPE 59362. (3) Chang, H.; Xingguang, S.; Long, X.; Heng, L.; Zhidong, G.; Yuming, Y.; Yuguo, X.; Gang, C.; Kaoping, S.; Mack, J. Successful Field Pilot of in-Depth Colloidal Dispersion Gel (CDG) Technology in Daqing Oil Field. Presented at SPE/DOE Symposium on Improved Oil Recovery, Tulsa, Oklahoma, U.S.A, 17−21 April, 2004; Paper SPE 89460.

Figure 9. The pressure contrast of two steps: microgel injection and 2% NaCl injection in the control experiment. 6635

dx.doi.org/10.1021/ef401881q | Energy Fuels 2013, 27, 6632−6636

Energy & Fuels

Article

(4) Bai, B.; Li, L.; Liu, Y.; Liu, H.; Wang, Z.; You, C. Preformed Particle Gel for Conformance Control: Factors Affecting Its Properties and Applications. SPE Reserv. Eval. Eng. 2007, 10 (4), 415−422. (5) Bai, B.; Huang, F.; Liu, Y.; Seright, R.; Wang, Y. Case Study on Prefromed Particle Gel for In-Depth Fluid Diversion. Presented at SPE/DOE Symposium on Improved Oil Recovery Tulsa, Oklahoma, U.S.A., 20−23 April, 2008; Paper SPE 113997. (6) Smith, J. Performance of 18 Polymers in Aluminum Citrate Colloidal Dispersion Gels. Presented at SPE International Symposium on Oilfield Chemistry, San Antonio, Texas, U.S.A., 14−17 February, 1995; Paper SPE 28989. (7) Chauveteau, G.; Tabary, R.; Renard, M.; Omari, A. Controlling in Situ Gelation of Polyacrylamides by Zirconium for Water Shutoff. Presented at SPE International Symposium on Oilfield Chemistry, Houston, Texas, U.S.A., 16−19 February, 1999; Paper SPE 50752. (8) Bai, B.; Zhang, H. Preformed-Particle-Gel Transport through Open Fractures and Its Effect on Water Flow. SPE J. 2011, 16 (2), 388−400. (9) Frampton, H.; Morgan, J.; Cheung, S.; Munson, L.; Chang, K.; Williams, D. Development of a Novel Waterflood Conformance Control System, Presented at SPE/DOE Symposium on Improved Oil Recovery, Tulsa, Oklahoma, U.S.A., 17−21 April, 2004; Paper SPE 89391. (10) Ghaddab, F.; Kaddour, K.; Tesconi, M.; Brancolini, A.; Carniani, C.; Galli, G. EI Borma-In-Depth Profile Modification: A Tertiary Method for Enhanced Oil Recovery for a Mature Field. Presented at SPE Production and Operations Conference and Exhibition, Tunis, Tunisia, 8−10 June, 2010; Paper SPE 136140. (11) Ohms, D.; McLeod, J.; Graff, C.; Frampton, H.; Morgan, J.; Cheung, S.; Chang, K.-T. Incremental-Oil Success From Waterflood Sweep Improvement in Alaska. SPE Prod. Oper. 2010, 25 (3), 247− 254. (12) Husband, M.; Ohms, D.; Frampton, H.; Carhart, S.; Carlson, B.; Chang, K.-T.; Morgan, J. Results of a Three-Well Waterflood Sweep Improvement Trial in the Prudhoe Bay Field Using a Thermally Activated Particle System. Presented at SPE Improved Oil Recovery Symposium, Tulsa, Oklahoma, U.S.A., 24−28 April, 2010; Paper SPE 129967. (13) Yañez, P. A. P.; Mustoni, J. L.; Relling, M. F.; Chang, K.-T.; Hopkinson, P.; Frampton, H. New Attempt in Improving Sweep Efficiency at the Mature Koluel Kaike and Piedra Clavada Waterflooding Projects of the S. Jorge Basin in Argentina. Presented at SPE Latin American and Carribbean Petroleum Engineering Conference, Buenos Aires, Argentina, 15−18 April, 2007; Paper SPE 107923. (14) Yao, C.; Lei, G.; Li, L.; Gao, X. Selectivity of Pore-scale Elastic Microspheres as a Novel Profile Control and Oil Displacement Agent. Energy Fuels 2012, 26, 5092−5101. (15) Adachi, Y.; Aoki, K. Restructuring of Small Flocs of Polystyrene Latex with Polyelectrolyte. Colloids Surf., A 2009, 342 (1−3), 24−29. (16) Barany, S.; Nagy, M.; Skvarla, J. Electrokinetic Potential of Polystyrene Particles in Polyelectrolyte and Polyelectrolyte Mixtures Solutions. Colloids Surf., A 2012, 413, 200−207. (17) Solberg, D.; Wågberg, L. Adsorption and Flocculation Behavior of Cationic Polyacrylamide and Colloidal Silica. Colloids Surf., A 2003, 219 (1), 161−172. (18) Zeng, Z.; Patel, J.; Lee, S.-H.; McCallum, M.; Tyagi, A.; Yan, M.; Shea, K. J. Synthetic Polymer Nanoparticle-Polysaccharide Interactions: A Systematic Study. J. Am. Chem. Soc. 2012, 134 (5), 2681− 2690. (19) Eriksson, L.; Alm, B.; Stenius, P. Formation and Structure of Polystyrene Latex Aggregates Obtained by Flocculation with Cationic Polyelectrolytes: 1. Adsorption and Optimum Flocculation Concentrations. Colloids Surf., A 1993, 70 (1), 47−60. (20) Hierrezuelo, J.; Vaccaro, A.; Borkovec, M. Stability of Negatively Charged Latex Particles in the Presence of a Strong Cationic Polyelectrolyte at Elevated Ionic Strengths. J. Colloid Interface Sci. 2010, 347 (2), 202−208.

(21) Li, H.; Hou, W.; Zhang, Y. Rheological Properties of Aqueous Solution of New Exopolysaccharide Secreted by a Deep-Sea Mesophilic Bacterium. Carbohydr. Polym. 2011, 84 (3), 1117−1125. (22) Papanagopoulos, D.; Pierri, E.; Dondos, A. Influence of the Shear Rate, Of the Molecular Architecture and of the Molecular Mass on the Critical Overlapping Concentration c*. Polymer 1998, 39 (11), 2195−2199. (23) Greenwood, R.; Luckham, P.; Gregory, T. The Effect of Diameter Ratio and Volume Ratio on the Viscosity of Bimodal Suspensions of Polymer Latices. J. Colloid Interface Sci. 1997, 191 (1), 11−21. (24) Ding, B.; Zhang, G.; Ge, J.; Liu, X. Research on Mechanisms of Alkaline Flooding for Heavy Oil. Energy Fuels 2010, 24 (12), 6346− 6352. (25) Sadeghpour, A.; Seyrek, E.; Szil gyi, I.; Hierrezuelo, J.; Borkovec, M. Influence of the Ionization Degree and Molecular Mass of Weak Polyelectrolytes on Charging and Stability Behavior of Oppositely Charged Colloidal Particles. Langmuir 2011, 27 (15), 9270−9276. (26) Peterson, C.; Kwei, T. The Kinetics of Polymer Adsorption onto Solid Surfaces. J. Phys. Chem. 1961, 65 (8), 1330−1333. (27) Bai, B.; Liu, Y.; Coste, J.-P.; Li, L. Preformed Particle Gel for Conformance Control: Transport Mechanism through Porous Media. Presented at SPE/DOE Symposium on Improved Oil Recovery, Tulsa, Oklahoma, U.S.A., 17−21 April, 2004; Paper SPE 89468.

6636

dx.doi.org/10.1021/ef401881q | Energy Fuels 2013, 27, 6632−6636