This paper presents a reverse-time migration algorithm as the conjugate of the forward time-extrapolation modeling. The conjugate operator is obtained by formulating the forward operator explicitly and by transposing the forward operator. The algorithm is tested by a numerical experiment for a 1-D layered model. A method to suppress artifacts, which are caused by the intrinsic non-unitary property of the operator, is introduced using the least-squares optimization technique.