GPU accelerated sparse curvelet-constrained wavefield reconstruction inversion with source estimation: Application to Chevron Benchmark 2014 blind test data set

Full-waveform inversion (FWI) is a pivotal tool in seismic exploration and seismology but encounters a persistent challenge known as “cycle-skipping,” causing it to frequently converge to local minima. Wavefield reconstruction inversion (WRI) has emerged as a potential solution to mitigate the risk of cycle skipping. Researchers have provided numerical examples indicating that WRI is less prone to being ensnared in local minima caused by the absence of low-frequency component data and the absence of a well-defined initial model compared with conventional FWI. Despite its potential, the computational demands of WRI have hindered its widespread application, especially in scenarios involving ocean towed streamer seismic acquisition and unknown sources, where the augmented systems differ from source to source. In our study, we introduce a novel approach — sparse curvelet-constrained WRI with source estimation (WRI-SE-CC) — accelerated by graphics processing unit (GPU). Real-time source function estimation is achieved through the variable projection method, and noise-related artifacts are suppressed using sparse curvelet constraints. By optimizing the utilization of hundreds of computation processors within a GPU for parallel computing of matrix-vector multiplications, we present a GPU-based grouped conjugate gradient method to accelerate the computation of WRI-SE-CC. Numerical experiments demonstrate a significant 240-fold acceleration compared with the preconditioned conjugate gradient using one CPU core for computations involving multiple sources. Inversion experiments with the overthrust model demonstrate the capability of our method in mitigating local minima and suppressing noise-related artifacts. Furthermore, we validate the framework on the Chevron 2014 blind test data set, showcasing its effectiveness in addressing practical challenges in the field.