The invention discloses a numerical 
simulation method and application of an 
enhanced geothermal system considering a working medium storage process, and the method comprises the following solving steps at any moment: solving a 
continuity equation and a 
momentum equation by taking an unsteady-state item in the 
continuity equation as a known value, and calculating a speed field and a pressure fieldof a working medium in a 
thermal reservoir at the moment; calculating a speed field, a pressure field and a temperature field of the working medium in the 
injection well shaft at the moment; calculating the temperature field of the 
thermal reservoir at the moment, and estimating the temperature field of the 
thermal reservoir at the next moment; calculating a speed field, a pressure field and a temperature field of the working medium in the production well shaft at the moment; calculating the average total 
stress field of the thermal reservoir at the moment, and estimating the average total 
stress field at the next moment; calculating the 
porosity and permeability distribution of the thermal reservoir at the moment, and estimating the 
porosity distribution at the next moment; and judging whether each 
physical field is converged or not, if not, repeating the steps, and if so, stopping. According to the invention, the 
thermal recovery performance and the working medium storage capacity inthe thermal reservoir can be simulated at the same time.