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.