Variational Quantum Physics-Informed Neural Networks for Hydrological PDE-Constrained Learning with Inherent Uncertainty Quantification

Solving reservoir seepage partial differential equations (PDEs) is critical for hydrocarbon development optimization and production forecasting. Conventional numerical solvers suffer from mesh errors and high computational cost. Standard Physics-Informed Neural Networks (PINNs) are limited by inefficient parameters. They also face weak high-dimensional representation and poor nonlinear fitting ability. We adapt the Discrete Variable (DV)-Circuit Quantum-Classical Physics-Informed Neural Network (QCPINN) to tackle these drawbacks. This work applies the method to four representative reservoir seepage models. The models cover heterogeneous single-phase pressure diffusion. They also include nonlinear Buckley-Leverett two-phase waterflooding. Another case is adsorption-coupled convection-diffusion compositional flow. The last one refers to fully coupled pressure-saturation two-phase flow with exponential permeability heterogeneity. The QCPINN combines classical preprocessing and postprocessing networks with a DV quantum core. It leverages quantum entanglement and superposition to boost high-dimensional feature mapping. Physical constraints are embedded within the framework to guarantee physically consistent solutions. We evaluate three quantum circuit topologies including Cascade, Cross-mesh and Alternate. Numerical tests prove QCPINNs outperform classical PINNs in prediction accuracy, and the circuit topologies show case-specific performance with clear applicable ranges. This work verifies QCPINN's practicability for reservoir engineering. It bridges quantum computing research with oilfield industrial applications.

Physics-informed neural networks (PINNs) have emerged as a versatile approach for solving nonlinear partial differential equations (PDEs), yet achieving high accuracy efficiently using these techniques remains challenging for high-dimensional or multiscale systems. Here, we present a hybrid quantum-classical framework that enhances Quantum PINNs (QPINNs) through adaptive collocation point sampling and loss-aware attention mechanisms. By dynamically prioritizing points in regions with large PDE residuals or steep solution gradients, our method mitigates the spectral bias inherent in conventional PINNs. Current Quantum Physics-Informed Neural Networks are commonly assumed to be limited by the expressive power of quantum circuits. In our work, we observed that, across diverse differential equations, optimization - not only expressivity - can be an important bottleneck. Furthermore, a trainable loss-weighting scheme balances contributions from physics residuals, boundary conditions, and data fidelity during training. Integrating these strategies with quantum computing techniques (including variational quantum circuits and quantum gradient estimation) can yield at least a 60% improvement in solution accuracy under specific regimes for benchmark fluid flows and reaction-diffusion systems. Finally, we argue that merely increasing model expressivity is insufficient for resolving complex PDEs via QPINNs, as they remain constrained by the structural optimization limitations of classical PINNs. This framework provides a scalable pathway for quantum-enhanced scientific machine learning, bridging physics-based modeling with emerging quantum computational capabilities.

Full waveform inversion (FWI) reconstructs heterogeneous material properties from receiver data but remains computationally demanding. Physics-informed neural networks (PINNs) and their domain-decomposed variants (FBPINNs) offer a mesh-free alternative but face convergence challenges when representing complex velocity fields. We present a hybrid quantum-classical FBPINN for acoustic FWI, bringing together quantum computing and classical machine learning, in which the decomposed wavefield network and the global velocity network are implemented as classical-to-quantum pipelines terminating in parameterized quantum circuits (PQCs). The PQCs are realized as differentiable JAX statevector simulators, enabling end-to-end automatic differentiation through the classical PINN, the quantum circuit, and the physics-informed loss. On a geophysical anomaly benchmark, the quantum hybrid reaches a lower L1 velocity error than the primary classical FBPINN baseline in approximately 8x fewer training iterations, despite using approximately 33% fewer trainable parameters, and it outperforms all 15 classical hyperparameter variants tested. A second benchmark (checkerboard) demonstrates the generality of the inversion pipeline, confirming that the quantum hybrid architecture can recover structured spatial variations beyond the localized anomaly benchmark. Our framework is broadly applicable to wave-based inverse problems beyond geophysics, including medical ultrasound tomography and non-destructive evaluation.