3-D elastic prestack, reverse‐time depth migration
By combining and extending previous algorithms for 2-D prestack elastic migration and 3-D prestack acoustic migration, a full 3-D elastic prestack depth migration algorithm is developed. Reverse‐time extrapolation of the recorded data is by 3-D elastic finite differences; computation of the image time for each point in the 3-D volume is by 3-D acoustic finite differences. The algorithm operates on three‐component, vector‐wavefield common‐source data and produces three‐component vector reflectivity distributions. Converted P‐to‐S reflections are automatically imaged with the primary P‐wave reflections. There are no dip restrictions as the full wave equation is used. The algorithm is illustrated by application to synthetic data from three models; a flat reflector, a dipping truncated wedge overlying a flat reflector, and the classical French double dome and fault model.