A high-temperature geothermal multi-source measurement data fusion prediction method and system in a complex structure area

Abstract

The application discloses a high-temperature geothermal multi-source measurement data fusion prediction method and system in a complex structure area, and relates to the technical field of geothermal resource exploration. The method comprises the following steps: obtaining multi-source measurement data in the complex structure area and classifying the data into measurement units, and marking deep geophysical anomaly units, heat-controlling structure units and drilling temperature anomaly units; taking the deep geophysical anomaly units as the starting point, connecting the heat-controlling structure units to the drilling temperature anomaly units, and generating candidate heat anomaly transmission paths; screening heat-controlling dominant paths based on anomaly continuity, thermal resistance and non-response intervals; combining thermal physical properties, heat generation data and temperature measurement residuals to correct the temperature field, and generating target isotherm burial depth and high-temperature geothermal target area prediction results. Through the construction of measurement units, heat anomaly transmission paths and residual attribution correction mechanisms, the multi-source data fusion, heat-controlling path identification and high-temperature geothermal target area prediction accuracy are improved.

Description

Technical Field

[0001] This invention relates to the field of geothermal resource exploration technology, specifically to a method and system for fusion and prediction of high-temperature geothermal multi-source measurement data in complex tectonic zones. Background Technology

[0002] High-temperature geothermal resources are an important type of clean energy resource. In complex tectonic zones, the distribution of their thermal anomalies is usually controlled by a combination of factors, including deep heat sources, fault structures, fold structures, lithological assemblage, thermal conductivity of sedimentary caprock, and local heat transfer disturbances. Current high-temperature geothermal prediction typically relies on data from geological surveys, gravity measurements, magnetotelluric sounding, well temperature measurements, rock thermophysical property testing, and radioactive heat generation testing. Favorable areas are delineated through comprehensive artificial interpretation or temperature field simulation.

[0003] However, there are significant differences in the detection scale, response depth, and physical meaning of different types of measurement data in complex geological zones. Gravity measurements mainly reflect differences in subsurface density, magnetotelluric sounding mainly reflects changes in electrical structure, well temperature data can only reflect the temperature response near the well site, and rock thermophysical property data and radioactive heat generation data mostly come from discrete sample points. If the above data are directly overlaid for analysis, it is easy to encounter problems with unclear correspondences between deep geophysical anomalies, heat-controlling structures, and well temperature anomalies, leading to isolated anomalies being mistakenly identified as high-temperature geothermal target areas.

[0004] Furthermore, existing temperature field predictions are mostly based on heat conduction models, which input thermophysical parameters, heat generation parameters, and boundary conditions into the model to generate temperature field results. This approach typically fails to adequately represent the path relationships of deep geophysical anomalies transmitted to drilling temperature anomalies through heat-controlling units such as faults, fold cores, low thermal conductivity strata, or basalt layers. It also struggles to attribute the residual between the measured drilling temperature and the initial predicted temperature to specific heat-controlling path segments, thereby failing to correct the corresponding heat conduction calculation conditions. Therefore, in areas with insufficient deep drilling data, limited distribution of temperature measurement wells, and complex structural and lithological conditions, the predicted depth of the target isotherm and the high-temperature geothermal target area are easily affected by single anomalies or local measurement points, and there is still room for improvement in prediction stability and target area delineation reliability. Summary of the Invention

[0005] Based on the shortcomings of the prior art described above, the purpose of this invention is to provide a method and system for fusion and prediction of high-temperature geothermal multi-source measurement data in complex structural areas, so as to solve the above-mentioned technical problems.

[0006] To achieve the above objectives, the present invention provides the following technical solution: a method for fusion and prediction of high-temperature geothermal multi-source measurement data in complex tectonic zones, comprising:

[0007] Acquire multi-source measurement data in complex structural areas, and categorize the multi-source measurement data into measurement units according to planar location, depth range, and stratigraphic interface to generate a measurement unit set;

[0008] Anomaly markers are applied to the measurement unit set to generate deep geophysical anomaly units, thermal control structure units, and drilling temperature anomaly units;

[0009] Starting from the deep geophysical anomaly unit, candidate thermal anomaly transmission paths are generated by connecting the adjacent measurement units through the thermal control structure unit segment by segment to the drilling temperature anomaly unit.

[0010] The consistency of the thermal anomaly path is calculated based on the continuous distribution of anomaly units in the candidate thermal anomaly propagation path, the thermal resistance of adjacent measurement units, and the interval of non-responding measurement units. The thermal anomaly path is then selected and generated according to the consistency of the thermal anomaly path.

[0011] The initial temperature field prediction results are generated by performing heat conduction calculations based on the thermophysical property data and radioactive heat generation data in the measurement unit corresponding to the dominant heat control path.

[0012] The temperature measurement residual is generated by calculating the difference between the actual drilling temperature and the initial temperature field prediction results in the same well location and the same depth range. The temperature measurement residual is then segmented and aggregated according to the distance relationship between the measurement unit where the temperature measurement residual is located and different path segments in the heat control dominant path and the layer attribution relationship to generate the temperature measurement residual attribution correction amount.

[0013] Based on the temperature measurement residual attribution correction, the heat conduction calculation conditions of the corresponding path segment are corrected, and the heat conduction calculation is re-performed to generate the corrected temperature field prediction results.

[0014] Extract the minimum burial depth surface that reaches the target temperature from the corrected temperature field prediction results to generate the target isothermal surface burial depth prediction results. Merge the measurement units whose target isothermal surface burial depth meets the preset exploration depth conditions and are connected to the dominant heat control path to generate the high temperature geothermal target area prediction results.

[0015] The present invention is further configured such that, in acquiring multi-source measurement data of complex structural areas, the multi-source measurement data is categorized into measurement units according to planar location, depth range, and stratigraphic interface to generate a measurement unit set, which includes:

[0016] Multi-source measurement data for complex tectonic zones include geological structural data, gravity measurement data, magnetotelluric sounding data, well temperature measurement data, rock thermal property data, and radioactive heat generation data;

[0017] Stratigraphic constraint subdivision data is established based on stratigraphic interfaces, fault boundaries, and well-exposed strata in geological structural data; planar subdivision data is established based on the location of measuring points, measuring lines, and well locations.

[0018] The planar subdivision data and the stratigraphically constrained subdivision data are matched, and measurement units are generated according to the depth range defined by the planar subdivision range and the adjacent stratigraphic interfaces.

[0019] Gravity measurement data and magnetotelluric sounding data are categorized into corresponding measurement units according to measurement point location, survey line coverage, inversion depth, and detection coverage to generate geophysical data items; well temperature measurement data are categorized into corresponding measurement units according to well location and temperature measurement depth to generate temperature data items; rock thermal property data are categorized into corresponding measurement units according to sample collection location, sample collection depth, and the strata to which the sample belongs to generate thermal property data items; radioactive heat generation data are categorized into corresponding measurement units according to sample collection location, sample collection depth, and the strata to which the sample belongs to generate heat generation data items; geological structural data are categorized into corresponding measurement units according to stratigraphic interfaces, fault boundaries, lithological contact interfaces, and well-exposed strata to generate geological attribute data items.

[0020] A measurement unit set is generated based on the measurement units containing geophysical data items, temperature data items, thermal property data items, heat generation data items, and geological attribute data items.

[0021] The present invention is further configured such that the step of generating deep geophysical anomaly units, thermal control structure units, and drilling temperature anomaly units by anomaly marking the measurement unit set includes:

[0022] Extract exploration data items, geological attribute data items, and temperature data items from the measurement unit;

[0023] Geophysical data items are layered according to the detection depth range or inversion depth layer, and measurement units corresponding to changes in deep density or deep electrical interface are identified. The corresponding measurement units are marked as deep geophysical anomaly units.

[0024] The geological attribute data items are used to identify thermally controlling geological elements. Measurement units containing fault cuts, fold cores, basalt layers, low thermal conductivity strata, or lithological contact interfaces are marked as thermally controlling structural units.

[0025] The drilling temperature curve is segmented and identified for the temperature data items. The measurement units containing the corresponding high temperature segments or temperature gradient abrupt change segments are marked as drilling temperature anomaly units.

[0026] The present invention is further configured such that the process of generating a candidate thermal anomaly transmission path, starting from a deep geophysical anomaly unit and connecting segment by segment along the connection direction of adjacent measurement units through heat-controlling structural units to a drilling temperature anomaly unit, includes:

[0027] The deep geophysical exploration anomaly unit is used as the starting unit of the path, and the drilling temperature anomaly unit is extracted as the ending unit of the path.

[0028] Establish an adjacent unit connection table based on the shared plane boundary, shared depth boundary and formation contact interface between measurement units;

[0029] Starting from the path starting unit, the adjacent measurement units with heat control structure markers or drilling temperature anomaly markers are connected segment by segment along the adjacent unit connection table to generate a unit connection sequence;

[0030] When a unit connection sequence passes through at least one thermal control construction unit and connects to a path termination unit, the unit connection sequence is recorded as a candidate thermal anomaly propagation path.

[0031] The present invention is further configured such that the step of calculating the thermal anomaly path consistency degree based on the continuous distribution of anomaly units in the candidate thermal anomaly propagation path, the thermal resistance of adjacent measurement units, and the interval of non-responding measurement units, and selecting and generating the dominant thermal control path according to the thermal anomaly path consistency degree includes:

[0032] Generate a path unit sequence according to the path connection order in the candidate thermal anomaly propagation path;

[0033] Abnormal continuous segment data is generated by merging the abnormal markers of measurement units in the path unit sequence;

[0034] Path thermal resistance data is generated based on the thermal property data and depth interval of adjacent measurement units in the path unit sequence, and the interval data of non-response measurement units between adjacent abnormal continuous segments is identified.

[0035] The thermal anomaly path consistency is generated based on the abnormal continuous segment data, path thermal resistance data, and non-response measurement unit interval data. Candidate thermal anomaly propagation paths are then selected based on the thermal anomaly path consistency to generate the dominant thermal control path.

[0036] The present invention is further configured such that the initial temperature field prediction result is generated by performing heat conduction calculations based on the thermophysical property data and radioactive heat generation data in the measurement unit corresponding to the dominant heat control path, including:

[0037] Extract the measurement units along the dominant heat control path to generate a sequence of temperature field calculation units;

[0038] The thermal property data items, radioactive heat generation data items, stratigraphic attribution identifiers, and heat-controlling structure markers are read from the temperature field calculation unit sequence to generate path heat conduction constraint data;

[0039] Based on the path heat conduction constraint data, the fracture penetration unit is configured as a heat conduction channel constraint, the low thermal conductivity stratum unit is configured as a heat insulation constraint, the basalt layer unit is configured as a local heat transfer disturbance constraint, and the radioactive heat generation data item is configured as a layer internal heat generation constraint.

[0040] Heat conduction calculations are performed based on constraints such as heat conduction channels, insulation, local heat transfer disturbances, and internal heat generation within the layer to generate initial temperature field prediction results.

[0041] The present invention is further configured such that the calculation of the temperature measurement residual and the segmentation and aggregation of the temperature measurement residual to generate the temperature measurement residual attribution correction amount include:

[0042] Extract the well location predicted temperature data corresponding to the drilling temperature measurement data from the initial temperature field prediction results, and match the well location predicted temperature data with the actual drilling temperature according to the drilling number and depth range to generate well temperature corresponding data.

[0043] The difference between the measured drilling temperature and the predicted well location temperature for the same drilling number and the same depth range in the well temperature data is processed to generate a temperature measurement residual record.

[0044] The heat control dominant path is divided into a heat conduction section, a heat insulation section, and a heat transfer disturbance section according to the heat control structure markings;

[0045] Based on the spatial proximity and layer-level affiliation between the measurement unit where the temperature measurement residual record is located and the heat conduction section, heat preservation section, and heat transfer disturbance section, the temperature measurement residual record is segmented and aggregated to generate the temperature measurement residual attribution correction amount.

[0046] The present invention is further configured such that, in correcting the heat conduction calculation conditions of the corresponding path segment based on the temperature measurement residual attribution correction amount, and re-performing the heat conduction calculation to generate the corrected temperature field prediction result, the following steps are taken:

[0047] The temperature measurement residual attribution correction amount is associated with the heat conduction section, heat preservation section or heat transfer disturbance section in the heat control dominant path according to the path segment type, and the path segment correction data is generated.

[0048] Based on the path segment correction data, the equivalent thermal conductivity or boundary heat flow input of the heat conduction segment, the equivalent thermal conductivity or interlayer thermal resistance of the insulation segment, and the local heat transfer correction terms of the heat transfer disturbance segment are corrected to generate the corrected heat conduction calculation conditions.

[0049] Based on the corrected heat conduction calculation conditions, the heat conduction of the measurement unit set is recalculated to generate the corrected predicted temperature of each measurement unit;

[0050] The corrected temperature field prediction results are generated by assembling the corrected predicted temperatures according to the planar location, depth range, and formation interface of each measurement unit.

[0051] The present invention is further configured such that the step of extracting the minimum burial depth surface that reaches the target temperature from the corrected temperature field prediction results to generate the target isothermal surface burial depth prediction results, and merging the measurement units whose target isothermal surface burial depth meets the preset exploration depth conditions and is connected to the dominant heat control path to generate the high-temperature geothermal target area prediction results, includes:

[0052] The corrected predicted temperature, plane location identifier, depth interval identifier, and stratigraphic affiliation identifier of each measurement unit are extracted from the corrected temperature field prediction results to generate corrected temperature unit data;

[0053] The calibration temperature unit data is retrieved in order from shallow to deep within the same plane position, and the measurement unit that first reaches the target temperature is taken as the target temperature arrival unit.

[0054] Based on the depth range of the target temperature reaching the unit and the location of the formation interface, isothermal surface burial points are generated. The isothermal surface burial points corresponding to each plane location are spatially aggregated to generate the target isothermal surface burial depth prediction results.

[0055] The predicted depth of the target isotherm surface is matched with the preset exploration depth conditions, and the measurement units whose depth of the target isotherm surface is within the preset exploration depth range are selected to generate candidate target area units.

[0056] Connectivity matching is performed between candidate target units and measurement units traversed by the dominant heat control path, and candidate target units that are connected to the dominant heat control path in adjacent plane positions, adjacent depth intervals, or within the same formation contact interface are retained.

[0057] The retained candidate target area units are merged according to their adjacency and stratigraphic continuity to generate target area connectivity unit groups;

[0058] The prediction results for the high-temperature geothermal target area are generated based on the outer edge measurement units of the target area connectivity unit group, the burial depth of the corresponding target isothermal surface, and the associated dominant heat control path.

[0059] This invention also provides a system for fusing and predicting high-temperature geothermal multi-source measurement data in complex tectonic zones, the system comprising:

[0060] Measurement Unit Module: Acquires multi-source measurement data in complex structural areas, and categorizes the multi-source measurement data into measurement units according to planar location, depth range, and stratigraphic interface to generate a measurement unit set;

[0061] Anomaly Unit Marking Module: Marks anomalies in the measurement unit set to generate deep geophysical anomaly units, thermal control structure units, and drilling temperature anomaly units;

[0062] Path generation module: Starting from the deep geophysical anomaly unit, it connects segment by segment along the connection direction of the adjacent measurement unit through the heat-controlling structural unit to the drilling temperature anomaly unit to generate candidate thermal anomaly transmission paths.

[0063] Path filtering module: Calculates the consistency of thermal anomaly paths based on the continuous distribution of anomaly units in candidate thermal anomaly propagation paths, the thermal resistance of adjacent measurement units, and the interval between non-responding measurement units, and generates the dominant thermal control path according to the consistency of thermal anomaly paths.

[0064] Temperature field prediction module: Based on the thermophysical property data and radioactive heat generation data in the measurement unit corresponding to the dominant heat control path, heat conduction calculation is performed to generate initial temperature field prediction results;

[0065] Residual Attribution Module: Calculates the difference between the measured drilling temperature and the initial temperature field prediction results in the same well location and depth range to generate temperature measurement residuals; Based on the distance relationship and layer attribution relationship between the measurement unit where the temperature measurement residuals are located and different path segments in the dominant heat control path, the temperature measurement residuals are segmented and aggregated to generate temperature measurement residual attribution correction quantities.

[0066] Temperature field correction module: Corrects the heat conduction calculation conditions of the corresponding path segment based on the temperature measurement residual attribution correction amount, and re-performs heat conduction calculation to generate the corrected temperature field prediction results;

[0067] Target area generation module: Extracts the minimum burial depth surface that reaches the target temperature from the corrected temperature field prediction results to generate the target isothermal surface burial depth prediction results. Merges the measurement units whose target isothermal surface burial depth meets the preset exploration depth conditions and is connected to the heat control dominant path to generate high temperature geothermal target area prediction results.

[0068] This invention provides a method and system for fusing and predicting high-temperature geothermal multi-source measurement data in complex structural areas. The method acquires multi-source measurement data from the complex structural area and categorizes the data into measurement units according to planar location, depth range, and formation interface to generate a measurement unit set. Anomaly markers are applied to the measurement unit set to generate deep geophysical anomaly units, thermal control structural units, and drilling temperature anomaly units. Starting from the deep geophysical anomaly unit, candidate thermal anomaly transmission paths are generated by connecting adjacent measurement units along the connection direction passing through the thermal control structural unit to the drilling temperature anomaly unit. The consistency of the thermal anomaly path is calculated based on the continuous distribution of anomaly units in the candidate thermal anomaly transmission path, the thermal resistance of adjacent measurement units, and the interval between unresponsive measurement units. Thermal control dominant paths are then selected based on the consistency of the thermal anomaly path. Finally, thermal property data are generated based on the thermal property data of the measurement units corresponding to the thermal control dominant path. The initial temperature field prediction result is generated by performing heat conduction calculations based on radioactive heat generation data; the difference between the measured drilling temperature and the initial temperature field prediction result is calculated in the same well location and depth range to generate the temperature measurement residual; according to the distance relationship and stratigraphic affiliation relationship between the measurement unit where the temperature measurement residual is located and different path segments in the heat control dominant path, the temperature measurement residual is segmented and collected to generate the temperature measurement residual attribution correction amount; the heat conduction calculation conditions of the corresponding path segments are corrected according to the temperature measurement residual attribution correction amount, and heat conduction calculation is re-performed to generate the corrected temperature field prediction result; the minimum burial depth surface that reaches the target temperature is extracted from the corrected temperature field prediction result to generate the target isothermal surface burial depth prediction result; the measurement units whose target isothermal surface burial depth meets the preset exploration depth conditions and are connected to the heat control dominant path are merged to generate the high-temperature geothermal target area prediction result. The beneficial effects include:

[0069] 1. Improve the continuity of multi-source measurement data fusion: By unifying geophysical, temperature, thermal properties and heat generation data through measurement unit sets, discrete measurement point, profile data and drilling data are transformed into serializable unit data, reducing fusion deviations caused by inconsistencies in different data scales and depths;

[0070] 2. Improve the accuracy of thermal control path identification: By comparing the candidate thermal anomaly transmission path and the consistency of the thermal anomaly path, the continuous distribution of anomaly units, unit thermal resistance, and non-response interval are used to screen the dominant thermal control path, so as to avoid isolated geophysical anomalies or local temperature anomalies directly affecting the target area judgment.

[0071] 3. Improve the reliability of temperature field prediction and target area delineation: Correct the heat conduction calculation conditions of the corresponding path segment by adjusting the temperature measurement residual attribution correction amount, and then generate the corrected temperature field and the target isothermal surface burial depth, so that the prediction results of high temperature geothermal target area are simultaneously constrained by the measured temperature and the heat control path.

[0072] The above description is only an overview of the technical solution of this application. In order to better understand the technical means of this application and to implement it in accordance with the contents of the specification, and to make the above and other objects, features and advantages of this application more obvious and understandable, the following are specific embodiments of this application. Attached Figure Description

[0073] To more clearly illustrate the technical solutions in the embodiments of the present invention, the accompanying drawings used in the description of the embodiments will be briefly introduced below. Obviously, the drawings described below are only some embodiments of the present invention. For those skilled in the art, other drawings can be obtained based on these drawings without creative effort. In the drawings:

[0074] Figure 1 A flowchart illustrating a method for fusing and predicting high-temperature geothermal multi-source measurement data in complex tectonic zones, as shown in an exemplary embodiment of the present invention;

[0075] Figure 2 A schematic diagram of the structure of a high-temperature geothermal multi-source measurement data fusion and prediction system in a complex tectonic zone is shown as an exemplary embodiment of the present invention.

[0076] Figure 3 A temperature measurement curve of a temperature measuring orifice is shown as an exemplary embodiment of the present invention;

[0077] Figure 4 Explanation diagram of deep electrical structure shown in an exemplary embodiment of the present invention;

[0078] Figure 5 Temperature profiles and isotherm diagrams are shown as an exemplary embodiment of the present invention;

[0079] Figure 6 The output diagram of a high-temperature geothermal prospect area is shown as an exemplary embodiment of the present invention. Detailed Implementation

[0080] The embodiments of the present invention will be described below with reference to the accompanying drawings and preferred embodiments. Those skilled in the art can easily understand other advantages and effects of the present invention from the content disclosed in this specification. The present invention can also be implemented or applied through other different specific embodiments, and various details in this specification can also be modified or changed based on different viewpoints and applications without departing from the spirit of the present invention. It should be understood that the preferred embodiments are only for illustrating the present invention and not for limiting the scope of protection of the present invention.

[0081] It should be noted that the illustrations provided in the following embodiments are only schematic representations of the basic concept of the present invention. Therefore, the drawings only show the components related to the present invention and are not drawn according to the actual number, shape and size of the components in the actual implementation. In the actual implementation, the form, quantity and proportion of each component can be arbitrarily changed, and the layout of the components may also be more complex.

[0082] In the following description, numerous details are explored to provide a more thorough explanation of embodiments of the invention. However, it will be apparent to those skilled in the art that embodiments of the invention may be practiced without these specific details. In other embodiments, well-known structures and devices are shown in block diagram form rather than in detail to avoid obscuring embodiments of the invention.

[0083] Example 1:

[0084] A method for fusion and prediction of high-temperature geothermal multi-source measurement data in complex tectonic regions, such as Figure 1 As shown, it includes:

[0085] Acquire multi-source measurement data in complex structural areas, and categorize the multi-source measurement data into measurement units according to planar location, depth range, and stratigraphic interface to generate a measurement unit set;

[0086] Anomaly markers are applied to the measurement unit set to generate deep geophysical anomaly units, thermal control structure units, and drilling temperature anomaly units;

[0087] Starting from the deep geophysical anomaly unit, candidate thermal anomaly transmission paths are generated by connecting the adjacent measurement units through the thermal control structure unit segment by segment to the drilling temperature anomaly unit.

[0088] The consistency of the thermal anomaly path is calculated based on the continuous distribution of anomaly units in the candidate thermal anomaly propagation path, the thermal resistance of adjacent measurement units, and the interval of non-responding measurement units. The thermal anomaly path is then selected and generated according to the consistency of the thermal anomaly path.

[0089] The initial temperature field prediction results are generated by performing heat conduction calculations based on the thermophysical property data and radioactive heat generation data in the measurement unit corresponding to the dominant heat control path.

[0090] The temperature measurement residual is generated by calculating the difference between the actual drilling temperature and the initial temperature field prediction results in the same well location and the same depth range. The temperature measurement residual is then segmented and aggregated according to the distance relationship between the measurement unit where the temperature measurement residual is located and different path segments in the heat control dominant path and the layer attribution relationship to generate the temperature measurement residual attribution correction amount.

[0091] Based on the temperature measurement residual attribution correction, the heat conduction calculation conditions of the corresponding path segment are corrected, and the heat conduction calculation is re-performed to generate the corrected temperature field prediction results.

[0092] The minimum burial depth surface reaching the target temperature is extracted from the corrected temperature field prediction results to generate the target isothermal surface burial depth prediction results. Measurement units whose target isothermal surface burial depth meets the preset exploration depth conditions and is connected to the dominant heat control path are merged to generate high-temperature geothermal target area prediction results. In this embodiment, dry hot rock resource prediction data from complex tectonic areas is used as the application data source. Input data includes geological structure data, gravity measurement data, magnetotelluric sounding data, well temperature measurement data, rock thermal property data, and radioactive heat generation data. Among these, the well temperature measurement data is obtained from... Figure 3 The temperature measurement curves shown are provided by the temperature measurement wells; the deep electrical structure and geophysical interpretation data are provided by [the relevant authority / organization]. Figure 4 The deep electrical structure interpretation diagram shown is provided, and the temperature field simulation and target isotherm results are provided by [the relevant authority / organization]. Figure 5 The temperature profiles and isotherm plots shown are provided, and the final target area output results are provided by [the relevant authority / organization]. Figure 6 The output map of the high-temperature geothermal prospect area is provided. Rock thermal property data and radioactive heat generation data are written into the embodiment using text data. The number of rock samples is 98, with thermal conductivity ranging from 1.516 W / (m·K) to 5.066 W / (m·K). Low thermal conductivity strata are used as thermal insulation constraints, and radioactive heat generation data are used as internal heat generation constraints within the strata.

[0093] The present invention is further configured such that, in acquiring multi-source measurement data of complex structural areas, the multi-source measurement data is categorized into measurement units according to planar location, depth range, and stratigraphic interface to generate a measurement unit set, which includes:

[0094] Multi-source measurement data for complex tectonic zones includes geological structural data, gravity measurement data, magnetotelluric sounding data, well temperature measurement data, rock thermal property data, and radiothermal generation data. Specifically, the high-temperature geothermal fusion prediction system for complex tectonic zones includes a data access terminal, a data preprocessing terminal, and a measurement unit construction terminal. Input objects include geological structural data, gravity measurement data, magnetotelluric sounding data, well temperature measurement data, rock thermal property data, and radiothermal generation data. Geological structural data is imported through geological maps, structural interpretation results, stratigraphic profiles, and well exposure records, and must at least include stratigraphic interfaces, fault boundaries, fold cores, lithological contact interfaces, and well exposure strata. Gravity measurement data is imported through gravimeter acquisition files or geophysical interpretation databases, and must at least include measurement point coordinates, Bouguer gravity anomalies, anomaly gradients, or gravity interpretation depth. Magnetotelluric sounding data is imported through MT or AMT inversion result files, and must at least include the location of the measuring line, measurement point coordinates, inversion depth, resistivity profile, and electrical interface. Drilling temperature data is imported via temperature well logging or geothermal monitoring interfaces, and must include at least the well location, measurement depth, measured temperature, and measurement time. Rock thermophysical property data is imported via laboratory testing databases, and must include at least the sample collection location, collection depth, formation, thermal conductivity, specific heat capacity, and thermal diffusivity. Radioactive heat generation data is imported via radioactive element testing results, and must include at least the sample collection location, collection depth, formation, uranium (U) content, thorium (Th) content, potassium (K) content, and heat generation rate. Gravity measurements and magnetotelluric sounding are batch-acquired data and are accessed according to measurement batch; drilling temperature data is depth sequence data and is accessed according to well number and measurement depth; thermophysical property data and radioactive heat generation data are sample point data and are accessed according to sample number. After access, an original multi-source measurement data package is formed. Each data entry in the data package retains the source type, spatial location, depth range, collection time, or test batch identifier, allowing for data source traceability when subsequently assigned to a measurement unit.

[0095] Stratigraphic constraint profile data is established based on stratigraphic interfaces, fault boundaries, and well-exposed strata in geological structural data. Planar profile data is established based on survey point locations, survey line locations, and well locations. Specifically, map coordinate checks, profile depth verification, and well-exposed strata verification are performed on the geological structural data to convert stratigraphic interfaces, fault boundaries, and well-exposed strata into stratigraphic constraint lines or surfaces usable for profile analysis. Stratigraphic interfaces define the upper and lower boundaries between adjacent strata, fault boundaries define the location of structural cuts and the fault relationships between adjacent strata, and well-exposed strata are used to correct the depth of the stratigraphic interface at the borehole location. When a depth difference exists between the stratigraphic interface and the well-exposed strata at the same location, the well-exposed strata are used as the depth verification point near the well location, and the stratigraphic interfaces in the adjacent area are locally adjusted to generate stratigraphic constraint profile data. A planar position check is performed on gravity measurement points, magnetotelluric survey lines, and well locations, projecting the measurement point locations, survey line coverage areas, and well locations onto the same planar reference frame. Planar profile data is generated based on measurement point density, survey line spacing, and well distribution. The boundaries of the planar profile data are jointly defined by the measurement point coverage area, survey line control area, and well influence area, enabling subsequent measurement units to simultaneously accommodate point, line, and well data. Formation-constrained profile data is a depth profile result jointly defined by formation interfaces, fault boundaries, and well-exposed strata; planar profile data is a planar profile result jointly defined by measurement points, survey lines, and well locations. Both provide unit boundaries in the depth and horizontal directions, respectively. This processing allows gravity measurements, magnetotelluric sounding, well temperature measurements, and rock sample test data to be integrated into the same spatial organization frame, providing a unified unit object for subsequent anomaly marking and path connection.

[0096] The planar profile data and stratigraphically constrained profile data are mapped to each other, and measurement units are generated according to the depth intervals defined by the planar profile range and adjacent stratigraphic interfaces. Specifically, each planar profile range is overlaid with the depth interval defined by the adjacent stratigraphic interface below it to form a measurement unit with a planar boundary, upper depth limit, lower depth limit, and stratigraphic attribution. The measurement unit is the basic carrier of multi-source measurement data in this method, used to uniformly accommodate data with different exploration scales, different depth ranges, and different physical meanings. When generating a measurement unit, a unit location identifier, a depth interval identifier, a stratigraphic attribution identifier, and a boundary relationship identifier are configured. Among them, the unit location identifier comes from the planar profile range, the depth interval identifier comes from the adjacent stratigraphic interface, the stratigraphic attribution identifier comes from the stratigraphic interpretation results in the geological structural data, and the boundary relationship identifier is used to record the planar contact, vertical contact, or stratigraphic interface contact relationship between adjacent measurement units. If a fault boundary crosses a measurement unit, the fault cutting state is preserved in that measurement unit; if the measurement unit is located near a lithological contact interface, the lithological contact state is preserved. The measurement unit output by this step is not a single grid data, but a spatial data unit that simultaneously has planar location, depth range and formation interface constraints. Subsequent anomaly marking, path connection and temperature field calculation all use this measurement unit as the data transfer object.

[0097] Gravity measurement data and magnetotelluric sounding data are categorized into geophysical data items by measurement point location, survey line coverage, inversion depth, and detection coverage area, and assigned to corresponding measurement units. Drilling temperature measurement data are categorized into temperature data items by drilling location and temperature measurement depth, and assigned to corresponding measurement units by rock thermal property data, based on sample collection location, sample collection depth, and the stratigraphic layer to which the sample belongs. Radioactive heat generation data are categorized into heat generation data items by sample collection location, sample collection depth, and the stratigraphic layer to which the sample belongs. Geological structural data are categorized into stratigraphic interfaces, fault boundaries, lithological contact interfaces, and drill-exposed layers. Geophysical data items are generated by assigning data to corresponding measurement units. Specifically, for gravity measurement data, the plane subdivision range is determined according to the location of the measurement point, and the corresponding depth interval is determined according to the interpretation depth or influence depth of the gravity anomaly. Gravity anomaly values, gravity anomaly gradients, and gravity anomaly boundaries are assigned to the corresponding measurement units. For magnetotelluric sounding data, the plane subdivision range is determined according to the coverage of the survey line and the location of the measurement point, and the corresponding depth interval is determined according to the inversion depth and electrical abrupt change interface. Resistivity values, electrical abrupt change interfaces, and inversion stratigraphic levels are assigned to the corresponding measurement units. Gravity measurement data and magnetotelluric sounding data together constitute geophysical data items. For well temperature measurement data, the plane subdivision range is determined according to the well location, and the depth interval is determined according to the temperature measurement depth. Measured temperature, temperature measurement depth, and temperature curve segments are assigned to the corresponding measurement units to form temperature data items. For rock thermal property data, the planar subdivision range is determined according to the sample collection location, and the depth interval is determined according to the sample collection depth and the strata to which the sample belongs. Thermal conductivity and its testing source are then assigned to the corresponding measurement unit to form a thermal property data item. For radiogenerated heat data, the sample is assigned to the corresponding measurement unit according to the sample collection location, sample collection depth, and the strata to which the sample belongs. The heat generation rate and its testing source are written into the heat generation data item. For geological structural data, the spatial relationship between stratigraphic interfaces, fault boundaries, lithological contact interfaces, and well-exposed strata and the measurement unit is used for assignment. Stratigraphic attribution, fault penetration status, lithological contact status, and well-exposed strata are written into the geological attribute data item. If multiple similar data records exist within the same measurement unit, the original record and its source identifier are retained, and an index is created according to the data type. If linear measurement data covers multiple measurement units, it is assigned to the covered measurement units according to the coverage area and inversion depth, while retaining the coverage relationship. Through this process, geophysical data items are used for subsequent deep geophysical anomaly unit marking, temperature data items are used for subsequent drilling temperature anomaly unit marking, thermal property data items are used for subsequent path thermal resistance and heat conduction calculations, heat generation data items are used for subsequent internal heat generation constraints within the strata, and geological attribute data items are used for subsequent heat-controlling structural unit marking.

[0098] A measurement unit set is generated based on measurement units containing geophysical data, temperature data, thermal property data, thermal generation data, and geological attribute data. Specifically, measurement units with completed data item inclusion are organized and arranged sequentially, and a measurement unit index is established according to planar location, depth range, and stratigraphic affiliation. Measurement units containing all three data items are marked as complete measurement units; for measurement units containing only partial data items, existing data items are retained and missing data types are marked, and supplemented or verified later in anomaly marking or heat conduction calculations based on existing data from the same stratigraphy, lithology, or adjacent depth ranges. The output measurement unit set includes measurement unit identifier, planar location, depth range, stratigraphic affiliation, geophysical data items, temperature data items, thermal property data items, thermal generation data items, and geological attribute data items. This is for application verification. Figure 3 The temperature measurement curves shown can be used to generate temperature data items based on drilling location and measurement depth. Figure 4 The interpretation results of the deep electrical structure shown can be used to form geophysical data items according to the coverage of the survey line and the inversion depth. Figure 5 The temperature profiles and isotherm results shown can be assembled using the measurement unit set in subsequent temperature field calculations. Figure 6 The output results of the high-temperature geothermal prospect area shown can be obtained by merging measurement units in the subsequent target area generation step. Thus, the measurement unit set can transform point temperature measurement data, linear geophysical profile data, sample test data, and geological structure data into a dataset with unified spatial boundaries and stratigraphic attributes, providing a continuous data foundation for subsequent anomaly unit labeling, candidate thermal anomaly propagation path generation, path consistency screening, temperature measurement residual attribution correction, and high-temperature geothermal target area prediction.

[0099] The present invention is further configured such that the step of generating deep geophysical anomaly units, thermal control structure units, and drilling temperature anomaly units by anomaly marking the measurement unit set includes:

[0100] Geophysical data items, geological attribute data items, and temperature data items are extracted from the measurement unit set. Specifically, at the data processing end of the high-temperature geothermal multi-source measurement data fusion and prediction system, the measurement unit set generated in the previous step is used as the input object. Each measurement unit contains at least a unit location identifier, depth interval identifier, stratigraphic attribution identifier, and data type index. When reading data items, gravity anomaly values, gravity anomaly gradients, resistivity values, resistivity interface depths, and geophysical inversion stratigraphic levels are extracted from the geophysical data items; fault penetration status, fold core location relationships, basalt layer attribution, low thermal conductivity strata attribution, and lithological contact interfaces are extracted from the geological attribute data items; and well number, temperature measurement depth, measured temperature, temperature measurement time, and temperature curve segment are extracted from the temperature data items. Before reading, the data items undergo integrity screening. Data items lacking spatial location or depth intervals are not included in the anomaly marking, and data items lacking measurement batch or source identifiers are retained as data to be verified. Anomaly peak removal and adjacent depth layer smoothing are performed on gravity and resistivity data; sequential verification and duplicate measurement point merging are performed on drilling temperature data; and consistency verification between formation interfaces and drill-revealed strata is performed on geological attribute data. The geophysical data items, geological attribute data items, and temperature data items generated in this step are then entered into three subsequent marking channels to avoid direct mixing of different data types before anomaly marking.

[0101] Geophysical data items are stratified according to the detection depth range or inversion depth layer, and measurement units corresponding to changes in deep density or deep electrical interfaces are identified and marked as deep geophysical anomaly units. Specifically, for gravity measurement data, the depth response range reflected by gravity anomalies is divided into shallow response layers, middle response layers, and deep response layers; for magnetotelluric sounding data, the electrical response layers corresponding to the depth intervals of the measurement units are divided according to the inversion depth layer. For measurement units within the same depth layer, the gravity anomaly gradient and resistivity abrupt change amplitude are extracted to form a deep geophysical response record. Deep geophysical response records are used to identify changes in deep density and deep electrical interfaces. The judgment logic is as follows: within the same depth layer, if the gravity anomaly gradient of a measurement unit is outside the background fluctuation range of that layer, or if the resistivity abrupt change amplitude of a measurement unit is outside the background fluctuation range of that layer, then a geophysical response anomaly record is generated. If this record is located in the middle or deep response layer and forms a continuous anomaly band with adjacent measurement units, then the corresponding measurement unit is marked as a deep geophysical anomaly unit. To make the threshold source transparent, let the first... The gravity anomaly gradient background threshold for each depth layer is: The background threshold for resistivity abrupt change is Both are determined by the stable background range of the gravity anomaly gradient and the resistivity abrupt change amplitude within that depth layer, respectively; measurement unit The deep geophysical anomaly marker is denoted as ,when Located in the middle or deep response layer, and satisfying the gravity anomaly gradient exceeding... Or the magnitude of the resistivity change exceeds At that time, This is designated as a deep geophysical anomaly. The above judgment rules are used to explain the triggering conditions for anomaly marking. and The data originates from the distribution of data within the same layer and geophysical interpretation results. After marking is completed, the system outputs deep geophysical anomaly unit data containing measurement unit identifiers, response depth layers, anomaly source types, and continuous anomaly band numbers, providing input for the selection of starting points for subsequent candidate thermal anomaly propagation paths;

[0102] The geological attribute data items are analyzed to identify heat-controlling geological elements. Measurement units containing fault penetration, fold cores, basalt layers, low thermal conductivity strata, or lithological contact interfaces are marked as heat-controlling structural units. Specifically, structural element analysis is performed on the geological attribute data items. First, the intersection status between the measurement unit and the fault boundary is read. Measurement units with fault penetration, fault adjacency, or fault influence zone coverage generate fault heat-controlling records. Next, the spatial relationship between the measurement unit and the fold core is read. Measurement units located in or near the fold core generate fold heat-controlling records. Then, the lithology and stratigraphic affiliation of the measurement unit are read. Measurement units belonging to basalt layers generate local heat transfer disturbance records, measurement units belonging to low thermal conductivity strata such as mudstone, shale, and coal seams generate insulation layer records, and measurement units located at lithological contact interfaces with high and low thermal conductivity generate lithological contact records. The identification of heat-controlling geological elements does not rely on a single structural name as the marking basis, but rather on whether the measurement unit has geological attributes such as heat conduction, heat accumulation, heat insulation, or heat transfer disturbance. Thermal conductivity corresponds to fracture penetration, thermal accumulation corresponds to fold cores, thermal insulation corresponds to low thermal conductivity strata, and thermal disturbance corresponds to basalt layers or lithological contact interfaces. When multiple thermally controlling geological elements exist simultaneously in the same measurement unit, the element type is retained and a composite thermally controlling structural marker is generated. After marking, the system outputs thermally controlling structural unit data, which includes the measurement unit identifier, thermally controlling element type, stratum, and adjacent geological boundaries. This output, along with deep geophysical anomaly units, enters the subsequent path connection step, enabling deep geophysical anomalies to continue to propagate along units with geological heat transfer significance.

[0103] The drilling temperature curves of the temperature data items are segmented for identification. Measurement units containing high-temperature segments or abrupt temperature gradient changes are marked as abnormal drilling temperature units. Specifically, a well temperature sequence is established for each temperature data item according to the drilling well number. Each well temperature sequence is sorted from shallow to deep according to the measurement depth. For repeated temperature measurements at the same depth, valid temperature values ​​are selected based on measurement time and measurement stability, and the valid temperature values ​​and their source records are retained. The well temperature sequence is then divided into several depth segments, and the temperature change trend with depth within each depth segment is calculated. Two types of temperature responses are identified: the first type is the high-temperature segment, i.e., the depth segment where the temperature in the same well or adjacent wells is significantly higher than the background temperature of the same layer; the second type is the abrupt temperature gradient segment, i.e., the depth segment where the temperature increase trend of adjacent depth segments changes significantly. To ensure the verification of the identification rules, let the... The measured temperature at each depth segment was Temperature gradients are formed by temperature changes in adjacent depth sections. The upper limit of the background temperature in the same layer is denoted as The upper limit of the background temperature gradient in the same layer is denoted as ;when Exceed When, the corresponding depth segment is recorded as the high temperature segment, when Exceed If the gradient change of adjacent depth segments exceeds the stable background range, the corresponding depth segment is recorded as the temperature gradient abrupt change segment. and Determined by stable temperature measurement sections within the same layer or well, the data originates from measured temperature data. The depth ranges containing high-temperature segments and abrupt temperature gradient changes are mapped back to the corresponding measurement units, generating drilling temperature anomaly markers. Upon completion, the system outputs the drilling temperature anomaly unit, containing the well number, anomaly depth segment, anomaly type, and corresponding measurement unit identifier. As an example, if a measurement unit is marked as a deep geophysical anomaly unit in the deep geophysical exploration channel, and adjacent units contain fracture penetration and low thermal conductivity formation markers, and high-temperature values ​​or abrupt temperature gradient changes occur within the same drilling depth segment, then this set of markers can form a continuous data basis from deep geophysical response to drilling temperature response in subsequent path generation. Compared to screening anomalies solely based on single-point high-temperature values, this embodiment simultaneously retains unit markers from three sources: geophysical response, thermally controlling geological elements, and well temperature response. This reduces the interference of isolated measurement points or single geophysical anomalies on subsequent path selection. Figure 3 As shown, the drilling temperature measurement curve forms a temperature change sequence along the depth direction. The system marks drilling temperature anomaly units according to the high temperature segment and the temperature gradient abrupt change segment. Figure 3 The middle part of the temperature measuring holes has a high heat flux value and obvious temperature gradient change, which can be used as a data source for temperature measurement residual calculation and temperature anomaly identification. Figure 4The deep electrical structure interpretation diagram shown is used to identify changes in deep electrical interfaces and to mark the corresponding measurement units as deep geophysical anomaly units. Figure 4 The regions corresponding to changes in the electrical structure and the structural interface are used as the basis for identifying thermal control structural units.

[0104] The present invention is further configured such that the process of generating a candidate thermal anomaly transmission path, starting from a deep geophysical anomaly unit and connecting segment by segment along the connection direction of adjacent measurement units through heat-controlling structural units to a drilling temperature anomaly unit, includes:

[0105] Using deep geophysical anomaly units as the starting units for the path and extracting drilling temperature anomaly units as the ending units, specifically, in the high-temperature geothermal multi-source measurement data fusion and prediction system for complex structural areas, the set of anomaly measurement units is used as the input object. The unit location identifier, depth interval identifier, formation attribution identifier, deep geophysical anomaly marker, heat-controlling structure marker, and drilling temperature anomaly marker for each measurement unit are read. Measurement units with deep geophysical anomaly markers are then assigned to the starting unit set. Measurement units marked with drilling temperature anomalies are grouped into the termination unit set. In this process, deep geophysical anomaly markers are derived from gravity response stratification results or magnetotelluric inversion depth layers, while drilling temperature anomaly markers are derived from high-temperature segments or temperature gradient abrupt change segments. Depth screening is performed on the starting unit set, retaining measurement units located in the mid-deep or deep response layers; well location validity screening is performed on the ending unit set, retaining data with complete well numbers, depth ranges, and measured temperature sources. This step establishes two-end constraints for the path search: the starting end corresponds to the deep geophysical response, and the ending end corresponds to the drilling temperature response. These two-end constraints limit the path generation range, preventing path searches from starting from arbitrary geophysical anomalies or arbitrary temperature measurement points, thus ensuring that subsequent path connections have clear input boundaries.

[0106] An adjacent unit connection table is established based on the shared planar boundaries, shared depth boundaries, and stratigraphic contact interfaces between measurement units. Specifically, spatial adjacency resolution is performed on the measurement unit set to read the planar subdivision range, depth interval, and stratigraphic interface information of each measurement unit. When two measurement units share a common planar boundary, they are recorded as lateral adjacency; when they share a common depth boundary, they are recorded as vertical adjacency; when two measurement units are located on opposite sides of the same stratigraphic contact interface, they are recorded as stratigraphic contact adjacency; when two measurement units are cut by the same fault boundary or located within the same fault influence zone, they are recorded as structural adjacency. The adjacent unit connection table is denoted as... ,in, Indicates measurement unit and measurement unit There exists at least one connectable relationship. This indicates that the two are not connected. This connection table only describes the spatial and stratigraphic connection basis between measurement units and does not directly determine the thermal anomaly path. Before establishing the connection table, overlapping boundaries, fault lines, and fault displacement positions at stratigraphic interfaces are preprocessed; when multiple connection relationships exist for the same pair of measurement units, a list of connection types is retained so that subsequent segment-by-segment connections can distinguish between lateral, vertical, and tectonic transfers. The adjacent unit connection table output in this step serves as the data carrier for path connections; subsequent steps only perform segment-by-segment access between units with connection relationships in the table.

[0107] Starting from the initial unit of the path, adjacent measurement units with heat-controlling structure markers or drilling temperature anomaly markers are connected segment by segment along the adjacent unit connection table to generate a unit connection sequence; specifically, starting from the initial unit set... any path starting unit For the current cell, read the adjacent cell connection table and... Adjacent measurement units are checked for thermal control structural markers or drilling temperature anomaly markers. Thermal control structural markers include thermally controlling geological elements such as fracture penetrations, fold cores, basalt layers, low thermal conductivity strata, or lithological contact interfaces. Drilling temperature anomaly markers include high temperature segments or temperature gradient abrupt change segments. Measurement units that satisfy the adjacent connection relationship and have the above markers are added to the current path, generating an initial unit connection sequence. The newly added measurement unit is then used as the current unit, and the process continues to read the next layer of adjacent measurement units along the adjacent unit connection table, repeating the connection relationship verification and marker verification to form a progressively extending unit connection sequence. To avoid path loops, measurement units already in the current unit connection sequence are not added again; to avoid the path deviating from the thermal control direction, several consecutive adjacent units without thermal control structural markers are not considered for extension. The path connection rule can be expressed as: if ,and or Then it is allowed by Access ,in, Indicates measurement unit It has heat-controlling structure markings. Indicates measurement unit It includes drilling temperature anomaly markers; this expression describes the triggering conditions for segment-by-segment access, which are derived from the anomaly marker results of the previous step and the adjacent cell connection table. The cell connection sequence output in this step includes the starting cell, the cells traversed, the connection type, and the current termination status, providing intermediate data for recording candidate thermal anomaly propagation paths;

[0108] When a unit connection sequence passes through at least one thermally controlled structural unit and connects to a path termination unit, the unit connection sequence is recorded as a candidate thermal anomaly propagation path. Specifically, termination determination and path recording are performed on each segment of the generated unit connection sequence. If the first end of a unit connection sequence is a deep geophysical anomaly unit, the last end is a drilling temperature anomaly unit, and the sequence contains at least one thermally controlled structural unit, then the unit connection sequence is recorded as a candidate thermal anomaly propagation path. The candidate thermal anomaly propagation path is denoted as... ,in, As the starting unit of the path, The intermediate units in the sequence are arranged in the order of connection, serving as the path termination unit. During recording, the path number, starting unit identifier, ending unit identifier, type of heat-controlling structure traversed, adjacent connection type, and path length are saved simultaneously. If the same starting unit can connect to the same or different drilling temperature anomaly units through different heat-controlling structure units, they are recorded as different candidate paths. If a unit connection sequence connects to a drilling temperature anomaly unit without passing through a heat-controlling structure unit, it is not recorded as a candidate thermal anomaly transmission path. If a unit connection sequence passes through a heat-controlling structure unit but cannot connect to a drilling temperature anomaly unit, it is recorded as an unclosed path and reserved as a path to be retested. As a verification example, if a deep geophysical anomaly unit connects to an adjacent fracture-cutting unit, the fracture-cutting unit continues to connect to a low thermal conductivity formation unit, and the low thermal conductivity formation unit finally connects to a unit containing a drilling temperature gradient abrupt change segment, then this unit sequence satisfies the closure condition of deep geophysical anomaly—heat-controlling structure—drilling temperature anomaly and is recorded as a candidate thermal anomaly transmission path. This processing prevents isolated deep geophysical anomalies, unconnected temperature anomalies, and short paths that do not pass through heat-controlling structures from being included in subsequent path consistency calculations. Subsequent screening steps can then be performed on a set of candidate paths with geological heat transfer data, such as... Figure 4 As shown, the interpretation results of deep electrical structures are used to constrain the positional relationship between deep geophysical anomaly units and thermal control structural units. The system uses the deep geophysical anomaly unit as the starting constraint for the path, and... Figure 3 The drilling temperature anomaly unit serves as a path termination constraint. By sequentially connecting adjacent measurement units to measurement units marked with heat-controlling structures, candidate thermal anomaly propagation paths are generated. If the unit connection sequence can be generated by… Figure 4 The location of the deep geophysical anomaly is connected to the thermal control structural unit. Figure 3 If the location of the drilling temperature anomaly is recorded, it will be recorded as a candidate thermal anomaly transmission path; if there is only an isolated geophysical anomaly or an isolated temperature anomaly, it will not be selected as a priority candidate for the thermal control dominant path.

[0109] The present invention is further configured such that the step of calculating the thermal anomaly path consistency degree based on the continuous distribution of anomaly units in the candidate thermal anomaly propagation path, the thermal resistance of adjacent measurement units, and the interval of non-responding measurement units, and selecting and generating the dominant thermal control path according to the thermal anomaly path consistency degree includes:

[0110] A path unit sequence is generated according to the path connection order in the candidate thermal anomaly transmission paths. Specifically, in the high-temperature geothermal multi-source measurement data fusion prediction system for complex structural areas, the set of candidate thermal anomaly transmission paths obtained in the previous step is used as the input object. Each candidate thermal anomaly transmission path includes a path starting unit, a path passing unit, and a path ending unit. The path starting unit originates from deep geophysical anomaly units, the path ending unit originates from drilling temperature anomaly units, and the path passing units originate from heat-controlling structural units or adjacent transition units sequentially connected from the adjacent unit connection table. A path integrity check is performed on the candidate thermal anomaly transmission paths. The check items include whether the starting unit has a deep geophysical anomaly marker, whether the ending unit has a drilling temperature anomaly marker, whether there are shared plane boundaries, shared depth boundaries, or formation contact interfaces between adjacent units, and whether there is at least one heat-controlling structural unit in the path. The candidate paths that pass the check are organized into a path unit sequence according to the connection order, denoted as . ,in, As the starting unit of the path, This is the path termination unit. For the first Each measurement unit is a separate measurement unit. Duplicate measurement units are deduplicated from the path. If deduplication results in a path break, the candidate path is added to the path verification set; if adjacent connections are maintained after deduplication, the path unit sequence is retained. The output of this step is path unit sequence data, including path number, unit sequence number, measurement unit identifier, unit connection type, and unit anomaly marker. This provides a unified path carrier for subsequent extraction of abnormal continuous segments, path thermal resistance, and no-response intervals.

[0111] Anomaly continuum data is generated by merging the anomaly markers of measurement units in the path unit sequence. Specifically, the anomaly marker field is read for each measurement unit in the path unit sequence. The anomaly marker field includes deep geophysical anomaly markers, thermal control structure markers, and drilling temperature anomaly markers. Measurement units with at least one anomaly marker are designated as response units, and measurement units without the above anomaly markers are designated as non-response measurement units. To clarify the source of the anomaly continuum, let the measurement unit... The response status is ,when When any one of the following is present: deep geophysical anomaly marker, heat-controlling structure marker, or drilling temperature anomaly marker, When none of the three types of markers exist, This response status is only used to indicate the identification of abnormal contiguous segments and is not considered an independent user-defined parameter. Adjacent segments are connected according to the connection order of the path unit sequence. The measurement units are grouped into the same abnormal continuous segment, and each abnormal continuous segment is recorded as... ,in, and These are the start and end measurement units for the continuous anomaly segment. For each continuous anomaly segment, the combination of anomaly types it contains is recorded synchronously. If a continuous segment simultaneously contains deep geophysical anomalies, thermal control structures, and drilling temperature anomalies, it is recorded as a composite anomaly segment; if it contains only one or two of these anomaly types, it is recorded as a single-type or double-type anomaly segment. This step outputs the anomaly segment data, including the segment number, start and end units, segment length, anomaly type combination, and corresponding path number. Through this processing, subsequent calculations can distinguish between multiple anomalies occurring consecutively along the same path and a single anomaly existing in isolation, providing a traceable data source for the consistency of thermal anomaly paths.

[0112] Based on the thermal property data and depth range of adjacent measurement units in the path unit sequence, path thermal resistance data is generated to identify the interval data of non-responding measurement units between adjacent abnormal continuous segments; specifically, for two adjacent measurement units in the path unit sequence... and Read the thermal property data items and depth interval data. The thermal property data items must include at least thermal conductivity. Measurement units lacking thermal conductivity data should have it supplemented using existing thermal conductivity data from the same stratum, lithology, or adjacent depth intervals, and the source of the supplementation should be retained in the path thermal resistance data. The thermal resistance between adjacent units is determined based on the center distance between adjacent units, the depth interval thickness, and the equivalent thermal conductivity. The equivalent thermal conductivity is determined according to the lithology and thermal conductivity of two adjacent measurement units. If the two units belong to the same stratum and have the same lithology, the same-layer thermal conductivity is used; if the two units are located in different strata or on opposite sides of a lithological interface, the interlayer equivalent value of the thermal conductivity on both sides is used. To clarify the source of thermal resistance, let the thermal resistance of adjacent units be denoted as . The formation rule is: the greater the heat transfer distance between adjacent units, the better. The larger the value, the higher the equivalent thermal conductivity of adjacent units. The smaller the value, the higher the thermal resistance of the path segment when low thermal conductivity formation units participate in the connection. Path thermal resistance data is generated by arranging the thermal resistances of all adjacent units in the order of path connection. Subsequently, non-response measurement units between abnormal contiguous segments are identified, when there is a non-response between two adjacent abnormal contiguous segments. When measuring units, these units are merged and recorded as non-response intervals, and non-response measurement unit interval data is generated based on their number, cumulative thickness, or cumulative path length. This step outputs path thermal resistance data and non-response measurement unit interval data, so that the consistency of thermal anomaly paths simultaneously reflects heat transfer resistance and the degree of evidence discontinuity;

[0113] A thermal anomaly path consistency score is generated based on anomaly continuum data, path thermal resistance data, and interval data of non-responding measurement units. Candidate thermal anomaly propagation paths are then selected based on this consistency score to generate the dominant heat-controlling path. Specifically, for each candidate thermal anomaly propagation path, anomaly continuum contribution is first calculated based on anomaly continuum data, then thermal resistance contribution is calculated based on path thermal resistance data, and finally interval penalty is calculated based on non-responding measurement unit interval data. The thermal anomaly path consistency score is denoted as... This is used to indicate whether deep geophysical anomalies, heat-controlling structures, and drilling temperature anomalies within a candidate thermal anomaly propagation path continuously correspond along a low thermal resistance path. One implementation is as follows: ,in, The contribution to anomaly continuity is formed by the proportion of anomaly continuity segments in the path and the combination of anomaly types; the longer the anomaly continuity segment and the more compound anomaly continuity segments, the better. The higher; The contribution to thermal resistance is formed from path thermal resistance data; the lower the total path thermal resistance, the better. The higher; The contribution to the no-response interval is formed by the interval data of the no-response measurement units; the shorter the no-response interval, The higher the value, the better; all three are normalized to the same range before being multiplied to avoid a single strong anomaly masking the problem of excessively high path thermal resistance or excessively long anomaly intervals. Trigger threshold for thermal anomaly path consistency. The consistency distribution of the candidate thermal anomaly propagation paths is determined. All candidate paths are sorted from highest to lowest consistency, and the lowest consistency among the top 25% of paths after sorting is taken as the baseline. When the number of candidate paths is less than four, the average consistency of the candidate paths is used as the initial threshold, and manual verification is performed in conjunction with whether the path passes through fracture penetration units and drilling temperature anomaly units. Sort the candidate paths from highest to lowest, Candidate thermal anomaly propagation paths are selected as the dominant heat-controlling paths. If multiple candidate paths correspond to the same well temperature anomaly unit, the path with the highest thermal anomaly path consistency and the most complete passage through the heat-controlling structural unit type is prioritized. After selection, the dominant heat-controlling path data is output, including path number, path unit sequence, thermal anomaly path consistency, continuous anomaly segments traversed, path thermal resistance summary, and no-response interval record. As a verification example, if the candidate path... Starting from a deep geophysical anomaly unit, passing through a fracture-cutting unit and a low thermal conductivity formation unit, and connecting to a drilling temperature anomaly unit, if the path has a long continuous anomaly segment, low thermal resistance, and short no-response intervals, then its If the value is higher than the threshold, it is output as the dominant heat control path; if the candidate path Although it includes deep geophysical anomaly units and drilling temperature anomaly units, if it lacks continuous support from heat-controlling structures and has a long interval of no response, its consistency is below the threshold and it is retained as a weakly supported path. This screening result provides a heat-controlling dominant path for subsequent temperature field prediction steps, allowing heat conduction calculations to preferentially load path constraints with continuous anomaly evidence and reasonable heat transfer conditions;

[0114] The present invention is further configured such that the initial temperature field prediction result is generated by performing heat conduction calculations based on the thermophysical property data and radioactive heat generation data in the measurement unit corresponding to the dominant heat control path, including:

[0115] The system extracts measurement units traversed by the dominant heat-controlling path to generate a sequence of temperature field calculation units. Specifically, in the high-temperature geothermal multi-source measurement data fusion prediction system for complex structural areas, the dominant heat-controlling path data and measurement unit sets are used as input objects. The system's data processing end reads the path number, path connection order, measurement unit identifier, and path segment type of the dominant heat-controlling path, and arranges the measurement units traversed by each dominant heat-controlling path in the connection order from deep geophysical anomaly units to drilling temperature anomaly units, generating a sequence of temperature field calculation units. Each temperature field calculation unit retains its planar location, depth range, formation affiliation, heat-controlling structural marker, geophysical anomaly marker, and drilling temperature anomaly marker. For the case where the same measurement unit is traversed by multiple dominant heat-controlling paths, the association number between the measurement unit and each path is retained; for measurement units lacking thermal property data or heat generation data in the path, they are marked as units to be supplemented and proceed to the next step of data supplementation processing. This step processes the selected dominant heat-controlling paths and outputs a sequence of temperature field calculation units, so that heat conduction calculations are not directly directed to discrete measurement points or single profiles, but to unit sequences with path order, layer location, and heat-controlling attributes.

[0116] The system reads thermal property data items, radioactive heat generation data items, stratigraphic attribution identifiers, and heat-controlling structure markers from the temperature field calculation unit sequence to generate path heat conduction constraint data. Specifically, the system reads thermal property data items unit by unit from the temperature field calculation unit sequence, including at least thermal conductivity, specific heat capacity, and thermal diffusivity; reads radioactive heat generation data items, including at least U content, Th content, K content, or heat generation rate converted from the above element contents; reads stratigraphic attribution identifiers to confirm the stratigraphic position of the measurement unit and its relationship with adjacent stratigraphic positions; and reads heat-controlling structure markers to confirm whether the measurement unit belongs to a fault-cutting unit, a low thermal conductivity stratigraphic unit, a basalt unit, or a common surrounding rock unit. For measurement units lacking thermal conductivity, thermal conductivity is supplemented according to the priority order of the same stratigraphy, the same lithology, and adjacent depth intervals, and the source of the supplementation is recorded. For measurement units lacking radioactive heat generation data, existing heat generation rate data of the same stratigraphy or the same lithology is supplemented; if supplementation is not possible, a blank heat generation constraint is set and a check mark is retained. After reading and completion, path heat conduction constraint data is generated, including element number, layer number, thermal conductivity, specific heat capacity, thermal diffusivity, heat generation rate, heat control structure type, and data source identifier. This step ensures that the data items required for heat conduction calculations have a clear source, avoiding gaps in the chain between thermal property data, heat generation data, and heat control structure markers;

[0117] Based on path heat conduction constraint data, fracture penetration units are configured as heat conduction channel constraints, low thermal conductivity strata units as heat insulation constraints, basalt strata units as local heat transfer disturbance constraints, and radioactive heat generation data items as intra-stratum heat generation constraints. Specifically, the system classifies and configures path heat conduction constraint data according to the heat-controlling structural markers of the temperature field calculation units. For fracture penetration units, heat conduction channel constraints are configured, represented by equivalent thermal conductivity adjustment along the fracture direction or heat flow input at the fracture boundary. For low thermal conductivity strata units, heat insulation constraints are configured, represented by a decrease in equivalent thermal conductivity between the unit and its adjacent strata or an increase in inter-stratum thermal resistance. For basalt strata units, local heat transfer disturbance constraints are configured, represented by an equivalent thermal conductivity correction term or a local heat transfer correction term for the basalt layer. For units with radioactive heat generation data items, intra-stratum heat generation constraints are configured, which apply the heat generation rate to the corresponding stratum measurement unit. The equivalent thermal conductivity is denoted as Its thermal conductivity is determined by actual measurement or supplementation. Determined together with the corresponding heat-controlling structure type; when When it is a fracture-cutting unit, Adjust according to the heat conduction channel constraint, when When it is a low thermal conductivity formation unit, Adjust according to insulation constraints, when When it is a basalt layer unit, Adjustments are made based on local heat transfer disturbance constraints, while ordinary units retain their original thermal property constraints. The aforementioned adjustment ranges are derived from measured thermal property data, stratigraphic lithology, fracture interpretation results, and temperature measurement fitting results, without setting separate empirical coefficients without a source. The initial heat conduction calculation conditions formed through this step include the equivalent thermal conductivity field, the internal heat generation rate field of the stratigraphic layer, the boundary of the fracture heat conduction channel, and local heat transfer disturbance units. This can transform the geothermal mechanism described in the report as "fracture heat conduction as the main factor and sedimentary cap layer insulation as a secondary factor" into loadable calculation conditions.

[0118] Heat conduction calculations are performed based on constraints such as heat conduction channels, insulation, local heat transfer disturbances, and internal heat generation within the formation. This generates an initial temperature field prediction. Specifically, the system loads the initial heat conduction calculation conditions into the measurement unit set, establishes a heat conduction calculation node for each measurement unit, and establishes heat transfer connections based on the planar adjacency, depth adjacency, and formation contact relationships between adjacent measurement units. Initial Temperature Field Satisfies steady-state heat conduction relationship: ,in, The equivalent thermal conductivity field is determined by the constraints of the heat conduction channel, the insulation constraint, and the local heat transfer disturbance constraint. This represents the heat generation rate field corresponding to the heat generation constraints within the strata. This is the initial temperature field prediction result. The surface boundary uses surface temperature or shallow thermal measurement boundary, the bottom boundary uses regional heat flow or existing heat flow interpretation results, and the lateral boundary is configured as an adiabatic boundary or a given heat flow boundary according to the boundary conditions of the study area. During the calculation, fracture penetration units participate in heat transfer as priority heat conduction channels, low thermal conductivity strata units limit upward heat loss through higher interlayer thermal resistance, basalt strata units adjust the heat distribution near the corresponding strata through local heat transfer perturbation constraints, and radioactive heat generation data terms are loaded as supplementary terms for intralayer heat sources. After the calculation is completed, the system outputs the initial predicted temperature of each measurement unit and assembles it according to planar location, depth range, and stratum interface to generate the initial temperature field prediction result. As a verification example, if a dominant heat-controlling path passes through fracture penetration units, low thermal conductivity strata units, and extends to drilling temperature anomaly units, the predicted temperature near this path in the initial temperature field should form a temperature distribution that continuously rises along the path; comparing this distribution with the drilling measured temperature in subsequent residuals can provide input for the temperature measurement residual attribution correction; such as Figure 5 As shown, the temperature profile and isotherm results are used to verify the initial temperature field prediction output. The system performs heat conduction calculations based on the thermophysical property data and radioactive heat generation data in the measurement unit corresponding to the dominant heat control path, generating the initial temperature field prediction results. Figure 5The burial depths corresponding to the 100℃ and 150℃ isotherms are used as references for extracting the target isothermal surface. The 100℃ temperature response is located at a depth of approximately 3100m, and the 150℃ temperature response is located at a depth of approximately 7600m. The burial depth of the 150℃ isotherm is approximately 7100m to 7900m. The initial temperature field prediction results output in this step are then fed into the subsequent temperature measurement residual calculation step, forming a closed-loop data transfer from the thermal control dominant path constraint to the temperature field prediction.

[0119] The present invention is further configured such that the calculation of the temperature measurement residual and the segmentation and aggregation of the temperature measurement residual to generate the temperature measurement residual attribution correction amount include:

[0120] The system extracts predicted well location temperature data corresponding to drilling temperature data from the initial temperature field prediction results. This predicted temperature data is then mapped to the actual measured temperature based on the well number and depth range, generating corresponding well temperature data. Specifically, in the high-temperature geothermal multi-source measurement data fusion prediction system for complex structural areas, the initial temperature field prediction results and drilling temperature data serve as inputs. The initial temperature field prediction results are output from the previous heat conduction calculation module, including the planar location, depth range, formation affiliation, and initial predicted temperature of each measurement unit. The drilling temperature data is imported from the drilling temperature record library or geothermal monitoring interface, including the well number, wellhead coordinates, measurement depth, actual temperature, measurement time, and measurement source. After data reading, a well location check is performed first, mapping the well number to the planar measurement unit; then a depth range check is performed, assigning the measurement depth to measurement units within the same depth range. For cases where multiple temperature records exist for the same well location and depth range, valid measured temperatures are generated based on measurement time, measurement stability, and abnormal peak filtering rules. Well location extraction of temperature field prediction results does not alter the original temperature field distribution; it only extracts the initial predicted temperature for the corresponding depth range from measurement units that coincide with or are adjacent to the well location. After mapping, well temperature mapping data is generated, which includes at least the well number, depth range, measurement unit identifier, measured well temperature, predicted well location temperature, and formation affiliation identifier. This step establishes a one-to-one correspondence between the measured well temperature and the initial temperature field prediction results within the same spatial location and depth range, providing traceable input for subsequent residual calculations.

[0121] For well temperature data with the same well number and depth range, the measured temperature and predicted temperature at the well location are processed to generate temperature residual records. Specifically, the well temperature data for the same well and depth range are processed to generate temperature residual records. Let the well number be... The depth range is The actual measured temperature during drilling was The predicted temperature at the well location is The temperature measurement residual is recorded as Then we have: ,in, For temperature measurement residuals, positive values ​​indicate that the measured temperature is higher than the initial predicted temperature, and negative values ​​indicate that the measured temperature is lower than the initial predicted temperature. After difference processing, the residual records undergo validity screening. The screening criteria include the allowable error of the temperature measuring instrument, the continuity of temperature changes in adjacent depths within the same well, and the reliability of the temperature measurement data source. Let the valid threshold for the residual be... , The accuracy is jointly determined by the accuracy of the temperature measuring instrument and the accuracy of the temperature field calculation; when Less than When the residual record is weak, it is marked as a weak residual record and used only for statistical verification; when Not less than At this point, the residual record is marked as a valid temperature measurement residual record and enters the segmented aggregation process. The temperature measurement residual record must include at least the well number, depth range, measurement unit identifier, residual value, residual direction, formation attribution, and validity identifier. This step transforms the deviation between the initial temperature field and the measured temperature data into a locationable and attributable residual record, ensuring that subsequent corrections are not directly applied to the entire temperature field, but rather to residual objects with well location and formation origin.

[0122] The dominant heat-controlling path is divided into heat-conducting, heat-insulating, and heat-transferring disturbance sections based on heat-controlling structural markers. Specifically, the path unit sequence, heat-controlling structural markers, stratigraphic affiliation identifiers, and adjacent unit connection relationships are extracted from the dominant heat-controlling path data. When a path unit has a fault penetration marker or a fault influence zone marker, the corresponding path segment is classified as a heat-conducting segment; when a path unit belongs to low-thermal-conductivity strata such as mudstone, shale, or coal seams, or is located at the interface between low-thermal-conductivity strata, the corresponding path segment is classified as a heat-insulating segment; when a path unit belongs to basalt layers, highly porosity igneous layers, or is near a lithological contact interface, the corresponding path segment is classified as a heat-transferring disturbance section. A path segment can consist of continuous measurement units of the same type, and the path segment boundary is determined by changes in heat-controlling structural markers, stratigraphic interfaces, or path connection types. Path segment data includes at least the path number, path segment number, path segment type, starting and ending measurement units, corresponding strata, spatial range, and adjacent connection relationships. This step further divides the dominant heat-controlling path into path segments capable of receiving residual attribution, allowing temperature measurement residuals to be mapped to different heat transfer control locations, such as fracture heat conduction, low thermal conductivity strata insulation, or basalt layer heat transfer disturbance. The research report indicates that the temperature field is controlled by fracture structures, and fracture heat conduction, basalt layer convective heat transfer, and low thermal conductivity strata insulation jointly influence local thermal anomalies. This path segment division is consistent with the aforementioned geothermal mechanism.

[0123] Based on the spatial proximity and stratigraphic affiliation between the measurement unit containing the temperature measurement residual record and the heat conduction, insulation, and heat transfer disturbance sections, the temperature measurement residual records are segmented and aggregated to generate temperature measurement residual attribution correction quantities. Specifically, for each valid temperature measurement residual record, the planar location, depth range, and stratigraphic affiliation of its measurement unit are read; for the heat conduction, insulation, and heat transfer disturbance sections in the dominant heat control path, the planar range, depth range, and stratigraphic affiliation of the path segment are read. First, stratigraphic matching is performed. If the measurement unit containing the residual record is located in the same stratigraphy, adjacent stratigraphic interfaces, or within the same fault connectivity range as a certain path segment, then that path segment is listed as a candidate aggregation object. Next, spatial proximity matching is performed. If the horizontal or vertical distance between the measurement unit containing the residual record and the candidate aggregation object is within the influence range of that path segment, then the residual record is assigned to the corresponding path segment. To make the aggregation rules transparent, the residual records are assumed to be... For path segments The aggregation weight is It is determined by both spatial proximity and stratigraphic affiliation; the closer the spatial distance, the more... When the size is larger, and the stratigraphic level is the same or within the same fracture connectivity range, Retained when there is no corresponding relationship between the layers. Set to zero. The specific calculation formula is as follows: , ,in, The path segment from the measurement unit where the temperature residual is located. Spatial distance, For the same measurement unit's planar dimensions or the center distance between adjacent measurement units, Distance weights , Layer-level attribution weight. Path segment. The temperature measurement residual attribution correction is denoted as The calculation formula is: ,in, To ensure effective recording of temperature measurement residuals, For this residual record, the path segment The aggregation weight, For path segment The attribution correction for the temperature measurement residual. If If the value is positive and the path segment type is a heat-conducting segment, then output the correction amount for the heat-conducting segment enhancement direction; if... If the value is positive and the path segment type is insulation segment, then output the correction amount for the reinforcement direction of the insulation segment; if... If the residuals are concentrated in the adjacent measurement units of the heat transfer disturbance section, the correction amount for the adjustment direction of the heat transfer disturbance section is output. If the number of effective residual records corresponding to a certain path segment is less than the preset minimum number of records, no correction amount is generated for that path segment, and a supplementary temperature measurement constraint mark is output. As a verification example, when the measured temperature at a certain well location is consistently higher than the initial predicted temperature in the depth segment near the fracture, and the measurement unit is within the same fracture connection range as the heat conduction segment in the dominant heat control path, the residual record is collected into the heat conduction segment and forms a positive temperature measurement residual attribution correction amount; when the high measured temperature value is concentrated in the depth segment under the low thermal conductivity formation, the residual record is collected into the insulation segment and forms an insulation constraint correction amount. This step outputs the temperature measurement residual attribution correction amounts corresponding to the heat conduction segment, insulation segment, and heat transfer disturbance segment, and passes them to the subsequent temperature field correction step to correct the heat conduction calculation conditions of the corresponding path segment, combined with Figure 3 and Figure 5 ,Will Figure 3 Measured temperature during drilling and Figure 5 The predicted temperature at the well location in the corresponding temperature field is correlated within the same well location and depth range to generate temperature measurement residual records. These residual records are then segmented and aggregated according to their spatial proximity and stratigraphic affiliation with the heat conduction, insulation, and heat transfer disturbance sections of the dominant heat control path. During aggregation, residual records within the same stratum are preferentially assigned to the corresponding path segment, while those on either side of adjacent strata are assigned secondary. Residual records lacking spatial and stratigraphic connectivity are not included in the path segment correction.

[0124] The present invention is further configured such that, in correcting the heat conduction calculation conditions of the corresponding path segment based on the temperature measurement residual attribution correction amount, and re-performing the heat conduction calculation to generate the corrected temperature field prediction result, the following steps are taken:

[0125] The temperature measurement residual attribution correction is associated with the heat-conducting segment, insulation segment, or heat transfer disturbance segment in the dominant heat-controlling path according to the path segment type, generating path segment correction data. Specifically, in the high-temperature geothermal multi-source measurement data fusion prediction system, the input objects include dominant heat-controlling path data, path segment data, initial heat conduction calculation conditions, and temperature measurement residual attribution correction. The dominant heat-controlling path data includes path number, path unit sequence, measurement unit identifier, and connection order; the path segment data includes three types of path segments: heat-conducting segment, insulation segment, and heat transfer disturbance segment. Among them, the heat-conducting segment originates from the fault penetration unit or fault influence zone unit, the insulation segment originates from the low thermal conductivity stratum unit, and the heat transfer disturbance segment originates from the basalt stratum unit or local heat transfer anomaly unit. The system first performs consistency verification on the path segment number, measurement unit number, and stratum attribution identifier, and then associates the temperature measurement residual attribution correction with the corresponding path segment according to the path segment number and path segment type to form path segment correction data. The path segment correction data must include at least the path number, path segment number, path segment type, associated measurement unit, temperature measurement residual attribution correction amount, residual direction, and the calculation conditions for the heat conduction to be corrected. The temperature measurement residual attribution correction amount is denoted as... ,in Indicates the corresponding path segment; The residuals are segmented and aggregated from the measured and initial predicted temperatures within the same or adjacent formations. Positive values ​​indicate that the measured temperature under the control of the corresponding path segment is higher than the initial predicted temperature, while negative values ​​indicate that the measured temperature is lower than the initial predicted temperature. The path segment correction data output in this step is used to clarify which type of path segment the residual correction should be applied to, avoiding the direct averaging of the temperature measurement residuals across the entire temperature field.

[0126] Based on the path segment correction data, the equivalent thermal conductivity or boundary heat flux input of the heat-conducting segment, the equivalent thermal conductivity or interlayer thermal resistance of the insulation segment, and the local heat transfer correction terms of the heat transfer disturbance segment are corrected to generate corrected heat conduction calculation conditions. Specifically, the system reads different heat conduction calculation conditions according to the path segment type. For the heat-conducting segment, the equivalent thermal conductivity and fracture boundary heat flux input are read; for the insulation segment, the equivalent thermal conductivity and interlayer thermal resistance are read; and for the heat transfer disturbance segment, the local heat transfer correction terms and the corresponding thermophysical property constraints of the measurement unit are read. To ensure the correction process has implementable boundaries, the path segment... The calculation conditions to be corrected are uniformly denoted as Among them, when the path segment is a heat-conducting segment, Corresponding equivalent thermal conductivity or boundary heat flux input; when the path segment is an insulation segment. Corresponding equivalent thermal conductivity or interlayer thermal resistance; when the path segment is a heat transfer disturbance segment. Corresponding to the local heat transfer correction term. The normalized residual is generated based on the temperature measurement residual attribution correction: ,in, For the bounded residual correction direction quantity, The reference temperature difference is determined jointly by the allowable error of the temperature measuring instrument, the initial temperature field fitting error, and the required survey accuracy; the hyperbolic tangent function is used to limit excessively large single-point residuals within a bounded range. The path segment calculation conditions are corrected according to the following formula: ,in, The calculation conditions for the path segment before correction. The calculation conditions for the corrected path segment. The allowable single correction range for this path segment is derived from the corresponding lithological thermophysical property test range, boundary heat flux value range, or existing temperature fitting results. and This represents the physical boundary of the calculation conditions for this path segment, derived from the measured thermal conductivity values ​​of the same lithology, the thermal resistance range of the same layer, the interpretation range of fracture heat flow, or the allowable range of heat transfer disturbance terms. In the heat-conducting segment, if... If positive, then increase the equivalent thermal conductivity or boundary heat flux input of the heat-conducting section; if... If the value is negative, then the corresponding calculation conditions should be lowered. In the insulation section, if... If the residual is positive and concentrated in the underlying unit of a low thermal conductivity stratum, then the interlayer thermal resistance increases or the equivalent thermal conductivity decreases; if A negative value reduces interlayer thermal resistance or increases equivalent thermal conductivity. In the heat transfer disturbance section, according to... The direction of the local heat transfer correction term is adjusted to align the local temperature field near the path segment with the direction of the measured residual. After the correction is completed, the direction of the local temperature field near the path segment is adjusted. Backfilling is performed to the corresponding measurement unit to generate corrected heat conduction calculation conditions. This step confines the residual correction to the corresponding path segment and corresponding physical conditions, ensuring that the correction results retain geological attribution relationships and thermal property boundaries.

[0127] Based on the corrected heat conduction calculation conditions, the heat conduction calculation is re-performed on the measurement unit set to generate the corrected predicted temperature for each measurement unit. Specifically, the system loads the corrected heat conduction calculation conditions into the measurement unit set, retaining the original planar position, depth range, and formation interface relationship of the measurement units, and only replacing the equivalent thermal conductivity, boundary heat flux input, interlayer thermal resistance, or local heat transfer correction terms corresponding to the path segment correction data. The heat conduction calculation adopts a steady-state discrete heat conduction calculation method, and the connection relationship between measurement units is determined by the shared planar boundary, shared depth boundary, and formation contact interface. The corrected temperature field is denoted as... The calculation relationship can be expressed as: ,in, This is the corrected equivalent thermal conductivity field. This represents the heat generation rate field corresponding to the heat generation constraints within the strata. This is the corrected predicted temperature field. For the heat conduction section, where boundary heat flow input correction is used, the heat flow conditions on the fracture boundary are synchronously replaced with the corrected boundary heat flow. For the insulation section, where interlayer thermal resistance correction is used, the vertical heat transfer connection conditions between adjacent measurement units are replaced with the corrected interlayer thermal resistance. For the heat transfer disturbance section, where local heat transfer correction is used, the local heat transfer or equivalent thermal conductivity correction result of the corresponding measurement unit is loaded into that section. After recalculation, the corrected predicted temperature of each measurement unit is output, and the calculation iteration status, boundary condition version, and corrected path segment number are retained. If the corrected residual at the same depth in the same well still exceeds the preset allowable range, the new residual record can be returned to the residual attribution step for one or more iterations. If the allowable range is reached or the preset number of iterations is reached, the corrected predicted temperature for this round is output. This step forms a closed-loop calculation process from the temperature measurement residual attribution correction amount to the corrected predicted temperature.

[0128] The system assembles the corrected predicted temperatures based on the planar location, depth range, and formation interface of each measurement unit to generate the corrected temperature field prediction results. Specifically, the system reads the corrected predicted temperatures, planar location identifiers, depth range identifiers, formation attribution identifiers, and corresponding corrected heat conduction calculation conditions for each measurement unit. The corrected predicted temperatures are then assembled according to the chronological order of planar location and the order of depth range from shallow to deep to form corrected temperature unit data. Formation interface continuity checks are performed on the corrected predicted temperatures of adjacent depth ranges at the same planar location; lateral continuity checks are performed on the corrected predicted temperatures of the same depth range at adjacent planar locations; and fault boundary records are retained for temperature differences on both sides of the fault penetration location. After assembly, the corrected temperature field prediction results are generated, including three categories: the corrected predicted temperatures of each measurement unit, corrected temperature profile data, and corrected temperature field spatial data. As a verification example, if the attribution correction for the temperature measurement residual corresponding to the heat conduction section of a certain dominant heat control path is positive, the recalculated well location predicted temperature after the bounded correction of the heat flow input at the boundary of the heat conduction section should be closer to the actual drilling temperature. If the residual of the underlying measurement unit of the low thermal conductivity formation is positive, the temperature rise of the lower part of the layer after the correction of the interlayer thermal resistance of the insulation section should vary within the allowable range of the thermal property boundary. This result is then passed to the target isotherm depth extraction step to generate the target isotherm depth prediction result and the high-temperature geothermal target area prediction result.

[0129] The present invention is further configured such that the step of extracting the minimum burial depth surface that reaches the target temperature from the corrected temperature field prediction results to generate the target isothermal surface burial depth prediction results, and merging the measurement units whose target isothermal surface burial depth meets the preset exploration depth conditions and is connected to the dominant heat control path to generate the high-temperature geothermal target area prediction results, includes:

[0130] The corrected temperature unit data is generated by extracting the corrected predicted temperature, planar location identifier, depth interval identifier, and stratigraphic affiliation identifier of each measurement unit from the corrected temperature field prediction results. Specifically, in the high-temperature geothermal multi-source measurement data fusion prediction system for complex structural areas, the corrected temperature field prediction results output by the temperature field correction module are used as input. This input includes the corrected predicted temperature, planar location identifier, depth interval identifier, stratigraphic affiliation identifier, associated heat-controlling dominant path number, and heat conduction calculation condition version for each measurement unit. During data extraction, the measurement units are first grouped according to their planar location identifiers, and then sorted from shallow to deep according to their depth intervals to form a vertical temperature sequence under the same planar location. For measurement units lacking corrected predicted temperatures, the temperature field calculation status of their adjacent measurement units at the same stratum is read. If the adjacent units have continuous stratigraphic interfaces and the same calculation version, they are marked as units to be supplemented; if they lack depth intervals or stratigraphic affiliation identifiers, they are marked as invalid units and are not included in the target isotherm extraction. After filtering and sorting, corrected temperature cell data is generated. This data includes at least the planar location number, upper and lower boundaries of the depth interval, corrected predicted temperature, formation affiliation, and the correlation status of the heat control path. This step organizes the three-dimensional temperature field results into a sequence of cell data that can be retrieved along the depth direction, enabling the subsequent target temperature arrival location to be extracted from a unified measurement cell structure.

[0131] The calibration temperature unit data is retrieved in ascending order from shallow to deep within the same planar location, and the measurement unit that first reaches the target temperature is designated as the target temperature arrival unit; specifically, the target temperature is set according to the exploration task. The target temperature can be set to 100℃, 150℃, or other target temperatures based on the high-temperature geothermal evaluation target. For the calibration temperature unit data at the same planar location, the calibration prediction temperature is retrieved sequentially from shallow to deep depth intervals. When the calibration prediction temperature of a measurement unit within a certain depth interval first reaches or exceeds the target temperature... Furthermore, the corrected predicted temperature for all shallower depth ranges at this planar location is lower than the target temperature. At that time, the measurement unit is recorded as the target temperature arrival unit. To make the arrival rule verifiable, the planar position can be... Next The corrected predicted temperature for each depth range is denoted as ,when and For all At the time of its establishment, the first The depth range represents the depth range at which the target temperature is reached. This formula is used to define the initial threshold for determining the temperature. Derived from the objectives of the exploration mission, It originates from the corrected temperature field prediction results. After the search is completed, the target temperature arrival cell is output, and the data includes the plane location number, target temperature, arrival depth range, corrected predicted temperature, and formation affiliation.

[0132] Isothermal surface depth points are generated based on the depth range of the target temperature arrival unit and the location of the formation interface. These isothermal surface depth points are then spatially aggregated to generate the target isothermal surface depth prediction result. Specifically, for each target temperature arrival unit, the upper and lower boundaries of its depth range and the location of adjacent formation interfaces are read. If both the upper and lower bound temperatures of the target temperature arrival unit can be determined by adjacent measurement units, linear interpolation is used to determine the depth point corresponding to the target temperature within that depth range. If only the unit center temperature is available within that depth range, the midpoint of that depth range or a preset representative depth is used as the isothermal surface depth point, and the source of the representative depth is recorded in the result. One implementation method for linear interpolation is as follows: ,in, Indicates planar position Target temperature The corresponding burial depth point, and These represent the representative depths of adjacent depth intervals. and These represent the corrected predicted temperatures at corresponding depths. This interpolation is only used when the temperature at adjacent depths increases with depth and the data is continuous; when a fracture interface causes a sudden change in local temperature, the location of the formation interface is retained as the boundary constraint for the isothermal surface depth points. All isothermal surface depth points corresponding to all planar locations are spatially aggregated according to their planar adjacency to form the target isothermal surface depth prediction result. This result includes the planar location, target temperature, isothermal surface depth, formation location, and data source status.

[0133] The predicted depth of the target isotherm is matched with preset exploration depth conditions to select measurement units whose depth falls within the preset exploration depth range, generating candidate target area units. Specifically, preset exploration depth conditions are set according to exploration engineering conditions, including the upper limit, lower limit, and workable depth range of the target isotherm. The system matches each isotherm depth point with the preset exploration depth conditions. When the target isotherm depth falls within the preset exploration depth range, the target temperature arrival unit corresponding to that plane position and its adjacent depth interval measurement units are included in the candidate target area units. When the target isotherm depth is shallower than the preset lower limit or deeper than the preset upper limit, the corresponding measurement unit is marked as not meeting the exploration depth conditions. The thresholds in this step are derived from exploration design conditions, drilling capabilities, and target reservoir temperature requirements, and are not arbitrarily generated by the model. After screening, the candidate target area unit data includes plane position, target isotherm depth, depth condition matching status, stratigraphic affiliation, and corresponding target temperature arrival unit.

[0134] Connectivity matching is performed between candidate target units and measurement units traversed by the dominant heat control path. Candidate target units connected to the dominant heat control path within adjacent planar locations, adjacent depth ranges, or the same stratigraphic interface are retained. Specifically, the measurement units traversed by the dominant heat control path and their path connection sequence are read, and spatial and stratigraphic connectivity matching is performed between candidate target units and measurement units along the dominant heat control path. If a candidate target unit shares a planar boundary, a depth boundary, or is located on both sides of the same stratigraphic interface with a measurement unit in the dominant heat control path, it is recorded as connected. If a candidate target unit is located within the same fault influence zone or the same fold core control area, it is also recorded as structurally connected. The connectivity matching results retain the connection type, including planar connectivity, vertical connectivity, stratigraphic interface connectivity, and structural connectivity. Candidate target units that do not form a connectivity relationship with any dominant heat control path are marked as isolated temperature units and are not included in target unit merging. This step constrains the temperature achievement results with the dominant heat control path, ensuring that the target unit simultaneously meets the temperature condition, depth condition, and heat control path connectivity condition.

[0135] The retained candidate target area units are merged according to their adjacency and stratigraphic continuity to generate target area connectivity unit groups. Specifically, unit merging is performed on candidate target area units that pass connectivity matching. During merging, laterally adjacent units are first identified based on shared planar boundaries, then vertically adjacent units are identified based on shared depth boundaries, and finally stratigraphic continuity units are identified based on the same stratum or the contact interface between adjacent strata. Candidate target area units that satisfy any adjacency relationship and whose target isothermal surface burial depth meets the preset exploration depth conditions are merged into the same target area connectivity unit group. During the merging process, the dominant heat control path number and connection type are retained. If multiple candidate target area units belong to different dominant heat control paths but form a connection within the same stratigraphic continuity area, they are recorded as a composite path target area connectivity unit group. Connectivity unit groups with an area or volume lower than the preset minimum target area size are marked as small-scale candidate areas; connectivity unit groups composed of multiple planar locations and multiple depth intervals and continuously connected to the dominant heat control path are retained as high-temperature geothermal target area output objects. The target area connectivity unit groups formed in this step have clear spatial boundaries, depth ranges, and sources of heat control paths.

[0136] High-temperature geothermal target area prediction results are generated based on the outer edge measurement units of the target area connected unit group, the corresponding target isothermal surface burial depth, and the associated dominant heat control path. Specifically, the system reads the outer edge measurement units of each target area connected unit group. The outer edge measurement units consist of measurement units adjacent to non-target area units in the connected unit group. Planar and depth boundaries are formed based on the outer edge measurement units. The target area burial depth range is formed based on the corresponding target isothermal surface burial depth. The target area heat control source record is formed based on the associated dominant heat control path. The high-temperature geothermal target area prediction results include target area number, planar range, burial depth range, target isothermal surface burial depth prediction results, associated dominant heat control path, main stratigraphic affiliation, and target area unit list. As a verification example, if the target isotherm depths of a group of candidate target area units are all within the preset exploration depth range, and these units form a continuous connection along fracture penetration units, low thermal conductivity formation units, and drilling temperature anomaly units, then this group of candidate target area units is merged into a single high-temperature geothermal target area prediction result. If another group of temperature-compliant units lacks a planar, depth, or formation interface connection with the dominant heat-controlling path, then this group of units is marked as an isolated temperature response area and is not output as a high-temperature geothermal target area. The high-temperature geothermal target area prediction result obtained through this step can be integrated with drilling deployment, supplementary measurement layout, and resource estimation modules to form a data closed loop from the corrected temperature field to the target area boundary output, such as... Figure 6 As shown, the final output prediction results for the high-temperature geothermal target area are jointly constrained by the target isothermal surface burial depth and the connectivity conditions of the dominant heat-controlling path. The system starts from... Figure 5 The minimum burial depth surface of the target temperature is extracted from the corresponding corrected temperature field. Measurement units that meet the preset exploration depth conditions are selected, and then connectivity matching is performed with the measurement units traversed by the dominant heat control path. The high-temperature geothermal target area prediction results generated after merging include the target area boundary, burial depth range, target isothermal surface, and associated heat control path. Its output results are consistent with... Figure 6 The output range of the distant scenic area is shown.

[0137] Example 2:

[0138] Please see Figure 2 This exemplary high-temperature geothermal multi-source measurement data fusion and prediction system for complex tectonic zones includes:

[0139] Measurement Unit Module: Acquires multi-source measurement data in complex structural areas, and categorizes the multi-source measurement data into measurement units according to planar location, depth range, and stratigraphic interface to generate a measurement unit set;

[0140] Anomaly Unit Marking Module: Marks anomalies in the measurement unit set to generate deep geophysical anomaly units, thermal control structure units, and drilling temperature anomaly units;

[0141] Path generation module: Starting from the deep geophysical anomaly unit, it connects segment by segment along the connection direction of the adjacent measurement unit through the heat-controlling structural unit to the drilling temperature anomaly unit to generate candidate thermal anomaly transmission paths.

[0142] Path filtering module: Calculates the consistency of thermal anomaly paths based on the continuous distribution of anomaly units in candidate thermal anomaly propagation paths, the thermal resistance of adjacent measurement units, and the interval between non-responding measurement units, and generates the dominant thermal control path according to the consistency of thermal anomaly paths.

[0143] Temperature field prediction module: Based on the thermophysical property data and radioactive heat generation data in the measurement unit corresponding to the dominant heat control path, heat conduction calculation is performed to generate initial temperature field prediction results;

[0144] Residual Attribution Module: Calculates the difference between the measured drilling temperature and the initial temperature field prediction results in the same well location and depth range to generate temperature measurement residuals; Based on the distance relationship and layer attribution relationship between the measurement unit where the temperature measurement residuals are located and different path segments in the dominant heat control path, the temperature measurement residuals are segmented and aggregated to generate temperature measurement residual attribution correction quantities.

[0145] Temperature field correction module: Corrects the heat conduction calculation conditions of the corresponding path segment based on the temperature measurement residual attribution correction amount, and re-performs heat conduction calculation to generate the corrected temperature field prediction results;

[0146] Target area generation module: Extracts the minimum burial depth surface that reaches the target temperature from the corrected temperature field prediction results to generate the target isothermal surface burial depth prediction results. Merges the measurement units whose target isothermal surface burial depth meets the preset exploration depth conditions and is connected to the heat control dominant path to generate high temperature geothermal target area prediction results.

[0147] It should be noted that the high-temperature geothermal multi-source measurement data fusion and prediction system for complex structural areas provided in the above embodiments and the high-temperature geothermal multi-source measurement data fusion and prediction method for complex structural areas provided in the above embodiments belong to the same concept. The specific operation methods of each module and unit have been described in detail in the method embodiments and will not be repeated here. In practical applications, the high-temperature geothermal multi-source measurement data fusion and prediction system for complex structural areas provided in the above embodiments can be assigned to different functional modules as needed, that is, the internal structure of the system can be divided into different functional modules to complete all or part of the functions described above. This is not a limitation here.

[0148] The above description is merely a specific embodiment of this application, but the scope of protection of this application is not limited thereto. Any variations or substitutions that can be easily conceived by those skilled in the art within the scope of the technology disclosed in this application should be included within the scope of protection of this application. Therefore, the scope of protection of this application should be determined by the scope of the claims.

Claims

1. A method for fusion and prediction of high-temperature geothermal multi-source measurement data in complex tectonic zones, characterized in that, include: Acquire multi-source measurement data in complex structural areas, and categorize the multi-source measurement data into measurement units according to planar location, depth range, and stratigraphic interface to generate a measurement unit set; Anomaly markers are applied to the measurement unit set to generate deep geophysical anomaly units, thermal control structure units, and drilling temperature anomaly units; Starting from the deep geophysical anomaly unit, candidate thermal anomaly transmission paths are generated by connecting the adjacent measurement units through the thermal control structure unit segment by segment to the drilling temperature anomaly unit. The consistency of the thermal anomaly path is calculated based on the continuous distribution of anomaly units in the candidate thermal anomaly propagation path, the thermal resistance of adjacent measurement units, and the interval of non-responding measurement units. The thermal anomaly path is then selected and generated according to the consistency of the thermal anomaly path. The initial temperature field prediction results are generated by performing heat conduction calculations based on the thermophysical property data and radioactive heat generation data in the measurement unit corresponding to the dominant heat control path. The temperature measurement residual is generated by calculating the difference between the actual drilling temperature and the initial temperature field prediction results in the same well location and the same depth range. The temperature measurement residual is then segmented and aggregated according to the distance relationship between the measurement unit where the temperature measurement residual is located and different path segments in the heat control dominant path and the layer attribution relationship to generate the temperature measurement residual attribution correction amount. Based on the temperature measurement residual attribution correction, the heat conduction calculation conditions of the corresponding path segment are corrected, and the heat conduction calculation is re-performed to generate the corrected temperature field prediction results. Extract the minimum burial depth surface that reaches the target temperature from the corrected temperature field prediction results to generate the target isothermal surface burial depth prediction results. Merge the measurement units whose target isothermal surface burial depth meets the preset exploration depth conditions and are connected to the dominant heat control path to generate the high temperature geothermal target area prediction results.

2. The method for fusion and prediction of high-temperature geothermal multi-source measurement data in complex tectonic zones according to claim 1, characterized in that, Acquire multi-source measurement data in complex structural areas, and then group the multi-source measurement data into measurement units according to planar location, depth range, and stratigraphic interface to generate a measurement unit set, including: Multi-source measurement data for complex tectonic zones include geological structural data, gravity measurement data, magnetotelluric sounding data, well temperature measurement data, rock thermal property data, and radioactive heat generation data; Stratigraphic constraint subdivision data is established based on stratigraphic interfaces, fault boundaries, and well-exposed strata in geological structural data; planar subdivision data is established based on the location of measuring points, measuring lines, and well locations. The planar subdivision data and the stratigraphically constrained subdivision data are matched, and measurement units are generated according to the depth range defined by the planar subdivision range and the adjacent stratigraphic interfaces. Gravity measurement data and magnetotelluric sounding data are categorized into corresponding measurement units according to measurement point location, survey line coverage, inversion depth, and detection coverage to generate geophysical data items; well temperature measurement data are categorized into corresponding measurement units according to well location and temperature measurement depth to generate temperature data items; rock thermal property data are categorized into corresponding measurement units according to sample collection location, sample collection depth, and the strata to which the sample belongs to generate thermal property data items; radioactive heat generation data are categorized into corresponding measurement units according to sample collection location, sample collection depth, and the strata to which the sample belongs to generate heat generation data items; geological structural data are categorized into corresponding measurement units according to stratigraphic interfaces, fault boundaries, lithological contact interfaces, and well-exposed strata to generate geological attribute data items. A measurement unit set is generated based on the measurement units containing geophysical data items, temperature data items, thermal property data items, heat generation data items, and geological attribute data items.

3. The method for fusion and prediction of high-temperature geothermal multi-source measurement data in complex tectonic zones according to claim 2, characterized in that, Anomaly markers are generated from the measurement unit set to create deep geophysical anomaly units, thermal control structure units, and drilling temperature anomaly units, including: Extract exploration data items, geological attribute data items, and temperature data items from the measurement unit; Geophysical data items are layered according to the detection depth range or inversion depth layer, and measurement units corresponding to changes in deep density or deep electrical interface are identified. The corresponding measurement units are marked as deep geophysical anomaly units. The geological attribute data items are used to identify thermally controlling geological elements. Measurement units containing fault cuts, fold cores, basalt layers, low thermal conductivity strata, or lithological contact interfaces are marked as thermally controlling structural units. The drilling temperature curve is segmented and identified for the temperature data items. The measurement units containing the corresponding high temperature segments or temperature gradient abrupt change segments are marked as drilling temperature anomaly units.

4. The method for fusion and prediction of high-temperature geothermal multi-source measurement data in complex tectonic zones according to claim 1, characterized in that, Starting from a deep geophysical anomaly unit, candidate thermal anomaly propagation paths are generated by connecting segment by segment to the drilling temperature anomaly unit along the connection direction of adjacent measurement units through the heat-controlling structural unit. These paths include: The deep geophysical exploration anomaly unit is used as the starting unit of the path, and the drilling temperature anomaly unit is extracted as the ending unit of the path. Establish an adjacent unit connection table based on the shared plane boundary, shared depth boundary and formation contact interface between measurement units; Starting from the path starting unit, the adjacent measurement units with heat control structure markers or drilling temperature anomaly markers are connected segment by segment along the adjacent unit connection table to generate a unit connection sequence; When a unit connection sequence passes through at least one thermal control construction unit and connects to a path termination unit, the unit connection sequence is recorded as a candidate thermal anomaly propagation path.

5. The method for fusion and prediction of high-temperature geothermal multi-source measurement data in complex tectonic zones according to claim 4, characterized in that, The consistency of the thermal anomaly path is calculated based on the continuous distribution of anomaly units in the candidate thermal anomaly propagation path, the thermal resistance of adjacent measurement units, and the interval between non-responding measurement units. The dominant thermal control path is then selected and generated according to this consistency, including: Generate a path unit sequence according to the path connection order in the candidate thermal anomaly propagation path; Abnormal continuous segment data is generated by merging the abnormal markers of measurement units in the path unit sequence; Path thermal resistance data is generated based on the thermal property data and depth interval of adjacent measurement units in the path unit sequence, and the interval data of non-response measurement units between adjacent abnormal continuous segments is identified. The thermal anomaly path consistency is generated based on the abnormal continuous segment data, path thermal resistance data, and non-response measurement unit interval data. Candidate thermal anomaly propagation paths are then selected based on the thermal anomaly path consistency to generate the dominant thermal control path.

6. The method for fusion and prediction of high-temperature geothermal multi-source measurement data in complex tectonic zones according to claim 3, characterized in that, The initial temperature field prediction results are generated based on the thermophysical property data and radioactive heat generation data in the measurement unit corresponding to the dominant heat control path, including: Extract the measurement units along the dominant heat control path to generate a sequence of temperature field calculation units; The thermal property data items, radioactive heat generation data items, stratigraphic attribution identifiers, and heat-controlling structure markers are read from the temperature field calculation unit sequence to generate path heat conduction constraint data; Based on the path heat conduction constraint data, the fracture penetration unit is configured as a heat conduction channel constraint, the low thermal conductivity stratum unit is configured as a heat insulation constraint, the basalt layer unit is configured as a local heat transfer disturbance constraint, and the radioactive heat generation data item is configured as a layer internal heat generation constraint. Heat conduction calculations are performed based on constraints such as heat conduction channels, insulation, local heat transfer disturbances, and internal heat generation within the layer to generate initial temperature field prediction results.

7. The method for fusion and prediction of high-temperature geothermal multi-source measurement data in complex tectonic zones according to claim 5, characterized in that, The temperature measurement residuals are calculated and generated. The temperature measurement residuals are then segmented and aggregated to generate the temperature measurement residual attribution correction amount, which includes: Extract the well location predicted temperature data corresponding to the drilling temperature measurement data from the initial temperature field prediction results, and match the well location predicted temperature data with the actual drilling temperature according to the drilling number and depth range to generate well temperature corresponding data. The difference between the measured drilling temperature and the predicted well location temperature for the same drilling number and the same depth range in the well temperature data is processed to generate a temperature measurement residual record. The heat control dominant path is divided into a heat conduction section, a heat insulation section, and a heat transfer disturbance section according to the heat control structure markings; Based on the spatial proximity and layer-level affiliation between the measurement unit where the temperature measurement residual record is located and the heat conduction section, heat preservation section, and heat transfer disturbance section, the temperature measurement residual record is segmented and aggregated to generate the temperature measurement residual attribution correction amount.

8. The method for fusion and prediction of high-temperature geothermal multi-source measurement data in complex tectonic zones according to claim 5, characterized in that, Based on the attribution correction of the temperature measurement residual, the heat conduction calculation conditions for the corresponding path segment are corrected, and the heat conduction calculation is re-performed to generate the corrected temperature field prediction results, including: The temperature measurement residual attribution correction amount is associated with the heat conduction section, heat preservation section or heat transfer disturbance section in the heat control dominant path according to the path segment type, and the path segment correction data is generated. Based on the path segment correction data, the equivalent thermal conductivity or boundary heat flow input of the heat conduction segment, the equivalent thermal conductivity or interlayer thermal resistance of the insulation segment, and the local heat transfer correction terms of the heat transfer disturbance segment are corrected to generate the corrected heat conduction calculation conditions. Based on the corrected heat conduction calculation conditions, the heat conduction of the measurement unit set is recalculated to generate the corrected predicted temperature of each measurement unit; The corrected temperature field prediction results are generated by assembling the corrected predicted temperatures according to the planar location, depth range, and formation interface of each measurement unit.

9. The method for fusion and prediction of high-temperature geothermal multi-source measurement data in complex tectonic zones according to claim 8, characterized in that, The minimum burial depth surface reaching the target temperature is extracted from the corrected temperature field prediction results to generate the target isothermal surface burial depth prediction results. Measurement units whose target isothermal surface burial depth meets the preset exploration depth conditions and are connected to the dominant heat control path are merged to generate high-temperature geothermal target area prediction results, including: The corrected predicted temperature, plane location identifier, depth interval identifier, and stratigraphic affiliation identifier of each measurement unit are extracted from the corrected temperature field prediction results to generate corrected temperature unit data; The calibration temperature unit data is retrieved in order from shallow to deep within the same plane position, and the measurement unit that first reaches the target temperature is taken as the target temperature arrival unit. Based on the depth range of the target temperature reaching the unit and the location of the formation interface, isothermal surface burial points are generated. The isothermal surface burial points corresponding to each plane location are spatially aggregated to generate the target isothermal surface burial depth prediction results. The predicted depth of the target isotherm surface is matched with the preset exploration depth conditions, and the measurement units whose depth of the target isotherm surface is within the preset exploration depth range are selected to generate candidate target area units. Connectivity matching is performed between candidate target units and measurement units traversed by the dominant heat control path, and candidate target units that are connected to the dominant heat control path in adjacent plane positions, adjacent depth intervals, or within the same formation contact interface are retained. The retained candidate target area units are merged according to their adjacency and stratigraphic continuity to generate target area connectivity unit groups; The prediction results for the high-temperature geothermal target area are generated based on the outer edge measurement units of the target area connectivity unit group, the burial depth of the corresponding target isothermal surface, and the associated dominant heat control path.

10. A system for fusing and predicting high-temperature geothermal multi-source measurement data in complex tectonic zones, used to implement the method for fusing and predicting high-temperature geothermal multi-source measurement data in complex tectonic zones as described in any one of claims 1-9, characterized in that, include: Measurement Unit Module: Acquires multi-source measurement data in complex structural areas, and categorizes the multi-source measurement data into measurement units according to planar location, depth range, and stratigraphic interface to generate a measurement unit set; Anomaly Unit Marking Module: Marks anomalies in the measurement unit set to generate deep geophysical anomaly units, thermal control structure units, and drilling temperature anomaly units; Path generation module: Starting from the deep geophysical anomaly unit, it connects segment by segment along the connection direction of the adjacent measurement unit through the heat-controlling structural unit to the drilling temperature anomaly unit to generate candidate thermal anomaly transmission paths. Path filtering module: Calculates the consistency of thermal anomaly paths based on the continuous distribution of anomaly units in candidate thermal anomaly propagation paths, the thermal resistance of adjacent measurement units, and the interval between non-responding measurement units, and generates the dominant thermal control path according to the consistency of thermal anomaly paths. Temperature field prediction module: Based on the thermophysical property data and radioactive heat generation data in the measurement unit corresponding to the dominant heat control path, heat conduction calculation is performed to generate initial temperature field prediction results; Residual Attribution Module: Calculates the difference between the measured drilling temperature and the initial temperature field prediction results in the same well location and depth range to generate temperature measurement residuals; Based on the distance relationship and layer attribution relationship between the measurement unit where the temperature measurement residuals are located and different path segments in the dominant heat control path, the temperature measurement residuals are segmented and aggregated to generate temperature measurement residual attribution correction quantities. Temperature field correction module: Corrects the heat conduction calculation conditions of the corresponding path segment based on the temperature measurement residual attribution correction amount, and re-performs heat conduction calculation to generate the corrected temperature field prediction results; Target area generation module: Extracts the minimum burial depth surface that reaches the target temperature from the corrected temperature field prediction results to generate the target isothermal surface burial depth prediction results. Merges the measurement units whose target isothermal surface burial depth meets the preset exploration depth conditions and is connected to the heat control dominant path to generate high temperature geothermal target area prediction results.

Images