The paper proposes a method for constructing the surrogate models for two-phase flow based on a combination of finite-difference solutions and fine-grid spatial approximation. The method provides the approximate models for solving optimal control the operating parameters of field development.