Temperature prediction method and system for geothermal fluid at outlet of geothermal field

Abstract

The invention provides a temperature prediction method and system for geothermal fluid at an outlet of a geothermal field, and belongs to the technical field of geothermal development. The method for predicting the temperature of the geothermal fluid at the outlet of the geothermal field comprises the following steps: dividing a development target interval to obtain a plurality of thermal storage fluid layers of the development target interval and interpretation data of each thermal storage fluid layer; according to the interpretation data of each thermal storage fluid layer in the development target interval, determining the heat obtained after geothermal fluid in each thermal storage fluid layer reaches a wellhead from the current thermal storage fluid layer in unit time; according to the heat which can be obtained after the geothermal fluid in each thermal storage fluid layer reaches the wellhead from the current thermal storage fluid layer in unit time, the heat which can be obtained after the geothermal fluid in the development target interval reaches the wellhead is determined; and according to the heat obtained after the geothermal fluid in the development target interval reaches the wellhead, the temperature of the geothermal fluid in the development target interval after the geothermal fluid reaches the wellhead is estimated.

Description

Technical Field

[0001] This invention relates to the field of geothermal development technology, specifically to a method for predicting the temperature of geothermal fluid at the outlet of a geothermal field, a system for predicting the temperature of geothermal fluid at the outlet of a geothermal field, a computer device, and a computer-readable storage medium. Background Technology

[0002] Currently, the method used in the geothermal field to estimate the wellhead water temperature of production wells is as follows: first, the mid-depth of the geothermal reservoir is given, and the reservoir temperature is calculated using the geothermal gradient formula. Then, the temperature loss of the geothermal fluid from the reservoir to the wellhead is subtracted. This method often results in a significant deviation between the estimated wellhead water temperature and the actual wellhead temperature measured during pilot tests, reducing the practicality and reliability of the approach.

[0003] Due to the strong heterogeneity within and between layers of the geothermal reservoir and the large differences in physical properties, the geothermal reservoir is not utilized evenly. The extracted hot water mainly comes from the main producing layer, but the middle depth of the main producing layer deviates from the middle depth of the total geothermal reservoir, resulting in temperature deviations. Summary of the Invention

[0004] To address the aforementioned technical deficiencies, this invention provides a method and system for predicting the temperature of geothermal fluids at the outlet of a geothermal field. The method involves stratifying the geothermal reservoir within the target development zone to obtain multiple reservoir fluid layers with different interpretation data. The method calculates the heat that the geothermal fluid in each of the multiple reservoir fluid layers can obtain upon reaching the wellhead per unit time, and then summarizes the calculations to obtain the total heat that the geothermal fluid in the target development zone can obtain upon reaching the wellhead. The temperature of the geothermal fluid in the target development zone upon reaching the wellhead is then calculated based on the total heat that the geothermal fluid in the target development zone can obtain upon reaching the wellhead, thus ensuring the accuracy of the calculation.

[0005] The first aspect of this invention provides a method for predicting the temperature of geothermal fluid at the outlet of a geothermal field, comprising:

[0006] Based on the layer data of the development target layer of the geothermal field reservoir area and the logging interpretation data of the completed wells in the geothermal field reservoir area, the development target layer of the geothermal field reservoir area is divided into multiple reservoir fluid layers and interpretation data of each reservoir fluid layer in the geothermal field reservoir area.

[0007] Based on the interpretation data of each geothermal reservoir, determine the amount of heat that geothermal fluids in each geothermal reservoir can obtain per unit time after reaching the wellhead from the current geothermal reservoir.

[0008] Based on the amount of heat that geothermal fluid in each geothermal reservoir can obtain after reaching the wellhead from the current geothermal reservoir per unit time, determine the amount of heat that geothermal fluid in the development target section can obtain after reaching the wellhead.

[0009] The temperature of the geothermal fluid at the wellhead is estimated based on the heat that can be obtained after the geothermal fluid in the target development zone reaches the wellhead.

[0010] In this embodiment of the invention, the step of dividing the geothermal field reservoir into development target zones based on the zone data and well logging interpretation data of completed wells in the geothermal field reservoir, to obtain multiple reservoir fluid layers within the development target zones and interpretation data for each of the multiple reservoir fluid layers, includes:

[0011] Based on the top and bottom depths of the development target section and the logging interpretation data of the completed wells in the geothermal reservoir area, the development target section is divided to obtain multiple thermal reservoir fluid layers in the development target section, as well as the top depth, bottom depth, effective thickness, and permeability of each thermal reservoir fluid layer.

[0012] In this embodiment of the invention, determining the amount of heat that geothermal fluid in each geothermal reservoir can obtain per unit time after reaching the wellhead from the current geothermal reservoir based on the interpreted data of each geothermal reservoir includes:

[0013] Calculate the exploitable heat of geothermal fluids in each geothermal reservoir per unit time based on the interpreted data of each reservoir.

[0014] Calculate the heat loss of geothermal fluid in each geothermal reservoir per unit time when it is extracted from the current geothermal reservoir to the wellhead based on the interpreted data of each geothermal reservoir.

[0015] The amount of heat that the geothermal fluid in each reservoir can obtain per unit time after reaching the wellhead is calculated based on the recoverable heat of the geothermal fluid per unit time in each reservoir and the heat loss when it is extracted from the current reservoir to the wellhead.

[0016] In this embodiment of the invention, calculating the exploitable heat of the geothermal fluid in each geothermal reservoir per unit time based on the interpreted data of each geothermal reservoir includes:

[0017] The distributable exploitable amount of geothermal fluid in each geothermal reservoir per unit time is calculated based on the interpreted data of each geothermal reservoir.

[0018] Calculate the temperature of the geothermal fluid in each geothermal reservoir layer based on the interpreted data of each geothermal reservoir layer;

[0019] The exploitable heat of the geothermal fluid in each reservoir per unit time is calculated based on the distributable exploitable amount of geothermal fluid in each reservoir per unit time and the temperature of the geothermal fluid in each reservoir per unit time.

[0020] In this embodiment of the invention, the calculation of the exploitable heat of the geothermal fluid in each geothermal reservoir per unit time based on the distributable exploitable amount of geothermal fluid in each reservoir and the temperature of the geothermal fluid in each reservoir specifically involves:

[0021] Q i =Q wi c w ρ w (T i -T0);

[0022] Among them, Q wi Let c be the distributable exploitable amount of geothermal fluid in the i-th geothermal reservoir layer per unit time. w ρ is the specific heat capacity of the geothermal fluid. w T represents the density of the geothermal fluid. i Ti represents the temperature of the geothermal fluid in the i-th thermal reservoir layer, and T0 represents the temperature of the isothermal layer.

[0023] In this embodiment of the invention, calculating the distributable exploitable amount of geothermal fluid in each geothermal reservoir per unit time based on the interpreted data of each geothermal reservoir includes:

[0024] Calculate the distribution coefficient of the recoverable geothermal fluid in each geothermal reservoir based on the interpreted data of each geothermal reservoir;

[0025] The allocatable amount of geothermal fluid in each geothermal reservoir layer per unit time is calculated based on the allocation coefficient of the exploitable amount of geothermal fluid in each geothermal reservoir layer.

[0026] In this embodiment of the invention, the calculation of the distributable extractable amount of geothermal fluid in each geothermal reservoir per unit time based on the distribution coefficient of the extractable amount of geothermal fluid in each geothermal reservoir layer specifically involves:

[0027] Q wi =n×g i ;

[0028] Among them, Q wi Let g be the distributable exploitable amount of geothermal fluid in the i-th geothermal reservoir layer per unit time. i is the distribution coefficient of the exploitable amount of geothermal fluid in the i-th geothermal reservoir layer, and n is the amount of geothermal fluid extracted per unit time.

[0029] In this embodiment of the invention, the step of calculating the distribution coefficient of the recoverable geothermal fluid in each geothermal reservoir based on the interpreted data of each geothermal reservoir includes:

[0030] The distribution coefficient of the exploitable geothermal fluid in each geothermal reservoir is calculated based on the permeability and effective thickness of each geothermal reservoir.

[0031] In this embodiment of the invention, the calculation of the distribution coefficient of the recoverable geothermal fluid in each geothermal reservoir layer based on the permeability and effective thickness of each geothermal reservoir layer is specifically as follows:

[0032]

[0033] Among them, g i K is the distribution coefficient of the recoverable geothermal fluid in the i-th geothermal reservoir layer. i h is the permeability of the i-th thermal reservoir fluid layer. i denoted as the effective thickness of the i-th thermal storage fluid layer, and n as the number of thermal storage fluid layers within the target development segment.

[0034] In this embodiment of the invention, calculating the temperature of the geothermal fluid in each geothermal reservoir layer based on the interpreted data of each geothermal reservoir layer includes:

[0035] Calculate the middle depth of each thermal reservoir fluid layer based on the top and bottom depths of each layer.

[0036] The temperature of the geothermal fluid in each geothermal reservoir layer is calculated based on the mid-depth of each reservoir layer, the geothermal gradient, and the depth of the isothermal layer.

[0037] In this embodiment of the invention, the calculation of the temperature of the geothermal fluid in each geothermal reservoir layer based on the midpoint depth, geothermal gradient, and isothermal layer depth specifically involves:

[0038]

[0039] Among them, T i h represents the temperature of the geothermal fluid in the i-th thermal reservoir layer. i denoted as the middle depth of the i-th thermal reservoir fluid layer, h0 as the depth of the isothermal layer, ΔT as the geothermal gradient, and T0 as the temperature of the isothermal layer.

[0040] In this embodiment of the invention, the calculation of the heat that the geothermal fluid in each geothermal reservoir can obtain per unit time after reaching the wellhead from the current geothermal reservoir, based on the exploitable heat of the geothermal fluid in each geothermal reservoir per unit time and the heat loss when it is extracted from the current geothermal reservoir to the wellhead, specifically involves:

[0041] Q si =Q i -Q hli ;

[0042] Among them, Q si Let Q be the heat generated by the geothermal fluid in the i-th thermal reservoir after it reaches the wellhead. i Let Q be the exploitable heat of the geothermal fluid in the i-th geothermal reservoir layer per unit time. hli This represents the heat loss of geothermal fluid in the i-th geothermal reservoir layer when it is extracted from the current geothermal reservoir layer to the wellhead.

[0043] In this embodiment of the invention, calculating the heat loss of geothermal fluid in each geothermal reservoir per unit time when it is extracted from the current geothermal reservoir to the wellhead based on the interpreted data of each geothermal reservoir includes:

[0044] Calculate the middle depth of each thermal reservoir fluid layer based on the top and bottom depths of each layer.

[0045] The heat loss of geothermal fluid in each geothermal reservoir per unit time is calculated based on the heat loss per unit length of wellbore within the middle depth and formation depth range of each geothermal reservoir.

[0046] In this embodiment of the invention, the calculation of the heat loss of geothermal fluid in each geothermal reservoir per unit time when it is extracted from the current geothermal reservoir to the wellhead, based on the heat loss per unit length of wellbore within the middle depth and formation depth range of each geothermal reservoir, specifically involves:

[0047] Q hli =q l h i ;

[0048] Wherein, Q hli h represents the heat loss of geothermal fluid in the i-th geothermal reservoir layer per unit time when it is extracted from the i-th geothermal reservoir layer to the wellhead. i Let q be the depth of the middle part of the i-th thermal reservoir fluid layer. l This represents the heat loss per unit length of wellbore within the formation depth range.

[0049] In this embodiment of the invention, determining the heat that the geothermal fluid in the target development section can obtain upon reaching the wellhead based on the heat that the geothermal fluid in each geothermal reservoir can obtain per unit time after reaching the wellhead from the current geothermal reservoir includes:

[0050] To obtain the amount of heat that geothermal fluid in each reservoir can obtain per unit time after it reaches the wellhead from the current reservoir.

[0051] The heat that the geothermal fluid can obtain after reaching the wellhead from the current geothermal fluid layer in each of the multiple geothermal fluid layers is summed per unit time to obtain the heat that the geothermal fluid in the development target section can obtain after reaching the wellhead.

[0052] In this embodiment of the invention, estimating the temperature of the geothermal fluid at the wellhead based on the heat that can be obtained after the geothermal fluid in the target development zone reaches the wellhead includes:

[0053] To obtain the isothermal layer temperature, specific heat capacity of geothermal fluids, and density of geothermal fluids from well logging interpretation data of geothermal reservoir areas;

[0054] The temperature of the geothermal fluid in the target development zone after reaching the wellhead is calculated based on the heat that the geothermal fluid can obtain after reaching the wellhead, the temperature of the isothermal layer, the specific heat capacity of the geothermal fluid, and the density of the geothermal fluid.

[0055] In this embodiment of the invention, the calculation of the temperature of the geothermal fluid in the target development zone after reaching the wellhead based on the heat available to the geothermal fluid in the target development zone, the temperature of the isothermal layer, the specific heat capacity of the geothermal fluid, and the density of the geothermal fluid specifically involves:

[0056]

[0057] Where T1 is the temperature of the geothermal fluid in the target development zone after it reaches the wellhead, and Q... t The heat that can be obtained by geothermal fluids in the target formation after reaching the wellhead, the geothermal fluid production rate per unit time, the geothermal reservoir, and c. w ρ is the specific heat capacity of the geothermal fluid. w T0 is the density of the geothermal fluid and T0 is the temperature of the isothermal layer.

[0058] A second aspect of the present invention provides a temperature prediction system for geothermal fluid at the outlet of a geothermal field, comprising:

[0059] The stratigraphic division module is used to divide the geothermal field reservoir area into development target stratigraphic segments based on the stratigraphic segment data of the development target stratigraphic segment and the well logging interpretation data of the completed wells in the geothermal field reservoir area, thereby obtaining multiple reservoir fluid layers of the development target stratigraphic segment of the geothermal field reservoir area and interpretation data of each reservoir fluid layer.

[0060] The first prediction module is used to predict, based on the interpretation data of each geothermal reservoir in the development target section, the amount of heat that geothermal fluids in each geothermal reservoir can obtain per unit time after reaching the wellhead from the current geothermal reservoir.

[0061] The second prediction module is used to predict the heat that geothermal fluid in the development target section can obtain after reaching the wellhead based on the heat that geothermal fluid in each geothermal reservoir can obtain after reaching the wellhead from the current geothermal reservoir per unit time.

[0062] The third prediction module is used to estimate the temperature of the geothermal fluid in the target development zone after it reaches the wellhead, based on the heat that can be obtained after the geothermal fluid in the target development zone reaches the wellhead.

[0063] A third aspect of the present invention provides a computer device, comprising:

[0064] Memory, which stores computer programs;

[0065] A processor is used to execute the computer program to implement the method for predicting the temperature of geothermal fluids at the outlet of a geothermal field as described above.

[0066] A fourth aspect of the present invention provides a computer-readable storage medium having a computer program stored thereon, the computer program being executed by a processor to implement the method for predicting the temperature of geothermal fluid at the outlet of a geothermal field as described above.

[0067] The method for predicting the temperature of geothermal fluids at the outlet of a geothermal field involves stratifying the geothermal reservoir within the target development zone to obtain multiple reservoir fluid layers with different interpretation data. The method then calculates the heat that the geothermal fluids in each of these layers can obtain per unit time upon reaching the wellhead. These calculations are then combined to obtain the total heat that the geothermal fluids in the target development zone can obtain upon reaching the wellhead. The temperature of the geothermal fluids in the target development zone upon reaching the wellhead is calculated based on this total heat, ensuring the accuracy of the calculation.

[0068] Other features and advantages of the technical solution of the present invention will be described in detail in the following detailed embodiments section. Attached Figure Description

[0069] The accompanying drawings, which are included to provide a further understanding of the invention and form part of this invention, illustrate exemplary embodiments of the invention and are used to explain the invention, but do not constitute an undue limitation of the invention. In the drawings:

[0070] Figure 1 This is a flowchart of a method for predicting the temperature of geothermal fluid at the outlet of a geothermal field, provided by an embodiment of the present invention.

[0071] Figure 2 This is a schematic diagram of the well location deployment in the H well area provided in an embodiment of the present invention;

[0072] Figure 3This is a structural block diagram of a temperature prediction system for geothermal fluid at the outlet of a geothermal field, provided in an embodiment of the present invention. Detailed Implementation

[0073] To make the technical solutions and advantages of the embodiments of the present invention clearer, the exemplary embodiments of the present invention will be further described in detail below with reference to the accompanying drawings. Obviously, the described embodiments are only a part of the embodiments of the present invention, and not an exhaustive list of all embodiments. It should be noted that, unless otherwise specified, the embodiments and features in the embodiments of the present invention can be combined with each other.

[0074] In the description of this invention, it should be understood that the terms "center," "longitudinal," "lateral," "length," "width," "thickness," "upper," "lower," "front," "rear," "left," "right," "vertical," "horizontal," "top," "bottom," "inner," and "outer," etc., indicate the orientation or positional relationship based on the orientation or positional relationship shown in the accompanying drawings. They are only for the convenience of describing this invention and simplifying the description, and do not indicate or imply that the device or element referred to must have a specific orientation, or be constructed and operated in a specific orientation. Therefore, they should not be construed as limitations on this invention.

[0075] Furthermore, the terms "first" and "second" are used for descriptive purposes only and should not be construed as indicating or implying relative importance or implicitly specifying the number of technical features indicated. Thus, a feature defined as "first" or "second" may explicitly or implicitly include one or more of that feature. In the description of this invention, "a plurality of" means at least two, such as two, three, etc., unless otherwise explicitly specified.

[0076] In this invention, unless otherwise explicitly specified and limited, terms such as "installation," "connection," "linking," and "fixing" should be interpreted broadly. For example, they can refer to a fixed connection, a detachable connection, or an integral part; they can refer to a mechanical connection, an electrical connection, or a connection that allows communication between them; they can refer to a direct connection or an indirect connection through an intermediate medium; they can refer to the internal communication of two components or the interaction between two components. Those skilled in the art can understand the specific meaning of the above terms in this invention according to the specific circumstances.

[0077] In the process of realizing this invention, the inventors discovered that if the water temperature of the extraction well can be determined relatively accurately during the geothermal scheme preparation stage, and combined with the designed water intake, the heat that the scheme can provide can be estimated.

[0078] Currently, the method used in the geothermal field to estimate the wellhead water temperature of production wells is as follows: first, the mid-depth of the geothermal reservoir is given, and the reservoir temperature is calculated using the geothermal gradient formula. Then, the temperature loss of the geothermal fluid from the reservoir to the wellhead is subtracted. This method often results in a significant deviation between the estimated wellhead water temperature and the actual wellhead temperature measured during pilot tests, reducing the practicality and reliability of the approach.

[0079] Due to the strong heterogeneity within and between layers of the geothermal reservoir and the large differences in physical properties, the geothermal reservoir is not utilized evenly. The extracted hot water mainly comes from the main producing layer, but the middle depth of the main producing layer deviates from the middle depth of the total geothermal reservoir, resulting in temperature deviations.

[0080] To address the aforementioned problems, this invention provides a method for predicting the temperature of geothermal fluids at the outlet of a geothermal field. The method includes: dividing the geothermal field reservoir into development target zones based on zone data and logging interpretation data from completed wells in the geothermal field reservoir, resulting in multiple reservoir fluid layers within the development target zones and interpretation data for each of these layers; determining the amount of heat that geothermal fluids in each reservoir fluid layer can obtain per unit time after reaching the wellhead from the current reservoir fluid layer based on the interpretation data of each layer; determining the amount of heat that geothermal fluids in the development target zones can obtain upon reaching the wellhead based on the amount of heat that geothermal fluids in each reservoir fluid layer can obtain per unit time after reaching the wellhead from the current reservoir fluid layer; and estimating the temperature of the geothermal fluids in the development target zones upon reaching the wellhead based on the amount of heat that geothermal fluids in the development target zones can obtain upon reaching the wellhead. The method for predicting the temperature of geothermal fluids at the outlet of a geothermal field involves stratifying the geothermal reservoir within the target development zone to obtain multiple reservoir fluid layers with different interpretation data. The method then calculates the heat that the geothermal fluids in each of these layers can obtain per unit time upon reaching the wellhead. These calculations are then combined to obtain the total heat that the geothermal fluids in the target development zone can obtain upon reaching the wellhead. The temperature of the geothermal fluids in the target development zone upon reaching the wellhead is calculated based on this total heat, ensuring the accuracy of the calculation.

[0081] Figure 1 This is a flowchart illustrating a method for predicting the temperature of geothermal fluid at the outlet of a geothermal field, provided by an embodiment of the present invention. Figure 1 As shown in the figure, this embodiment provides a method for predicting the temperature of geothermal fluid at the outlet of a geothermal field, which includes the following steps:

[0082] S1. Based on the layer data of the development target layer of the geothermal field reservoir area and the logging interpretation data of the completed wells in the geothermal field reservoir area, the development target layer of the geothermal field reservoir area is divided to obtain multiple reservoir fluid layers of the development target layer of the geothermal field reservoir area and interpretation data of each reservoir fluid layer.

[0083] S2. Determine the amount of heat that geothermal fluids in each geothermal reservoir can obtain per unit time after reaching the wellhead from the current geothermal reservoir, based on the interpretation data of each geothermal reservoir layer;

[0084] S3. Determine the amount of heat that the geothermal fluid in the target development section can obtain after reaching the wellhead within a unit time, based on the amount of heat that the geothermal fluid in each geothermal reservoir can obtain after reaching the wellhead from the current geothermal reservoir.

[0085] S4. Calculate the temperature of the geothermal fluid in the target development zone after it reaches the wellhead based on the heat that can be obtained after the geothermal fluid in the target development zone reaches the wellhead.

[0086] In step S1, the process of dividing the geothermal field reservoir area into development target zones based on the zone data and well logging interpretation data from completed wells in the geothermal field reservoir area, to obtain multiple reservoir fluid layers within the development target zones and interpretation data for each of these reservoir fluid layers, includes:

[0087] Based on the top and bottom depths of the development target section and the logging interpretation data of the completed wells in the geothermal reservoir area, the development target section is divided to obtain multiple geothermal fluid layers of the development target section, as well as the top depth, bottom depth, effective thickness, and permeability of each geothermal fluid layer.

[0088] Specifically, S11. Determine the development target layer of the geothermal reservoir area and obtain the layer data of the development target layer;

[0089] S12. Identify completed wells within the geothermal reservoir area that meet preset standards, and obtain logging interpretation data of the completed wells;

[0090] S13. The development target section is divided according to the section data and the logging interpretation data of the completed wells to obtain multiple thermal reservoir fluid layers and interpretation data of each thermal reservoir fluid layer.

[0091] Specifically, the top and bottom depths of the target development zone are matched with the depth data in the logging interpretation data of the completed wells. The depth region that is the same as the depth region of the target development zone is determined in the logging interpretation data, and the thermal reservoir fluid layer in the depth region of the logging interpretation data is determined.

[0092] Furthermore, the layer data of the development target layer includes: the top depth of the layer and the bottom depth of the layer;

[0093] The logging interpretation data of the completed wells includes: interpretation data of multiple reservoirs, isothermal layer depth, isothermal layer temperature, geothermal gradient, and heat loss per unit length of wellbore within the formation depth range;

[0094] The interpreted data for each thermal storage fluid layer includes: the top depth of each thermal storage fluid layer, the bottom depth of each thermal storage fluid layer, the effective thickness of each thermal storage fluid layer, and the permeability of each thermal storage fluid layer.

[0095] In step S2, calculating the heat that geothermal fluid in each geothermal reservoir can obtain per unit time after reaching the wellhead from the current geothermal reservoir based on the interpreted data of each geothermal reservoir includes:

[0096] S21. Calculate the exploitable heat of geothermal fluids in each geothermal reservoir per unit time based on the interpreted data of each geothermal reservoir;

[0097] S22. Calculate the heat loss of geothermal fluid in each geothermal reservoir layer per unit time when it is extracted from the current geothermal reservoir layer to the wellhead, based on the interpreted data of each geothermal reservoir layer;

[0098] S23. Calculate the amount of heat that the geothermal fluid in each geothermal reservoir can obtain per unit time after reaching the wellhead from the current geothermal reservoir, based on the recoverable heat of the geothermal fluid in each geothermal reservoir per unit time and the heat loss when it is extracted from the current geothermal reservoir to the wellhead.

[0099] In step S23, the formula for calculating the amount of heat that geothermal fluid in each thermal reservoir can obtain per unit time after reaching the wellhead from the current thermal reservoir is as follows:

[0100] Q si =Q i -Q hli ;

[0101] Among them, Q si Let Q be the heat generated by the geothermal fluid in the i-th thermal reservoir after it reaches the wellhead. i Let Q be the exploitable heat of the geothermal fluid in the i-th geothermal reservoir layer per unit time. hliThis represents the heat loss of geothermal fluid in the i-th geothermal reservoir layer when it is extracted from the current geothermal reservoir layer to the wellhead.

[0102] In step S21, calculating the exploitable heat of the geothermal fluid in each geothermal reservoir per unit time based on the interpreted data of each reservoir includes:

[0103] S211. Calculate the distributable exploitable amount of geothermal fluid in each geothermal reservoir per unit time based on the interpreted data of each geothermal reservoir.

[0104] S212. Calculate the temperature of the geothermal fluid in each geothermal reservoir layer based on the interpreted data of each geothermal reservoir layer;

[0105] S213. Calculate the recoverable heat of the geothermal fluid in each geothermal reservoir per unit time based on the distributable exploitable amount of geothermal fluid in each geothermal reservoir per unit time and the temperature of the geothermal fluid in each geothermal reservoir per unit time.

[0106] In step S213, the formula for calculating the exploitable heat of the geothermal fluid in each geothermal reservoir per unit time is as follows:

[0107] Q i =Q wi c w ρ w (T i -T0);

[0108] Among them, Q wi The unit is the distributable exploitable amount of geothermal fluid in the i-th geothermal reservoir layer per unit time, expressed in m³. 3 / h,c w ρ is the specific heat capacity of the geothermal fluid, expressed in kJ / (kg·K). w The density of the geothermal fluid is expressed in kg / m³. 3 T i Ti represents the temperature of the geothermal fluid in the i-th thermal reservoir layer, in K, and T0 represents the temperature of the isothermal layer, in K.

[0109] In step S211, calculating the distributable exploitable amount of geothermal fluid in each geothermal reservoir per unit time based on the interpreted data of each geothermal reservoir includes:

[0110] S2111. Calculate the distribution coefficient of the recoverable geothermal fluid in each geothermal reservoir based on the interpreted data of each geothermal reservoir;

[0111] S2112. Calculate the distributable amount of geothermal fluid in each geothermal reservoir per unit time based on the distribution coefficient of the exploitable amount of geothermal fluid in each geothermal reservoir.

[0112] Furthermore, the formula for calculating the distributable exploitable amount of geothermal fluid in each geothermal reservoir per unit time is as follows:

[0113] Q wi =n×g i ;

[0114] Among them, Q wi Let g be the distributable exploitable amount of geothermal fluid in the i-th geothermal reservoir layer per unit time. i Let n be the distribution coefficient of the recoverable geothermal fluid in the i-th geothermal reservoir layer, and n be the geothermal fluid recovery rate per unit time, in m³. 3 / h.

[0115] In step S2111, the calculation of the distribution coefficient of the recoverable geothermal fluid in each geothermal reservoir based on the interpretation data of each geothermal reservoir includes:

[0116] The distribution coefficient of the exploitable geothermal fluid in each geothermal reservoir is calculated based on the permeability and effective thickness of each geothermal reservoir.

[0117] In this embodiment of the invention, the formula for calculating the distribution coefficient of the recoverable geothermal fluid in each geothermal reservoir layer is as follows:

[0118]

[0119] Among them, g i K is the distribution coefficient of the recoverable geothermal fluid in the i-th geothermal reservoir layer. i Let be the permeability of the i-th thermal reservoir fluid layer, in mD, h. i denoted as the effective thickness of the i-th thermal storage fluid layer, in meters (m), and n as the number of thermal storage fluid layers within the target development segment.

[0120] In step S212, calculating the temperature of the geothermal fluid in each geothermal reservoir layer based on the interpreted data of each reservoir layer includes:

[0121] S2121. Calculate the middle depth of each thermal storage fluid layer based on the top and bottom depths of each thermal storage fluid layer;

[0122] S2122. Calculate the temperature of the geothermal fluid in each geothermal reservoir layer based on the mid-depth of each geothermal reservoir layer, the geothermal gradient, and the depth of the isothermal layer.

[0123] Furthermore, the formula for calculating the temperature of the geothermal fluid in each geothermal reservoir layer is as follows:

[0124]

[0125] Among them, T i The temperature of the geothermal fluid in the i-th thermal reservoir layer is expressed in K and h. i denoted as the midpoint depth of the i-th thermal reservoir fluid layer in meters, h0 as the depth of the isothermal layer in meters, ΔT as the geothermal gradient in ℃ / 100m, and T0 as the temperature of the isothermal layer in K.

[0126] In step S22, calculating the heat loss of geothermal fluid in each geothermal reservoir per unit time when it is extracted from the current geothermal reservoir to the wellhead based on the interpreted data of each geothermal reservoir includes:

[0127] S221. Calculate the middle depth of each thermal fluid reservoir layer based on the top and bottom depths of each thermal fluid reservoir layer;

[0128] S222. Calculate the heat loss of geothermal fluid in each geothermal reservoir per unit time when it is extracted from the current geothermal reservoir to the wellhead, based on the heat loss per unit length of the wellbore within the middle depth and formation depth range of each geothermal reservoir.

[0129] In step S222, the formula for calculating the heat loss of geothermal fluid in each geothermal reservoir layer when it is extracted from the current geothermal reservoir layer to the wellhead per unit time is as follows:

[0130] Q hli =q l h i ;

[0131] Wherein, Q hli h represents the heat loss of geothermal fluid in the i-th geothermal reservoir layer per unit time when it is extracted from the i-th geothermal reservoir layer to the wellhead. i Let q be the depth of the middle part of the i-th thermal reservoir fluid layer. l This represents the heat loss per unit length of wellbore within the formation depth range.

[0132] In step S3, calculating the heat that the geothermal fluid in the target development zone can obtain upon reaching the wellhead based on the heat that the geothermal fluid in each geothermal reservoir can obtain per unit time after reaching the wellhead from the current geothermal reservoir includes:

[0133] S31. Obtain the amount of heat that geothermal fluid in each thermal reservoir can obtain after reaching the wellhead from the current thermal reservoir per unit time;

[0134] S32. Sum the heat that the geothermal fluid in each of the multiple geothermal reservoir layers can obtain per unit time after reaching the wellhead from the current geothermal reservoir layer to obtain the heat that the geothermal fluid in the development target section can obtain after reaching the wellhead.

[0135] Furthermore, the formula for calculating the heat that can be obtained by the geothermal fluid in the target development zone after reaching the wellhead is as follows:

[0136]

[0137] Among them, Q t Q is the amount of heat that can be obtained by the geothermal fluid in the target formation after it reaches the wellhead. si Let be the amount of heat that geothermal fluid in the i-th geothermal reservoir can obtain per unit time after reaching the wellhead from the current geothermal reservoir, and n be the number of geothermal reservoirs in the development target section.

[0138] In step S4, calculating the temperature of the geothermal fluid at the wellhead based on the heat that can be obtained from the geothermal fluid in the target development zone after reaching the wellhead includes:

[0139] S41. Obtain the isothermal layer temperature, specific heat capacity of geothermal fluids, and density of geothermal fluids from the well logging interpretation data of the geothermal reservoir area;

[0140] S42. Calculate the temperature of the geothermal fluid in the target development zone after it reaches the wellhead based on the heat that can be obtained after the geothermal fluid reaches the wellhead, the temperature of the isothermal layer, the specific heat capacity of the geothermal fluid, and the density of the geothermal fluid.

[0141] Furthermore, the formula for calculating the temperature of the geothermal fluid in the target development zone after reaching the wellhead is as follows:

[0142] Where T1 is the temperature of the geothermal fluid in the target development zone after it reaches the wellhead, and Q... t The heat that can be obtained by the geothermal fluid in the target formation after reaching the wellhead is denoted by n, which is the amount of geothermal fluid extracted per unit time, in m³. 3 / h thermal storage fluid layer, c w ρ is the specific heat capacity of the geothermal fluid. w T0 is the density of the geothermal fluid and T0 is the temperature of the isothermal layer.

[0143] This embodiment uses permeability and effective thickness to determine the recoverable amount of geothermal fluid in each reservoir layer, ensuring the accuracy of the recoverable amount for each reservoir layer. By weighting the recoverable amount of geothermal fluid in each reservoir layer, the accurate calculation of the recoverable amount of geothermal fluid in the target development section is guaranteed. This addresses the problem of deviation between the actual reservoir temperature and the predicted reservoir temperature. This invention considers the influence of reservoir heterogeneity. The uneven utilization of reservoirs may lead to large errors in predicted temperature during the planned development stage. For the first time, this invention proposes a method using permeability and effective thickness weighting to calculate the recoverable amount of geothermal fluid per unit time in each small reservoir layer, thereby estimating the wellhead fluid temperature. This method considers more factors than the original method and is closer to reality.

[0144] This embodiment uses a specific mining area as an example, as follows:

[0145] The M mining area has a heating area of ​​760,000 square meters and a maximum heating load of 30.6 MW. It is planned to construct a 17.5 MW geothermal heating station within the M mining area's boiler room to replace the basic heating load of the M mining area. The H well area, located near the M mining area's coal-fired boiler room, is preferred as the geothermal reservoir. Structurally, the H well area is located in the northeastern part of a rift basin, with the main geothermal reservoir being the S3 section of the Shahejie Formation, belonging to a sedimentary basin-type geothermal reservoir. The target layer for this development is the S3 section. 2 The sub-section has a porosity of 18.9% and a permeability of 531.1 mD. The reservoir is at a depth of 1801 m, with single-layer thicknesses ranging from 3 to 19 m and a total thickness of 175 to 325 m, averaging 210 m. It belongs to a medium-porosity to high-permeability layered geothermal reservoir. The geothermal gradient is 3.0℃ / 100 m. Based on the principles of "production determined by injection, full utilization of old wells, and reinjection within the same layer," in the H well area, S3... 2 For the target layer, a reverse seven-point well pattern was adopted, with a well spacing of 560m along the long axis and 400m along the short axis of the sand body. A layout of 3 extraction wells and 5 irrigation wells was deployed, including 3 older wells with a single well extraction capacity of 100m³. 3 / h, single well reinjection volume 45-70m³ 3 / h, irrigation scale 300m 3 / h. For example... Figure 2 As shown, Figure 2 This is a schematic diagram of the well location deployment in the H well area provided in an embodiment of the present invention.

[0146] To predict the water intake temperature of Well 3-1, we can obtain the well logging interpretation data of the geothermal reservoir area based on the logging data of the completed well A, and then use this data to predict the water intake temperature of Well 3-1. The target interval of the dominant geothermal reservoir in Well A (i.e., the development target interval) is 1736.3-1865.8m, and a total of 7 water layers were interpreted by well logging, which are 1, 2, 3, 4, 5, 6, and 7 from top to bottom.

[0147] In a certain well area, the geothermal reservoir has a central depth of 1801m. The target development layers are sub-layers 1-7, with effective thicknesses of 1.6m, 6.6m, 2.5m, 4.2m, 2.6m, 4.4m, and 6.1m, respectively, and permeabilities of 432mD, 782mD, 916mD, 677mD, 646mD, 1004mD, and 925mD, respectively. The top depths are 1736.3m, 1775.7m, 1789.2m, and 1801m, respectively. The wellbore depths are 0.7m, 1820.4m, 1827.3m, and 1859.7m, with bottom depths of 1737.9m, 1782.3m, 1791.7m, 1811.2m, 1823m, 1831.7m, and 1865.8m, and middle depths of 1737.1m, 1779m, 1790.45m, 1809.1m, 1821.7m, 1829.5m, and 1862.75m, respectively. The isothermal layer has a depth of 30m and a temperature of 10.5℃. The geothermal gradient is 3℃ / 100m. The specific heat capacity of water is 4.180J / (Kg·K). The heat loss per unit length of wellbore within the depth range of 0-1800m is 0.28KJ / (m·h).

[0148] 1. Calculate the distribution coefficient g of the exploitable geothermal fluid in layers 1-7. i

[0149]

[0150]

[0151] 2. Calculate the distributable geothermal fluid extraction rate Q per unit time for layers 1-7. wi

[0152] Q w1 =n×g1=100×0.03=3m 3 / h

[0153] Q w2 =n×g2=100×0.23=23m 3 / h

[0154] Q w3 =n×g3=100×0.1=10m 3 / h

[0155] Q w4 =n×g4=100×0.13=13m 3 / h

[0156] Q w5 =n×g5=100×0.07=7m 3 / h

[0157] Q w6=n×g6=100×0.19=19m 3 / h

[0158] Q w7 =n×g7=100×0.25=25m 3 / h

[0159] 3. Calculate the thermal storage temperature T of layers 1-7. i

[0160]

[0161]

[0162] 4. Calculate the exploitable heat Q of geothermal fluid per unit time for layers 1-7. i

[0163] Q1 = Q w1 c w ρ w (T1-T0)=3×0.00418×1000×(61.713+273.15-10.5-273.15)=642.2KJ / h

[0164] Q2 = Q w2 c w ρ w (T2-T0)=23×0.00418×1.0×(62.97+273.15-10.5-273.15)=5044.5KJ / h

[0165] Q3 = Q w3 c w ρ w (T3-T0)=10×0.00418×1.0×(63.31+273.15-10.5-273.15)=2207.6KJ / h

[0166] Q4 = Q w4 c w ρ w (T4-T0)=13×0.00418×1.0×(63.873+273.15-10.5-273.15)=2900.3KJ / h

[0167] Q5 = Q w5 c w ρ w (T5-T0)=7×0.00418×1.0×(64.251+273.15-10.5-273.15)=1572.8KJ / h

[0168] Q6 = Q w6 cw ρ w (T6-T0)=19×0.00418×1.0×(64.485+273.15-10.5-273.15)=4287.5KJ / h

[0169] Q7 = Q w7 c w ρ w (T7-T0)=25×0.00418×1.0×(65.48+273.15-10.5-273.15)=5745.7KJ / h

[0170] 5. Calculate the heat loss Q per unit time of geothermal fluid extraction to the wellhead in layers 1-7. hli

[0171] Q hl1 =q l h1 = 0.28 × 1737.1 = 486.4 kJ / h

[0172] Q hl2 =q l h2 = 0.28 × 1779 = 498.1 kJ / h

[0173] Q hl3 =q l h3=0.28×1790.45=501.3KJ / h

[0174] Q hl4 =q l h4 = 0.28 × 1809.1 = 506.5 kJ / h

[0175] Q hl5 =q l h5 = 0.28 × 1821.7 = 510.1 kJ / h

[0176] Q hl6 =q l h6 = 0.28 × 1829.5 = 512.3 kJ / h

[0177] Q hl7 =q l h7=0.28×1862.75=521.6KJ / h

[0178] 6. Calculate the heat Q that can be extracted per unit time after the geothermal fluid reaches the wellhead in layers 1-7. si

[0179] Q s1 =Q1-Q hl1 =642.2-486.4=155.8KJ / h

[0180] Q s2 =Q2-Q hl2 =5044.5-498.1=4546.4KJ / h

[0181] Q s3 =Q3-Q hl3 =2207.6-501.3=1706.3KJ / h

[0182] Q s4 =Q4-Q hl4 =2900.3-506.5=2393.8KJ / h

[0183] Q s5 =Q5-Q hl5 =1572.8-510.1=1062.7KJ / h

[0184] Q s6 =Q6-Q hl6 =4287.5-512.3=3775.2KJ / h

[0185] Q s7 =Q7-Q hl7 =5745.7-521.6=5224.1KJ / h

[0186] 7. Calculate the total heat Q extracted from the geothermal fluid at the wellhead per unit time. t

[0187]

[0188] 8. Calculate the temperature T1 of the geothermal fluid extracted from the wellhead.

[0189]

[0190] Figure 3 This is a structural block diagram of a geothermal field outlet geothermal fluid temperature prediction system provided in an embodiment of the present invention. Figure 3 As shown in the figure, this embodiment provides a temperature prediction system for geothermal fluid at the outlet of a geothermal field, comprising:

[0191] The stratigraphic division module is used to divide the geothermal field reservoir area into development target stratigraphic segments based on the stratigraphic segment data of the development target stratigraphic segment and the well logging interpretation data of the completed wells in the geothermal field reservoir area, thereby obtaining multiple reservoir fluid layers in the development target stratigraphic segment of the geothermal field reservoir area and interpretation data of each reservoir fluid layer.

[0192] The first prediction module is used to predict, based on the interpretation data of each geothermal reservoir, the amount of heat that geothermal fluids in each geothermal reservoir can obtain per unit time after reaching the wellhead from the current geothermal reservoir.

[0193] The second prediction module is used to predict the heat that the geothermal fluid in the development target section can obtain after reaching the wellhead based on the heat that the geothermal fluid in each thermal reservoir can obtain after reaching the wellhead from the current thermal reservoir per unit time.

[0194] The third prediction module is used to estimate the temperature of the geothermal fluid after it reaches the wellhead based on the heat that the geothermal fluid in the target development zone can obtain after reaching the wellhead.

[0195] Specifically, the stratigraphic division module is used to determine the development target stratigraphic segment of the geothermal reservoir area and obtain the stratigraphic segment data of the development target stratigraphic segment; determine the completed wells in the geothermal reservoir area that meet the preset standards and obtain the logging interpretation data of the completed wells; divide the development target stratigraphic segment according to the stratigraphic segment data of the development target stratigraphic segment and the logging interpretation data of the completed wells to obtain multiple geothermal reservoir fluid layers and interpretation data of each geothermal reservoir fluid layer.

[0196] Specifically, the top and bottom depths of the target development zone are matched with the depth data in the logging interpretation data of the completed wells. The depth region that is the same as the depth region of the target development zone is determined in the logging interpretation data, and the thermal reservoir fluid layer in the depth region of the logging interpretation data is determined.

[0197] Furthermore, the layer data for the target development layer includes: the top depth and bottom depth of the layer; the logging interpretation data from the completed wells includes: interpretation data for multiple reservoirs, isothermal layer depth, isothermal layer temperature, geothermal gradient, and heat loss per unit length of wellbore within the formation depth range; the interpretation data for each thermal reservoir fluid layer includes: the top depth of each thermal reservoir fluid layer, the bottom depth of each thermal reservoir fluid layer, the effective thickness of each thermal reservoir fluid layer, and the permeability of each thermal reservoir fluid layer.

[0198] The first prediction module is specifically used for:

[0199] Calculate the exploitable heat of geothermal fluids in each geothermal reservoir per unit time based on the interpreted data of each reservoir.

[0200] Calculate the heat loss of geothermal fluid in each geothermal reservoir per unit time when it is extracted from the current geothermal reservoir to the wellhead based on the interpreted data of each geothermal reservoir.

[0201] The amount of heat that the geothermal fluid in each reservoir can obtain per unit time after reaching the wellhead is calculated based on the recoverable heat of the geothermal fluid per unit time in each reservoir and the heat loss when it is extracted from the current reservoir to the wellhead.

[0202] The formula for calculating the amount of heat that geothermal fluid can obtain per unit time after reaching the wellhead from the current geothermal reservoir is as follows:

[0203] Q si =Q i -Q hli ;

[0204] Among them, Q si Let Q be the heat generated by the geothermal fluid in the i-th thermal reservoir after it reaches the wellhead. i Let Q be the exploitable heat of the geothermal fluid in the i-th geothermal reservoir layer per unit time. hli This represents the heat loss of geothermal fluid in the i-th geothermal reservoir layer when it is extracted from the current geothermal reservoir layer to the wellhead.

[0205] The calculation of the exploitable heat of geothermal fluids in each geothermal reservoir per unit time based on the interpreted data of each reservoir includes:

[0206] The distributable exploitable amount of geothermal fluid in each geothermal reservoir per unit time is calculated based on the interpreted data of each geothermal reservoir.

[0207] Calculate the temperature of the geothermal fluid in each geothermal reservoir layer based on the interpreted data of each geothermal reservoir layer;

[0208] The exploitable heat of the geothermal fluid in each reservoir per unit time is calculated based on the distributable exploitable amount of geothermal fluid in each reservoir per unit time and the temperature of the geothermal fluid in each reservoir per unit time.

[0209] The formula for calculating the exploitable heat of geothermal fluid in each geothermal reservoir per unit time is as follows:

[0210] Q i =Q wi c w ρ w (T i -T0);

[0211] Among them, Q wi The unit is the distributable exploitable amount of geothermal fluid in the i-th geothermal reservoir layer per unit time, expressed in m³. 3 / h,c w ρ is the specific heat capacity of the geothermal fluid, expressed in kJ / (kg·K). w The density of the geothermal fluid is expressed in kg / m³. 3 T i Ti represents the temperature of the geothermal fluid in the i-th thermal reservoir layer, in K, and T0 represents the temperature of the isothermal layer, in K.

[0212] The calculation of the distributable exploitable amount of geothermal fluid in each geothermal reservoir per unit time based on the interpreted data of each geothermal reservoir includes:

[0213] Calculate the distribution coefficient of the recoverable geothermal fluid in each geothermal reservoir based on the interpreted data of each geothermal reservoir;

[0214] The allocatable amount of geothermal fluid in each geothermal reservoir layer per unit time is calculated based on the allocation coefficient of the exploitable amount of geothermal fluid in each geothermal reservoir layer.

[0215] Furthermore, the formula for calculating the distributable exploitable amount of geothermal fluid in each geothermal reservoir per unit time is as follows:

[0216] Q wi =n×g i ;

[0217] Among them, Q wi Let g be the distributable exploitable amount of geothermal fluid in the i-th geothermal reservoir layer per unit time. i Let n be the distribution coefficient of the recoverable geothermal fluid in the i-th geothermal reservoir layer, and n be the geothermal fluid recovery rate per unit time, in m³. 3 / h.

[0218] The calculation of the distribution coefficient of the recoverable geothermal fluid in each geothermal reservoir based on the interpreted data of each geothermal reservoir includes:

[0219] The distribution coefficient of the exploitable geothermal fluid in each geothermal reservoir is calculated based on the permeability and effective thickness of each geothermal reservoir.

[0220] In this embodiment, the formula for calculating the distribution coefficient of recoverable geothermal fluid in each geothermal reservoir layer is as follows:

[0221]

[0222] Among them, g i K is the distribution coefficient of the recoverable geothermal fluid in the i-th geothermal reservoir layer. i Let be the permeability of the i-th thermal reservoir fluid layer, in mD, h. i The effective thickness of the i-th thermal storage fluid layer is expressed in meters (m).

[0223] The calculation of the temperature of the geothermal fluid in each geothermal reservoir layer based on the interpreted data of each geothermal reservoir layer includes:

[0224] Calculate the middle depth of each thermal reservoir fluid layer based on the top and bottom depths of each layer.

[0225] The temperature of the geothermal fluid in each geothermal reservoir layer is calculated based on the mid-depth of each reservoir layer, the geothermal gradient, and the depth of the isothermal layer.

[0226] Furthermore, the formula for calculating the temperature of the geothermal fluid in each geothermal reservoir layer is as follows:

[0227]

[0228] Among them, T i The temperature of the geothermal fluid in the i-th thermal reservoir layer is expressed in K and h. i denoted as the midpoint depth of the i-th thermal reservoir fluid layer in meters, h0 as the depth of the isothermal layer in meters, ΔT as the geothermal gradient in ℃ / 100m, and T0 as the temperature of the isothermal layer in K.

[0229] The calculation of the heat loss of geothermal fluid in each geothermal reservoir per unit time when it is extracted from the current geothermal reservoir to the wellhead, based on the interpreted data of each geothermal reservoir, includes:

[0230] Calculate the middle depth of each thermal reservoir fluid layer based on the top and bottom depths of each layer.

[0231] The heat loss of geothermal fluid in each geothermal reservoir per unit time is calculated based on the heat loss per unit length of wellbore within the middle depth and formation depth range of each geothermal reservoir.

[0232] The formula for calculating the heat loss of geothermal fluid in each reservoir layer per unit time when it is extracted from the current reservoir layer to the wellhead is as follows:

[0233] Q hli =q l h i ;

[0234] Wherein, Q hli h represents the heat loss of geothermal fluid in the i-th geothermal reservoir layer per unit time when it is extracted from the i-th geothermal reservoir layer to the wellhead. i Let q be the depth of the middle part of the i-th thermal reservoir fluid layer. l This represents the heat loss per unit length of wellbore within the formation depth range.

[0235] The second prediction module is specifically used for:

[0236] To obtain the amount of heat that geothermal fluid in each reservoir can obtain per unit time after it reaches the wellhead from the current reservoir.

[0237] The heat that geothermal fluid can obtain from the current geothermal reservoir layer to the wellhead per unit time is summed to obtain the heat that geothermal fluid in the development target section can obtain after reaching the wellhead.

[0238] Furthermore, the formula for calculating the heat that can be obtained by the geothermal fluid in the target development zone after reaching the wellhead is as follows:

[0239]

[0240] Among them, Q t Q is the amount of heat that can be obtained by the geothermal fluid in the target formation after it reaches the wellhead. si Let be the amount of heat that geothermal fluid in the i-th geothermal reservoir can obtain per unit time after reaching the wellhead from the current geothermal reservoir, and n be the number of geothermal reservoirs in the development target section.

[0241] The third prediction module is specifically used for:

[0242] Obtain the isothermal layer temperature, specific heat capacity of geothermal fluids, and density of geothermal fluids from the well logging interpretation data of the geothermal reservoir area;

[0243] The temperature of the geothermal fluid at the wellhead in the target development zone is calculated based on the heat that the geothermal fluid can obtain after reaching the wellhead, the temperature of the isothermal layer, the specific heat capacity of the geothermal fluid, and the density of the geothermal fluid.

[0244] Furthermore, the formula for calculating the temperature of the geothermal fluid in the target development zone after reaching the wellhead is as follows:

[0245] Where T1 is the temperature of the geothermal fluid in the target development zone after it reaches the wellhead, and Q... t The amount of heat that geothermal fluids in the target formation can obtain after reaching the wellhead, expressed in m³ / s. 3 / h thermal storage fluid layer, c w ρ is the specific heat capacity of the geothermal fluid. w T0 is the density of the geothermal fluid and T0 is the temperature of the isothermal layer.

[0246] The present invention also provides a computer device, including: a memory, a processor, and a computer program, the computer program being stored in the memory and configured to be executed by the processor to implement the above-described method for predicting the temperature of geothermal fluid at the outlet of a geothermal field.

[0247] The present invention also provides a machine-readable storage medium storing computer program instructions thereon, which, when executed by a processor, implement the above-described method for predicting the temperature of geothermal fluid at the outlet of a geothermal field.

[0248] Those skilled in the art will understand that embodiments of the present invention can be provided as methods, systems, or computer program products. Therefore, the present invention can take the form of a completely hardware embodiment, a completely software embodiment, or an embodiment combining software and hardware aspects. Furthermore, the present invention can take the form of a computer program product implemented on one or more computer-usable storage media (including but not limited to disk storage, CD-ROM, optical storage, etc.) containing computer-usable program code. The solutions in the embodiments of the present invention can be implemented using various computer languages, such as the object-oriented programming language Java and the interpreted scripting language JavaScript.

[0249] This invention is described with reference to flowchart illustrations and / or block diagrams of methods, apparatus (systems), and computer program products according to embodiments of the invention. It will be understood that each block of the flowchart illustrations and / or block diagrams, and combinations of blocks in the flowchart illustrations and / or block diagrams, can be implemented by computer program instructions. These computer program instructions can be provided to a processor of a general-purpose computer, special-purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, generate instructions for implementing the flowchart illustrations and / or block diagrams. Figure 1 One or more processes and / or boxes Figure 1 A device that provides the functions specified in one or more boxes.

[0250] These computer program instructions may also be stored in a computer-readable storage medium that can direct a computer or other programmable data processing device to function in a particular manner, such that the instructions stored in the computer-readable storage medium produce an article of manufacture including instruction means, which are implemented in a process Figure 1 One or more processes and / or boxes Figure 1 The function specified in one or more boxes.

[0251] These computer program instructions may also be loaded onto a computer or other programmable data processing equipment to cause a series of operational steps to be performed on the computer or other programmable equipment to produce a computer-implemented process, thereby providing instructions that execute on the computer or other programmable equipment for implementing the process. Figure 1 One or more processes and / or boxes Figure 1 The steps of the function specified in one or more boxes.

[0252] Although preferred embodiments of the invention have been described, those skilled in the art, upon learning the basic inventive concept, can make other changes and modifications to these embodiments. Therefore, the appended claims are intended to be interpreted as including the preferred embodiments as well as all changes and modifications falling within the scope of the invention.

[0253] Obviously, those skilled in the art can make various modifications and variations to this invention without departing from its spirit and scope. Therefore, if these modifications and variations fall within the scope of the claims of this invention and their equivalents, this invention also intends to include these modifications and variations.

Claims

1. A method for predicting the temperature of geothermal fluid at the outlet of a geothermal field, characterized in that, include: Based on the layer data of the development target layer of the geothermal field reservoir area and the logging interpretation data of the completed wells in the geothermal field reservoir area, the development target layer of the geothermal field reservoir area is divided into multiple reservoir fluid layers and interpretation data of each reservoir fluid layer in the geothermal field reservoir area. Based on the interpretation data of each geothermal reservoir, determine the amount of heat that geothermal fluids in each geothermal reservoir can obtain per unit time after reaching the wellhead from the current geothermal reservoir. Based on the amount of heat that geothermal fluid in each geothermal reservoir can obtain after reaching the wellhead from the current geothermal reservoir per unit time, determine the amount of heat that geothermal fluid in the development target section can obtain after reaching the wellhead. The temperature of the geothermal fluid at the wellhead is estimated based on the heat that can be obtained after the geothermal fluid in the target development zone reaches the wellhead.

2. The method for predicting the temperature of geothermal fluid at the outlet of a geothermal field according to claim 1, characterized in that, The process involves dividing the geothermal field reservoir into development target zones based on the zone data and well logging interpretation data from completed wells. This results in multiple reservoir fluid layers within the development target zones, as well as interpretation data for each of these layers, including: Based on the top and bottom depths of the development target section and the logging interpretation data of the completed wells in the geothermal reservoir area, the development target section is divided to obtain multiple thermal reservoir fluid layers in the development target section, as well as the top depth, bottom depth, effective thickness, and permeability of each thermal reservoir fluid layer.

3. The method for predicting the temperature of geothermal fluid at the outlet of a geothermal field according to claim 2, characterized in that, The step of determining the amount of heat that geothermal fluids in each geothermal reservoir can obtain per unit time after reaching the wellhead from the current geothermal reservoir, based on the interpreted data of each geothermal reservoir, includes: Calculate the exploitable heat of geothermal fluids in each geothermal reservoir per unit time based on the interpreted data of each reservoir. Calculate the heat loss of geothermal fluid in each geothermal reservoir per unit time when it is extracted from the current geothermal reservoir to the wellhead based on the interpreted data of each geothermal reservoir. The amount of heat that the geothermal fluid in each reservoir can obtain per unit time after reaching the wellhead is calculated based on the recoverable heat of the geothermal fluid per unit time in each reservoir and the heat loss when it is extracted from the current reservoir to the wellhead.

4. The method for predicting the temperature of geothermal fluid at the outlet of a geothermal field according to claim 3, characterized in that, The calculation of the exploitable heat of geothermal fluids in each geothermal reservoir per unit time based on the interpreted data of each reservoir includes: The distributable exploitable amount of geothermal fluid in each geothermal reservoir per unit time is calculated based on the interpreted data of each geothermal reservoir. Calculate the temperature of the geothermal fluid in each geothermal reservoir layer based on the interpreted data of each geothermal reservoir layer; The exploitable heat of the geothermal fluid in each reservoir per unit time is calculated based on the distributable exploitable amount of geothermal fluid in each reservoir per unit time and the temperature of the geothermal fluid in each reservoir per unit time.

5. The method for predicting the temperature of geothermal fluid at the outlet of a geothermal field according to claim 4, characterized in that, The calculation of the exploitable heat capacity of the geothermal fluid in each geothermal reservoir per unit time, based on the distributable exploitable amount of geothermal fluid in each reservoir and the temperature of the geothermal fluid in each reservoir, specifically involves: Q i =Q wi c w ρ w (T i -T0); Among them, Q wi Let c be the distributable exploitable amount of geothermal fluid in the i-th geothermal reservoir layer per unit time. w ρ is the specific heat capacity of the geothermal fluid. w T represents the density of the geothermal fluid. i Ti represents the temperature of the geothermal fluid in the i-th thermal reservoir layer, and T0 represents the temperature of the isothermal layer.

6. The method for predicting the temperature of geothermal fluid at the outlet of a geothermal field according to claim 4, characterized in that, The calculation of the distributable exploitable amount of geothermal fluid in each geothermal reservoir per unit time based on the interpreted data of each geothermal reservoir includes: Calculate the distribution coefficient of the recoverable geothermal fluid in each geothermal reservoir based on the interpreted data of each geothermal reservoir; The allocatable amount of geothermal fluid in each geothermal reservoir layer per unit time is calculated based on the allocation coefficient of the exploitable amount of geothermal fluid in each geothermal reservoir layer.

7. The method for predicting the temperature of geothermal fluid at the outlet of a geothermal field according to claim 6, characterized in that, The calculation of the distributable extractable amount of geothermal fluid in each geothermal reservoir per unit time based on the distribution coefficient of the extractable amount of geothermal fluid in each reservoir is specifically as follows: Q wi =n×g i ; Among them, Q wi Let g be the distributable exploitable amount of geothermal fluid in the i-th geothermal reservoir layer per unit time. i is the distribution coefficient of the exploitable amount of geothermal fluid in the i-th geothermal reservoir layer, and n is the amount of geothermal fluid extracted per unit time.

8. The method for predicting the temperature of geothermal fluid at the outlet of a geothermal field according to claim 6, characterized in that, The calculation of the distribution coefficient of the recoverable geothermal fluid in each geothermal reservoir based on the interpreted data of each geothermal reservoir includes: The distribution coefficient of the exploitable geothermal fluid in each geothermal reservoir is calculated based on the permeability and effective thickness of each geothermal reservoir.

9. The method for predicting the temperature of geothermal fluid at the outlet of a geothermal field according to claim 8, characterized in that, The allocation coefficient for the recoverable geothermal fluid in each geothermal reservoir layer is calculated based on the permeability and effective thickness of each geothermal reservoir layer. Specifically: Among them, g i K is the distribution coefficient of the recoverable geothermal fluid in the i-th geothermal reservoir layer. i h is the permeability of the i-th thermal reservoir fluid layer. i denoted as the effective thickness of the i-th thermal storage fluid layer, and n as the number of thermal storage fluid layers within the target development segment.

10. The method for predicting the temperature of geothermal fluid at the outlet of a geothermal field according to claim 4, characterized in that, The calculation of the temperature of the geothermal fluid in each geothermal reservoir layer based on the interpreted data of each geothermal reservoir layer includes: Calculate the middle depth of each thermal reservoir fluid layer based on the top and bottom depths of each layer. The temperature of the geothermal fluid in each geothermal reservoir layer is calculated based on the mid-depth of each reservoir layer, the geothermal gradient, and the depth of the isothermal layer.

11. The method for predicting the temperature of geothermal fluid at the outlet of a geothermal field according to claim 10, characterized in that, The calculation of the geothermal fluid temperature in each geothermal reservoir layer based on the mid-depth, geothermal gradient, and isothermal layer depth is specifically as follows: Among them, T i h represents the temperature of the geothermal fluid in the i-th thermal reservoir layer. i denoted as the middle depth of the i-th thermal reservoir fluid layer, h0 as the depth of the isothermal layer, ΔT as the geothermal gradient, and T0 as the temperature of the isothermal layer.

12. The method for predicting the temperature of geothermal fluid at the outlet of a geothermal field according to claim 3, characterized in that, The calculation of the amount of heat that the geothermal fluid in each reservoir can obtain per unit time after reaching the wellhead from the current reservoir, based on the recoverable heat of the geothermal fluid in each reservoir per unit time and the heat loss during extraction from the current reservoir to the wellhead, specifically involves: Q si =Q i -Q hli ; Among them, Q si Let Q be the heat generated by the geothermal fluid in the i-th thermal reservoir after it reaches the wellhead. i Let Q be the exploitable heat of the geothermal fluid in the i-th geothermal reservoir layer per unit time. hli This represents the heat loss of geothermal fluid in the i-th geothermal reservoir layer when it is extracted from the current geothermal reservoir layer to the wellhead.

13. The method for predicting the temperature of geothermal fluid at the outlet of a geothermal field according to claim 3, characterized in that, The calculation of the heat loss of geothermal fluid in each geothermal reservoir per unit time when it is extracted from the current geothermal reservoir to the wellhead, based on the interpreted data of each geothermal reservoir, includes: Calculate the middle depth of each thermal reservoir fluid layer based on the top and bottom depths of each layer. The heat loss of geothermal fluid in each geothermal reservoir per unit time is calculated based on the heat loss per unit length of wellbore within the middle depth and formation depth range of each geothermal reservoir.

14. The method for predicting the temperature of geothermal fluid at the outlet of a geothermal field according to claim 13, characterized in that, The calculation of the heat loss of geothermal fluid in each geothermal reservoir per unit time when it is extracted from the current geothermal reservoir to the wellhead is based on the heat loss per unit length of wellbore within the middle depth and formation depth range of each geothermal reservoir. Specifically: Q hli =q l h i ; Wherein, Q hli h represents the heat loss of geothermal fluid in the i-th geothermal reservoir layer per unit time when it is extracted from the i-th geothermal reservoir layer to the wellhead. i Let q be the depth of the middle part of the i-th thermal reservoir fluid layer. l This represents the heat loss per unit length of wellbore within the formation depth range.

15. The method for predicting the temperature of geothermal fluid at the outlet of a geothermal field according to claim 2, characterized in that, The step of determining the heat that the geothermal fluid in the target development zone can obtain upon reaching the wellhead based on the heat that the geothermal fluid in each reservoir can obtain per unit time after reaching the wellhead from the current reservoir includes: To obtain the amount of heat that geothermal fluid in each reservoir can obtain per unit time after it reaches the wellhead from the current reservoir. The heat that the geothermal fluid can obtain after reaching the wellhead from the current geothermal fluid layer in each of the multiple geothermal fluid layers is summed per unit time to obtain the heat that the geothermal fluid in the development target section can obtain after reaching the wellhead.

16. The method for predicting the temperature of geothermal fluid at the outlet of a geothermal field according to claim 2, characterized in that, The method of estimating the temperature of the geothermal fluid at the wellhead based on the heat that can be obtained after the geothermal fluid in the target development zone reaches the wellhead includes: To obtain the isothermal layer temperature, specific heat capacity of geothermal fluids, and density of geothermal fluids from well logging interpretation data of geothermal reservoir areas; The temperature of the geothermal fluid in the target development zone after reaching the wellhead is calculated based on the heat that the geothermal fluid can obtain after reaching the wellhead, the temperature of the isothermal layer, the specific heat capacity of the geothermal fluid, and the density of the geothermal fluid.

17. The method for predicting the temperature of geothermal fluid at the outlet of a geothermal field according to claim 16, characterized in that, The calculation of the temperature of the geothermal fluid at the wellhead in the target development zone based on the heat available to the geothermal fluid upon reaching the wellhead, the temperature of the isothermal layer, the specific heat capacity of the geothermal fluid, and the density of the geothermal fluid is specifically as follows: Where T1 is the temperature of the geothermal fluid in the target development zone after it reaches the wellhead, and Q... t The heat that can be obtained by the geothermal fluid in the target formation after reaching the wellhead, where n is the amount of geothermal fluid extracted per unit time, and c is the geothermal reservoir layer. w ρ is the specific heat capacity of the geothermal fluid. w T0 is the density of the geothermal fluid and T0 is the temperature of the isothermal layer.

18. A temperature prediction system for geothermal fluid at the outlet of a geothermal field, characterized in that, include: The stratigraphic division module is used to divide the geothermal field reservoir area into development target stratigraphic segments based on the stratigraphic segment data of the development target stratigraphic segment and the logging interpretation data of the completed wells in the geothermal field reservoir area, thereby obtaining multiple reservoir fluid layers of the development target stratigraphic segment of the geothermal field reservoir area and interpretation data of each reservoir fluid layer. The first prediction module is used to predict, based on the interpretation data of each geothermal reservoir, the amount of heat that geothermal fluids in each geothermal reservoir can obtain per unit time after reaching the wellhead from the current geothermal reservoir. The second prediction module is used to predict the heat that the geothermal fluid in the development target section can obtain after reaching the wellhead based on the heat that the geothermal fluid in each thermal reservoir can obtain after reaching the wellhead from the current thermal reservoir in a unit time. The third prediction module is used to estimate the temperature of the geothermal fluid after it reaches the wellhead based on the heat that the geothermal fluid in the target development zone can obtain after reaching the wellhead.

19. A computer device, characterized in that, include: Memory, which stores computer programs; A processor for executing the computer program to implement the method for predicting the temperature of geothermal fluid at the outlet of a geothermal field, as described in any one of claims 1 to 17.

20. A computer-readable storage medium having a computer program stored thereon, characterized in that, The computer program is executed by a processor to implement the method for predicting the temperature of geothermal fluid at the outlet of a geothermal field as described in any one of claims 1 to 17.

Images