The application discloses a VTI medium double-way wave migration imaging method based on De Wolf approximation, and comprises the following steps: firstly, dividing a model into multiple thin plates; then, based on a seismic source, sequentially calculating and obtaining a positive transmission down-going migration
wave field and a positive transmission first up-going migration
wave field of each thin plate from a minimum depth to a
maximum depth; then, calculating a positive transmission up-going migration
wave field of each thin plate at different positions from the
maximum depth to the minimum depth and adding the reserved positive transmission first up-going migration wave field, so as to finally obtain the positive transmission up-going migration wave field of all the thin plates; then, based on seismic records, obtaining reverse transmission down-going and up-going migration wave fields according to the above calculation process; finally, obtaining all imaging results by using a cross-correlation
imaging condition. Not only the down-going migration wave field information can be used, but also the up-going migration wave field information can be used; a double-way wave migration operator is formed by two one-way migration operators, the influence of anisotropic parameters on migration imaging results can be corrected, and finally the migration imaging precision of steeply dipping structures is effectively improved.