Geothermal resource evaluation method and system

Abstract

According to the geothermal resource evaluation method and system, on the basis of an analytic hierarchy process, geothermal elements such as a heat source, a heat storage layer, a cover layer, a water source and a channel serve as criterion layers, and a set of geothermal resource potential evaluation system used for layered heat storage and strip heat storage is established; according to the evaluation system provided by the invention, the geothermal resource reservoir forming conditions are systematically analyzed and quantified, a geothermal resource potential evaluation map is formed, and geothermal exploration and development planning blocks and exploration and development prospects are determined.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of computer, in particular to a geothermal resource evaluation method and system. BACKGROUND

[0002] At present, the geothermal resource is related to the geothermal flow value, the geothermal gradient, the surface heat display, the buried depth of the Moho surface, the fracture cutting depth, the fracture extension length, the fracture activity period, the earthquake activity period, the magnitude, the stratum lithology, the fracture density, the plane distance from the supply area, the caprock thickness, the heat reservoir thickness, the fracture width, the heat reservoir lithology, the heat reservoir temperature, the infiltration coefficient, and the annual average precipitation. The geothermal resource can be used for urban heating, hot spring tourism and physiotherapy, and bathing. However, how to establish an evaluation system to evaluate the geothermal resource at each position is a problem. SUMMARY

[0003] The present application aims to provide a geothermal resource evaluation method and system to solve the above problems in the prior art.

[0004] In a first aspect, the present application provides a geothermal resource evaluation method, comprising: obtaining a geothermal resource data set at a plurality of positions; the geothermal resource data set comprises a plurality of geothermal resource data; constructing evaluation indexes for the plurality of geothermal resource data; a geothermal evaluation index system is established by a layered heat reservoir and a strip-shaped heat reservoir; the geothermal evaluation index system comprises a target layer, a criterion layer, a first index layer, and a second index layer; the second index layer comprises a layered heat reservoir evaluation index and a strip-shaped heat reservoir evaluation index; detecting the importance proportion of the evaluation indexes of adjacent levels in the geothermal evaluation index system to obtain a first judgment matrix, a second judgment matrix, and a third judgment matrix; detecting the weights of the plurality of evaluation indexes of the second index layer in the geothermal evaluation index system based on the first judgment matrix, the second judgment matrix, and the third judgment matrix to obtain an index weight set; the index weight set comprises a plurality of index weights; one layered heat reservoir index or strip-shaped heat reservoir evaluation index corresponds to one index weight; obtaining a geothermal resource potential evaluation map based on the geothermal resource data set and the index weight set; the geothermal resource potential evaluation map is used to represent the development potential of the geothermal resource in different regions.

[0005] Optionally, the detection of the weights of the plurality of evaluation indexes of the second index layer in the geothermal evaluation index system based on the first judgment matrix, the second judgment matrix, and the third judgment matrix to obtain the index weight set comprises: The first geothermal index weight vector is obtained based on judging the relationship between the evaluation indexes through the first judgment matrix; and the value in the first geothermal index weight vector represents the weight of the evaluation indexes in the criterion layer corresponding to the evaluation indexes in the target layer. The second judgment matrix corresponds to obtain a second geothermal index weight vector; the third judgment matrix corresponds to obtain a third geothermal index weight vector; the second geothermal index weight vector represents the weight of the evaluation indexes in the first index layer corresponding to the evaluation indexes in the criterion layer; and the third geothermal index weight vector represents the weight of the evaluation indexes in the second index layer corresponding to the evaluation indexes in the first index layer. Based on the first geothermal index weight vector, the second geothermal index weight vector and the third geothermal index weight vector, and the accuracy of the detected weight between the evaluation indexes of the geothermal resource potential, the first index weight accurate value, the second index weight accurate value and the third index weight accurate value are obtained. The first index weight accurate value, the second index weight accurate value and the third index weight accurate value are judged to obtain the index weight set.

[0006] Optionally, the importance proportion of the evaluation indexes of the adjacent levels in the geothermal evaluation index system is detected to obtain the first judgment matrix, the second judgment matrix and the third judgment matrix, which comprises: The multiple evaluation indexes in the criterion layer corresponding to one target index in the target layer are taken as the row and column to construct the first initial judgment matrix; and all the values in the first initial judgment matrix are 0. The importance of the evaluation indexes in the target layer and the multiple evaluation indexes in the criterion layer is respectively judged to fill in the corresponding positions of the first initial judgment matrix to obtain the first judgment matrix. The n second judgment matrices are obtained corresponding to the criterion layer and the first index layer; and the m third judgment matrices are obtained corresponding to the first index layer and the second index layer; n represents the number of the evaluation indexes in the first index layer; and m represents the number of the evaluation indexes in the second index layer.

[0007] Optionally, the geothermal resource potential evaluation map is obtained based on the geothermal resource data set and the index weight set, which comprises: The multiple geothermal resource scores are obtained based on the geothermal resource data set; one geothermal resource score corresponds to one geothermal resource data. The geothermal resource potential value is obtained by multiplying and adding the multiple geothermal resource scores and the corresponding index weight; the geothermal resource potential value represents the development potential of one position; and the multiple positions correspond to obtain the multiple geothermal resource potential values. The geothermal resource potential evaluation map is obtained by marking the multiple geothermal resource potential values in the multiple positions through gridding; and the geothermal resource potential evaluation map represents the development potential of each position.

[0008] Optionally, the second determination of the accuracy of the weights detected between the evaluation indicators of geothermal resource potential based on the first geothermal indicator weight vector, the second geothermal indicator weight vector, and the third geothermal indicator weight vector to obtain the first indicator weight accurate value, the second indicator weight accurate value, and the third indicator weight accurate value includes: Multiplying the first judgment matrix by the corresponding first geothermal index weight vector to obtain a first maximum eigenvalue vector; Find the value that is equal to the first maximum eigenvalue after multiplying the first geothermal index weight vector as the first maximum eigenvalue; Based on the first maximum eigenvalue and the first judgment matrix, the accuracy of the weights detected between the evaluation indicators of geothermal resource potential is secondarily judged to obtain the accurate value of the first indicator weight; The second judgment matrix corresponds to obtaining the accurate value of the second indicator weight; the third judgment matrix corresponds to obtaining the accurate value of the third indicator weight.

[0009] Optionally, the step of obtaining a plurality of geothermal resource scores based on the geothermal resource data set includes: Obtaining a geothermal standard table; the geothermal standard table includes standard scoring values ​​for multiple locations; the standard scoring values ​​represent scores corresponding to values ​​of multiple evaluation indicators of the second indicator layer in different ranges; The scores corresponding to the geothermal resource data are found in the geothermal standard set to obtain a geothermal resource score set; the geothermal resource score set includes multiple geothermal resource scores; the geothermal resource score represents the score of the geothermal resource corresponding to the layered heat storage index or the strip heat storage evaluation index.

[0010] Optionally, determining the accurate value of the first indicator weight, the accurate value of the second indicator weight, and the accurate value of the third indicator weight to obtain an indicator weight set includes: If the first indicator weight accuracy value is less than the geothermal relationship judgment threshold, and the second indicator weight accuracy value is less than the geothermal relationship judgment threshold, and the third indicator weight accuracy value is less than the geothermal relationship judgment threshold, it is assumed that the geothermal relationship judgment is accurate; If the geothermal relationship is accurately determined, the value in the third geothermal index weight vector is multiplied by the value in the corresponding second geothermal index weight vector to obtain an index weight set; the index weight set includes multiple index weights; the number of elements in the index weight set is equal to the number of evaluation indicators of the second index layer; If the geothermal relationship judgment is inaccurate, the first judgment matrix, the second judgment matrix, and the third judgment matrix are adjusted.

[0011] Optionally, the secondarily judging the accuracy of the weights detected between the evaluation indicators of geothermal resource potential based on the first maximum eigenvalue and the first judgment matrix to obtain the accurate value of the first indicator weight includes: Wherein, the first index weight accurate value is the first index weight accurate value, the first index weight accurate value is obtained by detecting the importance proportion of the evaluation index of the adjacent level in the geothermal evaluation index system, and the first index weight accurate value is obtained by detecting the weight of the evaluation index of the second index layer in the geothermal evaluation index system based on the first judgment matrix, the second judgment matrix and the third judgment matrix. The first maximum eigenvalue is represented by λmax, and the first minimum eigenvalue is represented by λmin. The number of values of the criterion layer is represented by n. The average random consistency index of the criterion layer is represented by CR.

[0012] Optionally, the relationship of the evaluation index is judged to obtain a first geothermal index weight vector, including: The row vectors of the first judgment matrix are geometrically averaged to obtain a plurality of first average values. The plurality of first average values are normalized to obtain the first geothermal index weight vector.

[0013] In a second aspect, an embodiment of the present application provides a geothermal resource evaluation system, including: An acquisition module is configured to acquire a plurality of geothermal resource data sets at different locations, wherein each geothermal resource data set includes a plurality of geothermal resource data. A hierarchical module is configured to construct evaluation indexes for the plurality of geothermal resource data, wherein a geothermal evaluation index system is established based on both stratiform thermal reservoirs and zonal thermal reservoirs, the geothermal evaluation index system includes a target layer, a criterion layer, a first index layer and a second index layer, and the second index layer includes stratiform thermal reservoir evaluation indexes and zonal thermal reservoir evaluation indexes. An association module is configured to detect the importance proportions of the evaluation indexes of adjacent levels in the geothermal evaluation index system to obtain a first judgment matrix, a second judgment matrix and a third judgment matrix. A weight module is configured to detect the weights of the plurality of evaluation indexes of the second index layer in the geothermal evaluation index system based on the first judgment matrix, the second judgment matrix and the third judgment matrix to obtain an index weight set, wherein the index weight set includes a plurality of index weights, one stratiform thermal reservoir index or zonal thermal reservoir evaluation index corresponds to one index weight. A potential module is configured to obtain a geothermal resource potential evaluation map based on the geothermal resource data set and the index weight set, wherein the geothermal resource potential evaluation map is used to represent the development potential of geothermal resources in different regions.

[0014] Compared with the prior art, the embodiment of the present application has the following beneficial effects: ​Embodiments of the present invention also provide a geothermal resource evaluation method and system. Based on geological characteristics, a geothermal evaluation index system can be used to determine the potential of geothermal areas in basins and peri-basin geothermal areas dominated by layered thermal reservoirs. The geothermal evaluation index system can identify basins dominated by layered thermal reservoirs and peri-basin geothermal areas with better potential. The Triassic thermal reservoir has slightly higher potential than the Permian thermal reservoir. Analysis suggests that this is related to the shallow burial depth and thickness of the Triassic thermal reservoir. Geothermal areas dominated by peri-zonal thermal reservoirs are generally better, reflecting the structural control of peri-zonal thermal reservoirs. The high-potential areas are consistent with actual understanding. The maximum score within the area is consistent with the rich geothermal resources in the above-mentioned areas. The distribution of existing mining rights and geothermal hotspots is consistent with the distribution of areas with medium to high potential. Secondary faults have a significant control effect on the replenishment, migration, and discharge of geothermal resources, which is very beneficial to geothermal formation. Generally, evaluation results with scores of medium or above are more consistent with reality. The geothermal resource potential evaluation system and typical geothermal hotspot evaluation can be verified through the geothermal evaluation index system. The evaluation results are consistent with the understanding of geothermal resource regularity and are consistent with the distribution and burial of geothermal resources. They can serve as a reference for subsequent geothermal resource evaluation. Based on the above potential evaluation results, comprehensive geothermal resource exploration in some areas can be considered in the future.

[0015] It can achieve the technical effect of scientifically supporting geothermal resource exploration and development by defining geothermal exploration and development planning blocks and exploration and development prospect areas based on evaluation results and combining them with plans at all levels. BRIEF DESCRIPTION OF THE DRAWINGS

[0016] Figure 1 This is a flow chart of a geothermal resource evaluation method provided by an embodiment of the present invention.

[0017] Figure 2 This is a framework diagram of a geothermal resource potential evaluation index system in a geothermal resource evaluation method provided by an embodiment of the present invention. DETAILED DESCRIPTION

[0018] The present invention will be described in detail below with reference to the accompanying drawings.

[0019] Example 1 like Figure 1 As shown, an embodiment of the present invention provides a geothermal resource evaluation method, the method comprising: S101: Obtain a set of geothermal resource data of multiple positions; the set of geothermal resource data comprises multiple geothermal resource data. In the embodiment, water-thermal geothermal resources in a region are taken as evaluation objects, and geothermal resource data of 15 positions are obtained. The geothermal resource data comprises specific values of terrestrial heat flow value, geothermal gradient, surface heat display, buried depth of living space, buried depth of Moho surface, fracture cutting depth, fracture extension length, fracture activity period, seismic activity period, magnitude, stratum lithology, fracture density, planar distance from recharge area, caprock thickness, heat reservoir thickness, fracture width, heat reservoir lithology, heat reservoir temperature, infiltration coefficient, and multi-year average precipitation. The geothermal resource data is data related to heat source reaction.

[0020] S102: Construct evaluation indexes for the multiple geothermal resource data; a geothermal evaluation index system is constructed for both layered heat reservoirs and zonal heat reservoirs; the geothermal evaluation index system comprises a target layer, a criterion layer, a first index layer, and a second index layer; the second index layer comprises layered heat reservoir evaluation indexes and zonal heat reservoir evaluation indexes.

[0021] In the embodiment, the layered heat reservoir evaluation indexes comprise terrestrial heat flow value (C1), geothermal gradient (C2), surface heat display (C3), buried depth of living space (C4), buried depth of Moho surface (C5), fracture cutting depth (C6), fracture extension length (C7), fracture activity period (C8), seismic activity period (C9), magnitude (C10), stratum lithology (C11), fracture density (C12), planar distance from recharge area (C13), caprock thickness (C14), heat reservoir thickness (C15), heat reservoir lithology (C17), heat reservoir temperature (C18), infiltration coefficient (C19), and multi-year average precipitation (C20); the zonal heat reservoir evaluation indexes comprise terrestrial heat flow value (C1), surface heat display (C3), buried depth of living space (C4), buried depth of Moho surface (C5), fracture cutting depth (C6), fracture extension length (C7), fracture activity period (C8), seismic activity period (C9), magnitude (C10), stratum lithology (C11), fracture density (C12), planar distance from recharge area (C13), caprock thickness (C14), fracture width (C16), heat reservoir lithology (C17), heat reservoir temperature (C18), infiltration coefficient (C19), and multi-year average precipitation (C20).

[0022] Wherein, as in this embodiment, the layered heat reservoir is a widely distributed Permian and Triassic carbonate rock heat reservoir, and the basin periphery and basin geothermal area dominated by the layered heat reservoir is evaluated using the layered heat reservoir evaluation index. The geothermal area dominated by the zonal heat reservoir is evaluated using the zonal heat reservoir evaluation index. In this embodiment, five criterion layers (A1-A5), ten secondary evaluation indexes (B1-B10), and twenty tertiary evaluation indexes (C1-C20) are selected. The evaluation index of the criterion layer is the first-level evaluation index, the evaluation index of the first index layer is the secondary evaluation index, and the evaluation index of the second index layer is the tertiary evaluation index. The tertiary evaluation index is a refinement of the secondary evaluation index.

[0023] Wherein, the evaluation index of the target layer is the geothermal resource potential. In combination with the geothermal resource accumulation conditions, five geothermal elements, i.e., heat source, heat reservoir, cap rock, water source, and channel, are taken as the criterion layer, so as to judge the geothermal resource potential based on the five criterion groups. In this embodiment, ten secondary indexes, i.e., geothermal field (B1), geophysical characteristics (B2), active faults (B3), seismic activity (B4), stratigraphic structure (B5), runoff distance (B6), cap rock thickness (B7), heat reservoir structure (B8), heat reservoir temperature (B9), and groundwater recharge (B10), are obtained. In this embodiment, around the ten secondary evaluation indexes, twenty tertiary indexes, i.e., terrestrial heat flow value (C1), geothermal gradient (C2), surface heat display (C3), crustal depth (C4), Moho depth (C5), fault cutting depth (C6), fault extension length (C7), fault activity period (C8), seismic activity period (C9), magnitude (C10), stratigraphic lithology (C11), fault density (C12), planar distance from recharge area (C13), cap rock thickness (C14), heat reservoir thickness (C15), fault width (C16), heat reservoir lithology (C17), heat reservoir temperature (C18), infiltration coefficient (C19), and annual average precipitation (C20), are further divided.

[0024] Wherein, the geothermal evaluation index system is shown in Table 1: Table 1 S103: The importance proportion of the evaluation indexes of the adjacent levels in the geothermal evaluation index system is detected to obtain a first judgment matrix, a second judgment matrix, and a third judgment matrix. Wherein, the first judgment matrix represents the importance proportion of the evaluation indexes of the criterion layer to the evaluation index of the target layer. The second judgment matrix represents the importance proportion of the evaluation indexes of the first index layer to the evaluation indexes of the corresponding criterion layer. The third judgment matrix represents the importance proportion of the evaluation indexes of the second index layer to the evaluation indexes of the corresponding first index layer.

[0025] S104: Based on the first judgment matrix, the second judgment matrix and the third judgment matrix, the weights of the multiple evaluation indexes of the second index layer in the geothermal evaluation index system are detected to obtain an index weight set; the index weight set includes multiple index weights; one stratiform heat reservoir index or one zonal heat reservoir evaluation index corresponds to one index weight.

[0026] S105: Based on the geothermal resource data set and the index weight set, a geothermal resource potential evaluation map is obtained; the geothermal resource potential evaluation map is used to represent the geothermal resource development potential of different regions.

[0027] Optionally, based on the first judgment matrix, the second judgment matrix and the third judgment matrix, the weights of the multiple evaluation indexes of the second index layer in the geothermal evaluation index system are detected to obtain an index weight set, including: Based on the first judgment matrix, the relationship between the evaluation indexes is determined to obtain a first geothermal index weight vector; the value in the first geothermal index weight vector represents the weight of the evaluation index corresponding to the target layer of the multiple evaluation indexes of the criterion layer. In this embodiment, the first judgment matrix is Z corresponding to Table 2. The first geothermal index weight vector is the normalized feature vector corresponding to Table 2.

[0028] The second judgment matrix corresponds to obtain a second geothermal index weight vector; the third judgment matrix corresponds to obtain a third geothermal index weight vector; the second geothermal index weight vector represents the weight of the evaluation index corresponding to the criterion layer of the multiple evaluation indexes of the first index layer; the third geothermal index weight vector represents the weight of the evaluation index corresponding to the first index layer of the multiple evaluation indexes of the second index layer. Wherein, the method of obtaining the second geothermal index weight vector corresponding to the second judgment matrix and the third geothermal index weight vector corresponding to the third judgment matrix is the same as that of obtaining the first geothermal index weight vector. As can be seen from Table 2, , , is the second judgment matrix, and the corresponding , , is the second geothermal index weight vector. And , is the third judgment matrix, and the corresponding , is the third geothermal index weight vector.

[0029] Table 2 The second judgment matrix corresponds to obtain a second geothermal index weight vector; the third judgment matrix corresponds to obtain a third geothermal index weight vector; the second geothermal index weight vector represents the weight of the first index layer corresponding to the criterion layer; and the third geothermal index weight vector represents the weight of the second index layer corresponding to the first index layer. Based on the first geothermal index weight vector, the second geothermal index weight vector and the third geothermal index weight vector, the accuracy of the detected weight between the evaluation indexes of the geothermal resource potential is determined, to obtain a first index weight accurate value, a second index weight accurate value and a third index weight accurate value.

[0030] The first index weight accurate value, the second index weight accurate value and the third index weight accurate value are determined to obtain an index weight set.

[0031] Optionally, the importance proportion of the adjacent level evaluation indexes in the geothermal evaluation index system is detected to obtain the first judgment matrix, the second judgment matrix and the third judgment matrix, including: The first initial judgment matrix is constructed by taking the multiple first-level evaluation indexes in the criterion layer corresponding to one target index in the target layer as rows and columns; all the values in the first initial judgment matrix are 0. In this embodiment, the heat source, the heat reservoir, the cap rock, the water source and the channel are taken as rows and columns in sequence.

[0032] The importance of the evaluation indexes in the target layer and the multiple evaluation indexes in the criterion layer is determined respectively, and the corresponding positions in the first initial judgment matrix are filled to obtain the first judgment matrix.

[0033] In this embodiment, the first judgment matrix is shown in Table 3. Table 3 In this embodiment, the row corresponding to A1 and the column corresponding to A2 are 3, which means that the importance of the evaluation index corresponding to A1 to the evaluation index corresponding to Z is 3 times the importance of the evaluation index corresponding to A2 to the evaluation index corresponding to Z. In this embodiment, the evaluation index corresponding to A1 is the heat source, the evaluation index corresponding to A2 is the channel, and the evaluation index corresponding to Z is the geothermal resource potential, so the importance of the heat source to the determination of the geothermal resource potential is 3 times the importance of the channel to the determination of the geothermal resource potential.

[0034] The criterion layer and the first index layer correspondingly obtain n second judgment matrices; the first index layer and the second index layer correspondingly obtain m third judgment matrices; n represents the number of evaluation indexes in the first index layer; m represents the number of evaluation indexes in the second index layer. Wherein, n and m are positive integers. In this embodiment, the number of evaluation indexes in the criterion layer is 5, the number of evaluation indexes in the first index layer is 10, and the number of evaluation indexes in the second index layer is 20, then n is 10 and m is 20. Wherein, Table 3-1 represents the second judgment matrix , Table 3-2 represents the second judgment matrix , Table 3-3 represents the second judgment matrix , Table 3-4 represents the third judgment matrix , Table 3-5 represents the third judgment matrix , Table 3-6 represents the third judgment matrix , Table 3-7 represents the third judgment matrix , Table 3-8 represents the third judgment matrix , Table 3-9 represents the third judgment matrix , Table 3-10 represents the third judgment matrix , Table 3-11 represents the third judgment matrix , Table 3-12 represents the third judgment matrix , Table 3-13 represents the third judgment matrix .

[0035] Table 3-1 Table 3-2 Table 3-3 Table 3-4 Table 3-5 Table 3-6 Table 3-7 Table 3-8 Table 3-9 Table 3-10 Table 3-11 Table 3-12 Table 3-13 Optionally, the geothermal resource potential evaluation map is obtained based on the set of geothermal resource data and the set of index weights, comprising: A plurality of geothermal resource scores are obtained based on the set of geothermal resource data; one geothermal resource score corresponds to one geothermal resource data.

[0036] The plurality of geothermal resource scores are multiplied by the corresponding index weights and then added to obtain a geothermal resource potential value; the geothermal resource potential value represents the development potential of a location; a plurality of locations correspond to a plurality of geothermal resource potential values.

[0037] The plurality of geothermal resource potential values are marked at a plurality of locations through gridding to obtain a geothermal resource potential evaluation map; the geothermal resource potential evaluation map represents the development potential of each location. In this embodiment, the geothermal resource potential evaluation map is obtained through the combination of Arcgis gridding, information mapping, and overlay analysis.

[0038] Optionally, the accuracy of the detected weights between the evaluation indexes of the geothermal resource potential is twice discriminated based on the first geothermal index weight vector, the second geothermal index weight vector, and the third geothermal index weight vector to obtain the first index weight accurate value, the second index weight accurate value, and the third index weight accurate value, comprising: The first judgment matrix is multiplied by the corresponding first geothermal index weight vector to obtain a first maximum eigenvalue vector; the value equal to the first maximum eigenvalue vector after multiplication with the first geothermal index weight vector is found as the first maximum eigenvalue. Wherein, A represents the first judgment matrix, w represents the corresponding first geothermal index weight vector, represents the first maximum eigenvalue. The accuracy of the detected weights between the evaluation indexes of the geothermal resource potential is twice discriminated based on the first maximum eigenvalue and the first judgment matrix to obtain the first index weight accurate value. The second judgment matrix corresponds to obtain the second index weight accurate value; the third judgment matrix corresponds to obtain the third index weight accurate value.

[0039] Wherein, the second judgment matrix corresponds to obtain the second index weight accurate value, and the method for obtaining the second index weight accurate value is as follows: the second judgment matrix is multiplied by the corresponding second geothermal index weight vector to obtain a second maximum eigenvalue vector; the value equal to the second maximum eigenvalue vector after multiplication with the second geothermal index weight vector is found as the second maximum eigenvalue; the accuracy of the detected weights between the evaluation indexes of the geothermal resource potential is twice discriminated based on the first maximum eigenvalue and the first judgment matrix to obtain the second index weight accurate value, comprising: wherein the second index weight accuracy value is obtained by multiplying the third judgment matrix and the third geothermal index weight vector, obtaining a third maximum eigenvalue vector; finding a value equal to the third maximum eigenvalue vector after multiplying the third geothermal index weight vector as a third maximum eigenvalue; based on the third maximum eigenvalue and the third judgment matrix, twice judging the accuracy of the weight detected between the evaluation indexes of the geothermal resource potential, obtaining a third index weight accuracy value, including: indicating the third maximum eigenvalue, indicating the number of values of the three index layer, indicating the average random consistency index of the third index layer. Wherein, the average random consistency index can be known by looking up the table, such as Z

[0040] The corresponding number of evaluation indexes is 5, such as The number of evaluation indexes is 4, such as The number of evaluation indexes is 3, . .

[0041] Optionally, the method further comprises: obtaining a geothermal standard table; the geothermal standard table includes a plurality of standard score values of positions; the standard score value indicates the score of the value corresponding to the plurality of evaluation indexes of the second index layer in different ranges. Wherein, the geothermal standard table is shown in Table 4 as follows: Table 4 finding the score corresponding to the geothermal resource data in the geothermal standard set to obtain a geothermal resource score set; the geothermal resource score set includes a plurality of geothermal resource scores; the geothermal resource score indicates the score of the geothermal resource corresponding to the stratified heat storage index or the strip-shaped heat storage evaluation index. Wherein, the geothermal resource score value indicates the score corresponding to the range where each value is located.

[0042] ​​​​​​Optionally, the first index weight accurate value, the second index weight accurate value and the third index weight accurate value are discriminated to obtain an index weight set, including: if the first index weight accurate value is less than a geothermal relationship judgment threshold value, and the second index weight accurate value is less than the geothermal relationship judgment threshold value, and the third index weight accurate value is less than the geothermal relationship judgment threshold value, it is determined that the geothermal relationship judgment is accurate. In the embodiment, the geothermal relationship judgment threshold value is 0.1. If the geothermal relationship judgment is accurate, the value in the third geothermal index weight vector is multiplied by the value in the corresponding second geothermal index weight vector, and the value in the corresponding second geothermal index weight vector is multiplied to obtain an index weight set; the index weight set includes a plurality of index weights; the number of elements in the index weight set is equal to the number of evaluation indexes in the second index layer. Because the geothermal heat flow value (C1) in the third geothermal index weight vector corresponds to the geothermal temperature field (B1) in the first index layer, and the geothermal temperature field (B1) corresponds to the heat source (A1) in the criterion layer, the weight of the geothermal heat flow value (C1) corresponding to the geothermal resource potential is . In the embodiment, the plurality of index weights are the weights of the evaluation indexes in the second index layer to the geothermal resource potential, as shown in Table 5.

[0043] Table 5 If the geothermal relationship judgment is not accurate, the first judgment matrix, the second judgment matrix and the third judgment matrix are adjusted. The importance of each evaluation index in the first judgment matrix, the second judgment matrix and the third judgment matrix is adjusted.

[0044] Optionally, the accuracy of the detected weight between the evaluation indexes of the geothermal resource potential is twice discriminated based on the first maximum eigenvalue and the first judgment matrix to obtain a first index weight accurate value, including: wherein the first index weight accurate value is the value of the first index weight accurate value, the first maximum eigenvalue is the value of the first maximum eigenvalue, the number of values in the criterion layer is the value of the number of values in the criterion layer, and the average random consistency index of the criterion layer is the value of the average random consistency index of the criterion layer. The table 6 indicates the accuracy of the detected weight between the evaluation indexes of the geothermal resource potential which is twice discriminated by the first judgment matrix, the second judgment matrix and the third judgment matrix.

[0045] Table 6 ​​​​​Optionally, the first geothermal index weight vector is obtained by judging the relationship of the evaluation indexes based on the first judgment matrix, comprising: geometrically averaging the row vectors of the first judgment matrix to obtain a plurality of first average values; and normalizing the plurality of first average values to obtain the first geothermal index weight vector.

[0046] Embodiment 2 Based on the above-mentioned geothermal resource evaluation method, the embodiment of the present application further provides a geothermal resource evaluation system, the system comprising an acquisition module, a hierarchy module, a correlation module, a weight module and a potential module.

[0047] The acquisition module is configured to acquire a plurality of geothermal resource data sets of a plurality of locations; the geothermal resource data set comprises a plurality of geothermal resource data; the hierarchy module is configured to construct evaluation indexes for the plurality of geothermal resource data; the geothermal evaluation index system is established by the stratiform heat reservoir and the zonal heat reservoir; the geothermal evaluation index system comprises a target layer, a criterion layer, a first index layer and a second index layer; the second index layer comprises stratiform heat reservoir evaluation indexes and zonal heat reservoir evaluation indexes; the correlation module is configured to detect the importance proportion of the evaluation indexes of adjacent levels in the geothermal evaluation index system to obtain a first judgment matrix, a second judgment matrix and a third judgment matrix; the weight module is configured to detect the weight of a plurality of evaluation indexes of the second index layer in the geothermal evaluation index system based on the first judgment matrix, the second judgment matrix and the third judgment matrix to obtain an index weight set; the index weight set comprises a plurality of index weights; one stratiform heat reservoir index or one zonal heat reservoir evaluation index corresponds to one index weight; and the potential module is configured to obtain a geothermal resource potential evaluation map based on the geothermal resource data set and the index weight set; the geothermal resource potential evaluation map is used to represent the geothermal resource development potential of different regions.

[0048] The specific implementation of the modules in the above-mentioned geothermal resource evaluation system can refer to the specific implementation of the corresponding steps in the above-mentioned geothermal resource evaluation method, which will not be described here again.

[0049] Embodiment 3 Taking the data of a certain region as an example, a plurality of geothermal resource data are determined by experts through multiple scoring based on various survey and exploration results, latest research results of scientific research units and the like.

[0050] 1, Heat source (A1): The geothermal field (B1) is further divided into three evaluation indexes of geothermal flux value (C1), geothermal gradient (C2), and surface thermal display anomaly (C3). The geophysical characteristics (B2) are further divided into two evaluation indexes of crustal thickness (C4) and Moho depth (C5). Active faults (B3): Sichuan Province is located at the junction of the Eurasian plate and the Indian Ocean plate, with strong tectonic activity. Compared with third-level and lower-level deep active faults, the contribution of heat source is dominant. Therefore, in this overall planning, the first and second-level faults in the whole province are considered, and 5km on both sides of the first-level fault and 4km on both sides of the second-level fault are taken as the evaluation range of the fault. The relevant information of the evaluation index is shown in Table 7. Seismic activity (B4) is further divided into two indexes of seismic activity period (C9) and magnitude (C10). The relevant information of the evaluation index is shown in Table 7.

[0051] Table 7 2, Channel (A2): Stratum structure (B5): The stratum lithology (C11) evaluation index distribution in the whole province is selected according to the geothermal geological map of the whole province. The fault density (C12) is taken as the evaluation block of 100km2 to obtain the fault density evaluation index distribution map of the whole province. The runoff distance (B6) is further divided into evaluation index and supply area plane distance (C13).

[0052] 3, Cap rock (A3): Further divided into evaluation index cap rock thickness (B7, C14).

[0053] 4, Heat reservoir (A4): Including heat reservoir structure (B8) including heat reservoir thickness (C15) stratified heat reservoir evaluation system, which is determined comprehensively according to the collected basin different stratum buried depth contour map, various drilling stratum exposure results, regional survey data, etc.; the evaluation index of zonal heat reservoir is not considered. The fault width (C16) considers the first, second and third-level faults in the whole province, and 5km on both sides of the first-level fault, 4km on both sides of the second-level fault, and 2km on both sides of the third-level fault as the evaluation of fault width. The heat reservoir lithology (C17) is confirmed according to the surface outcropped lithology within the evaluation range of the first, second and third-level faults, and the area with buried depth less than 4000m is determined according to the geothermal well and buried depth contour line. Heat reservoir temperature (B9) can be calculated by geothermal gradient and geochemical geothermal thermometer.

[0054] 5, Water source (A5): namely, groundwater recharge, further divided into two evaluation indexes of infiltration coefficient (C19) and multi-year average precipitation (C20). Infiltration coefficient (C19) is selected according to the results of 200,000, 500,000 and other hydrogeological survey work, various types of exploration and evaluation reports. Finally, the infiltration coefficient evaluation index distribution map of the whole province is obtained. Multi-year average precipitation (C20) is obtained according to the collected multi-year average precipitation data, and the multi-year average precipitation evaluation index distribution map of the whole province is obtained. The geothermal standard set is shown in Table 8. Among them, the scale of the third-order fault is less than the first and second order, the extension length is uniformly less than 100km, the fault activity is in the Pleistocene and before, the paleoseismic activity in the earthquake period, and the magnitude is less than 4. The score is evaluated.

[0055] Table 8 According to the score value of each evaluation index and the corresponding weight, the geothermal resource potential evaluation map of the whole province is obtained by using Arcgis rasterization and weighted superposition. As shown in Table 9, the area with score > 80 is defined as resource potential large area, the area with score 60-80 is defined as resource potential medium area, and the area with score < 60 is defined as resource potential small area.

[0056] Table 9 From the evaluation results: The basin and the surrounding geothermal area of the layered heat reservoir are good, especially near the Longmen Mountain, Micang Mountain, Dabashan Mountain, Wumeng Mountain and Huayun Mountain, the supply and discharge conditions are better than those in the basin, and the maximum score reaches 85; in general, the Triassic heat reservoir potential is slightly higher than that of the Permian heat reservoir, which is related to the shallow depth and large thickness of the Triassic heat reservoir.

[0057] The overall belt-shaped heat reservoir in the west, northwest and southwest of Sichuan is better, and the overall west geothermal area is better than the southwest geothermal area, which is better than the northwest geothermal area; it reflects the characteristics of the belt-shaped heat reservoir controlled by structure, and the potential large area is located in the Xianshuihe geothermal belt, Ganzi-Litang geothermal belt, Dege-Xiangcheng geothermal belt, Jinshajiang geothermal belt and Anninghe geothermal belt, which is consistent with the actual understanding; the maximum score in the area reaches 88, involving Kangding, Litang and Batang, which is consistent with the fact that the above-mentioned areas are rich in geothermal resources.

[0058] The distribution of existing mining rights and geothermal points is consistent with the distribution of potential medium to large areas.

[0059] The secondary fault has a significant control effect on the supply, migration and discharge of geothermal resources, and is very beneficial to the formation of geothermal resources. The evaluation result score is medium and above, and the evaluation result is more consistent with the actual situation.

[0060] According to the geothermal resource potential evaluation system and the typical geothermal point evaluation verification, the evaluation results are consistent with the understanding of the law of geothermal resources in the province, and basically coincide with the distribution and burial of geothermal resources, which can be used as a reference for later geothermal resource evaluation. Based on the above potential evaluation results, the Kangding Yulin River-Yala River, Tianquan-Xingjing in Ya'an, Litang Basin area, Guangyuan and other areas can be comprehensively considered for geothermal resource exploration in the future.

[0061] The algorithms and displays presented herein are not inherently related to any particular computer, virtual system, or other apparatus. Various general purpose systems can be used with these teachings, based upon the description as provided herein. Those skilled in the art will recognize that structures described above as being needed for the implementation can be structured in various ways. Furthermore, the present application is not intended to be limited to any particular programming language. It will be appreciated that there are many programming languages that can be used to implement the teachings described herein, and any such programming language can be used.

Claims

1. A geothermal resource evaluation method, characterized in that: include: Acquire a geothermal resource data set of multiple locations; the geothermal resource data set includes multiple geothermal resource data; Construct evaluation indicators for multiple geothermal resource data; A geothermal evaluation index system is established based on both layered heat reservoirs and strip heat reservoirs; the geothermal evaluation index system includes a target layer, a criterion layer, a first index layer, and a second index layer; the second index layer includes layered heat reservoir evaluation indicators and strip heat reservoir evaluation indicators; Detect the importance ratios of adjacent hierarchical evaluation indicators in the geothermal evaluation indicator system to obtain a first judgment matrix, a second judgment matrix, and a third judgment matrix; Based on the first judgment matrix, the second judgment matrix and the third judgment matrix, the weights of multiple evaluation indicators of the second indicator layer in the geothermal evaluation indicator system are detected to obtain an indicator weight set; The indicator weight set includes multiple indicator weights; one layered heat storage indicator or strip heat storage evaluation indicator corresponds to one indicator weight; Based on the geothermal resource data set and the indicator weight set, a geothermal resource potential evaluation map is obtained; the geothermal resource potential evaluation map is used to characterize the geothermal resource development potential of different regions.

2. The geothermal resource evaluation method according to claim 1, characterized in that: The method of detecting the weights of multiple evaluation indicators of the second indicator layer in the geothermal evaluation indicator system based on the first judgment matrix, the second judgment matrix, and the third judgment matrix to obtain an indicator weight set includes: Based on the first judgment matrix, the relationship between the evaluation indicators is judged to obtain a first geothermal indicator weight vector; the value in the first geothermal indicator weight vector represents the weight of the multiple evaluation indicators of the criterion layer corresponding to the evaluation indicators of the target layer; The second judgment matrix corresponds to obtaining a second geothermal index weight vector; the third judgment matrix corresponds to obtaining a third geothermal index weight vector; the second geothermal index weight vector represents the weights of the multiple evaluation indicators of the first indicator layer corresponding to the evaluation indicators of the criterion layer; the third geothermal index weight vector represents the weights of the multiple evaluation indicators of the second indicator layer corresponding to the evaluation indicators of the first indicator layer; Based on the first geothermal indicator weight vector, the second geothermal indicator weight vector and the third geothermal indicator weight vector, the accuracy of the weights detected between the evaluation indicators of the geothermal resource potential is secondary judged to obtain the first indicator weight accurate value, the second indicator weight accurate value and the third indicator weight accurate value; The accurate value of the first indicator weight, the accurate value of the second indicator weight and the accurate value of the third indicator weight are determined to obtain an indicator weight set.

3. The geothermal resource evaluation method according to claim 1, characterized in that: The detection of the importance ratios of the evaluation indicators of adjacent levels in the geothermal evaluation index system to obtain a first judgment matrix, a second judgment matrix, and a third judgment matrix includes: A first initial judgment matrix is ​​constructed by using multiple evaluation indicators in the criterion layer corresponding to one target indicator in the target layer as rows and columns; all values ​​in the first initial judgment matrix are 0; The importance of the evaluation indicators in the target layer and the multiple evaluation indicators in the criterion layer are respectively determined, and the corresponding positions in the first initial judgment matrix are filled in to obtain the first judgment matrix; The criterion layer and the first indicator layer obtain n second judgment matrices correspondingly; the first indicator layer and the second indicator layer obtain m third judgment matrices correspondingly; n represents the number of evaluation indicators in the first indicator layer; m represents the number of evaluation indicators in the second indicator layer.

4. The geothermal resource evaluation method according to claim 1, characterized in that: The step of obtaining a geothermal resource potential evaluation map based on the geothermal resource data set and the indicator weight set includes: Based on the geothermal resource data set, multiple geothermal resource scores are obtained; one geothermal resource score corresponds to one geothermal resource data; The geothermal resource potential value is obtained by multiplying the multiple geothermal resource scores by the corresponding indicator weights and then adding them together; the geothermal resource potential value represents the development potential of a location; multiple geothermal resource potential values ​​are obtained for multiple locations; By rasterizing, multiple geothermal resource potential values ​​are marked at multiple locations to obtain a geothermal resource potential evaluation map; the geothermal resource potential evaluation map represents the development potential of each location.

5. The geothermal resource evaluation method according to claim 2, characterized in that: The method of secondarily determining the accuracy of the weights detected between the evaluation indicators of geothermal resource potential based on the first geothermal indicator weight vector, the second geothermal indicator weight vector, and the third geothermal indicator weight vector to obtain the first indicator weight accurate value, the second indicator weight accurate value, and the third indicator weight accurate value includes: Multiplying the first judgment matrix by the corresponding first geothermal index weight vector to obtain a first maximum eigenvalue vector; Find the value that is equal to the first maximum eigenvalue after multiplying the first geothermal index weight vector as the first maximum eigenvalue; Based on the first maximum eigenvalue and the first judgment matrix, the accuracy of the weights detected between the evaluation indicators of geothermal resource potential is secondarily judged to obtain the accurate value of the first indicator weight; The second judgment matrix corresponds to obtaining the accurate value of the second indicator weight; the third judgment matrix corresponds to obtaining the accurate value of the third indicator weight.

6. The geothermal resource evaluation method according to claim 4, characterized in that: The method of obtaining a plurality of geothermal resource scores based on the geothermal resource data set includes: Obtaining a geothermal standard table; the geothermal standard table includes standard scoring values ​​for multiple locations; the standard scoring values ​​represent scores corresponding to values ​​of multiple evaluation indicators of the second indicator layer in different ranges; The scores corresponding to the geothermal resource data are found in the geothermal standard set to obtain a geothermal resource score set; the geothermal resource score set includes multiple geothermal resource scores; the geothermal resource score represents the score of the geothermal resource corresponding to the layered heat storage index or the strip heat storage evaluation index.

7. The geothermal resource evaluation method according to claim 2, characterized in that: The determining of the first indicator weight accurate value, the second indicator weight accurate value, and the third indicator weight accurate value to obtain an indicator weight set includes: If the first indicator weight accuracy value is less than the geothermal relationship judgment threshold, and the second indicator weight accuracy value is less than the geothermal relationship judgment threshold, and the third indicator weight accuracy value is less than the geothermal relationship judgment threshold, it is assumed that the geothermal relationship judgment is accurate; If the geothermal relationship is accurately determined, the value in the third geothermal index weight vector is multiplied by the value in the corresponding second geothermal index weight vector to obtain an index weight set; the index weight set includes multiple index weights; the number of elements in the index weight set is equal to the number of evaluation indicators of the second index layer; If the geothermal relationship judgment is inaccurate, the first judgment matrix, the second judgment matrix, and the third judgment matrix are adjusted.

8. The geothermal resource evaluation method according to claim 5, characterized in that: The method of secondarily judging the accuracy of the weights detected between the evaluation indicators of geothermal resource potential based on the first maximum eigenvalue and the first judgment matrix to obtain the accurate value of the first indicator weight includes: Among them, the is the accurate value of the first indicator weight, represents the first largest eigenvalue, represents the number of values ​​of the criterion layer, represents the average random consistency index of the criterion layer.

9. The geothermal resource evaluation method according to claim 2, characterized in that: The method of determining the relationship between the evaluation indicators based on the first judgment matrix to obtain the first geothermal indicator weight vector includes: Performing geometric averaging on the row vectors of the first judgment matrix to obtain multiple first average values; The multiple first average values ​​are normalized to obtain a first geothermal index weight vector.

10. A geothermal resource evaluation system, characterized in that: include: An acquisition module is used to acquire a geothermal resource data set of multiple locations; the geothermal resource data set includes multiple geothermal resource data; A hierarchical module is used to construct evaluation indicators for multiple geothermal resource data; a geothermal evaluation indicator system is established based on both layered heat reservoirs and strip heat reservoirs; the geothermal evaluation indicator system includes a target layer, a criterion layer, a first indicator layer, and a second indicator layer; the second indicator layer includes layered heat reservoir evaluation indicators and strip heat reservoir evaluation indicators; An association module is used to detect the importance ratios of evaluation indicators at adjacent levels in the geothermal evaluation indicator system to obtain a first judgment matrix, a second judgment matrix, and a third judgment matrix; A weight module is used to detect the weights of multiple evaluation indicators of the second indicator layer in the geothermal evaluation indicator system based on the first judgment matrix, the second judgment matrix, and the third judgment matrix to obtain an indicator weight set; the indicator weight set includes multiple indicator weights; one layered heat storage indicator or strip heat storage evaluation indicator corresponds to one indicator weight; The potential module is used to obtain a geothermal resource potential evaluation map based on the geothermal resource data set and the indicator weight set; the geothermal resource potential evaluation map is used to characterize the geothermal resource development potential size in different regions.

Images