Geothermal favorable area prediction method, device and equipment based on numerical simulation and medium

Abstract

The invention discloses a geothermal favorable area prediction method and device based on numerical simulation, equipment and a medium. The method comprises the following steps: determining a substrate earth heat flow in a research area; establishing a two-dimensional geologic model, and counting the thermal conductivity value of each lithology of each stratum above the bedrock in the research area; performing thermal conductivity assignment on the two-dimensional geologic model, performing thermal conduction numerical simulation, and storing a temperature curve with a fixed depth; according to the temperature curve of the fixed depth, drawing a ground temperature contour plane graph of the fixed depth in the research area; analyzing geothermal conditions in the research area according to the basic geological data and the geochemical data; and determining the correlation between the ground temperature isoline plane graph and the geothermal conditions, and predicting the geothermal favorable area of the research area. According to the method for predicting the geothermal resource favorable area through numerical simulation, the temperature information in the research area can be accurately obtained, the precision is high, the efficiency is good, the method is more economical and practical, and convenience is provided for detection of the ground surface geothermal abnormal area.

Description

Technical Field

[0001] The present invention relates to the field of calibration of favorable areas for geothermal anomalies, and more specifically, to a method, device, equipment and medium for predicting favorable areas for geothermal energy based on numerical simulation. Background Art

[0002] With the prosperity of the economy, the external dependence on fossil energy is getting higher and higher, which has constituted a potential strategic risk in the energy field. Geothermal resources are attracting wide attention due to their advantages of being clean and renewable.

[0003] Currently, there are few means and methods for predicting favorable areas for geothermal resources, which are time-consuming, costly, and have a high risk. The underground temperature data in the study area is scarce.

[0004] Currently, there is a need to develop a method, device, equipment and medium for predicting favorable areas for geothermal energy based on numerical simulation.

[0005] The information disclosed in the background art section of the present invention is only intended to deepen the understanding of the general background art of the present invention, and should not be regarded as an admission or any form of implication that this information constitutes the prior art known to those skilled in the art. Summary of the Invention

[0006] The present invention provides a method, device, equipment and medium for predicting favorable areas for geothermal energy based on numerical simulation. The method for predicting favorable areas for geothermal resources by using numerical simulation can obtain the temperature information in the study area more accurately, with high precision and good efficiency, and is more economical and practical, providing convenience for the detection of surface geothermal anomaly areas.

[0007] In a first aspect, an embodiment of the present disclosure provides a method for predicting favorable areas for geothermal energy based on numerical simulation, including:

[0008] Determine the basal terrestrial heat flow in the study area;

[0009] Establish a two-dimensional geological model and statistically analyze the thermal conductivity values of each lithology of each stratum above the bedrock in the study area;

[0010] After assigning thermal conductivity values to the two-dimensional geological model, perform numerical simulation of heat conduction and save the temperature curve at a fixed depth;

[0011] According to the temperature curve at the fixed depth, draw a plan of geothermal isotherms at the fixed depth in the study area;

[0012] Analyze the geothermal conditions in the study area based on basic geological data and geochemical data;

[0013] Determine the correlation between the plan of geothermal isotherms and the geothermal conditions, and predict the favorable areas for geothermal energy in the study area.

[0014] Preferably, the basal terrestrial heat flow is calculated by formula (1):

[0015] Q = -K·dT / dZ (1)

[0016] Wherein, Q is the basal terrestrial heat flow, K is the thermal conductivity of the basal rock, and dT / dZ is the geothermal gradient of the basement layer in the study area.

[0017] Preferably, if the well is shallow and the temperature data of the basement layer is missing, the basal terrestrial heat flow is calculated by formula (2):

[0018] Q = Qb - A·Z (2)

[0019] Wherein, Q is the basal terrestrial heat flow, Qb is the surface heat flow, Z is the thickness of the sedimentary cover, and A is the heat generation rate of the sedimentary cover.

[0020] Preferably, a two-dimensional geological model is established, and the thermal conductivity values of each lithology of each stratum above the bedrock in the study area are statistically included:[[]]

[0021] Through seismic interpretation, the seismic profile and the thickness distribution of each stratum in the study area are obtained, and a two-dimensional geological model is established through AutoCAD. The two-dimensional geological model is imported into COMSOL software, and each stratum is meshed and thermal conductivity data is assigned.

[0022] Preferably, after the thermal conductivity is assigned to the two-dimensional geological model, a numerical simulation of heat conduction is performed to save the temperature curve at a fixed depth:

[0023] The simulation is performed using the heat conduction simulation module in COMSOL software, with the density set to 2500 kg / m 3 , the specific heat is 1000 J / (kg·°C), the surface temperature is the average surface temperature of the study area, and the initial geothermal gradient is 3 °C / km. A numerical simulation of heat conduction is performed to save the temperature curve at a fixed depth.

[0024] Preferably, drawing the geothermal isotherm plan at a fixed depth in the study area includes:

[0025] According to the temperature curve at the fixed depth and combined with the temperature of the actual well drilling, the geothermal isotherm plan at the fixed depth in the study area is drawn using the Kriging interpolation method on Surfer software.

[0026] Preferably, the geothermal conditions include the tectonic distribution characteristics, lithological characteristics, heat source location, and fluid movement direction.

[0027] In a second aspect, the embodiments of the present disclosure further provide a geothermal favorable area prediction device based on numerical simulation, including:

[0028] A calculation module that determines the basal terrestrial heat flow in the study area;

[0029] A statistical module that establishes a two-dimensional geological model and statistically calculates the thermal conductivity values of each lithology in each stratum above the bedrock in the study area;

[0030] A simulation module that performs a thermal conductivity assignment on the two-dimensional geological model, conducts a numerical simulation of heat conduction, and saves the temperature curve at a fixed depth;

[0031] A plotting module that plots a planar isogeothermal line map at a fixed depth in the study area based on the temperature curve at the fixed depth;

[0032] An analysis module that analyzes the geothermal conditions in the study area based on basic geological data and geochemical data;

[0033] A prediction module that determines the correlation between the planar isogeothermal line map and the geothermal conditions, and predicts the favorable geothermal areas in the study area.

[0034] Preferably, the basal terrestrial heat flow is calculated by formula (1):

[0035] Q = -K·dT / dZ (1)

[0036] Where Q is the basal terrestrial heat flow, K is the thermal conductivity of the basal rock, and dT / dZ is the geothermal gradient of the basement layer in the study area.

[0037] Preferably, if the well is shallow and the temperature data of the basement layer is missing, the basal terrestrial heat flow is calculated by formula (2):

[0038] Q = Qb - A·Z (2)

[0039] Where Q is the basal terrestrial heat flow, Qb is the surface heat flow, Z is the thickness of the sedimentary cover, and A is the heat production rate of the sedimentary cover.

[0040] Preferably, establishing a two-dimensional geological model and statistically calculating the thermal conductivity values of each lithology in each stratum above the bedrock in the study area includes:

[0041] Through seismic interpretation, obtaining the seismic profile and the thickness distribution of each stratum in the study area, establishing a two-dimensional geological model through AutoCAD, importing the two-dimensional geological model into COMSOL software, and performing meshing and thermal conductivity data assignment on each stratum.

[0042] Preferably, after performing a thermal conductivity assignment on the two-dimensional geological model, conducting a numerical simulation of heat conduction, and saving the temperature curve at a fixed depth:

[0043] Using the heat conduction simulation module in COMSOL software for simulation, setting the density to 2500 kg / m 3, with a specific heat capacity of 1000 J / (kg·℃), the surface temperature being the average surface temperature of the study area, and the initial geothermal gradient being 3℃ / km, a numerical simulation of heat conduction is carried out, and the temperature curve at a fixed depth is saved.

[0044] Preferably, drawing the isogeothermal contour plan at a fixed depth in the study area includes:

[0045] According to the temperature curve at the fixed depth and combined with the temperature of the actual drilling, using the Kriging interpolation method on the Surfer software, draw the isogeothermal contour plan at a fixed depth in the study area.

[0046] Preferably, the geothermal conditions include tectonic distribution characteristics, lithological characteristics, heat source location, and fluid movement direction.

[0047] In a third aspect, an embodiment of the present disclosure also provides an electronic device, which includes:

[0048] A memory storing executable instructions;

[0049] A processor, and the processor runs the executable instructions in the memory to implement the above-mentioned method for predicting favorable geothermal areas based on numerical simulation.

[0050] In a fourth aspect, an embodiment of the present disclosure also provides a computer-readable storage medium, which stores a computer program, and when the computer program is executed by a processor, it implements the above-mentioned method for predicting favorable geothermal areas based on numerical simulation.

[0051] The method and device of the present invention have other characteristics and advantages, which will be obvious from the accompanying drawings incorporated herein and the subsequent specific embodiments, or will be described in detail in the accompanying drawings incorporated herein and the subsequent specific embodiments. These drawings and specific embodiments are jointly used to explain the specific principles of the present invention. BRIEF DESCRIPTION OF THE DRAWINGS

[0052] By describing the exemplary embodiments of the present invention in more detail in conjunction with the accompanying drawings, the above-mentioned and other objects, features, and advantages of the present invention will become more obvious. Among them, in the exemplary embodiments of the present invention, the same reference numerals generally represent the same components.

[0053] Figure 1 Shows a flowchart of the steps of the method for predicting favorable geothermal areas based on numerical simulation according to the present invention.

[0054] Figure 2 Shows a flowchart of the steps of the method for predicting favorable geothermal areas based on numerical simulation according to an embodiment of the present invention.

[0055] Figure 3Shows a schematic diagram of a seismic profile of a study area according to an embodiment of the present invention.

[0056] Figure 4a Shows a schematic diagram of a geological profile of the Heze Uplift - Yutai Sag according to an embodiment of the present invention. Figure 4b Shows a schematic diagram of a geological profile of the Yutai Sag - Dongmeng Uplift in the southwestern Shandong region according to an embodiment of the present invention. Figure 4c Shows a schematic diagram of a geological profile of the LQ85 - 378 survey line in the Shenxian - Yanggu - Shouzhang area according to an embodiment of the present invention. Figure 4d Shows a schematic diagram of a geological profile of the 99 - 436 survey line in the Shouzhang Sag according to an embodiment of the present invention.

[0057] Figure 5a Shows a schematic diagram of a temperature curve at a depth of 1 km in the geological profile of the Heze Uplift - Yutai Sag according to an embodiment of the present invention. Figure 5b Shows a schematic diagram of a temperature curve at a depth of 2 km in the geological profile of the Heze Uplift - Yutai Sag according to an embodiment of the present invention. Figure 5c Shows a schematic diagram of a temperature curve at a depth of 1 km in the geological profile of the Yutai Sag - Dongmeng Uplift in the southwestern Shandong region according to an embodiment of the present invention. Figure 5d Shows a schematic diagram of a temperature curve at a depth of 2 km in the geological profile of the Yutai Sag - Dongmeng Uplift in the southwestern Shandong region according to an embodiment of the present invention. Figure 5e Shows a schematic diagram of a temperature curve at a depth of 1 km in the LQ85 - 378 survey line in the Shenxian - Yanggu - Shouzhang area according to an embodiment of the present invention. Figure 5f Shows a schematic diagram of a temperature curve at a depth of 2 km in the LQ85 - 378 survey line in the Shenxian - Yanggu - Shouzhang area according to an embodiment of the present invention. Figure 5g Shows a schematic diagram of a temperature curve at a depth of 1 km in the geological profile of the 99 - 436 survey line in the Shouzhang Sag according to an embodiment of the present invention. Figure 5h Shows a schematic diagram of a temperature curve at a depth of 2 km in the geological profile of the 99 - 436 survey line in the Shouzhang Sag according to an embodiment of the present invention.

[0058] Figure 6a 、 Figure 6b Respectively show schematic diagrams of geothermal temperature plane maps at depths of 1 km and 2 km according to an embodiment of the present invention.

[0059] Figure 7a 、 Figure 7b Respectively show schematic diagrams of evaluation maps of favorable areas for geothermal resources at depths of 1 km and 2 km according to an embodiment of the present invention.

[0060] Figure 8The block diagram of a numerical simulation-based geothermal favorable area prediction device according to an embodiment of the present invention is shown.

[0061] Explanation of reference numerals:

[0062] 201, calculation module; 202, statistics module; 203, simulation module; 204, drawing module; 205, analysis module; 206, prediction module. Detailed implementation manners

[0063] The preferred embodiments of the present invention will be described in more detail below. Although the preferred embodiments of the present invention are described below, it should be understood that the present invention can be implemented in various forms and should not be limited by the embodiments set forth herein.

[0064] Figure 1 The flowchart showing the steps of a numerical simulation-based geothermal favorable area prediction method according to the present invention is shown.

[0065] As Figure 1 shown, the present invention provides a numerical simulation-based geothermal favorable area prediction method, including:

[0066] Step 101, determining the basement terrestrial heat flow in the study area; Step 102, establishing a two-dimensional geological model and statistically analyzing the thermal conductivity values of each lithology of each formation above the bedrock in the study area; Step 103, after assigning thermal conductivity to the two-dimensional geological model, performing a numerical simulation of heat conduction and saving the temperature curve at a fixed depth; Step 104, according to the temperature curve at the fixed depth, drawing a geothermal isotherm plan at the fixed depth in the study area; Step 105, analyzing the geothermal conditions in the study area based on basic geological data and geochemical data; Step 106, determining the correlation between the geothermal isotherm plan and the geothermal conditions, and predicting the geothermal favorable area in the study area.

[0067] In one example, the basement terrestrial heat flow is calculated by formula (1):

[0068] Q = -K·dT / dZ (1)

[0069] Wherein, Q is the basement terrestrial heat flow, K is the thermal conductivity of the basement rock, and dT / dZ is the geothermal gradient of the basement layer in the study area.

[0070] In one example, if the well is shallow and the temperature data of the basement layer is missing, the basement terrestrial heat flow is calculated by formula (2):

[0071] Q = Qb - A·Z (2)

[0072] Wherein, Q is the basement terrestrial heat flow, Qb is the surface heat flow, Z is the thickness of the sedimentary cover layer, and A is the heat generation rate of the sedimentary cover layer.

[0073] In one example, a two-dimensional geological model is established, and the thermal conductivity values of each lithology of each stratum above the bedrock in the study area are statistically obtained, including:

[0074] Through seismic interpretation, the seismic profiles and the thickness distributions of each stratum in the study area are obtained, and a two-dimensional geological model is established through AutoCAD. The two-dimensional geological model is imported into COMSOL software, and each stratum is meshed and assigned thermal conductivity data.

[0075] In one example, after the thermal conductivity of the two-dimensional geological model is assigned, a numerical simulation of heat conduction is carried out, and the fixed-depth temperature curve is saved:

[0076] The simulation is carried out using the heat conduction simulation module in COMSOL software, and the density is set to 2500 kg / m 3 , the specific heat is 1000 J / (kg·°C), the surface temperature is the average surface temperature of the study area, the initial geothermal gradient is 3 °C / km, a numerical simulation of heat conduction is carried out, and the temperature curve at a fixed depth is saved.

[0077] In one example, the plan of geothermal isotherms at a fixed depth in the study area includes:

[0078] According to the temperature curve at a fixed depth and combined with the temperature of the actual drilled wells, the plan of geothermal isotherms at a fixed depth in the study area is drawn using the Kriging interpolation method on Surfer software.

[0079] In one example, the geothermal conditions include the tectonic distribution characteristics, lithological characteristics, heat source location, and fluid movement direction.

[0080] Specifically, the basal terrestrial heat flow in the study area is obtained by calculation methods or literature research methods. The basal terrestrial heat flow can be calculated by formula (1). For the case where the basal stratum drilled is thin and temperature data is missing, the basal terrestrial heat flow can be calculated by formula (2).

[0081] The thermal conductivity values of each lithology of each stratum above the bedrock in the study area are statistically obtained through experiments or literature research. Through seismic interpretation, the seismic profiles and the thickness distributions of each stratum in the study area are obtained, and a two-dimensional geological model is established through AutoCAD. The model is imported into COMSOL software, and each stratum is meshed and assigned thermal conductivity data.

[0082] After the thermal conductivity of the model is assigned, the heat conduction simulation module in COMSOL software is used for simulation, and the density is set to 2500 kg / m 3 , the specific heat is 1000 J / (kg·°C), the surface temperature can be the average surface temperature of the study area, the initial geothermal gradient is 3 °C / km, the simulation time can be set to 400,000 years to achieve a steady-state effect, a numerical simulation of heat conduction is carried out, and the temperature curve at a fixed depth is obtained.

[0083] Based on the obtained temperature curve at a fixed depth and combined with the actual drilling temperature, use the Kriging interpolation method on Surfer software to draw the temperature contour plan of the fixed depth in the study area.

[0084] Combined with the basic geological data such as the crustal thermal structure, Moho depth, Curie surface data and the seismic profiles used in the study area, as well as geochemical data, analyze the geothermal conditions such as the tectonic distribution characteristics, lithological characteristics, heat source location, and fluid movement direction in the study area.

[0085] Combine the temperature contour plan and the geothermal conditions, determine the correlation between the two, and comprehensively predict the geothermal anomaly area in the study area.

[0086] The present invention also provides a device for predicting geothermal favorable areas based on numerical simulation, including:

[0087] A calculation module to determine the basal terrestrial heat flow in the study area;

[0088] A statistics module to establish a two-dimensional geological model and statistically calculate the thermal conductivity values of each lithology of each stratum above the bedrock in the study area;

[0089] A simulation module to perform thermal conductivity assignment on the two-dimensional geological model, then conduct numerical simulation of heat conduction, and save the temperature curve at a fixed depth;

[0090] A drawing module to draw the temperature contour plan of the fixed depth in the study area according to the temperature curve at a fixed depth;

[0091] An analysis module to analyze the geothermal conditions in the study area according to the basic geological data and geochemical data;

[0092] A prediction module to determine the correlation between the temperature contour plan and the geothermal conditions, and predict the geothermal favorable areas in the study area.

[0093] In one example, the basal terrestrial heat flow is calculated by formula (1):

[0094] Q = -K·dT / dZ (1)

[0095] Wherein, Q is the basal terrestrial heat flow, K is the thermal conductivity of the basal rock, and dT / dZ is the geothermal gradient of the basal layer in the study area.

[0096] In one example, if temperature data is missing, the basal terrestrial heat flow is calculated by formula (2):

[0097] Q = Qb - A·Z (2)

[0098] Wherein, Q is the basal terrestrial heat flow, Qb is the surface heat flow, Z is the thickness of the sedimentary cover, and A is the heat production rate of the sedimentary cover.

[0099] In one example, a two-dimensional geological model is established, and the thermal conductivity values of each lithology of each stratum above the bedrock in the study area are statistically analyzed, including:

[0100] Through seismic interpretation, the seismic profiles and the thickness distributions of each stratum in the study area are obtained. A two-dimensional geological model is established using AutoCAD, and the two-dimensional geological model is imported into COMSOL software to perform meshing and assign thermal conductivity data to each stratum.

[0101] In one example, after assigning thermal conductivity to the two-dimensional geological model, a numerical simulation of heat conduction is performed, and the temperature curve at a fixed depth is saved:

[0102] The simulation is carried out using the heat conduction simulation module in COMSOL software, with the density set to 2500 kg / m 3 , the specific heat is 1000 J / (kg·°C), the surface temperature is the average surface temperature of the study area, and the initial geothermal gradient is 3°C / km. A numerical simulation of heat conduction is performed to save the temperature curve at a fixed depth.

[0103] In one example, the isothermal contour plan at a fixed depth in the study area includes:

[0104] Based on the temperature curve at a fixed depth and combined with the temperature of the actual drill holes, the isothermal contour plan at a fixed depth in the study area is drawn using the Kriging interpolation method in Surfer software.

[0105] In one example, the geothermal conditions include the tectonic distribution characteristics, lithological characteristics, heat source location, and fluid movement direction.

[0106] Specifically, the basal terrestrial heat flow in the study area is obtained through calculation methods or literature research. The basal terrestrial heat flow can be calculated by formula (1). For the case where the basal stratum drilled through is thin and temperature data is missing, the basal terrestrial heat flow can be calculated by formula (2).

[0107] The thermal conductivity values of each lithology of each stratum above the bedrock in the study area are statistically analyzed through experiments or literature research. Through seismic interpretation, the seismic profiles and the thickness distributions of each stratum in the study area are obtained. A two-dimensional geological model is established using AutoCAD, and the model is imported into COMSOL software to perform meshing and assign thermal conductivity data to each stratum.

[0108] After assigning thermal conductivity to the model, the heat conduction simulation module in COMSOL software is used for simulation, with the density set to 2500 kg / m 3, with a specific heat capacity of 1000 J / (kg·℃), the surface temperature can be the average surface temperature of the study area, the initial geothermal gradient is 3℃ / km, and the simulation time can be set to 400,000 years to achieve a steady-state effect. Conduct a numerical simulation of heat conduction to obtain the temperature curve at a fixed depth.

[0109] Based on the obtained temperature curve at a fixed depth and combined with the actual drilling temperature, use the Kriging interpolation method in Surfer software to draw the temperature contour plan at a fixed depth in the study area.

[0110] Combined with the basic geological data such as the crustal heat structure, Moho depth, Curie depth data, and the seismic profiles used in the study area, as well as geochemical data, analyze the geothermal conditions such as the tectonic distribution characteristics, lithological characteristics, heat source location, and fluid movement direction in the study area.

[0111] Combine the temperature contour plan and geothermal conditions to determine the correlation between the two, clarify the mutual coupling relationship between the current geothermal field distribution characteristics, groundwater recharge sources, and tectonic characteristics in the region, clarify the main controlling factors of regional geothermal anomalies, and comprehensively predict the geothermal anomaly areas in the study area.

[0112] The present invention also provides an electronic device, which includes: a memory storing executable instructions; a processor that runs the executable instructions in the memory to implement the above-mentioned method for predicting favorable geothermal areas based on numerical simulation.

[0113] The present invention also provides a computer-readable storage medium, which stores a computer program that, when executed by a processor, implements the above-mentioned method for predicting favorable geothermal areas based on numerical simulation.

[0114] To facilitate understanding of the solutions and effects of the embodiments of the present invention, the following gives four specific application examples. Those skilled in the art should understand that this example is only for facilitating the understanding of the present invention, and any specific details are not intended to limit the present invention in any way.

[0115] Example 1

[0116] The specific example implemented this time is the regional geothermal numerical simulation under a large tectonic unit, and the experimental work area is the southern part of the Luxi Uplift Tectonic Area.

[0117] Figure 2 Shows a flowchart of the steps of a method for predicting favorable geothermal areas based on numerical simulation according to an embodiment of the present invention.

[0118] As Figure 2As shown in the figure, the method for predicting favorable geothermal areas based on numerical simulation includes: Step 101, determining the basement terrestrial heat flow in the study area; Step 102, establishing a two-dimensional geological model and statistically obtaining the thermal conductivity values of each lithology of each stratum above the bedrock in the study area; Step 103, after assigning thermal conductivity values to the two-dimensional geological model, conducting numerical simulation of heat conduction and saving the temperature curves at fixed depths; Step 104, according to the temperature curves at fixed depths, drawing the isogeothermal line plan of the study area at fixed depths; Step 105, analyzing the geothermal conditions in the study area based on basic geological data and geochemical data; Step 106, determining the correlation between the isogeothermal line plan and the geothermal conditions and predicting the favorable geothermal areas in the study area.

[0119] Through literature research, the terrestrial heat flow value at the bottom of sedimentary rocks in the Luxi Uplift area is 0.054 W / m 2 , and the thermal conductivity values of each lithology of each stratum above the bedrock in the study area are statistically obtained through experiments or literature research, as shown in Table 1.

[0120] Table 1

[0121]

[0122] Figure 3 Shows a schematic diagram of the seismic profile of the study area according to an embodiment of the present invention.

[0123] Figure 4a Shows a schematic diagram of the geological profile of the Heze Uplift - Yutai Sag according to an embodiment of the present invention, Figure 4b Shows a schematic diagram of the geological profile of the Yutai Sag - Dongmeng Uplift in the southwestern Shandong region according to an embodiment of the present invention, Figure 4c Shows a schematic diagram of the geological profile of the LQ85 - 378 survey line in the Shenxian - Yanggu - Shouzhang area according to an embodiment of the present invention, Figure 4d Shows a schematic diagram of the geological profile of the 99 - 436 survey line in the Shouzhang Sag according to an embodiment of the present invention.

[0124] Through seismic interpretation, the seismic profile of the study area and the thickness distribution of each stratum are obtained, as Figure 3 shown, and a two-dimensional geological model is established through AutoCAD, as Figures 4a - 4d shown, and the model is imported into the COMSOL software to perform meshing and assign thermal conductivity data to each stratum.

[0125] After assigning the thermal conductivity to the model, use the heat conduction simulation module in COMSOL software for simulation. Set the specific rock density and specific heat capacity. The surface temperature can be the average surface temperature of the study area, which is 25°C, and the initial geothermal gradient is 3°C / km. Simulate the steady-state effect, conduct numerical simulation of heat conduction, obtain the geothermal temperatures at depths of 1 km and 2 km for different wells, and draw the geothermal temperature curves at depths of 1 km and 2 km respectively.

[0126] Table 2

[0127]

[0128]

[0129] Figure 5a The schematic diagram shows the temperature curve graph at a depth of 1 km of the geological profile of Heze Uplift - Yutai Sag according to an embodiment of the present invention. Figure 5b The schematic diagram shows the temperature curve graph at a depth of 2 km of the geological profile of Heze Uplift - Yutai Sag according to an embodiment of the present invention. Figure 5c The schematic diagram shows the temperature curve graph at a depth of 1 km of the geological profile of Yutai Sag - Dongmeng Uplift in southwestern Shandong according to an embodiment of the present invention. Figure 5d The schematic diagram shows the temperature curve graph at a depth of 2 km of the geological profile of Yutai Sag - Dongmeng Uplift in southwestern Shandong according to an embodiment of the present invention. Figure 5e The schematic diagram shows the temperature curve graph at a depth of 1 km of the LQ85 - 378 survey line in Shenxian - Yanggu - Shouzhang area according to an embodiment of the present invention. Figure 5f The schematic diagram shows the temperature curve graph at a depth of 2 km of the LQ85 - 378 survey line in Shenxian - Yanggu - Shouzhang area according to an embodiment of the present invention. Figure 5g The schematic diagram shows the temperature curve graph at a depth of 1 km of the geological profile of the 99 - 436 survey line in Shouzhang Sag according to an embodiment of the present invention. Figure 5h The schematic diagram shows the temperature curve graph at a depth of 2 km of the geological profile of the 99 - 436 survey line in Shouzhang Sag according to an embodiment of the present invention.

[0130] Figure 6a 、 Figure 6b The schematic diagrams respectively show the geothermal temperature plan views at depths of 1 km and 2 km according to an embodiment of the present invention.

[0131] As Figures 5a - 5h shown, according to the geothermal temperature curves at depths of 1 km and 2 km, combined with the temperatures of actual drillings, use the Kriging interpolation method on the Surfer software to draw the geothermal temperature contour plan views at depths of 1 km and 2 km in the study area, as shown in Figure 6a 、 Figure 6b shown respectively.

[0132] Combined with the basic geological data such as the crustal thermal structure, Moho depth, Curie depth data, and the seismic profiles used in the study area, the sedimentary facies belt distribution data, and the groundwater chemical data, analyze the geothermal conditions such as the tectonic distribution characteristics, lithological characteristics, heat source location, and fluid movement direction in the study area.

[0133] Figure 7a 、 Figure 7b Schematically show the evaluation maps of favorable areas for geothermal resources at depths of 1 km and 2 km according to an embodiment of the present invention.

[0134] Combine the geothermal isotherm plane maps at depths of 1 km and 2 km with the geothermal conditions, determine the correlation between them, and comprehensively predict the geothermal anomaly areas in the study area, as Figure 7a 、 Figure 7b shown.

[0135] According to the seismic data, this method conducts seismic interpretation to obtain seismic profiles, establishes a geological model in AutoCAD, imports the geological model into COMSOL, obtains the main lithologies and their thermal conductivities of each stratigraphic age in the study area through experimental methods, and assigns values to the geological model in COMSOL to obtain a heat conduction model. Constraints are imposed on the heat conduction model according to the conditions such as the basal terrestrial heat flow value, density, specific heat capacity, initial geothermal gradient, and surface temperature in the study area, and two-dimensional heat conduction simulation is carried out. Finally, the obtained heat conduction results are combined with the actual temperature measurement data, and the Kriging interpolation method is used to draw the geothermal temperature field distribution map in the study area, and the favorable geothermal areas in the study area are comprehensively delineated in combination with data such as regional tectonics and groundwater chemistry, and the distribution law of the geothermal temperature field is studied to provide reference value for the later exploration and development of geothermal resources.

[0136] Example 2

[0137] Figure 8 Schematically show a block diagram of a device for predicting favorable geothermal areas based on numerical simulation according to an embodiment of the present invention.

[0138] As Figure 8 shown, the device for predicting favorable geothermal areas based on numerical simulation includes:

[0139] A calculation module 201 for determining the basal terrestrial heat flow in the study area;

[0140] A statistics module 202 for establishing a two-dimensional geological model and statistically obtaining the thermal conductivity values of each lithology of each stratigraphic layer above the bedrock in the study area;

[0141] A simulation module 203 for performing heat conduction numerical simulation after assigning thermal conductivity values to the two-dimensional geological model and saving the temperature curves at fixed depths;

[0142] The drawing module 204 draws the geothermal isotherm plane map at a fixed depth in the study area according to the temperature curve at the fixed depth.

[0143] The analysis module 205 analyzes the geothermal conditions in the study area according to the basic geological data and geochemical data.

[0144] The prediction module 206 determines the correlation between the geothermal isotherm plane map and the geothermal conditions, and predicts the favorable geothermal areas in the study area.

[0145] As an alternative, the basal terrestrial heat flow is calculated by formula (1):

[0146] Q = -K·dT / dZ (1)

[0147] where Q is the basal terrestrial heat flow, K is the thermal conductivity of the basal rock, and dT / dZ is the geothermal gradient of the basement layer in the study area.

[0148] As an alternative, if the well is shallow and the temperature data of the basement layer is missing, the basal terrestrial heat flow is calculated by formula (2):

[0149] Q = Qb - A·Z (2)

[0150] where Q is the basal terrestrial heat flow, Qb is the surface heat flow, Z is the thickness of the sedimentary cover, and A is the heat production rate of the sedimentary cover.

[0151] As an alternative, a two-dimensional geological model is established, and the thermal conductivity values of each lithology of each stratum above the bedrock in the study area are statistically included:

[0152] Through seismic interpretation, the seismic profiles and the thickness distributions of each stratum in the study area are obtained, and a two-dimensional geological model is established by AutoCAD. The two-dimensional geological model is imported into the COMSOL software, and the grids are generated for each stratum and the thermal conductivity data are assigned.

[0153] As an alternative, after the thermal conductivity is assigned to the two-dimensional geological model, a numerical simulation of heat conduction is carried out to save the temperature curve at a fixed depth:

[0154] The simulation is carried out using the heat conduction simulation module in the COMSOL software, with the density set to 2500 kg / m 3 , the specific heat is 1000 J / (kg·℃), the surface temperature is the average surface temperature of the study area, and the initial geothermal gradient is 3℃ / km. A numerical simulation of heat conduction is carried out to save the temperature curve at a fixed depth.

[0155] As an alternative, the drawing of the geothermal isotherm plane map at a fixed depth in the study area includes:

[0156] According to the temperature curve at a fixed depth and combined with the actual drilling temperature, the geothermal isotherm plane map at a fixed depth in the study area is drawn using Kriging interpolation method on Surfer software.

[0157] As an alternative, the geothermal conditions include tectonic distribution characteristics, lithological characteristics, heat source location, and fluid movement direction.

[0158] Example 3

[0159] The present disclosure provides an electronic device including: a memory storing executable instructions; and a processor that runs the executable instructions in the memory to implement the above-mentioned geothermal favorable area prediction method based on numerical simulation.

[0160] The electronic device according to an embodiment of the present disclosure includes a memory and a processor.

[0161] The memory is used to store non-transitory computer-readable instructions. Specifically, the memory may include one or more computer program products, and the computer program products may include various forms of computer-readable storage media, such as volatile memory and / or non-volatile memory. The volatile memory may include, for example, random access memory (RAM) and / or cache memory, etc. The non-volatile memory may include, for example, read-only memory (ROM), hard disk, flash memory, etc.

[0162] The processor may be a central processing unit (CPU) or other forms of processing units having data processing capabilities and / or instruction execution capabilities, and may control other components in the electronic device to perform desired functions. In an embodiment of the present disclosure, the processor is used to run the computer-readable instructions stored in the memory.

[0163] Those skilled in the art should understand that in order to solve the technical problem of how to obtain good user experience effects, this embodiment may also include well-known structures such as communication buses and interfaces, and these well-known structures should also be included in the protection scope of the present disclosure.

[0164] For the detailed description of this embodiment, reference may be made to the corresponding descriptions in the foregoing embodiments, and details will not be repeated here.

[0165] Example 4

[0166] An embodiment of the present disclosure provides a computer-readable storage medium storing a computer program, and when the computer program is executed by a processor, it implements the above-mentioned geothermal favorable area prediction method based on numerical simulation.

[0167] A computer-readable storage medium according to an embodiment of the present disclosure stores non-transitory computer-readable instructions. When the non-transitory computer-readable instructions are run by a processor, all or part of the steps of the methods of the various embodiments of the present disclosure described above are executed.

[0168] The above computer-readable storage medium includes, but is not limited to: optical storage media (such as CD-ROMs and DVDs), magneto-optical storage media (such as MOs), magnetic storage media (such as tapes or external hard drives), media with built-in rewritable non-volatile memories (such as memory cards), and media with built-in ROMs (such as ROM cartridges).

[0169] Those skilled in the art should understand that the purpose of the above description of the embodiments of the present invention is only to exemplarily illustrate the beneficial effects of the embodiments of the present invention, and is not intended to limit the embodiments of the present invention to any of the examples given.

[0170] The various embodiments of the present invention have been described above. The above description is exemplary and not exhaustive, and is also not limited to the disclosed embodiments. Many modifications and variations are obvious to those of ordinary skill in the art without departing from the scope and spirit of the described embodiments.

Claims

1. A method for predicting favorable geothermal areas based on numerical simulation, characterized in that Including: Determine the basal terrestrial heat flow in the study area; Establish a two-dimensional geological model and statistically obtain the thermal conductivity values of each lithology of each stratum above the bedrock in the study area; After assigning thermal conductivity values to the two-dimensional geological model, conduct numerical simulation of heat conduction and save the temperature curve at a fixed depth; Based on the temperature curve at the fixed depth, draw a plan of the isogeothermal lines at the fixed depth in the study area; Analyze the geothermal conditions in the study area according to the basic geological data and geochemical data; Determine the correlation between the plan of the isogeothermal lines and the geothermal conditions, and predict the favorable geothermal areas in the study area.

2. The method for predicting favorable geothermal areas based on numerical simulation according to claim 1, wherein, Calculate the basal terrestrial heat flow through formula (1): Q = -K·dT / dZ (1) Wherein, Q is the basal terrestrial heat flow, K is the thermal conductivity of the basal rock, and dT / dZ is the geothermal gradient of the basal layer in the study area.

3. The method for predicting geothermal favorable areas based on numerical simulation according to claim 1, wherein, If the well is shallow and the temperature data of the basal layer is missing, calculate the basal terrestrial heat flow through formula (2): Q = Qb - A·Z (2) Wherein, Q is the basal terrestrial heat flow, Qb is the surface heat flow, Z is the thickness of the sedimentary cover, and A is the heat generation rate of the sedimentary cover.

4. The numerical simulation-based geothermal favorable area prediction method according to claim 1, wherein, Establishing a two-dimensional geological model and statistically obtaining the thermal conductivity values of each lithology of each stratum above the bedrock in the study area includes: Through seismic interpretation, obtain the seismic profiles and the thickness distribution of each stratum in the study area, establish a two-dimensional geological model through AutoCAD, import the two-dimensional geological model into COMSOL software, and perform meshing and thermal conductivity data assignment for each stratum.

5. The method for predicting geothermal favorable areas based on numerical simulation according to claim 1, wherein, After assigning thermal conductivity values to the two-dimensional geological model, conduct numerical simulation of heat conduction and save the temperature curve at the fixed depth: The simulation is carried out using the heat conduction simulation module in COMSOL software, with the density set to 2500 kg / m 3 , the specific heat is 1000 J / (kg·℃), the surface temperature is the average surface temperature of the study area, and the initial geothermal gradient is 3℃ / km. A numerical simulation of heat conduction is performed, and the temperature curve at a fixed depth is saved.

6. The method for predicting geothermal favorable areas based on numerical simulation according to claim 1, wherein, Drawing a plan of the isogeothermal lines at the fixed depth in the study area includes: Based on the temperature curve at the fixed depth and combined with the temperature of the actual well, use the Kriging interpolation method on Surfer software to draw a plan of the isogeothermal lines at the fixed depth in the study area.

7. The method for predicting geothermal favorable areas based on numerical simulation according to claim 1, wherein, The geothermal conditions include the tectonic distribution characteristics, lithological characteristics, heat source location, and fluid movement direction.

8. A geothermal favorable area prediction device based on numerical simulation, characterized in that, Including: A calculation module to determine the basal terrestrial heat flow in the study area; A statistics module to establish a two-dimensional geological model and statistically obtain the thermal conductivity values of each lithology of each stratum above the bedrock in the study area; A simulation module to conduct numerical simulation of heat conduction after assigning thermal conductivity values to the two-dimensional geological model and save the temperature curve at a fixed depth; A drawing module to draw a plan of the isogeothermal lines at the fixed depth in the study area based on the temperature curve at the fixed depth; An analysis module to analyze the geothermal conditions in the study area according to the basic geological data and geochemical data; A prediction module to determine the correlation between the plan of the isogeothermal lines and the geothermal conditions and predict the favorable geothermal areas in the study area.

9. An electronic device, characterized in that, The electronic device includes: A memory storing executable instructions; A processor that runs the executable instructions in the memory to implement the method for predicting favorable geothermal areas based on numerical simulation 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, it implements the method for predicting favorable geothermal areas based on numerical simulation according to any one of claims 1-7.

Images