IP Library Granted Patent US 12,530,595
Granted Patent B2
US 12,530,595 · App. 17/165,714 · Granted Jan 20, 2026

Identification of a characteristic of a physical system based on collaborative sensor networks

Inventors: Ryan Alan Goldhahn (Livermore, CA); Priyadip Ray (Livermore, CA); Braden C. Soper (Livermore, CA); Hao Chen (Boise, ID); Deepak Rajan (Livermore, CA)
Assignee: LAWRENCE LIVERMORE NATIONAL SECURITY, LLC
G06N3/126G06N3/086G06N7/01
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,530,595
App. No.
17/165,714
Granted
Jan 20, 2026
Kind
B2
Abstract

A node having a sensor and a computing device is provided for identifying a location of a source of a physical process. The node collects via the sensor a measurement of the physical process. The node repeatedly recalculates parameters until termination criterion is satisfied. The node receives parameters of a probability distribution and a gradient from other nodes. The node generates parameters based on the parameters and the gradients received from the other nodes. The node samples from a distribution of source locations based on the generated parameters. The node calculates a gradient derived from the sampled source locations, the generated parameters, and a joint probability of the sampled source locations and the measurement. The node sends to a subset of other nodes the generated parameters and the calculated gradient. When the termination criterion is satisfied, the generated parameters represent the probability distribution of the source location.

Claims (43)

1 . A method performed by a first node for identifying a source location of a potentially hazardous physical event, the method comprising:

collecting, via a sensor in the first node, a measurement of a physical quantity and is associated with the potentially hazardous physical event, wherein the first node is one of a plurality of nodes that are spatially distributed across a geographic area, each node of the plurality of nodes including a sensor and a computing device for identifying the source location of the potentially hazardous physical event, the computing device in each of the plurality of nodes including a processing unit and a memory;

repeatedly re-calculating parameters representing a probability distribution of the source location of the potentially hazardous physical event, by the computing device in the first node, until a termination criterion is satisfied, to iteratively increase an accuracy of an approximation of the source location of the potentially hazardous physical event, including,

receiving, from each other node of the plurality of nodes, a probability distribution description and a gradient, each said probability distribution description representing a probability distribution of the source location as determined by a corresponding other node, wherein each received probability distribution description is based on a model of measurements of the potentially hazardous physical event over time and space, the model including either of: i) a Black Box Variational Inference (BBVI) algorithm and an Evidence Lower Bound (ELBO), or ii) a Distributed Markov chain Monte Carlo (“MCMC”) algorithm,

wherein, for a message routed through intermediate nodes, probability distribution descriptions and gradients generated by said intermediate nodes are integrated into the message to increase an amount of received data contributing to the approximation of the source location without increasing network communication overhead between the first node and each other node;

executing, by the first node, either of: i) a BBVI algorithm and an ELBO, or ii) a MCMC algorithm with an output of an additional probability distribution description as the approximation of the source location of the potentially hazardous physical event, based on the probability distribution descriptions received from the other nodes and the gradients received from the other nodes;

sampling the additional probability distribution description of the first node to produce a plurality of sampled probability distribution values of the probability distribution description generated by the first node;

calculating, by the first node, a gradient from the sampled probability distribution values of the probability distribution description generated by the first node, based on an ELBO relating to the probability distribution values;

sending, from the first node to a subset of the other nodes, the probability distribution description generated by the first node and the gradient calculated by the first node;

reducing propagation of unhelpful data to conserve network bandwidth in response to determining, based on a signal-to-noise ratio of the measurement or a magnitude of the calculated gradient, that the generated probability distribution description and calculated gradient are not informative;

filtering defective data in response to determining that the generated probability distribution description, the calculated gradient, or both are inconsistent with those received from other nodes; and

updating the parameters based on the probability distribution descriptions and gradients generated by said intermediate nodes,

wherein when the termination criterion is satisfied, the probability distribution generated by the first node is indicative of the approximation of the source location of the potentially hazardous physical event; and

outputting, by the first node, an indication of the approximation of the source location of the potentially hazardous physical event to an external entity.

2 . The method of claim 1 wherein the physical quantity is a gas that is released at the source location.

3 . The method of claim 1 wherein the physical quantity is radiation that is released at the source location.

4 . The method of claim 1 wherein the probability distribution is a normal distribution and the parameters are mean and variance of the normal distribution.

5 . The method of claim 1 wherein the termination criterion is based on the calculated gradient.

6 . The method of claim 1 wherein the termination criterion is based on a number of iterations of re-calculating the parameters.

7 . A node for identifying a source location of a potentially hazardous physical event based on a measurement relating to the potentially hazardous physical event and based on probability distribution descriptions received from other nodes, each probability distribution description indicating a probability distribution of the source location, the node comprising:

a sensor to collect a measurement of a physical quantity and is relating to the potentially hazardous physical event, wherein the sensor comprises at least one of a gas sensor or a radiation sensor and the potentially hazardous physical event comprises at least one of a gas release or a radiation release;

a communication component using which the node can communicate with other nodes; and

a computing device that includes:

one or more processors; and

one or more non-transitory computer-readable storage mediums storing computer-executable instructions, execution of which by the one or more processors causes the computing device to perform a process of iteratively increasing an accuracy of an approximation of the source location of the potentially hazardous physical event, wherein the process includes to:

receive via the communication component from each of a plurality of other nodes a probability distribution description and a gradient derived from samplings of the probability distributions by the other nodes, each node of the plurality of other nodes including a sensor and a computing device to identify the source location of the potentially hazardous physical event, the computing device in each of the plurality of other nodes including a processor and a memory, wherein each received probability distribution description is based on a model of measurements of the potentially hazardous physical event over time and space, the model including either of: i) a Black Box Variational Inference (BBVI) algorithm and an Evidence Lower Bound (ELBO), or ii) a Distributed Markov chain Monte Carlo (“MCMC”) algorithm,

wherein, for a message routed through intermediate nodes, probability distribution descriptions and gradients generated by said intermediate nodes are integrated into the message to increase an amount of received data contributing to the approximation of the source location without increasing network communication overhead between the node and each other node;

execute either of: i) a BBVI algorithm and an ELBO, or ii) a MCMC algorithm with an output of an additional another probability distribution description as the approximation of the source location of the potentially hazardous physical event, based on the received probability distribution descriptions and gradients, wherein the probability distribution description generated by the node is based on a model of measurements of the potentially hazardous physical event over time and space

sample the additional probability distribution description of the node to produce a plurality of sampled probability distribution values of the probability distribution description generated by the node;

calculate a gradient derived from the sampled probability distribution values of the probability distribution description generated by the node, based on an ELBO relating to the probability distribution values;

send via the communication component to a subset of the other nodes, the generated probability distribution description and the calculated gradient;

reduce propagation of unhelpful data to conserve network bandwidth in response to determining, based on a signal-to-noise ratio of the measurement or a magnitude of the calculated gradient, that the generated probability distribution description and calculated gradient are not informative;

filter defective data in response to determining that the generated probability distribution description, the calculated gradient, or both are inconsistent with those received from other nodes;

update the parameters based on the probability distribution descriptions and gradients generated by said intermediate nodes,

wherein when a termination criterion is satisfied, the generated probability distribution represents the approximation of the source location of the potentially hazardous physical event; and

output an indication of the approximation of the source location of the potentially hazardous physical event to an external entity.

8 . The node of claim 7 wherein the other nodes are a subset of nodes that are selected randomly.

9 . The node of claim 7 wherein the computer-executable instructions further include instructions to determine whether the generated probability distribution description and calculated gradient are informative and when not informative, suppress the sending of the generated probability distribution description and the calculated gradient.

10 . The node of claim 7 wherein the computer-executable instructions further include instructions to receive from a sensor an indication to forward the calculated gradient to another node and forward the generated probability distribution description and calculated gradient to the other node.

11 . The node of claim 7 wherein the computer-executable instructions further include instructions to, prior to sending the generated probability distribution description and calculated gradient to other nodes, quantize the generated probability distribution description, the calculated gradient, or both to reduce bandwidth needed for the sending.

12 . The node of claim 7 wherein the computer-executable instructions to determine whether the generated probability distribution description, the calculated gradient, or both are inconsistent with those received from other nodes, and when inconsistent, suppress the sending of the generated probability distribution description and the calculated gradient.

13 . The method of claim 1 , further comprising randomly selecting the subset of the other nodes, to which to send the generated parameters and the calculated gradient.

14 . The node of claim 7 wherein the node is at least one of an unmanned aerial vehicle (“UAV”), an unmanned ground vehicle (“UGV”), an unmanned underwater vehicle (“UUV”) or an unmanned space vehicle (“USV”).

Assignments (2)
CONFIRMATORY LICENSE Recorded Mar 19, 2021
From: LAWRENCE LIVERMORE NATIONAL SECURITY, LLC.
To: U.S. DEPARTMENT OF ENERGY
Reel/Frame 055657/0593 →
ASSIGNMENT OF ASSIGNOR'S INTEREST Recorded Feb 2, 2021
From: GOLDHAHN, RYAN ALAN; CHEN, HAO; RAY, PRIYADIP; RAJAN, DEEPAK; SOPER, BRADEN C.
To: LAWRENCE LIVERMORE NATIONAL SECURITY, LLC
Reel/Frame 055118/0953 →
Continuity (2)
Provisional Application 62969352 · Feb 3, 2020
Related Publication 20210241124A1 · Aug 5, 2021
References Cited (15)
US 20090125274A1 · Waldock · 2009 [cited by examiner]
US 20190148246A1 · Zhan · 2019 [cited by examiner]
US 20190259474A1 · Wang · 2019 [cited by examiner]
US 20200110175A1 · Chen · 2020 [cited by examiner]
US 20210158028A1 · Wu · 2021 [cited by examiner]
Ahn, S. et al., “Distributed stochastic gradient momc,” in International conference on machine learning, 2014, pp. 1044-1052. [cited by applicant]
Cattivelli, F. S., et al. (2009). Diffusion LMS strategies for distributed estimation. IEEE Transactions on Signal Processing, 58 (3), 1035-1048. [cited by applicant]
Draxler, R., R., “Forty-eight-hour atmospheric dispersion forecasts at selected locations in the united states,” 1981. [cited by applicant]
Erseghe, T. (2012). A distributed and scalable processing method based upon Admm. IEEE Signal Processing Letters, 19 (9), 563-566. [cited by applicant]
Ihler, A. T., et al. (2005). Nonparametric belief propagation for self-localization of sensor networks. IEEE Journal on Selected Areas in Communications, 23(4), 809-819. [cited by applicant]
Neiswanger, C. et al., “Asymptotically exact, embarrassingly parallel mcmc,” arXiv preprint arXiv:1311.4780, 2014. [cited by applicant]
Ranganath, R., et al. “Black Box Variational Inference,” Proceedings of the 17th International Conference on Artificial Intelligence and Statistics (AISTATS) 2014, Reykjavik, Iceland. Jmlr: W&Cp vol. 33, pp. 814-822. [cited by applicant]
Robbins, H., et al., “A stochastic approximation method,” The annals of mathematical statistics, pp. 400-407, 1951. [cited by applicant]
Scott, S. L. et al., “Bayes and big data: The consensus monte carlo algorithm,” International Journal of Management Science and Engineering Management, vol. 11, No. 2, 2016. [cited by applicant]
Welling, M., et al., “Bayesian learning via stochastic gradient langevin dynamics,” in Proceedings of the 28th International Conference on Machine Learning (ICML-11), 2011, pp. 681-688. [cited by applicant]