CN112632875B - Enhanced geothermal system heat-taking working medium and determination method thereof - Google Patents

Enhanced geothermal system heat-taking working medium and determination method thereof Download PDF

Info

Publication number
CN112632875B
CN112632875B CN202110009001.2A CN202110009001A CN112632875B CN 112632875 B CN112632875 B CN 112632875B CN 202110009001 A CN202110009001 A CN 202110009001A CN 112632875 B CN112632875 B CN 112632875B
Authority
CN
China
Prior art keywords
heat
working medium
temperature
rock
helium
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Expired - Fee Related
Application number
CN202110009001.2A
Other languages
Chinese (zh)
Other versions
CN112632875A (en
Inventor
刘雪玲
王源铭
汪健生
牛锦涛
付伟娟
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Tianjin University
Original Assignee
Tianjin University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Tianjin University filed Critical Tianjin University
Priority to CN202110009001.2A priority Critical patent/CN112632875B/en
Publication of CN112632875A publication Critical patent/CN112632875A/en
Application granted granted Critical
Publication of CN112632875B publication Critical patent/CN112632875B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • G06F30/28Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00Details relating to the application field
    • G06F2113/08Fluids
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/08Thermal analysis or thermal optimisation

Abstract

The invention provides a method for determining a heat-taking working medium of an enhanced geothermal system, which comprises the following steps: establishing a thermal-hydraulic coupling model of the hot dry rock based on physical property parameters of the deep hot dry rock thermal storage obtained by a rock physical experiment; constructing a cubic domain of a rock matrix, constructing a fracture channel with sine fluctuation morphology in the middle of the rock matrix by adopting a parametric surface method, determining the height of the fracture channel, and allowing fluid to flow into the fracture channel along the horizontal direction for heat exchange; simulation is carried out by using Commol Multiphysics to obtain different production temperatures of different working media under the same recharge flow and temperature. The invention determines that the heat-taking working medium with higher production temperature can be obtained according to the method as follows: supercritical carbon dioxide and helium are mixed according to the mass fraction of 9: 1 the mixed working medium obtained by mixing is a heat-taking working medium.

Description

Enhanced geothermal system heat-taking working medium and determination method thereof
Technical Field
The invention relates to a heat-taking working medium of an enhanced geothermal system for developing and utilizing geothermal energy of hot dry rocks, which is formed by mixing two working media, can realize higher production temperature and belongs to the field of application of geothermal energy of hot dry rocks.
Background
The utilization of geothermal energy has great significance for making up the uneven distribution of the traditional energy sources, relieving the problem of environmental pollution caused in the use process of the traditional energy sources and adjusting the energy source structure. Geothermal energy is a clean renewable energy source and has abundant reserves. The total available energy of geothermal heat is 5.4 x 10 estimated from the World Energy Council (WEC)27J, if only 1% of geothermal energy is utilized, the energy demand of 20000 years can be met. Compared with hydrothermal geothermal resources, the dry hot rock has low porosity, weak permeability, high temperature and rich reserves and is usually positioned3-10km underground is a clean energy which is currently researched and developed intensively in countries all over the world.
An enhanced geothermal system (EGS for short) is an effective means for developing and utilizing dry-hot rocks at present. A hot dry rock zone with a high temperature is selected and two or more parallel or near parallel wells are drilled in the zone as production wells and recharge wells. And constructing an artificial fracture network structure in a selected high-temperature rock region by means of hydraulic fracturing, wherein the artificial fracture network structure is called artificial heat storage. In the actual production process, heat exchange fluid is refilled into underground artificial heat storage through a refilling well for heat exchange, high-temperature fluid after heat exchange is pumped back to the ground surface through the production well and conveyed to a power plant for power generation and heat utilization, and low-temperature fluid after heat exchange is used for being refilled to the underground through the refilling well for heat exchange.
In order to efficiently extract heat in the dry-hot rock stratum, a working medium with excellent heat extraction property can be searched besides the improvement of a recharging mode. The enhanced geothermal system firstly adopts water as a heat-taking working medium. In 2000, enhanced geothermal systems with supercritical carbon dioxide as the heat-extracting working medium began to emerge. Compared with a water-based enhanced geothermal system, the supercritical carbon dioxide can reduce the pump power consumption required for pumping the supercritical carbon dioxide, and the interconnected fracture network of a plurality of secondary branches is easier to generate by using the supercritical carbon dioxide as a heat-taking working medium under the condition of relatively low fracturing pressure, so that a rock matrix generates a plurality of flow paths. Supercritical carbon dioxide has physical properties similar to the density and gas-like viscosity of liquids. Helium has small density, low viscosity and good heat conductivity, and is a working medium with very excellent heat exchange performance. Within the range known at present, helium is not found as a heat taking working medium of the enhanced geothermal system, and most of the current research is uniform working medium, and the excellent heat exchange characteristic of mixing two working media is not provided.
Reference documents:
[1]WANG J S,JIAO Y,LIU X L.Heat transfer and flow characteristics in a rectangular channel with small scale vortex generators.International journal of heat and mass transfer,2019,138:208–25.
[2]ZHANG L,JIANG P X,WANG Z,XU R N.Convective heat transfer of supercritical CO2 in a rock fracture for enhanced geothermal systems[J].Applied thermal engineering,2017,115:923-936.
disclosure of Invention
The invention aims to provide a method for determining a heat-taking working medium of an enhanced geothermal system and the determined heat-taking working medium. The invention mixes the supercritical carbon dioxide and the helium gas according to different mass fractions, provides a method for simulating the production temperature of different mixed working media and obtaining a new working medium, and the new working medium obtained by the invention can realize higher production temperature. The technical scheme is as follows:
a method for determining a heat-taking working medium of an enhanced geothermal system comprises the following steps:
1) establishing a thermal-hydraulic coupling model of the hot dry rock based on physical parameters of the deep hot dry rock thermal storage obtained by a rock physical experiment, wherein the physical parameters comprise temperature, permeability, density, thermal conductivity and specific heat;
2) constructing a cubic domain of a rock matrix, constructing a fracture channel with sine fluctuation morphology in the middle of the rock matrix by adopting a parametric surface method, determining the height of the fracture channel, and allowing fluid to flow into the fracture channel along the horizontal direction for heat exchange;
3) setting the density, specific heat, heat conductivity coefficient, initial temperature and initial pressure of the rock matrix; setting the temperature, the recharge flow and the outlet pressure of the recharge fluid; the upper wall surface and the lower wall surface of the cubic area adopt constant temperature boundary conditions; the peripheral wall surfaces adopt heat insulation boundary conditions, and the fluid adopts non-slip speed and non-step temperature boundaries on the upper wall surface and the lower wall surface of the fracture channel;
4) in the working range of underground heat extraction of the hot dry rock, the helium and the supercritical carbon dioxide are both in a supercritical state, an interpolation function changing along with the temperature and the pressure is generated according to physical parameters of a heat extraction working medium, and the density, the viscosity, the specific heat capacity and the heat conductivity coefficient of the supercritical carbon dioxide and the helium at different temperatures and pressures are respectively obtained under mixed working media with different helium mass fractions;
5) carrying out simulation by using a laminar flow module and a fluid and solid heat transfer module in Commol Multiphysics to obtain different production temperatures of different working media under the same recharge flow and temperature;
6) and determining that the heat-taking working medium with higher production temperature can be obtained.
The enhanced geothermal system heat-taking working medium determined by the method is characterized in that supercritical carbon dioxide and helium are mixed according to the mass fraction of 9: 1 the mixed working medium obtained by mixing is a heat-taking working medium.
The novel heat-taking working medium of the enhanced geothermal system provided by the invention has a better heat exchange effect, can obviously improve the production temperature of the system, improves the heat extraction capability of the working medium, and obviously improves the heat-taking effect when only supercritical carbon dioxide is used as the heat-taking working medium.
Drawings
FIG. 1 is a physical model of the present invention
FIG. 2 is a diagram of model validation
FIG. 3 is a trellis diagram
FIG. 4 shows production temperatures (10%, 20%, 30%, 40%, 50%, 60%, 70%, 80%, 90% in the legend represents the mass fraction of helium gas)
FIG. 5 physical Property parameters (10%, 20%, 30% in the drawing represent the mass fraction of helium gas)
Detailed Description
The method is based on the physical parameters of the deep dry hot rock thermal reservoir obtained by rock physics experiments, including temperature, density, thermal conductivity, specific heat and the like, and a dry hot rock thermal-hydraulic coupling model is established. And determining the flowing and spreading rule of heat between the hot dry rock fracture and the heat taking fluid through simulation.
Firstly, a cube of 50mm multiplied by 50mm is constructed in Commol Multiphysics as a rock matrix, and a parametric surface method is adopted in the middle of the rock matrix to construct a fracture channel with sine wave morphology, wherein the height of the fracture channel is 0.2 mm. The fluid flows into the fracture channel along the horizontal direction to exchange heat. The specific morphological parameters of the fracture channel are as follows:
Figure BDA0002884558180000021
(z1=24.9mm,z2=25.1mm)
Figure BDA0002884558180000022
the physical model of the underground heat-extraction part of the enhanced geothermal system adopted by the invention is shown in figure 1.
The present invention makes the following assumptions: the reservoir rock mass is homogeneous, and the physical property of the reservoir rock mass is irrelevant to temperature and pressure; secondly, deformation of reservoir rock mass is not considered; thirdly, viscous dissipation in the fluid flow process is ignored; and fourthly, adopting a non-slip boundary condition for the fracture channel boundary. Based on the above assumptions, the mathematical model constructed by the invention is as follows:
continuity equation of the fluid:
Figure BDA0002884558180000031
where ρ isfFluid density, kg/m3(ii) a τ is time, s; u is the fluid velocity vector, m/s.
Momentum governing equation of fluid:
Figure BDA0002884558180000032
wherein, mufIs hydrodynamic viscosity, Pa · s; p is the fluid pressure, Pa.
Energy control equation of fluid:
Figure BDA0002884558180000033
wherein, TfIs the fluid temperature, K; lambda [ alpha ]fIs the thermal conductivity of the fluid, W/(m.K); c. CfJ/(kg. K) is the specific heat capacity of the fluid.
Energy equation in rock matrix:
Figure BDA0002884558180000034
wherein, TsIs the rock matrix temperature, K; lambda [ alpha ]sThe thermal conductivity of the rock matrix is W/(m.K); rhosDensity of rock matrix, kg/m3;csThe specific heat capacity of the rock matrix is J/(kg. K).
The initial parameters selected by the invention are as follows: the density of the rock matrix is 2700kg/m3The specific heat is 1250J/(kg.K), the heat conductivity coefficient is 3W/(m.K), the initial temperature of the rock matrix is 473.15K, and the initial pressure is 20 MPa; the temperature of the recharging fluid is 343.15K, the recharging flow rate is 0.35kg/h, and the outlet pressure is 20 MPa; calculating the upper wall surface and the lower wall surface of the domain by adopting a constant temperature boundary condition; the peripheral wall surface adopts an adiabatic boundary condition, and the fluid adopts a non-slip speed and non-step temperature boundary on the upper wall surface and the lower wall surface of the fracture channel.
In order to ensure the calculation accuracy, the invention divides five grids with different sizes and numbers, and performs Richardson extrapolation method[1]The verification is carried out, the number of the finally adopted grids is 190259, the grids of the crack part are locally encrypted, and a physical model meshing schematic diagram is shown in fig. 2.
In order to ensure the accuracy and reliability of the calculation result, the simulation result and the reference document of the invention[2]The experimental results in (1) are compared: when the recharging parameters and the physical parameters are consistent with the literature, the production temperature is 0.35kg/h and the maximum deviation of the experimental result is 0.1% when the fluid pressure is 8MPa and the inlet temperature is 340K, which shows that the calculation result of the research is accurate and reliable.
The mixed working medium selected in the invention is helium and carbon dioxide, and both are in a supercritical state. The mass fractions of helium gas in the mixed working medium are respectively 10%, 20%, 30%, 40%, 50%, 60%, 70%, 80% and 90%, and the physical parameters of the 9 mixed working media are generated by REFPROP to generate interpolation functions changing with temperature and pressure, so as to respectively obtain the physical parameters of the mixed working medium such as density, viscosity, specific heat capacity, thermal conductivity coefficient and the like at different temperatures and pressures. The interpolation functions of the physical property parameters of the 9 mixed working media and the supercritical carbon dioxide and helium are respectively led into Comsol Multiphysics for simulation, so as to obtain different production temperatures of different working media under the same recharge flow and temperature, namely the temperature of the high-temperature working media extracted through the production well, as shown in FIG. 4.
Based on the enhanced geothermal system heat storage model and various mixed working media provided by the invention, the production temperature of helium under different mass fractions is researched to evaluate the heat extraction capability of different mixed working media, and the result shows that: when the mass fraction of helium is 10%, the production temperature of the mixed working medium is highest; the production temperature of pure helium is the lowest; in the mixed working medium, the production temperature is continuously reduced along with the increase of the helium proportion; when the mass fraction of the helium is 10% and 20%, the production temperature of the mixed working medium is higher than that of the supercritical carbon dioxide. As can be seen from fig. 5, in a certain temperature and pressure range, the specific heat capacity of helium is significantly higher than that of supercritical carbon dioxide, the specific heat capacity of the mixed working medium with helium mass fractions of 10% and 20% is lower than that of supercritical carbon dioxide, and when the helium mass fraction is 10%, the specific heat capacity of the mixed working medium is the lowest, so when the supercritical carbon dioxide and helium have a mass fraction of 9: 1, the highest production temperature can be realized, and a better heat exchange effect is achieved.
Therefore, the novel heat-taking working medium of the enhanced geothermal system provided by the invention has a better heat-exchanging effect, can obviously improve the production temperature of the system, improve the heat-extracting capacity of the working medium and obviously improve the heat-taking effect when only supercritical carbon dioxide is used as the heat-taking working medium.

Claims (2)

1. A method for determining a heat-taking working medium of an enhanced geothermal system comprises the following steps:
1) establishing a thermal-hydraulic coupling model of the hot dry rock based on physical parameters of the deep hot dry rock thermal storage obtained by a rock physical experiment, wherein the physical parameters comprise temperature, permeability, density, thermal conductivity and specific heat;
2) constructing a cubic domain of a rock matrix, constructing a fracture channel with sine fluctuation morphology in the middle of the rock matrix by adopting a parametric surface method, determining the height of the fracture channel, and allowing fluid to flow into the fracture channel along the horizontal direction for heat exchange;
3) setting the density, specific heat, heat conductivity coefficient, initial temperature and initial pressure of the rock matrix; setting the temperature, flow and outlet pressure of the recharging fluid; the upper wall surface and the lower wall surface of the cubic area adopt constant temperature boundary conditions; the peripheral wall surfaces adopt heat insulation boundary conditions, and the fluid adopts non-slip speed and non-step temperature boundaries on the upper wall surface and the lower wall surface of the fracture channel;
4) in the working range of underground heat extraction of the dry hot rock, helium and supercritical carbon dioxide are both in a supercritical state, an interpolation function changing along with temperature and pressure is generated according to physical property parameters of a heat extraction working medium, and the density, viscosity, specific heat capacity and heat conductivity coefficient of the supercritical carbon dioxide and the helium at different temperatures and pressures are respectively obtained under mixed working media with different helium mass fractions;
5) carrying out simulation by using a laminar flow module and a fluid and solid heat transfer module in Commol Multiphysics to obtain different production temperatures of different working media under the same recharge flow and temperature;
6) and determining that the heat-taking working medium with higher production temperature can be obtained.
2. The enhanced geothermal system heat extraction working medium determined by the method of claim 1, wherein supercritical carbon dioxide and helium are mixed according to the mass fraction of 9: 1 the mixed working medium obtained by mixing is a heat-taking working medium.
CN202110009001.2A 2021-01-05 2021-01-05 Enhanced geothermal system heat-taking working medium and determination method thereof Expired - Fee Related CN112632875B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110009001.2A CN112632875B (en) 2021-01-05 2021-01-05 Enhanced geothermal system heat-taking working medium and determination method thereof

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110009001.2A CN112632875B (en) 2021-01-05 2021-01-05 Enhanced geothermal system heat-taking working medium and determination method thereof

Publications (2)

Publication Number Publication Date
CN112632875A CN112632875A (en) 2021-04-09
CN112632875B true CN112632875B (en) 2022-06-10

Family

ID=75290767

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110009001.2A Expired - Fee Related CN112632875B (en) 2021-01-05 2021-01-05 Enhanced geothermal system heat-taking working medium and determination method thereof

Country Status (1)

Country Link
CN (1) CN112632875B (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112065521A (en) 2020-09-16 2020-12-11 天津大学 Based on CO2Mixed working medium supercharging heat absorption transcritical circulation hot dry rock geothermal power generation model

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102384046A (en) * 2011-06-24 2012-03-21 清华大学 Energy conversion system used in intensified geothermal system with CO2 as working medium
CN106870043A (en) * 2017-04-18 2017-06-20 长沙紫宸科技开发有限公司 The change system and method for carbon dioxide recycle generating are realized using geothermal energy
CN107100605A (en) * 2017-04-21 2017-08-29 中国石油大学(北京) A kind of method that dual horizontal well circulation supercritical carbon dioxide develops hot dry rock
CN107905778A (en) * 2017-10-19 2018-04-13 中国石油大学(华东) Supercritical CO2The enhanced geothermal system experimental provision of fluid fracturing and method
CN112065521A (en) * 2020-09-16 2020-12-11 天津大学 Based on CO2Mixed working medium supercharging heat absorption transcritical circulation hot dry rock geothermal power generation model

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102384046A (en) * 2011-06-24 2012-03-21 清华大学 Energy conversion system used in intensified geothermal system with CO2 as working medium
CN106870043A (en) * 2017-04-18 2017-06-20 长沙紫宸科技开发有限公司 The change system and method for carbon dioxide recycle generating are realized using geothermal energy
CN107100605A (en) * 2017-04-21 2017-08-29 中国石油大学(北京) A kind of method that dual horizontal well circulation supercritical carbon dioxide develops hot dry rock
CN107905778A (en) * 2017-10-19 2018-04-13 中国石油大学(华东) Supercritical CO2The enhanced geothermal system experimental provision of fluid fracturing and method
CN112065521A (en) * 2020-09-16 2020-12-11 天津大学 Based on CO2Mixed working medium supercharging heat absorption transcritical circulation hot dry rock geothermal power generation model

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
基于超临界二氧化碳-氦气混合工质的适应性流道回热器性能实验研究;程文龙等;《中国科学:技术科学》;20200925;第1-8页 *

Also Published As

Publication number Publication date
CN112632875A (en) 2021-04-09

Similar Documents

Publication Publication Date Title
Wang et al. Simulation of heat extraction from CO2-based enhanced geothermal systems considering CO2 sequestration
Zhang et al. System thermodynamic performance comparison of CO2-EGS and water-EGS systems
Ma et al. Study on the effect of well layout schemes and fracture parameters on the heat extraction performance of enhanced geothermal system in fractured reservoir
Wang et al. Improvement of energy efficiency for an open-loop surface water source heat pump system via optimal design of water-intake
Fard et al. Numerical simulation of a supercritical CO2 geothermosiphon
Zhang et al. Thermal characteristics analysis with local thermal non-equilibrium model during liquid nitrogen jet fracturing for HDR reservoirs
Ghoreishi-Madiseh et al. A transient natural convection heat transfer model for geothermal borehole heat exchangers
Liu et al. Comparative investigation on the heat extraction performance of an enhanced geothermal system with N2O, CO2 and H2O as working fluids
CN112632875B (en) Enhanced geothermal system heat-taking working medium and determination method thereof
Xie et al. The role of fracture networks randomness in thermal utilization of enhanced geothermal system
Liu et al. Convection heat transfer of supercritical CO2 in a single fracture in enhanced geothermal systems
Xiao et al. An investigation on spontaneous imbibition in tree-like branching network with fractal roughened surfaces
Wang et al. Influences of reservoir heterogeneity and anisotropy on CO 2 sequestration and heat extraction for CO 2-based enhanced geothermal system
Jiayan et al. Effects of variable thermophysical properties of water on the heat extraction of an enhanced geothermal system: A numerical case study
Śliwa et al. Theoretical model of borehole heat exchanger
Liu et al. The investigation on a hot dry rock compressed air energy storage system
Marzouk et al. Evaluating the effects of bifurcation angle on the performance of a novel heat exchanger based on contractual theory
Zhou et al. Simulation of thermal transport in aquifer: a GWHP system in Chengdu, China
Chávez et al. Polygeneration study of low-to-medium enthalpy geothermal reservoirs in Mexico
Zhang et al. Effects of heat transfer fluid and boundary conditions on temperature field of enhanced geothermal system
CN207556007U (en) For the taper heat exchange sleeve structure of ground rock heat absorption
Kusaiynov et al. Study of the heat-transfer processes of tubular elements of ground heat exchangers
CN109989736A (en) A method of design EGS heat storage pressure break and modification scheme
Nithyadevi et al. Convection of Cu-water nanofluid in a partially active porous cavity with internal heat generation
CN103352769A (en) Method for acquiring energy in pressurized and reversed internal energy power system

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20220610

CF01 Termination of patent right due to non-payment of annual fee