Organizations: Department of Electrical and Computer Engineering, University of New Mexico, Albuquerque, NM 87131 USA · Department of Mechanical Engineering, University of California, Merced, CA 95343 USA · Department of Electrical Engineering, Colorado School of Mines, Golden, CO 80401 USA · Department of Electrical Engineering, New Mexico Tech, Socorro, NM 87801 USA · Autonomous Solutions, Inc., Mendon, UT 84325 USA · Department of Mechanical Engineering, New Mexico Tech, Socorro, NM 87801 USA
Undocumented orphaned wells pose significant health and environmental risks to nearby communities by releasing toxic gases and contaminating water sources, with methane emissions being a primary concern. Traditional survey methods such as magnetometry often fail to detect older wells effectively. In contrast, aerial in-situ sensing using unmanned aerial vehicles (UAVs) offers a promising alternative for methane emission detection and source localization. This study presents a robust and efficient framework based on a multi-agent deep reinforcement learning (MARL) algorithm for the chemical plume source localization (CPSL) problem. The proposed approach leverages virtual anchor nodes to coordinate UAV navigation, enabling collaborative sensing of gas concentrations and wind velocities through onboard and shared measurements. Source identification is achieved by analyzing the historical trajectory of anchor node placements within the plume. Comparative evaluations against the fluxotaxis method demonstrate that the MARL framework achieves superior performance in both localization accuracy and operational efficiency.
Figures & tables
Figure 1: In-situ gas concentration measurement results and two UAV-sensor products [ 3 , 4 ] .
Figure 2: Envisioned Multi-UAV CPSL scenario.
Figure 3: Topology of the considered CPSL Problem.
Figure 4: Search area of the defined CPSL problem.
Figure 5: Concentration plot when a UAV is at h=2 m (left) and h=5 m (right) at a point 6 m downwind from the emitter location along the plume’s center-line. The time resolution is 0.05 sec. and the filament release rate is Nf=50 filaments/sec.
Notation
Value
Emitter Source Location (pes)
(80, 60) m
Mean Direction of Wind (Vw)
[1,0] m/s
Colored Noise Constants ([a,b,G])
[0.005,0.02,1]
Concentration Sensor Filter Coeffs. (B,ch)
(0.1,10−4)
Diffusivity Constants (Kx,Ky)
(1000,1000) m2 /s
Relative Diffusion Spectral Density (σvm,1,σvm,2)
(2,2) m/s/ Hz
Table I: Plume Parameter Values
Figure 6: Plume dynamics due to relative diffusion and constant advection, i.e. the mean flow field has velocity Vw=[1,0] m/s, and the emitter source location is at (x,y)=(80,60) m.
Figure 7: (a) Relative angle information of UAV i′ observed by UAV i (b) relative formation information observed by UAV i .
Figure 8: (a) Example of conditions where the anchor information is updated and (b) Example of conditions where the anchor information is NOT updated.
Description
Value
Geographical Area ( X×Y )
200m×200m
Grid element size ( Δx×Δy )
2m×2m
Number of UAVs ( N )
3
Number of Obstacles ( M )
5
UAV Object Sensing Radius ( rs )
20 m
Collision Safety Distance ( ϵ )
1 m
Table II: Simulation Parameter Values
Notation
Value
Emitter Location
(80, 60) m
(60, 120) m
Noise Constants ([a,b,G])
[0.005,0.02,1]
[0.005,0.02,3]
[0.005,0.02,5]
Table III: Test Plume Parameter Values
Figure 9: Simulation results for Plume Example 1: (a) Mean Total Reward Rui received per agent, (b) Mean Target Found Rate, (c) Mean Final Centroid Error, (d) Mean Agent-to-Agent and Agent-to-Obstacle Collision Rate.
Figure 10: Reward difference during training.
Figure 11: Plume distribution with different wind conditions
Scenario Setup
Succ. Rate
Mean Location Error (m)
No Offset
With Offset
Wind Meander
Emitter Location
All
Succ.
All
Succ.
No
[80, 60]
1.00
2.34
2.34
1.55
1.55
[60, 120]
1.00
2.17
2.17
1.52
1.52
Small
[80, 60]
1.00
2.16
2.16
1.76
1.76
[60, 120]
0.96
2.61
2.40
1.7
1.48
Table IV: Success rate and location error (i.e, final distance to emitter) under various environment conditions
Figure 12: Example of different UAVs behaviors during a successful episode
Figure 13: Trajectory and location error evolution for MARL and fluxotaxis under varying wind meandering conditions
Figure 14: CDF of final distance to emitter for MARL and fluxotaxis (100 tests per case)
Electrical Engineering & Computer Science Graduate Program, University of California, Merced, 5200 N. Lake Rd, Merced, CA 95343, USA · Dept. of Mechanical Engineering, University of California, Merced, 5200 N. Lake Rd, Merced, CA 95343, USA