Analysis and Optimization of Seismic Monitoring Networks with Bayesian Optimal Experiment Design

TL;DR

This paper develops a Bayesian optimal experimental design framework to optimize seismic sensor networks for improved earthquake detection and localization, considering sensor placement, types, and fidelity.

Contribution

It introduces a comprehensive Bayesian OED framework tailored for seismic monitoring, enabling optimized sensor configurations based on information gain about seismic events.

Findings

Optimized sensor placement improves seismic source localization accuracy.

Trade-offs between sensor fidelity and model uncertainty are quantified.

Sensor network configurations significantly influence uncertainty reduction.

Abstract

Monitoring networks increasingly aim to assimilate data from a large number of diverse sensors covering many sensing modalities. Bayesian optimal experimental design (OED) seeks to identify data, sensor configurations, or experiments which can optimally reduce uncertainty and hence increase the performance of a monitoring network. Information theory guides OED by formulating the choice of experiment or sensor placement as an optimization problem that maximizes the expected information gain (EIG) about quantities of interest given prior knowledge and models of expected observation data. Therefore, within the context of seismo-acoustic monitoring, we can use Bayesian OED to configure sensor networks by choosing sensor locations, types, and fidelity in order to improve our ability to identify and locate seismic sources. In this work, we develop the framework necessary to use Bayesian OED to optimize a sensor network's ability to locate seismic events from arrival time data of detected seismic phases at the regional-scale. Bayesian OED requires four elements: 1) A likelihood function that describes the distribution of detection and travel time data from the sensor network, 2) A Bayesian solver that uses a prior and likelihood to identify the posterior distribution of seismic events given the data, 3) An algorithm to compute EIG about seismic events over a dataset of hypothetical prior events, 4) An optimizer that finds a sensor network which maximizes EIG. Once we have developed this framework, we explore many relevant questions to monitoring such as: how to trade off sensor fidelity and earth model uncertainty; how sensor types, number, and locations influence uncertainty; and how prior models and constraints influence sensor placement.