Method and device for calculating geothermal resource quantity

Abstract

The present specification relates to the technical field of geothermal, and particularly relates to a geothermal resource quantity calculation method and device. The geothermal resource quantity calculation method comprises the following steps: obtaining geological parameters of a geothermal evaluation area; constructing a three-dimensional geological model of the geothermal evaluation area according to the geological parameters; selecting a thermal reservoir of the geothermal evaluation area according to the three-dimensional geological model; performing temperature numerical simulation of a heat transfer process of the geothermal evaluation area to obtain a temperature field distribution of the geothermal evaluation area, wherein the temperature field distribution is used to represent temperatures of each position point in the geothermal evaluation area; and performing integral calculation on the thermal reservoir according to the temperatures of each position point in the thermal reservoir to obtain a geothermal resource quantity of the thermal reservoir. The embodiment of the present specification can improve the accuracy of geothermal resource quantity calculation.

Description

TECHNICAL FIELD

[0001] The present specification relates to the technical field of geothermal, and particularly relates to a geothermal resource quantity calculation method and device. BACKGROUND

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

[0003] However, the temperature and the thickness are different everywhere in the heat reservoir. For a region with large changes in the thickness and the temperature, the calculated geothermal resource quantity has a large error, and thus the reliability of the evaluation result of the geothermal resource potential in the geothermal evaluation area is low. SUMMARY

[0004] The embodiments of the present specification provide a geothermal resource quantity calculation method and device, which are used for the accuracy of the geothermal resource quantity calculation.

[0005] The embodiments of the present specification provide a geothermal resource quantity calculation method, which comprises the following steps.

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

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

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

[0009] Performing temperature numerical simulation of a heat transfer process of the geothermal evaluation area to obtain a temperature field distribution of the geothermal evaluation area, wherein 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 heat reservoir layer according to the temperature of each position point in the heat reservoir layer to obtain the geothermal resource quantity of the heat reservoir layer.

[0011] The embodiments of the present specification also provide a geothermal resource quantity calculation device, which comprises the following.

[0012] An obtaining unit, configured to obtain a geological parameter 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 parameter;

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

[0015] The simulation unit is used to perform temperature numerical simulation of the heat transfer process in the geothermal evaluation area to obtain the temperature field distribution of the geothermal evaluation area. The temperature field distribution is used to represent the temperature at each location point in the geothermal evaluation area.

[0016] The calculation unit is used to perform integral calculations on the geothermal reservoir based on the temperature at various locations within the reservoir, thereby obtaining the geothermal resource quantity of the reservoir.

[0017] This specification also provides a computer device, including a memory, a processor, and a computer program stored in the memory and executable on the processor. When the processor executes the computer program, it implements the above-described method for calculating geothermal resources.

[0018] This specification also provides a computer-readable storage medium storing a computer program that, when executed by a processor, implements the above-described method for calculating geothermal resources.

[0019] This specification also provides a computer program product, which includes a computer program that, when executed by a processor, implements the above-described geothermal resource calculation method.

[0020] The technical solutions of the embodiments in this specification can obtain geological parameters of a geothermal evaluation area; construct a three-dimensional geological model of the geothermal evaluation area based on the geological parameters; select geothermal reservoirs in the geothermal evaluation area based on the three-dimensional geological model; perform temperature numerical simulation of the heat transfer process in the geothermal evaluation area to obtain the temperature field distribution of the geothermal evaluation area, which represents the temperature at each location point in the geothermal evaluation area; and perform integral calculations on the geothermal reservoir based on the temperature at each location point in the geothermal reservoir to obtain the geothermal resource quantity of the geothermal reservoir. Therefore, the embodiments of this specification, through temperature numerical simulation of the heat transfer process in the geothermal evaluation area, can obtain a high-precision temperature field distribution of the geothermal evaluation area. The temperature field distribution reflects the spatial distribution of temperature in the geothermal evaluation area. By utilizing the temperature field distribution of the geothermal evaluation area and performing integral calculations on the temperature of the geothermal reservoir, the geothermal resource quantity of the geothermal reservoir can be obtained. This approach fully considers the variations in reservoir temperature and thickness, avoiding errors caused by using the average reservoir thickness to represent the overall reservoir thickness and the average reservoir temperature to represent the overall reservoir temperature. This improves the accuracy of the calculated geothermal resource quantity, thus facilitating a more accurate assessment of the geothermal resource potential of the geothermal evaluation area. Attached Figure Description

[0021] To more clearly illustrate the technical solutions in the embodiments or prior art of this specification, the drawings used in the description of the embodiments or prior art will be briefly introduced below. The drawings described below are only some embodiments recorded in this specification. For those skilled in the art, other drawings can be obtained based on these drawings without creative effort.

[0022] Figure 1 This is a flowchart illustrating the geothermal resource calculation method in the embodiments of this specification;

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

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

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

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

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

[0028] Figure 7 This is a schematic diagram of the geothermal resource calculation device in the embodiments of this specification. Detailed Implementation

[0029] The technical solutions in the embodiments of this specification will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some embodiments of this specification, and not all embodiments. The specific embodiments described herein are only used to explain this disclosure, and not to limit this disclosure. All other embodiments obtained by those skilled in the art based on the described embodiments of this disclosure are within the scope of protection of this 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 such actual relationship or order between these entities or operations.

[0030] Please see Figure 1 This specification provides a method for calculating geothermal resources. The method can be applied to computer equipment and may specifically 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 region used to evaluate the potential of geothermal resources. The geological parameters may include geological structural data, temperature data, and rock property data. The geological structural data represents the geological structure of the geothermal evaluation area. This data may include high-altitude data, stratigraphic distribution data, and fault distribution data. The high-altitude data represents the topographic relief and altitude of the geothermal evaluation area. The stratigraphic distribution data represents the stratigraphic structure, lithological distribution, and stratigraphic thickness of the geothermal evaluation area. The fault distribution data represents the location, strike, dip angle, and other structural features of fault zones. The temperature data includes Curie depth data, Curie temperature data, and surface temperature data. The Curie isotherm is located in the Earth's crust. The Curie depth data represents the depth of the Curie isotherm. The Curie temperature data represents the temperature of the Curie isotherm. The surface temperature data represents the surface temperature. The rock property data represents the thermophysical properties of the rocks in the geothermal evaluation area, such as thermal conductivity, specific heat capacity, density, porosity, permeability, etc.

[0033] In some embodiments, an area requiring geothermal resource evaluation can be selected as the geothermal evaluation zone; geological parameters of the geothermal evaluation zone can be obtained. For example, an area within a defined latitude and longitude range that requires geothermal resource evaluation can be selected as the geothermal evaluation zone. Figure 2 This is a topographic map of the geothermal assessment area. The area includes multiple volcanoes and numerous hot springs and geothermal wells. For example, it contains hot spring clusters 1, 2, 3, and 4, as well as geothermal wells 1, 2, 3, and 4, indicating good geothermal resources. However, due to complex geological conditions and a coarse current geothermal field, the exact amount of geothermal resources in the assessment area is unclear.

[0034] Geological parameters of the geothermal evaluation area can be obtained through satellite remote sensing, seismic exploration, meteorological station data, and field measurements. For example, elevation data can be obtained through satellite remote sensing or field measurements. Surface temperature data can be obtained through meteorological station data, satellite remote sensing, or field measurements. Stratigraphic distribution data and fault distribution data can be obtained through seismic exploration. Curie depth data and Curie temperature data can be obtained through geothermal gradient measurements or geothermal flow measurements. Rock physical property data can be obtained through laboratory testing.

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

[0036] In some embodiments, a three-dimensional geological model framework with the Curie line as the base can be constructed for the geothermal evaluation area based on high-altitude data, stratigraphic distribution data, fault distribution data, and Curie depth data.

[0037] The three-dimensional geological model framework is used to represent the surface, stratigraphic interfaces, and curie zone of the geothermal evaluation area.

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

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

[0040] In some embodiments, surface temperature data, Curie temperature data, and rock property data can be integrated into the three-dimensional geological model framework to obtain a three-dimensional geological model. For example, surface temperature data can be assigned to the top layer of the three-dimensional geological model framework; Curie temperature data can be assigned to the bottom layer of the three-dimensional geological model framework; and rock property data can be assigned to the corresponding strata 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 temperature, and rock property data of each stratum in the geothermal evaluation area.

[0041] Step 13: Select the geothermal reservoir in the geothermal evaluation area based on the three-dimensional geological model.

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

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

[0044] The set physical property conditions may include a first physical property data interval. The physical property data within the first physical property data interval 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 interval is greater than or equal to the first porosity threshold, the permeability of the physical property data within the first physical property data interval 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 interval is greater than or equal to the first thermal conductivity threshold. Strata in the three-dimensional geological model whose corresponding physical property data are located within the first physical property data interval can be selected as the geothermal reservoir of the geothermal evaluation area. For example, strata in the three-dimensional geological model whose porosity is greater than or equal to the first porosity threshold, permeability is greater than or equal to the first permeability threshold, and thermal conductivity is greater than or equal to the first thermal conductivity threshold can be selected as the geothermal reservoir of the geothermal evaluation area.

[0045] Furthermore, the set physical property conditions may also 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. Additionally, 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 strata whose physical property data fall within the first physical property data interval can be selected from the three-dimensional geological model as one or more candidate geothermal reservoirs. For each candidate geothermal reservoir, the strata overlying it can be obtained as the overlying strata. From the one or more candidate geothermal reservoirs, the candidate geothermal reservoir whose corresponding overlying strata's physical property data falls within the second physical property data interval can be selected as the geothermal reservoir for the geothermal evaluation area. For example, candidate geothermal reservoirs with porosity less than or equal to a second porosity threshold, permeability less than or equal to a second permeability threshold, and thermal conductivity less than or equal to a second thermal conductivity threshold can be selected as the geothermal reservoir for the geothermal evaluation area.

[0047] Therefore, the geothermal reservoir in the geothermal assessment area should be a stratum with good physical properties, such as high porosity, permeability, and thermal conductivity; and it should be overlain by strata with poor physical properties, such as strata with low porosity, low permeability, and low thermal conductivity. This way, the selected reservoir is more representative of the geothermal assessment area, and its selection can more accurately characterize the geothermal resource potential of the area.

[0048] Please see Figure 3 The Permian and Triassic strata can be selected as the geothermal reservoirs for the geothermal evaluation area.

[0049] Step 14: Perform a temperature numerical simulation of the heat transfer process in 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. Based on the top and bottom temperature data of the three-dimensional geological model, and following the mesh within the model, the heat transfer process of the three-dimensional geological model can be calculated from bottom to top to obtain the temperature field distribution of the geothermal assessment area. The temperature field distribution can provide a data basis for calculating geothermal resources. The temperature field distribution represents the spatial distribution of temperature in the geothermal assessment area. For example, the temperature field distribution can represent the surface temperature, Curie temperature, and temperature data of various strata in the geothermal assessment area. The temperature field distribution can include the current temperature field distribution. The current temperature field distribution represents the current spatial distribution of temperature in the geothermal assessment area.

[0051] Figure 4 The three-dimensional geological model after gridding is shown. Figure 5 The temperature field distribution of the geothermal evaluation area is shown.

[0052] In some embodiments, the grid can be viewed as the granularity for numerical temperature simulation during heat transfer. The grid size affects the accuracy of the temperature field distribution. A smaller grid represents a smaller spatial area, allowing for more precise capture of temperature variations across different spatial locations and improving the overall accuracy of the temperature field distribution. However, this results in a higher computational load and longer simulation time. Conversely, a larger grid represents a larger spatial area, leading to lower overall accuracy of the temperature field distribution. However, this results in a lower computational load and shorter simulation time.

[0053] The grid size can be determined based on the complexity of the geothermal evaluation area and the required accuracy of the temperature numerical simulation.

[0054] A first set of grid sizes can be provided. This first set of grid sizes may include multiple grid sizes. Each grid size corresponds to a specific computational precision. The computational precision corresponding to each grid size can be selected from multiple pre-defined precision levels. For example, the multiple precision levels may include high, medium, and low precision levels. A target computational precision input by the user can be received. The target computational precision can be selected from the multiple pre-defined precision levels. A corresponding first target grid size can be selected from the first set of grid sizes based on the target computational precision. For example, the first target grid size could be 20 meters.

[0055] The three-dimensional geological model can be meshed according to the first target mesh size. Therefore, in the meshed three-dimensional geological model, the mesh size is always the first target mesh size. Based on the top and bottom temperature data of the three-dimensional geological model, and following the mesh within the model, the heat transfer process of the three-dimensional geological model can be calculated from bottom to top to obtain the temperature field distribution of the geothermal evaluation area.

[0056] Optionally, the geological conditions of the geothermal evaluation area are relatively complex. If a uniform first target grid size is used to mesh the three-dimensional geological model, it is impossible to balance the accuracy of the temperature field distribution with the computational load. 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. Based on the top and bottom temperature data of the three-dimensional geological model, and according to the optimized grid in the three-dimensional geological model, the heat transfer process of the three-dimensional geological model can be calculated from bottom to top to obtain the temperature field distribution of the geothermal evaluation area.

[0057] The target area is a region with complex geological conditions within the geothermal evaluation zone. For example, the target area may include areas where fault zones are located. The geological complexity of the target area can be calculated; the first target grid size of the target area can be optimized and adjusted based on the geological complexity. 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 range. The complexity range of the geological complexity of the target area can be obtained; based on the obtained complexity range, 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 meshed using the second target grid size. The magnitude of the geological complexity of the target area is inversely correlated with the size 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. Therefore, within the geothermal evaluation zone, a smaller second target grid size can be used for meshing the target area, while a larger first target grid size can be used for meshing other areas besides the target area. This allows for a better balance between the overall accuracy of the temperature field distribution and the computational load.

[0058] It can acquire characteristic data of a target area and calculate the geological complexity of that area using a predefined formula based on this data. For example, the target area may include the region containing a fault zone. The characteristic data of a fault zone may include the number of branches and intersections, curvature, and gradients in physical property changes. The geological complexity of the fault zone can be calculated using a predefined formula based on these factors.

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

[0060] Heat transfer is a physical phenomenon in physics, referring to the transfer of heat energy caused by a temperature difference. The change in the internal energy of an object can be measured in heat during heat transfer. Heat transfer can include three forms: conduction, radiation, and convection. As long as a temperature difference exists within or between objects, heat energy can be transferred from a high-temperature to a low-temperature area through one or more of these three methods. Within terrestrial areas, conduction and convection are the primary forms of heat transfer. Therefore, based on the top and bottom temperature data of a three-dimensional geological model and the rock 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 grid in the model, thus obtaining the current temperature field distribution of the geothermal assessment area.

[0061] Numerical simulation, also known as computer simulation, relies on electronic computers and combines the concepts of finite element or finite volume to achieve the purpose of studying engineering problems, physical problems, and even various problems in nature through numerical calculation or graphical display methods.

[0062] The geothermal assessment area possesses abundant basic geological data (such as fault distribution characteristics, stratigraphic distribution characteristics, rock property data, Curie depth and temperature, surface temperature, and elevation data). While this data cannot directly reflect the temperature field distribution of the geothermal assessment area, by substituting it into the heat conduction equation, the temperature distribution and geothermal heat flow distribution of the geothermal assessment area can be calculated, thus obtaining the temperature field distribution. Considering the enormous computational burden, a three-dimensional geological model of the geothermal assessment area can be established based on this data, and the temperature field distribution can be calculated using this model.

[0063] For example, the heat conduction equation can include ρ represents density; c represents specific heat capacity; Let T be the partial derivative of temperature T with respect to time t, representing the rate of change of temperature over time; The gradient operator is represented by λ; thermal conductivity is represented by λ. Q is the thermal conductivity term, representing the propagation of heat in space through thermal conduction, 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 and Newton's iteration method) can be used to solve the heat conduction equation, progressively updating the temperature values ​​of each grid in the three-dimensional geological model. 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 ending 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 stratum where the grid is located can be obtained. Based on the heat source term, thermal conductivity, specific heat capacity, and density, the heat conduction equation can be solved numerically to obtain the temperature value of that grid. By progressively calculating the temperature values ​​of each grid using the heat conduction equation in a bottom-up order, the temperature field distribution of the geothermal evaluation area can be obtained. For example, the temperature field distribution of the geothermal evaluation area may include: surface temperature 15℃, temperature at 100m depth 80℃, temperature at 300m depth 150℃, temperature at 600m depth 300℃, and Curie temperature 450℃.

[0065] Step 15: Based on the temperature at each location point in the geothermal reservoir, perform an integral calculation on the geothermal reservoir to obtain the geothermal resource quantity of the geothermal reservoir.

[0066] In some embodiments, the temperature field distribution can be used to represent the temperature at various locations within the geothermal assessment area. These locations may, for example, include grids in a three-dimensional geological model. The geothermal resource quantity of the reservoir can be obtained by integrating the temperatures at each location within the reservoir. This fully considers the variations in reservoir temperature and thickness, avoiding errors caused by using the average reservoir thickness to represent the overall reservoir thickness and the average reservoir temperature to represent the overall reservoir temperature. This improves the accuracy of geothermal resource quantity estimation and facilitates accurate assessment of the geothermal resource potential of the geothermal assessment area.

[0067] In some embodiments, temperature data at various points within the geothermal reservoir can be extracted from the temperature field distribution; the geothermal resource quantity of the reservoir can be calculated using a triple integral formula based on the temperature data at various points within the reservoir. Thus, by utilizing the temperature field distribution of the geothermal evaluation area and integrating the temperatures at various points within the reservoir, the geothermal resource quantity of the reservoir can be obtained.

[0068] For example, the formula Q=∫∫∫ can be used. Ω Rdv calculates the geothermal resource volume of a geothermal reservoir; where Q represents the geothermal resource volume of the reservoir, and Ω represents the volume of the reservoir. ρ r c represents the rock density of a thermal reservoir. r This indicates the specific heat capacity of the rock in a thermal reservoir. t represents the rock porosity of a thermal reservoir. r ρ represents the temperature at a point within the thermal reservoir, t0 represents the surface temperature, and ρ represents the temperature at that point. w c represents the density of geothermal water. w ρ represents the specific heat capacity of water, and dv represents the volumetric element of the thermal reservoir region. r c r , t r , t0, ρ w c w The values ​​of , etc., can be preset. dv = dxdydz. x, y, and z represent the coordinate data of a point within the thermal reservoir in a three-dimensional Cartesian coordinate system. The three-dimensional Cartesian coordinate system is a coordinate system established specifically for the thermal reservoir.

[0069] For example, in step 14, the three-dimensional geological model can be meshed. Location points within the reservoir can then include the mesh within the reservoir. Temperature values ​​of each mesh within the thermal reservoir can be extracted from the temperature field distribution; the formula Q=∫∫∫ can be used based on the temperature values ​​of each mesh within the reservoir. Ω Rdv calculates the geothermal resources of the reservoir. t r This represents the temperature value of the grid within the geothermal reservoir. Specifically, for example, the geothermal resource quantity of the reservoir could be 5.42 × 10⁻⁶.18 J.

[0070] In some embodiments, the geothermal resource potential of the geothermal assessment area can be calculated based on the geothermal resource quantity of the reservoir. For example, the geothermal resource quantity of the reservoir can be directly used to characterize the geothermal resource potential of the assessment area. Alternatively, the geothermal resource potential of the assessment area can be calculated using the formula E = Q × η, where Q represents the geothermal resource quantity of the reservoir and η represents the recoverability coefficient. The recoverability coefficient η can be an empirical value, or it can be calculated based on the physical property data of the rocks in the reservoir. For example, if the permeability of the rocks in the reservoir is high, the recoverability coefficient η can be taken as 0.2.

[0071] For example, the geothermal resources of the thermal reservoir could be 5.42 × 10⁻⁶. 18 J, the geothermal power generation potential is approximately 2.11 × 10⁻⁶. 11 kWh. This indicates that the geothermal evaluation area has good geothermal power generation potential. The geothermal reservoir is an area with advantageous geothermal resource potential. Figure 6 The location of the geothermal resource potential advantage area in the geothermal evaluation area is shown, and it is speculated that it is related to the underlying magma chamber system.

[0072] In some embodiments, multiple geothermal reservoirs can be selected through step 13. The geothermal resource quantities of the multiple reservoirs can be obtained through step 15. The geothermal resource potential of the geothermal assessment area can then be calculated based on the geothermal resource quantities of the multiple reservoirs. For example, the geothermal resource quantities of the multiple reservoirs can be added together, and the geothermal resource potential of the geothermal assessment area can be determined based on the sum. Alternatively, the geothermal resource potential of the geothermal assessment area can be directly characterized using the sum. Or, the geothermal resource potential of the geothermal assessment area can be calculated using the formula E = Q × η based on the sum. Here, Q represents the sum of the geothermal resource quantities of the multiple reservoirs, and η is the recoverability coefficient.

[0073] In some embodiments, temperature change data and thickness change data of the geothermal reservoir can be obtained. If the temperature change data is greater than or equal to a temperature change threshold, and / or the thickness change data is greater than or equal to a thickness change threshold, it indicates that the thickness and / or temperature change of the geothermal reservoir is relatively large. Therefore, step 14 can be used to perform a temperature numerical simulation of the heat transfer process in the geothermal evaluation area to obtain the temperature field distribution of the geothermal evaluation area; step 15 can be used to perform an integral calculation on the geothermal reservoir based on the temperature at each location point in the geothermal reservoir to obtain the geothermal resource quantity of the geothermal reservoir. If the temperature change data is less than a temperature change threshold, and the thickness change data is less than a thickness change threshold, it indicates that the thickness and temperature change of the geothermal reservoir are relatively small, and thus can be obtained using the formula... Calculate the geothermal resource quantity of the reservoir. Q represents the geothermal heat of the reservoir; A represents the area of ​​the reservoir; Z represents the average thickness of the reservoir; ρ r c represents the rock density of the thermal reservoir. r This indicates the specific heat capacity of the rock in a thermal reservoir; The ρ represents the rock porosity of the geothermal reservoir; t represents the average temperature of the geothermal reservoir; t0 represents the surface temperature of the geothermal evaluation area; ρ w c represents the density of geothermal water. w This indicates the specific heat capacity of water.

[0074] Compared to calculating the geothermal resources of the reservoir through steps 14 and 15, the formula...

[0075] Calculating the geothermal resources of a geothermal reservoir is relatively simple. Therefore, when the thickness and temperature variations of the reservoir are small, a simplified method can be used to calculate the geothermal resources, improving calculation speed.

[0076] In practical applications, the highest and lowest temperature values ​​within the thermal reservoir can be obtained. Subtracting the highest and lowest temperature values ​​yields the temperature variation data. Similarly, the maximum and minimum depth values ​​of the thermal reservoir can be obtained. Subtracting the maximum and minimum depth values ​​yields the thickness variation data. For example, the depth values ​​at various points on the top interface of the thermal reservoir can be obtained, and the minimum depth value can be selected from these values. Likewise, the depth values ​​at various points on the bottom interface of the thermal reservoir can be obtained, and the maximum depth value can be selected from these values.

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

[0078] The technical solutions of the embodiments in this specification can obtain geological parameters of a geothermal evaluation area; construct a three-dimensional geological model of the geothermal evaluation area based on the geological parameters; select geothermal reservoirs in the geothermal evaluation area based on the three-dimensional geological model; perform temperature numerical simulation of the heat transfer process in the geothermal evaluation area to obtain the temperature field distribution of the geothermal evaluation area, which represents the temperature at each location point in the geothermal evaluation area; and perform integral calculations on the geothermal reservoir based on the temperature at each location point in the geothermal reservoir to obtain the geothermal resource quantity of the geothermal reservoir. Therefore, the embodiments of this specification, through temperature numerical simulation of the heat transfer process in the geothermal evaluation area, can obtain a high-precision temperature field distribution of the geothermal evaluation area. The temperature field distribution reflects the spatial distribution of temperature in the geothermal evaluation area. By utilizing the temperature field distribution of the geothermal evaluation area and performing integral calculations on the temperature of the geothermal reservoir, the geothermal resource quantity of the geothermal reservoir can be obtained. This approach fully considers the variations in reservoir temperature and thickness, avoiding errors caused by using the average reservoir thickness to represent the overall reservoir thickness and the average reservoir temperature to represent the overall reservoir temperature. This improves the accuracy of the calculated geothermal resource quantity, thus facilitating a more accurate assessment of the geothermal resource potential of the geothermal evaluation area.

[0079] Please see Figure 7 This specification also provides a geothermal resource calculation device, comprising:

[0080] Acquisition unit 71 is used to acquire geological parameters of the geothermal evaluation area;

[0081] Construction unit 72 is used to construct a three-dimensional geological model of the geothermal evaluation area based on the geological parameters;

[0082] Selecting unit 73 is used to select the geothermal reservoir in the geothermal evaluation area based on the three-dimensional geological model;

[0083] Simulation unit 74 is used to perform temperature numerical simulation of the heat transfer process in the geothermal evaluation area to obtain the temperature field distribution of the geothermal evaluation area. The temperature field distribution is used to represent the temperature at each location point in the geothermal evaluation area.

[0084] The calculation unit 75 is used to perform integral calculations on the geothermal reservoir based on the temperature at each location point in the geothermal reservoir to obtain the geothermal resource quantity of the geothermal reservoir.

[0085] This specification also provides a computer device, including a memory, a processor, and a computer program stored in the memory and executable on the processor. When the processor executes the computer program, it implements the above-described method for calculating geothermal resources.

[0086] This specification also provides a computer-readable storage medium storing a computer program that, when executed by a processor, implements the above-described method for calculating geothermal resources.

[0087] This specification also provides a computer program product, which includes a computer program that, when executed by a processor, implements the above-described geothermal resource calculation method.

[0088] Those skilled in the art will understand that this specification can be provided as a method, system, or computer program product. Therefore, this specification may take the form of a completely hardware embodiment, a completely software embodiment, or an embodiment combining software and hardware. Furthermore, this specification may take the form of a computer program product embodied 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 flowchart illustrations and / or block diagrams of methods, apparatus (systems), and computer program products according to embodiments thereof. It should be understood that each block of the flowchart illustrations and / or block diagrams, and combinations of blocks in the flowchart illustrations and / or block diagrams, can be implemented by computer program instructions. The computer may be a personal computer, laptop computer, cellular phone, camera phone, smartphone, personal digital assistant, media player, navigation device, email device, game console, tablet computer, wearable device, or any combination of these devices.

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

[0091] Those skilled in the art will understand that the descriptions of the various embodiments in this specification have different focuses, and parts not described in detail in a certain embodiment can be referred to in the relevant descriptions of other embodiments. Furthermore, it is understood that those skilled in the art, after reading this specification, can conceive of any combination of some or all of the embodiments listed in this specification without creative effort, and such combinations are also within the scope of disclosure and protection of this specification.

[0092] Although this specification has been described through embodiments, those skilled in the art will understand that the above embodiments are merely illustrative of the core ideas of this specification. Those skilled in the art will appreciate that many variations and modifications are possible with this specification. It is intended that the appended claims encompass these variations and modifications without departing from the spirit of this specification.

Claims

1. A method for calculating geothermal resources, characterized in that, include: Obtain geological parameters of the geothermal evaluation area; Based on the geological parameters, a three-dimensional geological model of the geothermal evaluation area is constructed; The geothermal reservoir in the geothermal evaluation area was selected based on the three-dimensional geological model. A temperature numerical simulation of the heat transfer process is performed on the geothermal evaluation area to obtain the temperature field distribution of the geothermal evaluation area. The temperature field distribution is used to represent the temperature at each location point in the geothermal evaluation area. Based on the temperature at various points in the geothermal reservoir, the geothermal resource quantity of the geothermal reservoir is obtained by integral calculation. The step of integrating the geothermal resources of the reservoir based on the temperature at various points within the reservoir to obtain the geothermal resource quantity includes: extracting the temperature at various points within the reservoir from the temperature field distribution; calculating the geothermal resource quantity of the reservoir using a triple integral formula based on the temperature at various points within the reservoir; and calculating the geothermal resource quantity of the reservoir using a triple integral formula based on the temperature data at various points within the reservoir, including: using the formula... Calculate the geothermal resource quantity of the thermal reservoir; where Q represents the geothermal resource quantity of the thermal reservoir, and Ω represents the volume of the thermal reservoir. , ρ r c represents the rock density of a thermal reservoir. r φ represents the specific heat capacity of the rock in the thermal reservoir, φ represents the porosity of the rock in the thermal reservoir, and t represents the specific heat capacity of the rock in the thermal reservoir. r ρ represents the temperature at a point within the thermal reservoir, t0 represents the surface temperature, and ρ represents the temperature at that point. w c represents the density of geothermal water. w dv represents the specific heat capacity of water, and dv represents the volumetric element of the thermal reservoir region.

2. The method according to claim 1, characterized in that, The geological parameters include high-altitude data, stratigraphic distribution data, fault distribution data, Curie depth data, surface temperature data, Curie temperature data, and rock physical property data. The construction of the three-dimensional geological model of the geothermal evaluation area includes: Based on high-altitude data, stratigraphic distribution data, fault distribution data, and Curie depth data, a three-dimensional geological model framework with Curie as the bottom boundary is constructed for the geothermal evaluation area. Surface temperature data is assigned to the top layer of the three-dimensional geological model framework, Curie temperature data is assigned to the bottom layer of the three-dimensional geological model framework, and rock property data is assigned to the corresponding strata in the three-dimensional geological model framework to obtain the three-dimensional geological model.

3. The method according to claim 1, characterized in that, The selection of geothermal reservoirs in the geothermal evaluation area based on the three-dimensional geological model includes: In the three-dimensional geological model, strata that meet the set physical property conditions are selected as the geothermal reservoirs of the geothermal evaluation area.

4. The method according to claim 1, characterized in that, The temperature numerical simulation of the heat transfer process in the geothermal evaluation area, to obtain the temperature field distribution of the geothermal evaluation area, includes: The three-dimensional geological model is then meshed. Based on the top and bottom temperature data of the three-dimensional geological model, the heat transfer process of the three-dimensional geological model is calculated from bottom to top according to the grid 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 number of thermal reservoirs is multiple; the method further includes: The geothermal resource potential of the geothermal evaluation area is calculated based on the geothermal resource quantities of multiple thermal reservoirs.

6. A geothermal resource calculation device, characterized in that, include: The acquisition unit is used to acquire geological parameters of the geothermal evaluation area; A construction unit is used to construct a three-dimensional geological model of the geothermal evaluation area based on the geological parameters. Selecting a unit is used to select the geothermal reservoir in the geothermal evaluation area based on the three-dimensional geological model; The simulation unit is used to perform temperature numerical simulation of the heat transfer process in the geothermal evaluation area to obtain the temperature field distribution of the geothermal evaluation area. The temperature field distribution is used to represent the temperature at each location point in the geothermal evaluation area. The calculation unit is used to perform integral calculations on the geothermal reservoir based on the temperature at each location point in the geothermal reservoir to obtain the geothermal resource quantity of the geothermal reservoir; The step of integrating the geothermal resources of the reservoir based on the temperature at various points within the reservoir to obtain the geothermal resource quantity includes: extracting the temperature at various points within the reservoir from the temperature field distribution; calculating the geothermal resource quantity of the reservoir using a triple integral formula based on the temperature at various points within the reservoir; and calculating the geothermal resource quantity of the reservoir using a triple integral formula based on the temperature data at various points within the reservoir, including: using the formula... Calculate the geothermal resource quantity of the thermal reservoir; where Q represents the geothermal resource quantity of the thermal reservoir, and Ω represents the volume of the thermal reservoir. , ρ r c represents the rock density of a thermal reservoir. r φ represents the specific heat capacity of the rock in the thermal reservoir, φ represents the porosity of the rock in the thermal reservoir, and t represents the specific heat capacity of the rock in the thermal reservoir. r ρ represents the temperature at a point within the thermal reservoir, t0 represents the surface temperature, and ρ represents the temperature at that point. w c represents the density of geothermal water. w dv represents the specific heat capacity of water, and dv represents the volumetric element of the thermal reservoir region.

7. A computer device, characterized in that, include: processor; Memory used to store processor-executable instructions; The processor executes the instructions to implement the method as described in any one of claims 1-5.

8. A computer-readable storage medium, characterized in that, The computer-readable storage medium stores a computer program that, when executed by a processor, implements the method of any one of claims 1-5.

Images