Methane (CH4) leakage from natural gas (NG) pipelines poses an environmental, safety, and economic threat to the public. While previous leak detection and quantification studies focus on the aboveground infrastructure, the analysis of underground NG pipeline leak scenarios is scarce. Furthermore, no data from controlled release experiments have been published on the accuracy of methods used to (1) quantify emissions from an area source and (2) use these emissions to quantify the size of a subsurface leak. This proof-of-concept work uses CH4 mole fraction, as measured by a single gas sensor, as an input to a simple dispersion-based model (WindTrax) under ideal conditions (i.e., in a field) and compares the calculated emissions to the known controlled NG release rates. The aboveground and surface CH4 mole fractions were measured for 5 days at a field testbed using controlled underground release rates ranging from 0.08 to 0.52 kg hr–1 (3.83–24.94 ft3 hr–1). Results confirmed that the mean normalized CH4 mole fraction increases as the atmosphere transitions from the Pasquill–Gifford (PG) stability class A (extremely unstable) to G (extremely stable). The estimated surface CH4 emissions showed large temporal variability, and for the emission rates tested, at least 6 h of data are needed to have a representative estimate from subsurface pipeline leaks (±27% of the controlled release rate on average). The probability that the emission estimate is within ±50% of the controlled release rate (P±50%) is approximately 50% when 1 h of data is collected; the probability approaches 100% with 3–4 h of data. Findings demonstrate the importance of providing enough data over time for accurate estimation of belowground leak scenarios. By adopting the estimation method described in this study, operators can better estimate leakage rates and identify and repair the largest leaks, thereby optimizing annual greenhouse gas emissions reductions and improving public safety.
Fugitive methane (CH4) emissions from natural gas (NG) pipelines pose an ongoing threat to the environment and public safety. Common causes of pipeline leakage include corrosion, incorrect operation due to human error, excavation or similar ground disturbances, natural force damage caused by ground shifts or aboveground loads, and equipment failure (Pipeline and Hazardous Materials Safety Administration [PHMSA], 2018). Pipeline leakage events result in significant safety issues. From 2000 to 2019, the PHMSA reported pipeline incidents resulting in 202 fatalities, 918 injuries, and an estimated total cost of $2.94 billion (PHMSA, 2022). It is estimated that CH4 emissions from NG distribution main leaks are 0.69 Tg per year, accounting for 7.6% of the total U.S. CH4 emissions (Weller et al., 2020). CH4 emissions from the NG distribution networks contribute to 56% of the total estimated CH4 emissions in Paris (Defratyka et al., 2021). Given that the global warming potential of CH4 is 84 times greater than that of CO2 over a 20-year horizon (Intergovernmental Panel on Climate Change, 2014), the environmental benefits of NG as a short-term approach for mitigating climate change would be negated if CH4 emissions in the NG supply chain exceed 3.2% of total NG production (Alvarez et al., 2012; United Nations Framework Convention on Climate Change, 2022). Typically, local NG distribution companies detect atmospheric CH4 mole fraction from underground leaks using vehicle and walking surveys (PHMSA, 2017). However, studies showed that NG pipeline leaks are difficult to locate and quantify from the aboveground observations due to limited understanding of (1) how CH4 accumulates or disperses in the atmosphere under different environmental conditions and (2) the underground CH4 transport and surface manifestation of such leaks (Okamoto and Gomi, 2011; Mitton, 2018; National Association of Regulatory Utility Commissioners, 2019; Cho et al., 2020). Therefore, a better understanding of the influence of atmospheric conditions on the aboveground CH4 plume from underground leaks is required to improve leak detection and quantification and risk assessment.
NG distribution pipelines are surveyed by distribution companies through routine inspection. Once a leak is discovered, it is classified based on the existing or probable hazard to the person and property (Hendrick et al., 2016; Washington Administrative Code, 2019). Generally, Grade 1 leaks are prioritized for immediate repair. Grade 2 and 3 leaks are classified as nonhazardous based on location, CH4 magnitude, and potential for migration of the leaked gas. Grade 2 and 3 leaks require repair within 12–15 months of detection and then revaluation after 15 months of discovery. Before repair, unresolved Grade 2 and 3 leaks require periodic inspection.
Despite repeated observation, leak classification does not consider the effects of environmental conditions on the magnitude of the measured CH4 mole fraction. A recent vehicle survey in Durham, NC, showed that leak detection density was 50% higher at night (0.33 leaks per mile) compared to daytime conditions (0.22 leaks per mile) when following the same route (Gallagher et al., 2015). Therefore, many underground leaks may go undetected using current survey criteria that fail to incorporate the effects of environmental conditions on leak detection.
Previous studies have used several methods to estimate emissions from underground NG pipeline leaks, including bottom-up, local level methods, such as the surface dynamic chamber method (Lamb et al., 2015; Hendrick et al., 2016; Zimmerle et al., 2017), downwind tracer-ratio methods (Lamb et al., 2015; Weller et al., 2018), and vehicle-based statistical methods (Von Fischer et al., 2017; Weller et al., 2018; Weller et al., 2020). Enclosing a surface area to measure the CH4 enhancement, the surface dynamic chamber method relies on knowledge of the flow of air through a controlled chamber volume (National Academies of Sciences, 2018; Riddick et al., 2020). The surface enclosure area is typically 1 m2 in size with measurement times in the range of minutes (National Academies of Sciences, 2018) or a series of gridded enclosure placements to overcome the spatial limitations (Lamb et al., 2015). The main shortcoming of this method is that the surface CH4 expression from underground pipeline leaks has been documented to migrate 2–6 m from the leak location, depending upon the leak size and environmental conditions, such as soil properties and moisture (Okamoto and Gomi, 2011; Ulrich et al., 2019; Cho et al., 2020; Cho et al., 2022), and the chamber may not capture all emitted gas. In addition, recent work by Riddick et al. (2021) documents the temporal variability in the surface CH4 expression above an NG pipeline leak. The authors identify variations in wind speed, roughness length, and atmospheric stability as the key environmental variables most significantly affecting the surface CH4 mole fraction. Cho et al. (2022) measured this change continuously over 4 days for 2 different known leakage rates, demonstrating how it varies by over 100%, depending on the time of day. Downwind tracer-ratio methods are often used as a supplementary quality assurance method to evaluate the performance of other methods, such as the surface chamber method (Lamb et al., 2015). Such methods require accurate knowledge of the known leak location and are sensitive to environmental conditions, such as wind and the accurate co-placement of the tracer release and leak site (Mønster et al., 2019). Vehicle-based methods estimate emissions using a statistically fitted linear model with the maximum CH4 mole fraction, the total plume peak area, and the ratio of maximum CH4 mole fraction to peak area (Von Fischer et al., 2017). Although this method shows promise for city or regional studies, it relies on statistical analysis and does not incorporate the impacts of atmospheric conditions on the dispersion of the CH4 plume. When compared with surface chamber estimates, vehicle-based statistical estimates show larger emission estimates with very few indications of ground truth estimates from known underground leak rates (Weller et al., 2018).
As few of the above approaches have been verified using controlled underground releases, evaluating their accuracy remains an open challenge, and therefore, CH4 emissions from pipeline systems remain poorly characterized. To the authors’ knowledge, there have been no studies to date that investigate the effects of time-varied weather conditions on emission estimate accuracy from subsurface pipeline leaks in controlled conditions. Comparing a method’s emission estimates to known, controlled underground release over a range of environmental and operating conditions is the first step toward method validation.
Atmospheric dispersion models have been used to inversely calculate gas emissions across a wide range of scales, spanning from the microscale to the continental scale (Flesch et al., 2004; Jeong, 2011; Bonifacio et al., 2013; Barkley et al., 2017). The application considered here relies on simple downwind measurements and could be readily applied to open field environments free from obstructions, such as buildings, parked vehicles, and trees. Such conditions are typically seen in upstream sectors, such as gathering and flowlines or in the downstream sector in rural environments or along open extents. For these scenarios that represent short-range dispersion (1 km or less), the inverse dispersion model WindTrax has reported an accuracy of ±10% for emission estimates from agricultural sources (i.e., feedlots, pastures, and farms), landfills, and natural sources (Laubach and Kelliher, 2005; Yang et al., 2017; Riddick et al., 2018). As this method is not limited to specific source sizes and shapes, it allows for flexibility in sensor placement and therefore has good potential for use in estimating CH4 emissions from NG pipeline leakage scenarios (Riddick et al., 2019).
This study explores the effectiveness of using limited surface and aboveground CH4 mole fraction data as inputs to a near-field dispersion model to calculate CH4 emissions from subsurface pipeline leaks under realistic conditions. Explicitly, this study aims to (1) link the variations in atmospheric CH4 mole fraction magnitude to variations in atmospheric stability under a range of controlled NG leak rates, (2) test the accuracy of an inverse modeling approach to estimate emissions by comparing estimates directly against known, controlled underground leak rates, and (3) understand and quantify how much data are needed to estimate emissions within ±50% of controlled leak rates. To our knowledge, this is the first study of its kind to test the emission estimate accuracy against the known rates for belowground pipeline scenarios. Currently, pipeline operators use surface CH4 measurements to infer the size of a leak and prioritize repair. By adopting the method described in this study to correctly estimate leak rate, operators can estimate emissions and therefore optimize repair and greenhouse gas (GHG) emission reductions.
2. Materials and methods
2.1. Experimental testbed/setup
Field experiments were conducted at the Methane Emission Technology Evaluation Center (METEC) field site, Colorado State University, Colorado. The METEC site includes an underground pipeline testbed that simulates the behavior of pipeline leaks using controlled NG releases under varied subsurface, surface, and atmospheric conditions. Details of the testbed design can be found in Mitton (2018), Ulrich et al. (2019), Cho et al. (2020), and Gao et al. (2021) and will be summarized here. For this work, a testbed simulating a rural right of way with grassy surface cover was selected. The testbed consists of stainless-steel tubing (0.635 cm, model SS-T4-S-035-20, Swagelok, Solon, OH, USA) installed at 0.9 m below the ground surface along the east-west direction of the site, configured to release compressed NG at specific rates and locations. The tubing runs adjacent to a 10-cm diameter polyvinylchloride pipe to simulate a realistic underground pipeline. The pipeline is backfilled with native soil, classified previously as sandy loam, according to the U.S. Department of Agriculture Soil Texture Classification. The trench is slightly less compact than the surrounding undisturbed soil as discussed in Ulrich et al. (2019). NG (87 ± 2 vol % CH4) was injected through the stainless-steel tubing into the testbed from an aboveground 145-L compressed NG tank. The gas release rate (reported at kg hr–1 of CH4) was controlled by pressure regulators, solenoid valves, and choked flow orifices and measured using a thermal mass flow meter (range: 0–0.66 kg hr–1, accuracy: ±1%, Omega FMA1700 series). Gas was released to the soil through a stainless-steel mud dauber (0.635 cm, model SS-MD-4, Swagelok, USA) with a 10-cm wire mesh cube filled with gravel to prevent clogging.
CH4 mole fractions were collected at two elevations, 1 m above the ground and at the ground surface. The aboveground CH4 mole fractions were measured at two separate locations, depending on the experiment (Figure 1a). For the first three experiments (June 1 8:00 AM to June 4 8:00 AM), a trace gas analyzer (Picarro Inc., Santa Clara, CA, USA, G2210, precision: <0.1 ppb) was positioned at location T1 to measure the CH4 mole fraction at a distance of 10 m from the leak center and 1 m above the ground surface. The measurement interval was approximately 0.8–1 Hz. The analyzer was moved to location T2, located 10 m from the leak center, for the last experiment (June 4 8:00 AM to June 5 8:00 AM) with the same measurement height as the predominant wind direction changed during this period.
Surface CH4 was measured at the ground surface at 80 locations (Figure 1b) to determine the plume shape and size using a handheld CH4 detector (Bascom Turner Instruments Inc., Norwood, MA, USA, VGI-201 Gas-Rover, resolution: 1.0 ppm, detection limit: approximately 5.0 ppm). The surface CH4 mole fraction was measured during the last 1 h of each experiment following the method used in a previous study (Cho et al., 2022). In accordance with the manufacturer’s procedures, at each measurement point (Figure 1b), the rubber cup of the detector was lightly pressed onto the ground surface for 30 s. CH4 mole fraction was averaged over this time. To ensure that most of the surface expression was recorded, measurements were carried out from the leak center in each direction, A1–A8, until the instrument’s detection limit (approximately 5.0 ppm) was reached.
Two meteorological sensors were deployed to continuously measure wind speed, wind direction, solar radiation, and air temperature with a time resolution of 1 min (Figure S1). The first meteorological sensor (METER Group Inc., Pullman, WA, USA, ATMOS41, accuracy: 3% and ±5° for wind speed and direction, respectively; solar radiation: ±5%; air temperature: ± 0.6°C; air pressure: ±0.1 kPa from –10°C to 50°C, ± 0.5 kPa from –40°C to 60°C) was co-placed with the trace gas analyzer at T1 (1 m height). The second sensor (YOUNG Model 81000 Ultrasonic Anemometer, accuracy: wind speed: ±1% ±0.05 m s–1 [0–30 m/s] and ±3% [30–40 m s–1]; wind direction: ±2° [1–30 m s–1] and ±5° [30–40 m s–1]) was a fixed station installed at 6 m above ground, located 76.2 m northwest of the leak center.
2.2. Experiment procedure
Four 24-h experiments were conducted with controlled release rates of 0.08, 0.18, 0.27, and 0.52 kg hr–1 of CH4 (3.83, 8.64, 12.94, and 24.94 standard cubic feet per hour), representing medium to large distribution pipeline leaks (Lamb et al., 2015). Experiments were conducted between of June 1 and 5, 2020.
A 0.08-kg hr–1 leak was initiated on the first day at approximately 8:00 AM mountain daylight time, and continuous measurements were conducted over the following 24 h, at which time, the leak rate was increased to 0.18 kg hr–1. This process/timing was repeated for 0.27 and 0.52 kg hr–1. Mitton (2018) showed that NG released in the METEC testbed took approximately 4 h to migrate to the surface. Therefore, only the last 20 h of data (12:00 PM–8:00 AM) were used for each of the four experiments in the analysis to allow equilibration of the subsurface and ensure steady-state conditions. To properly represent atmospheric processes at the microscale, 10-min averaged data were calculated. For every 10 min, atmospheric stability was quantified as the Obukhov length (L, m) to express the relative roles of shear and buoyancy (Foken, 2006). The Obukhov length was calculated from the surface friction velocity (u*, m s–1), the mean absolute air temperature (T, K), the von Kármán constant (kv = 0.41), the gravitational acceleration (g = 9.81 m s–2), and the 3-D horizontal and vertical wind vectors (u, v, w, m s–1; Equations 1 and 2; Flesch et al., 2004; Foken, 2006).
The Obukhov length was converted to the PG stability class to provide practical insight for industry in situations, where on-site 3D sonic anemometer data are not available (Ulrich et al., 2019; Riddick et al., 2021). In this work, –100 ≤ L < 0, –500 ≤ L < –100, |L| > 500, 50 ≤ L < 100, and 0 < L < 50 are considered as extremely unstable (PG: A), unstable (PG: B and C), neutral (PG: D), stable (PG: E and F), and extremely stable conditions (PG: G), respectively (Gryning et al., 2007).
2.3. Model setup
The surface CH4 emissions were inversely calculated using WindTrax version 2.0 (Crenna, 2020). WindTrax is a publicly available backward Lagrangian stochastic-based model that can simulate the trajectories of thousands of tracer particles in the atmosphere surface layer (Flesch et al., 1995). The surface CH4 emissions Q(kg hr–1) were calculated inversely using Equation 3 by assuming the simulated normalized CH4 mole fraction (C/Q)simequals the actual normalized CH4 mole fraction (Flesch et al., 2004)
where Cobs is the measured CH4 mole fraction (ppm) at the downwind location, and Cb is the measured local CH4 background mole fraction (ppm). Additional details about the calculation of (C/Q)simin this model can be found in Flesch et al. (2004).
Modeling inputs included CH4 concentration, 10-min averaged wind speed and direction, air temperature, roughness length, surface expression area, and Obukhov length. The surface roughness length was estimated at 0.02 m, approximated as one-tenth of grass height (Gutmann and Small, 2005). The surface expression area was approximately 7–20 m2 (Table 1), which was estimated using the measured surface CH4 mole fractions (Figure S2). Atmospheric CH4 mole fractions were filtered for two criteria: (1) the measurement was downwind of the known leak source (±90° of leak center) and (2) the CH4 mole fraction was higher than the background CH4 mole fraction. The background CH4 mole fraction ranged from 1.88 to 1.92 ppm (Table 1), which was calculated by averaging the mole fraction within 25% of the cumulative distribution frequency (Figure S3).
|Experiment .||Background CH4 Mole Fraction at 1 m (ppm) .||Surface Expression Area (m2) .|
|Experiment .||Background CH4 Mole Fraction at 1 m (ppm) .||Surface Expression Area (m2) .|
Standard deviations from the background CH4 mole fraction are noted in parentheses.
Sensitivity analysis on seven stability cases (Srun1-7) was used to evaluate the effects of atmospheric stability on the surface CH4 emission estimate by changing L from –75% to 100% relative to the base case (Baserun; Table 2). The base case was the controlled case in WindTrax to calculate the surface emission using 10-min averaged CH4 mole fraction and associated meteorological variables as the model input parameters.
|Baserun .||Srun1 .||Srun2 .||Srun3 .||Srun4 .||Srun5 .||Srun6 .||Srun7 .|
|L||L × (1 – 75%)||L × (1 – 50%)||L × (1 – 25%)||L × (1 + 25%)||L × (1 + 50%)||L × (1 + 75%)||L × (1 + 100%)|
|Baserun .||Srun1 .||Srun2 .||Srun3 .||Srun4 .||Srun5 .||Srun6 .||Srun7 .|
|L||L × (1 – 75%)||L × (1 – 50%)||L × (1 – 25%)||L × (1 + 25%)||L × (1 + 50%)||L × (1 + 75%)||L × (1 + 100%)|
L in the base case (Baserun) was calculated using Equations 1 and 2. The Obukhov length for the 7 sensitivity cases (Srun1–7) was −75% to 100% of L, respectively (L ≠ 0 in WindTrax; Flesch et al., 2004); other input model parameters were the same as the Baserun.
3.1. Impact of atmospheric stability on CH4 mole fraction
Figure 2 shows the mean normalized CH4 mole fraction as a function of the PG stability class for the four experiments (Experiments 1–4). The mean normalized CH4 mole fractions consistently increase with the increase in atmospheric stability, ranging from 12.31 to 23.41s m–3 when the atmosphere transitions from the PG stability classes A–G. The large variability in mean normalized CH4 mole fraction may result from only measuring CH4 mole fraction at one location within the plume for each measurement interval as well as the variations in CH4 mole fraction within the plume for the same stability condition. The normalized CH4 mole fraction at the plume center was larger than that at the plume edge for the same stability condition, and the normalized mole fraction near the plume centerline during the unstable condition (e.g., PG stability class A) was larger than that at the edge of the plume during the stable condition (e.g., PG stability classes E and F) (Figure S4). In all, the results confirmed that the mean normalized CH4 mole fraction increases as the atmosphere transitions from PG stability A–G. For the conditions tested, our results show that measured CH4 mole fraction cannot be used to infer the size of a leak without considering the impacts of atmospheric stability conditions.
3.2. Methane emissions calculated by WindTrax
As seen in Figure 3a, the model estimated surface CH4 emissions showed large temporal variability, with the maximum estimates 3–5 times higher than the observed release rates. The absolute percentage difference between the estimated surface CH4 emissions and the controlled release rate ranged from 1% to 48% for the four experiments (27% on average). The large absolute percentage difference in Experiment 4 may be due to the lower sample size, compared to Experiments 1–3. While Experiments 1–3 had a sample size of 6–9 h, only 3.5 h were available for Experiment 4. It is possible that the shortened time (3.5 h) is too short to accurately represent the highly temporally variable surface CH4 emissions from a subsurface leak. Thus, the results demonstrate that for the emission rates tested, the average surface CH4 emissions from WindTrax can provide promising estimates (within 27% on average) of the leak rate from an underground NG pipeline release based on 6.5 h of downwind measurements from a single sensor location. The average ratio of the calculated surface CH4 emissions by WindTrax to the controlled release rate tends to decrease with the increase in atmospheric stability except for PG stability class G, that is, a larger emission estimate during unstable conditions than in stable conditions (Figure 3b). The larger average ratio in PG stability class G than E and F suggests that other factors not considered in the atmospheric model, such as soil conditions, impact the surface CH4 emissions (Cho et al., 2020; Gao et al., 2021). Previous studies have shown that temporal variability of surface CH4 emissions is associated with both soil and atmospheric conditions, such as soil moisture and fluctuations in atmospheric pressure (Czepiel et al., 2003; Poulsen and Møldrup, 2006; Forde et al., 2019; Kim et al., 2020).
3.3. Effect of atmospheric stability on CH4 emission estimates
Seven sensitivity cases were further used to evaluate the effect of atmospheric stability on the calculated emission rate by altering the Obukhov length in the model from –75% to 100% of L (Table 2). As shown in Figure 4a, the ratio of estimated surface CH4 emissions to the controlled release rate decreases with the increase in atmospheric stability, where the PG stability class changes from PG stability class A (extremely unstable) to G (extremely stable). As the atmosphere transitions from A to G, the average dispersion coefficients of the along-wind, crosswind, and vertical direction decrease rapidly (Figure 4b). The dispersion coefficient in the 3-D direction is about 0.4–0.8 m s–1 when the atmosphere is in PG stability class A (extremely unstable), but it reduces to below 0.2 in PG stability class G (extremely stable, a 50%–70% decrease). Thus, the strong dispersion during unstable atmospheric conditions results in enhanced vertical and horizontal mixing compared to the inhibited dispersion and hence mixing during stable conditions (Riddle et al., 2004; Leelőssy et al., 2014). This results in broad yet diluted CH4 plume during unstable conditions (e.g., PG stability class A) and narrow yet concentrated CH4 plume during stable conditions (e.g., PG stability classes E and F; Figure S4).
3.4. Effect of time averaging on CH4 emission estimates
Using the estimated emissions from the base case (Baserun), a rolling mean method was used to understand the amount of downwind data needed over time to accurately estimate the emission rate. Although our experiments collected data for 24 h, in field scenarios, this is often not possible. We therefore wanted to understand the effect of time averaging on the emission estimate. Figure 5 shows the probability that the emission estimate is within ±50% of the controlled release rate (P±50%) for data averaging time of 0.33–9.33 h. As expected, the probability increases with the increase in the rolling mean time in Experiments 1 through 3 (Figure 5). A probability of 50% is achieved when approximately 1 h of data is collected. The probability increases to 80% with about 2 h of data and reaches 100% with 3–4 h of data. However, the probability was not observed to increase with an increase in time in Experiment 4 (Figure 5). This may be due to the smaller sample size (3.5 h of downwind data in Experiment 4) and its distribution characteristics. As shown in Figure 3a, the distribution was normal for Experiment 4 but was right-skewed for Experiments 1–3. Despite this, it should be acknowledged that the low sample size for Experiment 4 may have affected the probability.
This work demonstrates the viability of linking a simple publicly available dispersion model with limited aboveground and surface CH4 mole fraction and micrometeorological data to estimate CH4 emissions from underground pipeline leaks. Previous downwind methods used for aboveground NG infrastructure leaks show that emission estimate accuracy is about ±70% when calculating based on 15–20 min of downwind stationary measurements (Heltzel et al., 2020). For belowground CH4 leaks, our method shows that downwind measurements measured over 4 h can be used to estimate a leak between 0.08 and 0.52 kg h–1 to within ±50% in an open field condition. Although the average CH4 emissions calculated by WindTrax showed good agreement with the controlled release rate, we also recognize that it might not be feasible to directly apply this method to estimate emissions from a large number of pipeline leaks (U.S. Energy Information Administration, 2021). However, this method could be modified to vehicle-based surveys using a gas sensor and associated 3-D anemometer data to extend its applicability. The method presented here is similar to the survey protocols outlined in other mobile methods for NG aboveground infrastructure leaks (Brantley et al., 2014; Caulton et al., 2018; Edie et al., 2020), that is, mobile monitoring and downwind stationary CH4 mole fraction measurements, yet demonstrates the importance of providing enough data over time for proper estimation of belowground leak scenarios. The length of time needed for mobile methods to achieve an equivalent accuracy remains unknown as there are currently no controlled belowground experiments results for mobile methods. The reported large emission estimate uncertainty (between –50% and +350%) of mobile methods on aboveground NG infrastructure leaks suggests more mobile transects are needed to improve the emission estimate (Lan et al., 2015; Rella et al., 2015; Yacovitch et al., 2015). A driving measurement campaign by Caulton et al. (2018) shows that at least 10 repeat transects are needed to reduce uncertainty, where the uncertainty is estimated at 0.05–6.5 times the emission rate for single-transect sites and 0.5–2.7 times the emission rate for sites with more than 10 transects. Therefore, Caulton’s emphasis on more transect time may also apply to the mobile methods for improving emission accuracy from pipeline leaks.
There is limited information regarding emissions from both flow and gathering lines (Zimmerle et al., 2017; Li et al., 2020). The emission factors of NG gathering lines in the current U.S. Environmental Protection Agency Greenhouse Gas Inventory (Hockstad and Hanel, 2018) are based on emission estimates from the distribution pipelines (Campbell et al., 1996). With enhanced mobility, the method presented here shows potential for estimating CH4 emissions from NG flow and gathering lines. Similarly, it is also useful for NG distribution and transmission pipeline operators when an accurate estimate of emission rate is required to inform the priority of leak repairs, such as finding and repairing the largest leaks and prioritizing work schedules.
This work tests an approach to estimate CH4 emission rates compared to known leak rates from subsurface pipeline leaks under realistic conditions. Our results confirm that the mean normalized CH4 mole fraction increases as the atmosphere transitions from PG stability classes A–G. The measured CH4 mole fraction cannot be used to infer the size of a leak without considering the impacts of atmospheric stability conditions. The estimated CH4 surface emissions from a simple publicly available dispersion model showed large temporal variability, and for the emission rates tested, at least 6 h of data are needed to have a representative estimate from subsurface pipeline leaks (±27% of the controlled release rate on average). The probability that the emission estimate is within ±50% of the controlled release rate (P±50%) is about 50% when approximately 1 h of data is collected. Probability increases to 80% with about 2 h of data and reaches 100% with 3–4 h of data. A sensitivity analysis also showed that the estimated surface CH4 emissions from the underground leaks tend to decrease when the atmosphere transitions from the PG stability class A (extremely unstable) to G (extremely stable) due to the inhibited dispersion and mixing in the atmosphere.
These findings of this study emphasize that a certain amount of data are necessary to have an estimate that can reflect the true subsurface leak rates. Atmospheric stability conditions should be incorporated in pipeline emission estimates based on the atmospheric CH4 measurements.
Data accessibility statement
Data sets for this research are available in the in-text data citation reference: Tian, Shanru; Smits, Kathleen M.; Cho, Younki; Riddick, Stuart; Zimmerle, Daniel; Duggan, Aidan. 2022. Replication data for estimating methane emissions from underground natural gas pipelines using an atmospheric dispersion-based method, https://doi.org/10.18738/T8/UAO5XX, Texas Data Repository.
The supplemental files for this article can be found as follows:
Figures S1–S4. (DOC)
Any opinion, findings, conclusions, or recommendations expressed here are those of the authors and do not necessarily reflect the views of those providing technical input or financial support. The trade names mentioned here are merely for identification purposes and do not constitute an endorsement by any entity involved in this study. The authors also thank Timothy Vaughn and Clay Bell for their experimental support.
This work is jointly supported by the U.S. Department of Transportation Pipeline and Hazardous Materials Safety Administration under Grant Nos. 693JK31810013 and 693JK32010011 and by the Colorado Oil and Gas Conservation Commission Mark Martinez &Joey Irwin Memorial Public Projects Fund.
The authors declare no competing financial interest.
Contributed to conception and design: ST, KS, YC, SR, DZ.
Contributed to the acquisition of data: SR, AD.
Contributed to analysis and interpretation of data: ST, KS, YC, SR.
Drafted and/or revised the article: ST, KS, YC, SR, DZ, AD.
Approved the submitted version for publication: ST, KS, YC, SR, DZ, AD.
How to cite this article: Tian, S, Smits, KM, Cho, Y, Riddick, SN, Zimmerle, DJ, Duggan, A. 2022. Estimating methane emissions from underground natural gas pipelines using an atmospheric dispersion-based method. Elementa: Science of the Anthropocene 10(1). DOI: https://doi.org/10.1525/elementa.2022.00045
Domain Editor-in-Chief: Detlev Helmig, Boulder AIR LLC, Boulder, CO, USA
Guest Editor: Brian Lamb, Washington State University, Pullman, WA, USA
Knowledge Domain: Atmospheric Science