ABSTRACTSeismic images of Earth's subsurface are essential for oil and gas exploration. Gaussian beam migration is popular for seismic imaging owing to its flexibility and efficiency of implementation. However, the practical use of classic Gaussian beam migration is limited for complicated structures such as subsalt because of poor seismic illumination in those areas. We propose an illumination‐compensated Gaussian beam migration under the framework of least‐squares migration, enhancing the subsalt imaging effectively. A novel scheme based on the Born modelling and migration of Gaussian beam is developed to estimate the point spreading function, with which the illumination compensation can be efficiently implemented in the local wavenumber domain. As the aim of least‐squares migration, the proposed illumination‐compensated Gaussian beam migration produces true‐amplitude images in subsalt areas, which facilitates the seismic interpretation of subsalt structures. The total computational cost of the proposed method includes one Born modelling process and two conventional Gaussian beam migrations, and thus it is much more efficient than the classic least‐squares migration which requires multiple iterations. Numerical examples have demonstrated the effectiveness of the adaptive strategies, manifesting applicable potential in oil and gas exploration of subsalt targets.