A non-local boundary value problem method for solving 2-D nonlinear heat equation backward in time is given. Some error estimates between the exact solution and its regularization approximation are provided and numerical examples show that the method works effectively.