The work is devoted to one of the approaches of wells modeling within numerical oil reservoir simulation. The approach can be consider as fictitious domain method at mixed finite element approximation, which is used for non-stationary filtration processes of two phase fluid in Bukley-Leverett problem. The numerical results are compared with the results for the problem with usual Neumann conditions at the wells boundaries.