<p style='text-indent:20px;'>In this paper, we establish the <inline-formula><tex-math id="M1">\begin{document}$ W^{1,p} $\end{document}</tex-math></inline-formula> estimates for solutions of second order elliptic problems with drift terms in bounded Lipschitz domains by using a real variable method. For scalar equations, we prove that the <inline-formula><tex-math id="M2">\begin{document}$ W^{1,p} $\end{document}</tex-math></inline-formula> estimates hold for <inline-formula><tex-math id="M3">\begin{document}$ \frac{3}{2}-\varepsilon<p<3+\varepsilon $\end{document}</tex-math></inline-formula> for <inline-formula><tex-math id="M4">\begin{document}$ d\geq3 $\end{document}</tex-math></inline-formula>, and the range for <inline-formula><tex-math id="M5">\begin{document}$ p $\end{document}</tex-math></inline-formula> is sharp. For elliptic systems, we prove that the <inline-formula><tex-math id="M6">\begin{document}$ W^{1,p} $\end{document}</tex-math></inline-formula> estimates hold for <inline-formula><tex-math id="M7">\begin{document}$ \frac{2d}{d+1}-\varepsilon<p<\frac{2d}{d-1}+\varepsilon $\end{document}</tex-math></inline-formula> under the assumption that the Lipschitz constant of the domain is small.</p>
$ W^{1, p} $ estimates for elliptic problems with drift terms in Lipschitz domains