Improving the Aquathermolysis Efficiency of Aromatics in Extra-Heavy

Nov 2, 2015 - †Faculty of Materials Science and Chemistry, and ‡Key Laboratory of Tectonics and Petroleum Resources of Ministry of Education, Chin...
2 downloads 13 Views 492KB Size
Subscriber access provided by University of Idaho Library

Article

Improving the Aquathermolysis Efficiency of Aromatics in Extraheavy Oil By Introducing Hydrogen-donor Ligands to Catalysts Ruilin Ren, Huachao Liu, Yu Chen, Jian Li, and YanLing Chen Energy Fuels, Just Accepted Manuscript • DOI: 10.1021/acs.energyfuels.5b01256 • Publication Date (Web): 02 Nov 2015 Downloaded from http://pubs.acs.org on November 2, 2015

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

Energy & Fuels is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 30

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

1

Improving the Aquathermolysis Efficiency of Aromatics in

2

Extra-heavy Oil By Introducing Hydrogen-donor Ligands to

3

Catalysts

4

Ruilin Rena, Huachao Liua, Yu Chena, Jian Lib ,Yanling Chena,b*

5

a

6

Resources of Ministry of Education, China University of Geosciences, Wuhan 430074, PR China

7

*Corresponding author: Tel: +86-13886113362; Fax: +86-027-87801763.

8

E-mail address: [email protected] (Y. Chen).

9

Postal address: Faculty of Materials Science and Chemistry, China University of Geosciences,

10

Faculty of Materials Science and Chemistry, bthe Key Laboratory of Tectonics and Petroleum

Wuhan 430074, PR China

11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27

ACS Paragon Plus Environment

Energy & Fuels

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

28

Abstract

29

We develop a new type of catalyst, dimethylbenzene sulfonic copper, which composed by

30

hydrogen-donative ligands and a Cu2+ center for catalytic aquathermolysis viscosity reduction of

31

extra-heavy oil. The component variations of extra-heavy oil before and after aquathermolysis

32

are designedly analyzed via elemental analysis (EL), 1H nuclear magnetic resonance (1H-NMR)

33

and gas chromatography-mass spectrometry (GC-MS). An optimal viscosity reduction value of

34

81.47% can be achieved using 0.3 wt% catalyst at 240℃for 24h, where 8.28 wt% of the heavy

35

components, namely as resins and asphaltenes, are converted to light saturates and aromatics.

36

The improved viscosity reduction efficiency can be ascribed to the hydrogen donating role from

37

the dimethylbenzene sulfonic ligands. The active hydrogens, not only generated from these

38

ligands but originated from the dehydrogenation of naphthene aromatics during the

39

aquathermolysis, are capable of tackling the heavy components or participating in the

40

hydrogenation of heavy components for heavy-light conversion, leading to enhanced visbreaking

41

effect and improved oil quality. The results enrich the catalytic mechanism of aquathermolysis

42

which can guide the design and development of new efficient catalysts.

43 44

Keywords: Dimethylbenzene sulfonic copper; Catalytic aquathermolysis; Hydrogen-donating

45

effect; Heavy-light conversion; Viscosity reduction.

46 47

1. Introduction

48

Catalytic aquathermolysis is a promising technology for chemically mining extra-heavy oil.

49

Developing low cost, efficient and universal catalysts as well as understanding their catalytic

50

mechanisms are of essential importance for reducing the viscosity of extra-heavy oil to facilitate

51

field mining. Since the first concept of aquathermolysis proposed by Hyne et al.1 in 1982,

52

numerous efforts have been attempted to develop appropriate catalysts, majorly transition metal

53

ion coordinates, and to explore the possible catalytic viscosity-reduced mechanisms in

54

aquathermolysis2-13. Many new type of catalysts have been reported for reducing the viscosity of

ACS Paragon Plus Environment

Page 2 of 30

Page 3 of 30

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

55

extra-heavy oil, including mineral, water-soluble catalyst, oil-soluble catalyst, dispersed catalyst 2.

56

Although the reported catalysts exhibit improved viscosity reduction value comparing with bare

57

aquathermolysis, however, most of them failed to meet the industrial requirements for field

58

mining of extra-heavy oil.

59

The viscosity reduction of heavy or extra-heavy oils in a typical aquathermolysis process is due to

60

the cleavage of the C-R (R =S, N, O, etc.) bonds14, which will generate large amounts of free

61

radicals. However, these free radicals may repolymerize and lead to the viscosity regress after

62

reaction. It has been recognized that, introducing hydrogen donors, such as tetralin15, can inhibit

63

the generated free radicals in coal liquefaction and residuum hydrocrackings. This concept has

64

been naturally employed to deal with the catalytic aquathermolysis of heavy or extra-heavy oil.

65

Many types of hydrogen donors, involving formic acid16,17, methylbenzene18,19, tetralin20-23, and

66

methane24-30, have been utilized to couple with different types of catalysts for aquathermolysis of

67

heavy oils. Chen et al.18 has proven that, using 0.8 wt% of toluene as additive and 0.2 wt% of

68

ammonium molybdate as catalyst, the viscosity reduction value could achieve 85%. Ovalles et

69

al.26 demonstrated that methane, as the hydrogen donor (4.8 MPa), in conjugation with a dispersed

70

molybdenum catalyst(250 ppm) at 410 °C, leading to an increase of 7 degrees in the API gravity,

71

16% of reduction in sulfur content, and 55% conversion of the >500°C fraction with respect to the

72

original crude oil, which was also testified by Li et al.25. The hydrogen donating mechanisms of

73

these hydrogen donors were sequentially discussed.18,22,29,31-32 In general, a hydrogen donor not

74

only serves as solvent to dilute the macromolecular radicals but mitigates the collision of free

75

radicals and the condensation of heavy components, which is beneficial for decreasing the

76

viscosity of heavy oils to a large extent. However, adding excessive hydrogen donors, either in

ACS Paragon Plus Environment

Energy & Fuels

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

77

fine or large amount, would of course result in technical complexity in field oil exploration, and

78

may cause health, environmental or safety issues (HSE).

79

To simultaneously take account of both the catalytic function of the transition metal center and the

80

hydrogen-donating function of the hydrogen donors, it is not hard to speculate that if we can

81

synergistically couple these two functions in one catalyst molecule. Inspired by this, here we

82

designed a new type of catalyst (dimethylbenzene sulfonic copper). Through the elemental

83

analysis (EL), 1H nuclear magnetic resonance (1H-NMR) and gas chromatography-mass

84

spectrometry (GC-MS) characterizations, we emphatically analyzes the changes of aromatic

85

hydrocarbons which contain the structure parameters and the compositions before and after

86

reaction. The results show that the naphthene aromatics which have the potential

87

hydrogen-donating capability could be generated in the catalytic aquathemolysis process. While in

88

the presence of the new catalyst, the naphthene aromatics could donate hydrogen by aromatization.

89

Thus, the hydrogen-donating effect achieved the goal to upgrade the heavy oil and to reduce the

90

viscosity. The current study enriches the mechanism of catalytic aquathermolysis which can be

91

used to guide screening more efficient catalyst.

92

2. Experimental

93

2.1. Preparation of catalysts

94

Preparation of catalyst: 0.025mol copper (II) sulfate was dissolved in distilled water and an

95

aqueous solution of sodium hydroxide was dropwisely added into the water until a pale blue

96

slurry was completely formed under stirring condition. When the pH of the mixture was equal to

97

7 approximately, 0.10 mol dimethylbenzene sulfonic acid was placed into a beaker in 100°C

98

oil-bath for 20 minutes. Then the prepared slurry was added into the beaker under the condition

ACS Paragon Plus Environment

Page 4 of 30

Page 5 of 30

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

99

of stirring and heated for 2h. The product was cooled to room temperature. Upon concentrating

100

and drying in vacuum at 60°C for 24 h, dimethylbenzene sulfonic copper was obtained.

101

2.2. Experiment of catalytic aquathermolysis

102

The extra-heavy oil with the viscosity of 178,000 mPa·s at 70°C and the density of 0.9915 kg/m3

103

at 20℃ extracted from Shengli Oilfield of China was selected in the experiment of catalytic

104

aquathermolysis. All experiments were conducted with 180g oil samples with the ratio of

105

oil/water of 8:2, and the water pH value of approximately 7~8. The quality of the catalyst is 0.3

106

percent of the crude oil quality (the optimized conditions were confirmed after a series of tests).

107

The aquathermolysis was carried out in a 500 mL FYX-0.5 high pressure reactor, and the initial

108

pressure was kept at about 3 MPa by aerating N2. Then the reaction proceeded at 240°C for 24h.

109

The viscosity of oil samples I (crude heavy oil prior to reaction), oil sample II (crude heavy oil

110

after aquathermolysis without catalyst), oil sample III (the post-reacted crude heavy oil sample

111

with the as-prepared catalyst) were tested using the BROOKFIELD DV-2+ programmable

112

Viscometer at 70°C. The viscosity reduction values (Table1) obtained were calculated by ∆ η =

113

((η0 – η)/ η0) × 100%, where η0 (mPa·s) is the viscosity of extra-heavy oil, and η (mPa·s) is the

114

viscosity of the oil samples after reaction9.

115

2.3. Separation of four group compositions

116

Four group compositions including the saturate hydrocarbons, aromatic hydrocarbons, resins, and

117

asphaltenes of the oil samples were separated by a chromatography column according to the

118

Industrial Specification of China Petroleum Standard (SY/T5119)33. Firstly, about 1g oil sample

119

was dissolved into 100ml n-hexane through ultrasonic, and stood for 12h to precipitate asphaltene.

120

After filtering out asphaltene, the filtrate was separated by chromatography column(1.5cm in the

ACS Paragon Plus Environment

Energy & Fuels

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

121

diameter, 40cm in the length) filled with the active alumina (activated at 450℃ for 6h, with 1

122

wt% water). Saturated hydrocarbons, aromatics and resins were eluted out respectively by 120ml

123

n-hexane solution, 120ml n-hexane mixed with dichloromethane in the volume ratio of 1:1, and

124

100ml dichloromethane mixed with methanol in the volume ratio of 1:1. (Fig. S1. Support

125

information). The water in oil samples was already removed before and after reaction. Firstly, the

126

sample oil was stirred and allowed to separate into layers after reaction, then the free water was

127

decanted. The oil sample was dehydrated by ultrasonic crude rapid dehydration instrument

128

(BT-350, Hubei Broad Electronics Co., Ltd., China). The water contents of the dehydrated

129

samples were determined via the Karl Fischer (Coulomb) SF101 Tester (Precision Instrument Co.,

130

Ltd., Taizhou, Jiangsu province, China) following the GB/T1146-2009 standard of People's

131

Republic of China. The measured water contents of the samples, either before or after the catalytic

132

aquathermolysis, are all below 0.002%. The content of four group compositions was obtained by

133

weight (Table1). Then, saturate hydrocarbons and aromatic hydrocarbons were prepared in a

134

solution of n-hexane, ready for the GC-MS analysis.

135

2.4. Measurement and Characterization

136

2.4.1. Elemental analysis

137

Organic elements (C, H, N, S) of the aromatics, resins and asphaltenes before and after reaction

138

were detected by an elemental analyzer (EL, VARIO EL-2).

139

2.4.2. 1H-NMR analysis

140

The aromatics in three oil samples before and after reaction were measured by 1H nuclear

141

magnetic resonance spectroscopy (1H−NMR, Bruker ARX300), using CDCl3 as solvent,

142

containing TMS as internal reference.

ACS Paragon Plus Environment

Page 6 of 30

Page 7 of 30

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

143

2.4.3. GC-MS analysis

144

Saturates and aromatics in three oil samples before and after reaction were analyzed by using gas

145

chromatography-mass spectrometer (GC-MS, Agilent 7890A/5975C). Saturates was injected in

146

the GC injection port (300°C) and operated with the splitless mode. GC oven was kept at the

147

temperature of 80°C for 3 minutes and then kept the incremental rate of 3°C /min until the

148

temperature became as high as 230°C. After that, the temperature increment was maintained at

149

2°C /min until it rose to 310°C, at which the reaction was kept for 15 minutes. The MS was

150

operated by using electron impact ionization with the default setting. While aromatic

151

hydrocarbons were injected into the GC injection port (290°C) and operated with the splitless

152

mode. The GC oven was programmed to retain the temperature at 60°C for 3 minutes and then

153

kept the incremental rate of 3°C /min until the temperature was up to 310°C. After that, the

154

reaction was maintained for 15 minutes sequentially. The MS was operated by using electron

155

impact ionization with the default setting.

156

3. Results and discussion

157

3.1. Aquathermolysis experiments

158

Table 1 shows the the changes of four components, viscosity, viscosity reduction value and

159

density and in extra-heavy oil. The results are the average of three parallel experiments. The

160

viscosity reduction value of oil sample Ⅱ was only 32.93%, compared with the oil sampleⅠ.

161

But in the participation of catalyst, the viscosity reduction value reached 81.47%. Obviously,

162

only a small part of heavy components cracked in aquathermolysis for the oil sample Ⅱ. For the

163

oil sample Ⅲ, the intensive conversion reaction occurred in the heavy component in the catalytic

164

aquathermolysis. This may indicate that light components have been generated through this

ACS Paragon Plus Environment

Energy & Fuels

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

165

cleavage reaction.

166

3.2. Elemental analysis before and after reaction

167

Organic elements (C, H, N, S) of aromatics, resins and asphaltenes in the three types of oil

168

samples were measured by EL as shown in Table 2. Table 2 also shows the molar H/C ratio of

169

each component before and after the reaction to compare with the crude oil after the catalytic

170

aquathermolysis reaction.

171

The mass balance of N and S were calculated according to the equation Nt= Wc∗Mc, where Wc is

172

the group composition mass fraction and Mc is the N, S elemental content in the group

173

compositions. Thus we obtained the total content (Nt) per 10000g oil samples. Table S1 and Table

174

S2 (Supporting Information) show the results of N and S elemental content (weight ratio) in

175

aromatics, resins, and asphaltenes of oil samples before and after the reaction. The saturates were

176

ignored, because the separation of the four group compositions depends on the polar property.

177

The decreased content of asphaltenes in oil sample Ⅱ indicates that a water-gas reaction may

178

have taken place and a part of asphaltenes may have been cracked into the light component (the

179

increased molar H/C ratio of aromatic hydrocarbon) based on the analysis of the data shown in

180

Table 1 and Table 2. The remarkable content decrease of resins and asphaltenes in oil sample Ⅲ

181

indicates that the catalyst has the ability to accelerate resins and asphaltenes to split into light

182

components, leading to an increase of molar H/C ratio of aromatic hydrocarbon and a decrease of

183

molar H/C ratio of asphaltenes and resins. It suggests that the catalyst is capable of accelerating

184

the catalytic dehydrogenation of aromatics to generate more free radicals, and to impede the

185

condensation polymerization of resins and asphaltenes.

186

The variation tendency of mass balance of element N and S is similar to the data discussed above.

ACS Paragon Plus Environment

Page 8 of 30

Page 9 of 30

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

187

This result further proves that the catalyst can catalyze the dehydrogenation reaction of aromatic

188

species and thereby promote the heavy components into light components34,35.

189

3.3. 1H-NMR analysis of the aromatic hydrocarbons before and after reaction

190

The distributions of different hydrogen (HA(9.0-6.0), Hα(4.0-2.0ppm), Hβ(2.0-1.0ppm), and

191

Hγ(1.0-0.5ppm)) (Table 3) were obtained based on different chemical shifts and integral area36.

192

Using the improved Brown-Ladner methods with the Equations (1), (2), (3) and (4),37 the

193

structure parameters of aromatic hydrocarbons such as aromaticity (fA), aromaticity condensation

194

(HAU/CA), branching index (BI) and substitution rate of aromatic ring (σ) (Table 3) were

195

calculated by the formula as shown in Supporting Information. 1H-NMR spectra of aromatic

196

compositions are in Fig.S2.(Support information).

197

Table 3 also shows the changes of fA , HAU/CA , BI and σ of all the samples before and after the

198

reaction. It demonstrates that the decreased fA value of aromatic hydrocarbon is attributed to the

199

aquathermolysis of the heavy components in the oil samples II. The aquathemolysis produced

200

small molecules with a high molar H/C ratio that enter into the light components. The increased fA

201

of aromatic hydrocarbon found only in the oil sample Ⅲ indicates that an aromatication reaction

202

occurred. This reaction may be viewed as a hydrogen-donating reaction. However, this

203

aromatication reaction takes place only in the presence of the prepared catalyst. The increment of

204

the HAU/CA ratio indicates that no aromatic hydrocarbon condensation reaction occurred. The

205

decrease of BI demonstrates that the branch of the side-chain of the aromatics may have fractured.

206

The increase of σ shows that a substitution reaction occurred in the aromatic ring. This substitution

207

reaction may be an effective way to terminate the active alkyl chains.

208

The above analyses lead to the conclusion that the aromatic hydrocarbons not only quench the

ACS Paragon Plus Environment

Energy & Fuels

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

209

alkyl chain radicals via the substitution reaction on the aromatic ring but also function as a

210

hydrogen donor via the aromatization reaction in the presence of the catalyst.

211

3.4. GC-MS analysis of aromatic hydrocarbons before and after reaction

212

Aromatization reaction of aromatic hydrocarbon provides an hydrogen source for the heavy

213

components in the catalytic aquathermolysis, while the naphthene aromatic hydrocarbon is one of

214

the ideal substances for hydrogen-donating . In order to find the source of the naphthene aromatics,

215

biomarked compounds monoaromatic sterane and triaromatic sterane were selected for a

216

comparative analysis using GC-MS because of their stable structures and their relationship to

217

hydrogenationship and dehydrogenation in content change.38

218

The biomarker monoaromatic steranes (m/z=253) (Fig.1) and triaromatic steranes (m/z =231)

219

(Fig. 2) were extracted from the total ion chromatorgraphy of the aromatic hydrocarbons in the

220

oil samples before and after the reaction. The compound identification of monoaromatic steranes

221

and triaromatic steranes are shown in Table 439-42.

222

As shown in Fig. 1 and Fig. 2, the content of monoaromatic sterane and triaromatic sterane in the

223

oil sample II is increased, compared with the oil sample I. They were produced through the

224

cracking and hydrogenation of heavy components in the extra-heavy oil. There was little

225

monoaromatic sterane left in the oil sample III, and the content of triaromatic sterane was increased.

226

This is attributed to the dehydrogenation reaction of monoaromatic sterane in the presence of the

227

catalyst, which provides active hydrogen for the cracking reaction of heavy components, leading to

228

the increase of triaromatic sterane.

229

These results indicate that the aromatic side chains of aromatic hydrocarbon in the oil sample II

230

were fractured during the catalyst aquathermolysis. For the oil sample III, the aromatic

ACS Paragon Plus Environment

Page 10 of 30

Page 11 of 30

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

231

hydrocarbon with unsaturated side chains and tetralin was disappeared. This indicates that the

232

catalyst not only catalyzed the cyclization but also promoted dehydrogenation aromatization

233

reaction, leading to the increase of fA of aromatics in the oil sample III. This is consistent with the

234

result shown in section 3.3.

235

A detail assessment of the aromatic hydrocarbons was made to better understand the

236

hydrogen-donating mechanism of aromatic hydrocarbon. The selected ion mass chromatogram

237

(m/z=132) of aromatic compositions in three oil samples was investigated. As shown in Fig. 3, the

238

aromatic substances with longer unsaturated side chains, such as 1-phenyl-1-nonyne,

239

octynyl-benzene ethanol, (3, 7-dimethyl-octa-2, 6-dienyl)-benzene, etc. (Fig. 3a and Table 5)

240

mainly exist in crude oil. The substances, such as 2-butenylbenzene and 1-methyl-1-propenyl-

241

benzene, etc. (Fig.3.b & Table 5), contain shorter unsaturated side chains detected in the oil sample

242

II. However, the above substances such as aromatics with unsaturated side-chain and tetralin were

243

disappeared in the oil sample III (Fig. 3c and Table 5). The results of the ion mass chromatograms

244

of naphthalenes (Fig. 4 and Table 6) obviously agree with the results above (Fig. 3). Furthermore,

245

the abundance of n-alkane was increased, which indicates that a dealkylation reaction took place in

246

the aromatic cycles. The iron mass chromatograms of n-alkane(m/z=85)

247

support the conclusions.

248

From the above analysis, it is clear that the naphthenic aromatic hydrocarbons come from three

249

sources: the crude oil, the thermolysis of heavy components and the hydrogenation of polycyclic

250

aromatic hydrocarbons. The naphthene aromatics donate hydrogen through aromatization in the

251

catalytic aquathermolysis.

ACS Paragon Plus Environment

43,44

(Fig. 5) could also

Energy & Fuels

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

252

In order to avoid the impact of the product derived from the catalysts after aquathermolysis,

253

experiments were performed in the absence of heavy oil. In the experiments, the hydrothermal

254

reactions of the as-prepared catalyst was carried out at 240°C for 24h. Subsequently, the products

255

were detected with GC-MS (Fig. S3, Supporting Information). The results (Table S3 which

256

referred to Fig. S3, Supporting Information) show that the hydrocarbons discussed in the paper

257

were not derived from these catalysts. This, in turn, proves that the hydrocarbons are derived from

258

the heavy oil. We also added the thermal stability characterization of the as-prepared catalyst in

259

the supporting information (Fig. S4). The TG-DSC curve shows that the catalyst is very stable at

260

the aquathermolysis temperature (240°C), which suggest that under our experimental condition,

261

the catalyst could maintain its structure and function. It can completely eliminate the influence of

262

the catalyst on the composition and structure changes of components in extra-heavy oil.

263

3.5 The transformation of various components and the hydrogen-donating mechanism in

264

extra-heavy oil

265

The above results show that, as the designed catalyst with hydrogen-donative ligands participating

266

in the catalytic aquathermolysis of heavy oil, a satisfactory viscosity reducing rate can be achieved,

267

which accompanies with the simultaneous dehydrogenation aromatization of naphthenic aromatics.

268

Interestingly, although the addition amount of hydrogen donors are relatively higher than the

269

introduced hydrogen-donative ligands in our experiments, the dehydrogenation aromatization

270

could not even be observed. This motivate us to speculate realistic role that the hydrogen-donative

271

ligands play during the aquathermolysis process. Considering the very few adding amount of the

272

designed catalyst (0.3 wt%), we can only expect the active hydrogen split from the ligand should

273

facilitate the initial stage of the catalytic viscosity reduction where partial of resins and

ACS Paragon Plus Environment

Page 12 of 30

Page 13 of 30

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

274

asphaltenes can be dissociated to saturates and aromatics with lower C/H ratio. In particular, under

275

our catalytic condition, the naphthenic aromatics such as polyaromatic steranes can successively

276

undergo dehydrogenation aromatization process which could generate more active hydrogens for

277

further attacking the intact heavy components in conjugation with the catalytic activity of the Cu2+

278

centers. Fig.6 shows the dehydrogenation aromatization reaction45,46 of naphthenic aromatics

279

occurred in catalytic aquathermolysis. Firstly, Cyclization reactions of aromatic hydrocarbons with

280

the unsaturated side chains occurred and generated into naphthenic aromatic hydrocarbons. Then

281

through further catalytic reaction, the naphthenic aromatics ’aromatization reaction happened and

282

produced active hydrogen. This dehydrogenation aromatization reaction provides a more active

283

hydrogen source for the resins and asphaltenes. At last, those active hydrogen entered into the

284

heavy components. On one hand, they can participate in the hydrogenation reaction of the heavy

285

components, on the other hand can inhibit the condensation of heavy components, such a

286

self-assistant process promising a continuous viscosity reduction efficiency of heavy oil.

287

4. Conclusion

288

Extra-heavy oil, one of the great potential hydrocarbon resources, contains many kinds of

289

substances with potential hydrogen-donating ability. Stimulated by the new catalyst, the quality

290

of the extra-heavy oil can be significantly improved and the viscosity can be remarkably reduced

291

by making full use of the abundant potential hydrogen donating source of the extra-heavy oil

292

itself. It is of tremendous practical significance to further study the aquathermolysis viscosity

293

reduction of extra-heavy oil through catalytic thermal cracking of autonomous hydrogen.

294 295 296 297

ACS Paragon Plus Environment

Energy & Fuels

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

298

Acknowledgements

299

This work was supported by the National Natural Science Foundation of People’s Republic of

300

China(No.51174179)and the Key Laboratory of Tectonics and Petroleum Resources of Ministry

301

of Education(China University of Geosciences).

302 303 304 305 306 307 308 309 310 311 312 313 314 315 316 317 318 319 320 321 322 323 324 325 326 327 328 329 330 331 332 333 334 335 336 337 338

ACS Paragon Plus Environment

Page 14 of 30

Page 15 of 30

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

339

References

340

[1] Hyne, J.; Greidanus, J.; Tyrer, J.; Verona, D.; Rizek, C. ; Clark, P. et al. Aquathermolysis of heavy oils, in

341

Proceedings of the 2nd International Conference on Heavy Crude and Tar Sands. Caracas, Venezuela, 1982; pp

342

25-30.

343

[2] Maity, S. K. ; Ancheyta, J. ; Marroquın, G. Energy & Fuels. 2010, 24(5), 2809-2816.

344

[3] Clark, P. D.; Hyne, J. B.; Tyrer, J. D. Fuel. 1983, 62 (5), 959-962.

345

[4] Clark, P. D.; Hyne, J. B.; Tyrer, J. D. Fuel. 1987, 66(5), 1353-1357.

346

[5] Clark, P. D.; Hyne, J. B.; Tyrer, J. D. Fuel. 1987;66(5):1699-1702.

347

[6] Fan, H.F.; Liu, Y.J.; Yang, F.L. Oilfield Chemistry. 2001, 18(1),13-16.

348

[7] Clark, P. D.; Kirk, M. J. Energy & Fuels. 1994, 8, 380 -387.

349

[8] Fan, H.F.; Liu, Y.J.; Zhao, X.F.; Zhong, L.G. Energy&Fuels. 2001, 29(5), 430-433.

350

[9] Chen, Y.L.; Wang, Y.Q.; Lu, J.Y.; Wu, C. Fuel. 2009, 88(8), 1426-1434.

351

[10] Chen, Y.L.; Wang, Y.Q.; Wu, C.; Xia, F. Energy & Fuels. 2008, 22(3), 1502-1508.

352

[11] Chen, Y.L.; Yang, C.; Wang, Y.Q. Journal of Analytical and Applied Pyrolysis. 2010, 89(2), 159-165.

353

[12] Chao, K.; Chen, Y.L; Li J.; Zhang, X.M.; Dong, B.Y. Fuel Processing Technology. 2012, 104, 174-180.

354

[13] Chao, K.; Chen, Y.L; Liu, H.C.; Zhang, X.M.; Li, J. Energy & Fuels. 2012, 26(2), 1152-1159.

355

[14] Wang, Y.Q.; Chen, Y.L.; He, J.; Li, P.; Yang, C. Energy& Fuels. 2010, 24, 1502−1510.

356

[15] Neavel, R. C. Fuel. 1976, 3(55), 237.

357

[16] Liu, Y.J; Zhao, F.J; Zhao, G.; Hu, S.B.; Lin, Z. Oilfield Chemistry. (Chin.).2008, 25(2), 133-136.

358

[17] Zhao, F.J.; Liu, Y.J.; Zhao, G.; Hu, S.B. Chemical industry and engineering progress. (Chin.). 2008, 27(9),

359

1453-1459.

360

[18] Chen, E.Y.; Liu, Y.J; Wen, S.B. Daqing Petroleum Institute. (Chin.).2005, 29(6), 38-39.

361

[19] Wu,C.;Chen,Y.L.;Wang,Y.;Yang,C. Journal of Oil and Gas Technology .(Chin.).2007, 29, 267-269.

ACS Paragon Plus Environment

Energy & Fuels

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

362

[20] Fan, H. F.; Liu, Y. J. Energy & Fuels. 2002, 16, 842−846.

363

[21] Fan, Z.X.; Zhao, F.L.; Wang, J.X.; Gu Y.G. Journal of Fuel Chemistry and Technology. (Chin.). 2006, 34(3),

364

315-318.

365

[22] Li, B. Daqing Petroleum Institute. (Chin.). 2004, 4(28), 24-26.

366

[23] Ovalles, C. Extra-heavy crude oil downhole upgrading process using hydrogen donor under steam injection

367

conditions [Z], SPE 69692, 2001.

368

[24] Ovalles, C.; Vallejos, C.; Vasquez, T.; Rojas, I.; Ehrman, U.; Benitez, J.L.; Martinez. R. Petroleum Science

369

and Technology. 2003,21 (1-2), 255-274.

370

[25] Li, S.F.;Mu, B.Q.;Liu, C.G. Chemicial engnierring of oil and gas. (Chin.). 2002, 31(3), 138-139.

371

[26] Ovalles, C.; Filgueiras, E.; Morales A, et al. Energy& Fuels, 1998,12(2), 379–385.

372

[27] Ovalles, C.; Antonia, H.; Rojas, I.; et al. Fuel.1995,74(8), 1162-1168.

373

[28] Ovalles, C.; Hamana, A.; Bolivar, R.; et al. Process for treating heavy crude oil [P], US 5269 909, 1993.

374

[29]Ovalles, C; Filgueiras, E; Morales, A; et al. Fuel. 2003, 8(82), 887 -892.

375

[30] Vallejos, C.;Vasquez. T.; Ovalles, C.; et al. Process for the downhole upgrading of extra heavy crude oil [P],

376

US 5 891 829, 1999.

377

[31] A.Del Bianco.; Garuti, G.; Pirovano C.; et al. Fuel. 1995, 74(5), 756 -760.

378

[32] Senthilnathan, V. P.; Stein, S. E. J Organic Chem. 1988, 53(13), 3000 -3007.

379

[33] Analysis of Family Composition of Rock Extract and Crude Oil by Coumn Chromatography, China Petroleum

380

Standard SY/T 5119-1995; China National Petroleum Corporation, Dec.25,1995.

381

[34] Muraza, O.; Galadima, A. Fuel. 2015, 157, 219-231.

382

[35] Wu, C.; Lei G.L; Yao, C.J; Gai, P.Y; Cao, Y.B.; Li, X.N. Journal of China University of

383

Petroleum.(Chin.).2011, 35(1), 164-168.

ACS Paragon Plus Environment

Page 16 of 30

Page 17 of 30

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

384

[36] W. Liang. Heavy oil chemistry, in:China University of Petroleum publishing house, Beijing, China University

385

of Petroleum Press (Dongying);2000; pp 52-54.

386

[37] Liang, W.J.; Que, G.H.; Chen, Y.Z. Chemcal composition and structure of vacuum residues of Chinese crudes

387

II. average structure of vacuum residues and their fractions.Acta Petrolei Sinica. 1991, 7 ,1-11.

388

[38] Bennett, B.; Larter, S.; Carbognanil; et al. Energy & Fuels. 2007, 22(1), 440-448.

389

[39] Isaac Kaplan; Lu, S.T.; Lee, R.P.; Warrick G. Environmental toxicology and chemistry. 1996, 15(5), 696-707.

390

[40] Fang, Y.X.; Liao, Y.H.; Wu, L.L.; Geng, A.S. Journal of Asian Earth Sciences. 2011; 41(2), 147-158.

391

[41] Barbanti, S.M.; Moldowan, J.M.; Watt, D.S.; Kolaczkowska, E. Organic geochemistry. 2011, 42(4). 409-424.

392

[42] Fustic, M.; Bennett, B.; Huang, H.; Larter, S. Marine and petroleum geology. 2012, 36(1), 50-69.

393

[43] Chen, Y.L.; He, J.; Wang, Y.Q. Energy. 2010, 35(8), 3454-3460.

394

[44] Elias, R.; Vieth, A.; Riva, A.; et al. Organic geochemistry. 2007, 38(12), 2111-2130.

395

[45] Lu, S.R.; Peng, B.; Bao, M. Process in chemistry.2007,19(10),1529-1536.

396

[46] Masahito M.; Naoki H.;Dsvid R.; et al. Organic Letters. 2014, 16, 4134-4137

397 398 399 400 401 402 403 404 405

ACS Paragon Plus Environment

Energy & Fuels

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

406 407 408 409 410

Page 18 of 30

Tables

Table 1 Four group compositions, viscosity, viscosity reduction value and density in different oil samples before and after reaction Oil samples Group compositions (wt %)

Oil sample Ⅲ

Saturates

18.82±0.05

20.48±0.03

22.53±0.03

25.77±0.06

26.97±0.05

30.34±0.05

resins

31.27±0.05

31.83±0.07

27.89±0.07

24.14±0.01

20.72±0.01

19.24±0.01

178,000

119,384

32900

0.9915

0.9558

0.9017



32.93%

81.47%

Asphaltane Viscosity(at 70 C) mPa·s ∆η

Oil sample Ⅱ

Aromatics

0

Density (at 200C)

Oil sample Ⅰ

Kg/m3

411

412 413

ACS Paragon Plus Environment

Page 19 of 30

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

414 415 416

Table 2 Percentage of hydrogen,carbon,nitrogen , sulfur elements and molar H/C ratio of aromatics, resins and asphaltenes in different oil samples Group

compositions

Aromatics

Resins

Asphaltenes

85.591

82.163

82.013

10.843

9.459

8.268

0.739

1.877

2.094

S

2.557

1.898

2.596

NH/NC

1.52

1.381

1.21

C

85.527

83.968

83.71

H

11.237

9.858

8.379

0.753

1.887

2.125

S

2.606

2.053

2.528

NH/NC

1.577

1.409

1.201

C

86.652

82.959

83.573

10.885

10.057

8.458

C H Oil sample Ⅰ

Oil sample Ⅱ

N

N

H Oil sample Ⅲ

N

(%)

(%)

(%)

0.814

1.841

2.187

S

2.720

2.039

2.496

NH/NC

1.507

1.455

1.215

417 418 419 420 421 422

ACS Paragon Plus Environment

Energy & Fuels

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

423 424 425

Page 20 of 30

Table3 1H-NMR analysis of aromatic compositions and aromatic structure parameters of the aromatic compositions in different oil samples Aromatics of Oil sample Ⅰ

Aromatics of oil sample Ⅱ

Aromatics of oil sample Ⅲ

13.89

16.07

16.79

54.85

53.89

53.64

22.92

21.56

20.44

HA

8.34

8.47

9.13

fA

0.303

0.279

0.315

HAU/CA

0.766

0.934

0.838

BI

0.222

0.205

0.193

σ

0.454

0.487

0.479

Aromatic hydrocarbons Hα composition

Hβ Hγ

(%)

426

ACS Paragon Plus Environment

Page 21 of 30

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

427 428 429

Table 4 The identification of monoaromatic sterance (m/z=253) and triaromatic sterance (m/z= 231) m/z

Peak

Chemical formula

Compounds

253

1

C22H32

C22 monoaromatic sterane (MAS)

2

C20H14

1,1-binaphthyl

3

C27H42

20R,5β C27 + 20R C27 dia MAS

4

C28H44

20S,5β C28 MAS

5

C29H46

20S,5β C29 + 20R,5β C29 MAS

6

C28H44

20S,5α C28 MAS

7

C29H46

20R,5β C29 + 20R,5α C29 MAS

8

C29H46

20R,5α C29 MAS

1

C20H20

C20 triaromatic sterane (TAS)

2

C21H22

C21 TAS

3

C26H32

20S,C26 TAS

4

C26H32 + C27H34

20R,C26 + 20S,C27 TAS

5

C28H36

20S,C28 TAS

6

C27H34

20R,C27 TAS

7

C28H36

20R,C28 TAS

231

430 431

ACS Paragon Plus Environment

Energy & Fuels

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

432 433 434

Table 5 The identification of ion mass chromatogram (m/z=132) of the aromatic hydrocarbons before and after reaction Peak

Chemical formula

Compounds

1

C10H12

2-butenylbenzene

2

C10H12

1-butenylbenzene

3

C10H12

o-isopropenyltoluene

4

C10H12

2-methylallylbenzene

5

C10H12

2-methyl-1-propenylbenzene

6

C10H12

2-butenylbenzene

7

C10H12

1-methyl-1-propenylbenzene

8

C10H12

1,2,3,4-tetrahydronaphthalene

9

C11H14

1,2,3,4-tetrahydro-1-methylnaphthalene

10

C11H14

2,3,-dihydro-4.7-dimethyl-1H-indene

11

C11H14

2-methyl-1-butenyl-benzene

12

C12H16

1,2,3,4-tetrahydro-1,4-dimethylnaphthalene

13

C12H12

1,4-cyclohexadien-1-ylbenzene

14

C13H16O

2-benzyllidenecyclohexanol

15

C12H15N

1,2,3,4-tetrahydro-2-naphthyl-1-aziridine

16

C14H20

1,4-dimethyl-2-cyclohexylbenzene

17

C13H14

2,3,6-trimethylnaphthalene

18

C15H24O

butylated hydroxytoluene

19

C15H20

1-phenyl-1-nonyne

20

C14H16O

1,2-bis(1-buten-3-yl)-benzene

21

C16H22O

octynyl-benzeneethanol

22

C14H20

1,2,3,4-tetrahydro-5,6,7,8-tetramethyl-naphthalene

23

C16H22

(3,7-dimethyl-octa-2,6-dienyl)-benzene

24

C15H18O

3a,9b-dimethyl-1,2,3a,4,5,9b-hexahydrocyclopenta-naphthalen-3-one

25

C15H20

9-methyl-S-octahydroanthrancene

435 436

ACS Paragon Plus Environment

Page 22 of 30

Page 23 of 30

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

437 438

Table 6 The identification of naphthalene and its homologous series (m/z=128+142+156) m/z

peak

Chemical formula

Compounds

128

1

C10H8

naphthalene

142

2

C11H10

1-methyl-naphthalene + 2-methyl-naphthalene

156

3

C12H12

1-ethyl-naphthalene + 2-ethyl-naphthalene

4

C12H12

2,3-dimethyl-naphthalene + 1,4-dimethyl-naphthalene

5

C12H12

2,3-dimethyl-naphthalene + 1,3-dimethyl-naphthalene

6

C12H12

2,6-dimethyl-naphthalene

7

C12H12

1,2-dimethyl-naphthalene

439 440 441 442 443 444 445 446 447 448 449 450 451 452 453 454 455 456 457 458 459 460 461 462 463 464 465 466 467

ACS Paragon Plus Environment

Energy & Fuels

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

468 469 470 471 472 473 474 475 476 477 478 479 480 481 482 483 484 485 486 487 488 489 490

Figure captions Fig.1 Monoaromatic steranes ion mass chromatogram (m/z=253) of aromatic hydrocarbons (a. Aromatic hydrocarbons of crude oil sampleⅠ b. Aromatic hydrocarbons of oil sampleⅡ c. Aromatic hydrocarbons of oil sample Ⅲ) Fig.2 Triaromatic steranes ion mass chromatogram (m/z=231) of aromatic hydrocarbons (a. Aromatic hydrocarbons of crude oil sampleⅠ b. Aromatic hydrocarbons of oil sampleⅡ c. Aromatic hydrocarbons of oil sampleⅢ) Fig.3 The selected ion mass chromatogram (m/z=132) of aromatic compositions (a. Aromatic hydrocarbons of crude oil sampleⅠ b. Aromatic hydrocarbons of oil sampleⅡ c. Aromatic hydrocarbons of oil sampleⅢ) Fig.4 Naphthalene and its homologous series ion mass chromatograms (m/z=128+142+156) of aromatic compositions (a. Aromatic hydrocarbons of crude oil sampleⅠ b. Aromatic hydrocarbons of oil sampleⅡ c. Aromatic hydrocarbons of oil sampleⅢ) Fig.5 The ion mass chromatograms of alkanes (a. Saturate hydrocarbons of crude oil sampleⅠ b. Saturate hydrocarbons of oil sampleⅡ c. Saturate hydrocarbons of oil sampleⅢ) Fig.6 The dehydrogenation aromatization of aromatic hydrocarbons with the unsaturated side chains in catalytic aquathermolysis of super-heavy oil with viscosity reducer

491 492 493 494 495 496 497 498 499 500 501 502 503 504 505 506 507

ACS Paragon Plus Environment

Page 24 of 30

Page 25 of 30

508 509

Fig. 1

20000

c

2 8 7 6 5 43

Abundance

15000

1

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

b

10000

5000

a

0 40

50

60

70

Time(min)

510

ACS Paragon Plus Environment

80

Energy & Fuels

511 512 513

Fig. 2

260000 240000 220000

5 4 7

200000

Abundance

1

6

160000

3

180000

2

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 26 of 30

c

140000 120000 100000

b

80000 60000 40000

a

20000 0 60

70

Time(min)

514 515

ACS Paragon Plus Environment

80

Page 27 of 30

516 517 518

Fig. 3

6000 18

5000

7

4000

Abundance

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

16

c

16 18

b

7

3000

6

8

12

9

10 11

13 14 15 17 20

2000

24

18 19 23 21

1000

a 25

34

5

16 15

12

8

0 10

15

20

25

30

Time(min)

519 520

ACS Paragon Plus Environment

35

40

Energy & Fuels

521 522 523

Fig. 4

200000

2

300000

7 6 5 4 3

150000

1

250000

100000 50000

Abundance

200000 0

2

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 28 of 30

23

24

25

26

27

28

29

30

31

150000

100000

c

50000

b a

0 10

20

30

40

50

60

Time(min)

524

ACS Paragon Plus Environment

70

80

Page 29 of 30

525 526 527

Fig. 5

C 25

80000 C 21

C 17

60000

Abundance

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

ph Pr

c

40000

b 20000

a 0

20

40

60

80

Time(min)

528 529 530 531 532 533 534 535 536 537 538 539 540 541 542 543 544 545 546 547 548 549 550 551

ACS Paragon Plus Environment

100

Energy & Fuels

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

552 553

Fig. 6

554 555 556

ACS Paragon Plus Environment

Page 30 of 30