Mathematical models of lunar and Earth shadow were established for lunar orbits. Also, the criterion was proposed to decide whether the satellite was in the shadow area. Golden section search optimization was used to search out the boundary values of the shadow area. In combination with specific cases, the optimal targeting lunar orbit was given by the genetic algorithm which was embedded in tabu search. The mathematical models of lunar orbit shadow given out in this article are simple enough to be used in onboard computer. Meanwhile, they can meet the engineering accuracy requirements. The embedded genetic algorithm has advantage of fast convergence without trapping in local optimum.