Subscriber access provided by UNIV OF NEWCASTLE
Article
Mineral Reactions in Shale Gas Reservoirs: Barite scale formation from reusing produced water as hydraulic fracturing fluid Amelia N. Paukert Vankeuren, J. Alexandra Hakala, Karl Jarvis, and Johnathan E. Moore Environ. Sci. Technol., Just Accepted Manuscript • DOI: 10.1021/acs.est.7b01979 • Publication Date (Web): 19 Jul 2017 Downloaded from http://pubs.acs.org on July 27, 2017
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.
Environmental Science & Technology 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 35
Environmental Science & Technology
1
Mineral reactions in shale gas reservoirs: barite scale formation from reusing produced water as
2
hydraulic fracturing fluid
3
Amelia N. Paukert Vankeuren* a,b, J. Alexandra Hakala b, Karl Jarvis c,d, Johnathan E. Moore c,d
4 5
a
Geology Department, California State University Sacramento, Sacramento, CA 95819
6
b
National Energy Technology Laboratory, U.S. Department of Energy, Pittsburgh, PA 15236
7
c
National Energy Technology Laboratory, U.S. Department of Energy, Morgantown, WV 26507
8
d
AECOM, Morgantown, WV 26507
9
*Corresponding author contact information: Mailing address: California State University,
10
Sacramento; 6000 J Street, PLR 1016; Sacramento, CA 95819-6043; Phone: (916) 278-7385;
11
Fax: (916) 278-4650; Email:
[email protected] 12 13
Abstract
14
Hydraulic fracturing for gas production is now ubiquitous in shale plays, but relatively little is
15
known about shale-hydraulic fracturing fluid (HFF) reactions within the reservoir. To
16
investigate reactions during the shut-in period of hydraulic fracturing, experiments were
17
conducted flowing different HFFs through fractured Marcellus Shale cores at reservoir
18
temperature and pressure (66oC, 20 MPa) for one week. Results indicate HFFs with
19
hydrochloric acid cause substantial dissolution of carbonate minerals, as expected, increasing
1
ACS Paragon Plus Environment
Environmental Science & Technology
20
effective fracture volume (fracture volume + near-fracture matrix porosity) by 56-65%. HFFs
21
with reused produced water composition cause precipitation of secondary minerals, particularly
22
barite, decreasing effective fracture volume by 1-3%. Barite precipitation occurs despite the
23
presence of antiscalants in experiments with and without shale contact, and is driven in part by
24
addition of dissolved sulfate from the decomposition of persulfate breakers in HFF at reservoir
25
conditions. The overall effect of mineral changes on the reservoir has yet to be quantified, but the
26
significant amount of barite scale formed by HFFs with reused produced water composition
27
could reduce effective fracture volume. Further study is required to extrapolate experimental
28
results to reservoir-scale, and to explore the effect that mineral changes from HFF interaction
29
with shale might have on gas production.
30
1. Introduction
31
Over the past decade, technological innovations in horizontal drilling and hydraulic fracturing
32
have spurred a shale gas boom and made the US the largest natural gas producer globally.1 The
33
Marcellus Shale, which underlies parts of Pennsylvania, New York, West Virginia, and Ohio, is
34
a key contributor to US natural gas production.2 In Pennsylvania alone, over 9,500
35
unconventional wells have been drilled in the past 10 years and annual natural gas production
36
has reached over 4.6 trillion cubic feet.3
37
The dramatic increase in gas production has led to a commensurate increase in produced
38
water – wastewater that is produced at the wellhead as a result of oil and gas production, either
39
through flowback of injected hydraulic fracturing fluid (HFF) or formation water from the
40
reservoir. There are now over 1.4 billion gallons of produced water generated in Pennsylvania
41
each year.4 This water contains high total dissolved solids (TDS) and naturally occurring
42
radioactive materials, making it unsuitable for disposal at wastewater treatment plants.5 Other
2
ACS Paragon Plus Environment
Page 2 of 35
Page 3 of 35
Environmental Science & Technology
43
states, e.g., Texas, dispose of the vast majority of their produced water through deep subsurface
44
injection wells, but Pennsylvania has only 8 active injection wells for the disposal of oil and gas
45
produced water.6-8 Some produced water from Pennsylvania is trucked to disposal wells in West
46
Virginia and Ohio, but that carries high costs and additional risks associated with the transport of
47
large volumes of wastewater.9-10 Due to the limited options for disposal of produced water in
48
Pennsylvania, most produced water is now diluted with freshwater and reused as HFF.9, 11-13
49
Reusing produced water as HFF introduces another complication – the high concentration of
50
alkali earth metals (calcium, barium, strontium) in produced water could lead to scale
51
formation.14-17 Precipitation of scale minerals could result in local reductions in reservoir
52
porosity and fracture aperture, preventing gas from reaching the wellbore.16, 18 The potential for
53
scale minerals to damage well productivity is well-recognized by the industry, and antiscalants
54
are usually included in the HFF chemical mixture.19 However, recent experiments have shown
55
antiscalants are ineffective at preventing precipitation of common scaling minerals (e.g., barite,
56
calcite) under supersaturated conditions.20-22 Given the high concentrations of alkali earth metals
57
in produced water,13, 23 mixing produced water with freshwater and HFF chemicals that generate
58
sulfate may cause supersaturation and precipitation of scaling minerals. Interaction between HFF
59
and shale minerals could also change the fluid composition (e.g., increasing dissolved sulfate
60
through pyrite oxidation) and lead to mineral scaling.
61
The most likely timeframe for scale mineral formation related to HFF is during injection–
62
when the fluid temperature rises from ambient to reservoir temperature– and the shut-in
63
(soaking) period– the days to month timeframe after hydraulic fracturing but before well
64
production when HFF sits in contact with the shale with little to no flow. What happens to the
65
fluid during the shut-in period is not well understood, and the relationship between well
3
ACS Paragon Plus Environment
Environmental Science & Technology
66
productivity and length of the shut-in period is debated.24-25 While fluid imbibition into shale
67
during the shut-in period has been examined as a mechanism for reduced well productivity,26
68
geochemical processes such as mineral reactions within the shale reservoir could also impact gas
69
production.
70
Though common mineral reactions have been well characterized in the lab, these
71
experiments are typically done in isolation with an emphasis on a particular reaction (e.g., calcite
72
dissolution or barite precipitation). Interactions between the complex mixture of solutes and
73
chemical additives present in HFF may alter mineral behavior and yield results that do not match
74
those predicted from laboratory studies. Thus it is important to conduct experiments in the lab
75
that use the complex mixtures employed in the field.
76
Prior experiments suggest that mineral reactions due to shale-fluid interaction could be
77
significant. Studies evaluating reactions between water with and without HFF chemicals and
78
various shales at ambient pressure and temperature have shown evidence of dissolution of
79
carbonate minerals and iron and manganese oxides, and precipitation of sulfate minerals.27-29 At
80
elevated temperature (80oC), under oxidizing conditions such as those created by the exposure to
81
the atmosphere and/or inclusion of an oxidant breaker in the HFF chemicals (typically
82
ammonium persulfate)30 these experiments showed evidence of pyrite oxidation, and
83
precipitation of iron oxides and sometimes barite.27, 29, 31 High temperature, high pressure rocking
84
autoclave experiments reacting Marcellus Shale with synthetic shale formation brine mixed with
85
HFF chemicals or synthetic diluted produced water also exhibited carbonate mineral dissolution
86
and minor anhydrite or gypsum precipitation, and predicted precipitation of barite and nontronite
87
clay.32-33 The formation of barite scale in relation to shale gas wells has been investigated in
4
ACS Paragon Plus Environment
Page 4 of 35
Page 5 of 35
Environmental Science & Technology
88
experiments performed at ambient pressure and temperature, including interactions with
89
antiscalants and using abandoned mine drainage to remove barium from HFF.34-35
90
All prior experiments have utilized shale powder and/or rock chips in benchtop or batch
91
reactor settings, precluding evaluation of spatial differences along a flow path and resulting in a
92
fluid:rock ratio much higher than expected in shale reservoirs. Most experiments were conducted
93
at ambient temperature and pressure, far from the conditions found in shale reservoirs. While
94
these experiments provide valuable preliminary data on mineral reactions in shale reservoirs,
95
there is a need for fractured core flow experiments at reservoir temperature and pressure to more
96
closely approximate conditions in the field.
97
To explore mineral reactions within the shale reservoir during the shut-in period and related
98
effects on gas production, fractured core flooding experiments were conducted exposing shale to
99
different compositions of HFF at reservoir temperature and pressure conditions. This study aims
100
to investigate:
101
•
What mineral reactions result from HFF injection into a shale reservoir?
102
•
What mineral reactions result from varying HFF composition?
103
•
What effect do these reactions have on fracture volume, and how might that affect well
104
productivity?
105
2. Experimental method
106
Experiments were conducted exposing Marcellus Shale cores to synthetic HFFs of varying
107
compositions for 7 days. Changes to fluid chemistry and shale fracture surfaces were used to
108
determine dissolution and secondary mineral precipitation reactions resulting from fluid-rock
109
interaction. Experiments were designed to closely approximate conditions in the field during
110
hydraulic fracturing. Pressure and temperature were set to those expected at a depth of 6,600 ft,
5
ACS Paragon Plus Environment
Environmental Science & Technology
111
where the Marcellus Shale occurs in southwestern Pennsylvania (Table 1). Oxygen was limited
112
in fluids by degassing fluids in a vacuum chamber and then flushing the headspace above HFF
113
with nitrogen gas to prevent addition of dissolved oxygen, but oxygen was not entirely removed
114
from the fluid prior to injection because the vast majority of operators do not incorporate oxygen
115
scavengers in their HFF.19 Because these experiments are intended to simulate the shut-in period,
116
the flow rate was set as slow as possible while still generating sufficient effluent for analyses.
117
The residence time of fluid in these experiments is approximately 100 minutes. While this is
118
likely shorter than would be expected during the shut-in period, when neither injection nor
119
production is occurring, it is a reasonable approximation of prolonged contact between shale and
120
HFF. For comparison, during hydraulic fracturing our industry collaborator reported average
121
fluid injection rates of 12,000-17,500 liters/min, and residence time is estimated at 0.5-0.8
122
minutes.36
123
Table 1: Experimental parameters Temperature 65.5oC Pore pressure
20 MPa
Confining pressure
21.4 MPa
Headspace
High purity N2 gas
Fluids tested
1) SW: local spring water (Table S1) 2) SWF: local spring water with HFF chemicals added (Table S2) 2) PW: synthetic reused produced water (Table S3) 3) PWF: synthetic reused produced water with HFF chemicals added 4) PWFNA: synthetic reused produced water with HFF chemicals added, excluding hydrochloric acid (HCl)
6
ACS Paragon Plus Environment
Page 6 of 35
Page 7 of 35
Environmental Science & Technology
Fluid flow rate
0.04 ml/min
Dimensions of material
Length 152.4 mm
in core holder
Diameter 38.1 mm Shale experiments
Control experiments
Duration
7 days
2 days
Fluid sampling interval
2 days, 7 days
2 days
Material in core holder a
Marcellus shale: 46.5% clay (32.5% illite, 7%
316 stainless
illite-smectite, 7% chlorite), 23.3% calcite,
steel cylindrical
18.6% quartz, 4.6% pyrite, 7% organic content
spacers
174 x 103 mm3
N/A
Shale volume
Average fracture volume 3927 mm3
124
Estimated fracture
>210 µm based on minimum proppant grain size, N/A
aperture b
500-1,500 µm based on x-ray CT imaging
Proppant composition c
40/70 mesh northern white sand (>99% SiO2)
a
b
127
Average shale mineralogical composition in percent mineral content from 37, average shale
It was not possible to measure fracture aperture at experimental pressure. Aperture is estimated from x-ray CT images taken at ambient pressure.
128
c
129
2.1 Fluid compositions
130 131
N/A
organic content from 38
125 126
N/A
Unimin Energy Solutions, The Woodlands, TX
All five experimental fluids (Table 1) used local spring water as the base fluid to ensure a common starting point and to match the fresh water used in industrial hydraulic fracturing
7
ACS Paragon Plus Environment
Environmental Science & Technology
132
operations. Synthetic reused produced water was generated by adding salts to spring water to
133
simulate the composition of HFF using diluted reused produced water from a well drilled in the
134
Marcellus Shale in Greene County, PA that was provided by an industry collaborator (Table
135
S3).39 The compositions of HFF from other wells targeting the same reservoir have been reported
136
and discussed in a previous publication.40 HFF chemicals were chosen to represent an average
137
HFF composition used in western Pennsylvania and eastern Ohio, collated from numerous
138
reports found on FracFocus.org.41-43 The selected HFF chemicals included a gelling agent, scale
139
inhibitor, corrosion inhibitor, clay stabilizer, iron control, friction reducer, cross linker, breaker,
140
surfactant, pH adjusters, and biocide. Two fluids also included HCl because HCl is usually
141
injected into a well prior to hydraulic fracturing to clean perforations and allow HFF access to
142
the shale (a step known as the acid pack). Because acid injection precedes injection of HFF,
143
PWFNA may be most representative of conditions within the shale reservoir at the time of
144
hydraulic fracturing. However, the acid neutralizing capacity of the shale will depend on the
145
amount of carbonate minerals accessible to dissolve and buffer the pH, will vary by shale
146
formation and by location within the formation. While it is possible that all acid is spent prior to
147
fracturing, it is also possible that some acid remains. Operators generally do not take samples of
148
fluid from the borehole after acid injection but prior to fracturing, so there are no data to
149
constrain the pH of the fluid in the well prior to fracturing. Due to this uncertainty, experiments
150
were run for both end-member cases: with HCl, as if no acid were spent resulting in pH ~2, and
151
without HCl, as if all acid were spent, resulting in pH ~8. The complicated mixture in HFF
152
contains many compounds that could– and in some cases are intended to– affect mineral
153
reactions, making it exceedingly difficult to predict mineral behavior within the reservoir.
154
2.2 Experimental set up
8
ACS Paragon Plus Environment
Page 8 of 35
Page 9 of 35
Environmental Science & Technology
155
A system was constructed to simulate HFF injection into a shale reservoir and subsequent
156
reaction between fluids and shale minerals (Table 1, Fig S1).
157
For each fluid, a control experiment was conducted in which the fluid flowed through a core
158
holder containing stainless steel spacers rather than shale at experimental temperature and
159
pressure. These control experiments allowed the differentiation between reactions that result
160
from heating and pressurizing the HFF and reactions that result from fluid-shale interaction. For
161
shale experiments, effluent was collected in the syringe pump reservoirs and sampled after 2
162
days (2-day samples), representing an average of effluent from 0-2 days, and again after 7 days
163
(7-day samples), representing an average of effluent from 2-7 days. Detailed experimental
164
protocols are provided in the supporting information.
165
2.3 Shale core preparation:
166
Cores were cut from the interior of large Marcellus Shale samples taken from outcrops away
167
from the weathered surface, artificially fractured using the Brazilian method 44 (Hydrasplit
168
Masonry Stone Splitter CM-10, Park Industries, St. Cloud, MN). Fractures were packed with
169
proppant to maintain a propped fracture at experimental pressure and simulate the nature of the
170
fracture once proppant has been emplaced and fluid pressure is reduced. The mineralogy of the
171
shale exposed along each fracture face varies slightly, but is assumed to be representative of an
172
average Marcellus Shale composition. Scanning electron microscopy with energy dispersive x-
173
ray spectroscopy (SEM/EDS; Quanta 600 FEG, FEI, Hillsboro, OR) supports this assumption
174
(Fig. S1).
175
2.4 Analytical methods:
176
Before and after experiments, SEM/EDS analyses were performed on shale fracture faces, and x-
177
ray computed tomography (CT) with 24 µm pixel resolution (M5000 Industrial Computed
9
ACS Paragon Plus Environment
Environmental Science & Technology
178
Tomography System, North Star Imaging, Rogers, MN) was performed on select core sections:
179
the quarter lengths of core nearest the fluid inlet and outlet (hereafter referred to as the inlet
180
38mm and outlet 38mm, respectively) (Fig S3). Images were processed with ImageJ.45
181
For experiments where CT scans showed significant mineral changes at the conclusion of the
182
experiments, analyses included segmentation of the CT images using ilastik, an interactive
183
segmentation toolkit.46 The pixel classification workflow was used to differentiate between
184
specific components: matrix shale, proppant grains, open fracture, mineral dissolution, and
185
secondary mineral precipitation, where applicable.
186
Fluid chemistry was measured for all fluid samples (see supporting information for detailed
187
analytical methods). The complicated composition of HFFs required special treatment to analyze
188
via instruments using inductively coupled plasma, which is discussed in the supporting
189
information. All but a few samples have charge balances within +4% electroneutrality, and all
190
samples are within +20%.
191
3. Results and discussion
192
The most significant mineral changes observed during experiments were calcite dissolution and
193
barite precipitation. Fluid chemistry data and plots of significant changes are provided in the
194
supporting information (Table S5a and S5b, Fig S2). Fluid chemistry data were speciated and
195
saturation indices calculated (SI = log Q/K, where Q is the ion activity product and K is the
196
equilibrium constant) with PHREEQC v3.2 using the LLNL thermodynamic database at 65.5oC
197
and 0.101 MPa (Table S6).47
198
Volume changes due to mineral reactions were calculated when quantifiable at CT scan
199
resolution (i.e., for the inlet 38 mm of core for experiments with HFF chemicals) (Table 2).
200
Spatial variations in volume changes are depicted in 3-D web enhanced objects highlighting
10
ACS Paragon Plus Environment
Page 10 of 35
Page 11 of 35
Environmental Science & Technology
201
initial fracture volume (yellow), matrix porosity added by dissolution (green), effective fracture
202
volume lost by secondary mineral precipitation (white), and quartz proppant (red) (WEO 1-3).
203
Effective fracture volume (defined for this paper as fracture volume + near-fracture matrix
204
porosity) changes are estimates; cores were depressurized and moved to the CT scanner, so shale
205
cores may have shifted slightly in the core sleeves, preventing the direct comparison of core
206
locations before and after the experiments. Additionally, initial fracture volumes are maximum
207
values because CT scans were performed at 0.101 MPa, while experiments were conducted at 20
208
MPa pore pressure, 21.4 MPa confining pressure. Mineral changes in the outlet 38 mm were
209
below pixel resolution and could not be quantified. All x-ray CT images are provided (Fig S3
210
and S4) but only results from the inlet 38 mm of core are discussed.
211
3.1 Experiments with natural spring water
212
Fluid chemistry
213
Both SW and SWF experiments show an increase in pH, Ca2+, and dissolved inorganic
214
carbon (DIC) in the effluent with respect to the influent (Table S5, Figure S2), though changes
215
are more pronounced in SWF experiments due to the much lower starting pH. In 2-day samples,
216
the increase in Ca2+ in the SWF experiments is over 12 times the increase in the SW experiments.
217
Dissolution of calcite likely is responsible for both pH buffering and the release of Ca2+ and DIC
218
from the shale into the fluid. In SWF, the increase in DIC is not stoichiometric with Ca2+
219
because DIC in the sample is likely an underestimate. Samples were exposed to the atmosphere
220
during sampling, and given the low pH of the solutions, much of the DIC should have been in the
221
form of carbonic acid and may have been lost by degassing of CO2. Nonetheless, calcite
222
dissolution is supported by saturation indices (Table S6): SW-influent and SWF-influent are
223
highly undersaturated with respect to calcite (SI of -6 and -11, respectively) and after 7 days of
11
ACS Paragon Plus Environment
Environmental Science & Technology
224
contact with shale, saturation indices in the effluents are much closer to calcite equilibrium (SI of
225
-1 and 1, respectively).
226
Page 12 of 35
Sulfate also increased in 2 day-samples with respect to the influents for both SW and
227
SWF experiments, with the increase 2.5 times larger in SWF experiments (Table S5b). With SW,
228
control experiments show little increase in SO42- compared to the experiments contacting shale.
229
This implies that the increase in SO42- when shale reacts with fresh water is likely due to
230
oxidation of pyrite, something observed in prior studies on HFF-shale reactions, even – to a
231
limited extent – under anoxic conditions.27, 29, 32 In SWF experiments, SO42- concentration
232
increases even in the control experiment due to thermal activation of ammonium persulfate in the
233
HFF chemicals to produce sulfate ions.48 In both the SW and SWF experiments, SO42- in 7-day
234
samples decreased slightly with respect to the 2-day samples. In SW experiments, though
235
gypsum remains undersaturated in fluid samples (SI -1.7 to -2.3), the increased Ca2+ from calcite
236
dissolution and SO42- from pyrite oxidation caused local oversaturation and precipitation of
237
gypsum, thus lowering dissolved SO42- between days 2 and 7. In SWF experiments, the larger
238
increase in SO42- from HFF chemicals may have allowed the formation of both barite and
239
gypsum and corresponding reductions in SO42-. Though one would expect the concentration of
240
sulfate to reach equilibrium with respect to barite and then remain at that concentration, there
241
may have been temporary oversaturation due to the rate of SO42- release by persulfate
242
decomposition exceeding the rate of removal by barite precipitation. Saturation indices for barite
243
are at equilibrium or moderately oversaturated throughout the experiments (SI of 0 to 0.4). Barite
244
precipitation dynamics are discussed in more detail following the experimental results.
245 246
In both SW and SWF experiments, small increases (up to 0.15 mmol/L) in Mg2+, Na+, SiO2, and Cl- are also observed. Na+ and Cl- increases may be due to dissolution of halite or
12
ACS Paragon Plus Environment
Page 13 of 35
Environmental Science & Technology
247
saline pore water from the shale, and Mg2+ and SiO2 from dissolution of dolomite and clay,
248
respectively.
249
Fracture surface imaging
250
SW experiments: SEM analysis reveals dissolution of calcium carbonate minerals near the fluid
251
inlet, and minor precipitation of iron oxides and gypsum near the outlet (Fig S3). Gypsum seems
252
to have preferentially precipitated on top of primary pyrite (Fig S3). Also near the outlet, new
253
Mg-Al-rich minerals appear, which could be due to precipitation or deflocculation and transport
254
of clays from closer to the inlet. SW experiments were the only ones showing new Mg-Al-rich
255
minerals; the low salinity of the spring water may be responsible for deflocculation of clays in
256
the shale.49
257
X-ray CT images do not show visible alteration to the shale (Fig S3) because changes in
258
mineral volumes are below the 24 µm pixel resolution. This prevents calculation of effective
259
fracture volume change using ilastik for these experiments.
260
SWF experiments: X-ray CT images from SWF experiments show a darkened reaction rim of
261
calcite dissolution extending along the surface of the main fracture from the fluid inlet to a
262
distance at least 30 mm along the core and penetrating a maximum of 0.8 mm into the shale (Fig
263
1). SEM images of the fracture surface confirm that extensive dissolution of calcite has occurred
264
in this area. When a dense grouping of proppant is present along part of the fracture, the reaction
265
rim is often narrower near the proppant than elsewhere along the fracture (Fig S3). Clusters of
266
proppant are thought to slow the flow of fluid, forcing most of the fluid through the open fracture.
267
Dissolution is faster along preferential pathways due to the greater availability of low pH
268
reactive fluid.50-52 Additionally, x-ray CT images reveal the presence of a secondary material less
269
dense than the shale, SEM images indicate that this is amorphous and may be residue of HFF
13
ACS Paragon Plus Environment
Environmental Science & Technology
270
chemicals, such as guar gum (Fig S5). This residue is not seen in experiments with reused
271
produced water and HFF chemicals, possibly due to decreased stability of guar gum with the
272
higher concentration of Na+. Studies testing the performance of guar based gelling agents in
273
reused produced water showed significantly decreased viscosity at levels of 85 mmol/L Na+; the
274
concentration of Na+ in HFF with reused produced water composition is 3.5 times that value.53
275 276
Figure 1: SWF experiments. A) and B) x-ray CT of inlet 38mm of shale parallel to flow, B) with
277
calcite dissolution (dark gray) along fracture. The decrease in fracture aperture in B is likely due
278
to the core halves shifting in the core sleeve during the experiment, not precipitation or
279
geomechanical alteration. C) and D) SEM of main fracture face near the fluid inlet, D) with
280
calcite dissolution (dark holes)
281
X-ray CT analysis indicates a significant increase in effective fracture volume (WEO 1).
282
In the inlet 38 mm, calcite dissolution in shale adjacent to the fracture creates matrix porosity,
283
increasing effective fracture volume by 65%. The amorphous secondary material thought to be
284
HFF chemical residue fills 4% of the fracture volume, yielding a net increase of 61%. For the
14
ACS Paragon Plus Environment
Page 14 of 35
Page 15 of 35
Environmental Science & Technology
285
outlet 38 mm, there was no measurable change in effective fracture volume on the x-ray CT, so
286
the impact of calcite dissolution on effective fracture volume appears to only apply to areas that
287
are exposed to HFF with residual acid.
288
3.2 Experiments with synthetic reused produced water
289
Fluid chemistry
290
All fluids with a reused produced water base have high TDS, so subtle changes in solute
291
concentration relative to the influents may not be perceptible. In PW and PWF experiments, pH
292
and DIC both increase due to dissolution of carbonate minerals. PW experiments have only
293
minor dissolution, causing a slight increase in DIC and no measurable change in Ca2+. PWF
294
experiments have a much more dramatic increase in pH, Ca2+, and DIC due to the lower starting
295
pH. Again, DIC measurements in these samples are likely underestimates due exposure to the
296
atmosphere during sampling. In both sets of experiments, PW and PWF, influents are
297
undersaturated with respect to calcite (SI of -0.03 and -13, respectively) and saturation index
298
increases in the 7-day samples (SI of 1 and -0.3, respectively) due to calcite dissolution.
299
PWFNA experiments show a decrease in pH and DIC due to precipitation of calcite.
300
PWFNA starts off oversaturated with respect to calcite (SI of 1.9) and draws closer to
301
equilibrium in the 7-day samples (SI of 1.3). Precipitation of calcite in the control experiment
302
may have occurred due to the increased temperature in the core holder lowering calcite solubility.
303
PW at 20oC is undersaturated with respect to calcite (SI of -0.65), while the same fluid at 66oC is
304
oversaturated (SI of 1.9).
305
All three experiments, PW, PWF, and PWFNA show decreases in Ba2+ in control
306
experiments, suggesting some barite precipitation occurs from the fluid even without interaction
307
with the shale. In experiments containing HFF chemicals (PWF and PWFNA) with or without
15
ACS Paragon Plus Environment
Environmental Science & Technology
308
shale, the magnitude of Ba2+ decrease in 2-day samples relative to the influents is 2-8 times the
309
decrease in the experiment without HFF chemicals. The greater decrease in Ba2+ is due to barite
310
precipitation driven by the addition of SO42- from activation of ammonium persulfate in the HFF
311
chemicals. In PW experiments, though Ba2+ decreases, SO42- increases in experiments with shale
312
cores so it is likely that SO42- removal from the fluid by barite precipitation is more than
313
balanced by addition of SO42- from pyrite oxidation.
314
Fracture surface imaging
315
PW experiments: Changes in fracture volume in x-ray CT images along the main fracture for the
316
PW experiment are below the 24 µm resolution, preventing quantification of fracture volume
317
change (Fig S4). However, a small side fracture does exhibit mineralization (Fig 2). SEM/EDS
318
of that location reveals the minerals filling the fracture to be a combination of barite and gypsum.
319
SEM images of the shale fracture surface near the inlet show dissolution of carbonate minerals,
320
but farther from the inlet, calcite dissolution ceases. SEM images also show precipitated barite
321
and gypsum at the inlet of the core. In some areas, the gypsum contains significant amounts of
322
strontium. These results are similar to previous experiments exposing samples of Marcellus
323
Shale to synthetic reused produced water at reservoir conditions.33
16
ACS Paragon Plus Environment
Page 16 of 35
Page 17 of 35
Environmental Science & Technology
324 325
Figure 2: PW experiments. A) and B) x-ray CT of shale perpendicular to flow, near fluid outlet,
326
B) with small fracture filled by secondary mineral precipitation. C) and D) SEM images of
327
fracture in B), partially filled with barite and gypsum.
328
PWF experiments: X-ray CT imaging and SEM analysis confirm the extensive dissolution of
329
calcite and precipitation of barite suggested by the fluid chemistry. X-ray CT imaging shows a
330
very bright, dense mineral lining the innermost surface of the main fracture from the inlet to a
331
distance of about 3 mm down the core, and SEM/EDS confirms this is a Ba-S mineral (hereafter,
332
Ba-S minerals are assumed to be barite) (Fig 3). X-ray CT also shows a darker, less dense rim
333
between the barite and unaltered shale extending for the inlet 38 mm; SEM/EDS identifies this as
334
calcite dissolution. In some instances, barite crystals grew into the space vacated by the dissolved
335
calcite (Fig 3). SEM/EDS imaging also shows an amorphous Ca-rich deposit that is interspersed
336
with the barite crystals near the fluid inlet. EDS analysis of the deposit rules out gypsum (no S)
337
and calcite is unlikely given the low pH of the fluid, so the exact composition has yet to be
338
identified. At the inlet, where the fluid first meets the core, x-ray CT and SEM/EDS show barite
17
ACS Paragon Plus Environment
Environmental Science & Technology
339
precipitated in the shape of the fluid distributor (Fig 3), suggesting considerable precipitation of
340
barite upon entry into the core holder. Similar to SWF, x-ray CT images show that when dense
341
groupings of proppant are present, mineral reactions (barite precipitation and calcite dissolution)
342
are more pronounced along open areas of the fracture, the preferred fluid pathway.
343
X-ray CT analysis indicates a net effect fracture volume increase of 55% due to extensive
344
dissolution of calcite in the shale matrix along the fracture (WEO 2). Barite precipitation near the
345
fluid inlet produces a slight decrease in effective fracture volume (1%), but that was more than
346
compensated for by calcite dissolution.
347
PWFNA experiments: There is a visible white coating on the shale fracture surface extending the
348
length of the shale core, though it is concentrated near the fluid inlet (Fig S6). SEM/EDS
349
confirms this is barite. X-ray CT scans show that in addition to precipitating along the main
350
fracture, secondary minerals (presumably barite) filled portions of smaller side fractures.
351
SEM/EDS and x-ray CT show barite precipitated extensively near the fluid inlet, again in the
352
shape of the flow distributor, and filled small fractures at the inlet (Fig 3). Barite also cemented
353
some proppant grains to the fracture face.
354
X-ray CT analysis for these experiments exhibited a net decrease in effective fracture volume
355
(WEO 3). The lack of calcite dissolution meant there was no mechanism for increasing effective
356
fracture volume, while barite precipitation decreased fracture volume by 2%.
357
The barite precipitated in this experiment exhibits many different morphologies (Fig S6). At the
358
inlet of the core, barite forms in both rounded and sharp tabular shapes, sometimes clustering
359
together in masses and sometimes as isolated crystals. In the inlet 38 mm, some of the barite has
360
a more porous look, while other barite crystals are much smaller and more angular. The reasons
361
for this are discussed in the next section.
18
ACS Paragon Plus Environment
Page 18 of 35
Page 19 of 35
362
Environmental Science & Technology
Table 2: Summary of experimental findings for the five different fluids tested. Experimental ID Fluid composition
SW
PW2
SWF2
PWF2
PWFNA2
Spring water
Diluted
Spring water
Diluted
Diluted
produced
with HFF
produced water
produced
water
chemicals
with HFF
water with
chemicals and
HFF
HCl
chemicals
Influent pH
5.1
7.2
1.7
2.0
8.2
Ending pH
7.0
7.9
7.2
6.0
7.4
Starting SI Calcite
-5.9
0.0
-12.8
-10.6
1.9
Ending SI Calcite
-1.0
0.9
1.0
-0.3
1.2
Starting SI Barite
-1.0
0.0
0.4
1.6
1.9
Ending SI Barite
-0.9
0.4
0.4
0.6
1.0
Starting SI Gypsum
-4.1
-3.4
-4.0
-1.8
-1.4
Ending SI Gypsum
-2.3
-2.9
-0.8
-1.9
-1.9
Changes apparent
Calcite
Minor
Calcite
Calcite
Barite and
dissolution,
calcite
dissolution,
dissolution,
minor calcite
gypsum
dissolution,
HFF
barite and
precipitation
precipitation,
barite and
chemical
amorphous Ca-
clay transport
gypsum
residue,
rich
precipitation
minor barite
precipitation
along shale fracture face
precipitation Initial fracture volume (mm3) a
661
1068
1216
Fracture volume lost to secondary mineral
N/A
12
29
19
ACS Paragon Plus Environment
Environmental Science & Technology
Page 20 of 35
precipitation (mm3) Fracture volume lost to secondary amorphous material (HFF chemical residue) (mm3)
24
N/A
N/A
(mm3)
430
601
N/A
Net effective fracture volume change (mm3)
406
589
-29
Net effective fracture volume change (%)
61%
55%
-2%
Matrix porosity added by mineral dissolution
363
a
364
volume changes were quantifiable at x-ray CT resolution (24 µm).
Volume changes were only calculated for the first 38mm of core from the fluid inlet, where
365
20
ACS Paragon Plus Environment
Page 21 of 35
Environmental Science & Technology
366
Figure 3: A-F) PWF and G-
367
L) PWFNA experiments. A-
368
B), and G-H) X-ray CT of
369
inlet 38mm of shale parallel
370
to flow. B) and H) bright
371
spots indicate barite
372
precipitation, and B)
373
darkening indicates calcite
374
dissolution along fracture.
375
C-D) and I-J) X-ray CT of
376
shale perpendicular to flow
377
near fluid inlet. D) and I)
378
with barite densely
379
precipitated in the shape of
380
the flow distributor (M). D) shows pitting from calcite dissolution (3 long dark lines are for core orientation). E) SEM of the fracture
381
parallel to flow with secondary barite in void space created by calcite vein dissolution, and F) close-up of minerals in E). K) and L)
382
SEM of secondary barite filling fractures in J).
21
ACS Paragon Plus Environment
Environmental Science & Technology
383
3.3 Controls on barite precipitation:
384
Barite precipitation occurred in all experiments with synthetic reused produced water. In
385
experiments with produced water and HFF chemicals, dense barite precipitation occurred near
386
the core inlet. At salinities up to those found in these experiments, barite solubility increases with
387
pressure and temperature up to about 150oC, so as the fluid enters the core holder, the solubility
388
should increase.54 There are several potentially complimentary hypotheses for elevated barite
389
precipitation near the fluid inlet: the thermal activation of persulfate in HFF chemicals may add
390
SO42- and increase saturation index with respect to barite, scale inhibitor performance may
391
degrade due to increased temperature and low pH, and the temperature increase may affect
392
nucleation and precipitation kinetics.
393
Effect of temperature
394
Scale inhibitors are less effective at preventing barite precipitation at higher temperature.
395
Previous investigations evaluated barite precipitation from supersaturated solutions with scale
396
inhibitors and found that induction time (time until barite begins to nucleate) in the presence of
397
phosphonate and polycarboxonate scale inhibitors is 10-100 times shorter at temperatures of 50-
398
70oC than at 25oC.55 That temperature swing is similar to the one in this study: when fluid enters
399
the heated core holder the temperature increases from ambient (20oC) to reservoir (65.5oC).
400
Though these experiments utilized a different scale inhibitor (ethylene glycol), it appears the
401
efficacy may have been similarly reduced by increased fluid temperature. Heating of the HFF
402
chemicals also increased dissolved SO42-. Persulfates, which are often included in HFF to
403
decompose gelling agents and decrease fluid viscosity, are activated by heat (50-70oC), acidic
404
pH, and transition metals – all of which may be present during hydraulic fracturing.48, 56 When
405
activated, the persulfate (S2O82-) degrades into sulfate radicals (SO4-) that then act as oxidizers
22
ACS Paragon Plus Environment
Page 22 of 35
Page 23 of 35
Environmental Science & Technology
406
and produce sulfate (SO42-) as a byproduct of the oxidation reaction. The addition of dissolved
407
sulfate from HFF chemicals is supported by the fact that SO42- increases 1.25 mmol/L in SWF
408
control experiments. The degradation of ammonium persulfate in SWF-influent (0.02 wt %)
409
could supply up to 1.8 mmol/L SO42-.
410
Further evidence of SO42- release from persulfate activation is found in the PWF and
411
PWFNA control experiments. In produced water experiments, the presence of Ba2+ and
412
consequent ability to remove SO42- from solution through barite precipitation makes it
413
impossible to directly measure SO42- increase. However, Ba2+ can be used as a proxy for SO42-.
414
In the PWF control experiment, dissolved Ba2+ decreased by 1.4 mmol/L. If this decrease is
415
driven solely by barite precipitation – which is likely since the low pH precludes formation of
416
barium carbonates – it requires a stoichiometric SO42- decrease. However, PWF influent only
417
contains 0.5 mmol/L SO42-. Given that there was no shale core present in the control, the 0.9
418
mmol/L SO42- difference between what was present in the influent and what was necessary for
419
barite precipitation must have come from HFF chemical decomposition, likely the ammonium
420
persulfate. The SO42- addition could be critical for barite precipitation in the presence of
421
antiscalants for PWF and PWFNA experiments. Influents for these experiments are already
422
supersaturated with respect to barite (SI of 1.6 and 1.9, respectively), and the addition of up to
423
1.8 mmol/L SO42- would increase supersaturation (SI of 2.2 and 2.4, respectively) to above the
424
threshold at which barite precipitates regardless of the presence of antiscalants.22, 57
425
Finally, barite precipitation near the inlet may be driven by the increase in temperature and
426
corresponding increase in barite nucleation and precipitation kinetics.58-59 Though the
427
temperature increase almost doubles barite solubility, from 0.04 to 0.075 mmol BaSO4/kg
428
water,54 the lowest concentration of Ba2+ and SO42- in all of the experiments with produced water
23
ACS Paragon Plus Environment
Environmental Science & Technology
429
(0.18 mmol/kg) remains greater than twice the solubility at 65.5oC. Thus, the increased
430
solubility at higher temperature should not impede barite precipitation. Indeed, increased barite
431
formation at the core inlets in PWF and PWFNA experiments is starkly illustrated by Fig 3. X-
432
ray CT images show barite precipitated in the shape of arcs and triangles matching the shape of
433
the flow distributors. The fluid does not even enter the shale core before barite precipitation
434
begins.
435
Effect of pH
436
The efficacy of scale inhibitors at preventing barite precipitation varies not only with
437
temperature, but also pH. In PWF experiments, barite precipitation was strongly concentrated at
438
the inlet of the core. In PWFNA experiments, while there is a concentrated barite deposit at the
439
core inlet, a thinner layer of barite is distributed throughout the length of the core. The
440
concentration of barite formation near the fluid inlet in the lower pH PWF experiments parallels
441
previous studies on the pH dependence of barite inhibition by phosphonate and polycarboxylate
442
antiscalants; in those experiments at pH