Improving explicit seismic depth migration with a stabilizing Wiener filter and spatial resampling
Bibliographic record
Abstract
Abstract We present a new approach to the design and implementation of explicit wavefield extrapolation for seismic depth migration in the space-frequency domain. Instability of the wavefield extrapolation operator is addressed by splitting the operator into two parts, one to control phase accuracy and a second to improve stability. The first partial operator is simply a windowed version of the exact operator for a half step. The second partial operator is designed, using the Wiener filter method, as a band-limited, least-squares inverse of the first. The final wavefield extrapolation operator for a full step is formed as a convolution of the first partial operator with the complex conjugate of the second. This resulting wavefield extrapolation operator can be designed to have any desired length and is generally more stable and more accurate than a simple windowed operator of similar length. Additional stability is gained by reducing the amount of evanescent filtering and by spatially downsampling the lower temporal frequencies. The amount of evanescent filtering is controlled by building two operator tables, one corresponding to significant evanescent filtering and the other to very little evanescent filtering. During the wavefield extrapolation process, most steps are taken with the second table while the first is invoked only for roughly every tenth step. Also, the data are divided into frequency partitions that are optimally resampled in the spatial coordinates to further enhance the performance of the extrapolation operator. Lower frequencies are downsampled to a larger spatial sample size. Testing of the algorithm shows accurate, high-angle impulse responses and run times comparable to the phase shift method of time migration. Images from trial depth migrations of the Marmousi model show very high resolution.
Fetched live from OpenAlex and de-inverted. Abstracts are not stored in this database: the inverted indexes are 8.6 GB of the frame’s 9.3 GB of text, and the host has 13 GB free.
How this classification was reachedexpand
Full frame machine prediction
Teacher imitationNot calibrated prevalence, not ground truth. Human validation pending. The Gemma side is a direct model label for every work in the frame, read from the title-only record. The Codex side is a classifier learned from the 10,348 direct Codex labels and calibrated to design-weighted sample rates; fields without enough sample support carry no Codex call. Candidate is the union of the two sides; consensus is their intersection. These outputs are machine_predicted_unvalidated and are not human labels.
Distilled classifier scores by category (both heads)
| Category | Codex | Gemma |
|---|---|---|
| Metaresearch | 0.001 | 0.002 |
| Meta-epidemiology (narrow) | 0.001 | 0.000 |
| Meta-epidemiology (broad) | 0.000 | 0.000 |
| Bibliometrics | 0.000 | 0.000 |
| Science and technology studies | 0.000 | 0.000 |
| Scholarly communication | 0.000 | 0.001 |
| Open science | 0.001 | 0.001 |
| Research integrity | 0.001 | 0.001 |
| Insufficient payload (model declined to judge) | 0.001 | 0.000 |
Machine scores (provisional)
The two teacher heads of the student model, read on this work. A score orders the frame for review; it never asserts a category, and the validation status ships verbatim with every row.
Baseline scores from an immature model (maturity gate not passed, 7 training rounds). Scores rank; they never assert a category.
score_only:v0-immature-baseline · verbatim from the scoring run: score_only means the number may rank works, and no category label ships from itClassification
machine, unvalidatedMachine predicted; a candidate call from one source (direct Gemma or distilled Codex), not a consensus.
How this classification was reached, model by model and score by score, is at the end of the page under "How this classification was reached".