The invention discloses a front fracturing well closing flowback numerical
simulation method considering reservoir damage, and belongs to the technical field of
petroleum development reservoir transformation. The method comprises the steps that initial parameters of a target oil reservoir are obtained, fracturing
fracture propagation simulation is conducted, and an embedded discrete fracture grid is constructed based on the fracture form obtained through
simulation; determining a function relationship between the reservoir damage coefficient and the
water saturation according to an experimental result; establishing a reservoir parameter dynamic calculation model; constructing a
mass conservation equation and carrying out differential
discretization to obtain a numerical
discretization solving model; and performing iterative solution on the model, dynamically updating reservoir
physical property parameters at each
time step, and finally obtaining a pressure field, a saturation field and
crude oil yield prediction in the whole fracturing process. The problem that reservoir damage caused by
fracturing fluid invasion cannot be dynamically reflected through traditional numerical simulation is solved, accurate simulation of the full-period process of preposed fracturing, well closing and flowback is achieved, and a reliable tool is provided for fracturing optimization and productivity prediction.