A method for dry hot rock identification of gravity, magnetic and electric integrated parameter exploration

Abstract

The application discloses a gravity-magnetic-electric comprehensive parameter exploration method for hot dry rock identification, and comprises the following steps: S1, acquiring gravity-magnetic-electric inversion parameters corresponding to a survey area in space, and performing grid processing on the gravity-magnetic-electric inversion parameters; S2, performing standardization processing on the grid inversion parameters; S3, calculating correlation coefficients of the density parameter and the magnetic parameter respectively by taking the resistivity as a main parameter; S4, calculating gravity-magnetic-electric comprehensive parameters; and S5, analyzing and interpreting the gravity-magnetic-electric comprehensive parameters. Through analyzing the gravity-magnetic-electric comprehensive parameters, the hot dry rock identification and delineation can be provided with higher exploration accuracy, so that the problems of low hot dry rock identification accuracy and strong randomness of identification are solved.

Description

Technical Field

[0001] This invention relates to the field of geophysical exploration, and in particular to a method for exploring comprehensive gravity, magnetic and electrical parameters for identifying hot dry rocks. Background Technology

[0002] Hot dry rocks are commonly found in metamorphic or crystalline rock formations. These rocks possess very high temperatures over a wide range and are in a hot-dry state. Due to current advanced drilling techniques and the fact that thermal energy systems do not require the rocks to be porous or contain fluids, hot dry rocks can be utilized as a thermal energy resource. Conservative estimates suggest that the energy contained in hot rocks (3–10 km deep) within the Earth's crust is equivalent to 30 times the energy contained in all of the world's oil, natural gas, and coal reserves.

[0003] Therefore, strengthening the exploration of hot dry rocks is particularly important for meeting my country's energy needs; how to find hot dry rocks more effectively, quickly and accurately is a topic worthy of study.

[0004] The main methods for exploring hot dry rocks currently include:

[0005] Geothermal fluidity method. This is the most direct method, and the study area is usually selected in sedimentary rock areas with normal thermal gradients, large depths, and relatively insulating conditions, or areas with high heat production (such as granite areas overlying sedimentary caps). The key points of this method are identifying high heat flow areas, determining basement temperature, measuring basement temperature, measuring rock cores for thermal conductivity and heat production, and studying the temperature at the top of the basement, basement heat production, heat flow on the basement surface, and surface heat flow. It primarily involves analysis from a planar perspective.

[0006] Seismic methods are used to determine crustal structure and geothermal conclusions that are thermally related but not directly related. This method is applicable to locating disturbed areas exhibiting characteristics of hydrothermal or magmatic heating.

[0007] Magnetotelluric methods play an important role in the exploration of hot and dry rocks because the conductivity of rocks is greatly affected by their water content and temperature.

[0008] Gravity and magnetic methods. Because they allow for large-area gravity and magnetic measurements, and the data shows a clear correlation with temperature, intrusions, and tectonic frameworks, they offer unique advantages for delineating large-scale hot dry rock target areas.

[0009] Existing geophysical exploration methods, especially gravity, magnetic, and electrical (GMO) methods for exploring hot dry rocks, either rely on single-parameter judgment and manually combine the characteristics of different parameters to identify target areas, or employ a combined GMO inversion approach to attempt to establish a relationship between hot dry rocks and their geophysical exploration parameters. The first approach is often highly subjective, with different people easily obtaining different results, making standardized operations difficult and reducing the accuracy of hot dry rock exploration. The second approach requires more manpower, financial resources, and time, although it can achieve relatively higher accuracy results. However, due to the ambiguity of geophysical data processing and inversion, the various GMO parameters are not strictly correlated, and the relationships between the various GMO parameters of the target body are not simple linear. Therefore, combined inversion of multiple GMO parameters can sometimes lead to a decrease in exploration accuracy. Summary of the Invention

[0010] In view of the above problems, the present invention is proposed to provide a comprehensive gravity, magnetic and electrical parameter exploration method for identifying hot dry rocks that overcomes or at least partially solves the above problems.

[0011] According to one aspect of the present invention, a method for exploring comprehensive gravity, magnetic, and electrical parameters for identifying hot dry rocks is provided, the method comprising:

[0012] Step S1: Obtain the spatially corresponding gravity, magnetoelectric inversion parameters of the survey area, and perform gridding processing on the gravity, magnetoelectric inversion parameters;

[0013] Step S2: Standardize the gridded inversion parameters;

[0014] Step S3: Calculate the correlation coefficients between resistivity and density parameters and magnetic parameters, respectively, using resistivity as the main parameter;

[0015] Step S4: Calculate the combined gravity, magnetism, and electrical parameters;

[0016] Step S5: Analyze and interpret the combined gravity, magnetism, and electrical parameters.

[0017] Optionally, step S1: obtaining the spatially corresponding gravity, magnetoelectric inversion parameters of the survey area, and performing gridding processing on the gravity, magnetoelectric inversion parameters specifically includes:

[0018] In the survey area where comprehensive gravity, magnetic and electrical exploration of hot dry rocks is carried out, the grid parameters for gravity, magnetic and electrical exploration are designed according to the requirements of comprehensive exploration work, comprehensive gravity, magnetic and electrical exploration work is carried out, the gravity, magnetic and electrical parameters of multiple measuring points are obtained, and the relevant parameters are inverted and calculated.

[0019] Obtain gravity, magnetic and electrical inversion parameters for each measuring point at different elevations, and perform gridded calculations on the gravity, magnetic and electrical inversion parameters;

[0020] Obtain the gridded inversion parameters W of gravity, magnetism, and electricity at different elevations for each gridded measuring point. k (x i ,y i ,z i );

[0021] When the planar measuring points of gravity, magnetic and electrical exploration do not completely overlap, the measuring points that overlap using three different methods are selected for gridding.

[0022] Where x i y i These are the planar coordinates of the measurement points after meshing; z i W represents the elevation corresponding to the gridded measurement point; where the subscript i is the number of the gridded measurement point, and i∈[1,n]; n is the total number of gridded measurement points, where n is a natural number not less than 1; the superscript k is the number of the gridded inversion parameter for gravity, magnetic and electrical exploration, and k∈[1,3]; where W 1 (x i ,y i ,z i ) represents the density gridded inversion parameters obtained from gravity exploration, W 2 (x i ,y i ,z i ) represents the magnetic gridding inversion parameters obtained from magnetic exploration, W 3 (x i ,y i ,z i ) represents the resistivity gridded inversion parameters obtained from electrical exploration.

[0023] Optionally, the inversion calculation of relevant parameters specifically includes: inverting the gravity, magnetic and electrical parameters obtained from gravity, magnetic and electrical exploration to obtain inversion results of relevant density, magnetic parameters and resistivity parameters.

[0024] Optionally, obtaining the gravity, magnetoelectric inversion parameters for each measuring point at different elevations specifically includes obtaining the relevant density, magnetic parameters, and resistivity parameters for different measuring points at different depths.

[0025] Optionally, when the planar measuring points of gravity, magnetic, and electrical exploration do not completely overlap, the measuring points overlapping by three methods are selected for gridding processing, specifically including:

[0026] When the minimum value of a certain gridded inversion parameter of gravity, magnetoelectricity is negative, then add twice the absolute value of the minimum value of the gridded parameter to all such gridded parameters.

[0027] When the minimum value of a certain gridded inversion parameter of gravity, magnetoelectricity is zero, the absolute value of the maximum value of the inversion parameter is added to all such inversion parameters.

[0028] If a certain gridded inversion parameter of gravity, magnetoelectricity is 0, then no further calculation of that parameter will be performed.

[0029] Optionally, step S2: standardizing the gridded inversion parameters specifically includes:

[0030] According to formula (1), the obtained elevation z i The gridded inversion parameters of all measuring points are standardized to unify all data into the interval [0,1], thus obtaining the elevation z of all measuring points. i Standardized parameters of gridded inversion parameters

[0031] Where p is the elevation z i The total number of a certain gridded inversion parameter; k is the gridded inversion parameter number of gravity, magnetic and electrical exploration, and k∈[1,3]; that is, W 1 (x i ,y i ,z i ) represents the density gridded inversion parameters obtained from gravity exploration, W 2 (x i ,y i ,z i ) represents the magnetic gridding inversion parameters obtained from magnetic exploration, W 3 (x i ,y i ,z i The resistivity gridded inversion parameters are obtained from electrical exploration.

[0032] in

[0033]

[0034] Repeat the previous step to calculate the standardized parameters of the gravity, magnetoelectric and electrical gridded inversion parameters for the remaining elevations, until the standardized parameters of the gravity, magnetoelectric and electrical gridded inversion parameters for all elevations are obtained.

[0035] Optionally, step S3: calculating the correlation coefficients with the density parameter and magnetic parameter, respectively, using resistivity as the main parameter, specifically includes:

[0036] Formula (2) is used to analyze a certain elevation z i The correlation between the normalized parameters of the resistivity gridded inversion parameters and the normalized parameters of the density gridded inversion parameters or the normalized parameters of the magnetic gridded inversion parameters at the same elevation;

[0037] Determine the elevation z i The correlation coefficients between the normalized parameters of the resistivity gridded inversion parameters and the normalized parameters of the density gridded inversion parameters or the normalized parameters of the magnetic gridded inversion parameters at the same elevation.

[0038] Where N is the total number of measurement points participating in the correlation analysis of standardized parameters; the subscript i is the number of the measurement point after gridding, and i∈[1,N]; The three-dimensional coordinates of the measuring point are (x i ,y i ,z i The standardized parameters of resistivity gridded inversion parameters involved in correlation analysis; The three-dimensional coordinates of the measuring point are (x i ,y i ,z i The parameter is a standardized parameter in the density and magnetic gridding inversion parameters involved in the correlation analysis; the subscript l is the standardized parameter number of the density and magnetic gridding inversion parameters of gravity, magnetic and electrical exploration, and l∈[1,2]. These are the standardized parameters for density gridded inversion parameters obtained from gravity exploration. These are the standardized parameters for magnetic gridding inversion parameters obtained from magnetic exploration.

[0039] Where σ1≠0, σ2≠0, σ3≠0;

[0040]

[0041] in:

[0042]

[0043] Repeat the previous step to obtain the correlation coefficients between the standardized parameters of the resistivity gridded inversion parameters at the remaining elevation and the standardized parameters of the density gridded inversion parameters and magnetic gridded inversion parameters at the same elevation. The correlation coefficients between the normalized parameters of the resistivity gridded inversion parameters at all elevations and the normalized parameters of the density gridded inversion parameters and the normalized parameters of the magnetic gridded inversion parameters at the same elevations were obtained.

[0044] Optionally, step S4: calculating the combined gravity, magnetism, and electrical parameters specifically includes:

[0045] Select a certain elevation z i The standardized parameters of the resistivity gridded inversion parameters are the primary standardized parameters;

[0046] Formula (3) is used to calculate all measuring points and a certain elevation z. i The comprehensive parameters of the standardized parameters

[0047] Repeat the previous step to obtain the comprehensive parameters of standardized parameters for all elevations and all measuring points.

[0048] Optionally, the selection of a certain elevation z i The standardized parameters of the resistivity gridded inversion parameters are the main standardized parameters, which specifically include:

[0049] Main standardized parameters; The standardized parameters are the density grid inversion parameters for the corresponding elevation. These are the standardized parameters for the magnetic grid inversion parameters at the corresponding elevation.

[0050] Optionally, the calculation of all measuring points and a certain elevation z using formula (3) is performed. i The comprehensive parameters of the standardized parameters Specifically, it includes:

[0051] in, For elevation z i The normalized parameters of the resistivity gridded inversion parameters and the elevation z i The correlation coefficients between the standardized parameters of the density gridded inversion parameters. For elevation z i The normalized parameters of the resistivity gridded inversion parameters and the elevation z i The correlation coefficients between the standardized parameters of the magnetic gridded inversion parameters;

[0052]

[0053] Optionally, step S5: analyzing and interpreting the combined gravity, magnetism, and electrical parameters specifically includes:

[0054] Comprehensive parameters for all measuring points and all elevations Analysis and interpretation were conducted to obtain the interpretation results of gravity, magnetic and electrical exploration of dry hot rocks in the survey area;

[0055] If a single profile exploration of gravity, magnetic, and electrical parameters is carried out, then two-dimensional contour lines are drawn for the comprehensive parameters of gravity, magnetic, and electrical parameters, and profile analysis and interpretation are performed.

[0056] If a multi-profile exploration of gravity, magnetoelectricity is carried out, then three-dimensional contour lines are drawn for the comprehensive parameters of gravity, magnetoelectricity, and electricity, and three-dimensional analysis and interpretation are performed.

[0057] This invention provides a method for exploring comprehensive gravity, magnetic, and electrical parameters for identifying hot dry rocks. The method includes: Step S1: Obtaining the spatially corresponding gravity, magnetic, and electrical inversion parameters of the survey area and performing gridding processing on these parameters; Step S2: Standardizing the gridded inversion parameters; Step S3: Calculating the correlation coefficients between resistivity and density and magnetic parameters, with resistivity as the main parameter; Step S4: Calculating the comprehensive gravity, magnetic, and electrical parameters; Step S5: Analyzing and interpreting the comprehensive gravity, magnetic, and electrical parameters. By analyzing the comprehensive gravity, magnetic, and electrical parameters, this method provides more accurate exploration results for the identification and delineation of hot dry rocks, thereby solving the problems of low accuracy and arbitrary identification of hot dry rocks.

[0058] The above description is merely an overview of the technical solution of the present invention. In order to better understand the technical means of the present invention and to implement it in accordance with the contents of the specification, and in order to make the above and other objects, features and advantages of the present invention more apparent and understandable, specific embodiments of the present invention are described below. Attached Figure Description

[0059] To more clearly illustrate the technical solutions of the embodiments of the present invention, the 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.

[0060] Figure 1 A flowchart of a comprehensive gravity, magnetic, and electrical parameter exploration method for identifying hot dry rocks provided in an embodiment of the present invention;

[0061] Figure 2 This is a schematic diagram of the field work layout provided in an embodiment of the present invention;

[0062] Figure 3 This is a density three-dimensional data volume obtained by gravity inversion provided in an embodiment of the present invention;

[0063] Figure 4 This is a three-dimensional data volume of magnetic susceptibility obtained by magnetic inversion according to an embodiment of the present invention;

[0064] Figure 5 This is a three-dimensional volume of resistivity logarithmic data obtained by magnetotelluric inversion provided in an embodiment of the present invention;

[0065] Figure 6 This is a three-dimensional data volume of density and resistivity correlation coefficients provided in an embodiment of the present invention;

[0066] Figure 7 This is a three-dimensional data volume of the correlation coefficient between magnetic susceptibility and resistivity provided in an embodiment of the present invention;

[0067] Figure 8 This is a three-dimensional data volume of gravity, magnetoelectric integrated parameters W3 (high density, high magnetic field, high resistance) provided in the embodiments of the present invention;

[0068] Figure 9 The three-dimensional data volume of gravity, magnetoelectric integrated parameters W3 (low density, high magnetic field, high resistance) provided in the embodiments of the present invention. Detailed Implementation

[0069] Exemplary embodiments of the present disclosure will now be described in more detail with reference to the accompanying drawings. While exemplary embodiments of the present disclosure are shown in the drawings, it should be understood that the present disclosure may be implemented in various forms and should not be limited to the embodiments set forth herein. Rather, these embodiments are provided so that this disclosure will be thorough and complete, and will fully convey the scope of the disclosure to those skilled in the art.

[0070] The terms "comprising" and "having," and any variations thereof, in the specification, embodiments, claims, and drawings of this invention are intended to cover non-exclusive inclusion, such as including a series of steps or units.

[0071] The technical solution of the present invention will be further described in detail below with reference to the accompanying drawings and embodiments.

[0072] Example 1

[0073] This invention standardizes the gravity, magnetic, and electrical parameters in geophysical exploration of hot dry rocks using relevant algorithms. Based on the correlation coefficients of these parameters, it obtains comprehensive gravity, magnetic, and electrical parameters more closely related to hot dry rocks. By analyzing these comprehensive parameters, it provides more accurate exploration results for the identification and delineation of hot dry rocks, solving problems such as low identification accuracy and high arbitrariness in hot dry rock identification. The steps include:

[0074] Step 1: By gridding the inversion parameters, obtain the spatially corresponding gravity, magnetoelectric and inversion parameters of the survey area;

[0075] In the survey area where comprehensive gravity, magnetic, and electrical (GMO) exploration of hot dry rocks is carried out, the GMO exploration grid parameters are designed according to the requirements of the comprehensive exploration work. Comprehensive GMO exploration is conducted to obtain GMO parameters from multiple measuring points, and inversion calculations are performed on relevant parameters (e.g., inversion of gravity, magnetic, and electrical parameters obtained from GMO exploration to obtain inversion results for parameters such as density, magnetic parameters, and resistivity). This yields the GMO inversion parameters for each measuring point at different elevations (i.e., obtaining relevant density, magnetic parameters, and resistivity parameters at different measuring points and depths). Because different inversion methods use different vertical grid partitioning methods, to obtain spatially one-to-one data, the GMO inversion parameters are calculated using a grid (to improve the comprehensive calculation between the three GMO parameters, the grid interval and grid range are taken to ensure comprehensive parameter calculation between relevant parameters). This yields the gridded GMO inversion parameters W for each gridded measuring point at different elevations. k (x i ,y i ,z i When the planar measuring points of gravity, magnetic, and electrical exploration do not completely overlap, the measuring points overlapping by the three methods are selected for gridding. To facilitate the subsequent calculation of comprehensive parameters, the gridded inversion parameters of gravity, magnetic, and electrical exploration must be positive, ensuring that all gridded inversion parameters of gravity, magnetic, and electrical exploration are positive and cannot have negative or zero values. When the minimum value of a certain gridded inversion parameter of gravity, magnetic, and electrical exploration is negative, then twice the absolute value of the minimum value of that gridded parameter is added to all such gridded parameters. When the minimum value of a certain gridded inversion parameter of gravity, magnetic, and electrical exploration is zero, then the absolute value of the maximum value of that inversion parameter is added to all such inversion parameters. When all gridded inversion parameters of gravity, magnetic, and electrical exploration are 0, then no further calculation of that parameter is performed. Where x i y i These are the planar coordinates of the measurement points after meshing; z i The elevation corresponding to the gridded measuring point (this parameter can also be defined as depth, and there is a one-to-one correspondence between depth and elevation, which can be converted to each other); where the subscript i is the number of the gridded measuring point, and i∈[1,n]; n is the total number of gridded measuring points, where n is a natural number not less than 1; the superscript k is the gridded inversion parameter number of gravity, magnetic and electrical exploration, and k∈[1,3]; where W 1 (x i ,y i ,z i ) represents the density gridded inversion parameters obtained from gravity exploration, W 2 (x i ,y i ,z i ) represents the magnetic gridding inversion parameters obtained from magnetic exploration, W 3 (x i ,yi ,z i The resistivity gridded inversion parameters are obtained from electrical exploration.

[0076] Step 2: Standardize the gridded inversion parameters;

[0077] According to formula (1), the elevation z obtained in step a) is... i (That is, only gridded inversion parameters with the same elevation are selected for subsequent standardized parameter calculations) A certain gridded inversion parameter of all measuring points is standardized to unify all data into the interval [0,1], eliminating the dimensional differences between different data, and obtaining the elevation z of all measuring points. i The standardized parameters of the gridded inversion parameters Where p is the elevation z i The total number of a certain gridded inversion parameter; k is the gridded inversion parameter number of gravity, magnetic and electrical exploration, and k∈[1,3]; that is, W 1 (x i ,y i ,z i ) represents the density gridded inversion parameters obtained from gravity exploration, W 2 (x i ,y i ,z i ) represents the magnetic gridding inversion parameters obtained from magnetic exploration, W 3 (x i ,y i ,z i ) represents the resistivity gridded inversion parameters obtained from electrical exploration; where

[0078]

[0079] Repeat the previous step to calculate the standardized parameters of the gravity, magnetoelectric, and electrical (MEO) gridded inversion parameters for other elevations until the standardized parameters of the MEO, MEO, and MEO gridded inversion parameters for all elevations are obtained. That is, MEO, MEO, and MEO gridded inversion parameters at different elevations do not calculate standardized parameters for each other; only MEO, MEO, and MEO gridded inversion parameters at the same elevation are standardized (i.e., if there are 101 different elevation values, then 101 sets of standardized parameters are calculated). For example, if we want to calculate the standardized parameters of the MEO, MEO, and MEO inversion parameters between 10000 and 0 meters (100-meter intervals), then we need to calculate the standardized parameters of the MEO, MEO, and MEO gridded inversion parameters at elevations of 10000 / 9900 / 9800 / 9700 / … / 300 / 200 / 100 / 0 meters, thus obtaining 101 sets of standardized parameters.

[0080] Step 3: Calculate the correlation coefficients between resistivity and density parameters and magnetic parameters, respectively, using resistivity as the main parameter;

[0081] Formula (2) is used to analyze a certain elevation z i The correlation between the normalized parameters of the resistivity gridded inversion parameters and the normalized parameters of the density gridded inversion parameters or the magnetic gridded inversion parameters at the same elevation is used to determine the elevation z. i The correlation coefficients between the normalized parameters of the resistivity gridded inversion parameters and the normalized parameters of the density gridded inversion parameters or the normalized parameters of the magnetic gridded inversion parameters at the same elevation. N represents the total number of measurement points participating in the correlation analysis of standardized parameters; the subscript i represents the number of the measurement point after gridding, and i∈[1,N]; The three-dimensional coordinates of the measuring point are (x i ,y i ,z i The standardized parameters of resistivity gridded inversion parameters involved in correlation analysis; The three-dimensional coordinates of the measuring point are (x i ,y i ,z i The standardized parameter in the density and magnetic gridded inversion parameters involved in the correlation analysis is a certain parameter; the subscript l is the standardized parameter number of the density and magnetic gridded inversion parameters of gravity, magnetic and electric exploration, and l∈[1,2]; These are the standardized parameters for density gridded inversion parameters obtained from gravity exploration. These are the standardized parameters for the magnetic gridding inversion parameters obtained from magnetic exploration; where σ1≠0, σ2≠0, σ3≠0;

[0082]

[0083] in:

[0084]

[0085] Repeat the previous step to obtain the correlation coefficients between the standardized parameters of the resistivity gridded inversion parameters at other elevations and the standardized parameters of the density gridded inversion parameters and magnetic gridded inversion parameters at the same elevation. The correlation coefficients between the normalized parameters of the resistivity gridded inversion parameters at all elevations and the normalized parameters of the density gridded inversion parameters and the normalized parameters of the magnetic gridded inversion parameters at the same elevations were obtained.

[0086] Step 4: Calculate the combined gravity, magnetic, and electrical parameters;

[0087] Select a certain elevation z i The standardized parameters of the resistivity gridded inversion parameters are the primary standardized parameters, i.e. Main standardized parameters; The standardized parameters are the density grid inversion parameters for the corresponding elevation. The standardized parameters for the magnetic grid inversion parameters at the corresponding elevation are used; formula (3) is used to calculate the parameters for all measuring points and a certain elevation z. i The comprehensive parameters of the standardized parameters in For elevation z i The normalized parameters of the resistivity gridded inversion parameters and the elevation z i The correlation coefficients between the standardized parameters of the density gridded inversion parameters. For elevation z i The normalized parameters of the resistivity gridded inversion parameters and the elevation z i The correlation coefficients between the standardized parameters of the magnetic gridded inversion parameters;

[0088]

[0089] Repeat the previous step to obtain the comprehensive parameters of standardized parameters for all elevations and all measuring points.

[0090] Step 5: Analyze and interpret the combined gravity, magnetism, and electrical parameters, including:

[0091] Comprehensive parameters for all measuring points and all elevations Analysis and interpretation are conducted to obtain interpretation results of gravity, magnetic and electrical exploration of dry hot rocks in the survey area; if a single-section gravity, magnetic and electrical exploration is carried out, two-dimensional contour lines are drawn for the comprehensive parameters of gravity, magnetic and electrical, and profile analysis and interpretation are carried out; if a multi-section gravity, magnetic and electrical exploration is carried out, three-dimensional contour lines are drawn for the comprehensive parameters of gravity, magnetic and electrical, and three-dimensional analysis and interpretation are carried out.

[0092] This invention is applied to the rapid delineation of regions with combined parameters of different parameters in the comprehensive exploration of hot, dry rocks using gravity, magnetic, and electrical systems, such as high-density, high-magnetic-resistivity regions and low-density, high-magnetic-resistivity regions. A novel algorithm is employed, which eliminates dimensional differences between different parameters through standardization, calculates correlation coefficients between different parameters, and obtains comprehensive parameters for analysis and interpretation, thereby improving the accuracy and effectiveness of target identification.

[0093] Example 2

[0094] like Figure 1 As shown, a comprehensive gravity, magnetic, and electrical parameter exploration method for identifying hot dry rocks includes:

[0095] a) In the survey area where comprehensive gravity, magnetic, and electrical exploration of hot dry rocks is carried out (e.g. Figure 2 As shown), the gravity, magnetic, and electrical exploration grid parameters are designed according to the requirements of comprehensive exploration work (e.g., Figure 2As shown, survey lines and points were arranged according to a 200m*200m grid to carry out comprehensive gravity, magnetic, and electrical exploration, obtaining multiple survey points (such as...). Figure 2 As shown, the gravity, magnetic, and electrical parameters of a total of 51*51=2601 measuring points were obtained, and inversion calculations were performed on the relevant parameters (such as inverting the gravity, magnetic, and electrical parameters obtained from gravity, magnetic, and electrical exploration to obtain inversion results for parameters such as density, magnetic parameters, and resistivity). Figure 3 , Figure 4 , Figure 5 As shown, gravity, magnetic, and electrical exploration data from 2601 measuring points were inverted to obtain gravity, magnetic, and electrical inversion results for elevations from 0 to 10000 meters. Gravity, magnetic, and electrical inversion parameters were obtained for each measuring point at different elevations (i.e., relevant density, magnetic parameters, resistivity, etc., for different measuring points and depths). The gravity, magnetic, and electrical inversion parameters were then calculated using a grid (to improve the comprehensive calculation of the three parameters, the grid interval and grid range were kept consistent during the gridding calculation to facilitate the comprehensive calculation of related parameters; assuming...). Figure 2 The gravity, magnetic, and electrical exploration conducted yielded gravity, magnetic, and electrical inversion parameters at elevations from 0 to 10000 meters with 200-meter intervals in the plane coordinate system. This also yielded the gridded inversion parameters W for gravity, magnetic, and electrical parameters at different elevations for each gridded measuring point. k (x i ,y i ,z i )(like Figure 3 , Figure 4 , Figure 5 As shown, a total of 262,701 density gridded inversion parameters from gravity exploration, 262,701 magnetic gridded inversion parameters from magnetic exploration, and 262,701 resistivity gridded inversion parameters from electrical exploration were obtained. When the planar measuring points of gravity, magnetic, and electrical exploration do not completely overlap, the measuring points overlapping by the three methods are selected for gridding (e.g., Figure 2 As shown, since the measurement points for gravity, magnetic, and electrical exploration completely overlap in this embodiment, all data from these points are used for subsequent calculations. To facilitate the calculation of subsequent comprehensive parameters, the gridded inversion parameters for gravity, magnetic, and electrical exploration need to be positiveized (assuming, for example...). Figure 3 , Figure 4 , Figure 5 The density, magnetic, and resistivity gridded inversion parameters shown are all greater than 0, therefore, as Figure 3 , Figure 4 , Figure 5 The gridded inversion parameters shown do not require positive value processing; where x i y i These are the planar coordinates of the measurement points after meshing (e.g.) Figure 2 As shown, x iThe value range of y is [0, 10000]. i The value range is [0, 10000], with intervals of 200 meters; z i The elevation corresponds to the gridded measuring point (as mentioned earlier, the elevation range is 0–10000 meters, with an elevation interval of 100 meters); where the subscript i is the number of the gridded measuring point, and i∈[1,n] (e.g., ...). Figure 2 As shown, i∈[1,51]); n is the total number of measurement points after gridding (e.g., Figure 2 As shown, n = 51); the superscript k is the gridded inversion parameter number of gravity, magnetic and electrical exploration, and k ∈ [1, 3]; where W 1 (x i ,y i ,z i ) represents the density gridded inversion parameters obtained from gravity exploration, W 2 (x i ,y i ,z i ) represents the magnetic gridding inversion parameters obtained from magnetic exploration, W 3 (x,y i ,z i The resistivity gridded inversion parameters are obtained from electrical exploration.

[0096] b) According to formula (1), the elevation z obtained in step a) is... i (As mentioned earlier, the selection is made from elevations between 0 and 10000 meters, specifically from 101 elevations: 10000 / 9900 / 9800 / 9700 / ... / 300 / 200 / 100 / 0 meters; assuming z is selected...) i The gridded inversion parameters of all measuring points (e.g., 10000 meters) are standardized to obtain the elevation z of all measuring points. i =10000 meters, the standardized parameters of this gridded inversion parameter Where p is the elevation z i The total number of a certain meshed inversion parameter (e.g.) Figure 2 As shown, p = 2601); k is the gridded inversion parameter number of gravity, magnetic and electrical exploration, and k∈[1,3]; that is, W 1 (x i ,y i ,z i ) represents the density gridded inversion parameters obtained from gravity exploration, W 2 (x i ,y i ,z i ) represents the magnetic gridding inversion parameters obtained from magnetic exploration, W 3 (x i ,y i ,zi The resistivity gridded inversion parameters are obtained from electrical exploration.

[0097] in

[0098]

[0099] Then obtain z i The standardized parameters of the gridded inversion parameters of density, magnetic field, and resistivity at an elevation of 10,000 meters are as follows:

[0100] c) Repeat step b) to calculate the standardized parameters of the gravity, magnetoelectric, and electrical gridded inversion parameters for other elevations until the standardized parameters of the gravity, magnetoelectric, and electrical gridded inversion parameters for all elevations are obtained; that is, the standardized parameters of gravity, magnetoelectric, and electrical gridded inversion parameters at different elevations are not calculated for each other, and the standardized parameters are only calculated for gravity, magnetoelectric, and electrical gridded inversion parameters at the same elevation (as mentioned above, if there are 101 different elevation values, then 101 sets of standardized parameters are calculated); in addition to the standardized parameters for the 10,000-meter elevation obtained in step b), 101 sets of standardized parameters can also be obtained for elevations of 9900 / 9800 / 9700 / … / 300 / 200 / 100 / 0 meters, for a total of 101 sets; which are as follows:

[0101]

[0102] d) Use formula (2) to analyze a certain elevation z i (As mentioned earlier, you can choose one of the 101 elevations mentioned above. Let's assume you choose elevation z.) i The correlation between the standardized parameters of the resistivity gridded inversion parameters (e.g., at elevation z = 10000 meters) and the standardized parameters of the density gridded inversion parameters or the magnetic gridded inversion parameters at the same elevation is calculated, and the elevation z is determined. i The correlation coefficient r between the standardized parameters of the resistivity gridded inversion parameters with a resistivity of 10000 and the standardized parameters of the density gridded inversion parameters or the magnetic gridded inversion parameters at the same elevation is given. 3,l,10000 N is the total number of measurement points participating in the correlation analysis of standardized parameters (e.g., ...). Figure 2 As shown, N = 2601); the subscript i is the number of the measurement point after gridding, and i ∈ [1, N] (as shown). Figure 2 As shown, i∈[1,2601]); The three-dimensional coordinates of the measuring point are (x i ,y i The standardized parameters of the resistivity gridded inversion parameters (10000) involved in the correlation analysis; The three-dimensional coordinates of the measuring point are (x i ,y iThe standardized parameter in the density and magnetic gridded inversion parameters involved in the correlation analysis is ,10000); the subscript l is the standardized parameter number of the density and magnetic gridded inversion parameters of gravity, magnetic and electric exploration, and l∈[1,2]; These are the standardized parameters for density gridded inversion parameters obtained from gravity exploration. These are the standardized parameters for the magnetic gridding inversion parameters obtained from magnetic exploration; where σ1≠0, σ2≠0, σ3≠0;

[0103]

[0104] in:

[0105]

[0106] e) Repeat step d) to obtain the correlation coefficients between the normalized parameters of the resistivity gridded inversion parameters at other elevations (as mentioned above, select one of the elevations of 9900 / 9800 / 9700 / ... / 300 / 200 / 100 / 0 meters) and the normalized parameters of the density gridded inversion parameters and the normalized parameters of the magnetic gridded inversion parameters at the same elevation. The correlation coefficients between the normalized parameters of the resistivity gridded inversion parameters for all elevations (as mentioned earlier, there are 101 elevations in total) and the normalized parameters of the density gridded inversion parameters and the normalized parameters of the magnetic gridded inversion parameters for the same elevations were obtained. Finally obtained Figure 6 , Figure 7 ;

[0107] f) Select an elevation z i (As mentioned earlier, you can choose one of the 101 elevations mentioned above. Let's assume you choose elevation z.) i The standardized parameters of the resistivity gridded inversion parameters (e.g., 10000 meters) are the primary standardized parameters, i.e. Main standardized parameters; The standardized parameters are the density grid inversion parameters for the corresponding elevation. The standardized parameters for the magnetic grid inversion parameters at the corresponding elevation are used; formula (3) is used to calculate the parameters for all measuring points and a certain elevation z. i The comprehensive parameters of the standardized parameters Where r 3,1,1000 For elevation z i Normalized parameters and elevation z of resistivity gridded inversion parameters = 10000 i The correlation coefficient between the standardized parameters of the density grid inversion parameters with a density of 10000. For elevation z i Normalized parameters and elevation z of resistivity gridded inversion parameters = 10000i The correlation coefficient between the standardized parameters of the magnetic gridded inversion parameters with a value of 10000;

[0108]

[0109] g) Repeat step f) to obtain the comprehensive parameters of standardized parameters for all elevations and all measuring points. (As mentioned earlier, standardized parameters for density, magnetic resistivity gridding can be obtained from 2601 measuring points and 101 elevations;) Finally, the following is obtained: Figure 8 , Figure 9 ;

[0110] Comprehensive parameters for all measuring points and all elevations (As mentioned earlier, the standardized parameters of density, magnetic, and resistivity gridded parameters obtained from 2601 measuring points and 101 elevations) are analyzed and interpreted to obtain interpretation results of gravity, magnetic, and electrical exploration of the hot dry rocks in the survey area; if a single-profile gravity, magnetic, and electrical exploration is carried out, two-dimensional contour lines are drawn for the comprehensive gravity, magnetic, and electrical parameters, and profile analysis and interpretation are conducted; if a multi-profile gravity, magnetic, and electrical exploration is carried out, three-dimensional contour lines are drawn for the comprehensive gravity, magnetic, and electrical parameters, and three-dimensional analysis and interpretation are conducted; if the hot dry rocks in the survey area generally exhibit low density, high magnetic properties, and high resistivity, then the overall comprehensive parameters will show high values ​​( Figure 9 This will primarily indicate a very strong correlation with hot dry rocks, interpreting the spatial locations of high-value comprehensive parameters as the presence of hot dry rocks, thus providing a more accurate identification effect for the exploration of hot dry rocks.

[0111] Beneficial effects: Improves the effectiveness of existing comprehensive exploration methods for hot dry rocks based on gravity, magnetics, and electricity, and enhances the speed and accuracy of hot dry rock identification. Reduces the difficulty of field data processing, data interpretation, and the arbitrariness of anomalies, providing a high-efficiency, low-cost, convenient, and rapid anomaly identification method and technology for comprehensive exploration of hot dry rocks.

[0112] The above specific embodiments further illustrate the purpose, technical solution, and beneficial effects of the present invention. It should be understood that the above are merely specific embodiments of the present invention and are not intended to limit the scope of protection of the present invention. Any modifications, equivalent substitutions, improvements, etc., made within the spirit and principles of the present invention should be included within the scope of protection of the present invention.

Claims

1. A method for exploring comprehensive gravity, magnetic, and electrical parameters for identifying hot dry rocks, characterized in that, The parameter exploration method includes: Step S1: Obtain the spatially corresponding gravity, magnetoelectric inversion parameters of the survey area, and perform gridding processing on the gravity, magnetoelectric inversion parameters; Step S2: Standardize the gridded inversion parameters; Step S3: Calculate the correlation coefficients between resistivity and density parameters and magnetic parameters, respectively, using resistivity as the main parameter; Step S4: Calculate the combined gravity, magnetism, and electrical parameters; Step S5: Analyze and interpret the combined gravity, magnetism, and electrical parameters.

2. The method for exploring comprehensive gravity, magnetic, and electrical parameters for identifying hot dry rocks according to claim 1, characterized in that, Step S1: Obtaining the spatially corresponding gravity, magnetoelectric inversion parameters of the survey area, and performing gridding processing on the gravity, magnetoelectric inversion parameters specifically includes: In the survey area where comprehensive gravity, magnetic and electrical exploration of hot dry rocks is carried out, the grid parameters for gravity, magnetic and electrical exploration are designed according to the requirements of comprehensive exploration work, comprehensive gravity, magnetic and electrical exploration work is carried out, the gravity, magnetic and electrical parameters of multiple measuring points are obtained, and the relevant parameters are inverted and calculated. Obtain gravity, magnetic and electrical inversion parameters for each measuring point at different elevations, and perform gridded calculations on the gravity, magnetic and electrical inversion parameters; Obtain the gridded inversion parameters W of gravity, magnetism, and electricity at different elevations for each gridded measuring point. k (x i ,y i ,z i ); When the planar measuring points of gravity, magnetic and electrical exploration do not completely overlap, the measuring points that overlap using three different methods are selected for gridding. Where x i y i is the plane coordinate of the measuring point after gridding; z_i is the elevation corresponding to the gridded measuring point; where the subscript i is the number of the measuring point after gridding, and i∈[1,n]; n is the total number of measuring points after gridding, where n is a natural number not less than 1; the superscript k is the gridding inversion parameter number of gravity, magnetic and electrical exploration, and k∈[1,3]; where W 1 (x i ,y i ,z i ) represents the density gridded inversion parameters obtained from gravity exploration, W 2 (x i ,y i ,z i ) represents the magnetic gridding inversion parameters obtained from magnetic exploration, W 3 (x i ,y i ,z i ) represents the resistivity gridded inversion parameters obtained from electrical exploration.

3. The method for exploring comprehensive gravity, magnetic, and electrical parameters for identifying hot dry rocks according to claim 2, characterized in that, The inversion calculation of relevant parameters specifically includes: inverting the gravity, magnetic and electrical parameters obtained from gravity, magnetic and electrical exploration to obtain inversion results of relevant density, magnetic parameters and resistivity parameters.

4. The method for exploring comprehensive gravity, magnetic, and electrical parameters for identifying hot dry rocks according to claim 2, characterized in that, The specific steps for obtaining gravity, magnetoelectric inversion parameters for each measuring point at different elevations include obtaining relevant density, magnetic parameters, and resistivity parameters for different measuring points at different depths.

5. The method for exploring comprehensive gravity, magnetic, and electrical parameters for identifying hot dry rocks according to claim 2, characterized in that, When the planar measuring points of gravity, magnetic, and electrical exploration do not completely overlap, the overlapping measuring points are selected using three methods for gridding, specifically including: When the minimum value of a certain gridded inversion parameter of gravity, magnetoelectricity is negative, then add twice the absolute value of the minimum value of the gridded parameter to all such gridded parameters. When the minimum value of a certain gridded inversion parameter of gravity, magnetoelectricity is zero, the absolute value of the maximum value of the inversion parameter is added to all such inversion parameters. If a certain gridded inversion parameter of gravity, magnetoelectricity is 0, then no further calculation of that parameter will be performed.

6. The method for exploring comprehensive gravity, magnetic, and electrical parameters for identifying hot dry rocks according to claim 1, characterized in that, Step S2: Standardizing the gridded inversion parameters specifically includes: According to formula (1), the obtained elevation z i The gridded inversion parameters of all measuring points are standardized to unify all data into the interval [0,1], thus obtaining the elevation z of all measuring points. i Standardized parameters of gridded inversion parameters Where p is the elevation z i The total number of a certain gridded inversion parameter; k is the gridded inversion parameter number of gravity, magnetic and electrical exploration, and k∈[1,3]; that is, W 1 (x i ,y i ,z i ) represents the density gridded inversion parameters obtained from gravity exploration, W 2 (x i ,y i ,z i ) represents the magnetic gridding inversion parameters obtained from magnetic exploration, W 3 (x i ,y i ,z i The resistivity gridded inversion parameters are obtained from electrical exploration. in Repeat the previous step to calculate the standardized parameters of the gravity, magnetoelectric and electrical gridded inversion parameters for the remaining elevations, until the standardized parameters of the gravity, magnetoelectric and electrical gridded inversion parameters for all elevations are obtained.

7. The method for exploring comprehensive gravity, magnetic, and electrical parameters for identifying hot dry rocks according to claim 1, characterized in that, Step S3, which involves calculating the correlation coefficients between resistivity and density and magnetic parameters, specifically includes: Formula (2) is used to analyze a certain elevation z i The correlation between the normalized parameters of the resistivity gridded inversion parameters and the normalized parameters of the density gridded inversion parameters or the normalized parameters of the magnetic gridded inversion parameters at the same elevation; Determine the elevation z i The correlation coefficients between the normalized parameters of the resistivity gridded inversion parameters and the normalized parameters of the density gridded inversion parameters or the normalized parameters of the magnetic gridded inversion parameters at the same elevation. Where N is the total number of measurement points participating in the correlation analysis of standardized parameters; the subscript i is the number of the measurement point after gridding, and i∈[1,N]; The three-dimensional coordinates of the measuring point are (x i ,y i ,z i The standardized parameters of resistivity gridded inversion parameters involved in correlation analysis; The three-dimensional coordinates of the measuring point are (x i ,y i ,z i The standardized parameter in the density and magnetic gridded inversion parameters involved in the correlation analysis is a certain parameter; the subscript l is the standardized parameter number of the density and magnetic gridded inversion parameters of gravity, magnetic and electric exploration, and l∈[1,2]. These are the standardized parameters for density gridded inversion parameters obtained from gravity exploration. These are the standardized parameters for magnetic gridding inversion parameters obtained from magnetic exploration. Where σ1≠0, σ2≠0, σ3≠0; in: Repeat the previous step to obtain the correlation coefficients between the standardized parameters of the resistivity gridded inversion parameters at the remaining elevation and the standardized parameters of the density gridded inversion parameters and magnetic gridded inversion parameters at the same elevation. The correlation coefficients between the normalized parameters of the resistivity gridded inversion parameters at all elevations and the normalized parameters of the density gridded inversion parameters and the normalized parameters of the magnetic gridded inversion parameters at the same elevations were obtained.

8. The method for exploring comprehensive gravity, magnetic, and electrical parameters for identifying hot dry rocks according to claim 1, characterized in that, Step S4: Calculating the combined gravity, magnetic, and electrical parameters specifically includes: Select a certain elevation z i The standardized parameters of the resistivity gridded inversion parameters are the primary standardized parameters; Formula (3) is used to calculate all measuring points and a certain elevation z. i The comprehensive parameters of the standardized parameters Repeat the previous step to obtain the comprehensive parameters of standardized parameters for all elevations and all measuring points.

9. A method for exploring comprehensive gravity, magnetic, and electrical parameters for identifying hot dry rocks according to claim 8, characterized in that, The selection of a certain elevation z i The standardized parameters of the resistivity gridded inversion parameters are the main standardized parameters, which specifically include: Main standardized parameters; The standardized parameters are the density grid inversion parameters for the corresponding elevation. These are the standardized parameters for the magnetic grid inversion parameters of the corresponding elevation.

10. The method for exploring comprehensive gravity, magnetic, and electrical parameters for identifying hot dry rocks according to claim 8, characterized in that, The formula (3) is used to calculate all measuring points and a certain elevation z. i The comprehensive parameters of the standardized parameters specific include: in, For elevation z i The normalized parameters of the resistivity gridded inversion parameters and the elevation z i The correlation coefficients between the standardized parameters of the density gridded inversion parameters. For elevation z i The normalized parameters of the resistivity gridded inversion parameters and the elevation z i The correlation coefficients between the standardized parameters of the magnetic gridded inversion parameters; 11. The method for exploring comprehensive gravity, magnetic, and electrical parameters for identifying hot dry rocks according to claim 1, characterized in that, Step S5: Analysis and interpretation of the combined gravity, magnetism, and electrical parameters specifically includes: Comprehensive parameters for all measuring points and all elevations Analysis and interpretation were conducted to obtain the interpretation results of gravity, magnetic and electrical exploration of dry hot rocks in the survey area; If a single profile exploration of gravity, magnetic, and electrical parameters is carried out, then two-dimensional contour lines are drawn for the comprehensive parameters of gravity, magnetic, and electrical parameters, and profile analysis and interpretation are performed. If a multi-profile exploration of gravity, magnetoelectricity is carried out, then three-dimensional contour lines are drawn for the comprehensive parameters of gravity, magnetoelectricity, and electricity, and three-dimensional analysis and interpretation are performed.

Images