Quantum computing has advanced enough that solving basic inverse problems in exploration seismology may be beneficial to show proof-of-concept results. The problems solved should be simple enough to work with the current early state of quantum computing hardware, which would be comparable to problems solved in the early stages of classical computing. We introduce an iterative algorithm for solving PDE-constrained optimization problems on a quantum annealer. We then use this method to invert for subsurface P-wave velocity—using a hybrid of classical and quantum computing—in a simple seismic experiment. Presentation Date: Tuesday, October 13, 2020 Session Start Time: 9:20 AM Presentation Time: 10:10 AM Location: Poster Station 11 Presentation Type: Poster
Quantum computing has advanced enough that solving basic inverse problems in exploration seismology may be beneficial to show proof-of-concept results. The problems solved should be simple enough to work with the current early state of quantum computing hardware, which would be comparable to problems solved in the early stages of classical computing. We introduce an iterative algorithm for solving PDE-constrained optimization problems on a quantum annealer. We then use this method to invert for subsurface P-wave velocity—using a hybrid of classical and quantum computing—in a simple seismic experiment. Presentation Date: Tuesday, October 13, 2020 Session Start Time: 9:20 AM Presentation Time: 10:10 AM Location: Poster Station 11 Presentation Type: Poster
Abstract ABSTRACT Quantum computing harnesses the principles of quantum mechanics to solve problems that are intractable for classical computers. Quantum annealing, a specialized approach within quantum computing, is particularly effective for optimization tasks, as it leverages quantum tunneling to escape local minima and efficiently explore complex energy landscapes. In geosciences, many problems are framed as high-dimensional optimization problems, including seismic inversion, which aims to estimate subsurface impedances from seismic data for accurate geological interpretation and resource exploration. A novel application of quantum computing for seismic inversion is presented, marking the first instance of inverting seismic data to estimate both P-wave and S-wave impedances using a quantum annealer. Building upon the prior work, which demonstrated the estimation of acoustic impedances from post-stack data using a two-step framework, we propose an enhanced workflow capable of inverting both post-stack and pre-stack seismic data in a single step. This advancement significantly reduces the number of qubits per model parameter (from 20 to 5) while improving computational speed (from 20 seconds to 6.3 seconds). The seismic inversion is implemented using the D-Wave Leap hybrid solver, achieving impedance estimation within 4-9 seconds, with the quantum processing unit (QPU) contributing just 0.043-0.085 seconds. Comparative analysis with simulated annealing reveals that quantum annealing produces impedance models closely matching true values in a single epoch, whereas simulated annealing requires 10 epochs for improved accuracy. These findings underscore the transformative potential of quantum computing for real-time, high-precision seismic inversion, marking a crucial step toward fully quantum-driven geophysical solutions.
Quantum computing exploits quantum mechanics to perform certain computations more efficiently than classical computers. Current quantum computers have performed carefully tailored computational tasks that would be difficult or impossible for even the fastest supercomputers in the world. This “quantum supremacy” result demonstrates that quantum computing is more powerful than classical computing in some computational regimes. At present, it is unknown if any computational problems related to the Earth's subsurface fall within these regimes. Here, we describe an approach to performing seismic inverse analysis that combines a type of quantum computer called a quantum annealer with classical computing. This approach improves upon past work on applying quantum computing to the subsurface (via subsurface hydrology) in two ways. First, the seismic inverse problem enables better performance from the quantum annealer because of the Earth's relatively narrow distribution of P-wave velocities compared to the broad distribution of hydraulic conductivities. Second, we develop an iterative approach to quantum-computational inverse analysis, which works with a realistic set of observations. By contrast, the previous method used an inverse method that depended on an impractically dense set of observations. In combination, these two advances significantly narrow the gap a quantum-computational advantage for a practical subsurface geoscience problem. Closing the gap completely requires more work, but has the potential to dramatically accelerate inverse analyses for subsurface geoscience.
Quantum computing exploits quantum mechanics to perform certain computations more efficiently than classical computers. Current quantum computers have performed carefully tailored computational tasks that would be difficult or impossible for even the fastest supercomputers in the world. This “quantum supremacy” result demonstrates that quantum computing is more powerful than classical computing in some computational regimes. At present, it is unknown if any computational problems related to the Earth's subsurface fall within these regimes. Here, we describe an approach to performing seismic inverse analysis that combines a type of quantum computer called a quantum annealer with classical computing. This approach improves upon past work on applying quantum computing to the subsurface (via subsurface hydrology) in two ways. First, the seismic inverse problem enables better performance from the quantum annealer because of the Earth's relatively narrow distribution of P-wave velocities compared to the broad distribution of hydraulic conductivities. Second, we develop an iterative approach to quantum-computational inverse analysis, which works with a realistic set of observations. By contrast, the previous method used an inverse method that depended on an impractically dense set of observations. In combination, these two advances significantly narrow the gap a quantum-computational advantage for a practical subsurface geoscience problem. Closing the gap completely requires more work, but has the potential to dramatically accelerate inverse analyses for subsurface geoscience.
This study demonstrates the application of quantum computing based quantum annealing to seismic traveltime inversion, a critical approach for inverting highly accurate velocity models. The seismic inversion problem is first converted into a Quadratic Unconstrained Binary Optimization problem, which the quantum annealer is specifically designed to solve. We then solve the problem via quantum annealing method. The inversion is applied on a synthetic velocity model, presenting a carbon storage scenario at depths of 1000-1300 m. As an application example, we also show the capacity of quantum computing to handle complex, noisy data environments. This work highlights the emerging potential of quantum computing in geophysical applications, providing a foundation for future developments in high-precision seismic imaging.
Accurate modeling of seismic wave propagation in the subsurface of the earth is essential for understanding earthquake dynamics, characterizing seismic hazards on global scales and hydrocarbon reservoir exploration and monitoring on local scales. These are among the most challenging computational problems in geoscience. Despite algorithmic advances and the increasingly powerful computational resources currently available, including fast CPUs, GPUs and large volumes of computer memory, there are still daunting computational challenges in simulating 3D seismic wave propagation in complex earth environments. Recent advances in quantum computing are suggestive that geoscience may soon begin to benefit from this promising field. For example, Finite Difference (FD) modeling is the most widely used method to simulate seismic wave propagation. In the frequency domain, FD methods reduce solutions of the wave equation into systems of linear equations; such systems are just the type that quantum algorithms may be capable of solving with exponential speedup, in comparison with classical algorithms. For the computational geophysicist, to prepare to take advantage of these speed-ups, which could arrive in as few as 5-10 years, the tasks at hand are (1) to become familiar with the logic and concepts associated with quantum computing, and (2) to map our key computational algorithms (e.g., frequency domain FD) to this domain. Presentation Date: Monday, October 15, 2018 Start Time: 1:50:00 PM Location: 204B (Anaheim Convention Center) Presentation Type: Oral
Quantum computing is increasingly recognized as a transformative technology in geophysics, offering the potential to significantly enhance computational power and efficiency. This advancement holds promise for the simulation and processing of complex geological data at unprecedented speeds. Recent studies have begun exploring the application of quantum computing methods to simplified versions of seismic inversion problems, underscoring the technology’s capability to address real-world inverse problems. The primary aim of this study is to tackle a realistic, scalable, and business-relevant problem by estimating acoustic impedances from seismic trace data using a quantum computer. To the best of our knowledge, this is the first instance of seismic impedances being predicted from seismic data through quantum computing, with a discussion on the advantages of solving the inverse problem on quantum processing units (QPUs). In this paper, we leverage the D-Wave quantum annealer to tackle the poststack seismic inversion problem, employing a novel two-step workflow. In the first step, we utilize the quantum annealer to estimate reflectivities from seismic data. Subsequently, these estimated normal-incidence reflectivities serve as the basis for predicting acoustic impedances with the same quantum technology. To validate the effectiveness of our approach, we present five examples, juxtaposing the impedance predictions from the D-Wave quantum annealer against those obtained through Simulated Annealing—a stochastic global optimizer traditionally employed for seismic inversion. Remarkably, the impedances derived from the quantum annealer matched closely with the true values in just a single epoch, whereas simulated annealing required 10 epochs to achieve a similar accuracy. Furthermore, the QPUs within our hybrid solver took only about 0.08 seconds to estimate these seismic impedances. This is remarkably efficient when compared to the time needed by both the classical component of the hybrid solver and simulated annealing, which each took over 10 seconds. This highlights the potential of QPUs to entirely resolve the seismic inverse problem in under a second, underscoring the transformative impact that quantum computing can have in the field of geophysics.