Addressing Uncertainty Pilot
Addressing Uncertainty Pilot
ADDRESSING UNCERTAINT Y IN
OIL AND NATUR AL G AS INDUSTRY
GREENHO USE G AS INVENTO RIE S
TECHNICAL CONSIDER ATIONS
AND CALCU L ATION METH ODS
SEPTEMBER 2009
PILOT TEST VERSION
IPIECA is the single global association representing both the upstream and downstream
oil and gas industry on key global environmental and social issues. IPIECA’s programme
takes full account of international developments in these issues, serving as a forum
for discussion and cooperation involving industry and international organizations.
IPIECA’s aims are to develop and promote scientifically-sound, cost-effective, practical,
socially and economically acceptable solutions to global environmental and social issues
pertaining to the oil and gas industry. IPIECA is not a lobbying organization, but provides
a forum for encouraging continuous improvement of industry performance.
5th Floor, 209–215 Blackfriars Road, London SE1 8NL, United Kingdom
The American Petroleum Institute is the primary trade association in the United States
representing the oil and natural gas industry, and the only one representing all segments
of the industry. Representing one of the most technologically advanced industries in the
world, API’s membership includes more than 400 corporations involved in all aspects
of the oil and gas industry, including exploration and production, refining and marketing,
marine and pipeline transportation and service and supply companies to the oil and
natural gas industry.
API is headquartered in Washington, D.C. and has offices in 27 state capitals and provides
its members with representation on state issues in 33 states. API provides a forum for
all segments of the oil and natural gas industry to pursue public policy objectives and
advance the interests of the industry. API undertakes in-depth scientific, technical and
economic research to assist in the development of its positions, and develops standards
and quality certification programs used throughout the world. As a major research institute,
API supports these public policy positions with scientific, technical and economic research.
For more information, please visit [Link].
Telephone: +1-202-682-8000
Website: [Link]
The Oil Companies' European Association for Environment, Health and Safety
in Refining and Distribution
Telephone: +32-2566-9160
Website: [Link]
© IPIECA 2009. All rights reserved. No part of this publication may be reproduced, stored
in a retrieval system, or transmitted in any form or by any means, electronic, mechanical,
photocopying, recording or otherwise, without the prior consent of IPIECA.
Table of Contents
SECTION Page
FOREWORD vi
DOCUMENT AT A GLANCE viii
1.0 INTRODUCTION 1-1
1.1 Importance of Accurate and Reliable GHG Accounting 1-1
1.2 Overview of Uncertainty Terminology 1-2
1.3 Types of Errors 1-3
1.4 Determination of Uncertainty Intervals 1-3
2.0 SOURCES OF UNCERTAINTY 2-1
2.1 Overview of Emissions Inventory Uncertainty 2-1
2.2 Emissions Inventory Uncertainty in the O&G Industry 2-3
2.3 Sources of Measurement Uncertainty 2-5
2.4 Emission Estimation Approaches 2-7
2.5 Inventory Steps and Data Aggregation 2-9
3.0 OVERVIEW OF MEASUREMENT PRACTICES 3-1
3.1 Flow Measurement Practices 3-3
3.1.1 Measurements by Orifice Meters 3-3
3.1.2 Measurement of Flow to Flares 3-5
3.1.3 Example: “Custody Transfer” Measurements 3-7
3.2 Uncertainties of Flow Measurements for GHG Inventories 3-9
3.3 Uncertainties of Sampling and Analysis for GHG Estimation 3-13
3.3.1 Gaseous Samples Collection and Handling 3-13
3.3.2 Quantifying Sampling Precision 3-15
3.4 Carbon Content Measurement Practices 3-15
3.4.1 Laboratory Based Measurements 3-16
3.4.2 On-line Measurements 3-18
3.5 Heat Content Determination 3-19
3.5.1 Direct Measurements 3-19
3.5.2 Computational Methods 3-21
3.6 Laboratory Management System 3-22
4.0 STATISTICAL CONCEPTS AND CALULATION METHODS 4-1
4.1 Measurement Uncertainty 4-1
4.1.1 Precision and Bias 4-1
4.1.2 Confidence Intervals 4-3
4.2 Overview of Uncertainty Propagation 4-4
4.2.1 Propagation Equations 4-6
4.2.2 Correlation Coefficient 4-7
List of Figures
SECTION Page
4-1 Measurement Error Over Time of an Unbiased Estimate 4-2
4-2 Measurement Error Over Time of a Biased Estimate 4-3
4-3 Decision Diagram for Emission Factor Uncertainty 4-15
4-4 Decision Diagram for Measurement Uncertainty 4-17
4-5 Step C – Decision Diagram for Uncertainty Aggregation 4-31
5-1 Onshore Oil Field: Summary of Emissions 5-20
5-2 Onshore Oil Field: Summary of CO2 Equivalent Emissions and 5-21
Uncertainties
The global oil and natural Gas (O&G) industry has been active in promoting consistency and
harmonization for industry greenhouse gas (GHG) emission inventories. Industry associations and their
members have been contributing to the development of guidance for accounting and reporting of GHG
emissions (API/IPIECA/OGP, 2003), and compiling methodologies that are appropriate for estimating
GHG emissions from industry operations (API, 2009). This guidance has been recently augmented with
guidelines to account for reductions associated with GHG projects (API/IPIECA, 2007).
The uncertainties inherent in the data used for emission inventories help inform and improve
understanding for the data’s use. The uncertainty of an O&G company’s GHG emission inventory, or of
its quantified emission reductions, is determined largely by the uncertainties in the estimates of its key
(largest) contributing sources. In turn, each of these uncertainties depends on the quality and availability
of sufficient data to estimate emissions. Our measured data robustness is receiving increased attention to
understanding determinations of GHG emissions and emission reductions.
The American Petroleum Institute (API), the oil companies’ association for the Conservation of Clean Air
and Water in Europe (CONCAWE) and the International Petroleum Industry Environmental Conservation
Association (IPIECA) convened an international workshop on the topic on January 16, 2007, in Brussels,
Belgium. The goals of this workshop were to:
• Develop an understanding of the relative importance of the key factors that contribute to
uncertainty;
• Review and discuss emerging techniques for quantitative assessment of the uncertainty and
accuracy of GHG emissions estimates;
• Identify emission sources and methods where O&G industry efforts are needed to improve
accuracy and reduce uncertainty to acceptable levels; and
• Create a prioritized list of topics to be addressed by the O&G industry to minimize emissions
estimation uncertainty and improve data accuracy.
A summary report, as well as all the workshop presentations, is posted on the IPIECA website
(API/IPIECA, 2007).
The workshop served as the first step in the process of addressing uncertainty and accuracy issues and in
the ensuing industry discussion, a list of priority issues was prepared. This list is comprised of items that
industry experts ought to address in a systematic fashion. As presented in the workshop summary report,
the issues listed by industry members fall into three thematic areas:
1. Measurement methods;
3. External communications.
Industry recognizes the need for meeting regulatory mandates and stakeholders’ expectations, and follow-
up activities will be designed to provide opportunities for continued dialogue and collaborative activities
with stakeholders.
The goal of these guidelines is to augment existing industry guidance and provide technically valid
approaches that would be applicable for use by the global O&G industry to improve GHG emissions
estimation robustness and data quality. The aim is to summarize in a single document overarching
guidance for meeting data needs of a range of initiatives as well as the requirements of diverse GHG
regimes.
These technical considerations and statistical calculation guidance are not designed to be an industry
standard. They will merely provide a brief summary and an overarching discussion of possible
approaches for addressing each of the topics covered. These include: clarification of the sources of
uncertainty in entity GHG inventories; information on measurement practices and their associated
uncertainties; and explanation of statistical procedures that can be used to quantify uncertainties. This
document also contains several case studies that illustrate how to implement the recommended
approaches.
These emission inventories of typical oil and natural gas (O&G) operations are quite complex; they are
based on a combination of measured and estimated emissions data, according to local requirements and
available information. The overall range of uncertainty associated with an entity GHG inventory is
determined primarily by the uncertainty associated with the largest (“key”) sources of emissions. In turn,
the confidence interval associated with each individual source depends on the availability of sufficient
data to estimate emissions, or on the quality of that data, in order to properly account for emission
variability.
Uncertainty analysis is a potential tool to not only assess confidence intervals, but more importantly, to
allow the targeting of specific areas for enhanced data collection. Such an analysis will enable a user to
rank order the importance of different emission sources in terms of their overall contribution to the
emissions inventory and its overall uncertainty range.
This document is a companion to the API Compendium of Greenhouse Gas Emission Methodologies for
the Oil & Gas Industry (2009 API Compendium). It provides a range of background information on
industry practices and specific calculation methods that would enable inventory developers to quantify
and better understand the uncertainty associated with the resultant GHG emissions.
Section 1.0 introduces some basic concepts and terms that provide a foundation for understanding GHG
emissions inventory uncertainty. This terminology is used throughout the document. This section covers:
the importance of reliable GHG accounting; a terminology overview; definition of error types; and a
description of the determination of uncertainty ranges (also known as confidence intervals).
Section 2.0 discusses the major sources of uncertainty in GHG inventories. It moves from general
concepts to issues that are germane to GHG inventories in the O&G industry. It also describes factors
that could introduce errors into the emission measurements process and contribute to the range of
uncertainties of estimated emissions. It introduces the categories of emission estimation approaches and
their uncertainty implications, and concludes with a short description of emission inventory steps and data
aggregation.
Section 3.0 provides an overview of measurement practices, focusing on gas flow measurements and the
determinations of carbon content and heating values of combusted fuels. The section recognizes industry
recommended practices and standards that have traditionally applied to “custody transfer”. This section
Section 4.0 provides calculation methods for quantifying GHG inventory uncertainty. It addresses the
statistical characterization of measurement uncertainty, and the uncertainty associated with published
emission factors, and concludes by providing specific statistical methods and industry relevant examples
for the propagation of uncertainty. This section links the statistical concepts introduced in earlier sections
to their application for measurement uncertainty and error propagation. Decision trees are provided to
guide the user through the steps needed to quantify uncertainty for each part of a GHG inventory.
Section 5.0 takes the guidelines, procedures, and equations for calculation of uncertainty that were
outlined in Section 4 and applies them to a hypothetical facility to demonstrate how to calculate total
uncertainty for a facility-level GHG inventory. The statistical approaches discussed previously are also
applied to two select refinery operations to examine and compare uncertainty estimates for different
emission estimation methods. The hypothetical example facilities and operations are taken directly from
the 2009 API Compendium containing details of how the inventory was generated or calculated. This
section also provides a detailed analysis of the resulting inventory data and its uncertainty in the context
of improving GHG data quality and reducing uncertainty.
Appendix D – Measurement method summaries for carbon content measurement methods, and heating
value measurement methods;
Appendix E – Unit conversions including energy units, common units of measure for fossil fuel heating
content values, and carbon content of selected fuels; and
Appendix F – Calculation details for uncertainty estimation for an example GHG inventory.
• Determination of the uncertainties associated with the individual measurements and factors used
in constructing the emissions inventory, and
• Propagation and aggregation of these individual terms to derive uncertainty intervals (at a pre-
designated probability level) for the whole inventory.
Clearly, the extent and scope of such analysis will depend on the likely uses of this information. For
example, the uncertainty analysis required for data that are merely used for relative ranking or comparison
of trends would be different than that required to demonstrate attainment of GHG emission limits, or
progress made towards meeting GHG emission reduction targets.
Higher quality GHG data leads to higher certainty of emission assessments, and improved confidence in
the data reported. This is true for national and government assessments, and is also important at the
entity, or facility level. To ensure that a company’s strategies and forward-looking actions are based on
the most robust data set and most appropriate computational methods, it is important that this data set and
method be based on four key factors (“The Four C’s”).
− Spurious errors;
− Systematic errors; and
− Random errors.
One or all of these errors could be associated with individual measurements or input variables used for
deriving an emissions inventory. However, such individual error determinations should not be confused
with the overall inventory uncertainty.
Indicators such as the range, confidence interval, or other error bounds typically are used to quantify an
emission estimate uncertainty. Errors may be due to the inherent variability of the emission processes and
the bias–or imprecision–in the terms typically used to define them. Bias is the result of a systematic error
in some aspect of the emissions inventory process. In contrast, imprecision is due to random errors or
fluctuations in the measurement process.
The limits of the confidence interval associated with GHG emissions from a source are directly dependent
on the probability distribution, or the probability function, used to represent that data set. For some
probability distributions, there are analytical relationships that relate the standard deviation to the required
confidence intervals. For example, when a normal distribution is assumed for the variable under
consideration, the confidence limits would be symmetric about the mean, and for a 95% confidence
interval the confidence limits are approximately 2 standard deviations above and below the mean.
Hence, the quantification of uncertainty intervals for calculated GHG emissions will depend both on the
accuracy and representativeness of measurement data used and the assumed distributions of other key
parameters used in the computations. The uncertainties associated with both emission factors and activity
data could be best described by probability density functions that are constructed from available data.
The applicable shapes of these probability density functions could be either determined empirically or by
expert judgment, following procedures described in many guideline documents and standards (ISO, 2005;
IPCC, 2000).
This section has introduced the basic concepts that are germane to estimating the uncertainty range of
GHG emissions. The applicable statistical calculation procedures are presented in greater detail in
Section 4.0.
− Spurious errors, which may be due to incomplete, unclear, or faulty definitions of emission
sources that result from human error or machine malfunction;
− Systematic errors, which may be due to the methods (or models) used to quantify emissions for
the process under consideration; and
− Random errors, which may be due to natural variability of the process that produces the
emissions.
When assessing the process or quantity under consideration, uncertainties might be associated with one or
more factors such as: sampling, measuring, incomplete reference data, or inconclusive expert judgment.
Uncertainties due to models or equations are related to the proper application of estimation methodologies
to the respective source categories. These errors are typically eliminated as far as possible in advance,
Adhering to appropriate sampling, measurement, and estimation procedures – with applicable quality
control and quality assurance measures – can help minimize uncertainties. Nonetheless, it is the nature of
the measurement process that makes it impossible to measure a physical quantity without error.
It is this collection of measurement errors and approximations that contribute largely to overall emission
inventory range of uncertainty. Emission inventory development and uncertainty analysis could be
viewed as components of a quality improved management process. Once estimates of uncertainty are
developed, the inventory preparer could review the inventory and target the most significant sources –
those exhibiting the largest range of uncertainty – for more research and refinement. Table 2-1 provides
an overview of selected methods recommended by the U.S. Environmental Protection Agency (EPA) for
qualitative and quantitative estimation of emissions ranges of uncertainty. The uncertainties associated
with natural variability inherent to the emission process and its underlying data can be assessed by
statistical analysis methods. Further discussion and specific procedures for statistical quantification of
uncertainty intervals are provided in Section 4.0.
LEVEL OF
METHODOLOGY DESCRIPTION OF METHOD
EFFORT
Qualitative − Sources of uncertainty are listed and discussed. Low
Discussion − 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 Low
Quality Ratings emission factor or parameter.
Data Attribute − Numerical values representing relative uncertainty are assigned through Medium
Rating System objective methods.
(DARS)
Expert − Experts estimate emission distribution parameters (i.e., mean, standard Medium
Estimation deviation, and distribution type).
Method − Simple analytical and graphical techniques are then used to estimate
confidence limits from the assumed distributional data.
− In the Delphi method, expert judgment is used to estimate uncertainty
directly.
LEVEL OF
METHODOLOGY DESCRIPTION OF METHOD
EFFORT
Propagation of − Emission parameter means and standard deviations are estimated using Medium
Errors Method expert judgment, measurements, or other methods.
− Standard statistical techniques of error propagation typically based upon
Taylor’s series expansions are then used to estimate the composite
uncertainty.
Direct − Monte Carlo, Latin hypercube, bootstrap (resampling), and other High
Simulation numerical methods are used to estimate directly the central value and
Method confidence intervals of individual emission estimates.
− 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 − Direct or indirect field measurements of emissions are used to compute High
Indirect emissions and emissions uncertainty directly.
Measurement − Methods include direct measurement such as stack sampling and indirect
(Validation) measurement such as tracer studies.
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,” U.S. EPA, Research Triangle Park, NC, July 1996
2.2 Emissions Inventory Uncertainty in the Oil & Natural Gas Industry
The API Compendium (API, 2004) and its forthcoming revision (API, 2009) provide an extensive
compilation and tabulation of methods that are used by companies in all the sectors of the O&G industry
to calculate their GHG emissions in a consistent manner. The API Compendium provides a wide range of
emission factors and other emission estimation methods that are directly applicable to all sectors of
industry operations. It also exhibits the ranges of uncertainty (at the 95% confidence level) associated
with each of the case study examples featured in Section 8 of the 2009 API Compendium. Related
industry guidelines include those of the Interstate Natural Gas Association of America (INGAA), which
provides supplemental guidance for estimating GHG emissions from key emission sources associated
with natural gas storage and transmission sector of the industry (INGAA, 2005) and the AGA which has
also published specific guidelines for estimating GHG emissions from gas distribution operations (AGA,
2008).
O&G industry operations vary widely among its operating sectors due to the nature of the operation, its
geographical locations and local practices. Operations in some of the industry sectors are highly
centralized in large and complex facilities while other extend over large geographical areas, with some of
The uncertainty associated with CO2 emissions from combustion would be primarily attributable to
variation in the composition of combusted fuels and their respective consumption rates (or total volumes).
For quantification of combustion emissions, quality data are typically available for industry facilities in
all sectors, though significant effort may be required in order to collect data for smaller operating
installations that are spread out in multiple locations.
Moreover, since a large fraction of industry operations rely on self-generated, if is the knowledge of the
carbon content of such fuels that is at the root of determining their associated CO2 emissions. For the
exploration and production sector, the composition of these self-generated fuels may vary with the nature
of the producing formations, while for refining it will depend on the composition of the crude oil
processed and the slate of products manufactured. On the other hand, for natural gas transmission and
distribution operations, gas quality and its composition are expected to adhere to contract requirements
and would vary only within a narrow specifications range. Hence, the use of average fuel compositions
data has to be evaluated when compiling an emission inventory since it might result in wide uncertainty
ranges for some sectors and operations, while they might be perfectly acceptable for others.
A different set of parameters is important for understanding emissions associated with process vents and
fugitive emissions. For many of the large process units that are found in refineries and natural gas
processing plants, numerical models (equations) are available for estimating these emissions. For high-
pressure pipelines transmitting natural gas over long distances, the main GHG emissions are due to
reciprocating engines and turbines, venting due to gas blow-down, and fugitive emissions associated with
leaking piping components. For low-pressure gas distribution, most of the GHG emissions come from
compressors and leaks from gas distribution mains and associated equipment. Quantifying emissions, and
their associated uncertainty ranges, for venting and fugitive emissions in the exploration and production
sector poses a real challenge. These emissions could be quite significant for high-pressure uncontrolled
natural gas production operations or miniscule for controlled oil and associated gas production.
For vented emissions, operators in the U.S., as well as other jurisdictions, typically maintain required
records for reporting (and archiving) gas-venting incidents. However, U.S. reports of “lost and
unaccounted for gas” from natural gas pipelines typically account from both vented and fugitive
emissions. Therefore, disaggregating the data would be required in order to derive a separate average
emission factor for venting incidents only. When it comes to fugitive emissions from equipment leaks,
In summary, since the most prevalent emissions from fuel combustion is CO2, and from venting and
fugitive emissions, CH4, the main contributors to the uncertainty ranges of these respective GHGs in an
inventory generally are:
− For estimating CH4 emissions – Uncertainty is primarily associated with estimates of vented and
fugitive emissions. The records that installations are required to keep (and report) on vented or
released gas might not be similar under all regimes globally. Fugitive emission estimates exhibit
the highest degree of uncertainty due to the use of average emission factors per component,
device or type of operation, and due to improper conversions of existing factors that are expressed
in terms of volatile organic compounds (VOCs) to CH4. Since the CH4 to VOC ratio varies
among installations, or even within different parts of a processing plant, these average emission
factors, coupled with generic conversions from VOC to CH4 may not be the best representation of
CH4 emissions.
The measurement process is comprised of different steps and each can introduce uncertainty into the final
results. These steps can be applied whether the measurement process involves measurements of activity
data (flow volumes), fuel carbon content speciation, screening for fugitive emission leaks, or direct
emissions testing.
The sources of uncertainty discussed below range from methods choice to physical constraints of the
measurement process itself and the processing of the data obtained.
b. Calibration uncertainty
Manufacturers and companies typically calibrate measurement instruments before they are used in the
field or in a plant. However, instrument calibration needs to be checked periodically to detect
instrument drift and thus reduce measurement uncertainty. The calibration process could achieve that
goal if the process is traceable to a known reference standard and allowance is made for adjustments
of the measurement instrument if a bias is detected during the calibration process. .
Quality control techniques and robust data management practices can minimize the effects of many of
these uncertainty sources. For example, comparing known input values with their measured or
computed results can provide an estimate of the data acquisition uncertainty. However, it is not
always possible to do this in practice. In these cases, it is necessary to evaluate each potential error
and aggregate these errors to assess the overall uncertainty.
Emissions information is typically obtained either through direct on-site measurement of emissions or by
using an emission estimation equation or model. Emission estimates are used for facility permitting,
development of control strategies, compliance review, and demonstration of attainment of environmental
goals. The four basic approaches for estimating emissions are:
a. Emissions Factors
Emissions factors are at the core of calculating GHG emissions when compiling O&G industry
inventories. An emissions factor relates the rate of emission of a specific compound with an activity
rate associated with its release.
The general equation for using emission factors for calculating emissions is:
In practice, for estimating GHG emissions from the O&G industry operations this means:
− For CO2
When calculating emissions and their associated uncertainties, it is important to note that the overall
uncertainty is based both on activity data (process flow, throughput, usage, or equipment count) and
on the emission factors used. Each contributor to uncertainty should be assessed independently and
then aggregated in the final analysis, as shown in Figure 4-3 in Section 4, and in the associated
examples provided.
Over the years, U.S. EPA and other emission factors repository databases have provided average
emission factors that can be used broadly across many industry source categories and operations. The
published emission factors are typically accompanied by a description of the group of processes and
the conditions they represent. Authoritative factors are generally published by EPA in AP-42 (U.S.
EPA, 1995 and further updates), or by the EU EMEP/CORINAIR Emission Inventory Guidebook
(EMEP/CORINAIR, 2007). The IPCC has also launched a new Emissions Factors Database (IPCC,
2006) for use with GHG emission calculations.
This 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 an
applicable monitor at the source. The error associated with these determinations varies for different
compounds and it is not necessarily a more reliable method for measuring emissions especially since
it is not available for monitoring all GHG emissions. Additionally, the acquisition and operating
costs of continuous emission monitoring are very high and thus this technique has limited application
to the evaluation of GHG emissions.
c. Source Testing
Like continuous emissions monitoring, source-testing data may be generated by either extracting a
sample or placing a monitor at a source, followed by analysis to characterize the emitted species. In
this application, the measurement campaign is limited to a specified number of hours. Facilities then
use the average emission rate calculated from source testing 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.
This measurement approach is generally useful when appropriate test methods are used, and emission
and process data are collected at a frequency that allows good characterization of emission
variability. However, the cost of increased measurement frequency should be balanced with the
contribution of the tested source to the overall inventory and the incremental improvement in the
range of uncertainty that is attainable by this increased testing.
For some emission sources, a material balance (based on fuel flow and carbon measurements) may be
an appropriate means of estimating GHG emissions. Use of a material balance requires knowledge of
total flows or throughput rates, and the corresponding compositions of those streams. Material
balance approaches can also be used for assessment of evaporative losses or for simulation of process
emissions under defined conditions.
Clearly, various methods are applicable to quantifying GHG emissions. Emission estimation methods
range from simple activity measurements multiplied by applicable emission factors to more sophisticated
estimation algorithms. The advanced methods represent an integrated approach that relies on the use of
factors and other data, including generic process simulation models, source specific models, and species
profiles databases.
In GHG emission inventories, emission estimates are obtained from many intermediate and independent
results, each of which is calculated from a separate set of data that is characterized by a different range of
uncertainties. The compilation of an entity-wide GHG emissions inventory typically follows a sequence
of steps:
Establishing boundaries – Where the organizational and operational boundaries are defined (for a
first-time inventory), or examined (for recurring cycles), this step will be largely dictated by local
requirements or corporate policies. It might involve facility-by-facility assessment prior to
aggregation, or it could use other pertinent entity indicators and information;
Collecting and inputting data – Where the activities data are collected and archived based on the
boundaries established above. The data are then incorporated into appropriate calculation tools for
emission calculations. The level of ‘granularity’ of the data collected and the details of the
calculation methods are dictated by local requirements with industry guidance (API Compendium) as
a resource to provide relevant technical details.
Validating data compiled – Where various techniques are used to compare the new data with earlier
versions (if available) to identify potential large errors. These errors could include: either large
changes or unchanged activity data for given facilities; operations that are not accounted for; lack of
supporting data for measurements or emission factors used; erroneous units or unit conversions; and
more.
Finalizing the inventory – Where the quality-checked and validated data are aggregated for
reporting based on company policy or local requirements. The preferable way of reporting the results
is in terms of the total emissions for each of the GHG species, along with the global warming
potential weighted sum of these emissions (also known as the CO2-E emission).
The overall uncertainty range for each GHG species, and CO2-E, should also be reported with the total
emissions in the format of:
Section 3.0 provides an overview of measurement practices that would result in high quality data when
properly implemented. It focuses on measurements that are applicable to the key sources that contribute
to the overall GHG emissions inventory, i.e., CO2 emissions from combustion devices.
Section 4.0 provides more detailed guidance on the calculation methods that are applicable for defining
single source uncertainties and for aggregating them at the facility (or entity) level. This section provides
a range of examples that are directly applicable to O&G facilities and their GHG reporting needs.
API publishes one of the more comprehensive sets of custody transfer measurement standards, but it is
neither complete nor the only widely recognized source for such industry practices. API’s MPMS
includes over 140 titles, and API publishes approximately eight new or revised measurement standards
each year. Appendix B presents a comprehensive list of specific measurement standards (and their
respective editions) from several standards setting organizations, which could be used to support the
Custody transfer measurements are typically expected to meet a set of performance specifications. For
other measurements, such as those performed to support the development of a GHG emissions inventory,
data quality objectives should be established prior to initiating any data collection in order to ensure that
the uncertainty ranges of the measured quantities are consistent with the intended use of the data.
Throughout this section we provide references and describe a select subset of industry standards that are
most typically used for the respective measurements discussed. Even so, the full list of measurement
methods provided in Appendix B could be used to provide equivalent measurements to meet company
practices and available instrumentation.
For many O&G industry installations, CO2 emissions from combustion and flaring are the largest
contributors to overall GHG emissions. Therefore, this section focuses on measurements and methods
typically used for quantification of these CO2 emissions. The subsections below provide details on flow
measurement practices, and their associated uncertainties, as well as methods for the measurement of
carbon content and heating values of combusted fuels.
For some other O&G industry sectors, such as exploration, production, processing, transmission and
distribution operations, emissions from other GHGs (such as methane) can contribute significantly to a
facility’s GHG emissions. Available emission factors exhibit large uncertainties and might not be
representative of current operating practices. Therefore, for inventories where CH4 constitutes a larger
fraction of the emissions, the overall uncertainty would be expected to be substantially higher.
Industry is now more fully internalizing the potential impact of measurement errors and bias through the
full chain of emission calculations including measurement equipment, software calculations, simulation
models, and the limitations of reliance on existing emission factors on the uncertainty range of resultant
GHG emissions inventories. Considerations of the need for more representative measurement data, and
the assessment of equipment design and age, are gaining more prominence in the process of assembling
high certainty emission inventories. This section focuses on measurements that are pertinent to improved
characterization of combusted fuels and the quantities used, but might not be directly applicable to
measurement of leakages and fugitive emissions. The measurement practices highlighted here represent
an initial step in what could end up being a multi-year effort to improve measurements of GHGs and
quantify the activities that cause their emissions.
Continuous handling of very large liquid and gas flow volumes is a characteristic of all the sectors of the
O&G industry. Therefore, an in-depth understanding of flow measurement is essential both for internal
process control and for transferring “custody” of intermediate streams or finished products. The accuracy
of measurements of “custody meters” is historically quite high, and practices follow rigorous industry
standards.
Industry has been instrumental in developing international voluntary standards such as ISO 5168 (see
Reference 8) establishing general principles and describing procedures for evaluating the uncertainty of
measuring fluid flow rate or quantity. Annex A of ISO 5168 provides a step-by-step procedure for
calculating and reporting these measurement uncertainties. Similarly, Chapter 14 of the API Manual of
Petroleum Measurement Standards (MPMS), contains detailed procedures and practices for all aspects of
natural gas (and similar gases) fluids measurement and calculation of their associated uncertainties, at the
point where oil or gas enters the marketplace (“custody transfer”) (API MPMS, 2006). Those same
practices are not as rigorously applied to internal accounting and process control during normal
operations.
Recommended practices for the installation, calibration and calculation of flows for these custody meters
are provided in Section 3 of Chapter 14 of API’s MPMS (API, 2005). This standard was developed by a
collaborative effort by members of API, AGA, and the Processing Gas Association and with contributions
from the Canadian Gas Association, American Chemical Council, the European Union, Norway, Japan
and others. It is designed to ensure global consistency for O&G transactions. The four-part standard for
square-edged, concentric orifice meters consists of:
The standard recognizes that many factors contribute to the overall measurement uncertainty associated
with many metering applications, as summarized in Exhibit 3-1.
In the reapproved 2006 version of the standard, several changes were incorporated to reduce the
uncertainty attributable to installation effects and to improve the rigor of the flow calculation routines.
The revised standard recognizes the lead-time necessary for upgrading existing installations, and leaves
this lead-time to the discretion of facility operators and their data quality targets for flow measurement
data.
However, it should be recognized that if orifice meter installations are not upgraded to conform to the
new recommendations, measurement bias error may occur. This bias might be due to improper upper and
lower distances from bends and points of flow turbulence that might lead to inadequate flow conditioning
prior to measurement. Additionally, even without changing equipment installations, the standard
recommends adopting new calculation procedures and techniques (explained in Part 1 and 3 of the
standard) that represent significant improvements over the previously adopted approach. It is important
to note that the expected uncertainty ranges for flow measurements quoted in Part 1 of the reaffirmed
standard may differ from those obtained in practice when the equipment installation differs.
The measurement of flow to flares is distinctly different than other flow measurements. Flares are
designed as safety relief systems and typically are capable of handling highly variable flow rates of
widely varying gas compositions. Therefore, some of the practices that are generally applicable to
custody transfer or process control flows have to be modified when addressing flows to flares. API
published a measurement standard addressing gas or vapor flare flow measurements, which also includes
cautionary details about the effects of fouling (due to entrained liquid droplets, aerosol mists, or other
contaminations) on the measurement (API MPMS, July 2007).
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 from
molecular weights approaching that of hydrogen to molecular weights of C5+.
As with other flow measurements, the accurate determination of flow to flares is dependent on many
parameters such as the ability to predict – or measure – mixture composition, pressure, temperature,
and/or density. The accuracy of measurements associated with highly variable flare gas mixtures will
depend largely on the meter technology type and the ability of the flare flow measurement system
(FFMS) to achieve the targeted response time and analytical accuracy levels. Exhibit 3-2 below lists the
Variability in flare composition may also be a significant factor in determining the measurement
uncertainty of an FFMS. Knowledge of flare composition may have a major effect on the calculation of
the actual volume, standard volume, or mass measured by the meters. For example:
− For a differential pressure meter – the output is a function of the square root of the flare gas
density.
− For thermal flow meters – knowledge of the actual volume (or standard volume) requires
consideration of the compositional effect on thermal conductivity and dynamic viscosity.
− For ultrasonic flow meters – sound speed is a function of gas composition.
Converting actual flow to mass flow requires the knowledge of gas composition in order to derive gas
density. When an analyzer is incorporated in the measurement system to correct for composition, care
must be taken to ensure that the response time of the system is short compared to the upset flow event
during flaring to ensure representative measurement during actual flaring.
An example of how all of these elements would be combined into the overall measurement uncertainty is
provided in Table 3-1.
“Custody transfer” measurements are defined as measurements that provide quantity and quality
information, which can be used as the basis for a change in ownership and/or a change in responsibility
for materials. In most O&G producing jurisdictions around the world national regulations and directives
have emerged to specify requirements for the expected accuracy and uncertainty ranges associated with
“custody transfer” and other critical measurements. For such precise metering applications, the
flowmeters and adjacent piping used in the measurement system are expected to meet the requirements of
the relevant, preferably the most stringent, specifications of the API and ISO standards that are cited in
many national regulations.
One such example of the measurement requirements promulgated for O&G operations are those of the
Alberta Energy Resources Conservation Board (ERCB), as shown in Exhibit 3-3.
a.) Directive 017 (EUB/Board, May 2007), spells out the measurement requirements for custody transfer within the Upstream
Oil and Gas operations in Alberta, Canada. The ERCB standards are stated as “maximum uncertainty of monthly volume”
and/or “single point measurement uncertainty”, as listed in Tables 3-2 for Oil Systems and Gas Systems measurements,
respectively.
b.) The uncertainties are to be applied as “plus/minus” (e.g., ±5%), and only measurements at the delivery or sales points are
required to meet the highest accuracy standards since they would have a direct impact on royalty determination.
The directive makes it clear that other measurement points that play a role in the overall control and accounting process would
be subject to less stringent accuracy standards. These less stringent accuracy standards are designed to accommodate physical
limitations and/or the overall economics of achieving very stringent accuracy standards for each volumetric measurement.
In Europe, the European Union Emissions Trading System (EU-ETS) specifies a tiered approach for
emission calculations together with required uncertainty ranges, according to the Monitoring and
Reporting Guidelines (EU-ETS MRG, 2007). It sets up a matrix of uncertainty requirements for different
facility sizes and measurement approaches used.
− For facilities emitting between 50,000 to 500,000 tonnes of fossil CO2 – Uncertainty ranges
specified are ±7.5%, +5%, and ±2.5% when the facilities employ tiers 1, 2, and 3 calculation
approaches, respectively;
− For facilities emitting over 500,000 tonnes of fossil CO2 – Uncertainty ranges are expected to
be as low as ±1.5% for facilities required to employ Tier 4 approaches.
It is important to note that the EU-ETS requirements are applicable to a limited set of O&G industry
installations and that facility inventories are being tracked only for CO2 emissions from fuel combustion
and flaring. The requirement to quantify these sources within such tight uncertainty ranges is a reflection
of the fact that these are the sources for which appropriate emission calculation methods are available,
while they are also the largest emission sources and the key contributors to most installations’ GHG
emissions.
When measuring fuel flow rate, or its total volume, and using the information to calculate GHG
emissions, it must be determined whether the flow meters used are properly installed and calibrated, and
that they are capable of providing data that are within the uncertainty ranges required by the governing
climate program. Differences must be considered between the manufacturers’ specifications of flow
meters’ expected measurement errors and those that are attained when using the flow meters in the field.
It is common practice to test flow meters in a laboratory setting under controlled conditions, prior to field
installations. However, these laboratory bench tests typically do not simulate ”real world” variations in
fluid flow and other possible fluctuations, and drift of the entire measurement system. For any given
As an example, Table 3-3 provides a listing of different meter types, their applicable fluid medium, and a
brief description of their operating principles. The table also lists manufacturers’ specified instrument
errors, as provided in a survey conducted by the EU-ETS Support Group (ETSG, 2007). The information
provided is an indication of potential error ranges and not the expected uncertainty ranges obtained in
practice. Appendix B provides additional details about the manufacturers’ recommendations for the
installation, calibration, and maintenance of these flow meters. The error levels cited refer primarily to
random errors that are observed under ‘ideal’ laboratory conditions and that decrease with repeated
measurements. They do not properly account for systematic errors (or bias) where the errors are due to
improper installations, inadequate calibrations, or devices drift, as discussed above.
However, these manufacturers’ specified measurement errors might not be attainable due to the practical
operational limitations of the facility. In most O&G industry facilities, detailed inspections, maintenance,
and recalibration of process control flowmeters are possible only once every few years when process units
are shut down for scheduled turnaround.
As indicated in Section 3.1.2 (Exhibit 3-2), a few key steps that should be followed to determine the
combined uncertainty associated with individual flow measurements. An expanded list of factors that
ought to be considered when evaluating the uncertainty of flow measurements that are used for GHG
emission calculations are provided in Exhibit 3-3.
The discussion above pertains to single flow measuring device uncertainty. Additional analyses are
required to convert this individual measurement uncertainty to an assessment of the uncertainty range of
annual measurements. Detailed guidance on the calculation steps is provided in Section 4.
The emission measurement process comprises either direct measurement at the source level of collection
of samples and their analysis in the laboratory to determine mass emissions. Sampling and analysis are
part of the same measurement process and their combined contribution to emissions estimation
uncertainty is obtained by their combined variances, as detailed in Section 4.0. Emission measurement
uncertainties for processes of sampling and analysis depend on random errors, measurement precision,
and systematic errors or bias.
The API MPMS provides specific details for collecting and handling natural gas samples for critical
measurements such as custody transfer (API, February 2006). Exhibit 3-4 provides guidance on general
considerations of inaccuracies that might be introduced in the measurement system when collecting
samples that are used for carbon content and/or heating values determinations for GHG emissions. At the
same time, designing a sampling and analysis system should – in all cases – take into account regulatory
requirements and contractual obligations.
EXHIBIT 3-4: KEY FACTORS IMPACTING 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, are typically 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 documented by
industry recommended practices for sample families
− Analyses are typically limited to the range of mixture concentrations and species previously identified
− For complicated sample matrices, the potential for interferences is usually noted.
Quantifying sampling precision requires that primary samples be collected according to a defined
protocol, but randomized in some way for each sample (in either space or time). For example, a pre-
selected percentage of the total number of samples can be collected in duplicate and the repeatability of
the measurement determined from these duplicate samples. Additionally, duplicate gaseous samples can
also be analyzed in duplicate and thus a full record of both sampling and analysis system variations can be
obtained. In this case, the sampling component of the variance represents the natural variability of the
sampling target as well as any errors in the sample collection and preparation.
When in situ measurement techniques are used (e.g., infrared gas analyzer), both the sampling and
analysis are addressed at the same time. The determination of sampling bias, or the difference between
the mean of several measurements and the true value, is more difficult. Biases could arise from several
causes such as sample fouling, inappropriate handling, or unrepresentative sampling. The “true” value for
the concentration of unknown species in a sample is never known since it is impossible to construct a true
reference laboratory standard for an unknown mixture of gases.
The importance of sampling uncertainty is a relatively new concept whose importance is slowly
beginning to be recognized. Recent research has shown that sampling uncertainty is often far greater than
analytical uncertainty. Therefore, combining sampling and analytical uncertainties to provide an estimate
of measurement uncertainty is an important component of quantifying the overall uncertainty of GHG
estimations that rely on sampling gaseous fuels of varying composition.
Different types of gas chromatography (GC) systems are normally used to analyze the carbon content of
gaseous streams. The GC systems might be laboratory based or set up as an online device for automated
collection of samples and their analysis. The systems are typically set up to analyze the individual
components in the sampled gas and provide detailed reports of properties including composition, calorific
value, and density.
The results of the determination of individual – or groups – of carbon-containing species are then used to
assess the total emissions of CO2 upon combustion of such a fuel. Several key considerations include:
− If the method is capable of determining CO2 content with the rest of the carbon containing
species, no further correction of the carbon content data is required in order to properly account
for all CO2 emissions;
− If the method is capable of a quantitative determination of CH4 content, these data can be used
separately for calculating evaporative and processing leaks along with venting losses; and
− All carbon content measurement data should be used in conjunction with the applicable fuel flow
measurements when calculating emissions.
Several ASTM and ISO methods are available for determining the composition and carbon content in the
natural gas range as well as for liquid and solid fuels. Although the terms “repeatability” and
“reproducibility” are applicable to all measurement situations (see Appendix A for terminology and
definitions), they are used here in the context in which they are defined in ASTM standards:
− Repeatability is defined as the difference between two successive results obtained by the same
operator with the same apparatus under constant operating conditions on identical test materials,
and
− Reproducibility is the difference between two results obtained by different operators in different
laboratories on identical test materials.
When using a natural gas measurement method for refinery fuel gas, care should be taken to ensure that the
range of compositions of individual components is within the bounds specified by the method. Table 3-4
lists selected commonly used laboratory measurement methods that are frequently used for the
determination of fuel carbon content. Additional details about these methods and further discussions of
their main features are provided in Appendix D.
ASTM D5291 - 02 Standard Test Methods − Applicable to samples such as − NOT recommended for
for Instrumental crude oils, fuel oils, additives, and the analysis of volatile
(2007)
Determination of residues for carbon and hydrogen materials such as
Carbon, Hydrogen, and and nitrogen analysis gasoline, gasoline-
Nitrogen in Petroleum − Tested in the concentration range oxygenate blends, or
Products and Lubricants of at least 75 - 87 wt% carbon, at gasoline type aviation
least 9 - 16 wt% hydrogen, and turbine fuels
0.1 - 2 wt% nitrogen
a
Measurement precision data is provided as an indication of attainable precision if all steps of the methodology are adhered to as described in the methods
procedures.
Instrumentation in this field has been developed to provide a measurement of stream components in order
to achieve optimum control and assure product quality. The configurations of such analyzers are
customized to accommodate typical site parameters and operating practices. Most such analyzers are
designed with ASTM and ISO standards in mind and their calibration routines are designed to provide
both reported data and its associated uncertainty. As mentioned previously, the two primary applications
are for refinery gas analyzers and natural gas analyzers, and these are discussed briefly below.
The heating value or calorific value of a substance is the amount of heat released during the combustion
of a specified amount of the substance. The calorific value is a characteristic for each substance and is
measured in units of energy per unit of the substance, such as: kcal/kg, kJ/kg, J/mole, Btu/m³. The heat
of combustion for fuels is expressed as:
− HHV – higher heating value (or gross calorific value or gross energy or upper heating value).
This value is determined by bringing all the products of combustion back to the original pre-
combustion temperature, and in particular condensing any water vapor produced.
− LHV – lower heating value (or net calorific value) is determined by subtracting the heat of
vaporization of the water vapor from the higher heating value and treating any water formed as a
vapor. The energy required to vaporize the water therefore is not realized as heat.
In the O&G industry, the two most prevalent modes of determining heating values of gaseous fuels is
either by measuring it directly, which can be accomplished either by stoichiometric combustion or by
calorimetric techniques, or by computational methods that are based on standardized calculation
procedures using gas sample composition data. A brief summary of these two types of practices is
provided below.
The heating value indicates the amount of energy that can be obtained as heat by burning a unit of gas.
The heating values of a gas depend not only upon the temperature and pressure, but also upon the degree
of saturation with water vapor.
As mentioned above, general practices for determining fuel gas heating values rely on either calorimetric
techniques or stoichiometric combustion practices. Table 3-5 provides a listing and a brief description of
some selected methods for heating value determination for gaseous, liquid, and solid fossil fuels. Further
discussion of the main features of these methods can be found in Appendix C.
The ideal gas heating value and ideal gas relative density are calculated from the molar composition and
the respective ideal gas values for the components; these values are then adjusted by means of a
calculated compressibility factor.
The ASTM approach recognizes that the calorific value of a mixture, such as in refinery fuel gas systems,
is a function of the mol% composition of the individual components of that mixture. For refinery fuel
gas, this mixture could contain both carbon containing species, which upon combustion will contribute
directly to CO2 emissions, as well as non-carbon containing species. Hydrogen, for example, is an
important contributor to refinery fuel gas heating values but it does not contribute to CO2 emissions when
combusted.
To implement this practice, the user would first determine the molar composition of the gas in accordance
with any applicable ASTM or GPA method. For a precise calculation, at least 98 % of the sample
constituents should be determined as individual components, with no more than a total of 2 % in terms of
groups of components (e.g. butanes, pentanes, hexanes, and so forth).
An ideal combustion reaction for fossil fuels (that may contain hydrogen) in the ideal gas state can be
generally represented as:
⎛ b c⎞ ⎛b ⎞
Ca H b + cH2 + ⎜ a + + ⎟ O2 = aCO2 + ⎜ + c ⎟ H 2 O (Equation 3-1)
⎝ 4 2⎠ ⎝2 ⎠
The ideal net heating value is the heating value that is observed when all the water remains in the ideal
gas state, while the ideal gross heating value is observed when all the water formed by the reaction
condenses to liquid. The difference between them is the enthalpy of vaporization of the water formed
during the combustion process. Therefore, the ideal gross heating value for a mixture can be expressed
as:
∑x
j =1
j × M j × H mj
Hm = n
(Equation 3-2)
∑x
j −1
j ×Mj
where
xj = the mole fraction of Component j;
Mj = the molar mass of Component j;
n = the total number of components; and
Hmj = the ideal gross heating value per unit mass for Component j (at 60°F or 15.6°C), as tabulated
in ASTM D3588. Values of Hm are independent of pressure, but they vary with temperature.
Errors that should be considered when computing heating values include errors in the heating values of
the components and in the composition data. The uncertainty ranges of the heating values for the
components cited in the ASTM practice are about 0.03%. These errors may affect the agreement between
calculated and measured heating values; however, those are negligible when compared to the errors
associated with the determination of the composition of individual species. Appendix D contains a listing
of common energy and fossil fuel unit conversion factors that could be used for these and other
calculations discussed in these document.
An additional consideration for uncertainty of GHG measurements is the credibility and technical veracity
of the laboratory performing the requested tests, as specified by different programs. For example, in the
EU-ETS program the MRGs require the demonstration of laboratories management systems (EU-ETS
MRGs, 2007; Section 13.5). Any laboratory used to determine an emission factor, calorific value, carbon
content, or composition data should be accredited according to ISO 17025:2005. If non-accredited
laboratories are used, the EU regulations provide specific requirements for the validation and testing of
such laboratories.
ISO/IEC 17025:2005 (ISO/IEC 17025, 2005) is not intended to be used for overall laboratory
certification, and it does not address compliance with regulatory and safety requirements. It merely
emphasizes the need for a well-developed and communicated laboratory management system that
addresses areas such as: quality, administrative procedures, and technical systems that govern the
operations of the laboratory. The ISO standard highlights the need for laboratories to ensure that their
personnel are aware of the relevance and importance of their activities to the achievement of the
management system’s objectives.
The standard specifies the general requirements for demonstrating competence to carry out specific tests
and calibrations, including field sampling. It covers testing and calibration procedures that are performed
ISO/IEC 17025:2005 is applicable to all laboratories regardless of the number of personnel or the extent
of the scope of testing and/or calibration activities. When a laboratory does not undertake one or more of
the activities covered by ISO/IEC 17025:2005, such as sampling and the design/development of new
methods, the requirements of those clauses would not apply.
As part of compliance with ISO/IEC 17025:2005, laboratories that want to be accredited to this standard
are mandated to seek feedback, both positive and negative, be it from in-house users or from external
customers. The information gathered is expected to help laboratories improve their management systems,
their testing and calibration activities, and customer service. The ISO standard seeks to improve
laboratory measurement proficiency and accuracy by continual improvement of the effectiveness of
laboratory management systems and the implementation of quality policy, quality objectives, internal
audits, data analysis, corrective and preventive actions, and periodic management review.
Activity factors are generally a measured quantity, such as a count of equipment or measure of fuel
consumed. Emission factors may be either based on site-specific measurements or based on published
values that were derived from averaging a variety of measurements. Where measurements are used for
either activity factors or emission factors, two components of uncertainty need to be considered: precision
and bias.
Bias, on the other hand, refers to how accurately a method estimates the true value. An example of bias
would be a meter that is not calibrated correctly and consistently overestimates the measurement. Ideally,
the data measurement scheme should be designed in a way to minimize bias. For example, if the
measuring device is well maintained and calibrated to minimize drift, there may be no bias in the
measurement.
Figure 4-1 shows an example of the error of an unbiased measurement over time, where the errors are
centered around zero. When these measurements are aggregated, some of the positive errors are offset by
some of the negative errors, resulting in a lower estimate of uncertainty for the aggregate measurement.
Figure 4-2 shows an example of the error in a biased measurement over time. Similar to the unbiased
case in Figure 4-1, some of the high errors are offset with some of the low errors so the uncertainty in the
precision of the estimate is lower for the aggregate. However, in Figure 4-2, the errors are not centered
on zero as in Figure 4-1; they are centered on five. This means that there is a bias in the estimate. Unlike
precision, the uncertainty due to the bias generally does not decrease as more data points are collected.
Thus, precision and bias should be considered separately, if possible.
Error 3
0
0 10 20 30
Time
If bias can be quantified, it should be corrected and thus eliminated from consideration in quantifying
uncertainty. For example, if a fuel stream has two types of measurement devices, data can be collected
from both devices to check for agreement. A bias would be indicated if the measurements differed.
However, quantifying and correcting for bias might not always be practical since it often requires
application of frequent calibration routines and implementation of quality control procedures that would
allow instrument adjustments and/or corrections under prescribed conditions. It will also require proper
quantification routines that account for drift or other causes of bias between calibrations on an ongoing
basis.
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.
and
1 n
s= ∑
n − 1 i =1
( xi − x ) 2 (Equation 4-2)
where
xi = the ith observation in the data set and
x = the mean of the data set.
∑x i
x= i =1
(Equation 4-3)
n
Tables for the Student’s t-distribution can be found in most basic statistics references. Most spreadsheet
software programs have a function that will calculate the necessary t-value. This is the preferred method
since the software generally retains more significant digits for the t-value than a look-up table would
display.
The following three assumptions are important when applying the uncertainty propagation technique for
overall uncertainty assessment (IPCC, Section A1.4.3.1, 2001):
1. The uncertainties are relatively small, which is defined as the standard deviation divided by the
mean value being less than 0.3;
2. The uncertainties have Gaussian (normal) distributions; and
3. The uncertainty values (i.e., the errors or uncertainties associated with the measured data or
reported values) are mutually independent.
In many cases, the first assumption may be difficult to meet. For example CH4 and N2O emissions often
have very sparse data and large associated uncertainties. Conducting a Monte Carlo simulation
(discussed further in Section 4.6.1) is an option if the standard deviation divided by the mean is greater
than 0.3. However, Monte Carlo simulations require a significant level of detail for the description data
The second assumption is based on the normality of the distribution of the underlying source data (i.e.,
symmetrical around the mean). According to the Central Limit Theorem, for a large enough sample size
(n>30), we can relax the normality assumption but still assume that the sampling distribution of the
sample means is normally distributed (Casella and Berger, 1990). Hence, if the calculated uncertainty is
based on statistical sampling of the population, one would need to obtain more samples to approach
normality. Alternatively, we might consider data transformation, i.e., mathematically transforming the
data to a different scale and using that transformed ‘normal’ distribution to derive the 95% confidence
interval. For example, in the case where the data distribution is skewed and the uncertainty is > 100% of
the mean (i.e., where the lower limit would be less than zero), the data could be transformed to a
lognormal distribution. This approach, however, requires the confidence interval to be transformed back
to the original scale to express the uncertainty in the original units, which can introduce error. As a result,
there is a trend away from using transformational approaches due to issues in transforming the data back
to their original scale.
The third assumption states that there is no significant covariance between the uncertainties that are to be
combined, which is equivalent to saying that the errors or uncertainties are independent or that there is no
correlation between the uncertainty terms. The uncertainties in two quantities would be considered
independent if they were estimated by entirely separate processes and there was no common source of
uncertainty. The uncertainties in two quantities would be dependent if they had a common source of
uncertainty (Williamson, 1996). Covariance between two uncertainty terms can be addressed through an
additional term in the uncertainty propagation equations (discussed further in Section 4.2.2). However,
the IPCC Good Practices document suggests avoiding the need for the covariance term in the equation by
“…stratifying the data or combining the categories where the covariance occurs” (IPCC, 2001).
Consider two quantities that can be measured: X and Y. The uncertainty for these values can be
expressed on an absolute basis as ±Ux and ±Uy, respectively, where U is calculated through statistical
analysis (as represented by Equation 4-1), determined through the Monte Carlo technique, or assigned by
expert judgment. Uncertainty may also be expressed on a relative basis, generally reported as percentage:
⎛U ⎞ ⎛U ⎞
± 100⎜ X ⎟ % or ± 100⎜ Y ⎟ %, respectively.
⎝ X ⎠ ⎝ Y ⎠
Depending on the uncertainty propagation equation, the absolute or relative uncertainty value may be
required. In addition, selection of the propagation equation also depends on whether the uncertainties
associated with the individual parameters are independent or correlated.
Two potential equations are used for computing the total uncertainty from the addition or
subtraction of two or more measured quantities. The selection between the two equations
depends on whether the uncertainties associated with the measured quantities, X and Y, are
correlated.
For uncertainties that are mutually independent, or uncorrelated (i.e., the uncertainty terms are not
related to each other), the aggregated uncertainty is calculated as the “square root of the sum of
the squares” using the absolute uncertainties, as shown in Equation 4-4.
where
U(abs) = the absolute uncertainty.
The absolute uncertainty values are used in the equations, and the resulting aggregated
uncertainty (UX+Y+…+N) is also on an absolute basis. Note that where a constant is also included in
the emission estimation calculation, the absolute uncertainty should include the constant. This is
demonstrated in the example provided in Section 4.4.1.
For two uncertainty parameters that are related to each other, the equation becomes:
However, the IPCC Good Practices guidance states, “Once the summation exceeds two terms and
the covariance occurs, the use of the Monte Carlo approach is preferable where data resources are
available” (IPCC, 2001).
The equation for propagating uncertainties from the product or quotient of two or more measured
and independent quantities is similar to Equation 4-4. However, in this case the relative
uncertainties are used, as shown in Equation 4-6. When multiplied by 100, the resulting
combined uncertainty (U(Rel)XxYxN) is expressed as a percentage.
2 2 2
⎛U ⎞ ⎛U ⎞ ⎛U ⎞
U (rel ) X ×Y ×...× N = U (rel ) X ÷Y ÷...÷ N = ⎜ X ⎟ + ⎜ Y ⎟ + ... + ⎜ N ⎟ (Equation 4-6)
⎝ X ⎠ ⎝ Y ⎠ ⎝ N ⎠
Equation 4-7 is used to estimate the uncertainty of a product or quotient of two parameters (X and
Y) where the uncertainties are correlated and positive values. Here also, relative uncertainty
values are used in the equation and the resulting combined uncertainty is on a relative basis.
2 2
⎛U ⎞ ⎛U ⎞ ⎛U U ⎞
U (rel )Correlated X ×Y = ⎜ X ⎟ + ⎜ Y ⎟ + 2r ⎜ X × Y ⎟ (Equation 4-7)
⎝ X ⎠ ⎝ Y ⎠ ⎝ X Y ⎠
c. Combining Uncertainties
If one suspects that the uncertainty parameters are correlated, but data are not available to plot or calculate
the correlation coefficient, the following rule-of-thumb values could be applied using expert judgment
(Franzblau, 1958):
Uncertainty Assessment:
Measurement Uncertainty:
o For this example, the activity data are based on a single point measurement. All of the devices were
accounted for, so there is no bias and the uncertainty of the activity value (i.e., the number of
pneumatic devices) is 0.
o The composition measurements are based on multiple measurements representing a sampling of the
composition. Equations 4-1, 4-2, and 4-3 are applied to calculate the standard deviation of the
composition samples collected throughout the year. The standard deviation accounts for the
measurement error and natural variability of the sampled values.
o For this example, the precision uncertainty of the composition data is assumed to be ± 1%. (A more
detailed demonstration of quantifying uncertainty for measured composition data is provided in a
separate example.) On an absolute basis, this equates to 0.8 for the CH4 composition, and 0.05 for
the CO2 composition.
o Bias associated with the sampling and analysis is assumed to be small due to equipment
maintenance and calibration practices. A value of 3 % is assigned for this assessment. On an
absolute basis, this equates to 2.40 for the CH4 composition, and 0.150 for the CO2 composition.
2
⎛ 2.53 ⎞
U (rel )CH 4 = 0 +⎜2
⎟ + 0.331 = 0.333 = 33.3%
2
⎝ 80 ⎠
2
⎛ 0.158 ⎞
U (rel )CO2 = 0 +⎜
2
⎟ + 0.331 = 0.333 = 33.3%
2
⎝ 5 ⎠
Next, the CH4 emissions are converted to a CO2 equivalent basis, since there are no additional emission
sources to sum for this example.
CO 2 e = ( 226 tonnes CH 4 /yr × 21) + 38.8 tonnes CO2 /yr
= 4, 740 + 38.8 = 4,780 tonnes CO 2 e/yr
Note that the global warming potential (GWP) of CH4 is treated as a constant.
The absolute uncertainty for the aggregated CO2 equivalent emissions are based on Equation 4-5 since the
uncertainties for both CH4 and CO2 are perfectly correlated (r = 1).
U (abs )Correlated X +Y = U X2 + U Y2 + 2r (U X × U Y )
Emissions Factor: Table 4-2 provides emission factors and their uncertainties for each of the components in
Table 4-1. Here also, bias is assumed to be eliminated due to a well designed sampling plan.
Table 4-2. Emission Factors and Uncertainties
Weight Fractions of Emissions: Table 4-3 contains the CH4 weight fractions for gas production and light
crude. The weight fractions are from Table C-7 of the 2009 API Compendium. Expert judgment was used
to estimate the uncertainty. Since the uncertainty of the components and the emission factor give lower and
upper uncertainty bounds, it is important to have lower and upper uncertainty levels for the weight fraction
as well.
Lower Upper
Uncertainty Uncertainty
Operations Wt. Fraction (%) (%)
Conventional Oil: 0.523 5% 5%
Gas Service
Conventional Oil: 0.0363 5% 5%
Light Liquid Service
Source: Picard, D. J., B. D. Ross, and D. W. H. Koon. A Detailed Inventory of CH4
and VOC Emissions from Upstream Oil and Gas Operations in Alberta, Volume II
Development of the Inventory, Canadian Petroleum Association, March 1992,
Tables 12 through 15.
Methane Emissions
Uncertainty Estimate:
o To estimate the uncertainty of the emissions for the individual components, use Equation 4-6, the
equation for the uncertainty of a product. This applies the relative uncertainties for uncertainties
that are independent. For this example, the equation is written as:
For example, the total methane emissions for compressor seals are:
0.713 kg 8760hr tonnes
Emissions = 2 seals × × 0.523 wt. fraction × × = 6.53 tonnes/yr
hr-source year 1000kg
o Note that in the emissions formula, there are two constants. Multiplication by a constant does not
change the relative uncertainty of an estimate. The total methane emissions are the sum of the
methane emissions for all of the components. The total uncertainty for methane emissions resulting
from fugitive sources associated with the electric reciprocating compressor are calculated using
Equation 4-4, which applies the absolute uncertainties for the uncertainty of a sum.
o Table 4-4 gives the methane emissions for the individual components and their uncertainties along
with the total methane emissions and its uncertainty. The 95% confidence interval bounds for the
total methane emissions is shown as (7.41, 30.6) tonnes/year.
Table 4-4. Estimated Fugitive CH4 Emission
Equipment Type Electric reciprocating compressor
Relative Absolute
Lower Upper
Component CH4 Uncertainty Uncertainty Lower Upper
Type Service Emissions (%) (%) Uncertainty Uncertainty
Connectors Gas/Vapor 0.889 36.2 36.2 0.322 0.322
Connectors Light Liquid 0.000350 100a 135 0.000350 0.000474
Controllers Fuel Gas 0.545 41.4 327 0.226 1.78
Compressor Gas/Vapor 6.53 83.3 83.3 5.44 5.44
Seals
Open-ended Gas/Vapor 7.83 85.0 171 6.66 13.4
Lines
Valves Gas/Vapor 0.225 22.5 22.5 0.0507 0.0507
Valves Light Liquid 0.00112 79.5 79.5 0.000889 0.000889
TOTAL 16.0 53.7 91.0 8.61 14.6
a
This value was truncated (set to -100%) since the emissions cannot be less than 0.
Expert judgment involves a person, or group of people, familiar with the systems assigning an
uncertainty to the estimate based on knowledge of the process. It is used when there is no
information to quantify the uncertainty based on data or manufacturer’s specifications of the
measured parameter. Expert judgment is covered in Sections 2.2, [Link], and Annex 2A.1 of the
2006 IPCC Guidelines for National Greenhouse Gas Inventories (IPCC, 2006). Section 7 of the
ISO-5168 explains how to deal with uncertainty estimates based on expert judgment (ISO,
2005).
Input Data:
o The CO2 emissions factor for natural gas in production (non-pipeline quality) is 0.0547 tonnes/106
Btu (HHV) (from Table 4-2 of the 2009 API Compendium). An estimate of uncertainty for this
value is not provided in the original reference, so an uncertainty of 10% at the 95% confidence
interval is assumed based on expert judgment. Bias is assumed to be accounted for in this value.
o Similarly, the natural gas heating values is 1020 Btu/scf (from Table 3-8 of the 2009 API
Compendium). An estimate of uncertainty for this value is not provided in the original reference so
an uncertainty of 10% at the 95% confidence interval is assumed based on expert judgment. Bias is
assumed to be accounted for in this value.
Emission Factor Estimate:
o The emission factor is quantified as the product of the carbon content and heating value, as shown
by the following formula:
Emission Factor Gas Carbon Content Heating Value 106 scf/MMscf
× × × 6
(tonnes CO2 /MMscf) (tonnes CO2 /MMBtu) (Btu/scf) 10 Btu/MMBtu
Uncertainty Estimate:
o For the purpose of this example, the uncertainties associated with the carbon content and heating
value are assumed to be independent because the values are cited from different literature
references and are not based on measured data. Therefore, Equation 4-6 is applied to estimate the
uncertainty for the emission factor. This applies the relative uncertainties for the carbon content
and heating value.
2 2
⎛U ⎞ ⎛U ⎞
U (Rel) X ×Y = ⎜ X ⎟ +⎜ Y ⎟
⎝ X ⎠ ⎝ Y ⎠
Default emission factor Do the default emission Yes Apply reported uncertainty and
factors have uncertainty progress to Uncertainty
specified or quantified? Aggregation
No
Assign uncertainty based on expert
judgment and document reasons
supporting the assignment.
Progress to Uncertainty Aggregation.
Although there is only a single measurement value in this application, such as the cumulative flow rate
through a totalizer meter, there may be more than one parameter that contributes to the uncertainty of the
measurement. For example, Section 3.2 indicates that the following items need to be considered when
estimating measurement uncertainty:
The aggregated uncertainty from these parameters is calculated by either applying Equation 4-4 for
independent parameters or Equation 4-5 for correlated parameters.
The following example demonstrates the uncertainty calculation for estimating CO2 emissions resulting
from the combustion of produced natural gas. The scenario presented is based on a single measurement of
the flow.
Do the measurements Apply Equations 4-1 through 4-3 to calculate the mean and standard
Multiple points represent a sampling of the Yes deviation, respectively.
measured parameter? Apply Equation 4-8 to estimate the uncertainty based on the measured
data, OR apply Equation 4-9 to estimate uncertainty for a single
observed estimate where other data are used to quantify the
No uncertainty.
Determine what (additional) parameters Apply Equation 4-4 to aggregate uncertainty for
contribute to uncertainty. Note that bias Yes addition/ subtraction of data.
may be one of the parameters. Apply Equation 4-6 to aggregate uncertainty for
Are the errors in the measurements multiplication/ division of data.
independent? (See text for guidance)
Input Data
Flow for a gaseous fuel is measured using a totalizer meter prior to being routed to a combustion device.
The total annual flow rate for the fuel is 18,361 MMscf/yr. Because the meter records a cumulative flow,
the uncertainty associated with this measurement is not reduced by taking daily or monthly readings of the
flow. The uncertainty of this measurement is determined by following the decision tree from Figure 4-4
through the measurement uncertainty. The uncertainty of the measurement value is then aggregated.
Uncertainty Assessment:
Measurement uncertainty
o From Section 3.2 Table 3-4, the random error that is expected for a properly installed and operated
orifice meter is 1.5% (95% confidence interval) when the meter is operating at 30-100% of the
measurement range. For this example, this equates to an absolute uncertainty value of:
MMscf
18,361 × 0.015 = 275
yr
o “Context-specific” factors and the uncertainty of pressure and temperature corrections for gas
meters can be based on expert data as summarized in industry standards referenced in Section 3.0.
In the absence of quantitative information about this uncertainty, this example provides a sensitivity
analysis of the impact of assuming a range of context-specific uncertainty values: 10%, 25%, and
100% (on a relative uncertainty basis). These values were selected to test the sensitivity of the
calculations and assess the effect over a wide range of variability.
o For this example, it is assumed that it is not possible to quantify the bias. This meter was installed
according to the manufacture’s requirements 17 years ago and has an expected life span of 25 years.
The meter was last calibrated 7 years ago. Due to the installation and calibration of the equipment,
expert judgment was used to estimate a 5% bias in the measurement. This equates to an absolute
uncertainty value of:
MMscf
18,361 × 0.05 = 918
yr
o As shown in the measurement uncertainty decision tree (Figure 4-4), the selection of an equation
depends on whether the uncertainties are related (correlated). For this example, the uncertainty
components (bias, meter error, and context-specific factors) are independent. To calculate the total
uncertainty of a single measurement for this example, Equation 4-4 is applied to account for the
uncertainties associated with the measurement, the other context-specific factors, and bias:
For example, if the context specific uncertainty is 10% on a relative basis (1,840 on an absolute basis), the
aggregated uncertainty is:
The following table shows the calculated aggregate uncertainty for the different values of the context-specific
factors due to pressure and temperature corrections.
Table 4-5. Uncertainty in a Single Measurement
Single
Meter Measurement Assigned Context Specific Bias, Single Measurement Measurement
Uncertainty, 106 scf Context Specific Uncertainty, 106 scf Total Uncertainty, Total
(1.5%) Uncertainty, % 106 scf (absolute) (5 %) 106 scf (absolute) Uncertainty, %
275 10 1,840 918 2,070 11.3%
275 25 4,590 918 4,690 25.5%
275 100 18,400 918 18,400 100%
The next example demonstrates the uncertainty calculation for estimating CO2 emissions resulting from the
combustion of produced natural gas. Two methods are compared. The first calculates an annual average
CO2 emission factor based on monthly natural gas composition measurements and an annual summation of
daily flow rates. The second determines a monthly CO2 emission factor based on a monthly natural gas
composition sample and applies a monthly summation of daily flow rates.
Table 4-5 shows flow measurement data for produced natural gas that is routed to a combustion device. The fuel
flow is measured continuously (scf/sec) with an orifice meter. A data acquisition system records the daily flow rate
measurements. The daily and monthly totals are shown in Table 4-6.
Uncertainty Assessment:
Measurement uncertainty
o The total volume of gas combusted for this example results in an annual flow rate of 1,187×106 scf.
Because the flow is measured continuously (with the total recorded daily), this is not considered a sampling
of data so Equations 4-1 and 4-3 do not apply.
o From Section 3.2 Table 3-4, the random error that is expected for a properly installed and operated orifice
meter is 1.5% (95% confidence interval) when the meter is operating at 30-100% of the measurement range.
o The context-specific factors discussed previously also apply for this example. Again, a sensitivity analysis
is applied for the context-specific uncertainty values of 10%, 25%, and 100%.
o As with the Single Flow Measurement example, expert judgment was used to assign a 5% bias in the
measurement.
Day Jan. Feb. Mar. April May June July Aug. Sept. Oct. Nov. Dec.
1 3,374 2,344 3,102 3,483 3,321 3,530 2,974 3,484 3,559 3,269 3,239 3,344
2 3,403 2,373 3,329 3,471 3,331 3,518 3,033 3,400 3,494 3,365 3,210 3,231
3 3,381 2,235 3,406 3,530 3,342 3,480 3,254 3,372 3,503 3,236 3,207 3,275
4 3,319 2,276 3,419 3,606 3,335 3,460 3,417 3,406 3,679 3,162 3,132 3,402
5 3,299 2,319 3,555 3,597 3,342 3,503 3,329 3,385 3,636 3,240 3,107 3,570
6 3,173 2,215 3,737 3,501 3,342 3,487 3,034 3,393 3,604 3,233 3,047 3,364
7 3,182 2,459 3,711 3,486 3,378 3,463 2,915 3,449 3,601 3,239 3,060 3,209
8 3,257 2,288 3,771 3,441 3,253 3,323 3,007 3,370 3,541 3,370 3,101 3,344
9 2,985 2,572 3,806 3,400 3,320 3,253 3,158 3,506 3,544 3,347 2,969 3,286
10 3,022 2,433 3,507 3,348 3,301 3,551 3,330 3,412 3,507 3,433 2,901 3,158
11 2,830 2,769 3,340 3,305 3,332 3,478 3,198 3,464 3,557 3,458 3,154 3,186
12 2,914 2,640 3,328 3,370 3,395 3,510 3,262 3,449 3,459 3,405 3,425 3,037
13 2,864 2,431 3,377 3,452 3,348 3,471 3,263 3,341 3,518 3,424 3,488 2,834
14 2,641 2,767 3,432 3,467 3,246 3,302 3,339 3,435 3,484 3,389 3,323 2,849
15 2,647 2,687 3,390 3,465 3,229 3,534 3,158 3,489 3,418 3,397 3,234 2,819
16 2,519 2,888 3,424 3,579 3,322 3,590 3,188 3,468 3,404 3,404 3,434 2,901
17 2,670 2,912 3,432 3,605 3,344 3,475 3,239 3,369 3,425 3,404 3,381 2,807
18 2,482 3,006 3,424 3,543 3,387 3,445 3,114 3,409 3,417 3,407 3,376 2,892
19 2,569 3,089 3,522 3,477 3,356 3,520 3,225 3,362 3,468 3,216 3,382 2,942
20 2,320 3,299 3,444 3,335 3,136 3,561 3,220 3,356 3,474 3,233 3,389 2,915
21 2,405 3,401 3,551 3,545 3,291 3,449 3,122 3,488 3,483 3,295 3,420 2,968
22 2,368 2,945 3,436 3,486 3,463 3,389 2,830 3,374 3,466 3,200 3,464 2,982
23 2,287 3,147 3,495 3,460 3,430 3,474 2,962 3,296 3,447 3,257 3,463 2,942
24 2,277 3,398 3,446 3,461 3,405 3,351 3,546 3,330 3,355 3,282 3,410 3,028
25 2,214 3,521 3,563 3,530 3,430 3,352 3,333 3,422 3,371 3,241 3,346 3,089
26 2,168 3,414 3,511 3,424 3,472 3,316 3,373 3,308 3,383 3,188 3,434 2,963
27 2,216 3,241 3,392 3,482 3,555 3,269 3,516 3,302 3,384 3,177 3,480 2,891
28 2,306 3,206 3,246 3,295 3,550 3,262 3,564 3,483 3,369 3,222 3,307 2,954
29 2,091 3,394 3,336 3,543 3,297 3,547 3,553 3,376 3,300 3,244 3,113
30 2,108 3,489 3,345 3,555 3,262 3,564 3,593 3,304 3,240 3,305 3,047
31 2,209 3,470 3,537 3,611 3,557 3,241 2,972
Total,
103 scf/ 83,500 78,275 107,449 103,825 104,590 102,874 100,624 106,027 104,232 102,277 98,432 95,313
month
o Since the uncertainties associated with the meter, the context-specific factors, and the bias are independent (there
is no autocorrelation), we can apply Equation 4-4 to estimate the uncertainty. The absolute uncertainties are used
in this equation. The results are shown in Table 4-7 below.
∑ (U )
December
U (abs )Year = 2
Measurementi + U context
2
specifici + U Biasi
2
i = January
As with the annual volume, the uncertainties associated with the meter (1.5%), the context-specific factors
(assigned 10%, 25%, and 100% to examine the sensitivity of this uncertainty), and bias (5%) apply to the
daily readings and monthly totals. Summing the daily flow rates for the monthly totals results in the same
relative uncertainties as shown in Table 4-6. Table 4-8 summarizes the absolute uncertainties for the monthly
flow rates.
Table 4-8. Uncertainty in Summation of Flow Measurements – Monthly
Monthly Measurement Flow Rate Uncertainty
Relative Uncertainty, % 11.28% 25.54% 100.14%
10% 25% 100%
Month MMscf (absolute uncertainties, 103 scf)
January 83.50 9,419 21,325 83,614
February 78.28 8,830 19,991 78,382
March 107.45 12,121 27,442 107,595
April 103.83 11,712 26,516 103,966
May 104.59 11,798 26,712 104,733
June 102.87 11,605 26,273 103,014
July 100.62 11,351 25,699 100,761
August 106.03 11,960 27,078 106,171
September 104.23 11,758 26,620 104,374
October 102.28 11,537 26,121 102,417
November 98.43 11,104 25,139 98,566
December 95.31 10,752 24,342 95,443
Do the measurements Apply Equations 4-1 through 4-3 to calculate the mean and
Multiple points Yes standard deviation, respectively.
represent a sampling of the
measured parameter? Apply Equation 4-8 to estimate the uncertainty based on the
measured data, OR apply Equation 4-9 to estimate uncertainty for a
single observed estimate where other data is used for to quantify
the uncertainty.
The sample standard deviation “includes contributions to the precision both from the measurement system
and from the material composition variation from sample to sample” (Coleman and Steele, 1989). In other
words, if the measured data for an emission source are derived from a statistical sample, the use of the
standard deviation as a measure of the spread of the data accounts for uncertainty of the measurement
instrument and the differences among the samples. For example, if three samples are taken every month
and the data for the 36 samples are used to calculate a yearly mean, the uncertainty calculated from these
samples accounts for both the uncertainty in the measurement instrument and the variability among
observations. Thus, this uncertainty will be larger than the uncertainty due to measurement error alone. As
the sample size increases, the uncertainty that combines the instrument error and context-specific factors
will decrease.
Next, we would calculate the standard deviation of the sample using Equation 4-1. The IPCC Good
Practices document and others recommend applying the Equation 4-8 to quantify the uncertainty of the
data set (IPCC, 2001), which is also referred to as the relative expanded uncertainty and can be thought of
as half the bandwidth of a 95% confidence interval (ISO, 2005).
t × s ( x) / n
U (rel )( x) = × 100% (Equation 4-8)
x
where
U(Rel)(x) = the relative expanded uncertainty in the data set (%);
t = a value based on Student’s t-distribution with n-1 degrees of freedom which gives a 95%
confidence interval;
x = is the mean for the set of data calculated in Equation 4-3;
s( x ) = the standard deviation of the data set calculated in Equation 4-2; and
n = the sample size for the set of data.
If one uses a single observation as an estimate and uses other data to calculate the uncertainty,
Section A1.2.3 of the IPCC Good Practices document explains that Equation 4-9 should be used instead of
Equation 4-8 (IPCC, 2006). IPCC’s example application of this approach is for the use of single emission
estimate for a particular year that has been calculated on more than one occasion. The recalculations have
t × s( x)
U ( rel )( x ) = × 100% (Equation 4-9)
x
where
U(Rel)(x) = the relative expanded uncertainty in the data set (%);
t = a value based on Student’s t-distribution with n-1 degrees of freedom which gives a 95%
confidence interval;
x = is the mean for the set of data calculated in Equation 4-3; and
s( x ) = the standard deviation of the data set calculated in Equation 4-2.
Continuing with the decision tree from Figure 4-4, uncertainty from other parameters may also require
consideration, despite the fact that the standard deviation of a statistical sample of observations accounts
for the measurement uncertainty and variability.
Do the measurements Apply Equations 4-1 through 4-3 to calculate the mean and
Multiple points Yes standard deviation, respectively.
represent a sampling of
the measured parameter? Apply Equation 4-8 to estimate the uncertainty based on the
measured data, OR apply Equation 4-9 to estimate uncertainty
for a single observed estimate where other data is used for to
No quantify the uncertainty.
Where multiple uncertainty parameters apply, selection of the equation to aggregate uncertainty parameters
depends on whether the uncertainties are independent. In addition, the activity value may be based on the
product or sum of multiple data or variables. The decision diagram references the appropriate equation
based on these criteria.
For this example, natural gas samples are collected to calculate the CO2 emission factor associated with
combusting the gas. (Note: this example would also apply to composition data for a flared gas stream.)
Although this example is examining uncertainty associated with an emission factor, the decision tree
provided in Figure 4-3 references the figure above (derived from Figure 4-4) for quantifying emission
factor uncertainty where the emission factor is based on multiple measurements from a statistical sample.
The molecular weight of the average composition is calculated by applying the following equation:
1 # compounds
MWMixture = × ∑ (Mole%i × MWi )
100 i =1
Mole% dry January February March April May June July August September October November December
Methane 90.73 91.3 91.37 91.21 94.86 93.97 93.91 94.3 94.45 94.33 94.08 92.85
Ethane 3.83 3.74 4.1 4.16 2.27 2.7 2.81 2.59 2.27 2.3 2.3 2.88
CO2 3.69 3.24 2.93 2.96 1.32 1.23 0.79 1.13 1.38 1.38 1.4 1.49
Propane 0.9 0.8 0.82 0.88 0.61 0.79 0.66 0.65 0.65 0.65 0.66 0.76
i-Butane 0.11 0.1 0.1 0.11 0.09 0.12 0.09 0.1 0.1 0.1 0.1 0.11
n-Butane 0.14 0.13 0.13 0.13 0.12 0.16 0.15 0.14 0.13 0.13 0.13 0.14
i-Pentane 0.03 0.04 0.03 0.04 0.04 0.05 0.04 0.05 0.04 0.04 0.04 0.05
n-Pentane 0.02 0.03 0.02 0.03 0.03 0.04 0.03 0.04 0.03 0.03 0.03 0.04
C6+ 0.03 0.05 0.04 0.05 0.05 0.07 0.06 0.24 0.05 0.05 0.04 0.08
Total 99.48 99.43 99.54 99.57 99.39 99.13 98.54 99.24 99.1 99.01 98.78 98.4
Mean, MW wt% C
mole% U(abs), Calculation, U(abs), calculation, U(abs),
dry Std dev t*s/sqrt(n) U(rel) mole % lb/lbmole lb/lbmole wt% wt%
Methane 93.11333 1.529369 0.971714736 1.04% 0.971715 14.93538 0.1559 64.7728 1.4797
Ethane 2.995833 0.747498 0.474937795 15.85% 0.474938 0.900847 0.1428 4.168007 0.6662
CO2 1.911667 0.988652 0.62815969 32.86% 0.62816 0.841325 0.2765 1.32982 0.4378
Propane 0.735833 0.100856 0.064080962 8.71% 0.064081 0.324503 0.0283 1.535612 0.1373
i-Butane 0.1025 0.00866 0.005502463 5.37% 0.005502 0.059573 0.0032 0.28521 0.0164
n-Butane 0.135833 0.010836 0.006885023 5.07% 0.006885 0.078946 0.0040 0.377961 0.0206
i-Pentane 0.040833 0.006686 0.004247814 10.40% 0.004248 0.029461 0.0031 0.142025 0.0151
n-Pentane 0.030833 0.006686 0.004247814 13.78% 0.004248 0.022246 0.0031 0.107244 0.0149
C6+ 0.0675 0.055942 0.035544067 52.66% 0.035544 0.058172 0.0306 0.281732 0.1485
Total 99.13417 17.25045 73.00041
Sum uncertainty (abs) 0.350567 1.693261
Sum uncertainty (rel) 2.03% 2.32%
For sample size of n=12, the Student’s t distribution value is 2.201.
Note, the resulting estimates are not rounded to show the comparison.
o The uncertainty associated with this value is calculated by applying Equation 4-4, using the absolute
uncertainties for each molar compound.
U (abs )
∑ ( mole %×MW )
= ∑U (abs) 2
mole %× MW
The weight percent carbon of the average composition is calculated by applying the following equation:
lbmolei x lbmole C 12 lb C 1
Wt% CTotal = ∑ × × ×
100 lbmole total lbmolei lbmole C MWTotal
This equates to 73.00041% C, as shown in Table 4-9. The uncertainty associated with this value is calculated
by applying Equation 4-1, using the absolute uncertainties for the weight percent carbon of each molar
compound.
o The uncertainty associated with this value is calculated by applying Equation 4-1, using the absolute
uncertainties for each molar compound:
U (abs )
∑ ( wt %)
= ∑U (abs) 2
wt %
o The facility collected one natural gas composition sample each month. For comparison with the annual average
gas composition example above, the same composition values shown in Table 4-8 will be used.
o For this example, only one data point is available each month. Table C-2 (Appendix C) provides reproducibility
uncertainty associated with natural gas samples. These values can be applied to account for the measurement
error in each monthly sample. Note that the reproducibility uncertainty varies based on the mole % of each gas
compound. The values are shown for January in Table 4-10.
o An additional 5% uncertainty is assigned by expert judgment to account for potential variability and bias in the
gas composition during the month. The combined uncertainty is calculated by applying Equation 4-4, using the
absolute uncertainties. The calculation is demonstrated for CH4. Results for January are shown in Table 4-11.
Combined wt%
Reproducibility Variability mole% MW Carbon
Uncertainty Uncertainty Uncertainty Calculation, Calculation, wt% C
mole% dry January (abs), mole% (Uabs=mole%x5%) (rel) lb/lbmole % U(rel)
Methane 90.73 0.15 4.5365 5.00% 14.5531 60.71 6.49%
Ethane 3.83 0.1 0.1915 5.64% 1.1517 5.13 6.99%
CO2 3.69 0.1 0.1845 5.69% 1.6240 2.47 7.03%
Propane 0.9 0.07 0.045 9.25% 0.3969 1.81 10.13%
i-Butane 0.11 0.07 0.0055 63.83% 0.0639 0.29 63.97%
n-Butane 0.14 0.07 0.007 50.25% 0.0814 0.37 50.42%
i-Pentane 0.03 0.02 0.0015 66.85% 0.0216 0.10 66.98%
n-Pentane 0.02 0.02 0.001 100.12% 0.0144 0.07 100.21%
C6+ 0.03 0.02 0.0015 66.85% 0.0259 0.12 66.98%
Total 99.48 17.9329 71.0717
Sum uncertainty (abs) 0.7404 3.9734
Sum uncertainty (rel) 4.13% 5.59%
o As with the annual composition example, the molecular weight of the monthly composition is calculated by
applying the following equation:
1 # compounds
MWMixture = × ∑ ( mole%i × MWi )
100 i =1
o For each individual gas compound, the relative uncertainty of the mole%i × MWi is equivalent to the
combined reproducibility and variability uncertainties since the molecular weight of each gas
compound is a constant. The aggregated uncertainty associated with the mixture’s MW is calculated
by applying Equation 4-1, using the absolute uncertainties for each molar compound.
U (abs )
∑ ( MW Total )
= ∑U (abs) 2
mole %× MW
o The weight percent carbon of the average composition is calculated by applying the following
equation:
lbmolei x lbmole C 12 lb C 1
Wt% CTotal = ∑ × × ×
100 lbmole total lbmolei lbmole C MWTotal
o For each compound in the gas mixture, the uncertainty from this calculation is determined by
applying Equation 4-6, using the relative uncertainties of the lbmolei and MWTotal. This is
demonstrated for the January CH4 composition.
o The uncertainty associated with the wt% C of the mixture is calculated by applying Equation 4-4,
using the absolute uncertainties for the weight percent carbon of each molar compound.
U (abs )
∑ ( wt %)
= ∑U (abs) 2
Wt %
o The monthly CO2 emissions factor is then calculated by applying the following equation:
o The uncertainty for this quantity is calculated by applying Equation 4-3 for the relative uncertainties
of MWTotal and wt% CTotal
2 2
⎛U ⎞ ⎛U ⎞
U (rel ) = ⎜ X ⎟ + ⎜ Y ⎟ = 4.132 + 5.592 =6.95%
⎝ X ⎠ ⎝ Y ⎠
o The same methods are applied to each of the monthly samples. Table 4-12 summarizes the resulting
emission factors (tonnes CO2/MMscf) and uncertainties.
Table 4-12. Measured Monthly Emission Factors
Tonnes U(abs), Tonnes
CO2/MMscf CO2/MMscf U(rel), %
January 55.8856 3.8840 6.95
February 55.7698 3.9276 7.04
March 55.9697 3.9414 7.04
April 56.1644 3.9304 7.00
May 54.7437 4.1879 7.65
June 55.2278 4.1441 7.50
July 54.7069 4.1606 7.61
August 55.4646 4.1751 7.53
September 54.6648 4.1705 7.63
October 54.6332 4.1647 7.62
November 54.4964 4.1527 7.62
December 54.8858 4.0798 7.43
Figure 4-5 provides the decision tree for aggregating uncertainties. This diagram steps through the
combination of activity data and emission factors to result in the final emission estimate. It also addresses
the aggregation of emissions from multiple sources, applying the GWP to convert non-CO2 emissions to a
CO2 equivalent basis, and the final summation of CO2 equivalent emissions.
An important point to note is that, although the GWP values reported by IPCC have corresponding
uncertainties, the GWP values are treated as constants in compiling a facility- or entity-wide GHG
inventory. Multiplying by a constant does not change the relative uncertainty.
Total CO2 emissions are calculated based on the product of the activity data (fuel consumption) and the
emission factor (tonnes CO2/volume fuel). Table 4-13 compares the emission estimate results from three
approaches:
− Annual emissions based on the measured gas consumption and a default emission factor;
− Annual emissions based on the measured gas consumption and an emission factor derived from an
annual average gas composition; and
− Annual emissions based on an aggregate of monthly gas composition and flow rate measurements.
Are the errors for the Independent errors: Apply Equation 4-4
following aggregations If two emission sources have the same to aggregate uncertainty for addition/
correlated? emission factor (for example, boilers subtraction of data.
Refer to the text for and heaters) combine the activity data Correlated errors: Apply Equation 4-5
guidance on testing for and aggregate the uncertainty before to aggregate uncertainty for addition/
autocorrelation and for applying the emission factor. subtraction of data.
methods to make the terms
independent.
Independent errors: Apply Equation 4-6
to aggregate uncertainty for
Aggregate emissions for multiplication/ division of data.
the source as AF x EF. Correlated errors: Apply Equation 4-7
to aggregate uncertainty for
multiplication/ division of data.
Uncertainty (rel)
Context specific factor sensitivity
Tonnes Emissions, analysis
Month CO2/MMscf U(rel) MMscf tonnes CO2 10% 25% 100%
Emissions based on January 55.8856 6.95% 83.50 4,666.44 13.25% 26.47% 100.38%
monthly flow and February 55.7698 7.04% 78.28 4,365.38 13.30% 26.49% 100.38%
carbon content March 55.9697 7.04% 107.45 6,013.89 13.30% 26.49% 100.38%
April 56.1644 7.00% 103.83 5,831.27 13.27% 26.48% 100.38%
May 54.7437 7.65% 104.59 5,725.67 13.63% 26.66% 100.43%
June 55.2278 7.50% 102.87 5,681.48 13.55% 26.62% 100.42%
July 54.7069 7.61% 100.62 5,504.85 13.60% 26.65% 100.42%
August 55.4646 7.53% 106.03 5,880.73 13.56% 26.63% 100.42%
September 54.6648 7.63% 104.23 5,697.83 13.62% 26.65% 100.43%
October 54.6332 7.62% 102.28 5,587.75 13.61% 26.65% 100.43%
November 54.4964 7.62% 98.43 5,364.20 13.61% 26.65% 100.43%
Annual Total 1,187.42 65,550.85 3.91% 7.71% 29.10%
Emissions based on annual flow and annual
average carbon content 55.2180 3.08% 1,187.42 65,566.90 11.69% 25.72% 100.18%
Emissions based on annual flow and default
emission factor 55.7940 14.14% 1,187.42 66,250.86 18.09% 29.19% 101.13%
Note, the resulting emission estimates are not rounded to show the comparison.
Inappropriate rounding of the data can lead to errors in the final estimate. Using computer software, such
as spreadsheets, helps the user avoid rounding during intermediate steps. The estimate of emissions should
be rounded to smallest unit of measure (API, [Link], 1985). The uncertainty should be rounded to the
same number of digits as the estimate.
As shown in Equations 4-5 and 4-7, the uncertainty propagation equation can be extended to account for
uncertainty terms that are correlated or not independent. A simple way to assess if the uncertainties are
correlated is to examine a graph of the uncertainties. If there is no pattern, the uncertainties are most likely
independent. If there is a pattern to the uncertainties, they are not independent.
The correlation between the uncertainties of two measured parameters can be calculated by applying the
following equation:
rU X ,UY =
∑ (U Xi ) (
− U X × U Yi − U Y ) (Equation 4-9)
(n − 1) × s (U X ) × s (U Y )
where
rUxUy = the correlation coefficient of the uncertainties for X and Y;
n= the sample size;
Uxi = the uncertainties associated with sample points from source X;
UX = the mean of the uncertainties from source X;
UY =the mean of the uncertainties from source y;
Uyi = the uncertainties associated with sample points from source Y;
s(UX) = the standard deviation for the uncertainties of source X; and
s(UY) = the standard deviation for the uncertainties of source Y.
Before combining the data it is important to eliminate the correlation of uncertainties, if possible. ISO
5168:2005(E) Annex F discusses methods for making measurement uncertainties independent, such as
calibrating instruments against different references and “redefining mathematical relationships to eliminate
correlations” (ISO, 2005).
Section A1.4.5 of the IPCC Good Practices document lists four sources for correlation:
• Mutual constraints on a group of emission estimates (such as a specified total fuel usage which
provides input to a number of processes);
• External drivers that affect a suite of emissions (economic, climatic, resource based) (IPCC, 2001).
To eliminate uncertainty correlations due to common activity data, IPCC suggests summing the emission
factor estimates and then multiplying the activity factor by the sum to obtain total emissions. For mutual
constraint, the IPCC Good Practices document suggests “to leave one of the proportions unspecified, and
to determine it by the difference between the other proportions and the total fraction.”
In addition, IPCC’s Good Practices document recommends that if there are more than two correlated
variables, a Monte Carlo simulation should be applied instead of uncertainty propagation.
Monte Carlo simulation is a more complex, model-based method for iteratively evaluating uncertainty
associated with individual parameters. It may be the preferred approach where more complex equations are
assessed. It is one of many methods for analyzing uncertainty propagation where the goal is to determine
how random variation, lack of knowledge, or error affects the sensitivity, performance, or reliability of the
resulting emission inventory.
Monte Carlo simulation is a “repeated sampling and calculation method” because the inputs are randomly
generated from probability distributions of the respective variables to simulate the process of sampling
from an actual population. Inputs are defined as probabilistic parameters rather than simple estimates. The
uncertainty model relies on repeated random sampling of all inputs and simultaneous recalculation of
emissions (outputs) to measure variation over the course of numerous model iterations. The data generated
from the Monte Carlo simulation can be represented as probability distributions (or histograms) or
converted to error bars and confidence intervals.
Monte Carlo simulation has an advantage over uncertainty propagation in that one can specify multivariate
distributions to account for correlations between different sources of uncertainty. The IPCC Good
Practices document recommends choosing one of the following distributions: normal, lognormal, uniform,
or triangular.1 The major difficulty in using the Monte Carlo simulation is the need to determine the
distribution of the data. If there are not enough data to assume normality by the Central Limit Theorem
(more than 30 data points), there are most likely not enough data to determine the underlying distribution of
the data. Consequently, such analyses are often forced to rely on subject matter expert opinion to
determine distributions rather than on observed data.
1 The IPCC Good Practice document discusses how to perform Monte Carlo simulation in Section 6.2 (IPCC, 2006).
It discusses choice of distribution in Section A1.2.5.
To perform a Monte Carlo simulation, one must first determine the distribution of the data for each of the
uncertain variables used in the emissions model estimate. Ideally, each of these distributions should be
derived from data and knowledge of the underlying process. It is helpful in many instances to first graph
the data, and using the shape of the graph to determine the underlying distributions. The parameters that
define the distribution could be derived from the data. For example, a normal distribution is defined by its
mean and variance. If data are limited, one may have to rely on expert judgment to determine the
underlying distribution.
It is then necessary to statistically test the hypothesis that the data follow a certain distribution. The test
will vary based on the hypothesized distribution. For example, the Shapiro-Wilks test is often used to test
if the data are normal or lognormal. Options to test for other distributions include Empirical Distribution
Functions (Singh and Singh, 2006).
Once distributions are determined for all of the data sources, the Monte Carlo simulation will proceed by
randomly sampling each of the distributions that describe the data used for estimating emissions. As many
as 10,000 replicate samples are typically taken, with the total emissions being estimated for each replicate.
These repeated determinations of emission are used to generate a distribution of the total emissions with its
mean being the estimate of total emissions, and its uncertainty determined by its variance.
Section 6.3.1 of the IPCC Good Practices document compares the uncertainty propagation method and the
Monte Carlo simulations (IPCC, 2006). It notes that the uncertainty propagation method’s assumption of
normality leads to symmetric 95% confidence intervals whereas the Monte Carlo method can take into
account the fact that emissions are bounded below by zero to fit an asymmetric (and thus narrower)
confidence interval. If the data are skewed and one transforms the data (discussed earlier), one could
achieve the asymmetric confidence intervals using uncertainty propagation, as well.
Since the Monte Carlo simulations can assume a truncated distribution, the lower confidence limits tend to
be closer to the mean than the upper confidence limits. The IPCC Good Practices document goes on to
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.
5.2.1 Background
Tables 5-1 and 5-2 summarize the emission sources and uncertainty values associated with this facility.
This example field is described in detail in the API Compendium, but is summarized here for context. A
discussion of the uncertainty values associated with the example facility parameters follows.
Facility Description: An onshore oil field in Texas consists of 320 producing oil wells.
Throughput: The average daily oil and gas production rates are 6,100 bbl/day and 30×106 scf/day,
respectively.
Operations: The facility operates approximately 343 days per year. The facility imports 917 MW-hr
annually from the eGRID subregion “ERCOT all”. The facility gas composition is presented in Table 5-3
and results in a heating value of 928 Btu/scf with an uncertainty of ± 4%, based on engineering judgment.
Table 5-3. Gas Composition for Onshore Oil Field (High CO2 Content)
(Reprinted from the API Compendium)
Most of the GHG emissions from this facility are from combustion sources. A small additional amount
comes from vented sources since the oil field gas contains a relatively high CO2 content (12 mole %).
In this example, emissions from some emission sources are calculated by use of emission factors, using the
following simple equation:
For these emission sources, the first task in the calculation of total inventory uncertainty is assigning
uncertainties to each of the activity factors and emission factors that were used in the inventory. Some
emission calculation equations have multiple terms, rather than only two simple terms. For example:
The same approaches discussed below for assigning uncertainties to simple activity factors will apply to
each term in a more complicated equation. Further detail is also provided in Section 4.0.
Sites may have a combination of directly measured emissions and emissions calculated using activity factors
and emission factors. For direct measurements of emissions, the equations and decision trees provided in
Section 4.0 can be used to quantify uncertainties.
This section discusses assignment of uncertainties to the activity factors. In some cases, measured data will
allow assignment of calculated confidence bounds, as covered in detail in the Section 4.0 examples.
However, there will be many cases where confidence bounds on activity factor data will be assigned by
expert judgment. This section covers how expert judgment assigned some of those bounds in this
hypothetical case.
Some assignment of expert judgment may require that the expert look at the activity value, say a total count
of 67 chemical injection pumps, and decide how many pumps a person might miss (undercount), or how
many devices might be non-operating, or incorrectly identified as chemical injection pumps (over counted).
If the expert, familiar with the count, decided that at most three pumps were under or over counted, the
uncertainty limit might be assigned simply to be 3/67 = 4%. In this case, the company expert decided that a
blanket 5% uncertainty was more appropriate. Some general categories of uncertainty assignments follows.
In this hypothetical example, some of the activity factors were simply counts of large equipment types at the
facility. These counts, such as the number of boilers at the facility, are known with absolute certainty since
they are discreet values and were produced by company personnel intimately familiar with the specific
facility. Therefore, the 95% confidence bounds (shown as “% uncertainty”) of the activity factor value is
zero (0), since the value is perfectly known. Table 5-1 shows uncertainties assumed for each activity factor.
The perfectly known activity factors are in the second column under “No. of Units” and the uncertainties
Counts of other devices on site, even if produced by company personnel, may not be perfectly known. It is
possible to know the exact number of amine units with perfect certainty, but not to know some less
significant equipment components. In this hypothetical case, counts were taken from the exact facility, such
as counts of pressure relief valves (PRVs), pneumatic devices, and chemical injection pumps. However,
some small uncertainty was assigned to these counts by expert judgment since the counts are not perfectly
known. In this hypothetical case, uncertainties between 1% and 5% were assigned to the counts of PRVs,
pneumatic devices, and chemical injection pumps. For unit capacities, the fifth column in Table 5-1, these
capacities are well known, and expert judgment assigned them to be ±5%. In this hypothetical example,
“days of operation” are supplied by company experts. In this example, these estimates were not based on
measured operating hours but were estimated based on operator logs. An uncertainty of 2% was assigned
by expert judgment.
For other sources, such as counts of valves, seals, flanges for fugitive emissions, there may be larger
uncertainties associated with the activity factor counts. This is true even if the count were produced for the
exact facility. The uncertainty is determined using expert judgment based on the quality of the component
inventory and the length of time since the last inventory of fugitive emission sources, since changes to the
facility and therefore component counts are possible. In this hypothetical example, the counts were not
taken for this exact facility, but were applied from average component count information from industry-
wide counts of fugitive components from this general type of facility. Therefore expert judgment assigned a
much wider uncertainty to each count of ±75%. In cases where exact component counts were taken for this
facility, the uncertainty percentage would be much lower. Table 5-2 shows the activity factor counts for
fugitive components (reprinted directly from the API Compendium), and the assigned uncertainties.
In this hypothetical example, some of the activity factors were flow rate information measured by on-site
totalizing meters that record data. These flow rate measurements were assigned a ±15% uncertainty by
Assignment of uncertainty to the most detailed level of activity data used in the inventory is the preferred
approach. Then statistical methods are used to propagate the uncertainty values.
For example, in Table 5-1, the sum of the Annual Activity Factors is calculated as:
Annual Activity Factor = Number of Units × Unit Capacity × Average Operations Time
Each term that is multiplied to result in the annual activity factor has its own uncertainties. For sources
such as the six boilers, there is no (0%) uncertainty associated with the count of equipment since it is
know that there are exactly six boilers at this facility. The uncertainty in the unit capacity is the
uncertainty in the unit’s average operational capacity.
Uncertainties for these combustion sources are calculated by applying Equation 4-6 and using the relative
uncertainty values.
For example, the uncertainty for the activity factor for the heaters/reboilers is calculated using the values
shown in Table 5-1:
Appendix F provides details on the uncertainty estimates for each emission source included in this example.
The resulting uncertainty estimates for this example facility are provided in Table 5-4 and Appendix F.
Source CO2 CH4 N2O and Other GHG Total Emissions, CO2 Eq
Uncertainty Uncertainty Uncertainty Uncertainty
Emissions % Emissions % Emissions % Emissions %
Source Type (tonnes/yr) Lower Upper (tonnes/yr) Lower Upper (tonnes/yr) Lower Upper (tonnes/yr) Lower Upper
Combustion Boiler/heaters 5,200 8.78 8.78 0.0865 26.1 26.1 0.0242 100 150 5,210 8.77 8.77
Sources Natural gas engines 13,900 15.7 15.7 0.904 29.4 29.4 0.325 100 151 14,100 15.6 15.6
Emergency generator IC
0.0108 27.8 27.8 0.00175 100 151
engine 219 15.6 15.6 220 15.5 15.5
Fire water pump IC engine 0.00112 100 106 0.0000467 100 151
Flares 27,400 23.4 23.4 153 25.3 25.3 0.223 100 200 30,700 21.1 21.1
Fleet vehicles 127 19.4 19.4 0.00643 100 151 0.00871 100 151 129 19.1 19.2
Combustion Total 46,800 14.5 14.5 154 25.2 25.2 0.582 68.9 114 49,900 13.7 13.7
Vented sources Dehydration and Kimray
105 77.5 77.5 254 77.5 77.5 NA NA NA 5,440 76.0 76.0
pump vents
Tanks – flashing losses 775 90.4 90.4 1,880 90.4 90.4 NA NA NA 40,300 88.7 88.7
Amine unit 62,600 6.97 6.97 193 100 119 NA NA NA 66,700 8.94 9.77
Pneumatic devices 64.6 50.2 50.2 157 50.2 50.2 NA NA NA 3,360 49.2 49.2
Chemical injection pumps 48.6 100 108 118 100 108 NA NA NA 2,530 98.1 106
Vessel blowdowns 0.0702 100 326 0.171 100 326 NA NA NA 3.65 98.1 319
Compressor starts 0.745 100 190 1.81 100 190 NA NA NA 38.7 98.1 187
Compressor blowdowns 0.333 100 179 0.808 100 179 NA NA NA 17.3 98.1 175
Well workovers 0.0181 100 300 0.0439 100 300 NA NA NA 0.939 98.1 294
Other non-routine (PRVs) 0.131 100 310 0.318 100 310 NA NA NA 6.81 98.1 319
Vented Total 63,600 6.95 6.95 2,610 66.3 66.5 NA NA NA 118,000 30.9 31.0
Source CO2 CH4 N2O and Other GHG Total Emissions, CO2 Eq
Uncertainty Uncertainty Uncertainty Uncertainty
Emissions % Emissions % Emissions % Emissions %
Source Type (tonnes/yr) Lower Upper (tonnes/yr) Lower Upper (tonnes/yr) Lower Upper (tonnes/yr) Lower Upper
Fugitive Sources Fugitive components NA NA NA 52.6 66.2 83.3 NA NA NA 1,100 66.2 83.3
Fleet vehicle
NA NA NA NA NA NA 0.00100 100 112 1.30 100 112
refrigeration, R-314a
Fugitive Total NA NA NA 52.6 66.2 83.3 0.00100 100 112 1,100 66.1 83.2
Indirect Sources Electricity consumed 551 10.2 10.2 0.00776 100 100 0.00628 100 100 553 10.2 10.2
Indirect Total 551 10.2 10.2 0.00776 100 100 0.00628 100 100 553 10.2 10.2
TOTAL (tonnes of each gas) 111,000 7.29 7.29 2,820 61.4 61.7 N2O:0.588 67.1 113
R134a:
0.00100 100 112
TOTAL (CO2e) 111,000 7.29 7.29 59,100 61.4 61.7 184 66.6 112 170,300 21.9 21.9
Appendix F demonstrates the use of an asymmetric uncertainty distribution. Where the uncertainty values for
emissions data were greater than 100%, the lower bound estimate was truncated at -100% and tracked
separately from the upper bound estimate. The upper and lower bound uncertainties were then propagated
through to the total CO2 equivalent emissions for the facility. For this example, the resulting difference
between the upper and lower uncertainty values was very small.
5.3.1 Background
Example 1 provides a very detailed, step-by-step demonstration of the uncertainty calculations for a crude oil
production facility. Because many of the source types are similar across all industry sectors, this example
focuses on a few distinct refinery process units. A comparison in the uncertainty calculations is provided for
methodologies presented in the API Compendium for fluid catalytic cracking unit (FCCU) and hydrogen plant
emissions. As with the previous example, much of the information is reprinted directly from the API
Compendium without change.
The following example applies uncertainty calculations to Exhibit 5.6 from the API Compendium for the
FCCU GHG emission estimation methods presented in the API Compendium. The same operating
parameters specified in the API Compendium are applied here, with the following assignment of uncertainties
added.
• The catalytic cracking unit has a coke burn rate of 119,750 tonnes per year ± 15% and a blower air
capacity of 2,150 m3/min ± 15% (assigned by expert judgment). The air blower is assumed to
operate continuously for the year (a ± 2% uncertainty is applied to this assumption).
• The carbon fraction of the coke is 0.93 ± 5.5% based on site-specific data (determined from
measured compositions).
• The flue gas concentrations are 11% for CO2 and 9% for CO exiting the regenerator. Table D-3 of
this uncertainty document provides reproducibility values for the precision of Reformed Gas
Samples based on ASTM 1946-90. For molar compositions between 5 and 25 percent, a
reproducibility factor of 0.5 applies. An additional 5% uncertainty is assigned by expert judgment
to account for potential variability in the composition.
• It is assumed that no CH4 is formed during the regeneration process.
• A CO boiler is used for control of the flue gas stream. Supplemental firing with natural gas is
employed at a rate of 100×106 ± 5% Btu/hr on a higher heating value basis.
The API Compendium presents three equations for estimating CO2 emissions from FCCUs. The following
demonstrates the uncertainty quantification for each of the three methods.
The uncertainty is calculated by applying Equation 4-6 and using the relative uncertainty values.
Applying the air rate to API Compendium Equation 5-5, the CO2 emission estimate is:
For this calculation, the uncertainty associated with the sum of the CO2 and CO is determined first. This
uncertainty is calculated by applying Equation 4-4, using the absolute uncertainties.
Equation 4-4, using the absolute uncertainties, is also applied to combine these compositions.
The CO2 uncertainty is then calculated by applying Equation 4-6 and using the relative uncertainty values.
From an uncertainty perspective, this calculation is the same as shown for the “K1, K2, K3” Approach (API
Compendium Equation 5-5). Both equations apply the summed composition of CO2 and CO, the air rate,
and an annual operational time. The uncertainty of 15.94%, calculated above, would apply for this
approach as well.
The emissions from the supplemental firing are in addition to the CO2 emissions from the FCCU
regenerator. Emissions from the supplemental firing of natural gas are estimated using the combustion
emission approaches presented in API Compendium Section 4. For CO2, the emission factor is taken from
API Compendium Table 4-3 for pipeline natural gas.
The CO2 uncertainty is calculated by applying Equation 4-6 and using relative uncertainty values. API
Compendium Table 4-3 does not provide an indication of the uncertainty associated with this emission
factor, so a value of 5% is assigned by expert judgment.
The CH4 emission factor is taken from API Compendium Table 4-7 for natural gas fired boilers.
The uncertainty is calculated by applying Equation 4-6 using the relative uncertainty values. The CH4
emission factor has a “B” quality rating assigned to it. Based on this, an uncertainty of 10% is applied to
the emission factor.
The N2O emission factor is also taken from API Compendium Table 4-7 for controlled natural gas fired
boilers.
The uncertainty is calculated by applying Equation 4-6 using the relative uncertainty values. The N2O
emission factor has an “E” quality rating assigned to it. Based on this, an uncertainty of 50% is applied to
the emission factor.
Table 5-5 summarizes the emission and uncertainty estimates for the different FCCU methodologies.
Equation 4-4, using the absolute uncertainty values, is applied where the emission estimates are summed.
Based on the assumptions applied for the FCCU methodologies, the uncertainty associated with each of the
methods provided in the API Compendium is comparable. In all three equations, the aggregated uncertainty is
influenced primarily by the ± 15% uncertainty values assigned to the coke burn rate (used in the first and third
methods) and the blower air capacity (used in the second method).
The API Compendium provides two rigorous calculation approaches for estimating the CO2 generation rate
from the hydrogen plant, both using a specific feed gas composition. The approaches are based on either the
volume of feedstock used or the hydrogen production rate. This section compares the uncertainty estimates
associated with these methods. For this comparison, the following operating parameters and uncertainties are
assigned:
• A hydrogen plant has a feedstock rate of 3×109 ± 15% standard cubic feet per year and produces
13×109 ± 15% standard cubic feet of hydrogen per year.
• The feed gas composition (molar basis) is CH4 = 85%, C2H6 = 8%, C4H10= 3%; the balance is inerts
(assume N2 for the inerts). Table D-2 of this Uncertainty document provides reproducibility
uncertainty associated with natural gas samples. These values can be applied to account for the
measurement error in the composition sample. An additional 5% uncertainty is assigned by expert
judgment to account for potential variability and bias in the gas composition during the month.
• It is assumed that no CH4 is entrained in the hydrogen product.
For both methods, the uncertainty associated with the feedstock composition is needed. The combined
uncertainty of the reproducibility and variability is calculated by applying Equation 4-4, using the absolute
uncertainties. The calculation is demonstrated for ethane below; results for all of the compounds are shown in
Table 5-6.
Combined
Uncertainty wt%
Reproducibility Variability (rel) MW Carbon wt%
mole% Uncertainty Uncertainty (applies to Calculation, Calculation, C
dry (abs), mole% (Uabs=mole%×5%) MW) lb/lbmole % U(rel)
Methane 85 0.15 4.25 5.00% 13.6340 53.96 6.24%
Ethane 8 0.12 0.4 5.22% 2.4056 10.16 6.41%
Butane 3 0.1 0.15 6.01% 1.7436 7.62 7.07%
N2 4 0.1 0.2 5.59% 1.1200 0.00 6.72%
Total 100 18.9032 71.7339
Sum uncertainty
(abs) 0.7042 3.4704
Sum uncertainty
(rel) 3.73% 4.84%
The molecular weight of the gas is calculated by applying the following equation:
1 # compounds
MWMixture = × ∑ (Mole% i × MWi )
100 i =1
For each individual gas compound, the relative uncertainty of the mole%i × MWi is equivalent to the
combined reproducibility and variability uncertainties, since the molecular weight of each gas compound is a
constant. The aggregated uncertainty associated with the mixture’s MW is calculated by applying
Equation 4-4, using the absolute uncertainties for each molar compound.
U (abs )
∑ ( MW Total )
= ∑U (abs) 2
mole %× MW
The weight percent carbon of the average composition is calculated by applying the following equation:
For each compound in the gas mixture, the uncertainty from this calculation is determined by applying
Equation 4-6, using the relative uncertainties of the lbmolei and MWTotal. This is demonstrated for ethane.
U ( R e l )Wt %Ci = U (rel ) 2mole %i × U (rel ) 2MW Total = 5.222 + 3.732 = 6.41%
The uncertainty associated with the wt% C of the mixture is calculated by applying Equation 4-4, using the
absolute uncertainties for the weight percent carbon of each molar compound.
U (abs )
∑ ( wt %)
= ∑U (abs) 2
Wt %
The uncertainties associated with the composition are used in the following comparison of the emission
estimation methodologies for H2 plants.
The CO2 uncertainty is then calculated by applying Equation 4-6 and using the relative uncertainty values.
Combined
#C # H2 Uncertainty (rel)
Compound Atoms Molecules Concentration (mole%) Moles C Moles H2
CH4 1 4 0.85 5.00% 0.85 3.4
C2H6 2 7 0.08 5.22% 0.16 0.56
C4H10 4 13 0.03 6.01% 0.12 0.39
Total Moles 1.13 4.35
Sum uncertainty
0.0439 0.0439
(abs)
Sum uncertainty
3.89% 3.89%
(rel)
The carbon-to-hydrogen ratio is calculated by dividing the total moles C by the total moles H2 (1.13/4.35 =
0.26). Since each mole of carbon produces 1 mole of CO2, the CO2/H2 ratio is the
same as the C/H2 ratio (0.26).
The uncertainty associated with the total moles of C and H2 is calculated by applying Equation 4-4, using the
absolute uncertainties of each molar compound.
U (abs )
∑ ( mole C)
= ∑U (abs) 2
mole C
U (abs )
∑ ( mole H2 )
= ∑U (abs) 2
mole H 2
The CO2 uncertainty is then calculated by applying Equation 4-6 and using the relative uncertainty values.
This section addresses the potential need to refine the emission inventory to reduce the uncertainty in the
overall emission estimate. There may be several reasons to do this. In some specific locations, there may be
state or national regulations, or guidelines for voluntary programs that suggest or require an emission
inventory to have uncertainties lower than a certain percentage level. A company may also independently
wish to refine its own uncertainty limits, if it considers the uncertainties too large.
The reader should note that none of the strategies mentioned here are aimed at reducing the actual emissions
of GHGs. That is a separate subject, and while emission reductions are achievable in some cases, they are
beyond the scope of this report. This section focuses on reducing the mathematical and statistical uncertainty
associated with an existing emission inventory.
Once a decision has been made to refine and reduce the uncertainty associated with a given inventory, some
strategic analysis of the major sources contributing to the uncertainty is in order. It is important to know
which sources significantly contribute to the overall inventory. It would make little sense to refine a term
with large confidence bounds, but that contributed very little to the overall inventory. It may also be useful to
have a target uncertainty in mind. For example, if the current total uncertainty is ±50.0 % and the company
wishes to reduce it to ±25.0%, then that target is useful in the analysis. As uncertainties of individual values
in the calculations are examined, values that have uncertainties that are already at or below the 25.0% target,
The largest contributors to uncertainty can be determined by multiplying the emission estimate by the
maximum error bound for each source. This would result in the upper bound emission estimate for the
particular source. The emission estimates can then be sorted by the largest contributors. This is demonstrated
in Table 5-8 for the example crude oil production facility.
Table 5-8. Emission Uncertainty Ranking for Onshore Oil Production Example
Maximum
Emissions Emissions,
(tonnes Maximum (tonnes
Source Type Source CO2e /yr) Uncertainty, % CO2e /yr) Ranking
Combustion Boiler/heaters 5,210 8.77 5,661 6
Sources Natural gas engines 14,100 15.6 16,244 4
Diesel engines 220 15.5 254
Flares 30,700 21.1 37,123 3
Fleet vehicles 129 19.2 154
Vented Sources Dehydration and Kimray pump vents 5,440 76.0 9,571 5
Tanks – flashing losses 40,300 88.7 76,039 1
Amine unit 66,700 9.77 73,216 2
Pneumatic devices 3,360 49.2 5,013 8
Chemical injection pumps 2,530 106 5,215 7
Vessel blowdowns 3.65 319 15.3
Compressor starts 38.7 187 111
Compressor blowdowns 17.3 175 47.6
Well workovers 0.939 294 3.70
Other non-routine (PRVs) 6.81 319 28.5
Fugitive Fugitive components
1,100 83.3 2,025 9
Sources
Fleet vehicle
1.30 112 2.75
refrigeration, R-314a
Indirect Electricity consumed
553 10.2 609 10
Emissions
Figure 5-1 provides an illustration of the emissions (which are not expressed in CO2e, but simply in tonnes of
each type of gas) for the onshore oil field which is a reprint from the API Compendium.
30,000
20,000
10,000
0
Combustion Vented Fugitive Indirect
Emission Source Category
5.4.1 Opportunities to Improve the Uncertainty Estimate for the Onshore Oil Field
Each facility should examine the major categories of emissions and emission uncertainty, and then examine
specific emission sources within the category. Prioritizing the largest sources of uncertainty can be done with
this simple approach.
Among the major types of emissions in Figure 5-2, there are three very significant emission categories:
1) combustion sources of CO2; 2) vented sources of CO2; and 3) vented sources of CH4. 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.
60,000
50,000
40,000
30,000
20,000
10,000
0
Combustion Vented Fugitive Indirect
Emission Source Category
Figure 5-2. Onshore Oil Field: Summary of CO2 Equivalent Emissions and Uncertainties
Examination of the uncertainties, as shown in the bars on Figure 5-2, or the absolute uncertainty values in
Table F-15, shows that “vented emissions of CH4” are the most significant source of uncertainty, contributing
36,400 CO2e tonnes of uncertainty (based on 2600 ±66.5% CO2e tonnes) in a total inventory that had only
37,300 CO2e tonnes of uncertainty (based on 170,000 ±21.9% CO2e tonnes). Carbon dioxide combustion
sources are the next largest contributor of uncertainty, contributing 6,780 CO2e tonnes of uncertainty. The
third largest source is CO2 vented, with 4,420 CO2e tonnes of uncertainty (based on 63,600 ±6.95% CO2e
tonnes).
Therefore, should the company wish to reduce uncertainty in the GHG emission inventory from this onshore
oil field facility, these categories would be the primary targets for uncertainty reduction. These appear to be
the sources where uncertainty reduction efforts could be effectively spent to refine the inventory and reduce
the uncertainty.
Within vented emissions, there are 10 sources listed in Table 5-8. The highest ranking source, tank flashing
losses, contributes a significant part of the vented emissions (33% of total vented emissions), as well as the
Examination of the detailed calculation for that category (presented earlier), shows that the largest uncertainty
is in the emission factor, which is a general industry-wide emission factor. The uncertainty in that factor is
±90.4%. This uncertainty can be reduced simply by using an improved estimation method to determine tank
flashing losses.
As elaborated in the API Compendium, other emission estimation methodologies can be used to estimate
emissions with lower uncertainties. For example, if the “EUB Rule of Thumb” approach (API Compendium
Section 5.4.1) were applied instead of the default tank flashing emission factor for the example onshore
production facility, the emissions for this category would be 20,500 ±49.7% tonnes CO2e.2 This one revision
would change the overall inventory emissions to 151,000 ±9.88% tonnes CO2e.
The company may also decide to take direct measurements of methane emissions from the tanks or other
vented sources with large uncertainties, such as amine unit emissions and dehydration unit emissions.
Repeated direct measurements or a site-specific emission factor would likely reduce uncertainty. In some
cases, the company might be able to produce an improved emission factor with lower uncertainty without
direct emission measurements. For example, dehydrator emission factors were produced by computer
simulations of dehydrators using national average input data; the company might elect to produce simulations
specifically for their dehydrators, using their dehydrator operating conditions, and thus produce a company
specific emission factor with less uncertainty.
• Boilers/heaters;
• Natural gas engines;
• Diesel engines;
• Flares; and
• Fleet vehicles.
2
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.
By examining the calculations used for flaring, the following general strategies could be selected by the
operating company to reduce uncertainty in flaring:
• Refine the measurement of total gas flared (the activity factor), to reduce the activity factor
uncertainty from the current value of 15% and/or;
• Refine the gas composition measurements to reduce the uncertainty from 4%.
Refining the measurement of total gas flared may result from many approaches, such as improving the meter
quality (even possibly replacing the meter), improving the quality control of the existing meter (such as
number of calibrations and inspections), or improving the number of measurements taken and recorded from
the meter that is used to calculate the total gas flared (this assumes measurements were not already
continuous). These approaches have varying costs, and some may be cost-prohibitive. The company would
have to determine which approach was the most cost effective.
Refining the gas composition data may also come from several methods, such as taking additional routine
samples, installing a continuous gas analyzer, or using a better analysis method. As with the gas flow rate
measurement, these approaches have varying costs, and some may be cost-prohibitive. The company would
have to determine which approach was the most cost effective.
If the company was effective in reducing the uncertainty in flare gas CO2 emissions, it might then elect to
proceed to the next largest CO2 uncertainty source. The end-user will have to recalculate total emissions for
the company or facility each time, and determine if the uncertainty goal or target has been reached.
If the company decided to take the next step of emission reductions, it may target the next largest category of
uncertainty, which is vented CO2 sources. Vented CO2 emissions result from seven sources listed in
Table F-15. The amine unit accounts for the over 98% of the emissions in this category.
Amine unit CO2 vented emissions have an uncertainty of 4,360 tonnes of CO2e in a category that only has an
uncertainty of 4,420 tonnes of CO2e. Therefore any uncertainty reduction efforts in this area should be
directed at the amine unit calculations. However, the calculations used for the amine units reveal that the
uncertainty is actually relatively low. The uncertainty in that emission category is less than 7%. It is possible
to reduce this uncertainty by improving the uncertainties for gas composition, and gas throughput. While it is
possible to reduce this uncertainty, given the low percentage level of the uncertainty, any improvement may
be small. Therefore this is likely a case where company efforts are best spent on the previous categories.
This Appendix provides a list of over 300 hydrocarbon measurement standards used to measure quantities of
petroleum products and natural gas liquids. Entire series of standards from particular standards developers are
listed as many of the standards are inter-related and have to be used in conjunction with one another for
measurement of quantities. This list is provided for users’ convenience to highlight the fact that companies may
use different sets of measurement standards that fit their operating environment and available instrumentation.
Each of the standards cited is designed to meet a variety of data needs and a defined level of accuracy, precision
and uncertainty ranges. The compilation provided here includes standards developed by the American
Petroleum institute (API); the American Society of Testing and Materials (ASTM)’ the American Gas
Association (AGA); the Gas Processors Association (GPA), and the international organization for
standardization (ISO). The list does not comprise an exhaustive reference to all standards; it does not cite
applicable standards from the EU, Japan, or other standards setting organizations.
Chapter 1 Vocabulary
Chapter 2.2A Measurement and Calibration of Upright Cylindrical Tanks by the Manual Strapping
Method
[Chapter 2.2A should be used in conjunction with Chapter 2.2B. These two standards
combined supersede the previous API Standard 2550, Measurement and Calibration of
Upright Cylindrical Tanks]
Chapter 2.2B Calibration of Upright Cylindrical Tanks Using the Optical Reference Line Method
Chapter 2.2C Calibration of Upright Cylindrical Tanks Using the Optical-Triangulation Method
Chapter 2.2D Calibration of Upright Cylindrical Tanks Using the Internal Electro-optical Distance
Ranging Method
Chapter 2.2E Petroleum and Liquid Petroleum Products—Calibration of Horizontal Cylindrical Tanks—
Part 1: Manual Methods
Chapter 2.2F Petroleum and Liquid Petroleum Products—Calibration of Horizontal Cylindrical Tanks—
Part 2: Internal Electro-optical Distance-ranging Method
Std 2552 Measurement and Calibration of Spheres and Spheroids
Std 2554 Measurement and Calibration of Tank Cars
Std 2555 Liquid Calibration of Tanks
RP 2556 Correcting Gauge Tables for Incrustation
Chapter 2.7 Calibration of Barge Tanks
Chapter 2.8A Calibration of Tanks on Ships and Oceangoing Barges
a. ASTM Test Method 1945-03, Analysis of Natural Gas by Gas Chromatography (July 2003)
This test method covers the determination of the chemical composition of natural gases and similar gaseous
mixtures within the range of composition shown in Table D-1. This test method may be abbreviated for the
analysis of lean natural gases containing negligible amounts of hexanes and higher hydrocarbons, or for the
determination of one or more components, as required.
COMPONENT MOL %
Helium 0.01 to 10
Hydrogen 0.01 to 10
Oxygen 0.01 to 20
Nitrogen 0.01 to 100
Carbon dioxide 0.01 to 20
Methane 0.01 to 100
Ethane 0.01 to 100
Hydrogen sulfide 0.3 to 30
Propane 0.01 to 100
Isobutane 0.01 to 10
n-Butane 0.01 to 10
Neopentane 0.01 to 2
Isopentane 0.01 to 2
n-Pentane 0.01 to 2
Hexane isomers 0.01 to 2
Heptanes+ 0.01 to 1
Components in a representative sample are physically separated by gas chromatography (GC) and compared
to calibration data obtained under identical operating conditions from a reference standard mixture of known
composition. The numerous heavy-end components of a sample can be grouped into irregular peaks by
reversing the direction of the carrier gas through the column at such time as to group the heavy ends either
as C5 and heavier, C6 and heavier, or C7 and heavier. The composition of the sample is calculated by
comparing either the peak heights, or the peak areas, or both, with the corresponding values obtained with
the reference standard. This test method is of significance for providing data for calculating physical
properties of the sample, such as heating value and relative density, or for monitoring the concentrations of
one or more of the components in a mixture.
COMPONENT
(MOLE %) REPEATABILITY REPRODUCIBILITY
0 to 0.09 0.01 0.02
0.1 to 0.9 0.04 0.07
1.0 to 4.9 0.07 0.10
5.0 to 10 0.08 0.12
Over 10 0.10 0.15
1
ASTM 1945-03, July 2003
This practice covers the determination of the chemical composition of reformed gases and similar gaseous
mixtures containing the following components: hydrogen, oxygen, nitrogen, carbon monoxide, carbon
dioxide, methane, ethane, and ethylene.
Components in a sample of reformed gas are physically separated by gas chromatography and compared to
corresponding components of a reference standard separated under identical operating conditions, using a
reference standard mixture of known composition. The composition of the reformed gas is calculated by
comparison of either the peak height or area response of each component with the corresponding value of
that component in the reference standard.
The chemical composition data can be used to calculate physical properties of the gas and its
interchangeability with other fuel gases. The quality of data obtained by this method depends on the
preparation of moisture-free and homogenous mixtures of known composition for comparison with the test
sample. The fraction of a component in the reference standard should not be < 0.5 mole%, nor differ by
more than 10 mole%, from the fraction of the corresponding component in the tested sample, and its
composition should be known to within 0.01 mole% for any component.
Method precision in terms of reproducibility and repeatability are provided in Table D-3 below.
COMPONENT
(MOLE %) REPEATABILITY REPRODUCIBILITY
0 to 1 0.05 0.1
1 to 5 0.1 0.2
5 to 25 0.3 0.5
Over 25 0.5 1.0
This method is for determining the composition of refinery gas samples or expanded liquefied petroleum gas
(LPG) samples that are obtained from refining processes or natural sources. It provides individual results
for non-condensable gases, hydrogen sulfide, C1 through C4 hydrocarbons and C5 paraffins, while C5
olefins and C6+ hydrocarbons are provided as a composite. The method yields quantitative results from 0.1
to 99.9 mole% for a single component or composite, except for hydrogen sulfide that yields quantitative
results between 0.1 and 25 mole%.
d. ISO 6974, Natural Gas - Determination of composition with defined uncertainty by gas chromatography
(October 2002)
Part 3: Determination of hydrogen, helium, oxygen, nitrogen, carbon dioxide and hydrocarbons up to C8
using two packed columns.
Part 4: Determination of nitrogen, carbon dioxide and C1 to C5 and C6+ hydrocarbons for a laboratory and
on-line measuring system using two columns.
Part 5: Determination of nitrogen, carbon dioxide and C1 to C5 and C6+ hydrocarbons for a laboratory and
on-line process application using three columns.
Part 6: Determination of hydrogen, helium, oxygen, nitrogen, carbon dioxide, and C1 to C8 hydrocarbons
using three capillary columns.
ISO 6974 describes a gas chromatographic methods for the quantitative determination of the content of
hydrogen, helium, oxygen, nitrogen, carbon dioxide, and C1 to C8 hydrocarbons in natural gas samples
using two or three packed or capillary columns combinations. They are applicable to the analysis of gases
containing constituents within the mole fraction ranges given in the respective parts of the standard, and is
commonly used for laboratory applications. These ranges do not represent the limits of detection, but the
limits within which the stated precision of the method applies. Although one or more components in a
sample may not be present at detectable levels, the method can still be applicable. This method can also be
applicable to the analysis of natural gas substitutes.
For example, the set-up for the determination of hydrogen, helium, oxygen, nitrogen, carbon dioxide, and
hydrocarbons from C1 to C8 by gas chromatography using three capillary columns (part 6 of the standard) is
as follows:
− A PLOT precolumn is used for the separation of carbon dioxide (CO2) and ethane (C2H6);
− A thick film WCOT2) column coated with an apolar phase is used for the separation of the C3 to C8
(and heavier) hydrocarbons.
The permanent gases helium (He), hydrogen (H2), oxygen (O2), nitrogen (N2), and methane (CH4) are
detected with a thermal conductivity detector (TCD). The C2 to C8 hydrocarbons are detected with a flame
ionization detector (FID).
For the analysis of natural gas substitutes, carbon monoxide (CO) and carbon dioxide (CO2) are detected
using an FID after reduction of the components to CH4 by a methanizer. Use of a methanizer makes it
possible to detect CO and CO2 at mole fractions greater than 0,001%. If the samples do not include CO or
CO2, or if the CO and/or the CO2 mole fraction exceeds 0,02%, a methanizer is not required. CO and CO2
may then alternatively be detected with the TCD.
When analyzing natural gas substitutes, the PLOT column described in 3.1 can also be used for the
separation of ethyne (C2H2) and ethene (C2H4) and the molecular sieve PLOT column can also be used for
the analysis of carbon monoxide (CO).
Typical precision values for this method are provided in Table A.1 of the standard. These values have been
obtained from practical experience and give an indication of the performance of the method. As such, they
cannot be compared with precision values mentioned in informative annexes of other parts of ISO 6974, as
they are very much dependent on the quality of the calibration gases used and the laboratory skills.
For specifies concentrations < 1.0 mole%, the relative repeatability and reproducibility is expected to be 2
and 4%, respectively. For higher concentrations, ranging from 1-50 mole%, the relative repeatability and
reproducibility are around 0.8 and 1.6% respectively.
e. ASTM D2650 – 04, Standard Test Method for Chemical Composition of Gases By Mass Spectrometry
(November 2004)
This test method is applicable for the quantitative analysis of gases containing specific combinations of the
following components: hydrogen; hydrocarbons with up to six carbon atoms per molecule; carbon
monoxide; carbon dioxide; mercaptans with one or two carbon atoms per molecule; hydrogen sulfide; and
air (nitrogen, oxygen, and argon). This test method is not applicable for the determination of constituents
that are present in amounts less than 0.1 mole %. This test method was developed on a specific type of
Mass Spectrometer, thus users of other instruments may have to modify operating parameters and the
calibration procedure and adapt it to their instrument.
The method sets out the experimental procedures, while the calculation procedures will depend on the
knowledge of qualitative mixture composition; errors due to components presumed absent; minimum cross
f. ASTM D5291 – 02, Standard Test Methods for Instrumental Determination of Carbon, Hydrogen, and
Nitrogen in Petroleum Products and Lubricants, (2007)
This ASTM standard covers the simultaneous determination of carbon, hydrogen, and nitrogen in petroleum
products and lubricants. The results, expressed as total carbon, total hydrogen, and total nitrogen, are useful
in determining the complex nature of sample types covered by this test method, and can also be used to
determine the total carbon content that could be converted to CO2 upon combustion of the product.
These test methods are applicable to samples such as crude oils, fuel oils, additives, and residues for carbon
and hydrogen and nitrogen analysis. These test methods were tested in the concentration range of at least 75
to 87-wt% for carbon, at least 9 to 16-wt% for hydrogen, and 0.1 to 2-wt% for nitrogen. These test methods
are not recommended for the analysis of volatile materials such as gasoline, gasoline-oxygenate blends, or
gasoline type aviation turbine fuels.
a. ASTM D4891 – 89, Test Method for Heating Value of Gases in Natural Gas Range by Stoichiometric
Combustion, (Reapproved 2006)
This test method covers the determination of the heating value of natural gases and similar gaseous mixtures
within the range of composition shown in Table D-4.
Table D-4. ASTM D4891-89 Range of Composition for Natural Gas Components
This test method provides an accurate and reliable procedure to measure the total heating value of a fuel gas,
on a continuous basis, which is used for regulatory compliance, custody transfer, and process control. Some
instruments which conform to the requirements set forth in this test method can have response times on the
order of 1 min or less and can be used for on-line measurement and control. The method is sensitive to the
presence of oxygen and unsaturated hydrocarbons. For components not listed and composition ranges that
fall outside those in Table 7, such as in process or refinery fuel gases, modifications in the method may be
required to obtain correct results.
In testing the precision and accuracy of this method, repeatability within a laboratory was shown to be 0.76
Btu/scf, with the corresponding 95% confidence of the repeatability interval being 2.1 Btu/scf.
Reproducibility between laboratories was determined to be 1.67 Btu/scf with the corresponding 95%
confidence reproducibility interval 5.1 Btu/scf. The average bias of all measurements agreed with the
average reference value to within 0.1%.
b. ASTM D1826 – 94, Standard Test Method for Calorific (Heating) Value of Gases in Natural Gas Range
by Continuous Recording Calorimeter, (Reapproved 2003)
This test method covers the determination – with the continuous recording calorimeter – of the total calorific
(heating) value of fuel gas produced or sold in the natural gas range from 900 to 1200 Btu/scf. The heating
value is determined by imparting the heat obtained from the combustion of the test gas to a stream of air and
measuring the rise of the air temperature. The streams of test gas and heat absorbing air are maintained in
fixed volumetric proportion to each other by metering devices similar to the ordinary wet test meters geared
together and driven from a common electric motor. The meters are mounted in a tank of water, the level of
which is maintained and the temperature of which determines the temperature of the entering gas and air.
The experimental set-up is such that the temperature rise produced in the heat-absorbing air is directly
proportional to the heating value of the gas. Since all the heat from the combustion of the test gas sample,
including the latent heat of vaporization of the water vapor formed in the combustion, is imparted to the
heat-absorbing air, the calorimeter makes a direct determination of total heating value. The temperature rise
is measured by nickel resistance thermometers and is translated into Btu/scf.
This test method provides an accurate and reliable method to measure the total calorific value of a fuel gas,
on a continuous basis, which is used for regulatory compliance, custody transfer, and process control. As
far as precision, the calorimeters tested were standardized with methane weekly, and a rigid control was
c. ASTM D7313 – 08, Standard Practice for Determination of the Heating Value of Gaseous Fuels using
Calorimetry and On-line/At-line Sampling (May 2008)
This practice is used for the determination of the heating value measurement of gaseous fuels using a
calorimeter with at-line and in-line instruments that are operated from time to time on a continuous basis.
This type of near-real time monitoring systems that measure fuel gas characteristics such as heating value
are prevalent in various gaseous fuel industries and in industries either producing or using gaseous fuel in
their industrial processes. The installation and operation of particular systems would vary depending on
process type, regulatory requirements, and the user’s objectives and performance requirements.
In operating the system, a representative sample of the gaseous fuel is extracted from a process pipe, a
pipeline, or other gaseous fuel stream and is transferred to an analyzer sampling system. After conditioning
that maintains the sample integrity, the sample is introduced into a calorimeter. Excess extracted process or
sample gas is vented to the atmosphere, a flare header, or is returned to the process in accordance with
applicable economic and environmental requirements and regulations. Post combustion gasses from the
calorimeter are typically vented to the atmosphere. The heating value is calculated based upon the
instrument’s response to changes in the heating value of the sample gas using an applicable computation
algorithm. This practice is providing guidance for standardized start-up procedures, operating procedures,
and quality assurance practices for calorimeter based on-line, at-line, in-line and other near-real time
heating-value monitoring systems.
d. ASTM D4809 – 06, Standard Test Method for Heat of Combustion of Liquid Hydrocarbon Fuels by
Bomb Calorimeter (Precision Method), (2006)
The heat of combustion is a measure of the energy available from a fuel, and it is essential when considering
the thermal efficiency of equipment for producing either power or heat. This procedure measures the mass
heat of combustion, namely, the heat of combustion per unit mass of fuel. The volumetric heat of
combustion, that is, the heat of combustion per unit volume of fuel, can be calculated by multiplying the
mass heat of combustion by the density of the fuel (mass per unit volume).
This test method covers the determination of the heat of combustion of hydrocarbon fuels, and is designed
for high precision with the difference between duplicate determinations in the order of 0.2%, for aviation or
turbine fuels, or greater differences for a wide range of volatile and nonvolatile materials.
In order to attain this high precision, strict adherence to all details of the procedure is essential since the
error contributed by each individual measurement that affects the precision ought to be kept below 0.04%,
e. ASTM D5865 – 07a, Standard Test Method for Gross Calorific Value of Coal and Coke, (2007)
The gross calorific value can be used to compute the total calorific content of the quantity of coal or coke,
and it can also be used for computing the calorific value versus sulfur content to determine whether the coal
meets regulatory requirements for industrial fuels.
This test method pertains to the determination of the gross calorific value of coal and coke by either an
isoperibol or adiabatic bomb calorimeter. The resulting values are to be regarded as standard, and no other
units of measurement are included in this standard.
Methods Cited
ASTM 1945 – 03, “Analysis of Natural Gas by Gas Chromatography”, current edition approved May 10, 2003.
Published July 2003. Originally approved in 1962. Last previous edition approved in 2001 as D1945–96(2001).
2ASTM D 1946 – 90 (Reapproved 2006), Standard Practice for Analysis of Reformed Gas by Gas
Chromatography, current edition approved June 1, 2006. Published June 2006 (Originally approved in 1962)
ASTM UOP UOP539-97, “Refinery Gas Analysis by Gas Chromatography”
ISO 6974, 2002, “Natural gas - Determination of composition with defined uncertainty by gas chromatography”
(in 6 parts); First edition 2002-10-15, Geneva, Switzerland
ASTM D2650, 2004, “Standard Test Method for Chemical Composition of Gases by Mass Spectrometry”,
Published November 1, 2004.
ASTM D5291 – 02, (2007) “Standard Test Methods for Instrumental Determination of Carbon, Hydrogen, and
Nitrogen in Petroleum Products and Lubricants”, Published 2007.
ASTM D 4891 – 89 (Reapproved 2006), Standard Test Method for Heating Value of Gases in Natural Gas
Range by Stoichiometric Combustion, Current edition approved June 1, 2006. Published June 2006. Originally
approved in 1989. Last previous edition approved in 2001 as D4891–89 (2001)
ASTM D 1826 – 94 (Reapproved 2003), “Standard Test Method for Calorific (Heating) Value of Gases in
Natural Gas Range by Continuous Recording Calorimeter”, Current edition approved May 10, 2003. Published
May 2003. Originally approved in 1961. Last previous edition approved in 1998 as D 1826 – 94 (1998).
ASTM D7313 – 08, “Standard Practice for Determination of the Heating Value of Gaseous Fuels using
Calorimetry and On-line/At-line Sampling”, Current edition approved May 1, 2008. Published May 2008.
ASTM D4809 – 06, “Standard Test Method for Heat of Combustion of Liquid Hydrocarbon Fuels by Bomb
Calorimeter (Precision Method)”, Published 2006.
ASTM D4809 – 06, “Standard Test Method for Heat of Combustion of Liquid Hydrocarbon Fuels by Bomb
Calorimeter (Precision Method)”, Published 2007.
UNITS CONVERSION
Appendix E
UNITS CONVERSION3
Energy units
Quantities
• 1.0 calorie = 4.187 J
• 1.0 gigajoule (GJ) = 109 joules = 0.948 million Btu = 239 million calories = 278 kWh
• 1.0 British thermal unit (Btu) = 1055 joules (1.055 kJ)
• 1.0 Quad = One quadrillion Btu (1015 Btu) = 1.055 exajoules (EJ)
− Approximately 172 million barrels of oil equivalent (boe)
• 1000 Btu/lb = 2.33 gigajoules per tonne (GJ/t)
• 1000 Btu/US gallon = 0.279 megajoules per liter (MJ/l)
Power
• 1.0 watt = 1.0 joule/second = 3.413 Btu/hr
• 1.0 kilowatt (kW) = 3413 Btu/hr = 1.341 horsepower
• 1.0 kilowatt-hour (kWh) = 3.6 MJ = 3413 Btu
• 1.0 horsepower (hp) = 550 foot-pounds per second = 2545 Btu per hour = 745.7watts = 0.746kW
Energy Costs
• $1.00 per million Btu = $0.948/GJ
• $1.00/GJ = $1.055 per million Btu
Fossil fuels4
• Barrel of oil equivalent (boe) = approx. 6.1 GJ (5.8 million Btu), equivalent to 1700 kWh
− "Petroleum barrel" is a liquid measure equal to 42 U.S. gallons (35 Imperial gallons or 159 liters);
about 7.2 barrels oil are equivalent to one tonne of oil (metric) = 42-45 GJ.
• Gasoline: LHV = 115,000 Btu/gallon = 121 MJ/gallon = 32 MJ/liter; HHV: 125,000 Btu/gallon = 132
MJ/gallon = 35 MJ/liter
• Metric tonne gasoline = 8.53 barrels = 1356 liter = 43.5 GJ/t (LHV); 47.3 GJ/t (HHV)
3
Based on data from the ORNL Bioenergy Feedstock Information Network; [Link]
4
The energy content (heating value) of petroleum products per unit mass is fairly constant, but their density differs
significantly, hence their energy content is different
5
The energy content (heating value) per unit mass varies greatly between different "ranks" of coal. "Typical" coal (rank
not specified) usually means bituminous coal, the most common fuel for power plants (27 GJ/t).
6
The average carbon content values above are indicative only. More accurate determinations should be sought from fuel
provider, or by testing fuels used.
As stated in Section 5.0, this Appendix provides the source-by-source calculation and aggregation of
uncertainty for a hypothetical facility.7
The following summarizes the characteristics of the onshore oil field example facility introduced in
Section 5. This example field is described in detail in the API Compendium (API, 2009). The calculation
exhibits that are reprinted directly from the API Compendium but have additional uncertainty information
added, are indicated by the highlighted text. Excerpts of Tables 5-1 and 5-2 are provided to reiterate
information used in the uncertainty calculations.
Facility Description: An onshore oil field in Texas consists of 320 producing oil wells.
Throughput: The average daily oil and gas production rates are 6,100 bbl/day and 30×106 scf/day,
respectively.
Operations: The facility operates approximately 343 days per year. The facility imports 917 MW-hr from
the eGRID subregion ERCOT all annually. The facility gas composition is presented in Table F-1 and
results in a heating value of 928 Btu/scf with an uncertainty of ±4%, based on engineering
judgment.
Uncertainty a
Gas Compound Produced Gas Mole % (±%)
CO2 12 4
N2 2.1 4
CH4 80 4
C2H6 4.2 4
C3H8 1.3 4
C4H10 0.4 4
Footnote:
a
Uncertainty is based on engineering judgment at a 95% confidence interval.
Note: the values shown above are for example only. They do not reflect average operations.
7
Note, the calculations shown in this Appendix are rounded to three significant figures. Actual calculations were
compiled in a spreadsheet with rounding reserved for the final results.
• Fuel consumption for boilers, turbines, and flares = ±15% (assigned by expert judgment in
Table 5-1).
• Heating Value of Gas Combusted = 928 Btu/scf ±4% (measured independent of gas composition;
assigned by expert judgment here).
• Gas Composition measurement = ±4% (determined by analysis of repeat samples using
techniques from Section 4.0; assigned by expert judgment here).
• Unit Capacity of Heaters and Reboilers = ±6.71% (Calculated for Table 5-1).
F.1 Uncertainty in Natural Gas Combustion Devices – CO2 Emissions
Carbon dioxide emissions are estimated based on the volume of fuel consumed and the fuel carbon
content. Boilers, heater/reboilers, and compressor engines-turbines use natural gas at this facility.
However, because CH4 and N2O emission factors are based on the type of equipment as well as the fuel
consumed, the CO2 emission estimates are also calculated for common equipment types.
We will first examine boilers and heaters/reboilers. The quantity of natural gas consumed by this
equipment is based on the parameters shown in Table F-2.
Annual
Average Activity
No. Unit Operation Factor (all
of Uncertainty Capacity Uncertainty (per unit Uncertainty units Uncertainty
Source Fuel Units (±%)a (per unit) (±%)a per year) (±%)a combined) (±%)a
Boilers Produced 40×106
6 ±0% units N/A N/A ±15% scf/yr
Gas scf/yr
Heaters/ Produced 2×106 343 53.2×106 ±6.71%
3 ±0% units ±5% Btu/hr ±2% days/yr
reboilers Gas Btu/hr days/yr scf/yr scf/yr
The volume of fuel combusted (V) for Boilers and Heaters/Reboilers is:
The uncertainty for the heaters/reboilers is examined first, then the combined natural gas uncertainty is
estimated.
The uncertainty for the boilers fuel consumption is 15% based on expert judgment.
Combined Uncertainty for Boilers and Heaters/Reboilers Activity Data: (Equation 4-4, Absolute
Uncertainties)
The natural gas usage is summed. Therefore, the uncertainty for the combined fuel consumption by the
boiler and heaters/reboiler is calculated by applying Equation 4-4 and using the absolute uncertainty values.
The uncertainty associated with the fuel consumed by turbines is assigned by expert judgment, as shown
in Table F-3. In this case, no further calculations are required for the turbine activity data.
The carbon content of the natural gas is needed to complete the estimation of CO2 emissions from natural
gas combustion. The gas composition data for this example are sown in Table F-4. Calculation of the
uncertainty values is described below.
For simplicity, we assume a 4% uncertainty on the mole% for all of the natural gas components. In reality,
the uncertainty of the gas composition should be determined through multiple samples, as shown in
Section 4.
U (abs )
∑ (mole%×MW)
= ∑ U(abs) 2
mole%×MW
The uncertainty for the wt% C for the individual components is calculated by applying Equation 4-6 and
using the relative uncertainties.
The uncertainty for the sum of the wt%C is calculated by applying Equation 4-4 and using the absolute
uncertainties.
Carbon dioxide emissions can then be calculated using mass balance approach. Emissions are calculated
below, by equipment type.
Boilers and Heaters:
93.2 ×106scf fuel lbmole 20.77 lb fuel 0.6124 lb C lbmole C
E CO2 = × × × ×
yr 379.3 scf fuel lbmole fuel lb fuel 12.01 lb C
lbmole CO 2 44.01 lb CO 2 tonne
× × × = 5, 200 tonnes CO 2 / yr
lbmole C lbmole CO 2 2204.62 lb
As shown earlier, 93.2×106 scf/yr has an uncertainty of ±7.49%. The uncertainty in the carbon content of
the fuel mixture is 3.71% and the uncertainty of the MW of the mixture is 2.69%. The uncertainty in the
CO2 emissions is calculated by applying Equation 4-6 and using the relative uncertainty values.
U (rel ) Emissions = U(rel ) 2Boilers and Heaters Gas Usage +U(rel ) Carbon
2 2
Content +U(rel ) MW mix
U ( r e l ) E m is s io n s = U ( r e l ) 2A F + U ( r e l ) C2 a rb o n C o n te n t + U ( r e l ) 2M W m ix
U ( r e l ) E m is s io n s = 1 5 + 3 .7 1 + 2 .6 7 = ± 1 5 .7 % (to n n e s C O 2 /yr)
2 2 2
The diesel firing rate is the sum of the activity factors for the generators and pumps as shown in
Table F-5.
Annual
Average Activity
No. Unit Operation Factor (all
of Uncertainty Capacity Uncertainty (per unit Uncertainty units Uncertainty
Source Fuel Units (±%)a (per unit) (±%)a per year) (±%)a combined) (±%)a
Emergency
±10% 2,912 ×106 ±12.3%
generator IC Diesel 1 ±0% units 1800 hp ±5% hr 200 hr/yr
(hr/yr) Btu/yr Btu/yr
engine
Fire water ±10%
24 hr/yr; 77.7 ×106 ±13.2%
pump IC Diesel 1 ±0% units 460 hp ±5% hr (hr/yr);
87% load Btu/yr Btu/yr
engine ±20% (load)
Because the diesel emission factor is provided on a heat input basis, the equipment ratings are converted
to volume of fuel consumed (V) on a heat input basis.
We assume there is no uncertainty in the number of units, which here is the count of diesel engines. The
uncertainty in the unit capacity is 5%, the uncertainty in the energy efficiency of the engine is 5%, and the
uncertainty in the average hours of operation is 10% based on expert judgment. The uncertainty in the CO2
emissions is calculated by applying Equation 4-6 and using the relative uncertainty values.
U ( rel ) Activity Factor = U(rel ) 2Number Of Units +U(rel ) 2Unity Capacity +U(rel ) 2Efficiency +U(rel ) 2Average Operations
The uncertainty for pumps includes one additional term: the uncertainty in the percent load is 20%,
2
U (abs ) Diesal Firing Rate = (U(abs )Generator +U(abs ) 2Pump
U (abs ) Diesal Firing Rate = (0.123 × 2912 × 106 ) 2 + (0.235 × 77.7 × 106 ) 2 = 356.85 ×106 Btu / yr
356.85×106 Btu/yr
U (rel ) = 100%× = ± 11.9% ( Btu/yr )
2989.7×106 Btu/yr
The uncertainty in the emissions factor is determined by engineering judgment to be 10%. The total
uncertainty for the CO2 emissions is then is calculated by applying Equation 4-6 and using the relative
uncertainty values.
2
U (rel )CO2 = U(rel ) AF 2
+U(rel ) EF = 11.92 + 10.02 = ±15.6% ( tonnes CO 2 /yr )
F.3 Uncertainty in Stationary Combustion Devices – CH4, N2O, and CO2e Emissions
In this hypothetical example, CH4 and N2O emissions from combustion are calculated using the fuel-based
emission factors presented in the API Compendium. The calculations for CH4 and N2O emissions from
combustion are reprinted directly from the API Compendium. Some uncertainty values were assigned as
follows:
Uncertainty Estimate of CH4 and N2O Emissions from Boilers and Heater/Reboilers:
Natural gas emission factors for CH4 and N2O are based on the energy consumed by the equipment. This
differs from the CO2 emission estimates above, which were based on the volume of natural gas consumed.
So, the fuel usage information from Table F-2 must first be converted to an energy input basis.
The uncertainty for the boiler activity factor is calculated by applying Equation 4-6 and using the relative
uncertainty values.
The uncertainty for the combined activity factor is calculated by applying Equation 4-4 and using the
absolute uncertainty values.
2 2
U (abs ) AF = (U(abs ) Boiler +U(abs ) Heater
U (abs ) AF = (0.155 × 37,120 × 106 ) 2 + (0.0539 × 49,392 × 106 ) 2 = 6,347 × 106 Btu/yr
6,347×106 Btu/yr
U (rel ) AF = 100% × = ±7.34% ( Btu/yr )
86,512×106 Btu/yr
As stated previously, the uncertainty in the emissions factor for N2O is determined by engineering judgment
to be 150%. The uncertainty in the emissions factor for methane is 25% based on expert judgment. The
uncertainty for emissions is calculated by applying Equation 4-6 and using the relative uncertainty values.
Because it is not possible to have negative emissions, the lower bound uncertainty will be truncated at zero.
As a result, the lower and upper bound uncertainties are estimated separately.
2
U (rel )emissions = U(rel ) AF 2
+U(rel )EF = 7.342 + 1502 = +150%, −100% ( tonnes N 2 O/yr )
The final step is to convert the emissions of CO2, CH4, and N2O to a CO2e. Note there are no uncertainties
for any of the GWPs, so the combined uncertainty is based on summing the individual emissions and
applying Equation 4-4 (using absolute uncertainty values). Also, because the N2O uncertainty estimate is
asymmetrical, separate upper and lower bound uncertainties are calculated for the CO2e. emissions.
U (abs ) ECO e = (0.0878 × 5, 200) 2 + (0.261× 21× 0.0865) 2 + (1.50 × 310 × 0.0242) 2
2
Upper:
= 456 tonnes CO2e/yr
456 tonnes CO2e/yr
U (rel ) ECO e = 100% × = +8.77% ( tonnes CO2e/yr )
2
5, 210 tonnes CO2e/yr
U (abs ) ECO e = (−0.0878 × 5, 200) 2 + (−0.261× 21× 0.0865) 2 + (−1.00 × 310 × 0.0242) 2
2
Lower:
= 456 tonnes CO2e/yr
456 tonnes CO2e/yr
U (rel ) ECO e = 100% × = −8.77% ( tonnes CO2e/yr )
2
5, 210 tonnes CO2e/yr
Here also, the natural gas emission factors for CH4 and N2O are based on the energy consumed by the
equipment. The fuel usage information from Table F-3 is first converted to an energy input basis.
⎝ yr scf ⎠
Then the emission factors are applied to estimated CH4 and N2O emissions:
The uncertainty for the turbines’ activity factor is 15%. The uncertainty in the heating values is 5% based
on expert judgment. The uncertainty in the CH4 emissions factor is determined by engineering judgment to
be 25%, and +150%, -100% for the N2O emission factor. The uncertainty of the emissions is calculated by
applying Equation 4-6 and using the relative uncertainty values.
Because it is not possible to have negative emissions, the lower bound uncertainty will be truncated at zero.
As a result, the lower and upper bound uncertainties are estimated separately.
The final step is to convert the emissions of CO2, CH4, and N2O to a CO2e basis. Equation 4-4 is applied
using absolute uncertainty values. Because the N2O uncertainty estimate is asymmetrical, separate upper
and lower bound uncertainties are calculated for the CO2e emissions.
U (abs ) ECO e = (0.157 ×13,900) 2 + (0.294 × 21× 0.905) 2 + (1.51× 310 × 0.325)2
2
Upper:
U (abs ) ECO e = 2,190 tonnes CO 2 e/yr
2
U (abs ) ECO e = (−0.157 ×13,900) 2 + (−0.294 × 21× 0.905) 2 + (−1.00 × 310 × 0.325) 2
2
Lower:
U (abs ) ECO e = 2,190 tonnes CO 2 e/yr
2
For diesel engines, the emission factor is provided on a heat input basis. The equipment ratings provided in
Table F-5 are converted to volume of fuel consumed on a heat input basis using the conversion factor 8,089
Btu/hp-hr (from API Compendium Table 4-2). An uncertainty of 5% is assumed based on engineering
judgment for this conversion factor.
As discussed earlier, we assume there is no uncertainty in the number of units, the count of diesel engines.
The uncertainty in the unit capacity is 5%, the uncertainty in the hp-hr to Btu conversion factor is 5%, and
the uncertainty in the average hours of operation is 10% based on expert judgment. The uncertainty in the
emissions factors are determined by engineering judgment to be 25% for CH4 and +150%, –100% for N2O.
The uncertainty of the emissions is calculated by applying Equation 4-6 and using the relative uncertainty
values.
Because it is not possible to have negative emissions, the lower bound uncertainty will be truncated at zero.
As a result, the lower and upper bound uncertainties are estimated separately.
U (rel )emissions = U(rel) 2NumberOfUnits +U(rel ) 2UnityCapacity +U(rel ) 2Btu/hp-hr conv +U(rel )OperationHours
2
+U(rel ) 2EF
As discussed earlier in Table 5-1, we assume there is no uncertainty in the number of units, the count of
diesel engines. The uncertainty in the unit capacity is 5%, the uncertainty in the percent load is 20%, the
uncertainty in the heating value is 5%, and the uncertainty in the average hours of operation is 10% based on
expert judgment. The uncertainty in the emissions factors are determined by engineering judgment to be
25% for CH4 and +150%,-100% for N2O. For CH4, the CH4 content of the TOC is an assumed value, so the
uncertainty with this composition is assigned by engineering judgment to be 100%. The uncertainty of the
emissions is calculated by applying Equation 4-6 and using the relative uncertainty values.
U (rel )emissions = 02 + 202 + 52 + 52 + 102 + 252 + 1002 = +106%, −100% ( tonnes CH 4 /yr )
The final step is to convert the emissions of CO2, CH4, and N2O to a CO2e. Equation 4-4 is applied using
absolute uncertainty values. Both CH4 and N2O uncertainty estimates are asymmetrical, so both upper and
lower bound uncertainties are calculated for the CO2e emissions.
Upper:
U (abs ) ECO e = (U(abs ) 2ECO +U(abs ) 2ECH +U(abs ) 2ECH +U(abs ) 2E N O-Diesel>600hp +U(abs ) 2E N O-Diesel<600hp
2 2 4 -Diesel>600hp 4 -Diesel<600hp 2 2
U (abs) ECO e = (U(abs ) 2ECO +U(abs ) 2ECH -Diesel>600hp +U(abs ) 2ECH -Diesel<600hp +U(abs ) 2E N O-Diesel>600hp +U(abs ) 2E N O-Diesel<600hp
2 2 4 4 2 2
Annual
Unit Activity Factor
No. of Uncertainty Capacity Uncertainty (all units Uncertainty
Source Fuel Units (±%)a (per unit) (±%)a combined) (±%)a
Emergency Produced
1 0 N/A 500×106 scf/yr ±15% scf/yr
flare Gas
For this example, we assume the flare gas composition is equal to the field gas composition (shown in
Table F-4).
For CH4, the uncertainty in the activity factor is 15%, and the uncertainty in the noncombusted methane is
20%, both assigned by expert judgment. The uncertainty in the mole % of CH4 is 4%. The uncertainty of
the emissions is calculated by applying Equation 4-6 and using the relative uncertainty values.
N2O emissions from flares are estimated based on the volume of crude produced.
6,100 bbl 365 days 1.0 ×10−4 tonnes N 2 O
E N 2O = × × = 0.223 tonnes N 2 O/yr
day yr 1, 000 bbl
The uncertainty in the activity factor is 5% based on expert judgment. There is no uncertainty in the number
of days of operation. The uncertainty in the emissions factor is +200%,–100% based on expert judgment.
The uncertainty of the emissions is calculated by applying Equation 4-6 and using the relative uncertainty
values.
Because it is not possible to have negative emissions, the lower bound uncertainty will be truncated at zero.
As a result, the lower and upper bound uncertainties are estimated separately.
CO2 emissions result from CO2 present in the gas and CO2 formed through combustion.
The flare emissions for CO2 in gas and CO2 formed are correlated, since they both use the same activity
factor. To calculate the uncertainty, first rearrange the CO2 flare emissions as follows to eliminate the
correlation.
For the uncertainty aggregation, start first with the uncertainty for the moles of carbon in the flared gas
stream (the terms in parenthesis) by applying Equation 4-4 and using the absolute uncertainties. The
uncertainty in the gas stream composition is 4%.
U (abs )
∑ lbmoleC
= ∑U (abs) 2
mole %
Then calculate the uncertainty of the product of the composition and the combustion efficiency by applying
Equation 4-4 and using the relative uncertainty.
Next, account for the uncertainty of the CO2 present in the gas by applying Equation 4-4 and using the
absolute uncertainty values. This will aggregate uncertainty for all the terms in the brackets.
U (abs ) 2 = U (abs )12 + U (abs ) 2mole % = (0.203 × 0.939 × 0.98) 2 + (0.04 × 0.120) 2 = 0.187
0.187
U (rel ) 2 = 100% × = 18.0%
( 0.939 × 0.98) + 0.120
Finally the uncertainty of the emissions is calculated by applying Equation 4-6 and using the relative
uncertainty values. For the CO2 that is in the flared gas, the uncertainty in the activity factor is 15% based
on expert judgment.
The final step is to convert the emissions of CO2, CH4, and N2O to a CO2e. Equation 4-4 is applied using
absolute uncertainty values. Because the N2O uncertainty estimate is asymmetrical, separate upper and
lower bound uncertainties are calculated for the CO2e emissions.
Upper:
U (abs ) ECO e = (0.234 × 27, 400) + (0.253 × 21×153) 2 + (2.00 × 310 × 0.233) 2 = 6, 460 tonnes CO 2 e/yr
2
2
Lower:
U (abs ) ECO e = (−0.234 × 27, 400)2 + (−0.253 × 21×153)2 + (−1.00 × 310 × 0.233) 2 = 6, 460 tonnes CO2e/yr
2
The calculations for fleet vehicle emissions are reprinted directly from the API Compendium. Some
uncertainty values were assigned as follows:
• Miles Driven per Vehicle = ±15% (assigned by expert judgment, based upon uncertainties in the
vehicle logs gathered by the company).
2 2 2
U (rel ) Energy of Fuel = U(rel ) AF +U(rel )fuel economy +U(rel ) Heat value of gasoline
Carbon dioxide emissions are then calculated by applying the emission factor for vehicles.
2 2 2
U (rel )CO2 Emissions = U(rel ) Vehicle +U(rel )Energy of Fuel +U(rel ) EF
Methane emissions are then calculated by applying the emission factor for vehicles.
2 2 2
U (rel )CH 4Emissions = U(rel )Vehicle +U(rel ) Fuel Economy +U(rel ) EF
Nitrous oxide emissions are then calculated by applying the emission factor for vehicles.
2
U (rel ) N2OEmissions = U(rel ) Vehicle +U(rel ) 2Fuel Economy +U(rel ) 2EF
The final step is to convert the emissions of CO2, CH4, and N2O to a CO2e. Equation 4-4 is applied using
absolute uncertainty values. Both CH4 and N2O uncertainty estimates are asymmetrical, so both upper and
lower bound uncertainties are calculated for the CO2e emissions.
U (abs ) ECO e = (0.194 ×127) + (1.51× 21× 0.00643) 2 + (1.51× 310 × 0.00871)2 = 24.8 tonnes CO 2 e/yr
2
2
U (abs ) ECO e = (−0.194 × 127) + (−1.00 × 21× 0.00643)2 + (−1.00 × 310 × 0.00871) 2 = 24.7 tonnes CO2 e/yr
2
2
Table F-8 provides the operating parameters for gas dehydration at the example facility.
Annual
Average Activity
No. Unit Operation Factor (all
of Uncertainty Capacity Uncertainty (per unit Uncertainty units Uncertainty
Source Fuel Units (±%)a (per unit) (±%)a per year) (±%)a combined) (±%)a
Dehydration
vents (also
Produced 30×106 343 10,290×106 ±5.39%
has Kimray 1 ±0% units ±5% scf/day ±2% days/yr
Gas scf/day days/yr scf/yr scf/yr
pump
emissions)
Methane emissions from the gas dehydrator vents and for Kimray pumps are estimated using the production
segment emission factors provided in the API Compendium. For both of these sources, emissions are based
on the quantity of gas processed. Carbon dioxide emissions are estimated using the facility CO2 and default
CH4 concentrations in the gas.
For the uncertainty estimate, the emission factors should be combined to eliminate the correlation associated
with the common activity data. The uncertainty of the emissions factor is calculated by applying
30 × 106 scf 343 day 0.80 tonne mole CH 4 (facility) 0.024326 tonne CH 4
E CH4 = × × ×
day gas processed yr 0.788 tonne mole CH 4 (default) 106scf
The uncertainty of the emissions can be calculated by applying Equation 4-6 and using the relative
uncertainty values. As noted in Table 5-1, there is no uncertainty associated with the count of equipment.
The uncertainty associated with the volume of gas is 5% and the uncertainty associated with the days of
operation is 2%. The uncertainty in the CH4 content for the facility gas is 4%. The uncertainty for the
default CH4 content is 5.53% from the GRI/EPA Methane Study (Shires and Harrison, Volume 6,
Tables A-1, A-2, A-3, 1996), which is provided in Table E-4 of the API Compendium.
30 ×106 scf 343 day 0.024326 tonne CH 4 tonne mole CH 4 0.80 tonne mole CH 4 (facility)
E CO2 = × × 6
× ×
day gas processed yr 10 scf 16.04 tonne CH 4 0.788 tonne mole CH 4 (default)
tonne mole gas 0.12 tonne mole CO 2 44.01 tonne CO 2
× × ×
0.8 tonne mole CH 4 tonne mole gas tonne mole CO 2
The uncertainty of the CO2 emissions is calculated by applying Equation 4-6 and using the relative
uncertainty values. As noted in Table 5-1, there is no uncertainty associated with the count of equipment.
The uncertainty associated with the volume of gas is 5% and the uncertainty associated with the days of
operation is 2%. The uncertainty in the CO2 content for the facility gas is 4%. The uncertainty for the
The uncertainty of the CO2e emissions is calculated by applying Equation 4-4, using the absolute
uncertainty values.
U (abs ) ECO e = (0.775 ×106) + (0.775 × 21× 254) 2 = 4,130 tonnes CO 2e/yr
2
2
Table F-9 summarizes the operating parameters for storage tanks at the example facility.
Annual
Average Activity
No. Unit Operation Factor
of Uncertainty Capacity Uncertainty (per unit Uncertainty (all units Uncertainty
Source Fuel Units (±%)a (per unit) (±%)a per year) (±%)a combined) (±%)a
Central
6,100 343 2,092,300 ±5.39%
tank Crude Oil 1 ±0% units ±5% bbl/day ±2% days/yr
bbl/day days/yr bbl/yr bbl/yr
battery
Storage
Chemical 1 ±0% units N/A N/A N/A
tanks
Naphtha 1 ±0% units N/A N/A N/A
Glycol 1 ±0% units N/A N/A N/A
Methane emissions from flashing losses are estimated using the simple emission factor provided in the API
Compendium. (Note that this simple emission factor approach is used in absence of the more detailed
information necessary for the other calculation approaches provided in the API Compendium).
6,100 bbl 343 day 8.86 × 10-4 tonne CH 4 ⎛ 0.80 mole CH 4 (facility) ⎞
E CH 4 = × × ×⎜ ⎟
day yr bbl ⎝ 0.788 mole CH 4 (default EF) ⎠
E CH 4 = 1,880 tonnes CH 4 /yr
As noted in Table 5-1, there is no uncertainty associated with the count of equipment. The uncertainty
associated with the volume of gas is 5% and the uncertainty associated with the days of operation is 2%.
The uncertainty in the emissions factor is 90% based on expert judgment and the uncertainty in the site
specific CH4 mole % is 4 % based on Exhibit 5-1. The uncertainty in the CH4 content is 5.53% from the
GRI/EPA Methane Study (Shires and Harrison), which is provided in Table E-4 of the API Compendium.
The uncertainty of the emissions is calculated by applying Equation 4-6 and using the relative uncertainty
values.
6,100 bbl 343 day 8.86 × 10-4 tonne CH 4 tonne mole CH 4 0.80 mole CH 4 (facility)
E CO2 = × × × ×
day yr bbl 16.04 tonne CH 4 0.788 mole CH 4 (default EF)
tonne mole gas 0.12 tonne mole CO 2 44.01 tonne CO 2
× × × = 775 tonnes CO 2 /yr
0.80 tonne mole CH 4 tonne mole gas tonne mole CO 2
The same approach applies to estimating the uncertainty for the CO2 emissions. The uncertainty associated
with the mole% of CO2 is 4%, so the results are the same as for CH4.
U (abs ) ECO e = (0.775 × 775) 2 + (0.775 × 21×1,880) 2 = 35, 700 tonnes CO 2 e/yr
2
F.8 Uncertainty in Vented Sources – Acid Gas Removal (Amine Unit Emissions)
Table F-10 summarizes the operating parameters for acid gas removal at the example facility.
Annual
Average Activity
Operation Factor (all
Unit Capacity Uncertainty (per unit per Uncertainty units Uncertainty
Source Fuel (per unit) (±%)a year) (±%)a combined) (±%)a
6
Inlet Produced Gas 30×10 ±5% scf/day 343 days/yr ±2% days/yr 10,290×106 ±5.39%
scf/day scf/yr scf/yr
Outlet Outlet Gas N/A N/A 8,997×106 ±5% scf/yr
(0.5 mole% CO2) scf/yr
Methane emissions from the amine unit are estimated based on the emission factor provided in the API
Compendium and the quantity of gas processed.
30 × 106 scf 343 day 0.0185 tonne CH 4 0.80 tonne mole CH 4 (facility)
E CH4 = × × ×
day processed yr 106scf 0.788 tonne mole CH 4 (default)
E CH4 = 193 tonnes CH 4 /yr
We assume there is no uncertainty in the number of units. The uncertainty in the daily volume of gas
processed is 5% and the uncertainty in the days of operation is 2% based on expert judgment. The
uncertainty in the emission factor is 119% based on Table 5-5 of the API Compendium. Because it is not
possible to have negative emissions, the lower bound uncertainty will be truncated at zero. As a result, the
lower and upper bound uncertainties are estimated separately.
The uncertainty in the facility mole% is 4%, the uncertainty in the default mole% is 5.53% from the
GRI/EPA Methane Study (Shires and Harrison, 1996), which is also provided in Table E-4 of the API
Upper:
U(rel ) 2Number of Units +U(rel ) Daily
2 2 2 2
Volume +U(rel ) Days of Operation +U(rel ) EF +U(rel ) Facility mole%
U (rel )emissions =
+U(rel ) 2Default mole%
Lower:
U(rel ) 2Number of Units +U(rel ) Daily
2 2 2 2
Volume +U(rel ) Days of Operation +U(rel ) EF +U(rel ) Facility mole%
U (rel )emissions =
+U(rel ) 2Default mole%
To calculate the uncertainty for the CO2 emissions, first calculate the uncertainty of the sour and sweet gas
components and then aggregate those uncertainties. For the uncertainty in the sour gas flow rate we assume
there is no uncertainty in the number of units, the uncertainty in the daily volume of gas processed is 5%,
and the uncertainty in the days of operation is 2% based on expert judgment. The uncertainty in the sweet
gas flow rate scf/yr is 5% based on expert judgment. The uncertainty in the CO2 concentrations for both the
sour and sweet gases are assumed to be 4%. The uncertainties of the first and second terms of the emissions
equation are calculated by applying Equation 4-6 and using the relative uncertainty values.
The aggregated uncertainty of the emissions is then calculated by applying Equation 4-4 and using the
absolute uncertainty values.
Upper: U (abs ) = (0.0697 × 62, 600) + (1.19 × 21× 190) 2 = 6,510 tonnes CO 2e/yr
2
E CO2e
F.9 Uncertainty in Vented Sources – Pneumatics Devices and Chemical Injection Pumps
Table F-11 summarizes the operating parameters for acid gas removal at the example facility.
Table F-11. Operating Parameters for Pneumatic Devices and Chemical Injection Pumps
Uncertainty Estimate for CH4 and CO2 Emissions from Pneumatic Devices:
A 5% uncertainty was assigned to the count of pneumatic devices based on expert judgment. The
uncertainty in the emission factor is 49.5% per Table 5-15 of the API Compendium. The uncertainty in the
default methane content is 5.53% from the GRI/EPA Methane Study (Shires and Harrison), which is also
provided in Table E-4 of the API Compendium. The uncertainty in the site specific mole % is 4%. The
uncertainty of both the CO2 and CH4 emissions are calculated by applying Equation 4-6 and using the
relative uncertainty values.
U (rel ) ECO ,E CH4 = U(rel ) 2AF +U(rel ) 2EF +U(rel )3Facility mole% +U(rel )3Default mole%
2
The uncertainty in the default methane content is 5.53% from the GRI/EPA Methane Study (Shires and
Harrison), which is provided in Table E-4 of the API Compendium. The uncertainty in the site specific
mole % is 4%. The uncertainty of both the CO2 and CH4 emissions are calculated by applying Equation 4-6
and using the relative uncertainty values.
= U(rel ) 2AF +U(rel ) 2EF +U ( rel )Facility mole% +U ( rel )Default mole%
2 2
U (rel ) ECO ,E CH4
Upper: 2
U(rel) ECO & E CH4 = U(rel ) 2AF + U(rel ) 2EF + U(rel ) 2Facility mole% + U(rel )Default
2
mole%
2
Lower:
( −5.00 ) + ( −100 ) + ( −4.00 ) + ( −5.52 ) = −100% ( tonnes/yr )
2 2 2 2
U(rel)ECO & E CH4 =
2
Table F-12 summarizes the equipment parameters for PRVs and other equipment contributing to
maintenance/turnaround emissions at the example facility.
Annual Activity
Factor (all units
Source Fuel combined) Uncertainty (±%)a
Vessel blowdowns (non-routine) Produced Gas 112 vessels ±0% vessels
Compressor starts (non-routine) Produced Gas 11 compressors ±0% compressors
Compressor blowdowns (non-routine) Produced Gas 11 compressors ±0% compressors
Well workovers (non-routine) Produced Gas 24 well workovers ±0% well workovers
Pressure relief valves Produced Gas 482 PRVs ±1% PRVs
Methane emissions from vessel blowdowns, compressor starts, compressor blowdowns, and oil well
workovers are estimated using the emission factors presented in API Compendium Table 5-23. Methane
emissions from PRVs are estimated using an emission factor from API Compendium Table 5-24. Carbon
dioxide emissions are estimated using the ratio of the facility gas CO2 to the default composition.
There is no uncertainty in the count of the major equipment (vessels, compressors, and wells). An
uncertainty of 1% was assigned to the count of PRVs based on expert judgment. The uncertainty in the
emissions factor for vessel blowdowns is 326% based on Table 5-23 of the API Compendium. The
uncertainty in the emissions factor for PRVs is 310% based on Table 5-24 of the Compendium. The
uncertainty in the emissions factor for compressor starts is 190% based on Table 5-23 of the API
Compendium. The uncertainty in the emissions factor of oil well workovers is 300% based on expert
judgment. The uncertainty in the emission factor of compressor blowdowns is 179% based on Table 5-23 of
the API Compendium. Because it is not possible to have negative emissions, the lower bound uncertainties
will be truncated at zero. As a result, the lower and upper bound uncertainties are estimated separately.
The uncertainty in the default methane content is 5.53% based on the GRI/EPA Methane Study (Shires and
Harrison, 1996), which is also provided in Table E-4 of the API Compendium. The uncertainty in the site
Uncertainty Estimate for CH4 and CO2 Emissions from Vessel Blowdowns:
U (rel ) ECO ,E CH4 = U(rel ) 2AF +U(rel ) 2EF +U(rel ) 2Facility mole% +U(rel ) 2Default mole%
2
Because the emissions cannot be less than zero, the lower bound uncertainty value for this emission source
is truncated at 100%. Therefore, separate upper and lower bound uncertainties are estimated for the CO2e
emissions.
Upper: U (abs ) = (3.26 × 0.0702) + (3.26 × 21× 0.171) 2 = 11.7 tonnes CO 2 e/yr
2
ECO2e
U (rel ) ECO , ECH 4 = U ( rel ) 2AF + U (rel ) 2EF + U (rel ) 2Facility mole % + U Default
2
mole %
2
Because the emissions cannot be less than zero, the lower bound uncertainty value for this emission source
is truncated at 100%. Therefore, separate upper and lower bound uncertainties are estimated for the CO2e
emissions.
Upper: U (abs) = (1.90 × 0.745) + (1.90 × 21×1.81)2 = 72.3 tonnes CO2 e/yr
2
ECO2e
U (rel ) ECO ,E CH4 = U(rel ) 2AF +U(rel ) 2EF +U(rel ) 2Facility mole% +U(rel ) 2Default mole%
2
Because the emissions cannot be less than zero, the lower bound uncertainty value for this emission source
is truncated at 100%. Therefore, separate upper and lower bound uncertainties are estimated for the CO2e
emissions.
Upper: U (abs ) = (1.79 × 0.333) 2 + (1.79 × 21× 0.808) 2 = 30.4 tonnes CO 2 e/yr
ECO2e
Uncertainty Estimate for CH4 and CO2 Emissions from Oil Well Workovers:
U (rel ) ECO ,E CH4 = U(rel ) 2AF +U(rel ) 2EF +U(rel ) 2Facility mole% +U(rel ) 2Default mole%
2
Because the emissions cannot be less than zero, the lower bound uncertainty value for this emission source
is truncated at 100%. Therefore, separate upper and lower bound uncertainties are estimated for the CO2e
emissions.
Upper: U (abs) = (3.00 × 0.0181)2 + (3.00 × 21× 0.0439)2 = 2.77 tonnes CO2e/yr
ECO2e
0.921tonnes CO 2 e/yr
U(rel ) ECO e = 100% × = −98.1% ( tonnes CO 2 e/yr )
2
0.939 tonnes CO 2 e/yr
Uncertainty Estimate for CH4 and CO2 Emissions from Pressure Relief Valves:
U (rel ) ECO ,E CH4 = 1.002 + 3102 + 4.002 + 5.532 = +310%, −100% ( tonnes/yr )
2
Because the emissions cannot be less than zero, the lower bound uncertainty value for this emission source
is truncated at 100%. Therefore, separate upper and lower bound uncertainties are estimated for the CO2e
emissions.
Upper: U (abs) = (3.10 × 0.131) + (3.10 × 21× 0.318)2 = 20.7 tonnes CO2e/yr
2
ECO2e
Table F-13 provides fugitive component counts associated with the high CO2 content onshore oil field
facility. The corresponding average component emission factors are also provided based on the API
average oil and gas production emission factors given in Table 6-14 of the API Compendium for light crude.
Because the emission factors are taken from API 4615, the weight fraction of CH4 is taken from the
“Generic” weight fraction for light crude give in the API Compendium, Table C-6, which assigns
composition by facility type for an aggregate of all streams associated with that facility. Note that the API
Compendium Table C-6 does not provide CO2 speciation data; therefore, these emissions are not calculated.
Table F-13. Onshore Oil Field (High CO2 Content) Fugitive Emission Factors
Average Uncertainty b Component EF, tonnes Uncertainty b
Component Service Component Count % TOC/comp./hr a %
Valves Liquid and Gas 2,740 75.0 1.32E-06 100
Pump Seals Liquid and Gas 185 75.0 3.18E-07 100
Connectors Gas 110 75.0 1.64E-07 100
Flanges Liquid and Gas 10,000 75.0 7.69E-08 100
OELs Gas 6 75.0 1.21E-06 100
Others Liquid and Gas 710 75.0 7.50E-06 100
Footnotes:
a
Note that for this example, the TOC weight fraction of the gas stream is not 100%. However, the TOC emissions are not adjusted here since
such an adjustment cancels out when calculating CH4 emissions as shown in Equation 6-9 of the API Compendium.
b
Uncertainty based on engineering judgment at a 95% confidence interval.
The uncertainty in the number of devices is 75.0%, the uncertainty in the emission factors is 100%, and the
uncertainty in the weight fraction of CH4 is 15% based on expert judgment. The uncertainty of the
emissions estimate is calculated by applying Equation 4-6 and using the relative uncertainty values. The
calculation is demonstrated for the liquid and gas pump seals.
U (rel ) ECH = U (rel ) 2Devices + U (rel ) 2EF + U (rel )2hr / yr + U (rel )Weight
2
%
4
Because the emissions cannot be less than zero, the lower bound uncertainty value for this emission source
is truncated at 100%. Therefore, separate upper and lower bound uncertainties are estimated for the CO2e
emissions.
The uncertainty of total fugitive emission estimate is calculated by applying Equation 4-4 using the absolute
uncertainties. This is demonstrated for the upper bound uncertainty of CH4 below.
CH 4 : U (abs )
∑ ( Fugitive )
= ∑U (abs) 2
Fugitive
Emissions associated with the air conditioning units in fleet vehicles for this example are calculated based
assuming the refrigerant is R-134a. As stated in Section F.5, the uncertainty associated with the count of
vehicles is 0%. The uncertainties associated with the refrigerant capacity and the emission factor are
assumed to be 100% and 50%, respectively. The calculations for fleet vehicle refrigeration emissions are
shown below.
The global warming potential for R-134a is treated as a constant, so the uncertainty in the CO2e emissions is
the same as for the refrigerant emissions.
Emissions associated with the electricity purchased by the facility are calculated using emission factors in
the API Compendium.
Based on expert judgment, the uncertainty in the activity factor is 2%, the uncertainty in the CO2 emission
factor is 10.0%, and the uncertainty in the CH4 and N2O emission factors is 100%. The uncertainty of the
emissions estimate is calculated by applying Equation 4-6 and using the relative uncertainty values.
U (rel )CO2 = U (rel ) 2AF + U (rel ) 2EF = 22 + 102 = ±10.2% ( tonnes CO 2 /yr )
U (rel )CH 4 and N2O = U (rel ) 2AF + U (rel ) 2EF = 22 + 1002 = ±100% ( tonnes/yr )
The final step is to convert the emissions of CO2, CH4, and N2O to a CO2e. Note there are no uncertainties
for any of the GWPs, so the combined uncertainty is based on summing the individual emissions and
applying Equation 4-4 (using absolute uncertainty values).
U (abs ) ECO e = (0.102 × 551) 2 + (1.00 × 21× 0.00776) 2 + (1.00 × 310 × 0.00628)2 = 56.2 tonnes CO 2 e/yr
2
Total emissions for this facility are summarized in Table F-15, along with the summary uncertainty
propagation. The uncertainty of the total emissions is calculated by applying Equation 4-4 and using the
absolute uncertainty values. In practice, uncertainty values are usually displayed in relative terms, which are
provided in terms of upper and lower % uncertainty values in Table F-15 for those emissions sources with
an asymmetrical uncertainty.
As an example, the summed upper bound uncertainties for N2O combustion emissions for the facility are
calculated as follows, using absolute uncertainties:
This calculation was repeated for the sum of each category, and for the totals in that GHG gas type. The
gases were each converted to CO2e using the GWPs. As was stated previously, for the purpose of compiling
a facility-level GHG inventory, the GWP’s are assumed to have no uncertainty. The bottom row of
Table F-15 shows the emissions in CO2e.