Elastic (seismic) wave stimulation is considered one of the unconventional enhanced oil recovery (EOR) methods. Increasing water quantity in the high permeability layer of a mature oil reservoir is highly challenging and can significantly decrease the ultimate recovery due to the reservoir heterogeneity. Using seismic waves can be considered low-cost, environmentally friendly, and illuminates the entire reservoir size compared to conventional EOR methods. A numerical model is developed by extending the Quintal approach for seismic attenuation due to wave-induced fluid flow (WIFF) to incorporate capillary pressure in partially saturated porous media and shift undrained boundary conditions to exclude external flow stress for drained boundary conditions. Therefore, the fluid distribution due to the capillary effect makes the developed finite element method (FEM) u-p model more widely applicable for oil recovery in mature reservoirs. A two-layer partially saturated media was subjected to compressive seismic stress at low frequency (3 Hz). The results indicated that the vertical displacement gradients of the bottom and upper layers decline with excitation time for both fully and partially saturated media. On the other hand, partially saturated pore pressure gradients of both the upper and bottom layers have higher amplitudes with excitation time than fully saturated pore pressure gradients due to the influence of capillar pressure. The cumulative crossflow oil volume for 180 days of continuous stimulation was 1176 bbl, 1032 bbl, and 648 bbl in low permeability layers: 200 md, 100 md, and 50 md, respectively. Therefore, the developed model has the potential for field-scale EOR applications. The study also suggests coupling elastic EOR with CO2 flooding to recover more oil due to increasing fluid mobility and relative permeability to oil in low-permeability reservoirs or tight formations.