The specialty of seismic imaging in foothill area originates from the complex surface conditions,which leads to low S/N and
serious moveout between gathers.Because of the rugged topography and near-surface lateral velocity variation,the conventional
root-mean-square velocity analysis based on CMP gathers stacking is no longer effective.Up to now,there does not exist an
effective imaging processing procedure widely accepted for foothill seismic data.After analyzing the seismic wave propagation
and the limitation of conventional processing methods,we proposed a new imaging workflow.The workflow includes three key
steps:①choosing a smoothed datum as close to the topography to remove the high wavenumber moveout.The prestack depth
migration and velocity model building are also based on the datum;②Kirchhoff PSDM based Monta Carlo velocity inversion is
used to estimate the background and smoothed initial velocity.③beam ray PSDM based velocity tomography in angle domain
is used to estimate the more fine and accurate migration velocity.In addition,the in-phase stacking the reflections from half
Fresnel band is used to do the noise attenuation,which is very important for the effective migration velocity estimation and
velocity model building.The near-surface velocity model building replaces the so-called static correction.The basic numerical
tests demonstrate the effectiveness of this procedure.