Method and device for predicting sweet section of uplift belt basement thermal reservoir, electronic equipment and medium

Abstract

The application discloses a raised zone base rock thermal reservoir sweet section prediction method and device, electronic equipment and medium. The method can include: respectively establishing a geological model of a raised zone base rock thermal reservoir and a top surface stratum temperature change model of the base rock thermal reservoir; establishing a numerical model of the base rock thermal reservoir according to the geological model of the raised zone base rock thermal reservoir; simulating a stratum temperature distribution field of the base rock raised zone according to the numerical model and the top surface stratum temperature change model of the base rock thermal reservoir; determining a geothermal sweet section and a geological sweet section of the base rock thermal reservoir, and further determining a development sweet section of the base rock thermal reservoir. The application analyzes the geothermal sweet section depth of the base rock raised zone in the shallow part of the sedimentary basin by the numerical simulation method, combines the geological sweet section analysis result of the base rock thermal reservoir, finds the depth sweet section in the geothermal exploration and development process, and effectively guides the geothermal field exploration and development.

Description

Technical Field

[0001] This invention relates to the field of geothermal exploration and development, and more specifically, to a method, apparatus, electronic device, and medium for predicting sweet sections of bedrock geothermal reservoirs in uplift zones. Background Technology

[0002] Geothermal resources can be classified into two types based on their tectonic setting and heat transfer mechanism: convection type in uplifted mountains and conduction type in sedimentary basins. Currently, most developed and utilized geothermal fields are of the conduction type in sedimentary basins. Among them, bedrock reservoirs (mainly carbonate rocks) in the uplifted zones within the basin have become the key targets for geothermal exploration and development due to their high temperature and ease of reinjection. Therefore, predicting the "sweet spot" (favorable depth) of bedrock reservoirs in uplifted zones within sedimentary basins is of significant production and economic importance for the efficient development and utilization of geothermal fields in the later stages.

[0003] Currently, exploration methods for geothermal resources in sedimentary basins mostly focus on analyzing the planar distribution of thermal anomalies using geophysical methods such as electrophysiology and magnetism, and further combining these with geological features such as structure, reservoir, and geothermal water migration to predict favorable zones. However, after delineating the thermal anomaly zone, there is little research on predicting the vertical "sweet spots" of the geothermal reservoir.

[0004] The low-to-medium temperature conductive geothermal resources within sedimentary basins are influenced by the depth of the Moho and the undulations of the basement. The higher the Moho, the greater the geothermal flow value of the basin basement, and the geothermal field in the basin caprock shows a significant positive correlation with the basement undulations. This is because, during the conduction of geothermal flow to the shallow surface, a "heat redistribution" phenomenon occurs due to the high thermal conductivity of the bedrock uplift zone. This involves "heat diversion" in the basin depression zone and "heat concentration" in the basin uplift zone, thus forming a positive geothermal anomaly in the shallow part of the bedrock uplift zone. However, the positive geothermal anomaly in the shallow part of the bedrock uplift zone gradually decreases with increasing depth until it disappears completely at a certain depth. Below this depth, it becomes a negative anomaly (i.e., the reservoir temperature is lower than the formation temperature at the same depth in the depression zone). Numerous geothermal well temperature measurements also confirm that the geothermal field in the bedrock uplift zone exhibits an "asymmetric mirror" distribution between the shallow positive anomaly and the deep negative anomaly in the vertical direction.

[0005] Therefore, it is necessary to develop a method, device, electronic equipment, and medium for predicting sweet sections of bedrock thermal reservoirs in uplift zones.

[0006] The information disclosed in the background section of this invention is intended only to enhance the understanding of the general background of this invention, and should not be construed as an admission or in any way implying that such information constitutes prior art known to those skilled in the art. Summary of the Invention

[0007] This invention proposes a method, device, electronic equipment, and medium for predicting the depth of geothermal sweet segments in bedrock geothermal reservoirs in uplift zones. It can analyze the depth of geothermal sweet segments in shallow bedrock uplift zones of sedimentary basins through numerical simulation, and combine the results of geological sweet segment analysis of bedrock geothermal reservoirs to find deep sweet segments in the process of geothermal exploration and development, effectively guiding the exploration and development of geothermal fields.

[0008] In a first aspect, embodiments of this disclosure provide a method for predicting sweet sections of bedrock geothermal reservoirs in uplift zones, including:

[0009] Geological models of the bedrock reservoir in the uplift zone and temperature variation models of the strata above the bedrock reservoir were established separately.

[0010] Based on the geological model of the bedrock thermal reservoir in the uplift zone, a numerical model of the bedrock thermal reservoir is established.

[0011] Based on the numerical model and the formation temperature change model at the top of the bedrock reservoir, the formation temperature distribution field of the bedrock uplift zone is simulated.

[0012] Identify the geothermal and geological sweet zones of the bedrock geothermal reservoir, and then determine the development sweet zones of the bedrock geothermal reservoir.

[0013] Preferably, establishing a numerical model of the bedrock reservoir based on the geological model of the uplift zone bedrock reservoir includes:

[0014] The geological model is then meshed.

[0015] The thermal conductivity range of different lithologies in the geological profile was determined by field outcrops or core sampling.

[0016] Thermal conductivity values ​​are assigned to different lithological combinations in the gridded formation to obtain a numerical model of the bedrock thermal reservoir.

[0017] Preferably, simulating the formation temperature distribution field of the bedrock uplift zone according to the numerical model includes:

[0018] Define the boundary conditions of the model and determine the numerical simulation formula;

[0019] For the numerical model, the implicit finite difference method is used for iterative solution, and forward modeling is used to obtain the formation temperature distribution at different times;

[0020] Determine whether the consistency rate between the output formation temperature lateral variation curve and the formation temperature variation model at the top of the bedrock reservoir reaches a set threshold. If yes, output the formation temperature field distribution map. If no, adjust the thermal conductivity parameters of the caprock and bedrock until the consistency rate reaches the set threshold.

[0021] Preferably, the numerical simulation formula is:

[0022]

[0023] Where ρ is the formation density, c is the specific heat capacity of the underground medium, T is the paleogeothermal temperature, i and j are the number of gridded cells in the x and z directions, respectively, x and z represent the horizontal and vertical coordinate axes, respectively, and n is the number of simulation time steps.

[0024] Preferably, the geothermal sweet section of the bedrock reservoir is identified as including:

[0025] Determine the depth of the maximum positive geothermal anomaly in the uplift zone;

[0026] Based on the formation temperature distribution field, geothermal wells in the uplift zone and the depression zone are selected respectively, and formation temperature-depth curves are plotted on the same coordinate axis. The depth of the intersection of the two curves is the heat conduction equilibrium depth line of the uplift zone and the depression zone.

[0027] The intersection point is located at the depth of the bedrock geothermal reservoir in the uplift zone, which is the high conductivity homogenization depth. The reservoir segment located between the maximum geothermal positive anomaly depth and the high conductivity homogenization depth is the geothermal sweet segment.

[0028] Preferably, the geological sweet sections of the bedrock geothermal reservoir include:

[0029] Identify the anhydrous argillaceous carbonate rock strata and the water-rich carbonate rock fractured strata.

[0030] Based on the porosity and permeability conditions of fractured sections in carbonate rocks, thermal reservoirs are classified into three types: Type I, Type II, and Type III. Type I reservoirs have an acoustic transit time greater than 200 μs / m, porosity greater than 9.5%, and permeability greater than 15 mD. Type II reservoirs have an acoustic transit time greater than 185 μs / m, porosity between 5.0% and 9.5%, and permeability between 4.0 mD and 15 mD. Type III reservoirs have an acoustic transit time less than 185 μs / m, porosity between 1.8% and 5%, and permeability between 0.1 mD and 4.0 mD.

[0031] The depth range where Class I and Class II reservoirs are located is the geological sweet section of the bedrock thermal reservoir.

[0032] Preferably, the overlapping section of the geothermal sweet section and the geological sweet section is the development sweet section of the bedrock thermal reservoir.

[0033] As one specific implementation of this disclosure,

[0034] Secondly, embodiments of this disclosure also provide a device for predicting sweet sections of bedrock geothermal reservoirs in uplift zones, comprising:

[0035] The model building module establishes a geological model of the bedrock reservoir in the uplift zone and a model of the temperature variation of the strata above the bedrock reservoir.

[0036] The numerical model building module establishes a numerical model of the bedrock thermal reservoir based on the geological model of the uplift zone bedrock thermal reservoir.

[0037] The simulation module simulates the formation temperature distribution field of the bedrock uplift zone based on the numerical model and the formation temperature change model of the bedrock reservoir top surface.

[0038] The sweet spot determination module identifies the geothermal sweet spots and geological sweet spots of the bedrock thermal reservoir, and then determines the development sweet spots of the bedrock thermal reservoir.

[0039] Preferably, establishing a numerical model of the bedrock reservoir based on the geological model of the uplift zone bedrock reservoir includes:

[0040] The geological model is then meshed.

[0041] The thermal conductivity range of different lithologies in the geological profile was determined by field outcrops or core sampling.

[0042] Thermal conductivity values ​​are assigned to different lithological combinations in the gridded formation to obtain a numerical model of the bedrock thermal reservoir.

[0043] Preferably, simulating the formation temperature distribution field of the bedrock uplift zone according to the numerical model includes:

[0044] Define the boundary conditions of the model and determine the numerical simulation formula;

[0045] For the numerical model, the implicit finite difference method is used for iterative solution, and forward modeling is used to obtain the formation temperature distribution at different times;

[0046] Determine whether the consistency rate between the output formation temperature lateral variation curve and the formation temperature variation model at the top of the bedrock reservoir reaches a set threshold. If yes, output the formation temperature field distribution map. If no, adjust the thermal conductivity parameters of the caprock and bedrock until the consistency rate reaches the set threshold.

[0047] Preferably, the numerical simulation formula is:

[0048]

[0049] Where ρ is the formation density, c is the specific heat capacity of the underground medium, T is the paleogeothermal temperature, i and j are the number of gridded cells in the x and z directions, respectively, x and z represent the horizontal and vertical coordinate axes, respectively, and n is the number of simulation time steps.

[0050] Preferably, the geothermal sweet section of the bedrock reservoir is identified as including:

[0051] Determine the depth of the maximum positive geothermal anomaly in the uplift zone;

[0052] Based on the formation temperature distribution field, geothermal wells in the uplift zone and the depression zone are selected respectively, and formation temperature-depth curves are plotted on the same coordinate axis. The depth of the intersection of the two curves is the heat conduction equilibrium depth line of the uplift zone and the depression zone.

[0053] The intersection point is located at the depth of the bedrock geothermal reservoir in the uplift zone, which is the high conductivity homogenization depth. The reservoir segment located between the maximum geothermal positive anomaly depth and the high conductivity homogenization depth is the geothermal sweet segment.

[0054] Preferably, the geological sweet sections of the bedrock geothermal reservoir include:

[0055] Identify the anhydrous argillaceous carbonate rock strata and the water-rich carbonate rock fractured strata.

[0056] Based on the porosity and permeability conditions of fractured sections in carbonate rocks, thermal reservoirs are classified into three types: Type I, Type II, and Type III. Type I reservoirs have an acoustic transit time greater than 200 μs / m, porosity greater than 9.5%, and permeability greater than 15 mD. Type II reservoirs have an acoustic transit time greater than 185 μs / m, porosity between 5.0% and 9.5%, and permeability between 4.0 mD and 15 mD. Type III reservoirs have an acoustic transit time less than 185 μs / m, porosity between 1.8% and 5%, and permeability between 0.1 mD and 4.0 mD.

[0057] The depth range where Class I and Class II reservoirs are located is the geological sweet section of the bedrock thermal reservoir.

[0058] Preferably, the overlapping section of the geothermal sweet section and the geological sweet section is the development sweet section of the bedrock thermal reservoir.

[0059] Thirdly, embodiments of this disclosure also provide an electronic device, the electronic device comprising:

[0060] Memory, which stores executable instructions;

[0061] A processor that executes the executable instructions in the memory to implement the method for predicting sweet segments of bedrock thermal reservoirs in uplift zones.

[0062] Fourthly, embodiments of this disclosure also provide a computer-readable storage medium storing a computer program that, when executed by a processor, implements the method for predicting sweet sections of bedrock thermal reservoirs in uplift zones.

[0063] Its beneficial effects are: to analyze the distribution pattern of geothermal field and heat transfer and heat accumulation mechanism in geothermal fields in sedimentary basins, thereby predicting the sweet spots in geothermal reservoirs, and to realize the integrated research from the selection of areas and zones in the early stage of geothermal resource exploration to the development of "sweet spots" in the later stage of geothermal field development, so as to obtain the relatively largest geothermal anomaly in the shallowest layer of the geothermal field and achieve the greatest economic benefits.

[0064] The methods and apparatus of the present invention have other features and advantages that will be apparent from or will be set forth in detail in the accompanying drawings and following detailed description, which together serve to explain the particular principles of the invention. Attached Figure Description

[0065] The above and other objects, features and advantages of the present invention will become more apparent from the more detailed description of exemplary embodiments of the invention in conjunction with the accompanying drawings, wherein the same reference numerals generally represent the same parts.

[0066] Figure 1 A flowchart illustrating the steps of a method for predicting sweet sections of bedrock thermal reservoirs in uplift zones according to an embodiment of the present invention is shown.

[0067] Figure 2a , Figure 2b and Figure 2c Schematic diagrams of a geological model, a formation temperature variation model, and a numerical model according to an embodiment of the present invention are shown respectively.

[0068] Figure 3a , Figure 3b , Figure 3c and Figure 3d The diagrams show the simulation results of the geothermal field for time periods t=1ka, t=100ka, t=200ka, and t=400ka according to an embodiment of the present invention, where k represents thousand and a represents year.

[0069] Figure 4 A schematic diagram of the maximum geothermal positive anomaly depth, high conductivity homogenization depth, and geothermal sweet spot is shown according to an embodiment of the present invention.

[0070] Figure 5 A schematic diagram of the superposition of geothermal sweet segments and geological sweet segments according to an embodiment of the present invention is shown.

[0071] Figure 6 A block diagram of a device for predicting sweet sections of bedrock thermal reservoirs in a protruding zone according to an embodiment of the present invention is shown.

[0072] Explanation of reference numerals in the attached figures:

[0073] 201. Model Building Module; 202. Numerical Model Building Module; 203. Simulation Module; 204. Sweet Segment Determination Module. Detailed Implementation

[0074] Preferred embodiments of the invention will now be described in more detail. While preferred embodiments of the invention are described below, it should be understood that the invention can be implemented in various forms and should not be limited to the embodiments set forth herein.

[0075] This invention provides a method for predicting sweet sections of bedrock geothermal reservoirs in uplift zones, comprising:

[0076] Geological models of the bedrock reservoir in the uplift zone and temperature variation models of the strata above the bedrock reservoir were established separately.

[0077] Specifically, a geological model of the bedrock geothermal reservoir in the uplift zone is established: by integrating the interpretation results of multiple seismic survey lines, a geological profile is compiled that traverses the basin uplift zone and reflects its main deformation characteristics, and the dividing boundary between the basin uplift zone and the depression zone is determined; the tectonic deformation patterns and dynamic mechanisms of the uplift zone are analyzed; and by combining well stratification data, the development depth of the bedrock geothermal reservoir and its spatial distribution characteristics are determined.

[0078] Establish a model for the temperature variation of the top layer of the bedrock geothermal reservoir: Using data such as well temperature measurement and combined with the regional caprock geothermal gradient variation law, establish the planar distribution characteristics of the temperature of the top layer of the bedrock in the basin uplift zone and adjacent areas, create typical geothermal geological profiles, and determine the lateral variation law of geothermal anomalies on the top layer of the bedrock in the uplift zone, which will serve as the verification condition for the geothermal field simulation of the uplift zone.

[0079] Based on the geological model of the bedrock reservoir in the uplift zone, a numerical model of the bedrock reservoir is established. In one example, the numerical model of the bedrock reservoir based on the geological model of the uplift zone includes:

[0080] The geological model was gridded.

[0081] The thermal conductivity range of different lithologies in the geological profile was determined by field outcrops or core sampling.

[0082] Thermal conductivity values ​​were assigned to different stratigraphic lithological combinations in a gridded manner to obtain a numerical model of the bedrock thermal reservoir.

[0083] Specifically, the geological model is gridded, that is, the grid nodes are divided into sections and blocks according to the thickness of strata of different ages revealed by the geological profile, and distinguished by different colors; the thermal conductivity range of different lithologies in the geological profile is determined by using field outcrops or well cores; thermal conductivity values ​​are assigned to different lithological combinations of gridded strata to obtain a numerical model of bedrock thermal reservoir.

[0084] Based on numerical models and formation temperature variation models at the top of bedrock reservoirs, the formation temperature distribution field of the bedrock uplift zone is simulated. In one example, the simulation of the formation temperature distribution field of the bedrock uplift zone based on the numerical model includes:

[0085] Define the boundary conditions of the model and determine the numerical simulation formula;

[0086] For the numerical model, the implicit finite difference method is used for iterative solution, and the formation temperature distribution at different times is obtained by forward modeling.

[0087] Determine whether the consistency rate between the output formation temperature lateral variation curve and the formation temperature variation model at the top of the bedrock reservoir reaches the set threshold. If yes, output the formation temperature field distribution map. If no, adjust the thermal conductivity parameters of the caprock and bedrock until the consistency rate reaches the set threshold.

[0088] In one example, the numerical simulation formula is:

[0089]

[0090] Where ρ is the formation density, c is the specific heat capacity of the underground medium, T is the paleogeothermal temperature, i and j are the number of gridded cells in the x and z directions, respectively, x and z represent the horizontal and vertical coordinate axes, respectively, and n is the number of simulation time steps.

[0091] Specifically, the model boundary conditions are set as follows: based on the geothermal geological characteristics of the study area, the initial surface temperature T0 and the initial geothermal gradient are set. Uniform heat flow q on the substrate.

[0092] Selection of numerical simulation formula: The two-dimensional heat conduction equation (i.e., Fourier heat conduction equation) without an additional heat source is adopted, and the formula is as follows:

[0093]

[0094] In the formula, T is the paleogeothermal temperature in °C, t is the time from the start of the simulation in seconds, and k is the kJ / ℃. x k z ρ is the thermal conductivity of the underground medium in the horizontal and vertical directions (W / (m·℃)), and ρ is the density of the formation (kg / m³). 3 ), c is the specific heat capacity of the underground medium / (J / (kg·℃)), x and z represent the horizontal and vertical coordinate axes respectively, and the boundary condition T| z＝0 =T0, q represents the homogenized heat flow value of the deep basement of the basin after the tectonic thermal event / (W / m) 2 ).

[0095] The solution to equation (2) can be obtained by using the finite difference method, which directly transforms the differential problem into an algebraic problem. According to the central second-order difference quotient method, equation (2) can be discretized into equation (1). Under the premise of setting the initial boundary conditions, the implicit finite difference method is used for iterative solution, which can obtain the formation temperature distribution at different times through forward modeling.

[0096] Output results verification: Output the formation temperature profile at the depth of the top surface of the bedrock in the uplift zone. Adjust the thermal conductivity parameters of the caprock and bedrock appropriately until the output formation temperature lateral variation curve and the measured variation characteristics of the formation temperature variation model at the top surface of the bedrock reservoir reach the set threshold. Output the formation temperature field distribution map of the uplift zone and its adjacent areas under given conditions.

[0097] Identify the geothermal and geological sweet zones of the bedrock geothermal reservoir, and then determine the development sweet zones of the bedrock geothermal reservoir.

[0098] In one example, identifying the geothermal sweet section of a bedrock reservoir includes:

[0099] Determine the depth of the maximum positive geothermal anomaly in the uplift zone;

[0100] Based on the formation temperature distribution field, geothermal wells in the uplift zone and the depression zone are selected respectively, and formation temperature-depth curves are plotted on the same coordinate axis. The depth of the intersection of the two curves is the heat conduction equilibrium depth line of the uplift zone and the depression zone.

[0101] The intersection point at the depth of the bedrock geothermal reservoir in the uplift zone is the high conductivity homogenization depth. The reservoir segment located between the maximum positive geothermal anomaly depth and the high conductivity homogenization depth is the geothermal sweet zone.

[0102] In one example, identifying the geological sweet section of a bedrock reservoir includes:

[0103] Identify the anhydrous argillaceous carbonate rock strata and the water-rich carbonate rock fractured strata.

[0104] Based on the porosity and permeability conditions of fractured sections in carbonate rocks, thermal reservoirs are classified into three types: Type I, Type II, and Type III. Type I reservoirs have an acoustic transit time greater than 200 μs / m, porosity greater than 9.5%, and permeability greater than 15 mD. Type II reservoirs have an acoustic transit time greater than 185 μs / m, porosity between 5.0% and 9.5%, and permeability between 4.0 mD and 15 mD. Type III reservoirs have an acoustic transit time less than 185 μs / m, porosity between 1.8% and 5%, and permeability between 0.1 mD and 4.0 mD.

[0105] The depth range where Class I and Class II reservoirs are located is the geological sweet section of the bedrock thermal reservoir.

[0106] In one example, the overlapping section of the geothermal sweet section and the geological sweet section is used as the bedrock geothermal reservoir development sweet section.

[0107] Specifically, the determination of the geothermal sweet zone of the bedrock reservoir: (1) Determination of the maximum positive geothermal anomaly depth in the uplift zone: The maximum positive geothermal anomaly depth refers to the depth at which the temperature difference between the bedrock strata in the uplift zone and the strata in the depression zone is the greatest. Under normal circumstances, the maximum positive geothermal anomaly depth is located at the top of the bedrock reservoir, that is, the lithological boundary between the overlying thermal cap and the bedrock reservoir in the uplift zone. (2) Identification of the thermal conduction equilibrium depth line: At a certain depth line in the uplift zone and the depression zone, the strata temperatures of the two are equal and the average geothermal gradient is consistent. This indicates that within the vertical interval from the depth line to the surface, the uplift zone and the depression zone have the same thermal conductivity. (3) Determination of the positive anomaly section of the bedrock reservoir, i.e., the geothermal sweet zone: Select typical geothermal wells (or artificial wells) in the uplift zone and the depression zone respectively, and draw formation temperature-depth curves on the same coordinate axis. The depth of the intersection of the two curves is the thermal conductivity equilibrium depth line of the uplift zone and the depression zone. The depth of this intersection point in the bedrock reservoir of the uplift zone is called the high conductivity homogenization depth. The reservoir section located between the maximum geothermal positive anomaly depth (or the top of the bedrock reservoir) and the high conductivity homogenization depth is the geothermal sweet zone.

[0108] Determination of the geological sweet section of bedrock thermal reservoir: (1) Determine the reservoir section and non-reservoir section: Based on the interpretation results of geothermal well logging curves, identify the waterless argillaceous carbonate rock section and the water-rich carbonate rock fracture section; classify the sections with natural gamma (GR) less than 10 API, porosity greater than 1.8%, and permeability greater than 0.1 mD as water-rich effective thermal reservoir sections; (2) Classification of reservoir performance: Based on the quality of the porosity and permeability conditions of the carbonate rock fracture section, classify the thermal reservoir into three types from good to poor; mainly based on sonic transit time, porosity and permeability. The reservoirs are classified according to three indicators: Class I reservoirs have a sonic transit time greater than 200 μs / m, porosity greater than 9.5%, and permeability greater than 15 mD; Class II reservoirs have a sonic transit time greater than 185 μs / m, porosity between 5.0% and 9.5%, and permeability between 4.0 mD and 15 mD; Class III reservoirs have a sonic transit time less than 185 μs / m, porosity between 1.8% and 5%, and permeability between 0.1 mD and 4.0 mD; (3) Determination of reservoir depth: The depth range where Class I and Class II reservoirs are located is the geological depth range of bedrock thermal reservoirs.

[0109] Determination of the sweet zone for bedrock geothermal reservoir development: The overlapping section of the geothermal sweet zone and the geological sweet zone is the sweet zone for bedrock geothermal reservoir development.

[0110] The present invention also provides a device for predicting sweet sections of bedrock thermal reservoirs in uplift zones, comprising:

[0111] The model building module establishes a geological model of the bedrock reservoir in the uplift zone and a model of the temperature variation of the strata above the bedrock reservoir.

[0112] Specifically, a geological model of the bedrock geothermal reservoir in the uplift zone is established: by integrating the interpretation results of multiple seismic survey lines, a geological profile is compiled that traverses the basin uplift zone and reflects its main deformation characteristics, and the dividing boundary between the basin uplift zone and the depression zone is determined; the tectonic deformation patterns and dynamic mechanisms of the uplift zone are analyzed; and by combining well stratification data, the development depth of the bedrock geothermal reservoir and its spatial distribution characteristics are determined.

[0113] Establish a model for the temperature variation of the top layer of the bedrock geothermal reservoir: Using data such as well temperature measurement and combined with the regional caprock geothermal gradient variation law, establish the planar distribution characteristics of the temperature of the top layer of the bedrock in the basin uplift zone and adjacent areas, create typical geothermal geological profiles, and determine the lateral variation law of geothermal anomalies on the top layer of the bedrock in the uplift zone, which will serve as the verification condition for the geothermal field simulation of the uplift zone.

[0114] The numerical model building module establishes a numerical model of the bedrock geothermal reservoir based on the geological model of the uplift zone bedrock geothermal reservoir. In one example, establishing the numerical model of the bedrock geothermal reservoir based on the geological model of the uplift zone bedrock geothermal reservoir includes:

[0115] The geological model was gridded.

[0116] The thermal conductivity range of different lithologies in the geological profile was determined by field outcrops or core sampling.

[0117] Thermal conductivity values ​​were assigned to different stratigraphic lithological combinations in a gridded manner to obtain a numerical model of the bedrock thermal reservoir.

[0118] Specifically, the geological model is gridded, that is, the grid nodes are divided into sections and blocks according to the thickness of strata of different ages revealed by the geological profile, and distinguished by different colors; the thermal conductivity range of different lithologies in the geological profile is determined by using field outcrops or well cores; thermal conductivity values ​​are assigned to different lithological combinations of gridded strata to obtain a numerical model of bedrock thermal reservoir.

[0119] The simulation module simulates the formation temperature distribution field of the bedrock uplift zone based on the numerical model and the formation temperature variation model of the bedrock reservoir top surface. In one example, the simulation of the formation temperature distribution field of the bedrock uplift zone based on the numerical model includes:

[0120] Define the boundary conditions of the model and determine the numerical simulation formula;

[0121] For the numerical model, the implicit finite difference method is used for iterative solution, and the formation temperature distribution at different times is obtained by forward modeling.

[0122] Determine whether the consistency rate between the output formation temperature lateral variation curve and the formation temperature variation model at the top of the bedrock reservoir reaches the set threshold. If yes, output the formation temperature field distribution map. If no, adjust the thermal conductivity parameters of the caprock and bedrock until the consistency rate reaches the set threshold.

[0123] In one example, the numerical simulation formula is:

[0124]

[0125] Where ρ is the formation density, c is the specific heat capacity of the underground medium, T is the paleogeothermal temperature, i and j are the number of gridded cells in the x and z directions, respectively, x and z represent the horizontal and vertical coordinate axes, respectively, and n is the number of simulation time steps.

[0126] Specifically, the model boundary conditions are set as follows: based on the geothermal geological characteristics of the study area, the initial surface temperature T0 and the initial geothermal gradient are set. Uniform heat flow q on the substrate.

[0127] Selection of numerical simulation formula: The two-dimensional heat conduction equation (i.e., Fourier heat conduction equation) without an additional heat source is adopted, and the formula is as follows:

[0128]

[0129] In the formula, T is the paleogeothermal temperature in °C, t is the time from the start of the simulation in seconds, and k is the kJ / ℃. x k z ρ is the thermal conductivity of the underground medium in the horizontal and vertical directions (W / (m·℃)), and ρ is the density of the formation (kg / m³). 3 ), c is the specific heat capacity of the underground medium / (J / (kg·℃)), x and z represent the horizontal and vertical coordinate axes respectively, and the boundary condition T| z＝0 =T0, q represents the homogenized heat flow value of the deep basement of the basin after the tectonic thermal event / (W / m) 2 ).

[0130] The solution to equation (2) can be obtained by using the finite difference method, which directly transforms the differential problem into an algebraic problem. According to the central second-order difference quotient method, equation (2) can be discretized into equation (1). Under the premise of setting the initial boundary conditions, the implicit finite difference method is used for iterative solution, which can obtain the formation temperature distribution at different times through forward modeling.

[0131] Output results verification: Output the formation temperature profile at the depth of the top surface of the bedrock in the uplift zone. Adjust the thermal conductivity parameters of the caprock and bedrock appropriately until the output formation temperature lateral variation curve and the measured variation characteristics of the formation temperature variation model at the top surface of the bedrock reservoir reach the set threshold. Output the formation temperature field distribution map of the uplift zone and its adjacent areas under given conditions.

[0132] The sweet spot determination module identifies the geothermal sweet spots and geological sweet spots of the bedrock reservoir, thereby determining the development sweet spots for the bedrock reservoir. In one example, determining the geothermal sweet spots of the bedrock reservoir includes:

[0133] Determine the depth of the maximum positive geothermal anomaly in the uplift zone;

[0134] Based on the formation temperature distribution field, geothermal wells in the uplift zone and the depression zone are selected respectively, and formation temperature-depth curves are plotted on the same coordinate axis. The depth of the intersection of the two curves is the heat conduction equilibrium depth line of the uplift zone and the depression zone.

[0135] The intersection point at the depth of the bedrock geothermal reservoir in the uplift zone is the high conductivity homogenization depth. The reservoir segment located between the maximum positive geothermal anomaly depth and the high conductivity homogenization depth is the geothermal sweet zone.

[0136] In one example, identifying the geological sweet section of a bedrock reservoir includes:

[0137] Identify the anhydrous argillaceous carbonate rock strata and the water-rich carbonate rock fractured strata.

[0138] Based on the porosity and permeability conditions of fractured sections in carbonate rocks, thermal reservoirs are classified into three types: Type I, Type II, and Type III. Type I reservoirs have an acoustic transit time greater than 200 μs / m, porosity greater than 9.5%, and permeability greater than 15 mD. Type II reservoirs have an acoustic transit time greater than 185 μs / m, porosity between 5.0% and 9.5%, and permeability between 4.0 mD and 15 mD. Type III reservoirs have an acoustic transit time less than 185 μs / m, porosity between 1.8% and 5%, and permeability between 0.1 mD and 4.0 mD.

[0139] The depth range where Class I and Class II reservoirs are located is the geological sweet section of the bedrock thermal reservoir.

[0140] In one example, the overlapping section of the geothermal sweet section and the geological sweet section is used as the bedrock geothermal reservoir development sweet section.

[0141] Specifically, the determination of the geothermal sweet zone of the bedrock reservoir: (1) Determination of the maximum positive geothermal anomaly depth in the uplift zone: The maximum positive geothermal anomaly depth refers to the depth at which the temperature difference between the bedrock strata in the uplift zone and the strata in the depression zone is the greatest. Under normal circumstances, the maximum positive geothermal anomaly depth is located at the top of the bedrock reservoir, that is, the lithological boundary between the overlying thermal cap and the bedrock reservoir in the uplift zone. (2) Identification of the thermal conduction equilibrium depth line: At a certain depth line in the uplift zone and the depression zone, the strata temperatures of the two are equal and the average geothermal gradient is consistent. This indicates that within the vertical interval from the depth line to the surface, the uplift zone and the depression zone have the same thermal conductivity. (3) Determination of the positive anomaly section of the bedrock reservoir, i.e., the geothermal sweet zone: Select typical geothermal wells (or artificial wells) in the uplift zone and the depression zone respectively, and draw formation temperature-depth curves on the same coordinate axis. The depth of the intersection of the two curves is the thermal conductivity equilibrium depth line of the uplift zone and the depression zone. The depth of this intersection point in the bedrock reservoir of the uplift zone is called the high conductivity homogenization depth. The reservoir section located between the maximum geothermal positive anomaly depth (or the top of the bedrock reservoir) and the high conductivity homogenization depth is the geothermal sweet zone.

[0142] Determination of the geological sweet section of bedrock thermal reservoir: (1) Determine the reservoir section and non-reservoir section: Based on the interpretation results of geothermal well logging curves, identify the waterless argillaceous carbonate rock section and the water-rich carbonate rock fracture section; classify the sections with natural gamma (GR) less than 10 API, porosity greater than 1.8%, and permeability greater than 0.1 mD as water-rich effective thermal reservoir sections; (2) Classification of reservoir performance: Based on the quality of the porosity and permeability conditions of the carbonate rock fracture section, classify the thermal reservoir into three types from good to poor; mainly based on sonic transit time, porosity and permeability. The reservoirs are classified according to three indicators: Class I reservoirs have a sonic transit time greater than 200 μs / m, porosity greater than 9.5%, and permeability greater than 15 mD; Class II reservoirs have a sonic transit time greater than 185 μs / m, porosity between 5.0% and 9.5%, and permeability between 4.0 mD and 15 mD; Class III reservoirs have a sonic transit time less than 185 μs / m, porosity between 1.8% and 5%, and permeability between 0.1 mD and 4.0 mD; (3) Determination of reservoir depth: The depth range where Class I and Class II reservoirs are located is the geological depth range of bedrock thermal reservoirs.

[0143] Determination of the sweet zone for bedrock geothermal reservoir development: The overlapping section of the geothermal sweet zone and the geological sweet zone is the sweet zone for bedrock geothermal reservoir development.

[0144] The present invention also provides an electronic device, comprising: a memory storing executable instructions; and a processor that executes the executable instructions in the memory to implement the above-described method for predicting sweet sections of bedrock thermal reservoirs in uplift zones.

[0145] The present invention also provides a computer-readable storage medium storing a computer program that, when executed by a processor, implements the above-described method for predicting sweet sections of bedrock thermal reservoirs in uplift zones.

[0146] To facilitate understanding of the solutions and effects of the embodiments of the present invention, four specific application examples are given below. Those skilled in the art should understand that these examples are merely for the purpose of understanding the present invention, and any specific details therein are not intended to limit the present invention in any way.

[0147] Example 1

[0148] Figure 1 A flowchart illustrating the steps of a method for predicting sweet sections of bedrock thermal reservoirs in uplift zones according to an embodiment of the present invention is shown.

[0149] like Figure 1 As shown, the method for predicting the sweet spot of the bedrock geothermal reservoir in the uplift zone includes: Step 101, establishing a geological model of the bedrock geothermal reservoir in the uplift zone and a temperature change model of the strata above the bedrock geothermal reservoir; Step 102, establishing a numerical model of the bedrock geothermal reservoir based on the geological model of the bedrock geothermal reservoir in the uplift zone; Step 103, simulating the strata temperature distribution field of the bedrock uplift zone based on the numerical model and the temperature change model of the strata above the bedrock geothermal reservoir; Step 104, determining the geothermal sweet spot and geological sweet spot of the bedrock geothermal reservoir, and then determining the development sweet spot of the bedrock geothermal reservoir.

[0150] Figure 2a , Figure 2b and Figure 2c Schematic diagrams of a geological model, a formation temperature variation model, and a numerical model according to an embodiment of the present invention are shown respectively.

[0151] like Figure 2a As shown, a geological model of the bedrock geothermal reservoir in the uplift zone is established: by integrating the interpretation results of multiple seismic survey lines, a geological profile is compiled that traverses the basin uplift zone and reflects its main deformation characteristics, and the dividing boundary between the basin uplift zone and the depression zone is determined; the tectonic deformation patterns and dynamic mechanisms of the uplift zone are analyzed; and by combining well stratification data, the development depth of the bedrock geothermal reservoir and its spatial distribution characteristics are determined.

[0152] like Figure 2b As shown, a model for the temperature variation of the top layer of the bedrock reservoir was established: using data such as well temperature measurement and combined with the regional caprock geothermal gradient variation law, the planar distribution characteristics of the temperature of the top layer of the bedrock in the basin uplift zone and adjacent areas were established, typical geothermal geological profiles were made, and the lateral variation law of geothermal anomalies on the top layer of the bedrock in the uplift zone was determined, which served as the verification condition for the geothermal field simulation of the uplift zone.

[0153] The geological model is gridded, meaning that based on the thickness of strata of different ages revealed by the geological profile, the grid nodes are divided into sections and blocks according to a rectangular step size ratio, and distinguished by different colors. Thermal conductivity ranges of different lithologies in the geological profile are determined using field outcrops or well core samples. Thermal conductivity values ​​are assigned to different lithological combinations within the gridded structure to obtain a numerical model of the bedrock reservoir, such as... Figure 2c As shown.

[0154] Model boundary conditions are set: Based on the geothermal geological characteristics of the study area, the initial surface temperature T0 and the initial geothermal gradient are set. Uniform heat flow q on the substrate.

[0155] Figure 3a , Figure 3b , Figure 3c and Figure 3d The diagrams show the simulation results of the geothermal field for time periods t=1ka, t=100ka, t=200ka, and t=400ka according to an embodiment of the present invention, where k represents thousand and a represents year.

[0156] The selection of numerical simulation formulas: For the solution of formula (2), the finite difference method can be used to directly transform the differential problem into an algebraic problem. According to the central second-order difference quotient method, formula (2) can be discretized into formula (1). Under the premise of setting the starting boundary conditions, the implicit finite difference method is used for iterative solution, that is, the formation temperature distribution at different times can be obtained by forward modeling, such as Figures 3a-3d As shown.

[0157] Output results verification: Output the formation temperature profile at the depth of the top surface of the bedrock in the uplift zone. Adjust the thermal conductivity parameters of the caprock and bedrock appropriately until the output formation temperature lateral variation curve and the measured variation characteristics of the formation temperature variation model at the top surface of the bedrock reservoir reach the set threshold. Output the formation temperature field distribution map of the uplift zone and its adjacent areas under given conditions.

[0158] Figure 4 A schematic diagram of the maximum geothermal positive anomaly depth, high conductivity homogenization depth, and geothermal sweet spot is shown according to an embodiment of the present invention.

[0159] like Figure 4As shown, the determination of the geothermal sweet section of the bedrock reservoir is as follows: (1) Determination of the maximum positive geothermal anomaly depth in the uplift zone: The maximum positive geothermal anomaly depth refers to the depth at which the temperature difference between the bedrock strata in the uplift zone and the strata in the depression zone is the greatest. Under normal circumstances, the maximum positive geothermal anomaly depth is located at the top of the bedrock reservoir, that is, the lithological boundary between the overlying thermal cap and the bedrock reservoir in the uplift zone. (2) Identification of the thermal conduction equilibrium depth line: At a certain depth line in the uplift zone and the depression zone, the strata temperatures of the two are equal and the average geothermal gradient is consistent. This indicates that within the vertical interval from the depth line to the surface, the uplift zone and the depression zone have the same thermal conductivity. (3) Determination of the positive anomaly section of the bedrock reservoir, i.e., the geothermal sweet zone: Select typical geothermal wells (or artificial wells) in the uplift zone and the depression zone respectively, and draw formation temperature-depth curves on the same coordinate axis. The depth of the intersection of the two curves is the thermal conductivity equilibrium depth line of the uplift zone and the depression zone. The depth of this intersection point in the bedrock reservoir of the uplift zone is called the high conductivity homogenization depth. The reservoir section located between the maximum geothermal positive anomaly depth (or the top of the bedrock reservoir) and the high conductivity homogenization depth is the geothermal sweet zone.

[0160] Determination of the geological sweet section of bedrock thermal reservoir: (1) Determine the reservoir section and non-reservoir section: Based on the interpretation results of geothermal well logging curves, identify the waterless argillaceous carbonate rock section and the water-rich carbonate rock fracture section; classify the sections with natural gamma (GR) less than 10 API, porosity greater than 1.8%, and permeability greater than 0.1 mD as water-rich effective thermal reservoir sections; (2) Classification of reservoir performance: Based on the quality of the porosity and permeability conditions of the carbonate rock fracture section, classify the thermal reservoir into three types from good to poor; mainly based on sonic transit time, porosity and permeability. The reservoirs are classified according to three indicators: Class I reservoirs have a sonic transit time greater than 200 μs / m, porosity greater than 9.5%, and permeability greater than 15 mD; Class II reservoirs have a sonic transit time greater than 185 μs / m, porosity between 5.0% and 9.5%, and permeability between 4.0 mD and 15 mD; Class III reservoirs have a sonic transit time less than 185 μs / m, porosity between 1.8% and 5%, and permeability between 0.1 mD and 4.0 mD; (3) Determination of reservoir depth: The depth range where Class I and Class II reservoirs are located is the geological depth range of bedrock thermal reservoirs.

[0161] Figure 5 A schematic diagram of the superposition of geothermal sweet segments and geological sweet segments according to an embodiment of the present invention is shown.

[0162] Determining the sweet zone for bedrock geothermal reservoir development: The overlapping section of the geothermal sweet zone and the geological sweet zone is the sweet zone for bedrock geothermal reservoir development, such as... Figure 5 As shown.

[0163] Example 2

[0164] Figure 6A block diagram of a device for predicting sweet sections of bedrock thermal reservoirs in a protruding zone according to an embodiment of the present invention is shown.

[0165] like Figure 6 As shown, the device for predicting sweet sections of bedrock geothermal reservoirs in the uplift zone includes:

[0166] Model building module 201 establishes a geological model of the bedrock reservoir in the uplift zone and a temperature variation model of the strata above the bedrock reservoir.

[0167] Numerical model building module 202 establishes a numerical model of the bedrock thermal reservoir based on the geological model of the bedrock thermal reservoir in the uplift zone.

[0168] Simulation module 203 simulates the formation temperature distribution field of the bedrock uplift zone based on the numerical model and the formation temperature change model of the bedrock reservoir top surface.

[0169] The sweet spot determination module 204 determines the geothermal sweet spot and geological sweet spot of the bedrock thermal reservoir, and then determines the development sweet spot of the bedrock thermal reservoir.

[0170] As an optional approach, based on the geological model of the bedrock reservoir in the uplift zone, a numerical model of the bedrock reservoir can be established, including:

[0171] The geological model was gridded.

[0172] The thermal conductivity range of different lithologies in the geological profile was determined by field outcrops or core sampling.

[0173] Thermal conductivity values ​​were assigned to different stratigraphic lithological combinations in a gridded manner to obtain a numerical model of the bedrock thermal reservoir.

[0174] As an optional approach, the formation temperature distribution field of the bedrock uplift zone can be simulated using a numerical model, including:

[0175] Define the boundary conditions of the model and determine the numerical simulation formula;

[0176] For the numerical model, the implicit finite difference method is used for iterative solution, and the formation temperature distribution at different times is obtained by forward modeling.

[0177] Determine whether the consistency rate between the output formation temperature lateral variation curve and the formation temperature variation model at the top of the bedrock reservoir reaches the set threshold. If yes, output the formation temperature field distribution map. If no, adjust the thermal conductivity parameters of the caprock and bedrock until the consistency rate reaches the set threshold.

[0178] As an alternative, the numerical simulation formula is:

[0179]

[0180] Where ρ is the formation density, c is the specific heat capacity of the underground medium, T is the paleogeothermal temperature, i and j are the number of gridded cells in the x and z directions, respectively, x and z represent the horizontal and vertical coordinate axes, respectively, and n is the number of simulation time steps.

[0181] As an alternative, the geothermal sweet sections of the bedrock reservoir include:

[0182] Determine the depth of the maximum positive geothermal anomaly in the uplift zone;

[0183] Based on the formation temperature distribution field, geothermal wells in the uplift zone and the depression zone are selected respectively, and formation temperature-depth curves are plotted on the same coordinate axis. The depth of the intersection of the two curves is the heat conduction equilibrium depth line of the uplift zone and the depression zone.

[0184] The intersection point at the depth of the bedrock geothermal reservoir in the uplift zone is the high conductivity homogenization depth. The reservoir segment located between the top of the bedrock geothermal reservoir and the high conductivity homogenization depth is the geothermal sweet zone.

[0185] As an alternative, identifying the geological sweet sections of bedrock geothermal reservoirs includes:

[0186] Identify the anhydrous argillaceous carbonate rock strata and the water-rich carbonate rock fractured strata.

[0187] Based on the thickness of the fractured section of carbonate rock and the porosity and permeability conditions, thermal reservoirs are classified into three types: good, medium, and poor.

[0188] The depth range where the good and medium reservoirs are located is the geological sweet section of the bedrock thermal reservoir.

[0189] As an alternative, the overlapping section of the geothermal sweet section and the geological sweet section is used as the bedrock thermal reservoir development sweet section.

[0190] Example 3

[0191] This disclosure provides an electronic device comprising: a memory storing executable instructions; and a processor executing the executable instructions in the memory to implement the above-described method for predicting sweet sections of bedrock thermal reservoirs in uplift zones.

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

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

[0194] The processor may be a central processing unit (CPU) or other form of processing unit with data processing capabilities and / or instruction execution capabilities, and may control other components in the electronic device to perform desired functions. In one embodiment of this disclosure, the processor is used to execute computer-readable instructions stored in the memory.

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

[0196] For a detailed description of this embodiment, please refer to the corresponding descriptions in the foregoing embodiments, which will not be repeated here.

[0197] Example 4

[0198] This disclosure provides a computer-readable storage medium storing a computer program that, when executed by a processor, implements the method for predicting sweet sections of bedrock thermal reservoirs in uplift zones.

[0199] A computer-readable storage medium according to embodiments of the present disclosure stores non-transitory computer-readable instructions. When these non-transitory computer-readable instructions are executed by a processor, all or part of the steps of the methods described in the foregoing embodiments of the present disclosure are performed.

[0200] The aforementioned computer-readable storage media include, but are not limited to: optical storage media (e.g., CD-ROM and DVD), magneto-optical storage media (e.g., MO), magnetic storage media (e.g., magnetic tape or portable hard drive), media with built-in rewritable non-volatile memory (e.g., memory card), and media with built-in ROM (e.g., ROM cartridge).

[0201] Those skilled in the art should understand that the above description of the embodiments of the present invention is only intended to 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.

[0202] The various embodiments of the present invention have been described above. These descriptions are exemplary and not exhaustive, nor are they limited to the disclosed embodiments. Many modifications and variations will be apparent to those skilled in the art without departing from the scope and spirit of the described embodiments.

Claims

1. A method for predicting sweet sections of bedrock thermal reservoirs in uplift zones, characterized in that, include: Geological models of the bedrock reservoir in the uplift zone and temperature variation models of the strata above the bedrock reservoir were established separately. Based on the geological model of the bedrock thermal reservoir in the uplift zone, a numerical model of the bedrock thermal reservoir is established. Based on the numerical model and the formation temperature change model at the top of the bedrock reservoir, the formation temperature distribution field of the bedrock uplift zone is simulated. Identify the geothermal and geological sweet zones of the bedrock thermal reservoir, and then determine the development sweet zones of the bedrock thermal reservoir; Among them, the geothermal sweet sections that have been identified as bedrock thermal reservoirs include: Determine the depth of the maximum positive geothermal anomaly in the uplift zone; Based on the formation temperature distribution field, geothermal wells in the uplift zone and the depression zone are selected respectively, and formation temperature-depth curves are plotted on the same coordinate axis. The depth of the intersection of the two curves is the heat conduction equilibrium depth line of the uplift zone and the depression zone. The intersection point is located at the depth of the bedrock thermal reservoir in the uplift zone, which is the high conductivity homogenization depth. The reservoir segment located between the maximum positive geothermal anomaly depth and the high conductivity homogenization depth is the geothermal sweet segment. Among them, the geological sweet sections that have been identified as bedrock geothermal reservoirs include: Identify anhydrous argillaceous carbonate rock sections and water-rich carbonate rock fracture sections; Based on the porosity and permeability conditions of fractured sections in carbonate rocks, thermal reservoirs are classified into three types: Type I, Type II, and Type III. Type I reservoirs have an acoustic transit time greater than 200 μs / m, porosity greater than 9.5%, and permeability greater than 15 mD. Type II reservoirs have an acoustic transit time greater than 185 μs / m, porosity between 5.0% and 9.5%, and permeability between 4.0 mD and 15 mD. Type III reservoirs have an acoustic transit time less than 185 μs / m, porosity between 1.8% and 5%, and permeability between 0.1 mD and 4.0 mD. The depth ranges where Class I and Class II reservoirs are located are considered the geological sweet sections of bedrock thermal reservoirs; The overlapping section of the geothermal sweet section and the geological sweet section is the bedrock thermal reservoir development sweet section.

2. The method for predicting sweet sections of bedrock geothermal reservoirs in uplift zones according to claim 1, wherein, Based on the geological model of the bedrock reservoir in the uplift zone, a numerical model of the bedrock reservoir is established, including: The geological model is then meshed. The thermal conductivity range of different lithologies in the geological profile was determined by field outcrops or core sampling. Thermal conductivity values ​​are assigned to different lithological combinations in the gridded formation to obtain a numerical model of the bedrock thermal reservoir.

3. The method for predicting sweet sections of bedrock geothermal reservoirs in uplift zones according to claim 1, wherein, The numerical model simulates the formation temperature distribution field of the bedrock uplift zone, including: Define the boundary conditions of the model and determine the numerical simulation formula; For the numerical model, the implicit finite difference method is used for iterative solution, and the formation temperature distribution at different times is obtained by forward modeling. Determine whether the consistency rate between the output formation temperature lateral variation curve and the formation temperature variation model at the top of the bedrock reservoir reaches a set threshold. If yes, output the formation temperature distribution field. If no, adjust the thermal conductivity parameters of the caprock and bedrock until the consistency rate reaches the set threshold.

4. The method for predicting sweet sections of bedrock geothermal reservoirs in uplift zones according to claim 3, wherein, The numerical simulation formula is as follows: （1） in, For the density of the formation, The specific heat capacity of the underground medium. It is an ancient geothermal area. , They are respectively , Number of meshed cells in the direction, , These represent the horizontal and vertical coordinate axes, respectively. This represents the number of time steps in the simulation.

5. A device for predicting sweet sections of bedrock thermal reservoirs in protruding zones, characterized in that, include: The model building module establishes a geological model of the bedrock reservoir in the uplift zone and a model of the temperature variation of the strata above the bedrock reservoir. The numerical model building module establishes a numerical model of the bedrock thermal reservoir based on the geological model of the uplift zone bedrock thermal reservoir. The simulation module simulates the formation temperature distribution field of the bedrock uplift zone based on the numerical model and the formation temperature change model of the bedrock reservoir top surface. The sweet spot determination module identifies the geothermal sweet spots and geological sweet spots of the bedrock thermal reservoir, and then determines the development sweet spots of the bedrock thermal reservoir; Among them, the geothermal sweet sections that have been identified as bedrock thermal reservoirs include: Determine the depth of the maximum positive geothermal anomaly in the uplift zone; Based on the formation temperature distribution field, geothermal wells in the uplift zone and the depression zone are selected respectively, and formation temperature-depth curves are plotted on the same coordinate axis. The depth of the intersection of the two curves is the heat conduction equilibrium depth line of the uplift zone and the depression zone. The intersection point is located at the depth of the bedrock thermal reservoir in the uplift zone, which is the high conductivity homogenization depth. The reservoir segment located between the maximum positive geothermal anomaly depth and the high conductivity homogenization depth is the geothermal sweet segment. Among them, the geological sweet sections that have been identified as bedrock geothermal reservoirs include: Identify anhydrous argillaceous carbonate rock sections and water-rich carbonate rock fracture sections; Based on the porosity and permeability conditions of fractured sections in carbonate rocks, thermal reservoirs are classified into three types: Type I, Type II, and Type III. Type I reservoirs have an acoustic transit time greater than 200 μs / m, porosity greater than 9.5%, and permeability greater than 15 mD. Type II reservoirs have an acoustic transit time greater than 185 μs / m, porosity between 5.0% and 9.5%, and permeability between 4.0 mD and 15 mD. Type III reservoirs have an acoustic transit time less than 185 μs / m, porosity between 1.8% and 5%, and permeability between 0.1 mD and 4.0 mD. The depth ranges where Class I and Class II reservoirs are located are considered the geological sweet sections of bedrock thermal reservoirs; The overlapping section of the geothermal sweet section and the geological sweet section is the bedrock thermal reservoir development sweet section.

6. An electronic device, characterized in that, The electronic device includes: Memory, which stores executable instructions; A processor that executes the executable instructions in the memory to implement the method for predicting sweet segments of bedrock thermal reservoirs in uplift zones as described in any one of claims 1-4.

7. 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 for predicting sweet sections of bedrock thermal reservoirs in any one of claims 1-4.

Images