We have constructed a drift-diffusion based finite element simulation model to quantitatively analyze the two-dimensional oxygen vacancy distribution behavior in four-terminal memristive devices. Our analysis revealed that consideration of temperature dependent thermal activation for oxygen vacancy migration is critical to reproduce threshold behavior of the resistive switching. Systematic simulation results under various conditions are compared with experimental results and dominant factors for resistive switching based on two-dimensional redistribution of oxygen vacancies are discussed.