There is an equation for Green function of a half-plane: $$ \Delta G = \delta (x - x_{0})\delta(y - y_{0}), \quad G(x, 0) = 0.$$ Of course, the Green function which satisfies the boundary condition is represented as $$ G(x, y, x_{0}, y_{0}) = G_{0}(x, y, x_{0}, y_{0}) - G_{0}(x, y, x_{0}, -y_{0}), $$ where $G_{0}(x, y, x_{0}, y_{0})$ is the Green function for a plane, so I need to find a Green function for a plane.
After using Fourier transform (and setting $x_{0}, y_{0}$ to zero) I got $$-(k_{x}^{2} + k_{y}^2)\tilde {G_{0}} = \frac{1}{2 \pi} \Rightarrow \tilde G_{0} = -\frac{1}{2 \pi}\frac{1}{k_{x}^{2} + k_{y}^{2}}.$$
But after that I got divergent integral. Did I make a mistake?