Professional Documents
Culture Documents
Fawole Etal Manuscript Revision Clean
Fawole Etal Manuscript Revision Clean
DOI:
10.1016/j.envpol.2016.05.075
License:
Creative Commons: Attribution-NonCommercial-NoDerivs (CC BY-NC-ND)
Document Version
Peer reviewed version
General rights
Unless a licence is specified above, all rights (including copyright and moral rights) in this document are retained by the authors and/or the
copyright holders. The express permission of the copyright holder must be obtained for any use of this material other than for purposes
permitted by law.
•Users may freely distribute the URL that is used to identify this publication.
•Users may download and/or print one copy of the publication from the University of Birmingham research portal for the purpose of private
study or non-commercial research.
•User may use extracts from the document in line with the concept of ‘fair dealing’ under the Copyright, Designs and Patents Act 1988 (?)
•Users may not further distribute the material nor use it for the purposes of commercial gain.
Where a licence is displayed above, please note the terms and conditions of the licence govern your use of this document.
If you believe that this is the case for this document, please contact [email protected] providing details and we will remove access to
the work immediately and investigate.
2 1 Introduction..................................................................................................... 5
3 1.1 Motivation ................................................................................................................... 5
4 1.2 Review summary ......................................................................................................... 7
5 1.3 Global oil and gas reserves.......................................................................................... 8
6 1.4 Temporal and geographical trends in gas flaring ........................................................ 9
7 2 Oil and Gas Production ................................................................................. 15
8 2.1 Exploration and exploitation ..................................................................................... 15
9 2.2 Composition of natural gas ....................................................................................... 16
10 2.3 Thermodynamic properties of varying fuel gas compositions .................................. 18
11 3 Overview of the gas flaring processes .......................................................... 19
12 3.1 Gas flaring emission and its environmental impact .................................................. 20
13 3.2 Emissions measurements around real-world gas flaring sites ................................... 23
14 3.3 Types of Flares .......................................................................................................... 26
15 3.3.1 Steam- Assisted Flare.................................................................................................... 27
16 3.3.2 Air-Assisted Flare ......................................................................................................... 28
17 3.3.3 Pressure-Assisted Flare ................................................................................................. 29
18 3.3.4 Non-Assisted Flare........................................................................................................ 29
30
31
1
32 Gas flaring and resultant air pollution: A review focusing on Black Carbon
1
34 School of Geography, Earth and Environmental Sciences, University of Birmingham,
35 Edgbaston, Birmingham, B15 2TT, UK
2
36 Birmingham Institute of Forest Research, University of Birmingham, Edgbaston,
37 Birmingham, B15 2TT, UK
3
38 Environmental Research Laboratory, Department of Physics, Obafemi Awolowo University,
39 Ile-Ife, Nigeria 220005
41 Abstract
42 Gas flaring is a prominent source of VOCs, CO, CO2, SO2, PAH, NOX and soot (black
43 carbon), all of which are important pollutants which interact, directly and indirectly, in the
44 Earth’s climatic processes. Globally, over 130 billion cubic metres of gas are flared annually.
45 We review the contribution of gas flaring to air pollution on local, regional and global scales,
46 with special emphasis on black carbon (BC, "soot"). The temporal and spatial characteristics
47 of gas flaring distinguishes it from mobile combustion sources (transport), while the open-
48 flame nature of gas flaring distinguishes it from industrial point-sources; the high
49 temperature, flame control, and spatial compactness distinguishes gas flaring from both
50 biomass burning and domestic fuel-use. All of these distinguishing factors influence the
52 consider this source separately in emissions inventories and environmental field studies.
53 Estimate of the yield of pollutants from gas flaring have, to date, paid little or no attention to
54 the emission of BC with the assumption often being made that flaring produces a smokeless
55 flame. In gas flares, soot yield is known to depend on a number of factors, and there is a need
2
56 to develop emission estimates and modelling frameworks that take these factors into
57 consideration. Hence, emission inventories, especially of the soot yield from gas flaring
58 should give adequate consideration to the variation of fuel gas composition, and to
59 combustion characteristics, which are strong determinants of the nature and quantity of
60 pollutants emitted. The buoyant nature of gas flaring plume, often at temperatures in the
61 range of 2000 K, coupled with the height of the stack enables some of the pollutants to escape
62 further into the free troposphere aiding their long-range transport, which is often not well-
64
65 Capsule Abstract
66 The review identified gaps in the estimation of emissions from the gas flaring process and
67 argues for explicit recognition of gas flaring emissions in emission inventories and global
68 models.
69
70 Keywords
71 Gas flaring; Oil and gas; Black Carbon; Incomplete combustion; Hydrocarbon; Emission
72 factors
73
74
75
76
77
78
3
79 Variables used and their meanings
82 ρe, ρ∞ the fuel gas and ambient densities respectively (kg m-3)
84 ΔTf mean peak flame temperature rise, K (taken as the difference between the adiabatic
91 ɸ equivalence ratio
92 T∞ ambient temperature
93
94
95
96
97
98
4
99 1 Introduction
102 living. Different forms of energy are in use and new ones are developed in order to meet the
103 increasing needs of society (MacKay, 2008). This quest for dependable, affordable and
104 environmentally benign energy sources has occurred throughout human history; for the last
105 century or so, crude oil has been the focus of world energy. In 2011, crude oil was estimated
106 to provide 52.8 % of the world’s energy (13113 Mtoe); with oil and natural gas accounting
108 Human reliance on oil and gas as an energy source is not without its attendant impact on the
109 environment. During production, detrimental impacts on the environment (air, water and soil)
110 include: oil spills and leakages; venting; sludge disposal; and flaring (Almanza et al., 2012;
111 Osuji and Adesiyan, 2005; Osuji and Onojake, 2004; Sonibare et al., 2010). The post-
112 production impact of oil and gas on the environment is also a major source of concern, but is
114 Gas flaring is often a routine daily activity in oil fields around the world, particularly in oil-
115 rich regions of the world where the infrastructure to capture, store and utilise the gas
116 produced is not available. Flaring is most often associated with Nigeria and the Russian
117 Arctic, but it does still occur in more developed economies: the North Dakota Bakken shale
119 accessed on 29 January 2016). Despite several calls by international bodies such as the World
120 Bank’s Global Gas Flaring Reduction (GGFR) initiative, the volume of gas flared globally
121 appears to have plateaued at around 130 billion cubic metres (bcm) since 2008, or may even
122 have increased (see section 1.4, below). According to GGFR (2013), there was an increase of
123 2 bcm in the volume of gas flared in 2011 compared to the previous year. The latest in this
5
124 series of initiatives by the World Bank - Zero Routine Flaring by 2030 - was launched by the
125 Secretary-General of the United Nations, Ban Ki-Moon and World Bank’s President, Jim
128 Gas flaring, a prominent source of black carbon (BC) and other pollutants, has been ignored
129 or underestimated in emission inventories and models, as a result of which models are
130 struggling to predict measurements of BC in regions of intense gas flaring. The intensity of
131 gas flaring and specifics of atmospheric transport can combine to enhance the role of gas
132 flaring emissions over very large areas (e.g., the Arctic) (Stohl et al., 2013). Presently,
133 treatment of emissions from gas flaring is rather rudimentary in most global emission
134 inventories. As at 2015, to the best of our knowledge, only two global pollutants inventories
135 explicitly accounted for emissions from gas flaring (see section 4.5).
136 Gas flaring is classified as a miscellaneous BC-rich source, a group which includes aviation
137 and shipping which together contribute about 9 % to global BC emission (Bond et al., 2013).
138 Gas flaring is a very different type of combustion compared to other BC sources in this
140 (see below) resulting in flames that can be 8-10 metres in length, with flame temperature as
141 high as 2000K. Gas flaring is a year-round activity in most of the intensive flaring regions of
142 the world, and so differs from transport sources in being stationary and differs from other
143 stationary sources (e.g. cooking and biomass burning) by being relatively constant in time.
144 The GAINS (Greenhouse gas Air pollution INteractions and Synergies) model estimated that
145 gas flaring contributes about 4 % of total anthropogenic BC emissions, majority of which are
146 from the leading gas-flaring nations; Russia, Nigeria, and countries in the Middle East (Bond
147 et al., 2013; Stohl et al., 2013). Gas flaring is estimated to contribute 260 Gg to global BC
148 estimates annually (Bond et al., 2013), while in Russia, it is estimated to have the largest
6
149 contribution of 36.2 % to anthropogenic BC emission (Huang et al., 2015). From a 3-year
150 model simulation, more than 40% of annual mean BC near the surface in the Arctic is
152 Considering the small number of nations that still flare gas, a contribution of 4 % to global
154 scale, the contribution of this ‘overlooked’ source of ambient aerosol loading is likely to be
155 significant.
158 contributions to ambient aerosol loadings. This review is the first to cover virtually every
159 aspects of gas flaring (process, trends, chemistry, flame dynamics and environmental impact).
160 The review starts with a brief assessment of the level and distribution of oil and gas reserves
161 around continents of the world. Compositional variation of the natural gas flared and their
162 thermodynamic properties are discussed because they are likely to have significant impact on
163 pollutant emission rates and overall amount. Next, we provide a discussion of the temporal
164 and geographical trends in gas flaring, including brief comments about how weather
165 conditions in regions where gas flaring is common will impact near-field dispersion and long-
166 range transport. This is followed by a broad discussion of the gas flaring process itself;
167 highlighting how engineering and technology decisions impact on the emission of air
168 pollutants. Finally, several techniques used to estimate emissions from gas flaring are
169 discussed. Throughout this review, special attention is paid to soot (i.e. carbonaceous aerosol
170 predominantly composed of BC) emission from gas flaring because of the now known
171 contribution of BC to global warming and the apparent neglect of the contribution of gas
172 flaring to ambient air aerosol loadings in inventories and global models.
7
173 1.3 Global oil and gas reserves
174 World-proven natural gas and oil reserves at the end of 2012 stood at 187.3 and 1668.9
175 trillion cubic metres (tcm) respectively, sufficient to meet 55.7 years of global production
176 (BP, 2013). The distribution of these reserves among regions of the world is shown in Figures
177 1(a) and (b). The natural gas collected during the exploration of crude oil from the Earth’s
178 crust can be a very good source of fuel; transported in pipes for industrial or domestic use and
179 also recycled back into the processing operation (Davoudi et al., 2013; Elvidge et al., 2009).
180
181
182 Figure 1: (a) The distribution of world natural gas reserves (b): The distribution of world oil
183 reserves (adapted from BP (2013).
184
185 In developing countries and oil-rich regions where the technology, infrastructure, and market
186 structure to put all of the natural gas to meaningful use are not available or inadequate, the
187 excess natural gas then becomes a waste stream and is flared or vented. Gas flaring has been
188 termed ‘gross waste’ by the World Bank’s initiative against gas flaring (Global gas flaring
189 reduction: GGFR), because flaring represents direct injection of fossil carbon into the
190 atmosphere without capture and utilization of the heat produced by combustion.
8
191 1.4 Temporal and geographical trends in gas flaring
192 The all-time peak of volume of gas flared, 172 bcm, was in 2005 (Elvidge et al., 2009).
193 Figure 2(a) shows the quantity of natural gas produced by the top 10 oil producing nations of
194 the world between 2000 and 2011, and for comparison, Figure 2(b) shows the estimated
195 quantity of gas flared by these major oil producing nations during the same period of time.
196 The estimated volume of gas flared globally in 2012 is also shown in Figure 2(b): the 2012
197 data is obtained from Elvidge et al. (2015). Data for the plots in figure 2(a) and (b) are from
198 IEA (2012) and GGFR (2012), respectively. Between 2008 and 2011, there was no
199 significant decrease in the amount of natural gas flared. In fact, there was a slight increase
200 between 2010 and 2011. According to (BP, 2015), the production excludes the quantity flared
201 or recycled. Therefore, total fossil fuel extracted and ultimately released to the environment is
203 Figure 3 gives the temporal variation of the fraction of total gas extracted that is flared
204 between 2000 and 2011. Nigeria, Libya and Kazakhstan still flare a sizeable amount of their
205 total production. As at 2011, Nigeria and Libya flared about 30 % of their total gas extracted
206 while Kazakhstan flares about 15 %. The 30 % (15.2 bcm) flared by Nigeria is more than
207 twice Libya’s total gas extracted (6.3 bcm) for the same period. Russia and Nigeria together
208 account for about 35 % of the gas flared globally (Elvidge et al., 2009). Although the fraction
209 of gas flared is decreasing for most countries and for the largest emitters, several countries
210 show flat or increasing fractions of gas flared. The estimated quantity of natural gas flared in
211 the US and Canada, as shown in Figure 2b, are just for the off-shore flaring which is
212 presumably responsible for the very low quantities recorded. After 2006, there is a factor of
213 22 and 10 increase in flaring from USA and Canada, respectively, resulting, presumably,
9
215
216
217
218 Figure 2: Trend of natural gas (a) production and (b) flaring in major oil producing nations
219 between 2000 and 2011 (adapted from IEA (2012) and (GGFR, 2012) respectively).
10
220
221
222 Figure 3: Trend of quantity flared compared to total production (data from Figure 2)
223
224
225 Table 1 gives the estimates, from satellite data, of the quantity of gas flaring by the top 20 gas
226 flaring nations of the world in 2008 and a brief summary of the climatic conditions of the
227 region based on the Köppen climate classification (see, e.g., Holden (2005))
228 Figure 4, reproduced, with permission, from Casadio et al. (2012), shows the global
229 geographic distribution of gas flaring sites obtained by remote sensing techniques. In Africa,
230 the flaring spots are clustered around the North - Algeria, Libya and Egypt - as well as around
231 the Atlantic Ocean especially the Niger Delta area of Nigeria. In Europe, they are around
232 Russia and Kazakhstan, and Iran, Iraq, Kuwait, and Saudi Arabia in the Middle East. Satellite
11
233 remote sensing not only gives a general picture of the spatial distribution of flares, but also a
235 Table 1: Year 2008 estimated gas flared by top 20 gas flaring nations with an estimated error
236 of ±2.11bcm.
Rank Country Gas Flared Regional Comments*
(bcm) Climate
1 Russia 40.5 (Dfa, Dfb, peri-Arctic emissions; pole-ward
Dwa, Dwb), atmospheric flow around Tibetan
(Dfc,Dfd, anticyclone in northern hemisphere
Dwc, Dwd) winter
and ET
2 Nigeria 15.1 Am and Aw Tropical monsoon and trade-wind
littoral and tropical wet and dry
climate resulting from West African
monsoon winds that changes
direction with season; climatic
condition is controlled by trade
wind and movement of the ITCZ
3 Iran 10.4 Bwh Tropical and subtropical desert
climate characterised by large
diurnal temperature range. Deep
turbulent boundary layer during the
day; shallow stable boundary layer
at night. Large-scale subsidence
(descending branch of the Hadley
circulation) above the boundary
layer.
4 Iraq 7.0 Bsh Mid-latitude steppe and desert
climate characterised by semiarid
annual rainfall distribution
5 Algeria 5.5 Bwh As for Iran
6 Kazakhstan 5.2 Bwk Mid-latitude arid wet and dry
climate
7 Libya 3.8 Bwh As for Iran
8 Saudi 3.5 Bwh As for Iran
Arabia
9 Angola 3.1 Cwa Humid subtropical climate;
equatorward and poleward
circulation during winter and
summer respectively cause changes
in the movement of air masses from
the cold polar and warm tropics
within this climate.
10 Qatar 3.0 Bwh As for Iran
11 Uzbekistan 2.7 Csa Mediterranean climate; it is
controlled by the variation between
subtropical high in summer and
12
polar westerlies in winter
12 Mexico 2.6 Af, Aw, Bsh, Complex climatic condition; two
Bsk, Bwh, tropical, two dry and three
Cwa, Cwb, temperate climates
Cfa, Cfb
13 Venezuela 2.6 Aw Tropical wet-dry climate
14 Indonesia 2.3 Aw, Am As for Nigeria (see comment above)
15 USA 2.3 Bsk, Bsh, Mediterranean/dry summer
Bwk, Csa, subtropical climate characterised by
Csb, Af moist winter and hot dry summer;
subtropical anticyclones are key
factor that control this climate
16 China 2.3 Cfa, Cwa Humid subtropical climate; during
winter the climate is influenced by
the Siberian cold and during
summer there is an inflow warm of
air from the west
17 Oman 1.9 Bwh As for Iran
18 Malaysia 1.9 Af Tropical rainforest characterised by
constant high temperatures and
evenly distributed high
precipitation; controlled by
movement of the ITCZ and rising
air along trade wind coast. Strongly
affected by el Nino Southern
Oscillation
19 Canada 1.8 Dfb, Dfd, Dsc, As for Russia (see above)
Af
20 Kuwait 1.8 Bwh As for Iran
TOTAL 119.3
237 Source: Quantity of gas flared is obtained from Elvidge et al. (2009). *See Holden (2005)
238
239
240 Elvidge et al. (2015), using data collected by National Aeronautics and Space
242 Radiometer Suite (VIIRS), identified more than 7000 active flare globally in 2012: the bulk
243 of which were found in the upstream sector of the oil and gas industries. In 2012, the
244 estimated volume of gas flared globally is 143±13.6 bcm (Elvidge et al., 2015). Compared to
245 119.3±2.11 bcm estimated for 2008 (Elvidge et al., 2009), this is an increment of about 20 %
246 in the central tendencies of the estimates for both years and is well outside the combined
247 uncertainty bounds of both estimates. An increment of about 20 % over four years, stands in
248 stark contrast to the decrease anticipated as a result of the World Bank's Global Gas Flaring
13
249 Reduction (GGFR) initiative. A new dimension to the problem is the inclusion of three new
250 countries – India, Egypt and Turkmenistan - in the list of top 20 gas flaring countries.
251 The geographical location of the sources of gas flaring emissions (as monitored from space,
252 see Figure 4), the atmospheric behaviour of emitted pollutants, and the pollutant “matrix”
253 from other sources into which emission are made, determine to a large extent the effect gas
254 flaring emissions will have on regional pollution loadings and on climate. BC emissions in
255 the peri-Arctic and West African Monsoon (WAM) regions have the potential to interact with
256 regional radiative energy budgets and atmospheric circulations, leading to impacts on their
257 respective regional climates. Gas flaring emissions in the tropics, especially the WAM and
258 South Asian Monsoon (SAM), could have significant regional impact as a result of the
260 Monsoon circulations are characterised by large-scale seasonal reversals of wind regimes.
261 Regions often referred to as ‘monsoonal’ include tropical and near-tropical regions which
262 experience a summer rainfall maximum and most of these regions have a double rainfall
263 maximum (Barry and Chorley, 2009). The annual climatic regime over West Africa has many
264 similarities to that of South Asia. Both are characterised by surface air-flow which is
265 determined by the position of the leading edge of a monsoon trough. Winds are south-
266 westerly to the south of the trough and north-easterly to the North. The lack of a large
267 mountain range in the north of West Africa strongly enhances the northward advance of the
268 WAM compared to its South Asian counterpart. The position of the leading edge of the
269 WAM trough may oscillate greatly from day-to-day through several degrees of latitude
270 (Barry and Chorley, 2009). In the WAM, deep convection occurs in organised systems
271 known as Mesoscale Convective Systems (MCS) (Mari et al., 2011; Mathon and Laurent,
272 2001). Deep convection in the tropics associated with the ITCZ is responsible for intense
14
273 mixing, venting of the atmospheric boundary layer, and long-range transport of aerosol (Law
274 et al., 2010; Reeves et al., 2010; Sultan and Janicot, 2003).
275
276
277 Figure 4: Flaring hot spot sites (1991-2009) as monitored from space are indicated as red
278 spots while grey spots represent position of active volcanoes during the same period Casadio
279 et al. (2012).
283 located in reservoirs buried far into the Earth. These locations can be in very remote and
284 inhospitable parts of the Earth and can be either on- or off-shore. The multi-staged process of
285 exploration, exploitation and processing of crude oil in its raw form from the Earth’s crust
286 can be broadly divided into upstream, midstream and downstream. Environmental
287 contaminants are expelled into the ambient environment - soil, air and water - at different
288 stages of the process. Crude oil is found in reservoirs which also contain gas. This gas is
289 known as “associated natural gas”, and it is separated from the crude oil at a Flow Station.
290 Natural gas includes both “non-associated” gas originating from fields producing only
291 gaseous hydrocarbons, and “associated” gas produced in association with crude oil (IEA,
15
292 2012). Natural gas comprises mainly hydrocarbons, predominantly short-chained alkanes. At
293 the separation stage, some of the natural gas is captured for domestic and industrial use while
294 the rest is disposed of, usually by flaring in open flames. Below, we will use the term ‘fuel
295 gas’ to refer to natural gas that is flared. The quantity of contaminants expelled at this stage
296 of processing depends on the technology employed, quantity of crude oil processed, flare
297 geometry and design, prevailing meteorological conditions and the composition of the flared
298 gas (E & P Forum, 1994; Ismail and Umukoro, 2014; Johnson and Coderre, 2011; Obanijesu
299 et al., 2009; Ouf et al., 2008; Sonibare and Akeredolu, 2004; Talebi et al., 2014).
300 Oil and gas are produced in many wells in different parts of the world at rates varying from
301 50 m3 to 700 m3 per day. As a result of the diverse nature of the geological features of the
302 area where these explorations take place, the composition of oil and gas varies from one
306 Fuel gas is a combination of C1 to C7+ hydrocarbons which are predominantly alkanes. A
307 typical fuel gas sample contains CH4, C2H6, C3H8, n-C4H10, i-C4H10, n-C5H12, i-C5H12, C6H14,
308 C7H16, H2S, CO2 and N2, where ‘n’ and ‘i’ stand for ‘normal’ i.e., straight chained, ‘iso’ or
309 branched-chained alkanes, respectively. The separation of gas and liquid is not perfect at the
310 Flow Stations and as such trace amounts of liquid can occur in the gas stream, enhancing the
16
312 Table 2: Composition (in mole %) and Some Properties of Fuel Gas from Field Stations and Literature.
Composition Flow Stations Lab-based
1 2 3 4 5 6 7 8 9 10 A B
CH4 74.3 79.85 56.9 55.5 82.23 78.41 68.14 68.42 72.32 69.58 85.24 74.54
C2H6 14.0 11.54 21.2 18.0 2.38 5.68 14.22 7.65 2.41 0.25 7.06 12.17
C3H8 5.8 2.25 6.0 9.8 4.24 0.23 10.27 11.27 6.24 12.54 3.11 5.37
nC4H10 2.0 2.58 3.7 4.5 0.94 0.70 3.23 4.39 8.12 2.35 1.44 2.49
iC4H10 - 0.14 - - 5.12 4.12 2.38 4.42 5.12 5.12 - -
nC5H12 0.9 3.24 1.6 1.6 2.25 9.12 0.75 0.94 3.14 5.20 - -
iC5H12 - - - - 2.14 0.25 1.01 1.55 2.48 2.54 - -
C6H14 - 0.14 - - 0.25 0.23 - 0.18 0.15 1.97 - -
N2 2.9 0.1 - 0.2 - 0.05 - 0.16 - 0.24 1.24 2.15
CO2 - 0.16 7.1 8.9 0.45 1.21 - 1.02 - 0.21 1.91 3.28
H2 S 0.1 - 3.5 1.5 - - - - 0.02 - - -
*C:H 0.2659 0.2715 0.2569 0.2602 0.2730 0.2751 0.2860 0.2852 0.2893 0.2924 0.2541 0.2570
*Molar mass (g mol-1) 21.4 21.4 25.8 26.9 22.9 24.7 24.3 25.8 27.0 28.6 19.2 21.5
*Molar mass (g mol-1) 20.5 21.3 22.7 22.9 22.7 24.1 24.3 25.3 27.0 28.4 18.0 19.4
without CO2 and N2
*HHV (MJ m-3) 44.8 46.4 47.8 48.7 49.1 52.0 52.2 54.2 57.2 60.3 39.8 42.5
313 2 and 5 – 10 are adapted from Sonibare and Akeredolu (2004), 1, 3 and 4 from Ismail and Umukoro (2014) and lab-based (A&B) from McEwen and Johnson (2012)
314 Compositions 1, 3 and 4 are from Saudi Arabia, Kuwait and Iraq respectively while 2 and 5 - 10 are from different flow stations in Nigeria.
315 C:H = mass-weighted carbon to hydrogen ratio
316 * evaluated in this work
17
317 2.3 Thermodynamic properties of varying fuel gas compositions
318 The higher heating value (HHV) is the total enthalpy of the complete combustion reaction for
319 the gas mixture, plus the heat of condensation of the water produced during combustion;
320 lower heating value (LHV) is the total enthalpy of the complete combustion when water
321 remains present in its gaseous form only (Flagan and Seinfeld, 2012). The Volumetric
322 heating value (VHV), measured in kJ mol-1, usually refer to the HHV unless otherwise stated.
323 In gas flaring, HHV, calculable using data from standard thermodynamic tables, defines the
324 total amount of energy available to provide buoyancy to the flare. Buoyancy is an important
325 parameter in determining the dispersion of flare emissions (Beychok, 1994; Leahey and
326 Davies, 1984). The heat content of a fuel gas depends on its molar mass and by extension the
327 density of the fuel. Figure 5 presents the best fit line relating HHV to the molar mass of the
328 fuel gases from the Flow Stations given in Table 2. This plot shows the extent of dependence
329 of HHV on the molar mass, and hence, density of the fuel gas. These properties of the fuel
330 gas, varies with the fuel composition. Note that the HHV values for the laboratory flares in
331 Table 1 are lower that the HHV for the Flow Station fuel gases.
18
332
333 Figure 5: Higher heating value (HHV) as a function of molar mass of fuel gas for the Flow
334 Station data reported in Table 1.
336 Gas flaring is carried out with the aim to convert its hydrocarbon content, especially methane,
337 to products that are less hazardous to the immediate vicinity of the flare site. Gas flaring is
338 classified as a stationary combustion source for the purpose of air pollution regulation
339 (USEPA, 2008). The combustion process involves the rapid oxidation of the fuel gas with the
340 release of heat, gaseous and particulate pollutants, whose nature and quantity depend on the
341 amount and composition of gas fuel burned, the combustion characteristics as well as the
342 flare geometry and design (Ouf et al., 2008; Torres et al., 2012a).
343 Gas flaring may be categorized as emergency, process and production flaring depending on
344 the basis of the flaring (Johnson and Coderre, 2011). Emergency flaring is unplanned: it is
345 carried out at large facilities for safety purposes for a short duration of time. During
19
346 emergency flaring, a large volume of gas is disposed of quickly and, hence, the flow rate of
347 the fuel gas is very high. Process flaring is an intermittent disposal of unwanted gas that may
348 last for a few hours or a couple of days at often low flow rate. It occurs during well-testing as
349 well as start-up and shutdown of process units. Production flaring may occur continuously
350 for years as long as the oil is being explored and exploited. The flow rate can be very high at
351 particular times especially during the initial development of a gas well (Johnson and Coderre,
352 2011). As a result of the length of time involved, which can be years, and the flow rate of the
353 gas flared, production flaring is the major process of concern for regional and global
357 combustion in an open flame in the open atmosphere, using a burner tip designed specifically
358 for that purpose, in the course of routine oil and gas production operations (McEwen and
359 Johnson, 2012; OGP, 2000; Stone et al., 1992). Gas flaring is a source of greenhouse gases,
360 precursor gases, VOCs, polycyclic aromatic hydrocarbon (PAH) and particulate matter (PM)
361 in the form of soot (Ana et al., 2012; Johnson et al., 2013; McEwen and Johnson, 2012;
362 USEPA, 2011; USEPA, 2012). These pollutants have been identified to have serious
363 detrimental impact on animals, vegetation and human (Burney and Ramanathan, 2014; Dung
364 et al., 2008; Pope III et al., 2002; USEPA, 2010). The soot emitted from the combustion of
365 natural gas is predominantly black carbon (BC) (Johnson et al., 2013; Smith and Chughtai,
366 1995).
367 BC is a strong climate forcer and plays a prominent role in the nature of the Earth’s climate
368 because of its ability to absorb solar radiation, and hence, to result in a changed vertical
369 gradient of warming by incoming solar radiation (Ramana et al., 2010; Ramanathan and
370 Carmichael, 2008). It also affects cloud processes as well as decreasing surface albedo on ice
20
371 and snow causing them to melt faster (IPCC, 2007). It can be short-lived in the atmosphere,
372 as it is removed from the atmosphere by dry and wet deposition (although insoluble in water,
373 it is wettable, particularly when ‘aged’ by atmospheric oxidation). When mixed with other
374 aerosol components in the atmosphere, BC can affect the climate for longer as its residence
375 time is increased in that state (Bond et al., 2013; Flanner et al., 2007; Quinn et al., 2008;
376 Ramana et al., 2010). Within current estimated uncertainties, BC is the second highest
377 contributor to global warming just after CO2 (IPCC, 2007), and also the main light-absorbing
378 component of atmospheric aerosols (Bond et al., 2013; Chung and Seinfeld, 2005; Jacobson,
380 The radiative forcing capacity of a cloud of aerosol particles depends on the ratio of the BC
381 to other components in the cloud (Ramana et al., 2010). The ubiquitous nature of BC coupled
382 with its several effects on the Earth’s climate makes the study of its sources and emission
383 rates important. If adequate measures are put in place to reduce BC emission, short-term
384 reduction of radiative forcing can be achieved in the Arctic and other oil rich regions of the
385 world (Arctic Council, 2013; Feichter and Stier, 2012; Tripathi et al., 2005).
386 Studies have shown that the contributions from oil and gas processing, especially gas flares,
387 to air pollution have been severely underestimated (Edwards et al., 2013; Johansson et al.,
388 2014; Schultz, 2014; Stohl et al., 2013). This is arguably due to the fact that emission factors
389 (EF) often used to estimate the emissions of pollutants from the oil and gas sector, especially
390 gas flaring, are too general and independent of site specifics such as fuel composition and
391 combustion characteristics. There is a pressing research need for more measurements and
392 development of EFs estimates that vary with the fuel gas composition and combustion
393 characteristics.
394 At the local scale, ground-level measurements within one kilometre downwind of the flaring
395 sites have indicated an elevated concentration of particulate and gaseous pollutants (Edwards
21
396 et al., 2013; Obanijesu et al., 2009; Sonibare et al., 2010). Likewise, on a regional scale,
397 flight campaign during ACCESS (Arctic Climate Change Economy and Society) around oil
398 and gas installations in Heidrun, Norway reported elevated concentrations of NO, SO2 and
400 NOX (NO+NO2), SO2, CO2, CO and unburned hydrocarbon are the major pollutant
401 constituents of gas flaring plumes (USEPA, 1995). The low amount of the nitrogen and
402 sulphur content of the fuel gas notwithstanding, the NOX and SO2 emission from gas flaring
403 remain significant because the ambient background level of both pollutants are usually low.
404 Hence, gas flaring can substantially enhance the local concentrations of these pollutants.
405 Complete combustion is often not achieved in most flaring conditions (Leahey et al., 2001).
406 During incomplete combustion, methane (CH4) and other unburned components of the fuel
407 gas are given off and some fuel components are partially oxidized to CO and soot rather than
408 completely oxidised to CO2 (RTI, 2011; Strosher, 2000; Villasenor et al., 2003). Equations 1a
409 and b, respectively, give simplified equations for the 'ideal' complete and partial oxidation of
410 the fuel gas. In equation (1b), 'rC' denotes soot production, ‘vO2’ excess oxygen in the case of
411 fuel-lean combustion and ‘CaHb’ are PAH and other semi-volatile organics resulting from the
412 pyrolysis of hydrocarbons in the fuel gas. The value of m in equation (1b) determines a fuel-
413 lean (over-fired; m > x + y/4) or fuel-rich (under-fired; m < x + y/4) combustion process. The
414 level of combustion (oxidation) of the fuel gas is dependent on a number of factors: the
415 nature of the flame during combustion, the level of mixing of the gas and air in the reaction
416 zone, the amount of oxidant (oxygen) available, the VHV of the fuel gas and the prevailing
417 condition of the ambient wind (Stone et al., 1992). Flare design and geometry are also key
𝑦 𝑦
419 𝐶𝑥 𝐻𝑦 + (𝑥 + 4) 𝑂2 → 𝑥𝐶𝑂2 + 2 𝐻2 𝑂 …………………………………………… (1a)
22
420 𝐶𝑥 𝐻𝑦 + 𝑚 𝑂2 → 𝑝𝐶𝑂2 + 𝑞𝐶𝑂 + 𝑟𝐶 + 𝑣𝑂2 + 𝑓𝐻2 𝑂 + ∑𝑎,𝑏 𝑑𝑎,𝑏 𝐶𝑎 𝐻𝑏 ……….... (1b)
421 For ‘ideal’ complete combustion, the oxidation of the hydrocarbon yields carbon dioxide and
422 water only (equation 1a), while the oxidation of the sulphur (as H2S) and nitrogen content of
423 the fuel gas gives SO2 and NOX respectively. In gas flaring, NOX is produced by thermal
424 cracking of the nitrogen content of the fuel gas and entrained atmospheric nitrogen. The
425 amount of carbon-containing emission (CO2, CO, CaHb, and BC) given off depends on the
426 molar mass of non-CO2 carbon per mole of the fuel gas (cf. Table 2). Smoking of a gas flare
427 does not necessarily imply that the combustion process is highly inefficient; because a small
428 amount of soot can absorb and scatter perceptible amounts of visible light, a flare with a
429 combustion efficiency of 99 % could still smoke visibly (Castineira and Edgar, 2006).
432 the ACCESS aircraft campaign experiment in the Arctic (Norway) and a couple of ground-
433 based in-situ measurements in Nigeria and the US show significant contributions from gas
435 In-situ ground measurements of air pollutants around typical oil and gas facilities where
436 varying degrees of gas flaring take place have been undertaken in the US (Edwards et al.,
437 2013; Edwards et al., 2014; Johansson et al., 2014), Mexico (Villasenor et al., 2003), Norway
438 (Roiger et al., 2015) and, Nigeria (Ana et al., 2012; Nwaichi and Uzazobona, 2011;
439 Obanijesu et al., 2009; Sonibare et al., 2010). Continuous noise levels higher than the WHO
440 limit of 70 dB were also observed around some gas flaring sites in Nigeria (Abdulkareem and
442 A 4-month sampling of three air pollutants (SO2, CO and NO2) around six flow stations in the
443 Niger Delta area of Nigeria was undertaken by Obanijesu et al. (2009) and Sonibare et al.
23
444 (2010). The measurements were made at 60, 200 and 500 m downwind of the flaring sites.
445 Although to a varying degree depending on the capacity of the station, gas flaring is a
446 prominent daily activity within the stations. Mean pollutants measurements around the six
447 flow stations are shown in Figures 6(a), (b) and (c). The variation bar on the bar-plots shows
448 the standard deviation of the measurements over the four-month period. The nature and
449 extent of dispersion of pollutants from a stationary source depend on the local meteorology
450 and topography of the area. As shown in Figures 6(a)-(c), the trend of the measurements from
451 the flaring sites are similar to observations from dispersion model studies where
452 concentrations of pollutants decreases exponentially with distance from the source (Hodgson
453 et al., 2007). Site 4 is the only place where a significant deviation from this trend, especially
454 for CO and NO2, was observed, suggesting the likelihood of contributions from other
455 source(s).
456
457 Figure 6a: Spatial variation of SO2 concentration downwind of six gas flaring sites
458 (adapted from Obanijesu et al. (2009))
459
460
461
24
462
463
464 Figure 6b: Spatial variation of CO concentration downwind of six gas flaring sites (adapted
465 from Sonibare et al. (2010))
466
467 Figure 6c: Spatial variation of NO2 concentration downwind of six gas flaring sites (adapted
468 from Sonibare et al. (2010))
469
470
25
471 A summary of the few available in-situ ground measurements downwind of gas flaring sites
473
474 Table 3: Pollutant measurements around several oil and gas facilities
BC O3 VOC PAH NO NO2 SO2 CO Ref
(ng kg-1) (ppbv) (ppbv) (ng m-3) (ppbv) (ppbv) (ppbv) (ppbv)
- > 120 100 -350 - > 3.5 > 7.5 - > 80 Edwards et
al. (2013)
> 40 > 25 - - > 1.2 - > 1.2 > 90 (Roiger et
al., 2015)
- - - 0.34 - - - - - (Ana et al.,
3.3x104 2012)
475
476
479 assisted, steam-assisted, non-assisted and pressure-assisted. Flares are assisted primarily to
480 enhance the turbulence and mixing of the fuel gas and air in the combustion zone, so as to
481 suppress smoking of the resulting flame (Castineira and Edgar, 2006; Enviroware, 2012;
482 Torres et al., 2012a). The choice of assistance, therefore, affects flame chemistry, as
483 discussed below. Air and pressure-assisted flares are not as efficient as steam-assisted flares
484 in terms of the carbon conversion efficiency (CCE) (Castineira and Edgar, 2006). Complete
485 combustion of the fuel gas requires sufficient air for combustion and adequate mixing of the
486 air and fuel gas. The efficiency of a gas flare at a given moment in time depends on the HHV
487 of fuel gas (see table 2), design of the burner, mixing of air and fuel gas in the combustion
488 zone, composition of the fuel gas, wind speed and direction, and ambient temperature and
489 pressure (Kostiuk et al., 2004; Stone et al., 1992; Torres et al., 2012b).
490 Equivalence ratio, ϕ, is a measure of the amount of oxygen available for the combustion of
491 the fuel gas (Flagan and Seinfeld, 2012; McAllister et al., 2011). It is defined as:
𝐴
( )
𝐹 𝑠
492 ϕ= 𝐴 ………………………………………………………………………. (2)
( )
𝐹 𝑎
26
𝐴 𝐴
493 where (𝐹 ) is the stoichiometric air-fuel mass ratio and (𝐹 ) is the actual air-fuel mass ratio.
𝑠 𝑎
494 ϕ < 1 implies fuel-lean mixture, that is, more oxygen than is needed for the combustion of
495 the fuel is available; ϕ = 1, a stoichiometric mixture, where the exact amount of oxygen
496 needed is made available, and ϕ > 1, a fuel-rich mixture, less oxygen than is needed is
497 available for the combustion of the fuel. Burners used in gas flaring and the entire gas flaring
498 set-up are designed to produce a flame operating at ϕ ~ 1 taking into account the economic
499 aspects such as cost effectiveness of applying a process like steam assistance.
502 flare to provide added momentum and turbulence to the fuel gas and air, enhancing mixing
503 and as such suppressing the tendency for smoke in the flare (Castineira and Edgar, 2006;
504 Müller-Dethlefs and Schlader, 1976; Stone et al., 1992). It is the most efficient assist given to
505 flare to suppress smoking during combustion because the steam affects flame chemistry as
506 well as mixing. The steam acts to break up long-chain hydrocarbons to smaller chains that
507 burn with less smoke (Castineira and Edgar, 2006; Fortner et al., 2012; Müller-Dethlefs and
508 Schlader, 1976; Torres et al., 2012a). Steam undergoes thermal dissociation in a flare flame
509 to give H and OH free radicals that react with carbon to give CH2- and –CHO radical
510 moieties. The steam-induced free radicals enhance the formation of C=O bonds rather than C-
511 C bonds, promoting completeness of combustion. The steam can also react with intermediate
512 products like CO, oxidizing it further to CO2 (Castineira and Edgar, 2006; Müller-Dethlefs
514 CO (g) + H2O (g) CO2 (g) + H2 (g) (ΔH = -41 kJ mol-1) ……… (3)
515 For a steam assisted flare, combustion efficiency starts to decrease when the steam-to-fuel
516 gas ratio goes beyond a threshold which depends on the heat content of the fuel gas and the
517 location at which the steam is injected into the combustion zone in the flare (USEPA, 2014).
27
518 In the US, a lower heating value (LHV) limit of 11.18 MJ/m3 is imposed on fuel gas in order
519 to be suitable for steam-assisted flares (Castineira and Edgar, 2006). Fuel gas with high heat
520 content promotes higher combustion efficiency in steam-assisted flares (McDaniel and
521 Tichenor, 1983; Torres et al., 2012b). Over-steaming results when too much steam is injected
522 into the combustion zone; over-steaming is analogous to over-aeration in air-assisted flares.
523 Over-steaming causes a decrease in the flame temperature by serving as a heat sink. A
524 decrease in combustion efficiency of the flare results, along with, increased noise caused by
525 cavitation created within the flame (Castineira and Edgar, 2006; Stone et al., 1992).
526 The formation of NOX in steam-assisted flares is reduced, compared to non-assisted flares,
527 and further reduced at large values of equivalence ratio (ϕ) due to a drop in the flame
528 temperature (Miyauchi et al., 1981; Müller-Dethlefs and Schlader, 1976). Steam-assisted
529 flares are rather expensive to maintain, especially for large gas facilities, as a large-scale
533 momentum and turbulence to the fuel gas in the combustion zone and, hence, enhances the
534 mixing of the fuel gas and air in the zone (Castineira and Edgar, 2006; Stone et al., 1992;
535 Torres et al., 2012a). Air-assisted flaring involves the installation of an air blower that
536 provides the forced air at the bottom of the stack. The major advantages of the air-assisted
537 flare are that it is less expensive to run, extends the life-span of the flare by cooling the tip of
538 the flare and is easier to maintain than other configurations (Castineira and Edgar, 2006). For
539 an air-assisted flare, the combustion efficiency of the flare decreases linearly above a
540 threshold limit of the air assist to fuel-gas ratio, but the rate of decrease is slow for a fuel gas
541 with a higher VHV (Torres et al., 2012b). Incomplete combustion can occur when air-fuel
28
542 gas ratios go beyond the optimum value, to the extent that the flame may be put out as a
546 flow of the fuel gas, and used to enhance the mixing of the fuel gas and air in the combustion
547 zone. A high-pressure burner is used to promote atomization of any liquid hydrocarbon and
548 enhance the mixing of the fuel gas with air to bring about a complete or near-complete
549 combustion. Pressure assistance often requires significant amount of space in a remote area
550 because of the burners arrangement at ground level. Fuel gas exit velocity increases with
551 pressure at the burner. Pressure-assisted flares usually have burners arranged on the ground
552 and as such must be located carefully within the oil and gas production plant (Enviroware,
556 the fuel gas and air. The method is often used for gases with low VHV, that is, fuel gas
557 having low C-to-H ratio (in alkanes, C-to-H ratio increases from a minimum value of 0.25 for
558 methane; see Table 2). The C-to-H ratio determines the smoking tendency of hydrocarbons,
559 with smokiness increasing with C-to-H (USEPA, 1995). Note that the experimental fuel gas
560 compositions listed in Table 2 are at or below the lowest C-to-H ratio of the Flow Station gas
561 compositions listed. Non-assisted flaring is used for gases that require smaller amounts of air
563 Table 4 summarises the features of the different types of flares discussed in Section 3.3.
564
565
29
566 Table 4: Summary of the properties of flare types
Steam-assisted Air-assisted Pressure-assisted Non-assisted
Method Steam is Air is introduced The vent pressure No assistance
introduced into from a blower to of the gas flow is is given to the
the combustion enhance the used to enhance combustion
zone to enhance mixing and mixing at the tip process
mixing. turbulence of the of the flare burner
fuel gas in the
combustion zone
Efficiency Most efficient in Less efficient Not as efficient as Only efficient
terms of than the steam steam and air- for non-
suppressing soot assisted flare but assisted but can sooting
formation. Some relatively equally suppress combustion
of the CO formed efficient that the sooting. especially in
can be oxidized other two types light
to CO2 hydrocarbons
Benefits Fuel with high Prolongs the life Enhance Can be used
heat value, and span of the flare combustion for occasional
hence, high tip. Less efficiency when emergency
sooting expensive than the gas flux flaring of near
propensity can be steam-assisted pressure is smokeless gas
disposed of with and easy to sufficiently high
relatively less maintain, hence, enough without
soot it is the most the additional cost
commonly used. of steam and air
generation
Relative size They are often Not as large as May be of same Often smaller
large flares as the steam size as air- in size
they include the assisted. assisted flare compared to
steam generator depending on the the other types
and are usually flow capacity of
employed in the facility
large gas
facilities.
Shortcoming Over-steaming Over- aeration The fluctuation of Cannot be used
can result in can also result in gas flow pressure for dense fuels
reduced less efficiency. A has a bearing with high
efficiency of limit of air assist consequence on sooting
flare. It is also to gas ratio must the efficiency of propensity
expensive to be maintained for the combustion. which are
maintain on a effectiveness of Requires large typical gas in
large-scale the flare. space in a remote oil and gas
area. processing
facilities
567
30
568 4 Estimating emissions from gas flaring
569 Emissions from a typical gas flare can be solids, liquids or gases. As a result of the
570 inaccessible nature of full-scale real-world gas flares, several techniques have been used to
571 quantify the amount of emissions from such flares. Such methods include measurement or
572 source monitoring often by lab-based, pilot-study-based or field-based study (Johnson et al.,
573 2011; McEwen and Johnson, 2012), application of emission factors obtained from
574 measurements and scaling calculations (Giwa et al., 2014; Huang et al., 2015; Sonibare and
575 Akeredolu, 2004; Talebi et al., 2014; USEPA, 1995) and, simulations, often by computational
577 In real-world flares, complete combustion cannot be achieved always and everywhere.
578 Incomplete combustion of the fuel gas can be due to poor efficiency of the flare system,
579 flame temperature (flame dynamical characteristics), insufficient oxygen resulting in poor
580 stoichiometric air/fuel gas mixing ratio, the condition of the fuel gas in the combustion zone
581 and prevailing ambient meteorological condition (Stone et al., 1992). Carbon monoxide (CO)
582 can represent 24 – 80 % (on carbon a molar basis) of emissions for an incomplete combustion
586 adequate estimation of the yield and transport of pollutant species from the combustion
587 process. Flames can be classified along a spectrum ranging from turbulent diffusion flames
588 (of the kinds discussed above) to pool fires (e.g. tar-pool fire) based on the nature and
589 dynamics of the fuel in the flame as well as the design of the burning process (Delichatsios,
590 1987, 1993a). In this review we are concerned with gas flares, which are classified as
591 turbulent jet-diffusion flames. They are so classified because of the high pressure associated
592 with the release of the fuel gas into the flame.
31
593 Jet-diffusion flames in the environment can be categorized based on the momentum flux
594 ratio, R, of the jet plume versus the horizontal momentum flux of the ambient wind (Huang
595 and Wang, 1999). Flares with high R (R>10) may be further categorized depending on
596 whether the flame characteristics are driven by the buoyancy of the hot plume or momentum
597 of the fuel gas (McEwen and Johnson, 2012). Both buoyancy and momentum are important
599 Richardson number, Richardson ratio, fire Froude number, gas Froude number, and Reynolds
600 number - have been used in studies to configure the regime of the flame in a jet-diffusion
601 flame. The Froude numbers measure the ratio of the inertia force on an element of the fluid
602 (in this case, gas or fire) to the weight of (i.e. gravitational force acting on) the fluid element.
603 The fire Froude number, 𝐹𝑟𝑓 , gas Froude number, 𝐹𝑟𝑔 and Reynolds number, Re, have
604 proved to be useful dimensionless parameter to define the flame regime (Becker and Liang,
605 1982; Delichatsios, 1993a; Delichatsios, 1993b; Delichatsios, 1987; Sivathanu and Faeth,
606 1990).
3
𝑢𝑒 𝑓𝑠 ⁄2
607 𝐹𝑟𝑔 = 1⁄ ………………………………………………………………… (4)
1 𝜌 4
(𝑔𝑑𝑒 ) ⁄2 ( 𝑒 )
𝜌∞
3
𝑢𝑒 𝑓𝑠 ⁄2
608 𝐹𝑟𝑓 = 1⁄ ……………………………………………………………... (5)
∆𝑇𝑓 2 𝜌 1⁄4
( 𝑔𝑑𝑒 ) ( 𝑒 )
𝑇∞ 𝜌∞
𝑢𝑒 𝑑𝑒
609 𝑅𝑒𝑠 = …………………………………………………………………………... (6)
ʋ𝑜
1
610 𝑓𝑠 = ………………………………………………………………………………. (7)
𝑆+1
∆𝑇
611 In the definition of 𝐹𝑟𝑓 in equation (5), ( 𝑇 𝑓 ) 𝑔 is the effective acceleration generated by
∞
612 individual hot eddies burning at the flame temperature (Delichatsios, 1987). Among the three
613 dimensionless numbers defined by equations (4) - (6), the Reynolds number is used to
32
614 determine the status of flow, either turbulent or laminar. The fire Froude number, Frf, is used
615 to identify the dominant mechanism between buoyant-generated turbulence and momentum-
616 generated turbulence. In practice, Frf can be used to parameterise soot yield from turbulent
620 atmosphere per unit activity or per unit raw material consumed. It can be obtained from
621 experimental measurements carried out on several sources which represent a particular
622 emission source type. For example, road transport emission factor can be compiled by
623 measuring the amount of each pollutant (CO2, CO, PM, NOx) given off by cars (petrol and
624 diesel), heavy duty vehicles and motorbikes per litre of fuel burned for every km travelled
625 (Gertler et al., 1998; USEPA, 1995; Zhang and Morawska, 2002) under given driving
626 conditions. It is often expressed in g m-3 (pollutant produced per unit volume of raw material
627 consumed), g kg-1 (pollutant produced per unit mass of raw material consumed), g km-1
628 (pollutant produced per unit distance travelled). Emission factors have been compiled by
629 several agencies, which include the US Environmental Protection Agency (USEPA), the
630 European Environmental Agency (EEA), the United Kingdom Department for Environment,
631 Food and Rural Affairs (Defra), and GAINS (Greenhouse gas Air Pollution Interactions and
632 Synergies), for several source categories based on technical sessions of lab-based studies,
633 pilot studies or actual field measurements. Of these agencies, ECLIPSE and USEPA have
634 EFs specifically for emissions from gas flares: the GAINS emission factor for BC from gas
635 flaring is 1.6 g m-3: the equivalent USEPA value has four discrete values between 0 and 6.4 g
636 m-3 depending on the smokiness of the flame (see below). Stohl et al. (2013) in their study
637 used emission factor of 1.6 g m-3 obtained from ECLIPSE (Evaluating the Climate and Air
33
638 Quality Impacts of Short-lived Pollutants) emissions data set to simulation BC emissions in
640 It is difficult to carry out an accurate estimate of emission from gas flares directly from field
641 measurements. Conventional experimental techniques are not suited due to the severe
642 operating conditions that occur in the flaring process and the almost uninhabitable nature of
643 the gas flaring area to both man and the field equipment during the process of gas flaring
644 (Ismail and Umukoro, 2014; McDaniel and Tichenor, 1983; McEwen and Johnson, 2012;
645 RTI, 2011; Talebi et al., 2014). As a result of the unsteady and opaque nature of gas flare
646 flames, remote sensing provides only a partial answer to the difficulties of in-situ monitoring.
647 In Chapter 13 of USEPA’s 5th edition of the compilation of air pollutant emission factors
648 known as AP-42, published in January 1995, emission factors (EF) for pollutant emitted from
649 industrial flaring of waste gas were given; these were recently updated in 2014 (though still
650 in the draft stage) except those for soot. The EFs published in 1995 were based on a study
651 conducted by McDaniel and Tichenor (1983) aimed at determining combustion efficiency
652 and hydrocarbon destruction efficiency for flares operated under different condition. The
653 recent updates of EFs in AP-42 give emission factors of 0.17 kg GJ-1, 1.43 kg GJ-1 and 0 - 6.4
654 g m-3 are given for CO, NOx and soot, respectively.
655 The fuel used in the pilot study to estimate EFs for industrial flare pollutants in the AP-42
656 compilation was predominantly propylene and inert diluents. As such, the EFs, especially for
657 soot, might not be an adequate representation of a typical flare in the oil and gas industry with
658 varying fuel gas composition (Table 2). For AP-42, the collection of soot at the experimental
659 stage was not done by a conventionally accurate method. For instance, the particulate matter
660 (soot) given off was not collected isokinetically in accordance with USEPA’s method 5 for
661 sampling particulate matter from stationary source (USEPA, 2000). An emission factor of 0.0
34
662 μg L-1 (microgram per litre of fuel gas at standard temperature and pressure) given for soot in
663 non-smoking flares of industrial flares in AP-42 is of limited utility and presumably denotes a
664 upper limit of 0.05 μg L-1, because the soot yield from flares is never likely to be precisely
665 zero. Even seemingly modest emission of soot (BC) can be significant, given what is now
666 known of the effects of BC on climate, weather and human health from findings in recent
667 studies (Bond et al., 2013; Ramanathan and Carmichael, 2008; Stocker et al., 2013; Tripathi
669 The mass/mole balancing technique for estimating the yield of pollutants from the
670 combustion of hydrocarbon is a widely used technique. Using this technique, Ismail and
671 Umukoro (2014) and Sonibare and Akeredolu (2004) estimated the yield for SO2, NOX, CO2
672 and CO from the combustion of fuel gas (hydrocarbon with inert diluents) at various levels of
673 combustion. Both studies did not account for unburned carbon (soot) even for their ‘severely’
674 incomplete combustion process reactions. Ismail and Umukoro (2014) varied combustion
675 efficiency (CE) and air available for combustion, and with CE as low as 0.5 (50 %) still did
679 of biomass, solid fuel and fossil fuel (Goldberg, 1985; Koch et al., 2009). Globally, fossil fuel
681 2004). Annually, the contribution of gas flaring to global BC concentration is estimated to be
682 260 Gg (Bond et al., 2013), that is approximately 0.1% of the total contribution from fossil-
683 fuel use, of which Russia is estimated to contribute 81.0 Gg (Huang et al., 2015).
684 The formation and quantification of soot from the combustion of hydrocarbon is a rather
685 complex thermo-chemical process that is not well understood, despite decades of research
35
686 (Castineira and Edgar, 2006; Haynes and Wagner, 1981; Johnson et al., 2011; Maricq, 2009).
687 Soot given off from the combustion of hydrocarbon is predominantly elemental carbon and
688 the amount given off depends on a number of factors, including the efficiency of the
690 The diameter of soot particles emitted ranges between 10 – 200 nm and most commonly lies
691 between 10 and 50 nm. Figure 7 shows a TEM micrograph of soot particle and agglomerates
692 from acetylene flames. The very fine particle of soot are from seemingly 'non-sooting' flames
693 while at the other extreme are those from heavily sooting flames (Flagan and Seinfeld, 2012;
694 Glassman, 1996; Haynes and Wagner, 1981). Combustion of hydrocarbon components of the
695 fuel gas has PAH and BC signatures which can be used as tracers for emissions from the
697
698 Figure 7: TEM micrograph of soot (a) microstructure (b) agglomerates (Tumolva et al.,
699 2010).
700
701 Soot is formed when the carbon particles are cooled below their ignition temperature and
702 there is a deficiency of oxygen (Stone et al., 1992). Considering the variation of the
36
703 composition of gas flared from one station to the other, EFs available for estimating
704 emissions from gas flares are overly generalised. The fuel gas used in most of the studies to
705 estimate the EFs is either propylene or propane or a mixture of both with nitrogen added to
706 alter the heat content as well as the use of predominantly methane-based fuels which are
708 Several attempts have been made to quantify and study the characteristics of soot yield from
709 gas flares using various approaches in lab-based, pilot-based and field-based studies as well
710 as simulation techniques. Emission factors from some of these studies are given in Table 5.
711 Other emission factors not highlighted in this table are based on the USEPA’s AP-42
713
714
715 Table 5: A Summary of Emission Factors/Emission rate for Soot from Industrial Flares
Study (year) Emission factor (g m-3 Emission rate Type of study Fuel
of fuel burned) (g s-1)
80% propylene
USEPA (1995) 0.0, 0.9, 4.2 6.4a - Pilot study 20% propane
Johnson et al. - 2.0 (±0.66) Field study associated
(2011) natural gas
Almanza et al. - 0.025 – 0.22 CFD simulation associated
(2012) natural gas
McEwen and 0.51 - Lab-scale flare hydrocarbon
Johnson (2012) (alkanes)
Johnson et al.
(2013) - 0.067 Field study -
GAINS 2011 1.6 - Modelling -
associated
IMP (2006)b - 3.37 Simulation natural gas
associated
CAPP (2007) 2.563 - - natural gas
a
716 0.0 is for non-smoking flame, 0.9 for light smoking flares, 4.2 for averagely smoking flare
717 and 6.4 for heavily smoking flares
b
718 cited in (Almanza et al., 2012)
719
720 The studies by Johnson et al. (2011) and Johnson et al. (2013) gave emission rates (g s-1) and
721 not emission factors (g m-3) because both were actually field studies that quantified the soot
37
722 given off per unit time by estimating the travel (speed) of the soot in space using a charge
723 couple device camera (CCD) camera viewing a real field flare.
724 Emission of soot from the combustion of hydrocarbon varies with the VHV of the gas and, as
725 such, the emission factor for soot can be estimated from its VHV (McDaniel and Tichenor,
726 1983; RTI, 2011; Torres et al., 2012b). McEwen and Johnson (2012), using a lab-based
727 experiment, studied the BC emission from the combustion of hydrocarbon. The study varied
728 the VHV of the fuel gas and measured the soot yield for each combustion process. The VHV
729 used here is the higher heating value (HHV), or gross heating value. The relationship between
730 the two variables (soot yield and HHV) obtained from the study by McEwen and Johnson
733 When the composition of the fuel gas is known, its VHV can be calculated from standard
734 thermochemical tables (cf. Table 2). It should, however, be noted that because the VHV
735 contains no information on flame dynamics, this relationship is only appropriate for the flame
736 dynamics conditions of the experiment. Although, equation (8) is a readily available
737 relationship to estimate soot emission from hydrocarbon combustion, its application is
738 restricted to a complete or near-complete combustion process. Inserting the estimated VHV
739 into the relation developed by McEwen and Johnson (2012) gives an estimate of the soot
38
747 been considered for such scaling purpose; these include Richardson ratio, 𝑅𝑖𝐿 (Becker and
748 Liang, 1982), fire Froude number, 𝐹𝑟𝑓 (McEwen and Johnson, 2012), and the first
749 Damköhler ratio, Da1 (Becker and Liang, 1982). Richardson ratio, 𝑅𝑖𝐿 , is defined as the ratio
750 of the buoyancy-generated turbulent kinetic energy (TKE) of the flame to the TKE of emitted
751 gas jet at the exit:
𝑔𝐿3 𝜌
752 𝑅𝑖𝐿 = (𝑢 2 ( 𝜌∞) .………………………………………. (9)
𝑒 𝑑𝑒 ) 𝑒
753 Richardson ratio is the basis to assess the turbulent regime of the flame: when 𝑅𝑖𝐿 ≪ 1,
754 buoyancy-induced mixing between emitted gas jet and background air is much weaker than
755 jet-induced mixing, and consequently the flame is dominated by forced convection; when
756 𝑅𝑖𝐿 ≫ 1, jet-induced mixing is much weaker than buoyancy-induced mixing and the flame is
757 dominated by natural convection. Fire Froude number, 𝐹𝑟𝑓 , is defined in Equation (5) and can
758 be interpreted as the ratio of the jet’s inertia to the buoyancy force acting on it. Fire Froude
759 number can be used to assess the dominating force to “stretch” the flame: when 𝐹𝑟𝑓 ≫ 1, the
760 jet’s momentum is the dominating factor and when 𝐹𝑟𝑓 ≪ 1, the flame’s buoyancy force is
761 the dominating factor. Comparing 𝐹𝑟𝑓 against 𝑅𝑖𝐿 , we understand that both are used to assess
762 the dominating factor between jet-related quantity and buoyancy-related quantity, but there
763 are two differences: (1) the quantity is momentum ( velocity) for 𝐹𝑟𝑓 , but TKE ( velocity
764 velocity) for 𝑅𝑖𝐿 ; (2) the ratios between jet-related quantity and buoyancy-related quantity
765 are reciprocal: 𝑅𝑖𝐿 (buoyancy-related quantity)/( jet-related quantity) whereas 𝐹𝑟𝑓 (jet-
766 related quantity)/( buoyancy -related quantity). Therefore, these two parameters are closely
767 related and mathematically their relationship should be 𝑅𝑖𝐿 ∝ 𝐹𝑟𝑓 −2. Precise relationship
768 between them can be derived from Equations (5) and (9):
3⁄ 3
𝜌 2 𝐿 𝑇
769 𝑅𝑖𝐿 = ( 𝜌∞) (𝑑 ) (∆𝑇∞ ) . 𝑓𝑠 3 . 𝐹𝑟𝑓 −2 …..………………………………......... (10)
𝑒 𝑒 𝑓
39
𝐿 3
770 It is noted that ( ) can be interpreted as the volume expansion ratio of the flaring gas due
𝑑𝑒
−1
𝑇∞
771 to burning and (∆𝑇 ) can be approximated as the temperature rising ratio of the flaring gas
𝑓
772 due to burning. Based on the gas law, these two are proportional to each other for an isobaric
773 process (from the exit to the flame tip) which is a good assumption for gas flaring, i.e.
3 −1 3
𝐿 𝑇 𝐿 𝑇
774 (𝑑 ) ∝ (∆𝑇∞ ) , or equivalently, (𝑑 ) (∆𝑇∞ ) ≈ const. Therefore, Equation (10) confirms
𝑒 𝑓 𝑒 𝑓
775 the relationship of 𝑅𝑖𝐿 ∝ 𝐹𝑟𝑓 −2 . This suggests a strong dependence between adopting 𝑅𝑖𝐿 and
776 adopting 𝐹𝑟𝑓 as the scaling parameter.
777 The first Damköhler ratio, Da1, however, is defined as the ratio of residence time of fuel in
778 flame (𝜏𝑟𝑒𝑠 ) to chemical time of burning process (𝜏𝑐ℎ𝑒𝑚 ):
𝐿
( )
= 𝑢𝑒 ⁄𝜏𝑐ℎ𝑒𝑚 .………………………… (11)
𝑟𝑒𝑠𝑖𝑑𝑒𝑛𝑐𝑒 𝑡𝑖𝑚𝑒 𝑖𝑛 𝑓𝑙𝑎𝑚𝑒, 𝜏𝑟𝑒𝑠
779 𝐷𝑎1 = 𝑐ℎ𝑒𝑚𝑖𝑐𝑎𝑙 𝑡𝑖𝑚𝑒, 𝜏𝑐ℎ𝑒𝑚
780 It describes the extent of the oxidation process within the flame in relation to the oxidant’s
781 feed rate. For a large Da1 (i.e. 𝜏𝑟𝑒𝑠 ≫ 𝜏𝑐ℎ𝑒𝑚 ), the velocity fluctuating component does not
782 have much influence on the chemistry of the flame. The chemical reaction is able to proceed
783 to completion within the residence time in the combustion zone, resulting in intensive
784 chemical reaction and hot diffusion flame. For a small Da1 (i.e. 𝜏𝑐ℎ𝑒𝑚 ≫ 𝜏𝑟𝑒𝑠 ), turbulence
785 can significantly affect the chemistry and structure of the flame. The rate of chemical reaction
786 and hence, heat release may be affected, causing combustion product to be mixed with
787 reactants within a time interval shorter than the chemical reaction time (Liberman, 2010;
788 William, 1985). From the perspective of processes, 𝐷𝑎1 involves an extra dimension (i.e.
789 chemical processes) which is not reflected by either 𝑅𝑖𝐿 or 𝐹𝑟𝑓 . In principle, we should
40
791 4.5 Soot modelling
792 Mathematical modelling is a technique that has been used by scientist and engineers over the
793 years to understand the relationship between sets of input and output parameters in a process,
794 especially where the 'real world' process is often remote or grossly complicated to assess. The
795 region of validity of such model outputs is often limited as several assumptions and constants
796 are applied in the modelling to further simplify the process being studied.
797 Soot formation and oxidation in pre-mixed and non-premixed (diffusion) hydrocarbon flames
798 have been studied using several modelling techniques including computational fluid
799 dynamics (CFD). The main problem with mathematical modelling of turbulent combustion of
800 hydrocarbon is the problem of modelling turbulent flow and chemical kinetics and the
801 interaction between flow and chemical reaction (Magnussen and Hjertager, 1977). Kennedy
802 (1997) classified the models for soot formation as empirical correlations, semi-empirical
803 correlation models and models with detailed chemistry and physics of soot formation. The
804 flame temperature, C:H ratio and number of carbon atoms in the fuel (hydrocarbon) are
805 important parameters considered to have strong influence on the sooting propensity of the
806 hydrocarbon (Harris et al., 1986). These parameters have been the basis of measurements
808 In their study to modelling soot formation and combustion, Magnussen and Hjertager (1977)
809 assumed that soot is formed from a gaseous hydrocarbon in two stages (a) formation of
810 radical nuclei, and (b) soot particle formation from these nuclei. Applying expressions for the
811 rate of formation of radical nuclei and rate of soot particle formation expressions developed
812 by Tesner et al. (1971a); (Tesner et al., 1971b) in their model, tested on both pre-mixed and
813 diffusion flame, they predicted soot concentrations that are in close agreement with
814 experimental data. They concluded that soot formed and was contained in eddies and burned
41
816 Moss et al. (1989), using a two-equation model for the evolution of soot volume fraction and
817 number density, simulated the formation of soot. They included the influences of nucleation,
818 surface growth and coagulation on the rate of soot formation. As given in equations 12(a) and
819 (b), their model contains simplified expressions to quantify particle nucleation, growth and
820 coagulation; using three empirical constants that are dependent on the fuel to control the rate
821 of these processes. A major finding from their study is that soot volume fraction is
𝑛
𝑑( ) 𝑛
𝑁𝑜
823 { } = 𝛼(𝜁) − 𝛽(𝜁)(𝑁 )2 ........................................................................... 12(a)
𝑑𝑡 0
𝑑𝑓
824 𝜌𝑠 { 𝑑𝑡𝑣} = 𝛾(𝜁)𝑛 + 𝛿(𝜁) .............................................................................. 12(b)
825 where ρs is the assumed density for solid carbon (1.8 x 103 kgm-3), No is Avogadro's number
826 and 𝜁 is the mixture fraction. Rates of the processes are expected to be a function of the
827 mixture fraction. The rates of the processes are defined explicitly in terms of the fuel density,
−𝑇𝛼⁄
829 𝛼 ≡ 𝐶𝛼 𝜌2 𝑇 1/2 𝑋𝑐 exp ( 𝑇) ; 𝛽 ≡ 𝐶𝛽 𝑇
1/2
−𝑇
830 𝛾 ≡ 𝐶𝛾 𝜌𝑇 1/2 𝑋𝑐 exp ( 𝛾⁄𝑇) ; 𝛿 ≡ 𝐶𝛿 𝛼 ...................................... 13
831 coefficients Cα,β,γ,δ and activation temperatures, Tα and Tγ are obtained from experimental
832 data. In equations 12(a), the first and second term on the left are nucleation and coagulation
833 processes, respectively. And, in equation 12(b), the first and second term on the left represent
42
835 Lautenberger et al. (2005) developed a CFD model to study the formation and oxidation of
836 soot. Since their model was to generate enough accurate predictions of soot emission
837 concentrations in order to estimate CFD simulations of fire radiation from turbulent flames,
838 considerations were given only to phenomena which were essential. Soot estimation in their
839 model was based on a further simplified form of the two-equation model of Moss et al.
840 (1989). There are no fuel-specific constants in their estimations, but rather, they used the
841 laminar smoke point height to account for the sooting propensity of different fuels. The
842 smoke point of a flame is its length just before the onset of the release of visible smoke.
843 Length of a flame is dependent on the extent of completeness of the combustion and heat
845 The use of CFD simulation soot has made it possible to make predictions about the
846 relationships between the various processes involved in soot formation and oxidation in both
847 pre-mixed and diffusion hydrocarbon flames as well as the sensitivity of the soot formation to
848 some of the complex phenomena. However, the veracity of CFD simulation results in limited
852 flaring emissions are EDGAR (Emissions Database for Global Atmospheric Research) and
853 ECLIPSE (Evaluating the Climate and Air Quality ImPacts of Short-livEd Pollutants)
855 In EDGAR v4.2, for gas venting and flaring, emissions were calculated for 1994 onwards
856 using the amount of gas flared estimated from satellite observations of intensities of light
857 from various gas flares. The estimated quantities of gas flared and emission factors obtained
858 from either inventory guidance documents or confidential information were used to generate
43
859 gridded annual emission datasets for different countries on a resolution of 0.1o x 0.1o. In
860 2005, the inclusion of new primary data sources for gas flaring in EDGAR v4.2 gives rise to
861 a change of +75 % of EDGAR 4.1 value (130 Tg) in global CO2 emission (European
863 ECLIPSE is provided by the IIASA (International Institute for Applied Systems Analysis).
864 Emission calculations for the historical years (2005 – 2010) were developed in a series of
865 regional and global projects. For gas flaring, emissions were calculated using data available
866 from NOAA, NASA and the World Bank collaborative work to estimate the volume of gas
867 flared globally. The volume of gas flared was estimated using NASA MODIS active fire
868 detection products (Elvidge et al., 2007; Elvidge et al., 2011). Emission factors and other
870 various technologies. The calculation was performed with the IIASA GAINS model (Klimont
871 et al., 2013). The ECLIPSE v4 global emission dataset is available on a 0.5o x 0.5o lon-lat
872 resolution.
873 Stohl et al. (2013), performing a 3-year black carbon (BC) simulation in the Arctic, using a
874 Lagrangian particle dispersion model FLEXPART driven by the ECLIPSE dataset, estimated
875 that more than 40 % of mean surface BC surface concentration in the Arctic is attributable to
876 gas flaring. Although, more studies on the sensitivity of global models to gas flaring emission
877 is needed to give a clear and precise quantification of the contribution of gas flaring emission
878 to global aerosol loadings, the change in global CO2 level in EDGAR4.2 underscores the
880 5 Conclusion
881 The World Bank has been at the fore-front of the campaign to reduce gas flaring through the
44
883 (http://www.worldbank.org/en/programs/gasflaringreduction). The “Zero Routine Flaring by
884 2030” initiative was launched in April 2015 by the World Bank, United Nations,
885 governments and oil companies. As at April 17 2015, a total of nine countries have agreed to
886 the “Zero Routine Flaring by 2030” initiative. At the time of writing, some major flaring
887 nations have yet to sign up to the initiative. Indeed, some major flaring countries still struggle
888 to meet targets for gas flaring set in the late 1980s. A steep increase in flaring since 2006 has
889 been reported for the USA and Canada, associated, presumably, with the exploitation of
891 Considering the wide spectrum, quantity and effects of pollutants emitted from gas flaring,
892 coupled with the estimated quantity of gas flared globally, it is surprising that so little effort
893 has been put into adequately understanding the yield of pollutants, especially BC, from the
894 process in real world field situations. The wide variation of fuel gas compositions from flow
895 stations around the world underpins the importance of developing strategies that take these
897 The steep decrease in the fraction of total gas production flared between 2000 and 2006
898 seems to have stabilised between 2007 and 2010, and the fraction has even increased in 2011.
899 The overall quantity of gas flared since 2000 has been steady between 93 and 110 bcm. An
900 increase in total production since 2009 has brought about a corresponding increase in the
901 quantity of gas flared. Incentives and stringent policies are not yet in place to encourage more
902 companies and countries to partner with the World Bank in their “Zero Routine Flaring by
904 Elevated concentrations of BC, CO, H2S, SO2, NO, NO2 and PAH measured around flaring
905 sites (ground based and aircraft) is indicative of the detrimental impact of gas flaring on the
906 environment. Clusters of gas flaring sites around the tropics and near-tropic regions of the
45
907 world, where there is the likelihood on enhanced atmospheric mixing of the emissions into
908 the lower and even mid-troposphere, coupled with the high temperature of the emitted plume
910 Emission factors used for BC emission from gas flaring are inadequate to estimate emission
911 from a typical real-world gas flare as most of the fuels used in the studies for such emission
912 factors are not representative of fuel gas from most Flow Stations around the world. In the
913 studies that employed the mass/mole balancing technique to estimate pollutants from the
914 combustion of hydrocarbons, estimated EFs for CO2 and CO have been given, but there is no
915 consideration of the amount of unburnt carbon given off. It should be noted that these studies
916 did not consider flame dynamics changes in their estimations. When estimating emissions
917 from gas flares, there is the need to ascertain the nature (regime) of the combustion flame as
918 the flame nature and temperature plays pivotal role in determining the pollutants yield, and
919 this has not yet been routinely considered in dispersion modelling and global inventories.
920 Global models need to update the sources of BC to include gas flaring, especially in regions
921 prone to long-range transfer of gas flaring emission from leading gas flaring nations
923
924 Acknowledgment
925 Olusegun G. Fawole is highly grateful to the UK government for funding his PhD studies
927 appreciates the anonymous reviewers for their helpful criticisms and comments.
928
929
46
930
931
932 6 Glossary
933 1. Heating value – the heat evolved per unit of a gas when it undergoes complete
934 combustion at standard temperature and pressure. It can be measured in unit mass,
936 2. Higher (gross) heating value (HHV) – The heating value when all of the water present
938 3. Lower (net) heating value (LHV) - The heating value when all of the water present in
940 4. Volumetric heating value (VHV) – Total heat evolved per unit volume of a gas.
941 5. Equivalent ratio, ϕ - This is the ratio of the stoichiometric air-fuel mass ratio to the
944 7. Laminar smoke point - The length of flame from a gas jet just before the onset of the
946
947
948
949
47
950 7 References
951 Abdulkareem, A. and Odigure, J., 2006. Deterministic Model for Noise Dispersion from gas
952 Flaring: A case study of Niger–Delta area of Nigeria. Chemical and biochemical
954 Almanza, V., Molina, L. and Sosa, G., 2012. Soot and SO 2 contribution to the supersites in
955 the MILAGRO campaign from elevated flares in the Tula Refinery. Atmospheric
957 Ana, G., Sridhar, M. and Emerole, G., 2012. Polycyclic aromatic hydrocarbon burden in
958 ambient air in selected Niger Delta communities in Nigeria. Journal of the Air &
960 Arctic Council, 2013. Recommendations to reduce black carbon and methane emissions to
961 slow Arctic climate change, Arctic task force on short-lived climate forcers.
962 Avwiri, G. and Nte, F., 2004. Environmental sound quality of some selected flow stations in
963 the Niger Delta of Nigeria. Journal of Applied Sciences and Environmental
965 Barry, R.G. and Chorley, R.J., 2009. Atmosphere, weather and climate. Routledge.
966 Becker, H. and Liang, D., 1982. Total emission of soot and thermal radiation by free
968 Beychok, M.R., 1994. Fundamentals of stack gas dispersion, 63. Milton R. Beychok Irvine.
969 Bond, T.C., Doherty, S.J., Fahey, D., Forster, P., Berntsen, T., DeAngelo, B., Flanner, M.,
970 Ghan, S., Kärcher, B. and Koch, D., 2013. Bounding the role of black carbon in the
48
973 Bond, T.C., Streets, D.G., Yarber, K.F., Nelson, S.M., Woo, J.H. and Klimont, Z., 2004. A
974 technology‐based global inventory of black and organic carbon emissions from
976 BP, 2013. Statistical review of world energy, British Petroleum, Great Britain
977 www.bp.com/statisticalreview.
979 Burney, J. and Ramanathan, V., 2014. Recent climate and air pollution impacts on Indian
981 CAPP, 2007. A recommended approach to completing the national pollutant release
982 inventory (NPRI) for the upstream oil and gas industry, Canadian Association of
984 Casadio, S., Arino, O. and Serpe, D., 2012. Gas flaring monitoring from space using the
986 Castineira, D. and Edgar, T.F., 2006. CFD for simulation of steam-assisted and air-assisted
988 Chung, S.H. and Seinfeld, J.H., 2005. Climate response of direct radiative forcing of
991 Davoudi, M., Rahimpour, M., Jokar, S., Nikbakht, F. and Abbasfard, H., 2013. The major
992 sources of gas flaring and air contamination in the natural gas processing plants: a
993 case study. Journal of Natural Gas Science and Engineering, 13: 7-19.
994 Delichatsios, M., 1993a. Smoke yields from turbulent buoyant jet flames. Fire safety journal,
49
996 Delichatsios, M., 1993b. Transition from momentum to buoyancy-controlled turbulent jet
997 diffusion flames and flame height relationships. Combustion and flame, 92(4): 349-
998 364.
999 Delichatsios, M.A., 1987. Air entrainment into buoyant jet flames and pool fires. Combustion
1001 Dung, E.J., Bombom, L.S. and Agusomu, T.D., 2008. The effects of gas flaring on crops in
1003 E & P Forum, 1994. Methods of estimating atmospheric emission from E&P operations;
1004 OGP Report No. 2.59/197 The Oil Industry International Exploration and Production
1006
1007 Edwards, P., Young, C., Aikin, K., deGouw, J., Dubé, W., Geiger, F., Gilman, J., Helmig, D.,
1008 Holloway, J. and Kercher, J., 2013. Ozone photochemistry in an oil and natural gas
1009 extraction region during winter: simulations of a snow-free season in the Uintah
1011 Edwards, P.M., Brown, S.S., Roberts, J.M., Ahmadov, R., Banta, R.M., deGouw, J.A., Dube,
1012 W.P., Field, R.A., Flynn, J.H., Gilman, J.B., Graus, M., Helmig, D., Koss, A.,
1013 Langford, A.O., Lefer, B.L., Lerner, B.M., Li, R., Li, S.-M., McKeen, S.A., Murphy,
1014 S.M., Parrish, D.D., Senff, C.J., Soltis, J., Stutz, J., Sweeney, C., Thompson, C.R.,
1015 Trainer, M.K., Tsai, C., Veres, P.R., Washenfelder, R.A., Warneke, C., Wild, R.J.,
1016 Young, C.J., Yuan, B. and Zamora, R., 2014. High winter ozone pollution from
1017 carbonyl photolysis in an oil and gas basin. Nature, advance online publication.
1018 Elvidge, C.D., Baugh, K., Pack, D., Milesi, C. and Erwin, E., 2007. A twelve year record of
1019 national and global gas flaring volumes estimated using satellite data. Boulder, CO,
50
1021 Elvidge, C.D., Baugh, K.E., Ziskin, D., Anderson, S. and Ghosh, T., 2011. Estimation of gas
1022 flaring volumes using NASA MODIS fire detection products. NOAA National
1024 Elvidge, C.D., Zhizhin, M., Baugh, K., Hsu, F.-C. and Ghosh, T., 2015. Methods for Global
1025 Survey of Natural Gas Flaring from Visible Infrared Imaging Radiometer Suite Data.
1027 Elvidge, C.D., Ziskin, D., Baugh, K.E., Tuttle, B.T., Ghosh, T., Pack, D.W., Erwin, E.H. and
1028 Zhizhin, M., 2009. A fifteen year record of global natural gas flaring derived from
1031 European Commission, 2009. Emission Database for Global Atmospheric Research
1034 Feichter, J. and Stier, P., 2012. Assessment of black carbon radiative effects in climate
1036 Flagan, R.C. and Seinfeld, J.H., 2012. Fundamentals of air pollution engineering. Courier
1038 Flanner, M.G., Zender, C.S., Randerson, J.T. and Rasch, P.J., 2007. Present‐day climate
1039 forcing and response from black carbon in snow. Journal of Geophysical Research:
1041 Fortner, E., Brooks, W., Onasch, T., Canagaratna, M., Massoli, P., Jayne, J., Franklin, J.,
1042 Knighton, W., Wormhoudt, J. and Worsnop, D., 2012. Particulate Emissions
1043 Measured During the TCEQ Comprehensive Flare Emission Study. Industrial &
51
1045 Gertler, A.W., Sagebiel, J.C., Dippel, W.A. and Farina, R.J., 1998. Measurements of dioxin
1046 and furan emission factors from heavy-duty diesel vehicles. Journal of the Air &
1048 GGFR, 2012. Estimated Flared Volumes from Satellite Data, 2007-2011., Global Gas Flaring
1050 2015).
1051 GGFR, 2013. The news flare (issue 14), Washington DC ; World Bank.
1052 http://documents.worldbank.org/curated/en/2013/06/18082251/news-flare-issue-14.
1053 Giwa, S.O., Adama, O.O. and Akinyemi, O.O., 2014. Baseline black carbon emissions for
1054 gas flaring in the Niger Delta region of Nigeria. Journal of Natural Gas Science and
1057 Goldberg, E.D., 1985. Black carbon in the environment: properties and distribution.
1059 Harris, M.M., King, G.B. and Laurendeau, N.M., 1986. Influence of temperature and
1062 Haynes, B. and Wagner, H.G., 1981. Soot formation. Progress in Energy and Combustion
1064 Hodgson, S., Nieuwenhuijsen, M.J., Colvile, R. and Jarup, L., 2007. Assessment of exposure
1065 to mercury from industrial emissions: comparing “distance as a proxy” and dispersion
1067 Holden, J., 2005. An introduction to physical geography and the environment. Pearson
1068 Education.
52
1069 Huang, K., Fu, J.S., Prikhodko, V.Y., Storey, J.M., Romanov, A., Hodson, E.L., Cresko, J.,
1070 Morozova, I., Ignatieva, Y. and Cabaniss, J., 2015. Russian anthropogenic black
1071 carbon: Emission reconstruction and Arctic black carbon simulation. Journal of
1073 Huang, R.F. and Wang, S.M., 1999. Characteristic flow modes of wake-stabilized jet flames
1075 IEA, 2012. IEA statistics: Natural gas information, International Energy Agency.
1076 IEA, 2013. Key World Energy Statistics International Energy Agency, France.
1077 IMP, 2006. Estudio de las emisiones de la zona industrial de Tula y su impacto en la calidad
1079 IPCC, 2007. Climate Change 2007: The Physical Science Basis. Contribution of Working
1080 Group I to the Fourth Assessment Report of the Intergovernmental Panel on Climate
1081 Change [Solomon, S., D. Qin, M. Manning, Z. Chen, M. Marquis, K.B. Averyt,
1082 M.Tignor and H.L. Miller (eds.)]. Cambridge University Press, Cambridge, United
1084 Ismail, O.S. and Umukoro, G.E., 2014. Modelling combustion reactions for gas flaring and
1086 http:dx.doi.org/10.1016/j.jksues.2014.02.003.
1087 Jacobson, M.Z., 2002. Control of fossil‐fuel particulate black carbon and organic matter,
1088 possibly the most effective method of slowing global warming. Journal of
1090 Johansson, J.K., Mellqvist, J., Samuelsson, J., Offerle, B., Lefer, B., Rappenglück, B., Flynn,
1091 J. and Yarwood, G., 2014. Emission measurements of alkenes, alkanes, SO2, and
1092 NO2 from stationary sources in Southeast Texas over a 5 year period using SOF and
53
1094 Johnson, M., Devillers, R. and Thomson, K., 2013. A Generalized Sky-LOSA Method to
1095 Quantify Soot/Black Carbon Emission Rates in Atmospheric Plumes of Gas Flares.
1097 Johnson, M.R. and Coderre, A.R., 2011. An analysis of flaring and venting activity in the
1098 Alberta upstream oil and gas industry. Journal of the Air & Waste management
1100 Johnson, M.R., Devillers, R.W. and Thomson, K.A., 2011. Quantitative field measurement of
1101 soot emission from a large gas flare using Sky-LOSA. Environmental science &
1103 Kennedy, I.M., 1997. Models of soot formation and oxidation. Progress in Energy and
1105 Klimont, Z., Smith, S.J. and Cofala, J., 2013. The last decade of global anthropogenic sulfur
1107 Koch, D., Schulz, M., Kinne, S., McNaughton, C., Spackman, J.R., Balkanski, Y., Bauer, S.,
1108 Berntsen, T., Bond, T.C. and Boucher, O., 2009. Evaluation of black carbon
1109 estimations in global aerosol models. Atmospheric Chemistry and Physics, 9(22):
1110 9001-9026.
1111 Kostiuk, L., Johnson, M. and Thomas, G., 2004. University Of Alberta Flare Research
1113 http://www.mece.ualberta.ca/groups/combustion/flare/papers/Final%20Report2004.pd
1114 f.
1115 Lautenberger, C.W., de Ris, J.L., Dembsey, N.A., Barnett, J.R. and Baum, H.R., 2005. A
1116 simplified model for soot formation and oxidation in CFD simulation of non-
54
1118 Law, K., Fierli, F., Cairo, F., Schlager, H., Borrmann, S., Streibel, M., Real, E., Kunkel, D.,
1119 Schiller, C. and Ravegnani, F., 2010. “Air mass origins influencing TTL chemical
1120 composition over West Africa during 2006 summer monsoon” published in Atmos.
1121 Chem. Phys., 10, 10753–10770, 2010. Atmos. Chem. Phys, 10: 10939-10940.
1122 Leahey, D. and Davies, M., 1984. Observations of plume rise from sour gas flares.
1124 Leahey, D.M., Preston, K. and Strosher, M., 2001. Theoretical and observational assessments
1125 of flare efficiencies. Journal of the Air & Waste Management Association, 51(12):
1126 1610-1616.
1127 Liberman, M.A., 2010. Introduction to physics and chemistry of combustion: explosion,
1129 MacKay, D., 2008. Sustainable Energy-without the hot air. UIT Cambridge.
1130 Magnussen, B.F. and Hjertager, B.H., 1977. On mathematical modeling of turbulent
1131 combustion with special emphasis on soot formation and combustion, Symposium
1133 Mari, C.H., Reeves, C.E., Law, K.S., Ancellet, G., Andrés‐Hernández, M.D., Barret, B.,
1134 Bechara, J., Borbon, A., Bouarar, I. and Cairo, F., 2011. Atmospheric composition of
1135 West Africa: highlights from the AMMA international program. Atmospheric Science
1137 Maricq, M.M., 2009. An examination of soot composition in premixed hydrocarbon flames
1138 via laser ablation particle mass spectrometry. Journal of Aerosol Science, 40(10):
1139 844-857.
1140 Mathon, V. and Laurent, H., 2001. Life cycle of Sahelian mesoscale convective cloud
1141 systems. Quarterly Journal of the Royal Meteorological Society, 127(572): 377-406.
55
1142 McAllister, S., Chen, J.-Y., Fernandez-Pello, A.C., Fernandez-Pello, A.C. and Fernandez-
1144 McDaniel, M. and Tichenor, B.A., 1983. Flare efficiency study, US Environmental
1146 http://www.tceq.state.tx.us/assets/public/implementation/air/rules/Flare/Resource_1.p
1147 df.
1148 McEwen, J.D. and Johnson, M.R., 2012. Black carbon particulate matter emission factors for
1149 buoyancy-driven associated gas flares. Journal of the Air & Waste Management
1151 Miyauchi, T., Mori, Y. and Yamaguchi, T., 1981. Effect of steam addition on NO formation,
1153 Moss, J., Stewart, C. and Syed, K., 1989. Flowfield modelling of soot formation at elevated
1155 Müller-Dethlefs, K. and Schlader, A., 1976. The effect of steam on flame temperature,
1156 burning velocity and carbon formation in hydrocarbon flames. Combustion and flame,
1158 Nwaichi, E.O. and Uzazobona, M.A., 2011. Estimation of the CO2 Level due to Gas Flaring
1159 in the Niger Delta. Research Journal of Environmental Sciences, 5(6): 565-572.
1160 Obanijesu, E., Adebiyi, F., Sonibare, J. and Okelana, O., 2009. Air-borne SO2 Pollution
1163 OGP, 2000. Flaring and venting in the oil and gas exploration and production industry: An
1164 overview of purpose, quantities, issues, practices and trends. Report No. 2.79/288,
1166 http://www.ogp.org.uk/pubs/288.pdf.
56
1167 Osuji, L.C. and Adesiyan, S.O., 2005. The Isiokpo Oil‐Pipeline Leakage: Total Organic
1168 Carbon/Organic Matter Contents of Affected Soils. Chemistry & biodiversity, 2(8):
1169 1079-1085.
1170 Osuji, L.C. and Onojake, C.M., 2004. Trace Heavy Metals Associated with Crude Oil: A
1171 Case Study of Ebocha‐8 Oil‐Spill‐Polluted Site in Niger Delta, Nigeria. Chemistry &
1173 Ouf, F.-X., Vendel, J., Coppalle, A., Weill, M. and Yon, J., 2008. Characterization of soot
1174 particles in the plumes of over-ventilated diffusion flames. Combustion Science and
1176 Pope III, C.A., Burnett, R.T., Thun, M.J., Calle, E.E., Krewski, D., Ito, K. and Thurston,
1177 G.D., 2002. Lung cancer, cardiopulmonary mortality, and long-term exposure to fine
1178 particulate air pollution. Journal of American Medical Association, 287(9): 1132-
1179 1141.
1180 Quinn, P., Bates, T., Baum, E., Doubleday, N., Fiore, A., Flanner, M., Fridlind, A., Garrett,
1181 T., Koch, D. and Menon, S., 2008. Short-lived pollutants in the Arctic: their climate
1182 impact and possible mitigation strategies. Atmospheric Chemistry and Physics, 8(6):
1183 1723-1735.
1184 Ramana, M., Ramanathan, V., Feng, Y., Yoon, S., Kim, S., Carmichael, G. and Schauer, J.,
1185 2010. Warming influenced by the ratio of black carbon to sulphate and the black-
1187 Ramanathan, V. and Carmichael, G., 2008. Global and regional climate changes due to black
1189 Reeves, C., Formenti, P., Afif, C., Ancellet, G., Attié, J.-L., Bechara, J., Borbon, A., Cairo,
1190 F., Coe, H. and Crumeyrolle, S., 2010. Chemical and aerosol characterisation of the
57
1191 troposphere over West Africa during the monsoon period as part of AMMA.
1193 Roiger, A., Thomas, J.-L., Schlager, H., Law, K.S., Kim, J., Schäfler, A., Weinzierl, B.,
1194 Dahlkötter, F., Krisch, I. and Marelle, L., 2015. Quantifying emerging local
1195 anthropogenic emissions in the Arctic region: the ACCESS aircraft campaign
1197 RTI, 2011. Emission protocol for petroleum refineries v 2.1.1, RTI International, Research
1199 http://www.epa.gov/ttnchie1/efpac/protocol/Emission_Estimation_Protocol_for_Petro
1200 leum_Refinerie_052011.pdf.
1201 Schultz, C., 2014. Texas refinery air pollution emissions are being severely underestimated.
1203 Seinfeld, J., 2008. Atmospheric science: Black carbon and brown clouds. Nature geoscience,
1205 Sivathanu, Y. and Faeth, G.M., 1990. Soot volume fractions in the overfire region of
1207 Smith, D. and Chughtai, A., 1995. The surface structure and reactivity of black carbon.
1208 Colloids and Surfaces A: Physicochemical and Engineering Aspects, 105(1): 47-77.
1209 Sonibare, J., Adebiyi, F., Obanijesu, E. and Okelana, O., 2010. Air quality index pattern
1212 Sonibare, J.A. and Akeredolu, F.A., 2004. A theoretical prediction of non-methane gaseous
1213 emissions from natural gas combustion. Energy Policy, 32(14): 1653-1665.
1214 Stocker, T., Qin, D., Plattner, G., Tignor, M., Allen, S., Boschung, J., Nauels, A., Xia, Y.,
1215 Bex, B. and Midgley, B., 2013. IPCC, 2013: climate change 2013: the physical
58
1216 science basis. Contribution of working group I to the fifth assessment report of the
1218 Stohl, A., Klimont, Z., Eckhardt, S., Kupiainen, K., Shevchenko, V., Kopeikin, V. and
1219 Novigatsky, A., 2013. Black carbon in the Arctic: the underestimated role of gas
1220 flaring and residential combustion emissions. Atmospheric Chemistry and Physics,
1222 Stone, D.K., Lynch, S.K., Pandullo, R.F., Evans, L.B. and Vatavuk, W.M., 1992. Flares. Part
1225 Strosher, M.T., 2000. Characterization of emissions from diffusion flare systems. Journal of
1227 Sultan, B. and Janicot, S., 2003. The West African monsoon dynamics. Part II: The
1228 “preonset” and “onset” of the summer monsoon. Journal of climate, 16(21): 3407-
1229 3427.
1230 Talebi, A., Fatehifar, E., Alizadeh, R. and Kahforoushan, D., 2014. The Estimation and
1231 Evaluation of New CO, CO2, and NOx Emission Factors for Gas Flares Using Pilot
1232 Scale Flare. Energy Sources, Part A: Recovery, Utilization, and Environmental
1234 Tesner, P., Smegiriova, T. and Knorre, V., 1971a. Kinetics of dispersed carbon formation.
1236 Tesner, P., Tsygankova, E., Guilazetdinov, L., Zuyev, V. and Loshakova, G., 1971b. The
59
1239 Torres, V.M., Herndon, S. and Allen, D.T., 2012a. Industrial flare performance at low flow
1242 Torres, V.M., Herndon, S., Kodesh, Z. and Allen, D.T., 2012b. Industrial flare performance
1243 at low flow conditions. 1. Study overview. Industrial & Engineering Chemistry
1245 Tripathi, S., Dey, S., Tare, V. and Satheesh, S., 2005. Aerosol black carbon radiative forcing
1247 Tumolva, L., Park, J.-Y., Kim, J.-s., Miller, A.L., Chow, J.C., Watson, J.G. and Park, K.,
1248 2010. Morphological and elemental classification of freshly emitted soot particles and
1249 atmospheric ultrafine particles using the TEM/EDS. Aerosol Science and Technology,
1251 USEPA, 1995. AP 42-Compilation of air pollutant emission factors. Section 13.5: Industrial
1252 flares, U.S Environmental Protection Agenecy, Office of Air quality planning and
1254 http://www.epa.gov/ttn/chief/ap42/ch13/final/c13s05.pdf.
1255 USEPA, 2000. Test Method 5-Determination of particulate matter emissions from stationary
1257 www.epa.gov/ttnemc01/methods/methods5.html.
1258 USEPA, 2008. Direct emissions from stationary combustion sources EPA430-K-08-003,
1260 http://www.epa.gov/climateleadership/documents/resources/stationarycombustionguid
1261 ance.pdf.
1262 USEPA, 2010. Integrated Science Assessment for Particulate Matter, US Environmental
60
1264 USEPA, 2011. Emission Estimation Protocol for Petroleum Refineries, US Environmental
1266 USEPA, 2012. Report to congress on Black Carbon. EPA-450/R-12-001, United States
1268 USEPA, 2014. AP 42-Compilation of air pollutant emission factors. Section 13.5: Industrial
1269 flares
1270 U.S Environmental Protection Agenecy, Office of Air quality planning and standards,
1272 http://www.epa.gov/ttn/chief/ap42/ch13/final/dc13s05_8-19-14.pdf.
1273 Villasenor, R., Magdaleno, M., Quintanar, A., Gallardo, J., López, M., Jurado, R., Miranda,
1274 A., Aguilar, M., Melgarejo, L. and Palmerın, E., 2003. An air quality emission
1275 inventory of offshore operations for the exploration and production of petroleum by
1277 Wang, X., Heald, C., Ridley, D., Schwarz, J., Spackman, J., Perring, A., Coe, H., Liu, D. and
1278 Clarke, A., 2014. Exploiting simultaneous observational constraints on mass and
1279 absorption to estimate the global direct radiative forcing of black carbon and brown
1281 William, F., 1985. Combustion Theory: The fundamental theory of chemically reacting flow
1283 Zhang, J.J. and Morawska, L., 2002. Combustion sources of particles: 2. Emission factors
1285
1286
61