Abstract:Inverting seismic data to build 3D geological structures is a challenging task due to the overwhelming amount of acquired seismic data, and the very-high computational load due to iterative numerical solutions of the wave equation, as required by industry-standard tools such as Full Waveform Inversion (FWI). For example, in an area with surface dimensions of 4.5km $\times$ 4.5km, hundreds of seismic shot-gather cubes are required for 3D model reconstruction, leading to Terabytes of recorded data. This paper presents a deep learning solution for the reconstruction of realistic 3D models in the presence of field noise recorded in seismic surveys. We implement and analyze a convolutional encoder-decoder architecture that efficiently processes the entire collection of hundreds of seismic shot-gather cubes. The proposed solution demonstrates that realistic 3D models can be reconstructed with a structural similarity index measure (SSIM) of 0.8554 (out of 1.0) in the presence of field noise at 10dB signal-to-noise ratio.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is to efficiently reconstruct 3D geological structure models from vast amounts of seismic data through deep - learning methods. Specifically, traditional iterative numerical solutions based on the wave equation (such as full - waveform inversion, Full Waveform Inversion, FWI) face the challenges of huge computational volume and large data volume when processing large - scale seismic data. For example, in a 4.5km×4.5km area, hundreds of seismic ray - beam cubes are required to reconstruct a 3D model, which will result in several terabytes of data.
To solve these problems, this paper proposes a deep - learning solution based on the convolutional encoder - decoder architecture to achieve the following goals:
1. **Efficiently reconstruct complex 3D models**: Through efficient input - data dimension reduction, this model can complete inference within an average of 0.165 seconds (using a single NVIDIA A100 GPU), which is much faster than any iterative global - optimization solver.
2. **Be robust to noise**: This method shows inherent robustness to noise in the recorded seismic data and can maintain high reconstruction quality even in the case of a 10 - dB signal - to - noise ratio.
3. **Evaluate performance in real - world scenarios**: This method has been evaluated on a real 3D geological model, verifying its effectiveness in practical applications.
### Mathematical Formulas
The main formulas involved in the paper include:
- Wave equation:
\[
\frac{\partial^2 u}{\partial t^2}-V_P^2\nabla(\nabla\cdot u)-V_S^2\Delta u = f
\]
where \(u\) is the seismic - wave displacement, \(V_P\) is the P - wave velocity (compression/rarefaction), \(V_S\) is the S - wave velocity (shear stress), and \(f\) is the source function.
- Acoustic - wave equation (assuming no density change in the medium):
\[
\frac{\partial^2 u}{\partial t^2}-V^2\nabla^2 u = f
\]
where \(u\) is the wave displacement, \(V\) is the P - wave velocity model, and \(f\) is the perturbation - source (i.e., seismic - source) function.
- Expression for the inversion problem:
\[
d = F(m)+\epsilon
\]
where \(d\) is the synthetic seismic data, \(F\) is the forward operator acting on the artificial model \(m\), and \(\epsilon\) is the environmental noise.
- Data cube after dimension reduction:
\[
\bar{d}(R_x,R_y,t)=b(t)\frac{1}{N_x\times N_y}\sum_{S_x = 1}^{N_x}\sum_{S_y = 1}^{N_y}d(S_x,S_y,R_x,R_y,t)
\]
where \(b(t)\) is a monotonically increasing time - enhancement function used to compensate for the attenuation of reflected waves from the lowest geological layer.
### Summary
This paper aims to solve the computational bottleneck and noise - sensitivity problems faced by traditional seismic - inversion methods when processing large - scale data through deep - learning techniques, especially the convolutional encoder - decoder architecture, thereby achieving efficient and high - quality 3D geological - structure reconstruction.