A wave equation migration method for receiver function imaging: 1. Theory