Well-to-Seismic Tie and Wavelet Extraction
Establish a precise time-depth relationship and extract a representative source wavelet.
Part 1/3 — Advanced Theory & Mechanics
The foundational efficacy of seismic inversion rests upon the fidelity of the well-to-seismic tie and the subsequent extraction of a representative wavelet. This process serves as the bridge between the depth-domain measurements of petrophysical boreholes and the time-domain observations of the seismic volume. By integrating sonic ($\Delta t$) and bulk density ($\rho_b$) logs, geophysicists derive the acoustic impedance ($Z_p = V_p \cdot \rho$), which is then converted into a reflectivity series via the Zoeppritz equations or their linear approximations, such as the Aki-Richards or Shuey formulations. The primary objective is to resolve the vertical misalignments induced by velocity anisotropy, dispersion, and the inherent differences between the high-frequency $(\sim 10^4 \text{ Hz})$ sonic tool measurements and the bandwidth-limited $(\sim 10^1–10^2 \text{ Hz})$ seismic signals.
Achieving a high-correlation coefficient between the synthetic seismogram and the extracted composite trace is the prerequisite for removing the source signature, thereby transitioning from an interface-based seismic interface to a layer-based rock property model.
The Physics of Time-Depth Conversion and Checkshot Calibration
The initial phase of any seismic tie involves the conversion of well logs from the depth domain ($z$) to the two-way travel time domain ($t$) using a velocity model derived from the sonic log ($DT$). However, sonic logs are susceptible to "cycle skipping" and are influenced by the disturbed zone around the borehole, leading to cumulative drift when integrated over thousands of feet. To mitigate this, checkshot surveys or Vertical Seismic Profiles (VSP) are employed as the ultimate reference for absolute travel time. The checkshot provides discrete time-depth pairs that account for the low-frequency velocity component often missed by the high-frequency sonic tool. A "drift curve" is calculated—representing the difference between the integrated sonic time and the checkshot time—and applied via a knee-point or spline interpolation to the sonic log to produce a calibrated velocity profile.
$$t(z) = 2 \int_{0}^{z} \frac{1}{V_{sonic}(z')} dz' + \Delta t_{drift}$$
> Expert Note: When reconciling sonic logs with seismic data, one must account for seismic dispersion. According to the Futterman (1962) model, higher frequency signals travel faster than lower frequency signals in anelastic media. Failure to apply a frequency-dependent dispersion correction (converting the kHz sonic velocity to the Hz seismic velocity) will result in a systematic depth-to-time mismatch that manifests as a progressive phase rotation in the correlation.
```mermaid
graph TD
A[Sonic & Density Logs] --> B[Acoustic Impedance Zp]
B --> C[Reflectivity Series R]
D[Checkshot/VSP Data] --> E[Drift Correction]
E --> F[Corrected Time-Depth Relationship]