Jump to ContentJump to Main Navigation
Show Summary Details
More options …

Open Chemistry

formerly Central European Journal of Chemistry


IMPACT FACTOR 2018: 1.512
5-year IMPACT FACTOR: 1.599

CiteScore 2018: 1.58

SCImago Journal Rank (SJR) 2018: 0.345
Source Normalized Impact per Paper (SNIP) 2018: 0.684

ICV 2017: 165.27

Open Access
Online
ISSN
2391-5420
See all formats and pricing
More options …
Volume 16, Issue 1

Issues

Volume 13 (2015)

Performance and thermal decomposition analysis of foaming agent NPL-10 for use in heavy oil recovery by steam injection

Fa-Jun Zhao
  • Corresponding author
  • Key Laboratory of Oil Recovery Enhance of Ministry of Education, Northeast Petroleum University, Daqing City, Heilongjiang province, China, 163318
  • Email
  • Other articles by this author:
  • De Gruyter OnlineGoogle Scholar
/ Yun-Long Wang
  • College of Chemistry and Chemical Engineering, Northeast Petroleum University, Daqing, 163318, Heilongjiang, China
  • Other articles by this author:
  • De Gruyter OnlineGoogle Scholar
/ Jun Song
  • College of Chemistry and Chemical Engineering, Northeast Petroleum University, Daqing, 163318, Heilongjiang, China
  • Other articles by this author:
  • De Gruyter OnlineGoogle Scholar
/ Hai-Cheng Ma
  • College of Chemistry and Chemical Engineering, Northeast Petroleum University, Daqing, 163318, Heilongjiang, China
  • Other articles by this author:
  • De Gruyter OnlineGoogle Scholar
/ Hao-Liang Liu
  • College of Chemistry and Chemical Engineering, Northeast Petroleum University, Daqing, 163318, Heilongjiang, China
  • Other articles by this author:
  • De Gruyter OnlineGoogle Scholar
Published Online: 2018-02-13 | DOI: https://doi.org/10.1515/chem-2018-0002

Abstract

Foaming agents, despite holding potential in steam injection technology for heavy oil recovery, are still poorly investigated. In this work, we analyzed the performance of the foaming agent NPL-10 in terms of foam height and half-life under various conditions of temperature, pH, salinity, and oil content by orthogonal experiments. The best conditions of use for NPL-10 among those tested are T=220°C, pH 7, salinity 10000 mg·L–1 and oil content 10 g·L–1. Thermal decomposition of NPL-10 was also studied by thermogravimetric and differential thermal analyses. NPL-10 decomposes above 220°C, and decomposition is a two-step process. The kinetic triplet (activation energy, kinetic function and pre-exponential factor) and the corresponding rate law were calculated for each step. Steps 1 and 2 follow kinetics of different order (n = 2 and ½, respectively). These findings provide some criteria for the selection of foaming agents for oil recovery by steam injection.

Keywords: thermal recovery; foam flooding; foaming agent; foamability; thermal decomposition

1 Introduction

The extraction of heavy crude oil from reservoirs is usually carried out by steam injection. There are several forms of the technology, with the major ones being steam stimulation, steam flooding and steam-assisted gravity drainage [1,2]. The main challenge in the use of steam for heavy oil recovery is steam channeling, i.e., the ability to make steam flow in a single direction (from the injector to the extractor) in such a way as to maximize the extraction yield. Another problem posed by steam injection is that steam can seriously damage the underground well structure. This phenomenon can endanger the lives of workers while resulting into thermal dispersion and poor mining effect [3, 4, 5]. Heavy oil reservoirs in China are quite heterogeneous in terms of layer permeability, thus requiring highly different conditions for optimal steam channeling [4, 5]. Foams can improve steam channeling and enhance oil recovery by reducing the surface tension and viscosity of fluids [6, 7, 8, 9]. However, the use of foaming agents in heavy oil recovery is still limited due to the small number of studies on the behavior of such agents at high temperature and salinity, which prevents the selection of the most appropriate agent.

In this work, we analyzed the performance of a high temperature foaming agent (NPL-10) in terms of foam height and half-life in order to establish the optimal conditions of its use. The thermal decomposition process of NPL-10 was analyzed by thermogravimetric and differential thermal analyses, and the kinetics of non-isothermal decomposition was studied to quantify the effect of temperature on decomposition rate.

2 Experimental

2.1 Chemicals

Test solutions for orthogonal experiments were prepared using deionized water. Three solutions with different total concentrations (5000, 10000 and 15000 mg·L–1) were used. Solutions had the following composition: 5000 mg·L–1 solution:NaHCO33000mg·L–1, NaCl700mg·L–1;10000mg·L–1 solution: NaHCO3 6000 mg·L–1, NaCl 1400 mg·L–1; 15000 mg·L–1 solution: NaHCO3 9000 mg·L–1, NaCl 2100 mg·L–1. Concentrations of other ions were simulated by a water preparation and salinity calculation software. Heavy oil was collected from Du-66 block area of Liaohe Oilfield and dehydrated before use. Foaming agent NPL-10 (long-chain alpha olefin sulfonate) was prepared according to literature [10, 11], and had purity > 97%.

2.2 Methods

Reactions were performed in a high temperature, high pressure reaction kettle WYF-1 supplied by Haian Jiangsu oil factory. Foamability and foam stability were determined by a 2151 Roche foam instrument (standard form) using a customized GB/T 7462-94 Ross-Miles method. Thermogravimetric and differential thermal analyses were carried out by a thermal analyzer DuPont 2100 (Perkin Elmer, USA) under nitrogen.

Ethical approval

This research is not related to human or animal use.

3 Results and discussion

3.1 Performance analysis of NPL-10

We studied the effect of four factors (temperature, salinity, pH, and oil content) on the foam height and half-life of NPL-10 by orthogonal experiments [12, 13, 14]. The height and half-life of NPL-10 foam were measured under various combinations of T (220, 260, and 300°C), pH (5, 7, and 9), salinity (5000, 10000, and 15000 mg·L–1) and oil content (0, 5, and 10 g·L–1; Table 1). The sensitivity of foam height and half-life to each factor was examined by a trend chart (Figure 1). Results showed that foam height and half-life depend on temperature (A), pH (B), salinity (C) and oil content (D) to a different extent (Figure 1). At 220°C, the order of factors observed for foam height is A>B>C>D; in other words, foam height mainly depends on temperature. Foam height decreases sharply on increasing temperature, and the most marked decrease is observed when temperature rises from 260°C to 300°C. Conversely, an increase in the other factors (pH, salinity and oil content) has almost no effect on foam height.

NPL-10 trend chart of factors at three different levels. Level 1 = 220°C, level 2 = 260°C; level 3 = 300°C. Top lines (dark blue, light blue, black and pink): effect on foam height. Bottom lines (violet, red, brown, green): effect on foam half-life.
Figure 1

NPL-10 trend chart of factors at three different levels. Level 1 = 220°C, level 2 = 260°C; level 3 = 300°C. Top lines (dark blue, light blue, black and pink): effect on foam height. Bottom lines (violet, red, brown, green): effect on foam half-life.

Table 1

Orthogonal experimental data for foam height and half-life.

The order of factors observed for foam half-life at 220°C is A>D>C>B. The most important factor for foam half-life is still temperature, but the second one is the oil content, not pH. Analogously to foam height, foam halflife decreases significantly on increasing temperature but is not significantly influenced by other factors.

3.2 Thermal decomposition of NPL-10

After determining the best conditions of use for NPL-10, we studied its thermal decomposition process by non-isothermal thermogravimetric analysis (TGA) and differential thermal analysis (DTA) [15, 16, 17, 18].

3.2.1 Non-isothermal thermogravimetric analysis (TGA) and differential thermal analysis (DTA)

Non-isothermal thermogravimetric analysis (TGA) and differential thermal analysis (DTA) were carried out at three different heating rates (5, 15, and 20 K·min–1). TGA indicated that, regardless of heating rate, NPL-10 starts decomposing at about 250°C, and decomposition is complete at about 500°C (Figure 2). DTA highlighted that thermal decomposition is a two-step process: the first decomposition step occurs at 200–300°C, with a mass loss of 4–26%; the second step occurs at 300–400°C, with 15–59% mass loss (Figure 3). Both steps are endothermic. Furthermore, both the rate of mass loss and the peak temperature of thermal decomposition gradually increase on increasing the heating rate.

Thermogravimetric analysis (TGA) curve of NPL-1.
Figure 2

Thermogravimetric analysis (TGA) curve of NPL-1.

Differential thermal analysis (DTA) curve of NPL-10.
Figure 3

Differential thermal analysis (DTA) curve of NPL-10.

3.2.2 Kinetics of thermal decomposition reaction

The thermal decomposition rate can be expressed by Eq. (1) as a function of two independent variables: the transformation rate α and temperature T [19]:

(dαdt)T=k(T)f(α)(1)

where T (K) is the absolute temperature of sample, α (%) is the decomposition rate, and k is the rate constant for the decomposition reaction. The heating rate β is constant in GTA, and DTA, and is expressed by Eq. (2):

β=dT/dt(2)

The rate constant k depends on temperature T according to the Arrhenius equation [20]:

k(T)=Aexp(EαRT)(3)

where A is the pre-exponential factor (s–1), Ea is the activation energy (J·mol–1), and R is the gas constant (8.314 J·mol–1·K–1).

Combining Eqs. (1), (2) and (3) provides Eq. (4), which is the rate law for the decomposition reaction:

dαdT=1βAexp(EαRT)f(α)(4)

We calculated the kinetic triplet (the activation energy Ea, the pre-exponential factor A, and the mechanism function g(α)) for both steps of NPL-10 decomposition process.

The activation energy Ea was obtained by Eq. (5) using KAS iterative method:

1nβh(x)T2=1n[AEag(α)R]EaRT(5)

where

h(x)=x4+18x3+88x2+96xx4+20x3+120x2+240x+120(6)

and

x=EaRT(7)

KAS iterative method is divided into three steps: ① Assume h(x)=1 to estimate the initial activation energy E1 for given β and g(α) values by Eq. (5). ② Bring E1 and the corresponding T into Eq. (7), then use the calculated x value to derive h(x) by Eq. (6); enter x and h(x) values into Eq. (5), and use the linear regression of ln[β/h(x)T2] vs 1/T to calculate the line slope (–E2/R) and hence E2. ③ Repeat step ② by replacing E1 by E2. When |EiEM| < 0.01 kJ·mol–1, Ei can be considered a real value.

The activation energy Ea was calculated for α = 0.2–0.8 at 0.05-unit intervals (Table 2). A plot of the activation energy Ea as a function of a for steps 1 and 2 (Figure 4) shows two curves of different shape, indicating that decomposition follows a different mechanism in steps 1 and 2. The non-linearity of curves suggests a complex mechanism for both steps.

Line chart for thermal decomposition of NPL-10.
Figure 4

Line chart for thermal decomposition of NPL-10.

Table 2

Activation energies calculated by KAS iterative method.

To determine the most probable kinetic functior g(α) for the thermal decomposition of NPL-10, we performed the linear regressen anaiyais of kinetic functions (Table 3) [18]. Eq. (5) was modified into Eq. (6):

Table 3

Linear regression analysis of 41 kinetic functions for steps 1 and 2.

1ng(α)=[1nAER+1nexx+1np(x)]1nβ(6)

where

p(x)=exx2h(x).

Linear regression analysis did not provide a straight line for four functions, namely 21 (dependent on pressure in nddhion to temperature, 25 (tew of Mampel Pewer), 39 and 40 (both functions of index law). The kinetic functions with the slope closest to –1 and best R2 values were g(α)= 1–(1–α)2 (Table 3, entry 33) for step 1 and g(α)= α1/2 for step 2 of thermal decomposition (Table 3, entry 31). These results indicate that the two steps of thermal decomposition follow a different order of reaction (n=2 for step 1 and n=1/2 for step 2).

Pre-exponential factors A were calculated by entering Ea values (Table 2) and g(α) values (Table 3) into Eq. (5). The pre-exponential factor of step 1 is A1=3.86·1011s–1, the pre-exponential factor of step 2 is A2=5.60·1013 s–1. Calculation of Ea, g(α), and A allowed to derive the rate law of the decomposition process. The rate law for the first step of thermal decomposition is:

dαdT=3.861011βexp(62.21103RT)[1(1α)2];

the rate law for the second step is:

dαdT=5.601013βexp(75.16103RT)α12.

4 Conclusions

  1. The optimal conditions (in terms of both foam height and half-life) for the use of foaming agent NPL-10 in heavy oil recovery by steam injection are the following: T=220°C, pH 7, 10000 mg·L–1 salinity, and 10 g·L–1 oil content;

  2. NPL-10 starts decomposing above 220°C;

  3. Decomposition of NPL-10 occurs in two steps that follow a different kinetic law, with the rate of the first step being more sensitive to temperature than the rate of the second step.

The derived kinetic model provides a guidance for the selection of the optimal quantity of NPL-10 in heavy oil recovery by steam injection.

Acknowledgment

this work was supported by the Natural Science Foundation of Heilongjiang Province of China (Project No. E2015036), the National Science and Technology Major Projects of China for Oil and Gas (Projects No. 2016ZX05055-006 and 2016ZX05012- 001), and the Cultivation Fund of Northeast Petroleum University of China.

References

  • [1]

    Koottungal L., Special Report 2010 worldwide EOR survey. Oil & Gas Journal, 2010, 108(14), 41-53. Google Scholar

  • [2]

    Kamari A., Nikookar M., Sahranavard L., et al., Efficient screening of enhanced oil recovery methods and predictive economic analysis. Neural Computing and Applications, 2014, 25(3-4), 815-824. CrossrefGoogle Scholar

  • [3]

    Hiraski G.J., The steam-foam process. Journal of Petroleum Technology, 1989, 41(05), 449-456. CrossrefGoogle Scholar

  • [4]

    Ma K., Lopez-Salinas J. L., Puerto M. C., et al., Estimation of parameters for the simulation of foam flow through porous media. Part 1: the dry-out effect. Energy & Fuels, 2013, 27(5), 2363-2375. Web of ScienceCrossrefGoogle Scholar

  • [5]

    Farajzadeh R., Krastev R., Zitha P.L.J., Foam films stabilized with alpha olefin sulfonate (AOS). Colloids and Surfaces A: Physicochemical and Engineering Aspects, 2008, 324(1), 35-40. CrossrefWeb of ScienceGoogle Scholar

  • [6]

    Farajzadeh R., Andrianov A., Krastev R., et al., Foam–oil interaction in porous media: Implications for foam assisted enhanced oil recovery. Advances in colloid and interface science, 2012, 183, 1-13. Web of SciencePubMedGoogle Scholar

  • [7]

    Ma K., Lopez-Salinas J.L., Puerto M.C., et al., Estimation of parameters for the simulation of foam flow through porous media. Part 1: the dry-out effect. Energy & Fuels, 2013, 27(5), 2363-2375. Web of ScienceCrossrefGoogle Scholar

  • [8]

    Worthen A.J., Bryant S.L., Huh C., et al., Carbon dioxide-in-water foams stabilized with nanoparticles and surfactant acting in synergy. AlChE Journal, 2013, 59(9), 3490-3501. CrossrefGoogle Scholar

  • [9]

    El-Amin M.F., Salama A., Sun S., Numerical and dimensional analysis of nanoparticles transport with two-phase flow in porous media. Journal of Petroleum Science and Engineering, 2015, 128, 53-64. Web of ScienceCrossrefGoogle Scholar

  • [10]

    Keijzer P.P.M., Muijs H.M., Janssen-van Rosmalen R., et al., Application of Steam Foam in the Tia Juana Field, Venezuela:Laboratory Tests and Field Results. SPE Enhanced Oil Recovery Symposium, 20-23 April, Tulsa, Oklahoma, 1986, SPE-14905. Google Scholar

  • [11]

    Maini B.B., Ma V., Thermal Stability of Surfactants for Steamflood Applications. SPE Oilfield and Geothermal Chemistry Symposium, 9-11 March, Phoenix, Arizona,1985, SPE-13572. Google Scholar

  • [12]

    Kleijnen J.P.C., Sensitivity analysis and related analyses: a review of some statistical techniques. Journal of Statistical Computation and Simulation, 1997, 57(1-4), 111-142. CrossrefGoogle Scholar

  • [13]

    Tang B., Orthogonal array-based Latin hypercubes. Journal of the American statistical association, 1993, 88(424), 1392-1397. CrossrefGoogle Scholar

  • [14]

    Kleijnen J.P.C., An overview of the design and analysis of simulation experiments for sensitivity analysis. European Journal of Operational Research, 2005, 164(2), 287-300. CrossrefGoogle Scholar

  • [15]

    Einaga H., Futamura S., Ibusuki T., Heterogeneous photocatalytic oxidation of benzene, toluene, cyclohexene and cyclohexane in humidified air: comparison of decomposition behavior on photoirradiated TiO2 catalyst. Applied Catalysis B: Environmental, 2002, 38(3), 215-225. CrossrefGoogle Scholar

  • [16]

    Patel P., Hull T.R., McCabe R.W., et al., Mechanism of thermal decomposition of poly(ether ether ketone)(PEEK) from a review of decomposition studies. Polymer Decomposition and Stability, 2010, 95(5), 709-718. CrossrefGoogle Scholar

  • [17]

    Zhao Z., Chaos M., Kazakov A., et al., Thermal decomposition reaction and a comprehensive kinetic model of dimethyl ether. International Journal of Chemical Kinetics, 2008, 40(1), 1-18. Web of ScienceCrossrefGoogle Scholar

  • [18]

    Robledo-Ortiz J.R., Zepeda C., Gomez C., et al., Non-isothermal decomposition kinetics of azodicarbonamide in high density polyethylene using a capillary rheometer. Polymer Testing, 2008, 27(6), 730-735. Web of ScienceCrossrefGoogle Scholar

  • [19]

    Kandelbauer A., Wuzella G., Mahendran A., et al., Model-free kinetic analysis of melamine–formaldehyde resin cure. Chemical Engineering Journal, 2009, 152(2), 556-565. Web of ScienceCrossrefGoogle Scholar

  • [20]

    Ptáček P., Bartoníčková E., Švec J., et al., The kinetics and mechanism of thermal decomposition of SrCO3 polymorphs. Ceramics International, 2015, 41(1), 115-126. Web of ScienceCrossrefGoogle Scholar

About the article

Received: 2017-08-18

Accepted: 2017-10-27

Published Online: 2018-02-13


Conflict of interest: Authors state no conflict of interest.


Citation Information: Open Chemistry, Volume 16, Issue 1, Pages 29–34, ISSN (Online) 2391-5420, DOI: https://doi.org/10.1515/chem-2018-0002.

Export Citation

© 2018 Fa-Jun Zhao et al., published by De Gruyter. This work is licensed under the Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 License. BY-NC-ND 4.0

Citing Articles

Here you can find all Crossref-listed publications in which this article is cited. If you would like to receive automatic email messages as soon as this article is cited in other publications, simply activate the “Citation Alert” on the top of this page.

[1]
Xiang Rao, Linsong Cheng, Lizhun Wang, and Mingqiang Chen
Arabian Journal of Geosciences, 2019, Volume 12, Number 9
[2]
Weitao Li, Falin Wei, Chunming Xiong, Jian Ouyang, Guang Zhao, Liming Shao, and Mingli Dai
Arabian Journal of Geosciences, 2019, Volume 12, Number 7

Comments (0)

Please log in or register to comment.
Log in