Addressing Uncertainty
Addressing Uncertainty
Acknowledgements
Table of Contents
INTRODUCTION 7
DOCUMENT AT A GLANCE 8
1. INTRODUCTION 10
1.1 Importance of accurate and reliable GHG accounting 10
1.2 Overview of uncertainty terminology 11
1.3 Types of errors 11
1.4 Numerical determination of uncertainty intervals 12
1.5 Emissions inventory uncertainty assessment 12
2. ASSESSMENT OF UNCERTAINTY RELEVANT
TO GREENHOUSE GAS EMISSION QUANTIFICATION 13
2.1 Overview of emissions inventory uncertainty 13
2.2 Sources of measurement uncertainty 15
2.3 Emission calculation approaches 15
2.4 Inventory steps and data aggregation 18
2.5 Emissions inventory and uncertainty
assessment in the petroleum and natural gas industry 18
3. OVERVIEW OF MEASUREMENT PRACTICES FOR
ESTIMATING GREENHOUSE GAS EMISSIONS 22
3.1 Flow measurement practices 23
3.1.1 Measurements by Orifice Meters 23
3.1.2. Measurement of flow to flares 25
3.1.3 Flow measurements uncertainty analysis 26
3.1.4 Uncertainty specifications for “custody transfer” measurements 26
3.2 Flow measurements for GHG inventories 28
3.2.1 Flow meter types 28
3.2.2 Fuel gas measurement 28
3.2.3. Integrated measurement systems 30
3.3 Sampling and analysis for quantifying GHG emissions 32
3.3.1 Gaseous samples collection and handling 32
3.3.2 Quantifying sampling precision 33
3.4 Carbon content measurement practices 34
3.4.1 Laboratory-based measurements 34
3.4.2 On-line measurements 34
3.4.3 Carbon content calculation 36
3.5 Heat content determination 37
3.5.1 Direct measurements 38
3.5.2 Computational methods 38
3.6 Venting and fugitive emissions measurements 40
3.6.1 Vented emission measurements 41
3.6.2 Fugitive emission measurements 41
3.7 Laboratory management system 42
2
ADDRESSING UNCERTAINTY IN OIL & NATURAL GAS INDUSTRY GREENHOUSE GAS INVENTORIES
3
IPIECA • API
List of Tables
Table 2-1. Overview of methods used to estimate emissions uncertainty 14
Table 2-2. Summary of typical sources of measurement uncertainty 16
Table 2-3. Common emission estimation approaches and their error sources 17
Table 3-1. Example of flare flow measurement systems combined uncertainty 26
Table 3-2. Summary of Alberta ERCB accuracy requirements 27
Table 3-3. Compilation of specifications for common flow meters 29 - 30
Table 3-4. Summary of selected carbon content measurement methods 35
Table 3-5. Summary of selected heating value measurement methods 39
Table 3-6. GHGRP measurement methods for vented and fugitive emission sources 40
Table 4-1. Inter-Country comparison of select mandatory GHG reporting programs 45 - 47
Table 4-2. Comparison of reporting scope between EPA’s GHGRP and CDP 47
Table 4-3. EU regulation maximum permissible uncertainty for activity data tiers 48
Table 4-4. Highlights of calibration requirements under the EPA GHGRP 49
Table 4-5. Estimated uncertainty levels based on Australia National
Greenhouse Gas and Energy Reporting (NGER) 51
Table 5-1. Example of Pedigree Matrix for determining uncertainty
scaling factors based on data quality ratings 55
Table 5.2. Basic Pedigree Matrix uncertainty factors for selected pollutants 56
Table 5-3. Adapted IPCC Template: Approach for Aggregating Uncertainty 59
Table 5-4. Example uncertainty calculation and reporting table 60
Table 5-5. Pedigree Matrix example 60
Table 5-6. Example of flow measurement total uncertainty 61
Table 5-7. Emission uncertainty ranking for onshore oil production example 63
List of Figures
Section
2-1 Hierarchy of emission estimation approaches 19
5-1 Measurement error over time of an unbiased estimate 54
5-2 Measurement error over time of a biased estimate 54
5-3 Onshore oil field: summary of CO2 equivalent emissions and uncertainties 63
List of Figures in Appendices document
F-1 Decision diagram for emission factor uncertainty F-5
F-2 Decision diagram for measurement uncertainty F-6
F-3 Step C – decision diagram for uncertainty aggregation F-22
4
ADDRESSING UNCERTAINTY IN OIL & NATURAL GAS INDUSTRY GREENHOUSE GAS INVENTORIES
Acronym List
ADEME French Environment and IPCC Intergovernmental Panel on
Energy Management Agency Climate Change
AF Activity Factor IPIECA Global oil and gas industry
AGA American Gas Association association for environmental
ANSI American National and social issues
Standards Institute ISO International Standards
API American Petroleum Institute Organization
ASTM American Society of kW Kilowatt
Testing and Materials kWh Kilowatt-hour
AWP Alternative Work Practice l Liter
BOE Barrels of Oil Equivalent LACT Lease Automatic Custody
BTU British Thermal Unit Transfer
CARB California Air Resources Board LHV Lower Heating Value
CDP Carbon Disclosure Project LPG Liquefied Petroleum Gas
CE Calibration Error MJ Megajoules
CEMS Continuous Emissions MRG Monitoring and Reporting
Monitoring System Guidelines
CER Clean Energy Regulator MPMS Manual of Petroleum
CH4 Methane Measurement Standards
CIP Chemical Injection Pumps NGER National Greenhouse and
CO2 Carbon Dioxide Energy Reporting
CORINAIR The Core Inventory of Air Nm3 Normal Cubic Meters
Emissions in Europe OBQ On-board Quantity
DP Differential Pressure OECD Organization for Economic
DQI Data Quality Indicators Cooperation
EDR Electronic Data Reporting IOGP International Oil and Gas
EF Emission Factor Producers Association
eGGRT Electronic GHG Reporting Tool OSCAR Online System
EIPP Emissions Inventory Comprehensive Activity
Improvement Program Reporting
EJ Exajoules OVA Organic Vapor Analyzer
EMEP European Monitoring and PD Positive Displacement
Evaluation Programme PEI Primary Element Inspection
EPA Environmental Protection PPM Parts Per Million
Agency ROB Remaining on Board
ERCB Energy Resources SAR Second Assessment Report
Conservation Board SI International System of Units
ETS European Union Emissions t tonne
Trading System TJ Terajoule
ETSG ETS Support Group TCD Thermal Conductivity
EU European Union Detector
EVCI Electronic Volume TVA Toxic Vapor Analyzer
Conversion Instrument UNFCCC United Nations Framework
FCCU Fluid Catalytic Cracking Unit Convention on Climate
FID Flame Ionization Detector Change
FFMS Flare Flow Measurement VEF Vessel Experience Factor
System VEFD Vessel Experience Factor on
GC Gas Chromatography Discharging
GHG Greenhouse Gas VEFL Vessel Experience Factor on
GHGRP GHG Reporting Program Loading
Gj Gigajoule VOC Volatile Organic Compounds
GOR Gas-to-Oil Ratio WBCSD World Business Council for
GPA Gas Processors Association Sustainable Development
GWP Global Warming Potential WCI Western Climate Initiative
HHV Higher Heating Value WRI World Resources Institute
Hp Horsepower
INGAA Interstate Natural Gas
Association of America 5
ADDRESSING UNCERTAINTY IN OIL & NATURAL GAS INDUSTRY GREENHOUSE GAS INVENTORIES
INTRODUCTION
The global petroleum and natural gas industry The uncertainties inherent in the data used for
has been active in promoting consistency emission inventories help inform and improve
and harmonization for industry greenhouse understanding of the data’s use. The uncertainty
gas (GHG) emission inventories. Industry of petroleum and natural gas companies’ GHG
associations and their members have emission inventories, or of its quantified emission
been contributing to the development of reductions, is determined largely by uncertainties
guidance for accounting and reporting of in the estimates of the key (largest) contributing
GHG emissions (API/IPIECA/OGP, 2011), and sources. In turn, each of these uncertainties
compiling methodologies that are appropriate depends on the quality and availability of
for estimating GHG emissions from industry sufficient data to estimate emissions. The
operations (API, 2009). This guidance has robustness of industry data disclosures is
been recently augmented with guidelines to receiving increased attention with special
account for reductions associated with GHG emphasis on the need to understand how GHG
projects (API/IPIECA, 2007a). emissions and emission reductions are quantified.
Background
The American Petroleum Institute (API), The 2007 workshop served as the first step in
Concawe (the scientific and technical addressing uncertainty and accuracy issues. In
organization of the European petroleum refining the ensuing industry discussion, a list of priority
industry) and the The global oil and gas industry issues was prepared. This list is comprised of
association for environmental and social issues items that industry experts ought to address
(IPIECA) convened an international workshop in a systematic fashion. As presented in the
on the topic on 16 January 2007, in Brussels, workshop summary report, the issues listed by
Belgium. The goals of this workshop were to: industry members fall into three thematic areas:
• develop an understanding of the relative 1. Measurement methods
importance of the key factors that
contribute to uncertainty; 2. Computational methods
• review and discuss emerging techniques 3. External communications
for quantitative assessment of the
uncertainty and accuracy of GHG emissions Because the industry recognizes the need
estimates; for meeting regulatory mandates and
stakeholders’ expectations, follow-up activities
• identify emission sources and methods will be designed to provide opportunities for
where petroleum and natural gas industry continued dialogue and collaborative activities
efforts are needed to improve accuracy with stakeholders.
and reduce uncertainty to acceptable
levels; and A pilot version of this document was published
• create a prioritized list of topics to be in September 2009 to allow broad review
addressed by the petroleum and natural gas and implementation of guidance provided
industry to minimize emissions estimation by all stakeholders. Comments received,
lessons learned, and new developments that
uncertainty and improve data accuracy. have occurred in the area since then are now
A summary report, as well as all the workshop summarized in this final guideline document.
presentations, is posted on the IPIECA website
(API/IPIECA, 2007b, http://www.ipieca.org/
publication/greenhouse-gas-emissions-
estimation-and-inventories-addressing-
uncertainty-and-accuracy ).
7
IPIECA • API
9
IPIECA • API
1. INTRODUCTION
Policymakers use entity GHG inventories and SECTION FOCUS
reported facility-level GHG emissions to develop
strategies and policies for emission reductions This is an introductory section that
and to track the progress of these policies. introduces some basic concepts and
Both regulatory agencies and corporations rely terms that are the foundation for
on inventories to better understand emission understanding GHG emissions inventory
sources and trends. GHG inventory data are uncertainty. This terminology will be
associated with varying degrees of uncertainty, further expanded throughout the next
and such actual uncertainties have both sections of the document.
technical and policy implications.
The subsections include:
“Uncertainty analysis” has been increasingly
recognized as an important tool for improving • Importance of accurate and reliable
national, sectoral, and corporate inventories of GHG accounting;
GHG emissions and removals (IPCC, 2000). This • Overview of uncertainty
increased attention on accurate inventories terminology;
has resulted in the need to provide guidance • Types of errors;
to industry on technical considerations and
calculation methods for consistent estimation
• Numerical determination of
uncertainty intervals; and
of GHG inventory uncertainty. This typically
would consist of: • Emissions inventory uncertainty
assessment.
• determination of the uncertainties
associated with the individual
measurements and factors used in
constructing the emissions inventory; and Since an understanding of the magnitude and
• propagation and aggregation of these sources of GHG emissions is critical for properly
individual terms to derive uncertainty managing these emissions, using a consistent
intervals (at a pre-designated probability approach can significantly improve industry-
level) for the whole inventory. wide, comparable estimates of emissions and
The extent and scope of such analysis will emission reductions.
depend on the likely uses of this information. Higher-quality GHG data lead to higher
For example, the uncertainty analysis required certainty of emission assessments and improved
for data that are merely used for relative confidence in the data reported. This is true for
ranking or comparison of trends would be national and government assessments, and is also
different from that required to demonstrate important at the entity or facility level. To ensure
attainment of GHG emission limits or progress that a company’s strategies and forward-looking
made towards meeting GHG emission actions are based on the most robust data set and
reduction targets. most appropriate computational methods, it is
important to address the following factors:
a. Comparability: Uncertainty analysis enables
1.1 Importance of accurate and comparison of data trends and between
reliable GHG accounting entities.
Key areas that benefit from reliable GHG a. Consistency: Science-based estimation and
accounting include: measurement methods should include
• focused GHG emissions management; consideration of accuracy and precision.
• reduced business risk and reputation a. Certainty: Emission inventories are estimates;
management; and uncertainty analysis provides a likely range of
• participation in GHG emissions mitigation those estimates.
programs.
a. Confidence: Users of the inventory need to
understand the reliability of the estimated
emissions, especially when they are used for
Comparability, Consistency, Certainty, Confidence policy development, or to target inventory
improvements.
10
ADDRESSING UNCERTAINTY IN OIL & NATURAL GAS INDUSTRY GREENHOUSE GAS INVENTORIES
11
IPIECA • API
12
ADDRESSING UNCERTAINTY IN OIL & NATURAL GAS INDUSTRY GREENHOUSE GAS INVENTORIES
2. ASSESSMENT OF
UNCERTAINTY RELEVANT TO
GREENHOUSE GAS EMISSION SECTION FOCUS
QUANTIFICATION
A myriad of sources contribute to the uncertainty
This section discusses the major
of an emission inventory. Whether at the national,
sources affecting the uncertainty of
entity, or facility level, the ability to quantify
GHG inventories. It moves from general
emissions and understand their associated
concepts to issues that are germane to
uncertainty hinges on two main factors:
GHG inventories in the petroleum and
• readily available data for emission natural gas industry. It also describes
quantification; and factors that could introduce errors into
• knowledge of input parameters for the emission measurements process and
statistical calculation of uncertainty. contribute to the range of uncertainties
of estimated emissions.
The overall range of uncertainty associated
with an entity GHG inventory is primarily The subsections address:
determined by the uncertainty associated
with the largest (“key”) sources of emissions.
• Overview of emissions inventory
uncertainty;
These key emission sources differ by industry
sector and characteristics of company • Sources of Measurement Uncertainty;
operations1. Although very large confidence • Emission Calculation Approaches;
intervals may be associated with the data • Inventory Steps and Data
used to characterize some small sources, the Aggregation; and
overall impact on the range of uncertainty
at the entity, or installation level, may often
• Emissions Inventory and Uncertainty
Assessment in the Petroleum and
be very small. In turn, the confidence interval
Natural Gas Industry.
associated with each individual source depends
on the availability of sufficient data to estimate
emissions, or on the quality of the data in order
to properly account for emission variability.
1 Emissions inventory examples are provided in Chapter 8 of the API Compendium of Greenhouse Gas Emission Methodologies
for the Oil and Natural Gas Industry (API, 2009) and the relative uncertainty of the various sources is addressed in Appendices G
and H of this document.
13
IPIECA • API
Uncertainties due to models (or equations) Uncertainty estimates are tools that enable
depend on the proper application of the inventory preparer to assess the major
quantification methods for the respective contributing factors to the emissions inventory
source categories. These errors typically can and target the most significant ones, i.e., those
be eliminated as far as possible in advance, exhibiting the largest range of uncertainty for
when planning the compilation of an emissions more research and refinement. Table 2-1 provides
inventory, and are often addressed as part of an overview of selected methods recommended
emission inventory assurance processes (API/ by the U.S. Environmental Protection Agency
IPIECA/OGP, 2011). (EPA) for qualitative and quantitative estimation
of the ranges of uncertainty for emissions. The
Parameter uncertainties can be evaluated
uncertainties associated with natural variability
through statistical analysis, measurement
inherent to the emission process and its
equipment precision determinations, and
underlying data can be assessed by the statistical
a variety of expert judgment techniques.
analysis methods introduced in Section 5 and
Quantifying parameter uncertainties,
discussed further in Appendix F.
and then estimating source category
uncertainties based on these parameter If measured parameters uncertainties are
uncertainties, is the focus of this document. unknown, a pedigree matrix approach can
Hence, adhering to appropriate sampling, be used to calculate uncertainties. In this
measurement and estimation procedures – approach, which is a method for structuring
with applicable quality control and quality the expert judgment, qualitative data quality
assurance measures – are all part of a quality assessment results are used to relate data
improvement management system and can quality indicators to uncertainty ranges for
help minimize uncertainties. individual parameters (WRI/WBCSD a, 2011).
Qualitative Discussion
-- Sources of uncertainty are listed and discussed.
Low
-- General direction of bias and relative magnitude of imprecision are given if known.
Subjective Data -- Subjective rankings based on professional judgment are assigned to each emission
Low
Quality Ratings factor or parameter.
Data Attribute Rating -- Numerical values representing relative uncertainty are assigned through objective
Medium
System (DARS) methods.
-- Experts estimate emission distribution parameters (i.e. mean, standard deviation,
and distribution type).
Expert Estimation
-- Simple analytical and graphical techniques are then used to estimate confidence
limits from the assumed distributional data. Medium
Method
-- In the Delphi method, expert judgment is used to estimate uncertainty directly.
-- In the Pedigree Matrix method, experts are used to set up the appropriate matrix for
estimating data quality
-- Emission parameter means and standard deviations are estimated using expert
Propagation of Errors judgment, measurements or other methods.
Medium
Method -- Standard statistical techniques of error propagation typically based upon Taylor’s
series expansions are then used to estimate the composite uncertainty.
-- Monte Carlo, Latin hypercube, bootstrap (resampling), and other numerical methods
are used to estimate directly the central value and confidence intervals of individual
Direct Simulation emission estimates.
High
Method -- In the Monte Carlo method, expert judgment is used to estimate the values of the
distribution parameters prior to performance of the Monte Carlo simulation.
-- Other methods require no such assumptions.
-- Direct or indirect field measurements of emissions are used to compute emissions
and emissions uncertainty directly.
Direct or Indirect
Measurement
-- Methods include direct measurement such as stack sampling and indirect
High
measurement such as tracer studies.
(Validation) Method
-- These methods also provide data for validating emission estimates and emission
models.
a
Extracted from Table 4.1-1 of the Emissions Inventory Improvement Program (EIIP), Chapter IV: “Evaluating the
Uncertainties of Emission Estimates,” EPA, Research Triangle Park, NC, July 1996
14
ADDRESSING UNCERTAINTY IN OIL & NATURAL GAS INDUSTRY GREENHOUSE GAS INVENTORIES
15
IPIECA • API
Data Quality -- Comparing known input values with their measured or computed results can provide an estimate of the data
acquisition uncertainty.
-- When it is not possible to do this in practice, it is advisable to evaluate potential individual and aggregated
errors when assessing uncertainty.
– or –
16
ADDRESSING UNCERTAINTY IN OIL & NATURAL GAS INDUSTRY GREENHOUSE GAS INVENTORIES
The four basic approaches for estimating emissions be used broadly across many industry source
and their potential error sources are summarized in categories and operations. The published
Table 2-3. Each contributor to uncertainty should emission factors are typically accompanied by
be assessed independently and then aggregated a description of the group of processes and the
in the final analysis, as discussed below. conditions they represent. Authoritative factors
are generally published by the EPA in AP-42
In the absence of direct emission measurements,
(EPA, 1995 and further updates), or by the EU
the use of emission factors has traditionally
EMEP/CORINAIR Emission Inventory Guidebook
been the simple, low cost method of choice for
(EEA, 2013). The IPCC has also launched a new
estimating emissions. Over the years, the EPA
Emissions Factors Database (IPCC, EFDB) for use
and other emission factor repository databases
with GHG emission calculations.
have provided average emission factors that can
Table 2-3. Common emission estimation approaches and their error sources
Emissions Estimation
Description
Approach
An emission factor relates the rate of emission of a specific compound to an activity
rate associated with its release.
The errors associated with the use of emissions factors are attributed to:
-- Variations in operating conditions during the collection of emission and activity
data that were used in the development of the emission factors.
Emissions Factors -- Variability in emissions that arises from differences in operating conditions
among different facilities where the factors are used.
-- Difference between the actual composition of the stream to which the
emission factor is applied and the default composition on which the emission
factor is based.
-- Uncertainty that is due to measurement errors, systematic errors, and random
sampling errors.
The technique involves continuously measuring flow and concentrations of species
directly emitted into the atmosphere from a specific source, such as a stack. It is
accomplished by placing applicable monitors at the source.
Uncertainty associated with CEMS generally includes:
Continuous -- Stack flow measurements.
Emissions Monitoring -- Concentration measurements.
System (CEMS) -- Measurements of stack temperature and pressure required to correct emissions
to standard conditions.
The error associated with these determinations varies for different compounds.
CEMS will not necessarily produce emission estimates with lower uncertainties than
alternative methods especially since CEMS are not available for monitoring all GHG
emissions. CEMS may not be practical for a large number of emission sources.
This technique involves either extracting a sample or placing a monitor at a
source, followed by analysis to characterize the emitted species. The measurement
campaign is typically limited to a specified number of hours, and the average
emission rates calculated are used to estimate total annual emissions. For
characterizing GHG emissions over a longer period of time (such as a year), periodic
sampling and analysis can be used to determine emission variability.
Source Testing The errors associated with these measurements are due to:
-- Sporadic testing of limited duration that is used and extrapolated to
characterize emissions over a longer period, i.e., a year, not allowing for a
robust estimate of variability.
-- Testing methods used are either improperly calibrated or do not have Version
2, September 2014 2-7 sufficient sensitivity to enable characterization of the
full gamut of emission rates.
This technique is based on the material balance or total quantity of material used,
e.g., fuel flow into a combustion device and its carbon speciation.
Material Balance The errors associated with this method include:
-- Improper calibration and measurement of applicable flow devices.
-- Inaccurate determination of material stream composition.
17
IPIECA • API
2.4 Inventory steps and data total emissions for each of the GHG species,
aggregation along with the global warming potential-
weighted sum of these emissions (also known
When developing GHG emission inventories, as the CO2e emission).
emission estimates are obtained from many
intermediate and independent results, each The overall uncertainty range for each GHG
of which is calculated from a separate set of species, and CO2e, should also be reported with
data that is characterized by a different range the total emissions in the format of:
of uncertainties. The compilation of an entity- Emissions = Average Value ± %
wide GHG emissions inventory typically follows (at the 95% confidence limit).
a sequence of steps:
(Equation 2-5)
a. Establishing boundaries – Where the
organizational and operational boundaries
are defined (for a first-time inventory), 2.5 Emissions inventory and
or examined (for recurring cycles), this uncertainty assessment in the
step will be largely dictated by local petroleum and natural gas industry
requirements or corporate policies. It
might involve facility-by-facility assessment Inventorying of GHG emissions by entities is
prior to aggregation, or it could use other a ‘bottom-up’ summation of emissions from
pertinent entity indicators and information. individual sources (or emissions from the total
consumption of different fuel types) at a report-
b. Collecting and inputting data – Where ing unit to create an inventory for that report-
the activity data are collected and archived ing unit. Emissions from disparate reporting
based on the boundaries established units may be aggregated to create an entity,
above. The data are then incorporated into or corporate, inventory. Reporting units may
appropriate calculation tools for emission be defined by the reporting entity to represent
calculations. The level of ‘granularity’ of logical groupings of activities and assets, or
the data collected and the details of the could be mandated by governments as part
calculation methods are dictated by local of GHG reporting regulations. Developing
requirements with industry guidance (API these GHG inventories entails both the proper
Compendium) as a resource to provide accounting of activity levels for operations that
relevant technical details. may lead to GHG emissions as well as the con-
c. Validating data compiled – Where various version of these activity levels to quantitative
techniques are used to compare the new GHG emissions using proper measurement and
data with earlier versions (if available) to calculation methods.
identify potential large errors. These errors
could include: either large changes or The petroleum and natural gas industry en-
unchanged activity data for given facilities; compasses a wide variety of activities, rang-
operations that are not accounted for; lack ing from the exploration and production of
of supporting data for measurements or petroleum and natural gas to the delivery of
emission factors used; erroneous units or products to consumers. As part of defining the
unit conversions among others. scope of GHG inventories, industry guidelines
(IPIECA/API/OGP, 2011) are available to assist
d. Assessing data uncertainty – Where the companies in determining which emissions
confidence intervals associated with the related to their activities should be included
data available for each of the emission within the organizational boundaries they
sources are characterized independently, have established. In all cases, mandatory GHG
as discussed later in this document. The reporting requirements would take precedence
uncertainty information could be based over industry guidelines, as may be the case in
on documentation of data repositories different jurisdictions.
(API, 2009), expert judgment, or on
measurements conducted during the For a complex and dynamic sector such as
inventory year. the oil and natural gas industry, a variety of
e. Finalizing the inventory – Where the methods are applicable to quantifying GHG
quality-checked and validated data are emissions, ranging from simple activity mea-
aggregated for reporting based on company surements multiplied by applicable emission
policy or local requirements. The preferable factors to more sophisticated quantification
way of reporting the results is in terms of the algorithms. Advanced engineering estimation
18
ADDRESSING UNCERTAINTY IN OIL & NATURAL GAS INDUSTRY GREENHOUSE GAS INVENTORIES
methods consist of an integrated approach Figure 2-1. Hierarchy of emission estimation approaches
that relies on the use of emissions and physical
activity factors in conjunction with generic Types of Approaches Hierarchy
process simulation models, source-specific
models, and species profiles databases. Each of
these estimation methods will have a different Published emission factors
level of uncertainty. For example, combustion
emissions are dependent on the amount of fuel
Equipment manufacturer
consumed. This can be estimated based on the
emission factors
hours of equipment operation, equipment rating, Improved accuracy
and thermal efficiency. The uncertainty associated
with this approach differs from the uncertainty Engineering calculations Additional data requirements
associated with directly measuring fuel use. Higher cost
Monitoring over a range
Selection of an appropriate emission estima-
of conditions and deriving
tion method must consider:
emission factors
• The objectives of the emission inventory;
• The contribution of the emission source Periodic monitoring of
to the overall accuracy of the emission
emissions or parameters for
inventory; and
calculating emissions
• The costs and practicality of the emission
estimation method.
Continuous emissions* or
Figure 2-1 provides a general hierarchy of accu-
parameters monitoring
racy associated with common emission estima-
tion approaches. Accuracy is generally improved *Continuous emissions monitoring may not be directly
as emission estimation progresses from default applicable to certain greenhouse gases or to all emission
emission factors to more direct estimation or sources.
measurement methods. For example, measuring from key emission sources associated with
CH4 and CO2 emissions from every dehydrator the natural gas storage and transmission
in a natural gas production field may produce sector of the industry (INGAA, 2005) and
a highly accurate emission rate, but at great the American Gas Association (AGA), which
expense for a source that has minimal contribu- has also published specific guidelines for
tion to the GHG inventory of the field and can estimating GHG emissions from gas distribution
be estimated using an engineering modeling operations (AGA, 2008).
approach, such as GLYCalc™ (GTI, 2000).
Industry guidelines also served as major input
to the formulation of regulations governing
The third edition of the API Compendium (API, the mandatory reporting of GHG, such as the
2009) provides an extensive compilation and U.S. GHG Reporting Program (EPA, 2009), the
tabulation of methods used by companies Western Climate Initiative (WCI, 2009), and the
in all segments of the petroleum and natural California Air Resources Board (CARB, 2012).
gas industry for consistent calculation of GHG
emissions. The API Compendium provides The categories of sources that contribute to
a wide range of emission factors and other overall GHG emissions as classified in the API
emission estimation methods that are directly Compendium (API, 2009) are:
applicable to all sectors of industry operations. • Stationary combustion sources – linked to
Summary tables and decision trees are emissions resulting from the combustion of
provided to illustrate the variety of available fuels in boilers, furnaces, burners, heaters,
emission estimation options and the associated and stationary turbines and engines. This
considerations. It also lists the uncertainty category may also include the combustion
ranges (at the 95% confidence level) associated of waste gases-or emergency releases-in
with multiple case study examples featured in incinerators and flares.
Section 8 of the 2009 API Compendium.
• Mobile combustion sources – linked to
Related industry guidelines include those emissions resulting from combustion
of the Interstate Natural Gas Association of of fuels in ships, barges, trains, trucks,
America (INGAA), which provides supplemental automobiles, aircrafts and other off-road
guidance for estimating GHG emissions devices such as drilling rigs.
19
IPIECA • API
20
ADDRESSING UNCERTAINTY IN OIL & NATURAL GAS INDUSTRY GREENHOUSE GAS INVENTORIES
21
IPIECA • API
3. OVERVIEW OF
SECTION FOCUS
MEASUREMENT PRACTICES FOR
ESTIMATING GREENHOUSE GAS
EMISSIONS This section provides an overview
of measurement practices focusing
In defining aggregated uncertainty of on gas flow measurements and the
measurement ensembles used for developing determination of carbon content or
emission inventories, the uncertainty for each the heating values of combusted
measurement stream must be assessed in a fuels. The section recognizes industry-
way that is applicable to that measurement recommended practices and standards
method and its implementation in practice. including those that traditionally applied
Random errors could be a major factor in to “custody transfer” of fuel products and
the uncertainty of an individual observation; new ones that have been more recently
however, their contribution to the overall developed for use in an emissions
emission inventory diminishes as more estimation context.
measurements are obtained during the
reporting period. This section goes on to link these
practices to the acquisition of GHG
Note: Random errors tend to average out emissions and related activity data,
during the year, whereas systematic errors (or representing a move from reliance on
measurement bias) become more important and available data or engineering judgment.
tend to accumulate rather than diminish over The subsections address:
longer periods of time such as a year.
• Flow Measurement Practices;
In fact, determining the true value of any
measured variable is not practical due to the
• Flow Measurements for GHG
Inventories;
limitations of measurement equipment and
procedures, and the possibility of human error. • Sampling and Analysis for
Hence, industry measurement procedures and Quantifying GHG Emissions;
standards have been developed to emphasize • Carbon Content Measurement
practices that lead to collecting better quality Practices;
data, especially for critical measurements. • Heat Content Determination;
Industry uses standards from several different
standards-developing organizations, resulting
• Venting and Fugitive Emissions
Measurements; and
in equivalent measurements, based on the
scope of the standard, company preference, • Laboratory Management System.
and type of devices used. Consensus industry
standards, such as those developed by ANSI,
API, ASTM, ISO and other standard-setting of specific measurement standards (and their
organizations, have rigor in their development respective editions) from several standards-
process, and measurement standards are setting organizations, which could be used
reviewed at least every five years to ensure to support the calculation of GHG emissions.
that standards are in step with technological It is up to the user of these measurement
changes and advancements. Most, if not all, standards to reference the specific standards/
of the measurement standards are developed editions used for a given measurement, and
for measurements associated with ‘custody to incorporate the updated measurement
transfer’ and to define quantities that are procedures, as applicable.
essential for robust financial transactions. Custody transfer measurements are typically
API publishes one of the more comprehensive are expected to meet a set of performance
sets of custody transfer measurement specifications. For other measurements,
standards, but it is neither complete nor the such as those performed to support the
only widely recognized source for such industry development of a GHG emissions inventory,
practices. API’s MPMS (Manual of Petroleum data quality objectives should be established
Measurement Standards) includes over 140 prior to initiating any data collection to
titles, and API publishes approximately eight ensure that the uncertainty ranges of the
new or revised measurement standards each measured quantities are consistent with the
year. Appendix B presents a comprehensive list intended use of the data. Throughout this
section, we provide references and describe
22
ADDRESSING UNCERTAINTY IN OIL & NATURAL GAS INDUSTRY GREENHOUSE GAS INVENTORIES
Industry is now more fully internalizing the 3.1.1 Measurements by Orifice Meters
potential impact of measurement errors
and bias through the full chain of emission Orifice meters are by far the most prevalent
calculations including measurement flow meter type used in the petroleum
equipment, software calculations, and natural gas industry, and are used for
simulation models, and the limitations of metering products during “custody transfer”
reliance on existing emission factors on the as well as for process control and internal
uncertainty range of resultant GHG emissions accounting. These flow meters are also used
inventories. Considerations of the need for to account for fuel volumes when estimating
more representative measurement data, CO2 emissions. Flow meters are designed
and the assessment of equipment design for long-term reliability and ruggedness
and age, are gaining more prominence in under a variety of component mixtures and
the process of assembling high certainty conditions that are essential for consistent
emission inventories. This section focuses on fluid blending and processing. Although for
measurements that are pertinent to improved this type of meter, temperature and pressure
characterization of calibrations can be done while the units
are operating; they generally have limited
combusted fuels and the quantities used, access to direct orifice plate inspections and
but might not be directly applicable to maintenance outside of planned shutdown
measurement of leakages and fugitive (‘turnaround’) cycles.
emissions. The measurement practices
highlighted here represent an initial step in Recommended practices for the installation,
what could end up being a multi-year effort to calibration, and calculation of flows for
improve measurements of GHGs and quantify these custody meters are provided in
the activities that cause their emissions. Section 3 of Chapter 14 of API’s MPMS (API,
23
IPIECA • API
2013). This standard was developed through • Part 1: General equations and uncertainty
a collaborative effort by members of API, guidelines.
AGA, and the Processing Gas Association, • Part 2:Specifications and installation
with contributions from the Canadian Gas requirements;
Association, American Chemistry Council, the
European Union, Norway, Japan, and others.
• Part 3: Natural gas applications.
It is designed to ensure global consistency for • Part 4: Background, development,
petroleum and natural gas transactions and is implementation procedures, and
recognized by the Version 2, September 2014 3-3 subroutine documentation.The standard
recognizes that many factors contribute to
American National Standards Institute (ANSI) as overall measurement uncertainty associated
an American National Standard. The four-part with many metering applications, as
standard for square-edged, concentric orifice summarized in Exibit 3-1
meters consists of:
24
ADDRESSING UNCERTAINTY IN OIL & NATURAL GAS INDUSTRY GREENHOUSE GAS INVENTORIES
All these factors should be assessed when from molecular weights approaching that of
estimating the overall range of uncertainty for hydrogen to molecular weights of C5+.
flow measurements using thin plate, concentric,
As with other flow measurements, the accurate
square-edged metering systems.
determination of flow to flares is dependent on
In the reapproved 2012 version of the standard, many parameters such as the ability to predict –
several changes were incorporated to reduce the or measure – the quantity of gas flared, mixture
uncertainty attributable to installation effects composition, pressure, temperature, and/
and to improve the rigor of the flow calculation or density. For example, gas-to-oil ratio (GOR)
routines. The revised standard recognizes the measurements can be used to estimate the
lead time necessary for upgrading existing volume of gas flared for isolated oil production
installations, and leaves this lead time to the operations, or engineering estimates based on
discretion of facility operators and their data purge rates or operating logs can be used to
quality targets for flow measurement data. estimate flare volumes for refineries or gas plants.
However, it should be recognized that if orifice The accuracy of measurements associated with
meter installations are not upgraded to conform highly variable flare gas mixtures will depend
to the new recommendations, measurement largely on the meter technology type and the
bias error may occur. This bias might be due ability of the flare flow measurement system
to improper upper and lower distances from (FFMS) to achieve the targeted response time
bends and points of flow turbulence that and analytical accuracy levels. Exhibit 3-2
might lead to inadequate flow conditioning below lists the basic steps needed to conduct
prior to measurement. Additionally, even a simplified analysis for determining the
without changing equipment installations, the uncertainty ranges of a given flow measurement
standard recommends adopting new calculation system. The actual approaches for the required
procedures and techniques (explained in calculations are provided in Section 5, and
Part 1 and 3 of the standard) that represent Appendix F provides examples demonstrating
significant improvements over the previously how to apply these methods.
adopted approach. It is important to note
that the expected uncertainty ranges for flow
measurements quoted in Part 1 of the reaffirmed EXHIBIT 3-2: STEPS FOR SIMPLIFIED
standard may differ from those obtained in UNCERTAINTY ANALYSIS
practice when the equipment installation differs.
a) Determination of the equation that defines the
meter output
3.1.2. Measurement of flow to flares • Governing equations are those applicable for the
meter technology type used.
The measurement of flow to flares is distinctly
different than other flow measurements. b) Determination of the combined sensitivity
Flares are designed as safety relief systems coefficient
and typically are capable of handling highly • Numerical values of uncertainty are associated with
variable flow rates of widely varying gas pressure, temperature, and composition.
compositions. Therefore, some of the practices • Meter accuracy is estimated from calibration data,
that are generally applicable to custody transfer pipe size, or other installation effects.
or process control flows have to be modified • Sensitivity coefficients are obtained by dividing the
when addressing flows to flares. API published calculated percent uncertainty for an input variable
a measurement standard addressing gas or by the percent change in that input variable.
vapour flare flow measurements, which also c) Derive the combined uncertainty range (extended
includes cautionary details about the effects standard deviation)
of fouling (due to entrained liquid droplets, • Combined uncertainty is calculated by summing
aerosol mists, or other contamination) on the the square of the errors for pressure, temperature,
measurement (API MPMS, July 2007). composition, meter calibration, and installation
effects.
Most flare headers are designed to operate
during both non-upset conditions at near
atmospheric pressure and ambient temperature,
and during flare episodes, at a wide range of
pressure, temperature, and flow velocities.
During such episodes, flare gas compositions
are also highly variable and could range
25
IPIECA • API
Variability in flare composition may also Each component in the measurement system
be a significant factor in determining the will exhibit variations that contribute to
measurement uncertainty of an FFMS. variations in the process result. Variability
Knowledge of flare composition may have a may be due to procedure(s), operator(s),
major effect on the calculation of the actual environment, instrument(s) and maintenance.
volume, standard volume, or mass measured by Consequently, repeated results from the same
the meters. For example: measurement process may vary, and there is
always the potential that many small variations
• For a differential pressure meter: the
will affect the final results.
output is a function of the square root of
the flare gas density. Therefore, it is advisable to analyze the overall
• For thermal flow meters: knowledge uncertainty of the flow measurement process
of the actual volume (or standard to affirm its adequacy for the intended
volume) requires consideration of application. As indicated in Exhibit 3-2, a few
the compositional effect on thermal key steps should be followed to determine
conductivity and dynamic viscosity. the combined uncertainty associated with
individual flow measurements. An expanded
• For ultrasonic flow meters: sound speed discussion of factors to be considered
is a function of gas composition.
when evaluating the uncertainty of flow
Converting volume flow to mass flow requires measurements used for GHG emission
knowledge of gas composition in order calculations is included in Section 3.2.
to derive gas density. When an analyzer is
incorporated in the measurement system to
correct for composition, care must be taken to 3.1.4 Uncertainty specifications for “custody
ensure that the response time of the system is transfer” measurements
short compared to the upset flow event during
flaring to ensure representative measurement “Custody transfer” measurements are defined
during actual flaring. as measurements that provide quantity and
quality information, which can be used as
An example of how all of these elements would the basis for a change in ownership and/or
be combined into the overall measurement a change in responsibility for materials. In
uncertainty is provided in Table 3-1. most petroleum and natural gas producing
jurisdictions around the world, national
regulations and directives have emerged
3.1.3 Flow measurements uncertainty to specify requirements for the expected
analysis accuracy and uncertainty ranges associated
The result of a measurement process is the with “custody transfer” and other critical
determined quantity of the parameters measurements. For such precise metering
measured. For flow measurements the applications, the flow meters and adjacent
quantities may include pressure, temperature piping used in the measurement system
and volumes, (i.e., multiple components) are expected to meet the requirements of
contribute to the accuracy, or uncertainty, of the relevant, preferably the most stringent,
the final result. specifications of the API and ISO standards
that are cited in many national regulations.
Table 3-1. Example of flare flow measurement systems combined
uncertainty a One such example of measurement
requirements promulgated for petroleum
Combined
Variable Sensitivity and Error (S x U95)2 and natural gas operations is Directive
(S x U95) 017 of the Alberta Energy Regulator (AER,
2013), which spells out the measurement
Pressure 2.0% 4.00
requirements for upstream oil and gas
Temperature 0.1% 0.01 operations in Alberta, Canada:
Flare Composition 2.0% 4.00
• The standards are stated as “maximum
Meter error (calibration 1.4% 1.96 uncertainty of monthly volume” and/or
Installation effects 0.5% 0.25 “single point measurement uncertainty,”
Sum of squares 10.22 as listed in Table 3-2 below for petroleum
(oil) Systems and Natural Gas Systems
Square root of sum of square 3.2%
measurements, respectively.
a Table 4 from API MPMS Section 14.10 (API, 2007)
26
ADDRESSING UNCERTAINTY IN OIL & NATURAL GAS INDUSTRY GREENHOUSE GAS INVENTORIES
Note: Extracted from Section 1.8.1 and 1.8.2, respectively, ERCB Directive 017, September 2012
a
27
IPIECA • API
29
IPIECA • API
30
ADDRESSING UNCERTAINTY IN OIL & NATURAL GAS INDUSTRY GREENHOUSE GAS INVENTORIES
31
IPIECA • API
32
ADDRESSING UNCERTAINTY IN OIL & NATURAL GAS INDUSTRY GREENHOUSE GAS INVENTORIES
EXHIBIT 3-4: KEY FACTORS AFFECTING GAS SAMPLING AND ANALYSIS UNCERTAINTY
a) Inappropriate sampling techniques or equipment
• All sampling methods would require the use of a sample container for transporting the
sample from the field location to the laboratory.
• Whenever practical, samples should be collected on a flow proportional or flow weighted
basis, since spot samples – by their nature – may not fully represent a gas stream of varying
composition.
• Gaseous samples of interest are a mixture of organic and inorganic gases, and their integrity
will be compromised.
b) Inappropriate sample conditioning and handling
• Bias could be introduced if any components of a sample are either depleted or augmented
during the sampling, transport, or laboratory handling phases prior to analyses.
• Condensation and revaporization of hydrocarbons can cause significant distortions in the
gas sample.
• Care should be taken to sample above the hydrocarbon dew point and/or to prevent
retrograde condensation when pressure is reduced during sampling.
c) Collection of samples from non-representative locations and/or under non-
representative operating conditions
• Sampling systems that are used in conjunction with on-line analyzers, such as
chromatographs or gravitometers, typically are designed to extract, condition, and deliver a
representative sample to the analyzer.
• Sampling lines are kept as short as possible in conjunction with proper heating and
insulation to avoid condensation.
• The flow rate of the sampling system is adjusted to allow for close to real time data, while at
the same time not increasing the flow to a level that might lead to turbulence.
d) Inappropriate analytical methods
• The threshold sensitivity of the analytical methods used are typically those that are well
3.3.2 Quantifying sampling precision same time. The determination of sampling bias,
or the difference between the mean of several
Quantifying sampling precision requires that
measurements and the true value, is more
primary samples be collected according to
difficult. Biases could arise from several causes
a defined protocol, but randomized in some
such as sample fouling, inappropriate handling,
way for each sample (in either space or time).
or unrepresentative sampling. The “true” value
For example, a pre-selected percentage of
for the concentration of unknown species in a
the total number of samples can be collected
sample is never known since it is impossible to
in duplicate and the repeatability of the
construct a true reference laboratory standard
measurement determined from these duplicate
for an unknown mixture of gases.
samples. Additionally, duplicate gaseous
samples can also be analyzed in duplicate Sampling uncertainty is a relatively new
and thus a full record of both sampling and concept whose importance is slowly
analysis system variations can be obtained. beginning to be recognized. Recent research
In this case, the sampling component of the has shown that sampling uncertainty is
variance represents the natural variability of often far greater than analytical uncertainty.
the sampling target as well as any errors in the Therefore, combining sampling and analytical
sample collection and preparation. uncertainties to provide an estimate of
measurement uncertainty is an important
When in situ measurement techniques are
component of quantifying the overall
used (e.g., infrared gas analyzer), both the
uncertainty of GHG estimations that rely on
sampling and analysis are addressed at the
sampling gaseous fuels of varying composition.
33
IPIECA • API
3.4 Carbon content measurement by the same operator with the same
practices apparatus under constant operating
conditions on identical test materials
Different types of gas chromatography (GC)
systems are normally used to analyze the
• Reproducibility is the difference between
two results obtained by different operators
carbon content of gaseous streams. The GC
in different laboratories on identical test
systems might be laboratory based or set up
materials.
as an online device for automated collection
of samples and their analysis. The systems When using a natural gas measurement
typically are set up to analyze the individual method for refinery fuel gas, care should be
components in the sampled gas and provide taken to ensure that the range of compositions
detailed reports of properties including of individual components is within the bounds
composition, calorific value, and density. specified by the method. Table 3-4 lists selected
commonly used laboratory measurement
The results of the determination of individual methods that are frequently used for the
or groups of carbon-containing species are determination of fuel carbon content.
then used to assess the total emissions of CO2 Additional details about these methods and
upon combustion of such a fuel. Several key further discussions of their main features are
considerations include: provided in Appendix D.
• If the method is capable of determining
CO2 content with the rest of the carbon
containing species, no further correction 3.4.2 On-line measurements
of the carbon content data is required to On-line determination of fluid stream
properly account for all CO2 emissions. compositions is quite challenging due to
• If the method is set up to provide possible variations in these compositions. This
information only on hydrocarbon species, is especially notable for self-generated fuel gas
the CO2 content should be obtained by an such as refinery fuel gas or other processing
independent measurement and added to plant gas. Conversely, for commercial products
the fuel carbon content data. such as natural gas, liquid fuels, coal and coke, or
for the analysis of associated gas in exploration
• If the method is capable of a quantitative and production operations, the challenges are
determination of CH4 content, these data more related to the ability to analyze multiple
can be used separately for calculating streams rapidly and ascertain that they all are
evaporative and processing leaks along within a desired property range.
with venting losses.
Instrumentation in this field has been
• All carbon content measurement data developed to provide a measurement of stream
should be used in conjunction with the components in order to achieve optimum
applicable fuel flow measurements when control and assure product quality. The
calculating emissions. configurations of such analyzers are customized
to accommodate typical site parameters and
operating practices. Most such analyzers
3.4.1 Laboratory-based measurements
are designed with ASTM and ISO standards
Several ASTM and ISO methods are available in mind, and their calibration routines are
for determining the composition and carbon designed to provide both reported data and
content in the natural gas range as well as for its associated uncertainty. As mentioned
liquid and solid fuels. The methods usually previously, the two primary applications are for
document measurement precision in terms refinery gas analyzers and natural gas analyzers,
of repeatability and reproducibility of the and these are discussed briefly below.
measurement when comparing multiple
sets of measurements. Although the terms • Refinery Gas Analyzer: Refinery gas
“repeatability” and “reproducibility” are samples are delivered to the sample inlet
applicable to all measurement situations (see of the GC after passing through a sample
Appendix A for terminology and definitions), conditioning system that selectively
they are used here in the context in which they removes any liquid fractions and
are defined in ASTM standards: particulate matter from the sample. This
ensures that only the gas phase sample
• Repeatability is defined as the difference is delivered to the analyzer. An internal
between two successive results obtained vacuum pump draws this conditioned
34
ADDRESSING UNCERTAINTY IN OIL & NATURAL GAS INDUSTRY GREENHOUSE GAS INVENTORIES
35
IPIECA • API
sample into each injector, which then injects (API, 3rd edition, August 2009). The carbon
the mixture onto each of the columns for content of a mixture is the weighted
analysis. A complete analysis of hydrogen, average of the individual component carbon
saturated and olefinic hydrocarbons (C1-C5, contents and can be calculated from the
and C6+ grouped peaks), plus fixed gases compositional analysis of the mixture (see
(O2, N2, CO, and CO2) typically is performed. API Compendium Table 3-6). Hydrocarbon
Precise “retention times” information and gas chromatographic analysis is generally
component areas translate into accurate reported in mole fraction or mole percent,
component identification and quantification and hydrocarbon liquid chromatographic
of the relative magnitude of individual analysis is generally reported in mass fraction
components present in refinery gas. or mass percent; therefore, calculations
• Natural Gas Analyzer: These types of for both types are included in API’s TR
analyzers are applicable to natural gas 2572. Molecular weights of hydrocarbon
samples from wellhead to pipeline- components can be found in the Table of
quality gas. Samples are introduced Physical Properties for Hydrocarbons and Other
using sample cylinders, Tedlar bags, or Compounds of Interest to the Natural Gas
by direct connection to the pipeline or Industry (GPA, 2145).
wellhead sampling points. Usually, two Carbon content can also be estimated from
chromatographic modules are used other process or fuel supply properties
to quickly separate and measure the such as heating values that are correlated
individual components in natural gas. to the fuel type and carbon content. The
The analyzer separates and measures API Compendium (3rd edition, August 2009)
the permanent gases and hydrocarbons summarizes a number of these correlations in
present via an optimized, dual-channel its Tables 3-8 and 3-9, which provide carbon
portable gas chromatograph. Wellhead content and heating value information
samples may often contain significant for petroleum products and natural gas
amounts of H2S. Many instruments heating value ranges, respectively. Use of
take this into account and there are no these correlation factors works well for
interferences, which mean that H2S can be compositions with limited variability or for
measured from 50 parts per million (ppm) fuel supplies where the heating value is
to 30-mole%. continuously monitored or controlled.
For process or fuel supplies with widely
3.4.3 Carbon content calculation varying compositions or that are subject to
periodic upset conditions, the frequency,
API Technical Report (API, 2013) provides
duration and carbon content of this
guidance and a methodology for determining
variability needs to be characterized to
the carbon content of hydrocarbon-based
accurately calculate the average carbon
petroleum and petrochemical products, and
content. Section 3.6 below provides further
the uncertainty of the average carbon content
details on heat content determination.
as calculated from multiple samples taken
during a sampling period. This method is API TR 2572 provides references and
intended to make use of industry-accepted supplemental information on applicable
mixture property data and test methods with industry practices based on published
no new or modified test methods introduced resources, existing industry standards,
in this document. The method is applicable to industry-accepted physical constants or
carbon-content-based reporting or trading for properties of hydrocarbons for measurement,
all gaseous and liquid hydrocarbons. sampling, sampling frequency, and analysis of
hydrocarbon samples. It provides guidance on
It is a companion to API TR 2571 on Fuel
estimating carbon content sampling period
Gas Measurement (see Section 3.2.2), which
uncertainty based on multiple periodic samples
provides guidance on measuring the volume/
taken at intervals that are independent of the
mass of process fuel gas or feedstock.
process or fuel supply operation.
Equations for the calculation of carbon
To determine the number of samples required
content for a pure component and a product
to meet a certain average carbon content
analysis typically are based on the definitions
uncertainty, one can use Equation 3-1 below:
found in Section 4.3 of the API Compendium
36
ADDRESSING UNCERTAINTY IN OIL & NATURAL GAS INDUSTRY GREENHOUSE GAS INVENTORIES
2
σ
k 95 × ×100
CCaverage
Number of Sampleas =
TargetPercentCCuncertainty
(Equation 3-1)
Where:
Number of samples = the number of samples to be collected
k95 = the 95 % confidence coverage factor (typically assumed to
be 3)
σ = the carbon content standard deviation of the samples
analyzed
CC = carbon content
Target Percent = the target percent uncertainty for the reporting period
CCUncertaint
This equation could be reversed to derive the level 3.5 Heat content determination
of uncertainty given a specific sampling frequency.
The heating value or calorific value of a
According to the discussion in Annex B of API substance is the amount of heat released during
TR 2572, the uncertainty of reporting period the combustion of a specified amount of the
fuel carbon content (i.e., its 95 % confidence substance. The calorific value is a characteristic
interval) is dependent on the number of for each substance and is measured in units
multiple fuel samples used. Increasing the of energy per unit of the substance, such as:
number of samples reduces the uncertainty of kcal/kg, kJ/kg, J/mole, Btu/m³. The heat of
the reporting period average carbon content combustion for fuels is expressed as:
by the square root of the number of samples.
• HHV: higher heating value (or gross calorific
For example, if the average annual carbon value or gross energy or upper heating value)
content is determined from monthly fuel is determined by bringing all the products
sampling and analysis (12 spot samples) and of combustion back to the original pre-
the standard deviation of these 12 samples combustion temperature, and in particular
is 3.9%, the 95% confidence interval for condensing any water vapor produced.
annually reported carbon content will be: • LHV: lower heating value (or net calorific
± (3×3.9% / √¯12), or ± 3.4 %. However, if value) is determined by subtracting the
quarterly samples are used (four samples heat of vaporization of the water vapor
per year with the same sampling standard from the higher heating value and treating
deviation of 3.9%), the 95% confidence any water formed as a vapor. The energy
interval for the annually reported carbon required to vaporize the water therefore is
content would be + 5.9%. not realized as heat.
37
IPIECA • API
In the petroleum and natural gas industry, the The ASTM approach recognizes that the
two most prevalent modes of determining calorific value of a mixture, such as in refinery
heating values of gaseous fuels are either fuel gas systems, is a function of the mol%
by measuring it directly, which can be composition of the individual components of
accomplished either by stoichiometric that mixture. For refinery fuel gas, this mixture
combustion or by calorimetric techniques, or could contain both carbon containing species,
using computational methods that are based on which upon combustion will contribute
standardized calculation procedures using gas directly to CO2 emissions, as well as non-carbon
sample composition data. A brief summary of containing species. Hydrogen, for example, is
these two types of practices is provided below. an important contributor to refinery fuel gas
heating values but it does not contribute to CO2
emissions when combusted.
3.5.1 Direct measurements
To implement this practice, the user would first
The heating value indicates the amount determine the molar composition of the gas in
of energy that can be obtained as heat by accordance with any applicable ASTM or GPA
burning a unit of gas. The heating values of a method. For a precise calculation, at least 98% of
gas depend not only upon the temperature the sample constituents should be determined
and pressure, but also upon the degree of as individual components, with no more than
saturation with water vapor. a total of 2% in terms of groups of components
(e.g. butanes, pentanes, hexanes, and so forth).
As mentioned above, general practices for
determining fuel gas heating values rely on An ideal combustion reaction for fossil fuels
either calorimetric techniques or stoichiometric (that may contain hydrogen) in the ideal gas
combustion practices. Table 3-5 provides a state can be generally represented as in the
listing and a brief description of some selected reaction depicted in Eq. 3-2:
methods for heating value determination for
gaseous, liquid, and solid fossil fuels. Further b c b
discussion of the main features of these
Ca H b + cH 2 + a+ + O 2 = aCO 2 + + c H 2 O
4 2 2
methods can be found in Appendix C.
(Equation 3-2)
The ideal net heating value is the heating value
3.5.2 Computational methods
that is observed when all the water remains in
Heating value may be determined from gas the ideal gas state, while the ideal gross heating
compositional analysis in accordance with a value is observed when all the water formed by
standard practice established by the ASTM for the reaction condenses to liquid. The difference
calculating heating values for natural gas and between them is the enthalpy of vaporization of
similar mixtures from compositional analysis. the water formed during the combustion process.
ASTM D3588-98 (Reapproved 2003), is the Therefore, the ideal gross heating value for a
standard recommended practice for calculating mixture can be expressed as shown in Eq. 3-3:
heating values, compressibility factors, and relative
densities of gaseous fuels (ASTM D3588-98). This n
practice covers procedures for calculating these
quantities at base conditions (14.696 psia and
∑x
j =1
j × M j × H mj
60°F or15.6°C) for natural gas mixtures from Hm = n
compositional analysis. It applies to all common
types of utility gaseous fuels, i.e., dry natural gas,
∑x
j −1
j ×M j
reformed gas, oil gas (both high and low Btu), (Equation 3-3)
propane-air, coke oven gas, and other gaseous Where:
fractions for which suitable methods of analysis
are designated.
xj= the mole fraction of Component j;
Mj = the molar mass of Component j;
The ideal gas heating value and ideal gas
relative density are calculated from the molar n = the total number of components; and
composition and the respective ideal gas Hmj = the ideal gross heating value per unit mass
values for the components; these values for Component j (at 60°F or 15.6°C), as
are then adjusted by means of a calculated tabulated in ASTM D3588 (ASTM, 2003).
compressibility factor. Values of Hm are independent of pressure,
but they vary with temperature.
38
ADDRESSING UNCERTAINTY IN OIL & NATURAL GAS INDUSTRY GREENHOUSE GAS INVENTORIES
Errors to be considered when computing heating values; however, these errors are
heating values include: errors in the heating negligible when compared to the errors
values of the mixture components and errors associated with the determination of the
in determining the mixture composition. The composition of individual species. Appendix D
uncertainty ranges of the heating values for contains a listing of common energy and fossil
the components cited in the ASTM practice fuel unit conversion factors that could be used
are about 0.03%. These errors may affect the for these and other calculations discussed in
agreement between calculated and measured these document.
39
IPIECA • API
3.6 Venting and fugitive emissions of instruments, drip pots, and scrubbers.
measurements Continuous applications may include
disposal of associated gas, treater off-gas
Vented emissions are releases to the and tank vapors at oil production facilities,
atmosphere that are due to process or casing gas at heavy oil wells, process waste,
equipment design, or operational practices. or byproduct streams. Depending on local
Vented emissions may come from a variety requirements, some of the vented streams
of emission points or non-fired stacks, or may be either captured for on-site use/sale
sometimes from flares as is commonly the or flared, depending on the economic value
approach used for categorizing emissions in of the stream and the availability of local
GHG emission inventory compilations that infrastructure.
adhere to the IPCC 2006 national inventory
guidance. In all cases, these are part of the These emission sources tend to be very specific
process design and are intended for pressure to the type of operation and detailed methods
relief purposes. Fugitive emissions are due to for quantifying these emissions are provided
equipment leaks from piping components and in Section 5.0 and 6.0 of the API Compendium
a variety of seals found in typical petroleum (3rd edition, August 2009).
and natural gas installations. Table 3-6 provides a summary of the vented
Vents may either be intermittent or continuous. and fugitive emission source types and
Intermittent applications may include the their measurement requirements under the
disposal of waste volumes from emergency EPA’s GHG Reporting Program (GHGRP). The
pressure relief episodes, operator initiated or emissions sources are categorized around the
instrumented depressurization events, plant subparts of the EPA regulation: Petroleum
or system upsets, well servicing and testing, and Natural Gas Systems (Subpart W) and
pigging events, and routine blowdown Petroleum Refining (Subpart Y).
Table 3-6. GHGRP measurement methods for vented and fugitive emission sources
Emission
Source Type GHGRP Measurement Method
Category
Petroleum and Natural Gas Systems
Vented -- Pneumatic device and pump vents, Engineering Estimate using
Emissions -- Gas well liquids unloading monitored process operating
-- Gas well completions parameters and either simulation
-- Gas well workovers models or emission factors
-- Storage tanks
-- Dehydrators
-- Acid gas removal
-- Compressor wet seal vents Direct Measurement where
-- Compressor rod packing vents meters depending on source
-- Transmission condensate tanks measured
Fugitive -- Compressor station leaks Leak Detection Surveys in
Emissions -- Processing plant equipment leaks conjunction with emission
-- Fugitives and non-compressor plant/station fugitives factors for identified “leaking”
components
-- Valves, connectors and pump seals from: Equipment Count in conjunction
• Offshore petroleum and natural gas production with applicable population
• Onshore petroleum and natural gas production emission factors
• Local distribution company mains,
• Service lines, vaults, and other inaccessible sources
Vented -- Fluid Catalytic Cracking Engineering Estimate using
Emissions -- Catalyst regeneration monitored process operating
-- Hydrogen production parameters and either simulation
-- Acid gas removal models or emission factors
-- Sulphur recovery
-- Cokers
40
ADDRESSING UNCERTAINTY IN OIL & NATURAL GAS INDUSTRY GREENHOUSE GAS INVENTORIES
41
IPIECA • API
42
ADDRESSING UNCERTAINTY IN OIL & NATURAL GAS INDUSTRY GREENHOUSE GAS INVENTORIES
43
IPIECA • API
4.1 Global approaches applicable to For example, the EU ETS now operates in 30
petroleum and natural gas sector countries and covers CO2 emissions from some
11,000 installations. In Japan, in 2009, over
Global trends indicate that a growing number 11,000 enterprises reported their CO2 emissions
of companies assess and address the potential under the mandatory GHG Accounting and
threats and opportunities of climate change Reporting system, accounting for about half of
for their business. At the same time, there is the total emissions of Japan nationwide. In the
also an increasing demand from governments U.S. around 8,000 entities reported data in 2012
and other stakeholders for corporate climate under the GHG Reporting Program, covering an
change-related information, therefore, estimated 85-90% of total U.S. GHG emissions.
encouraging or mandating companies to
disclose GHG emissions data. This has given Table 4-1 provides an inter-country comparison
rise to a myriad of sector-specific guidance for of mandatory reporting programs in selected
voluntary reporting of GHG emissions, such countries with significant petroleum and
as for the petroleum and natural gas industry natural gas industry operations. As indicated in
(IPIECA, 2011), along with increased requests the table, these programs have varying policies
from investors for data disclosure (CDP, 2010). for implementing their quality assurance,
At the same time, a wide range of regulatory verification, and uncertainty assessment
schemes are being developed that mandate requirements but they all share the common
GHG emissions reporting by installations and goal of ensuring the reporting of high quality
business entities at the national and regional data to be used as the basis of policy decisions.
levels (OECD, 2012).
2 http://www.ghgprotocol.org/programs-and-registries
44
ADDRESSING UNCERTAINTY IN OIL & NATURAL GAS INDUSTRY GREENHOUSE GAS INVENTORIES
Program
Countries
Elements
Australia Canada EU France UK USA
Program National Environment European Union Bilan démissions Greenhouse Greenhouse
Name Greenhouse Canada GHG Emissions Trading de GES Gas Emissions Gas Reporting
and Energy Emissions System Reporting Program
Reporting Reporting (GHGRP)
(EU-ETS)
(NGER) Program
(GHGRP)
Inception 2009 2004 2005 2011 2009 (Mandatory 2009
Date April 2013)
Legal National Canadian EU Regulation “Grenelle II” Climate Change FY2008
Authority Greenhouse Environmental 601/2012 Act of 2008 Consolidated
n°2010-788,
and Energy Protection Act and pending Appropriations
(21-JUN-2012),
Reporting (CEPA, 1999) – 12 July 2010 regulations Act (H.R. 2764;
Monitoring and
(NGER) Act of Section 46 (Article 75) PL 110– 161),
Reporting of
2007
GHG emissions
pursuant to
Directives
2009/29/EC and
2003/87/EC
Geographical Australia Canada EU- 28?, Norway, France Global emissions U.S. entities
Coverage Lichtenstein and of UK companies
Switzerland
Gases 6 Kyoto GHGs 6 Kyoto GHGs Primarily CO2, 6 Kyoto GHGs 6 Kyoto GHGs 6 Kyoto GHGs +
Reported with provision to NF3
include CH4 and
N2O in the future
Reporting Corporate Facility level Installation - Companies over All London Stock Facilities
Boundaries Control 500 employees; Exchange- listed
Fossil fuels or
- Sub-national firms
industrial GHGs
government
Consistent with suppliers
over 50,000
GHG Protocol
inhabitants; Vehicles
- Public bodies and engines
over 250 manufacturers
employees
Reporting 25 KtCOe 50 KtCO2e 25 KtCO2e per No No 25 KtCOe per
threshold (facility) (facility) year (excluding year
emissions from
50 KtCOe
biomass);
(company)
Combustion
Production /
of fuels in
consumption
installations with
> 200TJ
a total rated
thermal input
exceeding 20 MW
45
IPIECA • API
Table 4-1. Inter-Country comparison of select mandatory GHG reporting programs (continued)
Program
Countries
Elements
Australia Canada EU France UK USA
consumption > 200TJ Combustion
of fuels in
installations with a
total rated thermal
input exceeding
20 MW
Scope of - Direct operational - Direct operational - Direct regular - Direct operational Direct Direct
emissions emissions emissions operational emissions operational operational
- Indirect emissions emissions - Indirect emissions emissions emissions
from energy - Abnormal events from energy
Indirect Potential
consumption including start-up consumption
emissions from emissions from
- Other indirect and shut-down
energy distribution of
emissions (optional) - Emergency
fossil fuels and
situations
industrial gases
- Exception of
emissions from
mobile machinery
for transportation
purposes.
Calculation - 2009 Technical - Gov’t of Canada - Commission - Bilan methods DEFRA Guidance - Detailed
methods Guidelines Technical Guide Regulation (EU) démissions de GES on how to calculation
- Australian EFs Based on IPCC No. 601/2012: - Inspired by ISO measure and methods for
provided Guidelines and Good - Subsection 2 – 14064-1, 2006 report GHG each of the
Practice Guidance Activity Data - Consistent with emissions 41 subparts
- Sector-specific - Subsection 3 – GHG protocol addressing
(September 2009)
manuals Calculation Factors - EFs from ADEME sectors of the
- National EFs and - Subsection 4 – database economy that
GWPS Specific calculation are required to
- Conversion factors factors report
and Electricity Industry
intensity consensus
standards for
measurements
Reporting - Annual report to - Annual report to - Annual verified - Every three years Annual report Annual
Ministry of Climate Environment Canada report with financial reporting
Change and Energy - Electronic Data - Information report through the
Efficiency Reporting (EDR) reported in EPA’s electronic
- Online System system on the GHG rounded-off tonnes GHG reporting
Comprehensive reporting website of CO2e tool (eGGRT)
Activity Reporting
(OSCAR)
Quality - Internal - Information should - Member States - Internal Not specified - QA/QC checks
Assurance be “verifiable” are required specified in
to ensure that each of the rule
emissions are subparts
monitored in - Guidance for
accordance with substitution of
the monitoring missing data
and reporting
guidelines.
46
ADDRESSING UNCERTAINTY IN OIL & NATURAL GAS INDUSTRY GREENHOUSE GAS INVENTORIES
Table 4-1. Inter-Country comparison of select mandatory GHG reporting programs (continued)
Program
Countries
Elements
Australia Canada EU France UK USA
Verification - 3rd party - No specific - Member States are required - Not required - No specific - Self-
verification not requirements to ensure that the reports by law requirement certification
required for a 3rd party submitted by operators are for 3rd party by facilities’
verification. verified verification of registered
GHG data designated
representatives
- Verification
checks by EPA
Uncertainty - Required - No specified - Maximum Permissible - Not specified - Not - Specifications
assessment - Uncertainty ranges assessment Uncertainty for each tier of beyond GHG addressed for calibrations
provided with EFs required activity data (Anx II, Table1) protocol in General
guidance Provisions
(Subpart A)
Since not all governments have reporting by major investors. Table 4-2 provides a high
programs, a few non-governmental initiatives, level comparison of the scope and coverage
require or encourage enterprises to quantify of mandatory reporting in the USA. under the
and report their GHG emissions. These voluntary GHGRP as compared to the voluntary disclosure
initiatives are part of environmental and other of GHG emissions under the CDP, which
non-financial disclosure requirements, with brings into the voluntary CDP reporting the
the most notable one being that of the Carbon measurements and QA/QC checks associated
Disclosure Project, which is backed globally with mandatory GHG reporting in the USA.
Table 4-2. Comparison of reporting scope between EPA’s GHGRP and CDP
EPA Greenhouse Gas
CDP Comment
Reporting Program
Reporting - Monitor and report GHG emissions - Monitor and report GHG emissions - EPA’s GHGRP supplements and
requirements - Include climate change related complements other US government
risks, opportunities, strategies & programs that may include CDP subject
performance) matter
Reporting - Facilities over 25,000 tonnes CO2e - At the corporate level (including - Information about GHG emissions from
responsibilities the facilities they own, manage or facilities is a subset of the information
- Suppliers report at the corporate
operate) requested by CDP
level.
Scope of - 41 industry sectors - Same as EPA GHGRP, including - CDP includes all of the activities covered
reporting combustion and fugitive emissions by the EPA GHG RP and a few more
- Examples: stationary fuel
combustion, electricity generation, - Mobile source emissions
industrial and chemical processing,
- Petroleum and Natural Gas
Suppliers of fuels, products and
industrial gases.
Reporting - Annually (calendar year) - Annually (latest or more recent - For CDP reporting could coincide with
frequency 12 month period for which data is financial reporting
reported)
47
IPIECA • API
Table 4-3. EU regulation maximum permissible uncertainty for activity data tiersa
Parameter to which uncertainty
Tier 1 Tier 2 Tier 3 Tier 4
is applied
Combustion of fuels and fuels used as process input
Commercial standard Amount of fuel [metric tons] or
7.5% 5.0% 2.5% 1.5%
fuels [Nm3]
Other gaseous and liquid Amount of fuel [metric tons] or
7.5% 5.0% 2.5% 1.5%
fuels [Nm3]
Solid fuels Amount of fuel [metric tons] 7.5% 5.0% 2.5% 1.5%
Flaring Amount of flare gas [Nm3] 17.5% 12.5% 7.5%
Petroleum Refining
Hydrocarbon feed [metric tons]
(mineral oil refining)
Uncertainty requirements apply
Catalytic cracker, catalyst separately for each emission source
regeneration and flexi- 10.0% 7.5% 5.0% 2.5%
cokers (Total uncertainty of all emissions
from a given source)
Hydrogen production 7.5% 2.5%
a
Excerpt from Table 1 of Annex II of the June 21, 2012 EU-ETS regulation.
3 For the EU-ETS: Category A installation, where average verified annual CO emissions are equal to or less than 50,000 tonnes
2
of CO2e; Category B installation, where the average verified annual CO2 emissions are more than 50,000 tonnes of CO2e and
equal to or less than 500,000 tonnes of CO2e; Category C installation, where the average verified annual CO2 emissions are
more than 500,000 tonnes of CO2e.
48
ADDRESSING UNCERTAINTY IN OIL & NATURAL GAS INDUSTRY GREENHOUSE GAS INVENTORIES
c. EPA GHGRP. The regulations governing date at which data collection is required
the requirements for the GHG Reporting to begin using a measurement device. All
Program (GHGRP) in the USA. are included measurement devices must be calibrated
in many subparts to Part 98 in Title according to a method specified either in a
40 of the Code of Federal Regulations relevant subpart of the regulations in Part
(40CFR98). The rules are specified under 98 or the manufacturer’s recommended
general provisions along with sector- procedures, or an appropriate industry
specific subparts that invoke the general consensus standard.
provisions and add more requirements that
Highlights of the calibration and
are pertinent to specific industry sectors.
maintenance requirements are provided
General Provision §98.3 specifies that
in Table 4-4 below for all equipment and
reporters must follow the procedures for
measurements used in association with
emission calculation, monitoring, quality
applicable subparts to provide data for the
assurance, missing data, recordkeeping,
GHG emission calculations. Subsequent
and reporting that are specified in each
recalibrations of the flow meters and
relevant subpart.
other measurement devices must also be
In accordance with these provisions, flow performed, as shown in Table 4-4. If no
meters and other devices that measure recalibration frequency is specified in the
data used to calculate GHG emissions applicable subpart, the recalibrations must be
(e.g., GC analyzers used to determine performed at the frequency recommended
composition data and flow meters used by the manufacturer or as specified by an
for calculations) must be calibrated by the industry consensus standard practice.
49
IPIECA • API
4.3 Tools for calculating uncertainty variables required for the Monte Carlo
simulation software.
A few reporting programs have developed
tools to assist reporters in quantifying GHG Monte Carlo simulation has an advantage
emission uncertainty. Alternatively, general over uncertainty propagation in that one
statistical methods or commercially available can specify multivariate distributions
tools can also be used to evaluate uncertainty. to account for correlations between
A brief description of some common methods different sources of uncertainty. Ideally,
are provided below. the distributions should be derived from
data and knowledge of the underlying
process. It is helpful in many instances
4.3.1 Numerical methods for calculating to first graph the data, and then use the
uncertainty shape of the graph to determine the
underlying distributions. The IPCC Good
a. Monte Carlo Simulation
Practices document recommends choosing
Monte Carlo simulation is a complex, one of the following distributions: normal,
model-based method for iteratively lognormal, uniform, or triangular.4
evaluating uncertainty associated with
It is then necessary to statistically test the
individual parameters. It applies to any
hypothesis that the data follow a certain
probability distribution and may be the
distribution. The test will vary based on the
preferred approach where more complex
hypothesized distribution. For example, the
equations are assessed, or where more
Shapiro-Wilks test is often used to test if
than two correlated variables exist. It is one
the data are normal or lognormal. Options
of many methods for analyzing uncertainty
to test for other distributions include
propagation where the goal is to determine
Empirical Distribution Functions.
how random variation, lack of knowledge,
or error affects the sensitivity, performance, Selecting the proper probability
or reliability of the resulting emission distribution function for the model
inventory. parameters may be difficult in light of
inadequate data or a lack of understanding
The principle of the Monte Carlo
of the underlying physical processes. If
analysis is to generate random inputs
there are not enough data to assume
from probability distributions of the
normality by the Central Limit Theorem
respective variables and to calculate
(more than 30 data points), there are
the corresponding emission values. This
most likely not enough data to determine
procedure is repeated many times using
the underlying distribution of the data.
computer software, such as @RISK or
Consequently, such analyses are often
Crystal Ball. The uncertainty model relies
forced to rely on subject matter expert
on repeated random sampling of all
opinion to determine distributions rather
inputs and simultaneous recalculation
than on observed data.
of emissions (outputs) to measure
variation over the course of numerous Once distributions are determined for
model iterations. The data generated all of the data sources, the Monte Carlo
from the Monte Carlo simulation can be simulation will proceed by randomly
represented as probability distributions sampling each of the distributions that
(or histograms) or converted to error bars describe the data used for estimating
and confidence intervals. emissions. As many as 10,000 replicate
samples are typically taken, with the
Monte Carlo simulation can also be
total emissions being estimated for each
computationally intensive, with 10,000
replicate. These repeated determinations
simulations being the norm. A significant
of emissions are used to generate a
difficulty/barrier is the need to build a
distribution of the total emissions
spreadsheet-based simulation model
with its mean being the estimate of
for uncertainty analysis that reflects the
total emissions, and its uncertainty
complexity of the emissions estimation
determined by its variance.
model and allows generation of the
4 The IPCC Good Practices document discusses how to perform Monte Carlo simulation in Section 6.2 (IPCC, 2000). It discusses
choice of distribution in Section A1.2.5.
50
ADDRESSING UNCERTAINTY IN OIL & NATURAL GAS INDUSTRY GREENHOUSE GAS INVENTORIES
51
IPIECA • API
52
ADDRESSING UNCERTAINTY IN OIL & NATURAL GAS INDUSTRY GREENHOUSE GAS INVENTORIES
5. UNCERTAINTY CALCULATION
Section Focus
OVERVIEW
Uncertainty is used to characterize the
dispersion of values that could be reasonably This section provides an overview
attributed to a measured quantity. Uncertainty of general uncertainty calculation
may be expressed as a qualitative ranking, such principles and more detailed guidance
as the letter ratings assigned in EPA’s AP-42 on the calculation methods that are
publication series (USA the EPA, 1currently applicable for defining single source
has995 with Supplements through 2000), uncertainties and for aggregating
or as a quantified value. For the purpose of them at the facility (or entity) level.
quantifying the uncertainty of a GHG inventory, Uncertainty aggregation equations and
this section first addresses measurement an example are provided. Methods for
uncertainty, then discusses uncertainty reducing uncertainty are also addressed.
associated with emission factors, and finally
The subsections address:
addresses the propagation of uncertainty.
A step-by-step tutorial with examples of • Measurement uncertainty
the uncertainty calculations is provided in • Uncertainty propagation
Appendix F. The aggregation of uncertainty
estimates for two example facility-based
• Uncertainty aggregation
GHG emission inventories are provided in • Reducing uncertainty
Appendices G and H.
Figure 5-1. Measurement error over time of an unbiased estimate The general approach for quantifying bias
would depend on prior experience in the
laboratory or from specifically designed field
measurement campaigns. Measurement bias can
vary from small to very large, depending on the
application, and can even change over time if
the measurement instrument is allowed to drift
without calibration. Therefore, in practice, bias
will most commonly be determined using expert
judgment and will be based on such parameters
as the length of time since the equipment was
calibrated and other factors that would cause the
measurement to systematically overestimate or
underestimate the true value.
s
±t ×
n (Equation 5-1)
Where:
s = standard deviation;
n = sample size; and
t = t-value for “n-1” degrees of
freedom
And
uncertainty due to the bias generally does not
decrease as more data points are collected. 1 n
Thus, precision and bias should be considered s= ∑
n − 1 i=1
( xi − x ) 2
(Equation 5-2)
separately, if possible.
Where:
If bias can be quantified, it should be corrected
and thus eliminated from consideration in
quantifying uncertainty. For example, if a
xi = the ith observation in the data
set and
fuel stream has two types of measurement
devices, data can be collected from both x = the mean of the data set.
devices to check for agreement. A bias would
n
be indicated if the measurements differed.
However, quantifying and correcting for ∑x i
54
ADDRESSING UNCERTAINTY IN OIL & NATURAL GAS INDUSTRY GREENHOUSE GAS INVENTORIES
5.1.3 Implementation of Pedigree Matrix Table 5-1. Example of Pedigree Matrix for
approach determining uncertainty scaling factors
based on data quality ratingsa
The Pedigree Matrix could be a useful tool when
measured single parameter uncertainties are Very
unknown. This approach may be used to elicit Indicator score Good Fair Poor
good
a structured expert’s judgment for determining
Precision 1.00 1.10 1.20 1.50
single parameter uncertainty values, which can
then be propagated using techniques such as Completeness 1.00 1.05 1.10 1.20
series expansion or Monte Carlo simulations, as Temporal
discussed below. 1.00 1.10 1.20 1.50
Representativeness
In the pedigree matrix approach, qualitative data Geographical
1.00 1.02 1.05 1.10
quality assessment is used to select appropriate Representativeness
uncertainty scaling factors for single parameters. Technological
Assessment of data quality should be conducted 1.00 1.20 1.50 2.00
Representativeness
separately for activity data and emission factors, a
Weidema and Wesnaes, 1996
which will be merged in the ensuing uncertainty
propagation analysis.
In this proposed adaptation, which is based on and represent the contribution to the square of
guidance provided by the GHG Protocol (WRI/ the geometric standard deviation.
WBCSD a, 2011), an uncertainty scaling factor is
Data sources used when compiling the
selected by assigning a data quality score (very
emissions inventories need to be assessed. The
good, good, fair and poor) to each of five data
total uncertainty, expressed as a 95% confidence
quality indicators (precision, completeness,
interval, SDg95 (the square of the geometric
temporal geographical and technological
standard deviation), is calculated using the
representativeness). The ranges of recommended
formula shown below (WRI/WBCSD a, 2011):
uncertainty scaling factors are shown in Table 5-1
55
IPIECA • API
Table 5.2. Basic Pedigree Matrix uncertainty factors for selected pollutants a
Indicator score Very good Good
CO2 1.05 1.05
CH4 1.50 Not available
N0x N20 1.50 1.40
NMVOC total 1.50 1.50
NH3 1.50 1.20
Individual Hydrocarbons 1.50 2.00
a
Overview and Methodology, Data v2.0 (2007). Ecoinvent Report No. 1, Dübendorf, December 2007.
http://www.ecoinvent.org/fileadmin/documents/en/01_OverviewAndMethodology.pdf
56
ADDRESSING UNCERTAINTY IN OIL & NATURAL GAS INDUSTRY GREENHOUSE GAS INVENTORIES
there were no common source of uncertainty. is calculated as the “square root of the
The uncertainties in two quantities would be sum of the squares” using the absolute
dependent if they had a common source of uncertainties, as shown in Equation 5-5.
uncertainty (Williamson et al, 1996).
Covariance between two uncertainty terms U ( abs ) X +Y +...=
N
U X2 + U Y2 + ... + U N2
can be addressed through an additional term
in the uncertainty propagation equations (Equation 5-5)
(discussed further in Section 4.2.2). However, Where:
the IPCC Good Practices document suggests
avoiding the need for the covariance term U(abs) = the absolute uncertainty.
in the equation by “…stratifying the data or
The absolute uncertainty values are
combining the categories where the covariance
used in the equations, and the resulting
occurs” (IPCC, 2001).
aggregated uncertainty (UX+Y+…+N) is also
on an absolute basis. Note that where a
constant is also included in the emission
5.2.1 Propagation equations
estimation calculation, the absolute
Four general equations for propagating uncertainty should include the constant.
uncertainty are used in this document and the This is demonstrated in the example
API Compendium for compiling the uncertainty provided in Appendix F.
associated with a GHG inventory. A general
For two uncertainty parameters that are
introduction to the equations is presented
related to each other, the equation becomes:
here, with example applications provided in
Appendix F.
Consider two quantities that can be measured: U ( abs )CorrelatedX +Y = U X2 + U Y2 + 2r (U X × U Y )
X and Y. The uncertainty for these values can
be expressed on an absolute basis as ±Ux and (Equation 5-6)
±Uy respectively, where U is calculated through Where:
statistical analysis (as represented by Equation
5-1), determined through the Monte Carlo r = the correlation coefficient between UX,
technique, or assigned by expert judgment. UY, (discussed further in Sections 5.2.3 and
Uncertainty may also be expressed on a relative Appendix F).
basis, generally reported as percentage: However, the IPCC Good Practices
guidance states, “Once the summation
U UY exceeds two terms and the covariance
±100 X % or ± 100 % respectively
X Y occurs, the use of the Monte Carlo
approach is preferable where data
Depending on the uncertainty propagation resources are available” (IPCC, 2000).
equation, the absolute or relative uncertainty
value may be required. In addition, selection of the
b. Uncertainty propagation for a product
propagation equation also depends on whether
(or quotient)
the uncertainties associated with the individual
parameters are independent or correlated. The equation for propagating uncertainties
from the product or quotient of two
a. Uncertainty propagation for a sum (or
or more measured and independent
difference)
quantities is similar to Equation 5-5.
Two potential equations are used for However, in this case the relative
computing the total uncertainty from the uncertainties are used, as shown in
addition or subtraction of two or more Equation 5-7. When multiplied by 100, the
measured quantities. The selection between resulting combined uncertainty (U(Rel)XxYxN)
the two equations depends on whether the is expressed as a percentage.
uncertainties associated with the measured
2 2 2
quantities, X and Y, are correlated. U X UY UN
U ( rel ) X ×= U (rel ) X ÷=
Y ÷...÷ N + + ... +
For uncertainties that are mutually
Y ×...× N
X Y N
independent, or uncorrelated (i.e., the
uncertainty terms are not related to (Equation 5-7)
each other), the aggregated uncertainty
57
IPIECA • API
Equation 5-7 is used to estimate the These guidelines concentrate solely on the
uncertainty of a product or quotient uncertainty propagation method due to the
of two parameters (X and Y) where the potential to introduce further errors in assigning
uncertainties are correlated and positive the probability distributions for the Monte Carlo
values. Here also, relative uncertainty simulations. As stated previously, applying the
values are used in the equation and the uncertainty propagation methods, even where
resulting combined uncertainty is on a the assumptions are not met, is advised in these
relative basis. guidelines, particularly for emission sources
with a small contribution to the overall GHG
2 2
U X UY U X UY inventory. As data collection methods improve
U ( rel )CorrelatedX ×Y = + + 2r × for GHG inventories, the ability to quantify
X Y X Y uncertainties will also improve.
(Equation 5-8)
5.2.3 Correlation coefficient
The correlation coefficient, r, used in Equations
c. Combining uncertainties
5-5 and 5-7, is a number between -1 and 1
It may be necessary to combine multiple that measures the linear relationship between
uncertainty parameters associated the errors or uncertainties of two measured
with a single measured value, such as parameters. The value of r is zero when the
combining uncertainties for precision and parameters are independent. As stated
bias. For uncertainty parameters that are previously, once the uncertainty propagation
independent, the combined uncertainty is exceeds two terms and covariance occurs, the
calculated using the absolute uncertainties use of the Monte Carlo approach (described
as shown in Equation 5-5. Similarly, for further in Section 4.0) is preferable (IPCC, 2001).
uncertainty parameters that are related to Additional details on calculating the correlation
each other, Equation 5-6 applies. coefficient are provided in Appendix F. A
simplified explanation follows.
5.2.2 Comparison of uncertainty For two terms that might be correlated, the
propagation and Monte Carlo errors or uncertainties are plotted against each
other. For the purpose of this discussion, Ux
Section 6.3.1 of the IPCC Good Practices represents the uncertainties of one variable
document compares the uncertainty plotted along the x-axis, and Uy represents the
propagation method and the Monte Carlo uncertainties of the second variable plotted
simulations (IPCC, 2000). It notes that on the y-axis. The correlation coefficient, r, is
the uncertainty propagation method’s determined by a linear regression of the Ux and
assumption of normality leads to symmetric Uy values.
95% confidence intervals whereas the Monte
Carlo method can take into account the fact If one suspects that the uncertainty parameters
that emissions are bounded below by zero are correlated, but data are not available to
to fit an asymmetric (and thus narrower) plot or calculate the correlation coefficient, the
confidence interval. If the data are skewed and following rule-of-thumb values could be applied
one transforms the data (discussed earlier), using expert judgment (Franzblau, 1958):
one could achieve the asymmetric confidence • r = 0: no correlation, the data are
intervals using uncertainty propagation, as well. independent
Since the Monte Carlo simulations can • r = ±0.2: weak correlation
assume a truncated distribution, the lower • r = ±0.5: medium correlation
confidence limits tend to be closer to the • r = ±0.8: strong correlation
mean than the upper confidence limits. The
IPCC Good Practices document goes on to
• r = ±1: perfect correlation, the data fall on a
straight line.
state that the two methods produce results
that are fairly comparable. It recommends that
countries report the results of the uncertainty
propagation method and those countries with
“sufficient resources and expertise” report
Monte Carlo results as well.
58
ADDRESSING UNCERTAINTY IN OIL & NATURAL GAS INDUSTRY GREENHOUSE GAS INVENTORIES
5.3.1 IPCC aggregation template Table 5-3 exhibits a generic template that
is an adaptation of the IPCC-recommended
IPCC provides an example calculation table approach for calculating aggregated
for estimating the aggregated uncertainty of uncertainties. This template could be tailored
GHG emissions, which is based on the error for specific oil and natural gas inventories by
propagation method for the multiplication expanding or collapsing the emission sources
of an emission factor and activity value (IPCC, and industry segments included based on the
2006, page 3.31). The approach entails the goals of the assessment and data availability.
Offshore
Production
Natural
Gas
Processing
CH4
Onshore
Production
Offshore
Production
Natural
Gas
Processing
2
ln
ΣH ∑
n
i =1
× Hn
100
TOTAL l= ×100%
∑ H
59
IPIECA • API
Table 5-4 provides an example using such a emissions from combusting natural gas.
table structure for calculating uncertainty, Equations 5-6 and 5-4 are applied, assuming
similar to the IPCC table and WRI tool. This that the variables are uncorrelated. More
example illustrates the aggregation of detailed examples are provided in Appendices
uncertainty for the combined CO2 and CH4 F and G.
A B C D E F G H I
Emission
Source Activity Emission Emis- Emissions Uncertainty of Calculated
Activity
Description: data Factor sions in in tonnes Emissions (expressed in
Data (e.g. Activity GHG Emission
Example for Quantity Data Uncertainty emission Factor Uncertainty tonnes CO2e ± Percent)
Natural Gas of fuel Units factor Units
Combustion used) ±% ±% A×D G × GWP
1 = C2 + F 2
tonnes
MMBtu/
CO2 Emission 3000 +/- 5.0% 0.0732 CO2/ +/- 10.0% 219.6 219.6 11.2%
yr
MMBtu
tonnes
MMBtu/
CO4 Emission 3000 +/- 5.0% 3.01x10-6 CO4/ +/- 15.0% 0.009 0.19 15.8%
yr
MMBtu
2
ln
∑ i=1 100
n
Σ
× Hn
H l=
×100%
TOTAL ∑ H
219.79 24.59
= ×100% = 11.2%
219.79
60
ADDRESSING UNCERTAINTY IN OIL & NATURAL GAS INDUSTRY GREENHOUSE GAS INVENTORIES
5.3.2 Pedigree Matrix example fuel combusted is metered, and the meter
is calibrated annually. The natural gas is
Table 5-5 provides an example application of
purchased from a local distribution company.
the Pedigree Matrix approach for determining
the uncertainty associated with CH4 emissions The uncertainty, expressed as a 95% confidence
from the combustion of natural gas in a interval, is calculated by applying Equation 5-4
compressor. For this example, the amount of and using the uncertainty factors from Table 5-5.
( lnU 1)2 + ( lnU 2)2 + ( lnU 3)2 + ( lnU 4)2 + ( lnU 5)2 + (ln Ub )2
SDg95 ≅ σ 2 =
exp
( ln1.1)2 + ( ln1.0)2 + ( ln1.0)2 + ( ln1.05)2 + ( ln1.0)2 + ( ln1.5)2
= exp
SDg95 = 1.52
The 95% probability range for this example 5.3.3 Uncertainty assessment from API
would be expressed as (median/SDg95) to Technical Report 2571
(median* SDg95). Hence, for this example, if the
Table 5-6 provides an illustration for performing
median CH4 emissions are 5,000 tonnes per
the uncertainty assessment in conjunction with the
year the 95% probability range wo uld
equations and processes set forth in the TR 2571
be 3,290 to 7,600 tonnes per year.
document, using average values as an example.
Table 5-6. Example of flow measurement total uncertaintya
Standard Sensitivity Total
Nominal Standard
Component Units Uncertainty Coefficient Uncertainty
Value Uncertainty
(Ux)% (Sx) (Ux×Sx)2
Unit Conversion Factor NA 1.899 0.0 0.0000 1.00 0.000
Discharge Coefficient Dimensionless 0.6008 0.00411 0.6834 1.00 0.467
Expansion Factor Dimensionless 0.9989 0.00014 0.0144 1.00 0.000
Bore Diameter Inches 0.5 0.0004 0.0750 2.04 0.023
Pipe Diameter Inches 2.067 0.0257 1.2425 0.04 0.002
Base Compressibility Dimensionless 0.9979 0.0010 0.1000 1.00 0.010
Base Temperature °Fahrenheit 60 0 0.0000 1.00 0.000
°Rankine 519.67
Flowing Pressure psia 100 0.46 0.4575 0.5 0.052
Universal Gas Constant ft-lbf/lbmol-°R 1545.33 0 0.0000 0.5 0.000
Gravitational Constant lbm-ft/lbf-s2 32.174 0 0.0000 1.00 0.000
Differential Pressure Inches-H2O 10 0.1120 1.1204 0.50 0.314
Base Pressure psia 14.73 0 0.0000 1.00 0.000
Molecular Weight lbm/lb-mol 16.799 0.1260 0.7500 0.50 0.141
Flowing Dimensionless 0.98664 0.007 0.7500 0.50 0.141
Compressibility
Flowing Temperature °Fahrenheit 60 2.598 0.5000 0.50 0.063
°Rankine 519.67
Installation Effect Dimensionless 1 0.02 2.0000 1.00 4.000
Factor
Meter Condition Factor Dimensionless 1 0.02 2.0000 1.00 4.000
Energy Content Btu/ft3 1050 7.88 0.7500 1.00 0.563
Carbon Content b
lbm/ft 3
0.035 0.0002625 0.7500 1.00 0.563
SUMMARY
Volume Flow Rate SCFH 848 Total Volume Uncertainty 3.04
Energy Flow Rate Btu per hour 890,320 Total Energy Uncertainty 3.13
Carbon Flow Rate Tons per hour 0.015 Total Carbon Uncertainty 3.13
a
Extracted from API TR 2571, Table A.3.5
b
Carbon Content is the fraction of carbon in the fluid expressed as percent by weight.
61
IPIECA • API
The example in Table 5-6 assumes a 2 in. line certain percentage level. A company may also
with a 0.5 in. orifice plate using a multivariable independently wish to refine its own uncertainty
transmitter for measuring temperature, limits, if it considers the uncertainties too large.
pressure, and pressure differential (DP). It also
The reader should note that none of the
assumes that the gas composition is measured
strategies mentioned here are aimed at reducing
at the fuel gas drum by an online analyzer;
the actual emissions of GHGs. That is a separate
hence lower molecular weight uncertainty
subject, and while emission reductions are
values are used.
achievable in some cases, they are beyond
The total uncertainty was derived by summing the scope of this report. This section focuses
the individual combined uncertainties (Ux × on reducing the mathematical and statistical
Sx)2 and taking the square root of the sum. uncertainty associated with an existing emission
The calculation for this example is done inventory.
three different ways to derive three different
Once a decision has been made to refine and
uncertainties (volume, energy, and carbon).
reduce the uncertainty associated with a given
• Volume Uncertainty: all components are emissions inventory, some strategic analysis
included with the exception of energy of the major sources contributing to the
content and carbon content; uncertainty is in order. It is important to know
• Energy Uncertainty: all components are which emission sources and which uncertainties
included with the exception of the carbon significantly contribute to the overall inventory
content; uncertainty. It would make little sense to refine
a term with large confidence bounds, but that
• Carbon Uncertainty: all components are
contributed very little to the overall inventory
included with the exception of the energy
uncertainty. It may also be useful to have a target
content.
uncertainty in mind. For example, if the current
total uncertainty is ±50.0 % and the company
5.4 Strategic reduction of wishes to reduce it to ±25.0%, then that target
uncertainty is useful in the analysis. As uncertainties
of individual values in the calculations are
examined, values that have uncertainties that
5.4.1 Periodic updates are already at or below the 25.0% target, are less
likely to be fruitful targets for reduction. The goal
Activities data that comprise a GHG emission is to find the largest contributors to emissions
inventory will likely change from year to year. that have the largest uncertainties.
Emission factor information may also vary
over time. However, aspects that impact the The largest contributors to uncertainty can
calculation of uncertainty such as measurement be determined by multiplying the emission
equipment or techniques, quantification estimate by the maximum error bound for each
methods, and emission factors, are relatively source. This would result in the upper bound
constant year-on-year. Therefore, it may emission estimate for the particular source. The
not be necessary to conduct an uncertainty emission estimates can then be sorted by the
assessment on an annual basis. A frequency of largest contributors. This is demonstrated in
every three to five years may be sufficient, but Table 5-7 on page 59 for the example crude oil
should ultimately be determined based on the production facility presented in Appendix G.
variability of the factors that contribute to the a. Onshore oil field
quantification of uncertainty.
Each facility should examine the major
categories of emissions and emission
5.4.2 Reducing uncertainty of emission uncertainty, and then examine specific
estimates emission sources within the category.
Prioritizing the largest sources of uncertainty
This section addresses the potential need to can be done with this simple approach.
refine the emission inventory to reduce the
uncertainty in the overall emission estimate. Figure 5-3 illustrates emissions for the
There may be several reasons to do this. In facility with each emission category
some specific locations, there may be state or and each gas converted to CO2e, using
national regulations, or guidelines for voluntary each gas’ GWP. It should be noted that
programs that suggest or require an emission although the GWP itself has uncertainty
inventory to have uncertainties lower than a associated with it, this report treats the
62
ADDRESSING UNCERTAINTY IN OIL & NATURAL GAS INDUSTRY GREENHOUSE GAS INVENTORIES
Table 5-7. Emission uncertainty ranking for onshore oil production example
Maximum
Emissions Maximum
Source Type Source Uncertainty, Ranking
(tonnes CO2e /yr) Uncertainty, %
(tonnes CO2e /yr)
Combustion Boiler/heaters 5,210 8.77 457 9
Sources
Natural gas engines 14,100 15.6 2,200 6
Diesel engines 220 15.5 34
Flares 30,700 21.1 6,478 3
Fleet vehicles 129 19.2 25
Vented Sources Dehydration and Kimray 5,440 76.0 4,134 4
pump vents
Tanks-flashing losses 40,300 88.7 35,746 1
Amine unit 66,700 9.77 6,517 2
Pneumatic devices 3,360 49.2 1,653 7
Chemical injection pumps 2,530 106 2,682 5
Vessel blowdowns 3.65 319 12
Compressor starts 38.7 187 72 10
Compressor blowdowns 17.3 175 30
Well workovers 0.939 294 3
Other non-routine (PRVs) 6.81 319 22
Fugitive Sources Fugitive components 1,100 83.3 916 8
Fleet vehicle refrigeration, 1.30 112 1.5
R-314a
Indirect Emissions Electricity consumed 553 10.2 56
GWP as a selected constant. Therefore, no Figure 5-3. Onshore oil field: summary of co2
uncertainty is associated or propagated equivalent emissions and uncertainties
from the GWP values. Figure 5-3 also
illustrates the bounds of uncertainty (at
95% confidence) for each emission source
category.
Among the major types of emissions
in Figure 5-3, are three very significant
emission categories: vented sources
of CO2; vented sources of CH4; and
combustion sources of CO2. These are the
most significant source of greenhouse
gas emissions for this facility. Together
they comprise almost 95% of all GHG
emissions from the facility. Examination
of the uncertainties, as shown in the bars
on Figure 5-3, or the absolute uncertainty
values in Table 5-7, shows that vented
emissions of CH4 from tank flashing
losses are the most significant source of Therefore, should the company wish to
uncertainty, contributing 35,746 CO2e reduce uncertainty in the GHG emission
tonnes of uncertainty (based on 40,300 inventory from this onshore oil field
±88.7% CO2e tonnes). Carbon dioxide facility, these categories would be the
vented from the amine unit is the next primary targets for uncertainty reduction.
largest source of uncertainty, contributing These appear to be the sources where
6,517 CO2e tonnes of uncertainty. The third uncertainty reduction efforts could be
largest source are emissions from flares, effectively undertaken to refine the
with 6,478 CO2e tonnes of uncertainty. inventory and reduce the uncertainty.
63
IPIECA • API
5
This emission estimate is based on an assumed separator pressure of 30 ±5% psi, and an assumed uncertainty of ±50% applied to
the correlation constant used in API Compendium, Equation 5-20.
64
ADDRESSING UNCERTAINTY IN OIL & NATURAL GAS INDUSTRY GREENHOUSE GAS INVENTORIES
65
IPIECA • API
6. REFERENCES
AER (2013). Directive 017: Measurement Requirements for Upstream Oil and Gas Operations. Alberta
Energy Regulator. Formerly released by the Energy Resource Conservation Board. Last released 15
May 2013. www.aer.ca/documents/directives/Directive017.pdf (Accessed 12 September 2014.)
AGA (2008). Greenhouse Gas Emission Estimation Guidelines for the Natural Gas Distribution Sector.
American Gas Association, Washington D.C., April 2008.
API (1985). Manual of Petroleum Measurement Standards. Chapter 13: ‘Statistical Aspects of Measuring
and Sampling’, Part 1: ‘Statistical Concepts and Procedures in Measurement’. First Edition, June 1985;
reaffirmed February 2011. American Petroleum Institute. API Publications, Washington D.C.
API (1999). Amine Unit Air Emissions Model, AMINECalc Version 1.0. API Publication 4679. Software by
American Petroleum Institute, 31 January 1999 (errata June 1999). http://www.techstreet.com/api/
searches/1088133 (Accessed 12 September 2014.)
API (2006a). Manual of Petroleum Measurement Systems. Chapter 14: ‘Natural Gas Fluids
Measurement’, Sections 1–10. 6th Edition, February 2006; reaffirmed September 2011. American
Petroleum Institute. API Publications, Washington D.C.
API (2006b). Manual of Petroleum Measurement Systems. Chapter 14.1: ‘Collecting and Handling of
Natural Gas Samples for Custody Transfer’. 6th Edition, February 2006; reaffirmed September 2011.
American Petroleum Institute. API Publications, Washington D.C.
API (2007). Manual of Petroleum Measurement Systems. Chapter 14.10: ‘Measurement of Flow to
Flares’. 1st Edition, July 2007; reaffirmed June 2012. American Petroleum Institute. API Publications,
Washington D.C.
API (2009). Compendium of Greenhouse Gas Emissions Estimation Methodologies for the Oil and Natural
Gas Industry. Version 3.0, August 2009. American Petroleum Institute, Washington D.C.
API (2011). Fuel Gas Measurement. API Technical Report 2571, 1st Edition, March 2011. American
Petroleum Institute, Washington D.C.
API (2013). Carbon Content, Sampling, and Calculation. API Technical Report 2572, 1st Edition, May
2013. American Petroleum Institute, Washington D.C.
API (2013). Manual of Petroleum Measurement Systems. Chapter 14.3: ‘Concentric, Square-Edged
Orifice Meters.’ 4th Edition, November 2013. American Petroleum Institute. API Publications,
Washington D.C.
API/IPIECA (2007a). Oil and Natural Gas Industry Guidelines for Greenhouse Gas Reduction Projects.
March 2007. IPIECA, London.
API/IPIECA (2007b). Greenhouse Gas Emissions Estimation and Inventories: Addressing Uncertainty and
Accuracy. Report of the IPIECA/API Workshop held in Brussels, Belgium on 16 January 2007. IPIECA,
London, January 2007. http://www.ipieca.org/publication/greenhouse-gas-emissions-estimation-
and-inventories-addressing-uncertainty-and-accuracy (Accessed 12 September 2014.)
API/IPIECA/OGP (2011). Petroleum industry guidelines for reporting greenhouse gas emissions. 2nd
Edition. American Petroleum Institute/IPIECA/International Association of Oil and Gas Producers.
IPIECA, London, November 2011.
ASTM D1826-94 (2010). Standard Test Method for Calorific (Heating) Value of Gases in Natural Gas Range
by Continuous Recording Calorimeter. Originally approved in 1961. Last previous edition approved in
2003 as D1826-94 (2003). ASTM International, West Conshohocken, PA, 2010, www.astm.org
66
ADDRESSING UNCERTAINTY IN OIL & NATURAL GAS INDUSTRY GREENHOUSE GAS INVENTORIES
ASTM D1945-03 (2010). Standard Test Method for Analysis of Natural Gas by Gas Chromatography.
Current edition approved 10 May 2003. Reapproved 2010. Originally approved in 1962. Last previous
edition approved in 2001 as D1945–96 (2001). ASTM International, West Conshohocken, PA, 2010,
www.astm.org
ASTM D1946-90 (2011). Standard Practice for Analysis of Reformed Gas by Gas Chromatography.
Current edition approved 1 June 2006. Reapproved 2011. Originally approved in 1962. ASTM
International, West Conshohocken, PA, 2011, www.astm.org
ASTM D2650-10 (2010). Standard Test Method for Chemical Composition of Gases by Mass Spectrometry.
Previous edition published 1 November 2004 as D2650-04. ASTM International, West Conshohocken,
PA, 2010, www.astm.org
ASTM D3588-98 (2011). Standard Practice for Calculating Heat Value, Compressibility Factor, and
Relative Density of Gaseous Fuels. Originally approved in 1998. Previous edition approved in 2003 as
D3588-98 (2003). ASTM International, West Conshohocken, PA, 2011, www.astm.org
ASTM D4809-13 (2013). Standard Test Method for Heat of Combustion of Liquid Hydrocarbon Fuels
by Bomb Calorimeter (Precision Method). Previous edition approved in 2009 as D4809-09ae1. ASTM
International, West Conshohocken, PA, 2013, www.astm.org
ASTM D4891-13 (2013). Standard Test Method for Heating Value of Gases in Natural Gas and Flare Gases
Range by Stoichiometric Combustion. Previous edition approved in 2006 as D4891–89 (2006). ASTM
International, West Conshohocken, PA, 2013, www.astm.org
ASTM D5720-95 (2009). Standard Practice for Static Calibration of Electronic Transducer-based Pressure
Measurement Systems for Geotechnical Purposes. Previous edition approved in 2002 as D5720-95
(2002). ASTM International, West Conshohocken, PA, 2009, www.astm.org
ASTM D5291-10 (2010). Standard Test Methods for Instrumental Determination of Carbon, Hydrogen,
and Nitrogen in Petroleum Products and Lubricants. Previous edition approved in 2009 as D5291-09.
ASTM International, West Conshohocken, PA, 2010, www.astm.org
ASTM D5865-13 (2013). Standard Test Method for Gross Calorific Value of Coal and Coke. Last previous
edition approved in 2012 as D5865-12. ASTM International, West Conshohocken, PA, 2013, www.
astm.org
ASTM D7314-10 (2010). Standard Practice for Determination of the Heating Value of Gaseous Fuels using
Calorimetry and On-line/At-line Sampling. Previous edition approved in 2008 as D7314-08. ASTM
International, West Conshohocken, PA, 2010, www.astm.org
ASTM UOP539-12 (2012). Refinery Gas Analysis by Gas Chromatography. Previous edition UOP539-97.
ASTM International, West Conshohocken, PA, 2012, www.astm.org
CARB (2012). Regulation for the Mandatory Reporting of Greenhouse Gas Emissions. Final Regulation.
California Air Resources Board, Sacramento, California, 19 December 2012.
Casella, G., and Berger, R. L. (1990). Statistical Inference. Wadsworth & Brooks/Cole, Belmont, California.
67
IPIECA • API
CER (2012). Uncertainty Calculator 2013-14 (Microsoft® Excel™ document). Australian Government
Clean Energy Regulator, National Greenhouse and Energy Reporting (NGER) Scheme. September
2014. www.cleanenergyregulator.gov.au/National-Greenhouse-and-Energy-Reporting/Forms-and-
calculators/Pages/default.aspx (Last accessed 12 September 2014.)
Cochran, W. G. (1977). Sampling Techniques. 3rd Edition. John Wiley & Sons, Inc., New York.
Coleman, H. W. and Steele, W. G. (1989). Experimentation and Uncertainty Analysis for Engineers. John
Wiley & Sons, New York.
EC (2012). Guidance Document: The Monitoring and Reporting Regulation – Guidance on Uncertainty
Assessment. MRR (Monitoring and Reporting Regulation) Guidance document No. 4, Final Version of 5
October 2012. European Commission, Directorate-General Climate Action.
http://ec.europa.eu/clima/policies/ets/monitoring/docs/gd4_guidance_uncertainty_en.pdf
(Last accessed 12 September 2014.)
EEA (2013). EMEP/EEA air pollutant emission inventory guidebook: Technical guidance to prepare
national emission inventories. European Environment Agency Technical report No. 12/2013.
Copenhagen, Denmark. www.eea.europa.eu//publications/emep-eea-guidebook-2013 (Last
accessed 12 September 2014.)
EPA (1971). ‘Method 21 Guidelines’. In Standards of Performance for New Stationary Sources – Title 40
of the Code of Federal Regulations (40 CFR), Chapter 1, Subchapter C, Part 60, Appendix A-7: ‘Test
Methods 19 through 25E’. www.ecfr.gov/cgi-bin/text-idx?SID=1d40301728f289bfce4a3216ebdf2db4
&node=pt40.8.60&rgn=div5#ap40.8.60.a
EPA (1995). Technology Transfer Network, Clearinghouse for Inventories & Emissions Factors
(website): ‘Emissions Factors & AP-42, Compilation of Air Pollutant Emission Factors’. United States
Environmental Protection Agency, Research Triangle Park, North Carolina. www.epa.gov/ttn/chief/
ap42/ (Last accessed 12 September 2014.)
EPA (2008). ‘Alternative work practice for equipment leaks’. In Consolidated Federal Air Rule – Title 40 of
the Code of Federal Regulations (40 CFR), Chapter 1, Subchapter C, Part 65, Subpart A, Section 65.7:
‘Monitoring, recordkeeping, and reporting waivers and alternatives, and alternative work practice for
equipment leaks. www.ecfr.gov/cgi-bin/text-idx?SID=e80816c4cda3834bac4e17f73b08cba4&node=
se40.16.65_17&rgn=div8
EPA (2009). Mandatory Greenhouse Gas Reporting – Title 40 of the Code of Federal Regulations (40
CFR), Chapter 1, Subchapter C, Part 98. Published in the Federal Register (www.regulations.gov) on
30 October 2009 under Docket ID No. EPA-HQ-OAR-2008-0508-2278, and became effective on 29
December 2009; implementation of 40 CFR Part 98 is referred to as the Greenhouse Gas Reporting
Program (GHGRP). www.ecfr.gov/cgi-bin/text-idx?SID=dba613069e26d03a053c2b80fe0e07f9&node
=40:21.0.1.1.3&rgn=div5 (Last accessed 12 September 2014).
ETSG (2007). Compendium of notes ‘establishing guidelines for the monitoring and reporting of
greenhouse gas emissions pursuant to Directive 2003/87/EC of the European Parliament and of the Council.’
Document ref. 2007/589/EC. European Union Emissions Trading System (EU-ETS), Emissions Trading
Technical Support Group (ETSG), Brussels, Belgium, 18 July 2007. http://www.emissions-euets.com/
attachments/232_M&R%20Guidelines%20589_2007.pdf (Last accessed 12 September 2014.)
Franzblau, A. (1958). A Primer of Statistics for Non-Statisticians. Chapter 7, Harcourt, Brace & World,
New York.
GTI (2000). GRI-GLYCalc™ Version 4.0, Software, Ref. GRI-00/0102. Gas Technology Institute, Illinois.
http://sales.gastechnology.org/000102.html (Last accessed 12 September 2014.)
68
ADDRESSING UNCERTAINTY IN OIL & NATURAL GAS INDUSTRY GREENHOUSE GAS INVENTORIES
INGAA (2005). Greenhouse Gas Emission Estimation Guidelines for Natural Gas Transmission and Storage.
Revision 2.0, 28 September 2005. Interstate Natural Gas Association of America, Washington, D.C.
http://www.ingaa.org/File.aspx?id=5485
IPCC (EFDB). Emissions Factors Database (website), National Greenhouse Gas Inventories Programme
(NGGIP) of the Intergovermental Panel on Climate Change.
www.ipcc-nggip.iges.or.jp/EFDB/main.php (Last accessed 12 September 2014.)
IPCC (1996). 1996 IPCC Guidelines for National Greenhouse Gas Inventories. Intergovernmental Panel on
Climate Change, Geneva. http://www.ipcc-nggip.iges.or.jp/public/gl/invs1.html
IPCC (2000). Good Practices Guidance and Uncertainty Management in National Greenhouse Gas
Inventories. Good practice guidance report accepted by the IPCC Plenary at its 16th session held
in Montreal on 1–8 May 2000; Corrigendum, 15 June 2001. Intergovernmental Panel on Climate
Change, Geneva. www.ipcc-nggip.iges.or.jp/public/gp/english
IPCC (2006). 2006 IPCC Guidelines for National Greenhouse Gas Inventories. Volume 1: ‘General
Guidance and Reporting’, Chapter 3: ‘Uncertainties’. Intergovernmental Panel on Climate Change,
Geneva. http://www.ipcc-nggip.iges.or.jp/public/2006gl/pdf/1_Volume1/V1_3_Ch3_Uncertainties.
pdf (Last accessed 12 September 2014.)
ISO (2002). ISO 6974-6:2002: Natural gas – Determination of composition with defined uncertainty
by gas chromatography. Part 6: ‘Determination of hydrogen, helium, oxygen, nitrogen, carbon
dioxide and C1 to C8 hydrocarbons using three capillary columns.’ First edition, 15 October 2002.
International Organization for Standardization, Geneva, Switzerland.
ISO (2005). ISO 5168:2005: Measurement of fluid flow – Procedures for the evaluation of uncertainties.
International Organization for Standardization, Geneva, Switzerland.
ISO (2006). ISO 14064-1:2006: Greenhouse gases. Part 1: ‘Specification with guidance at the
organization level for quantification and reporting of greenhouse gas emissions and removals.’
International Organization for Standardization, Geneva, Switzerland.
ISO/IEC (2005). ISO/IEC 17025:2005: General requirements for the competence of testing and calibration
laboratories. Current version effective as of 12 May 2005. Revises ISO/IEC 17025:1999. International
Organization for Standardization, Geneva, Switzerland
ISO/IEC (2008). ISO/IEC Guide 98-3:2008: Uncertainty of measurement – Part 3: Guide to the expression
of uncertainty in measurement (GUM:1995). International Organization for Standardization, Geneva,
Switzerland. http://www.iso.org/iso/catalogue_detail.htm?csnumber=50461
Lloyd, S. M and Reis, R. (2007). ‘Characterizing, propagating, and analyzing uncertainty in life-cycle
assessment – a survey of quantitative approaches.’ In Journal of Industrial Ecology, Vol. 11, No. 1,
pp.161–179, January 2007.
Shires, T. M., and Harrison, M. R. (1996). Methane Emissions from the Natural Gas Industry—Volume 6:
Vented and Combustion Source Summary. Final Report, refs. GRI-94/0257.23 and EPA-600/R-96-080f.
Gas Research Institute and US Environmental Protection Agency, June 1996. http://www.epa.gov/
gasstar/documents/emissions_report/6_vented.pdf
Weidema, B. P. and Wesnaes, M.S. (1996). ‘Data quality management for life cycle inventories—an
example of using data quality indicators’. In Journal of Cleaner Production, Vol. 4, No. 3-4, pp. 167–174.
69
IPIECA • API
WCI (2009). Final Essential Requirements for Mandatory Reporting. Western Climate Initiative, 15 July
2009. www.westernclimateinitiative.org/component/remository/Reporting-Committee-Documents/
Final-Essential-Requirements-for-Mandatory-Reporting (Last accessed 12 September 2014).
Williamson, H. J., Hall, M. B. and Harrison, M. R. (1996). Methane Emissions from the Natural Gas
Industry—Volume 4: Statistical Methodology. Final Report, refs. GRI-94/0257.21 and EPA-600/R-96-
080d, Gas Research Institute and US Environmental Protection Agency, June 1996. http://www.epa.
gov/gasstar/documents/emissions_report/4_statisticameth.pdf (Last accessed 12 September 2014.)
WRI/WBCSD (2004). The Greenhouse Gas Protocol: A Corporate Accounting and Reporting Standard.
Revised Edition, March 2004. World Resources Institute and World Business Council for Sustainable
Development. http://www.ghgprotocol.org/standards/corporate-standard (Last accessed 12
September 2014).
WRI/WBCSD (2011a). Greenhouse Gas Protocol: Quantitative Inventory Uncertainty. Guidance for
companies wishing to quantify their inventory uncertainty. World Resources Institute and World
Business Council for Sustainable Development, October 2011. http://www.ghgprotocol.org/files/
ghgp/tools/Quantitative%20Uncertainty%20Guidance.pdf (Last accessed 12 September 2014.)
70
ADDRESSING UNCERTAINTY IN OIL & NATURAL GAS INDUSTRY GREENHOUSE GAS INVENTORIES
71
BG Group Mærsk Oil Santos
BP Marathon Oil Saudi Aramco
Chevron Nexen Shell
CNOOC Noble Energy SNH
ConocoPhillips NOC Libya Statoil
EDF OMV Talisman
eni Petrobras Total
ExxonMobil Petronas Tullow Oil
Hess Petrotrin Woodside Energy
Hunt Oil PTT EP
IPIECA also has an active global network of oil and gas industry association members.
Please refer to our website for a full list.
IPIECA is the global oil and gas industry association for environmental and social issues. It develops,
shares and promotes good practices and knowledge to help the industry improve its environmental and
social performance, and is the industry’s principal channel of communication with the United Nations.
Through its member-led working groups and executive leadership, IPIECA brings together the collective
expertise of oil and gas companies and associations. Its unique position within the industry enables its
members to respond effectively to key environmental and social issues.
IPIECA
5th Floor, 209–215 Blackfriars Road, London SE1 8NL, United Kingdom
Telephone: +44 (0)20 7633 2388 Facsimile: +44 (0)20 7633 2389
E-mail:
Members [email protected] Internet: www.ipieca.org
Addax
© IPIECAPetroleum Husky Energy
2014 All rights reserved. Qatargas
Anadarko INPEX RasGas
Bashneft KPC Repsol
BG Group Mærsk Oil Santos
BP Marathon Oil Saudi Aramco
Chevron Nexen Shell
The American Petroleum Institute is the primary trade association in the United States representing
the oil and natural gas industry,
CNOOC and the only oneSNH
Noble Energy representing all segments of the industry.
ConocoPhillips NOC Libya Statoil
Representing one of the most technologically advanced industries in the world, API’s membership
includes more than 400 corporations
EDF OMV involved in all Talisman
aspects of the oil and gas industry, including
exploration and production, refining and marketing, marine and pipeline transportation and service
eni
and supply companies to the Petrobras Total API is headquartered in Washington, D.C.
oil and natural gas industry.
and has offices in 27 state Petronas
ExxonMobil capitals and provides its members
Tullow Oilwith representation on state issues in
33 states. API provides a forum for all segments of the oil and natural gas industry to pursue public
Hess Petrotrin
policy objectives and advance the interests of the Woodside
industry. Energy
API undertakes in-depth scientific,
technical
Hunt Oil and economic research
PTT EP to assist in the development of its positions, and develops
standards and quality certification programmes used throughout the world. As a major research
institute, API supports these public policy positions with scientific, technical and economic research.
IPIECA
5th Floor, 209–215 Blackfriars Road, London SE1 8NL, United Kingdom
Telephone: +44 (0)20 7633 2388 Facsimile: +44 (0)20 7633 2389
E-mail: [email protected] Internet: www.ipieca.org