IP Library Granted Patent US 12,254,622
Granted Patent B2
US 12,254,622 · App. 18/480,279 · Granted Mar 18, 2025

Computing emission rate from gas density images

Inventors: Andrew J. Speck (Milton, MA); Manasi Doshi (Cambridge, MA); Lukasz Zielinski (Arlington, MA)
Assignee: SCHLUMBERGER TECHNOLOGY CORPORATION
G06T7/001G01M3/16G06T7/62G06T2207/20221G06T2207/30232H04N23/695
View Patent ↗
Loading inventors, assignments & file history…
Monitor This Case
Get email alerts when status or documents change.
Order Certified Copies
Most orders are placed with the USPTO same day — all within 24 business hours.
Order via The Patent Place →
Pre-filled with this patent's details
Quick Facts
Patent No.
US 12,254,622
App. No.
18/480,279
Granted
Mar 18, 2025
Kind
B2
Abstract

Systems and methods are described for calculating an emission rate of a fugitive gas based on a gas density image of the fugitive gas. In an example, a computing device receives a gas density image of a fugitive gas from a camera. The computing device determines how to optimize the fugitive gas in the camera's field of view and instructions the camera to adjust its bearing and zoom accordingly. The camera captures one or more additional images of the fugitive gas, and the computing device stitches the images together where appropriate. The computing device then calculates the emission rate by delineating the fugitive gas in the image and determining a flux of the gas using one of various calculation methods.

Claims (95)

1. A method for determining an emission rate of a fugitive gas, comprising:

receiving, from an imaging device, a first gas density image;

detecting the fugitive gas in the first gas density image;

delineating the fugitive gas in the first gas density image;

sending, to the imaging device, instructions for adjusting a bearing of the imaging device bearing to optimize a position of the fugitive gas in the imaging device's field of view;

sending, to the imaging device, instructions for capturing a second gas density image;

receiving the second gas density image from the imaging device;

calculating, based on the second gas density image, the emission rate of the fugitive gas;

prior to sending the instructions for capturing the second gas density image, determining that a percentage of the fugitive gas is present in less than a threshold percentage of a vertical field of view of the first gas density image; and

in response to the percentage of the fugitive gas being below the threshold percentage, sending, to the imaging device, instructions for modifying a zoom of the imaging device so that the fugitive gas is present in more than the threshold percentage of a vertical field of view of the imaging device.

2. The method of claim 1 , further comprising:

prior to sending the instruction for capturing the second gas density image, determining that the second gas density image includes a first portion of the fugitive gas that does not include all of the fugitive gas;

sending instructions to the imaging device for capturing a third gas density image that includes a second portion of the fugitive gas that is horizontally adjacent the first portion of the fugitive gas;

receiving, from the imaging device, the third gas density image; and

stitching the first gas density image and the second gas density image together.

3. The method of claim 1 , wherein calculating the emission rate of the gas density image comprises:

determining a volume of the fugitive gas based on the second density image and a measured distance between the imaging device and the fugitive gas;

receiving wind measurement data, the wind measurement data including measured wind speed and direction of wind in the vicinity of the fugitive gas;

inputting the wind measurement data and the volume of the fugitive gas into a closed surface integral equation for calculating flux; and

receiving an output from the closed surface integral equation, the output being a calculated amount of the fugitive gas.

4. The method of claim 1 , wherein calculating the emission rate of the gas density image comprises:

creating a cross-section in the fugitive gas in the second gas density image;

dividing the cross-section into a plurality of subsections based on spacing between image pixels in the second gas density image, each of the plurality of subsections having an estimated length;

determining a concentration of pixels for each of the plurality of subsections;

determining a fugitive gas plume travel velocity orthogonal to the cross-section;

inputting the determined concentration of pixels for each of the plurality of subsections, the fugitive gas plume travel velocity, and a weighting function into an equation; and

receiving, as an output from the equation, a flux rate of the fugitive gas for the cross-section.

5. The method of claim 4 , wherein calculating the emission rate of the gas density image further comprises:

calculating the flux rate for additional cross-sections of the fugitive gas;

determining that total flux of the fugitive gas by calculating an average of the calculated flux rate for the cross-section and the additional cross-sections of the fugitive gas.

6. The method of claim 4 , wherein the estimated length of the plurality of subsections is determined by converting pixels to length using degrees of the imaging device's bearing and a measured length between the imaging device and the fugitive gas.

7. A non-transitory, computer-readable medium containing instructions that, when executed by a hardware-based processor, causes the processor to perform stages for determining an emission rate of a fugitive gas, the stages comprising:

receiving, from an imaging device, a first gas density image;

detecting the fugitive gas in the first gas density image;

delineating the fugitive gas in the first gas density image;

sending, to the imaging device, instructions for adjusting a bearing of the imaging device bearing to optimize a position of the fugitive gas in the imaging device's field of view;

sending, to the imaging device, instructions for capturing a second gas density image;

receiving the second gas density image from the imaging device;

calculating, based on the second gas density image, the emission rate of the fugitive gas;

prior to sending the instructions for capturing the second gas density image, determining that a percentage of the fugitive gas is present in less than a threshold percentage of a vertical field of view of the first gas density image; and

in response to the percentage of the fugitive gas being below the threshold percentage, sending, to the imaging device, instructions for modifying a zoom of the imaging device so that the fugitive gas is present in more than the threshold percentage of a vertical field of view of the imaging device.

8. The non-transitory, computer-readable medium of claim 7 , the stages further comprising:

prior to sending the instruction for capturing the second gas density image, determining that the second gas density image includes a first portion of the fugitive gas that does not include all of the fugitive gas;

sending instructions to the imaging device for capturing a third gas density image that includes a second portion of the fugitive gas that is horizontally adjacent the first portion of the fugitive gas;

receiving, from the imaging device, the third gas density image; and

stitching the first gas density image and the second gas density image together.

9. The non-transitory, computer-readable medium of claim 7 , wherein calculating the emission rate of the gas density image comprises:

determining a volume of the fugitive gas based on the second density image and a measured distance between the imaging device and the fugitive gas;

receiving wind measurement data, the wind measurement data including measured wind speed and direction of wind in the vicinity of the fugitive gas;

inputting the wind measurement data and the volume of the fugitive gas into a closed surface integral equation for calculating flux; and

receiving an output from the closed surface integral equation, the output being a calculated amount of the fugitive gas.

10. The non-transitory, computer-readable medium of claim 7 , wherein calculating the emission rate of the gas density image comprises:

creating a cross-section in the fugitive gas in the second gas density image;

dividing the cross-section into a plurality of subsections based on spacing between image pixels in the second gas density image, each of the plurality of subsections having an estimated length;

determining a concentration of pixels for each of the plurality of subsections;

determining a fugitive gas plume travel velocity orthogonal to the cross-section;

inputting the determined concentration of pixels for each of the plurality of subsections, the fugitive gas plume travel velocity, and a weighting function into an equation; and

receiving, as an output from the equation, a flux rate of the fugitive gas for the cross-section.

11. The non-transitory, computer-readable medium of claim 10 , wherein calculating the emission rate of the gas density image further comprises:

calculating the flux rate for additional cross-sections of the fugitive gas;

determining that total flux of the fugitive gas by calculating an average of the calculated flux rate for the cross-section and the additional cross-sections of the fugitive gas.

12. The non-transitory, computer-readable medium of claim 10 , wherein the estimated length of the plurality of subsections is determined by converting pixels to length using degrees of the imaging device's bearing and a measured length between the imaging device and the fugitive gas.

13. A system for determining an emission rate of a fugitive gas, comprising:

an imaging device;

a memory storage including a non-transitory, computer-readable medium comprising instructions; and

a hardware-based processor that executes the instructions to carry out stages comprising:

receiving, from the imaging device, a first gas density image;

detecting the fugitive gas in the first gas density image;

delineating the fugitive gas in the first gas density image;

sending, to the imaging device, instructions for adjusting a bearing of the imaging device bearing to optimize a position of the fugitive gas in the imaging device's field of view;

sending, to the imaging device, instructions for capturing a second gas density image;

receiving the second gas density image from the imaging device;

calculating, based on the second gas density image, the emission rate of the fugitive gas;

prior to sending the instructions for capturing the second gas density image, determining that a percentage of the fugitive gas is present in less than a threshold percentage of a vertical field of view of the first gas density image; and

in response to the percentage of the fugitive gas being below the threshold percentage, sending, to the imaging device, instructions for modifying a zoom of the imaging device so that the fugitive gas is present in more than the threshold percentage of a vertical field of view of the imaging device.

14. The system of claim 13 , the stages further comprising:

prior to sending the instruction for capturing the second gas density image, determining that the second gas density image includes a first portion of the fugitive gas that does not include all of the fugitive gas;

sending instructions to the imaging device for capturing a third gas density image that includes a second portion of the fugitive gas that is horizontally adjacent the first portion of the fugitive gas;

receiving, from the imaging device, the third gas density image; and

stitching the first gas density image and the second gas density image together.

15. The system of claim 13 , wherein calculating the emission rate of the gas density image comprises:

determining a volume of the fugitive gas based on the second density image and a measured distance between the imaging device and the fugitive gas;

receiving wind measurement data, the wind measurement data including measured wind speed and direction of wind in the vicinity of the fugitive gas;

inputting the wind measurement data and the volume of the fugitive gas into a closed surface integral equation for calculating flux; and

receiving an output from the closed surface integral equation, the output being a calculated amount of the fugitive gas.

16. The system of claim 13 , wherein calculating the emission rate of the gas density image comprises:

creating a cross-section in the fugitive gas in the second gas density image;

dividing the cross-section into a plurality of subsections based on spacing between image pixels in the second gas density image, each of the plurality of subsections having an estimated length;

determining a concentration of pixels for each of the plurality of subsections;

determining a fugitive gas plume travel velocity orthogonal to the cross-section;

inputting the determined concentration of pixels for each of the plurality of subsections, the fugitive gas plume travel velocity, and a weighting function into an equation; and

receiving, as an output from the equation, a flux rate of the fugitive gas for the cross-section.

17. The system of claim 16 , wherein calculating the emission rate of the gas density image further comprises:

calculating the flux rate for additional cross-sections of the fugitive gas; and

determining that total flux of the fugitive gas by calculating an average of the calculated flux rate for the cross-section and the additional cross-sections of the fugitive gas.

Assignments (1)
ASSIGNMENT OF ASSIGNOR'S INTEREST Recorded Jan 9, 2024
From: SPECK, ANDREW J.; DOSHI, MANASI; ZIELINSKI, LUKASZ
To: SCHLUMBERGER TECHNOLOGY CORPORATION
Reel/Frame 066066/0543 →
Continuity (2)
Provisional Application 63508669 · Jun 16, 2023
Related Publication 20240420311A1 · Dec 19, 2024
References Cited (105)
US 11525764B2 · Bryning · 2022 [cited by examiner]
US 11630022B2 · Hedberg · 2023 [cited by applicant]
US 20150103178A1 · Itoh · 2015 [cited by applicant]
US 20150213697A1 · Knox · 2015 [cited by examiner]
US 20150369730A1 · Schmidt · 2015 [cited by examiner]
US 20160146696A1 · Steele · 2016 [cited by applicant]
US 20170024642A1 · Xiong · 2017 [cited by applicant]
US 20170364818A1 · Wu · 2017 [cited by applicant]
US 20180039885A1 · Albrecht · 2018 [cited by applicant]
US 20180292286A1 · Dittberner · 2018 [cited by applicant]
US 20180365522A1 · Jiang · 2018 [cited by applicant]
US 20190196892A1 · Matei · 2019 [cited by applicant]
US 20190285504A1 · Muralidhar · 2019 [cited by applicant]
US 20190285600A1 · Klein · 2019 [cited by applicant]
US 20190302013A1 · Wang · 2019 [cited by examiner]
US 20200011789A1 · Sandsten · 2020 [cited by examiner]
US 20200176088A1 · Kanamarlapudi · 2020 [cited by applicant]
US 20200182779A1 · Kasten · 2020 [cited by applicant]
US 20200226426A1 · Jarquin Arroyo · 2020 [cited by applicant]
US 20200393434A1 · Li · 2020 [cited by applicant]
US 20200398959A1 · Pomerantz · 2020 [cited by applicant]
US 20210140934A1 · Smith · 2021 [cited by applicant]
US 20210156793A1 · Leen · 2021 [cited by applicant]
US 20210255157A1 · Rashid · 2021 [cited by applicant]
US 20220244127A1 · Aubrey · 2022 [cited by applicant]
US 20220327708A1 · Satish · 2022 [cited by examiner]
US 20230111401A1 · Machover · 2023 [cited by applicant]
US 20230175914A1 · Diven · 2023 [cited by applicant]
US 20230177726A1 · Spears · 2023 [cited by applicant]
US 20230194483A1 · Some · 2023 [cited by examiner]
US 20230222271A1 · Spackova · 2023 [cited by applicant]
US 20230324430A1 · Reed · 2023 [cited by applicant]
US 20230326201A1 · Rashid · 2023 [cited by applicant]
US 20240067919A1 · Kamino · 2024 [cited by examiner]
US 20240200991A1 · Salman · 2024 [cited by applicant]
CN 109195132A · 2019 [cited by applicant]
CN 109614526A · 2019 [cited by applicant]
CN 113887044A · 2022 [cited by applicant]
CN 114329245A · 2022 [cited by applicant]
CN 114416692A · 2022 [cited by applicant]
CN 115018327A · 2022 [cited by applicant]
CN 115187767A · 2022 [cited by applicant]
EP 0943085B1 · 2004 [cited by applicant]
JP H06294500A · 1994 [cited by applicant]
JP H07083786A · 1995 [cited by applicant]
JP H07140031A · 1995 [cited by applicant]
JP H07190879A · 1995 [cited by applicant]
JP H07198523 · 1995 [cited by applicant]
JP 2996349B2 · 1999 [cited by applicant]
KR 101473956B1 · 2014 [cited by applicant]
KR 102015003470A · 2015 [cited by applicant]
KR 1020160123467A · 2016 [cited by applicant]
KR 20180041828A · 2018 [cited by applicant]
KR 1020190058230A · 2019 [cited by applicant]
WO 2019213280A1 · 2019 [cited by applicant]
WO 2020018867A1 · 2020 [cited by applicant]
WO 2020206008A1 · 2020 [cited by applicant]
WO 2020263693A1 · 2020 [cited by applicant]
WO 2021067844A1 · 2021 [cited by applicant]
WO 2021156864A1 · 2021 [cited by applicant]
WO 2022023226A1 · 2022 [cited by applicant]
WO 2022051572A1 · 2022 [cited by applicant]
WO 2022056152A1 · 2022 [cited by applicant]
WO 2023108041A1 · 2023 [cited by applicant]
WO 2023133345A1 · 2023 [cited by applicant]
WO 2024030523A1 · 2024 [cited by applicant]
WO 2024030525A1 · 2024 [cited by applicant]
WO 2024054640A1 · 2024 [cited by applicant]
WO 2024064113A1 · 2024 [cited by applicant]
WO 2024155857A1 · 2024 [cited by applicant]
Titchener et al. “Single photon Lidar gas imagers for practical and widespread continuous methane monitoring.” Applied Energy 306 (2022): 118086 (11 pages). [cited by applicant]
Saunois, M., et al. (2020), The Global Methane Budget 2000-2017, Earth Syst. Sci. Data, 12, 1561-1623. [cited by applicant]
Pomerantz A.E. et al., 2022. Present global warming: a justifiable and stable metric for evaluating short-lived climate pollutants. Environmental Research Letters, 17(11), p. 114052. (6 pages). [cited by applicant]
Cusworth , D.H. et al., 2021. Intermittency of large methane emitters in the Permian Basin. Environmental Science Technology Letters, 8(7), pp. 567-573. [cited by applicant]
Chakrabarti et al., Sep. 2022. Rapid Detection of Super-Emitters Utilizing an IoT-Enabled Continuous Methane Emissions Monitoring System. In SPE Annual Technical Conference and Exhibition. OnePetro. (6 pages). [cited by applicant]
Chakrabarti et al., Oct. 2022. Rapid Detection of Methane Super-Emitters Through Advanced Interpretation. In ADIPEC. OnePetro (6 pages). [cited by applicant]
Potyrailo, R.A. et al., 2020. Extraordinary performance of semiconducting metal oxide gas sensors using dielectric excitation. Nature Electronics, 3(5), pp. 280-289. [cited by applicant]
Zimmerle, D. Metec Controlled Test Protocol: Continuous Monitoring Emission Detection and Quantification, Energy Institute, Colorado State University. https://energy.colostate.edu/wp-content/uploads/sites/28/2021/03/Con… [cited by applicant]
Take control of your emissions management programme, downloaded on Dec. 15, 2023 from link https://sensorup.com/methane-emissions-management/ (12 pages). [cited by applicant]
BlueSky resources makes sense out of sensors, downloaded on Dec. 15, 2023 from linke https://blueskyresources.com/solutions/ (3 pages). [cited by applicant]
International Search Report and Written Opinion issued in PCT Application PCT/US2023/010475, dated May 2, 2023 (11 pages). [cited by applicant]
S. R. Hanna, G. A. Briggs, R. P. Hosker, 1982, Handbook on Atmospheric Diffusion, DOE/TIC-11223 (7 pages). [cited by applicant]
International Search Report and Written Opinion issued in PCT Application PCT/US2023/032307, dated Dec. 26, 2023 (10 pages). [cited by applicant]
International Search Report and Written Opinion issued in PCT Application PCT/US2023/029361, dated Nov. 14, 2023 (11 pages). [cited by applicant]
Yu et al., “Methane leakage source location based on a near-infrared off-axis integrated cavity output spectroscopic ppbv-level sensor and an optimized inverse model”, Infrared Physics & Technology, 2022, vol. 121, pp. … [cited by applicant]
Weidmann et al., “Locating and Quantifying Methane Emissions by Inverse Analysis of Path-Integrated Concentration Data Using a Markov-Chain Monte Carlo Approach”, ACS Earth and Space Chemistry, 2022, vol. 6, Issue 9, pp… [cited by applicant]
Ye et al., “Leakage source location based on Gaussain plume diffusion model using a near-infrared sensor”, Infrared Physics & Technology, 2020, vol. 109, pp. 1-5 (Jun. 19, 2020). [cited by applicant]
International Search Report and Written Opinion issued in PCT Application PCT/US2023/033106, dated Jan. 23, 2024 (7 pages). [cited by applicant]
Search Report issued in Norwegian Patent Application No. 20230720 on Jan. 23, 2024; 10 pages. [cited by applicant]
Majumder et al., “Development and Evaluation of Ensemble Learning-based Environmental Methane Detection and Intensity Prediction Models”, Dec. 17, 2023, pp. 1-16 (Year: 2023). [cited by applicant]
Andrews et al., “Quantitative Mapping of Methane Emissions in Oil & Gas Facilities”, Oct. 16-18, 2023, Society of Petroleum engineers, SPE-214909-MS, pp. 1-10 (Year: 2023). [cited by applicant]
Imbiriba et al., “Augmented physics-based machine learning for navigation and tracking”, 2023, IEEE Transactions on Aerospace and electronic Systems, pp. 1-13 (Year: 2023). [cited by applicant]
Chakrabarti et al., “Rapid Detection of Methane Super-Emitters Through Advanced Interpretation”, Oct. 31-Nov. 3, 2022, Society of Petroleum Engineers, SPE-210970-MS, pp. 1-6 (Year: 2022). [cited by applicant]
International Search Report and Written Opinion issued in PCT Application PCT/US2023/083724, dated Apr. 30, 2024 (9 pages). [cited by applicant]
Combined Search and Exam Report issued in United Kingdom patent Application No. GB2409015.1 dated Oct. 8, 2024, 5 pages. [cited by applicant]
International Search Report and Written Opinion issued in the PCT Application No. PCT/US2024/033102 dated Sep. 24, 2024, 10 pages. [cited by applicant]
Notice of Allowance issued in U.S. Appl. No. 18/537,978 dated Aug. 19, 2024, 10 pages. [cited by applicant]
Yu, L. et al., “Methane leakage source location based on a near-infrared off-axis integrated cavity output spectroscopic ppbv-level sensor and an optimized inverse model”, Infrared Physics Technology, 2022, vol. 121, pp… [cited by applicant]
Ye, W. et al., “Leakage source location based on Gaussain plume diffusion model using a near-infrared sensor”, Infrared Physics Technology, 2020, vol. 109, pp. 1-5. [cited by applicant]
“AERMOD” downloaded from Wikipedia on Jul. 31, 2024 from [https://en.wikipedia.org/wiki/AERMOD], Last Edited Mar. 4, 2022, 2 pages. [cited by applicant]
International Search Report and Written Opinion issued in the PCT Application No. PCT/US2023/029361 dated Nov. 14, 2023, 11 pages. [cited by applicant]
International Search Report and Written Opinion issued in the PCT Application No. PCT/US2024/012087 dated May 23, 2024, 11 pages. [cited by applicant]
International Preliminary Report on Patentability issued in the PCT Application No. PCT/US2023/010475 dated Jul. 25, 2024, 8 pages. [cited by applicant]
International Search Report and Written Opinion issued in PCT Application PCT/US2023/029363, dated Nov. 21, 2023 (10 pages). [cited by applicant]
Extended Search Report issued in European Patent Application No. 24182656.9 dated Nov. 14, 2024, 6 pages. [cited by applicant]