Geothermal resource quantity calculation method and device

Abstract

The invention relates to the technical field of geothermal science, in particular to a geothermal resource quantity calculation method and device. The geothermal resource quantity calculation method comprises the following steps: acquiring geological parameters of a geothermal evaluation area; according to the geological parameters, constructing a three-dimensional geological model of the geothermal evaluation area; according to the three-dimensional geologic model, a heat reservoir of a geothermal evaluation area is selected; temperature numerical simulation in the heat transfer process is carried out on the geothermal evaluation area, temperature field distribution of the geothermal evaluation area is obtained, and the temperature field distribution is used for representing the temperature of each position point in the geothermal evaluation area; and according to the temperature of each position point in the thermal reservoir, performing integral calculation on the thermal reservoir to obtain the geothermal resource quantity of the thermal reservoir. According to the embodiment of the invention, the accuracy of geothermal resource quantity calculation can be improved.

Description

Technical Field

[0001] This specification relates to the technical field of geothermal science, and particularly to a method and device for calculating geothermal resource quantity. Background Art

[0002] The geothermal resource quantity is an important data for evaluating the potential of geothermal resources and is the basis for developing geothermal resources in a geothermal evaluation area. Currently, the methods for obtaining the geothermal resource quantity include methods such as the reservoir method. Taking the reservoir method as an example, it represents the overall thickness of the reservoir with the average thickness of the reservoir and the overall temperature of the reservoir with the average temperature of the reservoir, so as to calculate the geothermal resource quantity.

[0003] However, the temperature and thickness vary everywhere in the reservoir. For areas with large variations in thickness and temperature, there are large errors in the calculated geothermal resource quantity, thus making the credibility of the evaluation result of the geothermal resource potential in the geothermal evaluation area relatively low. Summary of the Invention

[0004] Embodiments of this specification provide a method and device for calculating geothermal resource quantity, for improving the accuracy of geothermal resource quantity calculation.

[0005] Embodiments of this specification provide a method for calculating geothermal resource quantity, including:

[0006] Obtaining geological parameters of a geothermal evaluation area;

[0007] Constructing a three-dimensional geological model of the geothermal evaluation area according to the geological parameters;

[0008] Selecting a reservoir layer of the geothermal evaluation area according to the three-dimensional geological model;

[0009] Performing a temperature numerical simulation on the heat transfer process of the geothermal evaluation area to obtain the temperature field distribution of the geothermal evaluation area, where the temperature field distribution is used to represent the temperature of each position point in the geothermal evaluation area;

[0010] Performing integral calculation on the reservoir layer according to the temperature of each position point in the reservoir layer to obtain the geothermal resource quantity of the reservoir layer.

[0011] Embodiments of this specification also provide a device for calculating geothermal resource quantity, including:

[0012] An obtaining unit, configured to obtain geological parameters of a geothermal evaluation area;

[0013] A constructing unit, configured to construct a three-dimensional geological model of the geothermal evaluation area according to the geological parameters;

[0014] A selecting unit, configured to select a reservoir layer of the geothermal evaluation area according to the three-dimensional geological model;

[0015] A simulation unit is configured to perform a numerical simulation of the temperature in the heat transfer process of the geothermal evaluation area, so as to obtain the temperature field distribution of the geothermal evaluation area, and the temperature field distribution is used to represent the temperature of each position point in the geothermal evaluation area.

[0016] A calculation unit is configured to perform an integral calculation on the heat reservoir according to the temperature of each position point in the heat reservoir, so as to obtain the geothermal resource amount of the heat reservoir.

[0017] An embodiment of this specification also provides a computer device, including a memory, a processor, and a computer program stored on the memory and executable on the processor. When the processor executes the computer program, the above-mentioned method for calculating the geothermal resource amount is implemented.

[0018] An embodiment of this specification also provides a computer-readable storage medium. The computer-readable storage medium stores a computer program, and when the computer program is executed by a processor, the above-mentioned method for calculating the geothermal resource amount is implemented.

[0019] An embodiment of this specification also provides a computer program product. The computer program product includes a computer program, and when the computer program is executed by a processor, the above-mentioned method for calculating the geothermal resource amount is implemented.

[0020] The technical solution of the embodiment of this specification can obtain the geological parameters of the geothermal evaluation area; can construct a three-dimensional geological model of the geothermal evaluation area according to the geological parameters; can select a heat reservoir in the geothermal evaluation area according to the three-dimensional geological model; can perform a numerical simulation of the temperature in the heat transfer process of the geothermal evaluation area to obtain the temperature field distribution of the geothermal evaluation area, and the temperature field distribution is used to represent the temperature of each position point in the geothermal evaluation area; can perform an integral calculation on the heat reservoir according to the temperature of each position point in the heat reservoir to obtain the geothermal resource amount of the heat reservoir. Thus, through the numerical simulation of the temperature in the heat transfer process of the geothermal evaluation area, the embodiment of this specification can obtain a high-precision temperature field distribution of the geothermal evaluation area. The temperature field distribution reflects the distribution of the temperature in the geothermal evaluation area with respect to the spatial position. By using the temperature field distribution of the geothermal evaluation area and performing an integral calculation on the temperature of the heat reservoir, the geothermal resource amount of the heat reservoir can be obtained. In this way, the changes in the temperature and thickness of the heat reservoir are fully considered, avoiding the errors caused by using the average thickness of the heat reservoir to represent the overall thickness of the heat reservoir and using the average temperature of the heat reservoir to represent the overall temperature of the heat reservoir, and being able to improve the accuracy of the calculated geothermal resource amount, thereby being conducive to more accurately evaluating the geothermal resource potential of the geothermal evaluation area. Description of the Drawings

[0021] To more clearly illustrate the technical solutions in the embodiments of this specification or the prior art, the following will briefly introduce the accompanying drawings required for the description of the embodiments or the prior art. The accompanying drawings in the following description are only some embodiments recorded in this specification. For those of ordinary skill in the art, without creative efforts, other drawings can also be obtained based on these drawings.

[0022] Figure 1 It is a schematic flowchart of the geothermal resource quantity calculation method in the embodiments of this specification;

[0023] Figure 2 It is a topographic map of the geothermal evaluation area in the embodiments of this specification;

[0024] Figure 3 It is a schematic diagram of the three-dimensional geological model of the geothermal evaluation area in the embodiments of this specification;

[0025] Figure 4 It is a schematic diagram of the three-dimensional geological model after grid processing in the embodiments of this specification;

[0026] Figure 5 It is a schematic diagram of the temperature field distribution of the geothermal evaluation area in the embodiments of this specification;

[0027] Figure 6 It is a schematic diagram of the location of the geothermal resource potential advantage area in the geothermal evaluation area in the embodiments of this specification;

[0028] Figure 7 It is a schematic diagram of the structure of the geothermal resource quantity calculation device in the embodiments of this specification. Detailed implementation manners

[0029] The following will clearly and completely describe the technical solutions in the embodiments of this specification in conjunction with the accompanying drawings in the embodiments of this specification. Obviously, the described embodiments are only a part of the embodiments of this specification, rather than all of them. The specific embodiments described here are only used to explain the present disclosure, rather than limiting the present disclosure. Based on the described embodiments of the present disclosure, all other embodiments obtained by those of ordinary skill in the art belong to the scope of protection of the present disclosure. In addition, relational terms such as "first" and "second" are only used to distinguish one entity or operation from another entity or operation, and do not necessarily require or imply any actual relationship or order between these entities or operations.

[0030] Please refer to Figure 1 . The embodiments of this specification provide a geothermal resource quantity calculation method. The geothermal resource quantity calculation method can be applied to a computer device and specifically may include the following steps.

[0031] Step 11: Obtain the geological parameters of the geothermal evaluation area.

[0032] In some embodiments, the geothermal evaluation area is a spatial area for evaluating the potential of geothermal resources. The geological parameters may include geological structure data, temperature data, and rock physical property data. The geological structure data is used to represent the geological structure of the geothermal evaluation area. The geological structure data may include elevation data, stratigraphic distribution data, and fault distribution data, etc. The elevation data is used to represent the terrain undulation and elevation of the geothermal evaluation area. The stratigraphic distribution data is used to represent the stratigraphic structure, lithology distribution, and stratigraphic thickness of the geothermal evaluation area. The fault distribution data is used to represent the structural characteristics such as the position, strike, and dip of the fault zone. The temperature data includes Curie depth data, Curie temperature data, and surface temperature data. The Curie surface is the Curie isothermal surface in the earth's crust. The Curie depth data is used to represent the depth of the Curie surface. The Curie temperature data is used to represent the temperature of the Curie surface. The surface temperature data is used to represent the temperature of the surface. The rock physical property data is used to represent the thermophysical properties of the rocks in the geothermal evaluation area, and may include, for example, thermal conductivity, specific heat capacity, density, porosity, permeability, etc.

[0033] In some embodiments, an area that needs to be evaluated for geothermal resources can be selected as the geothermal evaluation area; the geological parameters of the geothermal evaluation area can be obtained. For example, an area within a set longitude and latitude range that needs to be evaluated for geothermal resources can be selected as the geothermal evaluation area. Figure 2 FIG. is a topographic map of the geothermal evaluation area. The geothermal evaluation area includes multiple volcanoes. At the same time, a large number of hot springs and geothermal wells are distributed in the geothermal evaluation area. For example, hot spring group 1, hot spring group 2, hot spring group 3, hot spring group 4, geothermal well 1, geothermal well 2, geothermal well 3, geothermal well 4, etc. are distributed in the geothermal evaluation area, indicating good geothermal resources. However, due to the complex geological conditions and the rough current geothermal field, the amount of geothermal resources in the geothermal evaluation area is not clear.

[0034] The geological parameters of the geothermal evaluation area can be obtained through satellite remote sensing, seismic exploration, meteorological station data, field measurement, etc. For example, the elevation data can be obtained through satellite remote sensing or field measurement. The surface temperature data can be obtained through meteorological station data, satellite remote sensing, or field measurement. The stratigraphic distribution data and the fault distribution data can be obtained through seismic exploration. The Curie depth data and the Curie temperature data can be obtained through geothermal gradient measurement or geothermal flow measurement. The rock physical property data can be obtained through laboratory tests.

[0035] Step 12: Construct a three-dimensional geological model of the geothermal evaluation area according to the geological parameters.

[0036] In some embodiments, a three-dimensional geological model framework with the Curie surface as the bottom boundary of the geothermal evaluation area can be constructed based on altitude data, formation distribution data, fault distribution data, and Curie surface depth data.

[0037] The three-dimensional geological model framework is used to represent the surface, each formation interface, and the Curie surface of the geothermal evaluation area.

[0038] For example, altitude data can be used to generate a surface model of the geothermal evaluation area as the upper boundary of the three-dimensional geological model framework; formation distribution data can be used to generate the top and bottom interfaces of each formation; the top and bottom interfaces of each formation can be connected to obtain a three-dimensional model framework of the formation; fault distribution data can be used to generate a three-dimensional model framework of the fault zone; the three-dimensional model framework of the fault zone can be embedded into the three-dimensional model framework of the formation. The Curie surface model can be generated based on the Curie surface depth data as the bottom boundary of the three-dimensional geological model framework. The surface model, the three-dimensional model framework of the formation, the three-dimensional model framework of the fault zone, and the Curie surface model can be fused to obtain the three-dimensional geological model framework.

[0039] Figure 3 A three-dimensional geological model framework of the geothermal evaluation area is shown.

[0040] In some embodiments, surface temperature data, Curie surface temperature data, and rock physical property data can be incorporated into the three-dimensional geological model framework to obtain a three-dimensional geological model. For example, the surface temperature data can be assigned to the top layer of the three-dimensional geological model framework; the Curie surface temperature data can be assigned to the bottom layer of the three-dimensional geological model framework; the rock physical property data can be assigned to the corresponding formations in the three-dimensional geological model framework, thereby obtaining the three-dimensional geological model. In this way, the three-dimensional geological model can be used to represent the surface temperature, Curie surface temperature, and rock physical property data of each formation in the geothermal evaluation area.

[0041] Step 13: Select a heat reservoir in the geothermal evaluation area according to the three-dimensional geological model.

[0042] In some embodiments, the heat reservoir is the carrier of geothermal resources in the geothermal evaluation area, capable of storing and conducting geothermal fluids (such as water or steam), and capable of characterizing the geothermal resource potential of the geothermal evaluation area. The heat reservoir is a formation with good physical property data, such as a formation with high porosity, high permeability, and high thermal conductivity. The number of the heat reservoirs can be one or more.

[0043] In some embodiments, a formation that meets the set physical property conditions can be selected in the three-dimensional geological model as the heat reservoir in the geothermal evaluation area.

[0044] The set physical property conditions may include a first physical property data range. The physical property data within the first physical property data range is greater than or equal to the first set physical property data. For example, the physical property data includes porosity, permeability, and thermal conductivity. The first set physical property data includes a first porosity threshold, a first permeability threshold, and a first thermal conductivity threshold. The porosity of the physical property data within the first physical property data range is greater than or equal to the first porosity threshold, the permeability of the physical property data within the first physical property data range is greater than or equal to the first permeability threshold, and the thermal conductivity of the physical property data within the first physical property data range is greater than or equal to the first thermal conductivity threshold. The formation corresponding to the physical property data within the first physical property data range can be selected in the three-dimensional geological model as the heat reservoir of the geothermal evaluation area. For example, the formation with porosity greater than or equal to the first porosity threshold, permeability greater than or equal to the first permeability threshold, and thermal conductivity greater than or equal to the first thermal conductivity threshold can be selected in the three-dimensional geological model as the heat reservoir of the geothermal evaluation area.

[0045] Further, the set physical property conditions may further include a second physical property data range. The physical property data within the second physical property data range is less than or equal to the second set physical property data. For example, the physical property data includes porosity, permeability, and thermal conductivity. The second set physical property data includes a second porosity threshold, a second permeability threshold, and a second thermal conductivity threshold. The porosity of the physical property data within the second physical property data range is less than or equal to the second porosity threshold, the permeability of the physical property data within the second physical property data range is less than or equal to the second permeability threshold, and the thermal conductivity of the physical property data within the second physical property data range is less than or equal to the second thermal conductivity threshold. In addition, the second set physical property data is less than the first set physical property data. For example, the physical property data includes porosity, permeability, and thermal conductivity. The first set physical property data includes a first porosity threshold, a first permeability threshold, and a first thermal conductivity threshold. The second set physical property data includes a second porosity threshold, a second permeability threshold, and a second thermal conductivity threshold. The second porosity threshold is less than the first porosity threshold. The second permeability threshold is less than the first permeability threshold. The second thermal conductivity threshold is less than the first thermal conductivity threshold.

[0046] One or more formations in the three-dimensional geological model whose corresponding physical property data is within the first physical property data interval can be selected as one or more candidate heat reservoirs. For each candidate heat reservoir, the formation covering the candidate heat reservoir can be obtained as the covering formation of the candidate heat reservoir. Among the one or more candidate heat reservoirs, the candidate heat reservoir whose physical property data of the corresponding covering formation is within the second physical property data interval can be selected as the heat reservoir of the geothermal evaluation area. For example, the candidate heat reservoir with a porosity less than or equal to the second porosity threshold, a permeability less than or equal to the second permeability threshold, and a thermal conductivity less than or equal to the second thermal conductivity threshold can be selected as the heat reservoir of the geothermal evaluation area.

[0047] Thus, the heat reservoir of the geothermal evaluation area is a formation with good physical property data, such as a formation with high porosity, permeability, and thermal conductivity; and there should be a formation with poor physical property data covering it, such as a formation with low porosity, low permeability, and low thermal conductivity covering it. The selected heat reservoir can better represent the geothermal evaluation area, and the selected heat reservoir can more accurately characterize the geothermal resource potential of the geothermal evaluation area.

[0048] Please refer to Figure 3 . The Permian and Triassic systems can be selected as the heat reservoirs of the geothermal evaluation area.

[0049] Step 14: Perform a temperature numerical simulation on the heat transfer process of the geothermal evaluation area to obtain the temperature field distribution of the geothermal evaluation area.

[0050] In some embodiments, the three-dimensional geological model can be meshed. According to the top and bottom temperature data of the three-dimensional geological model, the heat transfer process of the three-dimensional geological model can be calculated from bottom to top according to the grids in the three-dimensional geological model to obtain the temperature field distribution of the geothermal evaluation area. The temperature field distribution can provide a data basis for the calculation of geothermal resources. The temperature field distribution is used to represent the distribution of the temperature in the geothermal evaluation area with respect to the spatial position. For example, the temperature field distribution can represent the surface temperature, Curie temperature, and temperature data of each formation in the geothermal evaluation area. The temperature field distribution can include the present-day temperature field distribution. The present-day temperature field distribution is used to represent the distribution of the temperature in the current geothermal evaluation area with respect to the spatial position.

[0051] Figure 4 Shows the three-dimensional geological model after meshing. Figure 5 Shows the temperature field distribution of the geothermal evaluation area.

[0052] In some embodiments, a grid can be regarded as the granularity for temperature numerical simulation during the heat transfer process. The grid size can affect the accuracy of the temperature field distribution. The smaller the grid, the smaller the range of spatial positions represented by the grid, so that the changes in temperature at different spatial positions can be captured more precisely, improving the overall accuracy of the temperature field distribution. However, the computational cost during temperature numerical simulation is relatively large and the time spent is relatively long. The larger the grid, the larger the range of spatial positions represented by the grid, and thus the overall accuracy of the temperature field distribution is relatively small. However, the computational cost during temperature numerical simulation is relatively small and the time spent is relatively short.

[0053] The grid size can be determined according to the complexity of the geothermal evaluation area and the computational accuracy requirements of the temperature numerical simulation.

[0054] A first set of grid sizes can be provided. The first set of grid sizes can include multiple grid sizes. Each grid size corresponds to a computational accuracy. The computational accuracy corresponding to the grid size can be selected from a plurality of preset accuracy levels. For example, the plurality of accuracy levels can include high level, medium level, low level, etc. A target computational accuracy input by the user can be received. The target computational accuracy can be selected from the plurality of preset accuracy levels. The corresponding first target grid size can be selected from the first set of grid sizes according to the target computational accuracy. The first target grid size can be, for example, 20 meters.

[0055] The three-dimensional geological model can be meshed according to the first target grid size. Thus, in the meshed three-dimensional geological model, the sizes of all grids are the first target grid size. According to the top and bottom temperature data of the three-dimensional geological model, the heat transfer process of the three-dimensional geological model can be calculated from bottom to top according to the grids in the three-dimensional geological model, and the temperature field distribution of the geothermal evaluation area can be obtained.

[0056] Optionally, the geological conditions of the geothermal evaluation area are relatively complex. If the three-dimensional geological model is meshed using a unified first target grid size, it is impossible to well balance the accuracy of the temperature field distribution and the computational cost. Therefore, a target area can be selected in the three-dimensional geological model; the first target grid size of the target area can be optimized and adjusted. According to the top and bottom temperature data of the three-dimensional geological model, the heat transfer process of the three-dimensional geological model can be calculated from bottom to top according to the optimized grids in the three-dimensional geological model, and the temperature field distribution of the geothermal evaluation area can be obtained.

[0057] The target area is an area with relatively complex geological conditions within the geothermal evaluation area. For example, the target area may include the area where the fault zone is located. The geological complexity of the target area can be calculated; based on the geological complexity, the first target grid size of the target area can be optimized and adjusted. For example, a second set of grid sizes can be provided. The second set of grid sizes may include multiple grid sizes. Each grid size corresponds to a complexity interval. The complexity interval in which the geological complexity of the target area is located can be obtained; according to the obtained complexity interval, the corresponding second target grid size can be selected from the second set of grid sizes. The first target grid size of the target area can be adjusted to the second target grid size. In this way, the target area can be gridded using the second target grid size. The magnitude of the geological complexity of the target area is inversely correlated with the magnitude of the second target grid size. The greater the geological complexity of the target area, the smaller the second target grid size. The smaller the geological complexity of the target area, the larger the second target grid size. The second target grid size is smaller than the first target grid size. Thus, within the geothermal evaluation area, a smaller second target grid size can be used for gridding the target area, and a larger first target grid size can be used for gridding the other areas except the target area. This can achieve a better balance between the overall accuracy of the temperature field distribution and the computational amount.

[0058] The characteristic data of the target area can be obtained, and the geological complexity of the target area can be calculated using a set formula based on the characteristic data. For example, the target area may include the area where the fault zone is located. The characteristic data of the fault zone may include the number of branches and intersections, curvature, physical property change gradient, etc. of the fault zone. The geological complexity of the fault zone can be calculated using a set formula based on the number of branches and intersections, curvature, physical property change gradient, etc. of the fault zone.

[0059] In some embodiments, based on the finite element method, a temperature numerical simulation of the three-dimensional geological model of the geothermal evaluation area can be performed based on the heat conduction equation and the heat convection equation to obtain the current temperature field distribution of the geothermal evaluation area.

[0060] Heat transfer is a physical phenomenon in physics, which refers to the phenomenon of heat energy transfer caused by a temperature difference. In heat transfer, the change in the internal energy of an object can be measured by heat. Heat transfer can include three forms: heat conduction, heat radiation, and heat convection. As long as there is a temperature difference within an object or between objects, heat energy can be transferred from a high temperature to a low temperature in one or more of the above three ways. Within the terrestrial range, heat conduction and heat convection are the main forms of heat transfer. Therefore, according to the top and bottom temperature data of the three-dimensional geological model and the rock physical property data of each stratum, the heat transfer process of the three-dimensional geological model can be calculated from bottom to top according to the grids in the three-dimensional geological model to obtain the current temperature field distribution of the geothermal evaluation area.

[0061] Numerical simulation is also called computer simulation. Relying on an electronic computer, combined with the concepts of finite element or finite volume, through numerical calculation or image display methods, it aims to study engineering problems, physical problems, and various problems in nature.

[0062] The geothermal evaluation area has abundant basic geological data (such as fracture distribution characteristics, stratigraphic distribution characteristics, rock physical property data, Curie depth and temperature, surface temperature, and elevation data). These data cannot directly reflect the temperature field distribution in the geothermal evaluation area. However, by substituting these data into the heat conduction equation, the temperature distribution and terrestrial heat flow distribution in the geothermal evaluation area can be calculated, and then the temperature field distribution in the geothermal evaluation area can be obtained. Considering the extremely large amount of calculation, a three-dimensional geological model of the geothermal evaluation area can be established based on these data, and the temperature field distribution in the geothermal evaluation area can be calculated using the three-dimensional geological model.

[0063] For example, the heat conduction equation can include ρ represents density; c represents specific heat capacity; is the partial derivative of temperature T with respect to time t, representing the rate of change of temperature with time; represents the gradient operator; λ represents thermal conductivity; is the heat conduction term, representing the propagation of heat through heat conduction in space, and Q is the heat source term, representing the heat generated per unit volume.

[0064] For example, numerical methods (such as the conjugate gradient method, Newton iteration method) can be used to solve the heat conduction equation, and the temperature values of each grid in the three-dimensional geological model are gradually updated. Specifically, for example, the Curie temperature in the three-dimensional geological model can be used as the starting condition, and the surface temperature in the three-dimensional geological model can be used as the termination condition. For each grid, the grid below it can be used as the heat source term, and the thermal conductivity, specific heat capacity, and density of the formation where the grid is located can be obtained; the heat conduction equation can be solved numerically according to the heat source term, thermal conductivity, specific heat capacity, and density to obtain the temperature value of the grid. By gradually calculating the temperature values of each grid in ascending order from bottom to top through the heat conduction equation, the temperature field distribution in the geothermal evaluation area can be obtained. For example, the temperature field distribution in the geothermal evaluation area can include: surface temperature of 15 °C, temperature of 80 °C at a depth of 100 m, temperature of 150 °C at a depth of 300 m, temperature of 300 °C at a depth of 600 m, and Curie temperature of 450 °C.

[0065] Step 15: Integrate the heat reservoir according to the temperatures at each position point in the heat reservoir to obtain the geothermal resource amount of the heat reservoir.

[0066] In some embodiments, the temperature field distribution can be used to represent the temperatures of various position points in the geothermal evaluation area. The position points can include, for example, the grids in a three-dimensional geological model. The geothermal resource amount of the thermal reservoir can be obtained by performing an integral calculation on the thermal reservoir according to the temperatures of the various position points in the thermal reservoir. In this way, the variations in the temperature and thickness of the thermal reservoir are fully considered, avoiding the errors caused by representing the overall thickness of the thermal reservoir with the average thickness of the thermal reservoir and representing the overall temperature of the thermal reservoir with the average temperature of the thermal reservoir, improving the accuracy of the geothermal resource amount, and being conducive to accurately evaluating the geothermal resource potential of the geothermal evaluation area.

[0067] In some embodiments, the temperature data of each position point in the thermal reservoir can be extracted from the temperature field distribution; the geothermal resource amount of the thermal reservoir can be calculated using the triple integral formula according to the temperature data of each position point in the thermal reservoir. Thus, by utilizing the temperature field distribution of the geothermal evaluation area and performing an integral calculation on the temperatures of the various position points in the thermal reservoir, the geothermal resource amount of the thermal reservoir can be obtained.

[0068] For example, the formula Q = ∫∫∫ Ω Rdv can be used to calculate the geothermal resource amount of the thermal reservoir; where Q represents the geothermal resource amount of the thermal reservoir, Ω represents the volume of the thermal reservoir, ρ r represents the rock density of the thermal reservoir, c r represents the specific heat capacity of the rock of the thermal reservoir, represents the porosity of the rock of the thermal reservoir, t r represents the temperature of the position point in the thermal reservoir, t0 represents the surface temperature, ρ w represents the density of geothermal water, c w represents the specific heat capacity of water, dv represents the volume element of the thermal reservoir area. Among them, the values of ρ r , c r , t r , t0, ρ w , c w etc. can be preset. dv = dxdydz. x, y, z represent the coordinate data of the points in the thermal reservoir in a three-dimensional rectangular coordinate system. The three-dimensional rectangular coordinate system is a coordinate system established for the thermal reservoir.

[0069] For example, in step 14, the three-dimensional geological model can be meshed. Then the position points in the reservoir can include the grids in the reservoir. The temperature values of each grid in the thermal reservoir can be extracted from the temperature field distribution; the geothermal resource amount of the thermal reservoir can be calculated using the formula Q = ∫∫∫ Ω Rdv. t r represents the temperature value of the grid in the thermal reservoir. Specifically, for example, the geothermal resource amount of the thermal reservoir can be 5.42×1018 J.

[0070] In some embodiments, the geothermal resource potential of the geothermal evaluation area can be calculated according to the geothermal resource amount of the heat reservoir. For example, the geothermal resource amount of the heat reservoir can be directly used to characterize the geothermal resource potential of the geothermal evaluation area. For another example, according to the geothermal resource amount of the heat reservoir, the formula E = Q×η can be used to calculate the geothermal resource potential of the geothermal evaluation area. Q represents the geothermal resource amount of the heat reservoir, and η represents the recoverable coefficient. The recoverable coefficient η can be an empirical value, or it can also be calculated based on the physical property data of the rock in the heat reservoir. For example, if the permeability of the rock in the heat reservoir is high, the recoverable coefficient η can be taken as 0.2.

[0071] For example, the geothermal resource amount of the heat reservoir can be 5.42×10 18 J, and the geothermal power generation potential is about 2.11×10 11 kw·h. This indicates that the geothermal evaluation area has good geothermal power generation potential. The heat reservoir is an area with advantages in geothermal resource potential. Figure 6 The position of the area with advantages in geothermal resource potential in the geothermal evaluation area is shown, and it is speculated to be related to the underlying magma chamber system.

[0072] In some embodiments, multiple heat reservoirs can be selected through step 13. The geothermal resource amounts of multiple heat reservoirs can be obtained through step 15. Then, the geothermal resource potential of the geothermal evaluation area can be calculated according to the geothermal resource amounts of the multiple heat reservoirs. For example, the geothermal resource amounts of the multiple heat reservoirs can be added together, and the geothermal resource potential of the geothermal evaluation area can be determined according to the addition result. For example, the addition result can be directly used to characterize the geothermal resource potential of the geothermal evaluation area. For another example, according to the addition result, the formula E = Q×η can be used to calculate the geothermal resource potential of the geothermal evaluation area. Q represents the addition result of the geothermal resource amounts of the multiple heat reservoirs, and η is the recoverable coefficient.

[0073] In some embodiments, the temperature change data and thickness change data of the heat reservoir can be obtained; under the condition that the temperature change data is greater than or equal to the temperature change threshold, and / or, the thickness change data is greater than or equal to the thickness change threshold, it indicates that the thickness and / or temperature of the heat reservoir change greatly. Thus, through step 14, a numerical simulation of the temperature in the heat transfer process of the geothermal evaluation area can be carried out to obtain the temperature field distribution of the geothermal evaluation area; through step 15, according to the temperature at each position point in the heat reservoir, an integral calculation of the heat reservoir can be carried out to obtain the geothermal resource amount of the heat reservoir. Under the condition that the temperature change data is less than the temperature change threshold and the thickness change data is less than the thickness change threshold, it indicates that the thickness and temperature of the heat reservoir change little, so the formula Calculate the geothermal resource quantity of the thermal reservoir. Q represents the geothermal resource heat of the thermal reservoir; A represents the area of the thermal reservoir; Z represents the average thickness of the thermal reservoir; ρ r represents the rock density of the thermal reservoir; c r represents the specific heat capacity of the rock of the thermal reservoir; represents the porosity of the rock of the thermal reservoir; t represents the average temperature of the thermal reservoir; t0 represents the surface temperature of the geothermal evaluation area; ρ w represents the density of geothermal water; c w represents the specific heat capacity of water.

[0074] Compared with calculating the geothermal resource quantity of the thermal reservoir through Steps 14 and 15, calculating the geothermal resource quantity of the thermal reservoir through the formula

[0075] is relatively simple. Thus, when the thickness and temperature of the thermal reservoir change little, a simple method can be adopted to calculate the geothermal resource quantity of the thermal reservoir to improve the calculation speed.

[0076] In practical applications, the maximum temperature value and the minimum temperature value in the thermal reservoir can be obtained. The maximum temperature value can be subtracted from the minimum temperature value to obtain the temperature change data. The maximum depth value and the minimum depth value of the thermal reservoir can be obtained. The maximum depth value of the thermal reservoir can be subtracted from the minimum depth value to obtain the thickness change data of the thermal reservoir. For example, the depth values of each position point on the top interface of the thermal reservoir can be obtained, and the minimum depth value can be selected from the depth values of each position point on the top interface. The depth values of each position point on the bottom interface of the thermal reservoir can be obtained, and the maximum depth value can be selected from the depth values of each position point on the bottom interface.

[0077] The temperature change threshold and the thickness change threshold can be empirical values respectively.

[0078] The technical solution of the embodiment of this specification can obtain the geological parameters of the geothermal evaluation area; can construct a three-dimensional geological model of the geothermal evaluation area according to the geological parameters; can select a heat reservoir of the geothermal evaluation area according to the three-dimensional geological model; can perform a numerical simulation of the temperature of the heat transfer process in the geothermal evaluation area to obtain the temperature field distribution of the geothermal evaluation area, and the temperature field distribution is used to represent the temperature of each position point in the geothermal evaluation area; can perform an integral calculation on the heat reservoir according to the temperature of each position point in the heat reservoir to obtain the geothermal resource amount of the heat reservoir. Thus, the embodiment of this specification can obtain a high-precision temperature field distribution of the geothermal evaluation area through the numerical simulation of the temperature of the heat transfer process in the geothermal evaluation area. The temperature field distribution reflects the distribution of the temperature in the geothermal evaluation area with respect to the spatial position. By using the temperature field distribution of the geothermal evaluation area and performing an integral calculation on the temperature of the heat reservoir, the geothermal resource amount of the heat reservoir can be obtained. In this way, the changes in the temperature and thickness of the heat reservoir are fully considered, avoiding the errors caused by using the average thickness of the heat reservoir to represent the overall thickness of the heat reservoir and using the average temperature of the heat reservoir to represent the overall temperature of the heat reservoir, and being able to improve the accuracy of the calculated geothermal resource amount, thereby being beneficial to more accurately evaluating the geothermal resource potential of the geothermal evaluation area.

[0079] Please refer to Figure 7 。 This embodiment of the specification also provides a geothermal resource amount calculation device, including:

[0080] An acquisition unit 71, configured to acquire the geological parameters of the geothermal evaluation area;

[0081] A construction unit 72, configured to construct a three-dimensional geological model of the geothermal evaluation area according to the geological parameters;

[0082] A selection unit 73, configured to select a heat reservoir of the geothermal evaluation area according to the three-dimensional geological model;

[0083] A simulation unit 74, configured to perform a numerical simulation of the temperature of the heat transfer process in the geothermal evaluation area to obtain the temperature field distribution of the geothermal evaluation area, and the temperature field distribution is used to represent the temperature of each position point in the geothermal evaluation area;

[0084] A calculation unit 75, configured to perform an integral calculation on the heat reservoir according to the temperature of each position point in the heat reservoir to obtain the geothermal resource amount of the heat reservoir.

[0085] This embodiment of the specification also provides a computer device, including a memory, a processor, and a computer program stored on the memory and executable on the processor. When the processor executes the computer program, the above-mentioned geothermal resource amount calculation method is implemented.

[0086] An embodiment of this specification also provides a computer-readable storage medium storing a computer program, which, when executed by a processor, implements the above-mentioned geothermal resource quantity calculation method.

[0087] An embodiment of this specification also provides a computer program product, which includes a computer program that, when executed by a processor, implements the above-mentioned geothermal resource quantity calculation method.

[0088] Those skilled in the art can understand that this specification can be provided as a method, a system, or a computer program product. Therefore, this specification can take the form of a complete hardware embodiment, a complete software embodiment, or an embodiment combining software and hardware. Moreover, this specification can take the form of a computer program product implemented on one or more computer-usable storage media (including but not limited to disk storage, CD-ROM, optical storage, etc.) containing computer-usable program code.

[0089] This specification is described with reference to the flowcharts and / or block diagrams of the methods, devices (systems), and computer program products according to the embodiments of this specification. It should be understood that each process and / or block in the flowcharts and / or block diagrams, and the combination of processes and / or blocks in the flowcharts and / or block diagrams, can be implemented by computer program instructions. The computer can be a personal computer, a laptop computer, a cellular phone, a camera phone, a smart phone, a personal digital assistant, a media player, a navigation device, an email device, a game console, a tablet computer, a wearable device, or any combination of these devices.

[0090] Each functional unit in the embodiments of this specification can be integrated into one processing unit, or each functional unit can exist physically alone, or two or more functional units can be integrated into one processing unit.

[0091] Those skilled in the art can understand that the descriptions of the embodiments in this specification each have their own emphases. For parts not detailed in a certain embodiment, reference can be made to the relevant descriptions of other embodiments. Additionally, it can be understood that after reading this specification document, those skilled in the art can, without creative labor, think of arbitrarily combining some or all of the embodiments listed in this specification, and these combinations are also within the scope of disclosure and protection of this specification.

[0092] Although this specification is depicted through embodiments, those of ordinary skill in the art know that the above embodiments are only used to help understand the core idea of this specification. Those skilled in the art can understand that this specification has many modifications and variations. It is hoped that the appended claims will cover these modifications and variations without departing from the spirit of this specification.

Claims

1. A method for calculating geothermal resource quantity, characterized in that, Including: Obtaining geological parameters of the geothermal evaluation area; Constructing a three-dimensional geological model of the geothermal evaluation area according to the geological parameters; Selecting a heat reservoir of the geothermal evaluation area according to the three-dimensional geological model; Performing a numerical simulation of the temperature in the heat transfer process of the geothermal evaluation area to obtain the temperature field distribution of the geothermal evaluation area, where the temperature field distribution is used to represent the temperature of each position point in the geothermal evaluation area; Performing an integral calculation on the heat reservoir according to the temperature of each position point in the heat reservoir to obtain the geothermal resource amount of the heat reservoir.

2. The method according to claim 1, wherein The geological parameters include altitude data, formation distribution data, fracture distribution data, Curie depth data, surface temperature data, Curie temperature data, and rock physical property data; The constructing of the three-dimensional geological model of the geothermal evaluation area includes: Constructing a three-dimensional geological model framework of the geothermal evaluation area with the Curie surface as the bottom boundary according to the altitude data, formation distribution data, fracture distribution data, and Curie depth data; assigning the surface temperature data to the top layer of the three-dimensional geological model framework, assigning the Curie temperature data to the bottom layer of the three-dimensional geological model framework, and assigning the rock physical property data to the corresponding formations in the three-dimensional geological model framework to obtain a three-dimensional geological model.

3. The method according to claim 1, wherein The selecting of the heat reservoir of the geothermal evaluation area according to the three-dimensional geological model includes: Selecting a formation that meets the set physical property conditions in the three-dimensional geological model as the heat reservoir of the geothermal evaluation area.

4. The method according to claim 3, characterized in that The performing of the numerical simulation of the temperature in the heat transfer process of the geothermal evaluation area to obtain the temperature field distribution of the geothermal evaluation area includes: Performing a meshing process on the three-dimensional geological model; Calculating the heat transfer process of the three-dimensional geological model from bottom to top according to the top and bottom temperature data of the three-dimensional geological model and the grids in the three-dimensional geological model to obtain the temperature field distribution of the geothermal evaluation area.

5. The method according to claim 1, characterized in that The performing of the integral calculation on the heat reservoir according to the temperature of each position point in the heat reservoir to obtain the geothermal resource amount of the heat reservoir includes: Extracting the temperature of each position point in the heat reservoir from the temperature field distribution; Calculating the geothermal resource amount of the heat reservoir using the triple integral formula according to the temperature of each position point in the heat reservoir.

6. The method according to claim 5, wherein The calculating of the geothermal resource amount of the heat reservoir using the triple integral formula according to the temperature data of each position point in the heat reservoir includes: Use the formula Q = ∫∫∫ Ω Rdv to calculate the geothermal resource quantity of the geothermal reservoir; where Q represents the geothermal resource quantity of the geothermal reservoir, Ω represents the volume of the geothermal reservoir, ρ r represents the rock density of the geothermal reservoir, c r represents the specific heat capacity of the rock of the geothermal reservoir, represents the porosity of the rock of the geothermal reservoir, t r represents the temperature at a position point within the geothermal reservoir, t0 represents the surface temperature, ρ w represents the density of geothermal water, c w represents the specific heat capacity of water, and dv represents the volume element of the geothermal reservoir area.

7. The method according to claim 1, characterized in that, The number of heat reservoirs is multiple; the method further includes: Calculating the geothermal resource potential of the geothermal evaluation area according to the geothermal resource amounts of multiple heat reservoirs.

8. A geothermal resource quantity calculation device, characterized in that, Including: An obtaining unit for obtaining geological parameters of the geothermal evaluation area; A constructing unit for constructing a three-dimensional geological model of the geothermal evaluation area according to the geological parameters; A selecting unit for selecting a heat reservoir of the geothermal evaluation area according to the three-dimensional geological model; A simulating unit for performing a numerical simulation of the temperature in the heat transfer process of the geothermal evaluation area to obtain the temperature field distribution of the geothermal evaluation area, where the temperature field distribution is used to represent the temperature of each position point in the geothermal evaluation area; A calculating unit for performing an integral calculation on the heat reservoir according to the temperature of each position point in the heat reservoir to obtain the geothermal resource amount of the heat reservoir.

9. A computer device, characterized in that, Including: A processor; A memory for storing processor-executable instructions; The processor executes the instructions to implement the method according to any one of claims 1-7.

10. A computer-readable storage medium, characterized in that, The computer-readable storage medium stores a computer program, and when the computer program is executed by a processor, the method according to any one of claims 1-7 is implemented.

Images