Method, device and equipment for predicting high-ground-temperature abnormal area through gravity and magnetic data and medium

Abstract

The invention provides a method, a device, equipment and a medium for predicting a high-ground-temperature abnormal area by gravity and magnetic data, and the method comprises the following steps: carrying out the geological interpretation of gravity data, and explaining the geological significance of gravity data abnormity; calculating gravity average data according to the gravity data, determining a gravity gradient zone by using the gravity average data, and analyzing a tectonic fracture zone according to the gravity gradient zone and a regional tectonic evolution event; geologically explaining the magnetic data, and explaining the geological significance of the magnetic data anomaly; calculating residual magnetic force data according to the magnetic force data, and determining a magnetic force abnormal area by utilizing the residual magnetic force data and the structural fracture zone; and determining a distribution range of a magma invasion area by using the magnetic force abnormal area and the area structure evolution event, and predicting a high ground temperature abnormal area according to the distribution range of the magma invasion body and ground temperature limit data.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of geothermal and geophysics, and particularly relates to a method and device for predicting a high-geothermal-anomaly area based on gravity and magnetic data, equipment and a medium. BACKGROUND

[0002] The geothermal field characteristics of a sedimentary basin are a comprehensive response of tectonic deformation of the lithosphere and regional tectonic evolution of the basin, and are mainly affected by factors such as magmatic activity and basin structure. The study of the geothermal field characteristics on a basin scale is of great scientific significance to deepening the understanding of the dynamic process of basin formation.

[0003] At present, the study of the geothermal field characteristics of a sedimentary basin mainly adopts two methods of borehole temperature measurement and calculation of terrestrial heat flow. The borehole temperature measurement method is to measure the temperature of drilling fluid to display the original temperature of the underground rock formation. The calculation of terrestrial heat flow needs to obtain the geothermal gradient value and the rock thermal conductivity in the corresponding depth section. However, both methods need to be supported by drilling data and cannot be applied to sedimentary basins with low exploration degree and lack of drilling data. SUMMARY

[0004] In view of the above problems, the present application provides a method and device for predicting a high-geothermal-anomaly area based on gravity and magnetic data, equipment and a medium.

[0005] In a first aspect, the present application provides a method for predicting a high-geothermal-anomaly area based on gravity and magnetic data, characterized in that the method comprises the following steps:

[0006] S1, geologically interpreting gravity data and explaining the geological significance of gravity data anomalies; calculating gravity average data from the gravity data, determining a gravity gradient zone by using the gravity average data, and analyzing a tectonic fracture zone according to the gravity gradient zone and regional tectonic evolution events;

[0007] S2, geologically interpreting magnetic data and explaining the geological significance of magnetic data anomalies; calculating magnetic residual data from the magnetic data, and determining a magnetic anomaly area by using the magnetic residual data and the tectonic fracture zone;

[0008] S3, determining the distribution range of a magmatic intrusion area by using the magnetic anomaly area and the regional tectonic evolution events, and predicting a high-geothermal-anomaly area according to the distribution range of the magmatic intrusion body and geothermal limiting data.

[0009] In some embodiments, the step S1 further comprises, before the step S1:

[0010] collecting regional tectonic evolution information and extracting the regional tectonic evolution events;

[0011] The regional tectonic evolution events include regional tectonic compression events and / or regional tectonic stretching events.

[0012] In some embodiments, the tectonic fracture zone is parsed from the regional tectonic evolution events according to the gravity gradient zone and regional tectonic evolution event, comprising:

[0013] Obtaining gravity gradient characteristic parameters of the gravity gradient zone; wherein the gravity gradient characteristic parameters comprise at least one of gravity gradient anomaly value, gravity gradient vector, gravity gradient modulus value, and gravity gradient vector rotation angle;

[0014] Parsing the tectonic fracture zone from the regional tectonic evolution events according to the gravity gradient characteristic parameters.

[0015] In some embodiments, the gravity average data is calculated according to the gravity data, comprising:

[0016] The gravity average data is calculated according to the gravity data based on a sliding average method.

[0017] In some embodiments, the magnetic force residual data is calculated according to the magnetic force data, comprising:

[0018] The magnetic force average data is calculated according to the magnetic force data based on a sliding average method;

[0019] The magnetic force residual data is determined according to the difference between the magnetic force data and the magnetic force average data.

[0020] In some embodiments, the magnetic force anomaly zone is determined by using the magnetic force residual data and the tectonic fracture zone, comprising:

[0021] The correlation between the magnetic force residual data and the tectonic fracture zone is established according to the geological interpretation of the magnetic force residual data;

[0022] The magnetic force anomaly zone is determined according to the correlation between the magnetic force residual data and the tectonic fracture zone.

[0023] In some embodiments, the distribution range of the magmatic intrusion zone is determined by using the magnetic force anomaly zone and the regional tectonic evolution events, and a high geothermal anomaly zone is predicted according to the distribution range of the magmatic intrusion body and geothermal limit data, comprising:

[0024] Obtaining magnetic force anomaly characteristic parameters of the magnetic force anomaly zone; wherein the magnetic force anomaly characteristic parameters comprise at least one of magnetic field strength anomaly value, magnetic deflection angle, magnetic gradient vector rotation angle, and magnetic polarity anomaly value;

[0025] The magmatic intrusion zone and its distribution range are determined from the regional tectonic evolution events according to the magnetic force anomaly characteristic parameters;

[0026] According to the distribution range of the magma intrusive body and the ground temperature limit data, a high ground temperature anomaly area is predicted.

[0027] In a second aspect, the application further provides a device for predicting a high ground temperature anomaly area based on gravity and magnetic data, characterized in that the device comprises:

[0028] A tectonic fracture zone analysis module is configured to perform geological interpretation on gravity data and explain the geological significance of the gravity data anomaly; gravity average data is calculated according to the gravity data, gravity gradient zones are determined using the gravity average data, and tectonic fracture zones are analyzed according to the gravity gradient zones and regional tectonic evolution information.

[0029] A magnetic anomaly area determination module is configured to perform geological interpretation on magnetic data and explain the geological significance of the magnetic data anomaly; magnetic residual data is calculated according to the magnetic data, and magnetic anomaly areas are determined using the magnetic residual data and the tectonic fracture zones.

[0030] A high ground temperature anomaly area prediction module is configured to determine the distribution range of a magma intrusive area using the magnetic anomaly areas and the regional tectonic evolution events, and to predict a high ground temperature anomaly area according to the distribution range of the magma intrusive body and the ground temperature limit data.

[0031] In a third aspect, the application further provides an electronic device comprising a memory and a processor, wherein the memory stores a computer program,

[0032] The processor executes the computer program to implement the steps of the method according to any one of the above embodiments.

[0033] In a fourth aspect, the application further provides a computer readable storage medium, characterized in that the computer readable storage medium stores at least one instruction, and the at least one instruction is executed by a processor to implement the steps of the method according to any one of the above embodiments.

[0034] The application has the following beneficial effects:

[0035] The application provides a method for predicting a high ground temperature anomaly area based on gravity and magnetic data. Based on the analysis of the geological background such as basin tectonic evolution events and basin structure, the tectonic fracture zones of the basin structure are analyzed using the gravity field characteristics (gravity data, gravity average data, gravity gradient zones), the magnetic anomaly areas and the distribution range of the magma intrusive body are analyzed using the magnetic field characteristics (magnetic data, magnetic residual data), and then the high ground temperature anomaly area of the basin is predicted by using the set ground temperature data as a constraint, and the ground temperature field characteristic analysis of the entire basin is realized. The application can not only break through the constraint of the lack of drilling data in low exploration basins, but also has important significance for the evaluation of the basin geothermal resources and the oil and gas potential of the system. BRIEF DESCRIPTION OF DRAWINGS

[0036] Figure 1 A flowchart of a method for predicting high-temperature anomaly areas based on gravity and magnetic data is provided for the embodiments of the present application.

[0037] Figure 2 An EGM2008 original free-air gravity map of the Powder River Basin is provided for the embodiments of the present application.

[0038] Figure 3 A 49*49 window sliding average field map of the EGM2008 free-air gravity data of the Powder River Basin is provided for the embodiments of the present application.

[0039] Figure 4 An EGM2008 original magnetic map of the Powder River Basin is provided for the embodiments of the present application.

[0040] Figure 5 A 39*39 window sliding average field map of the EGM2008 free-air gravity data of the Powder River Basin is provided for the embodiments of the present application.

[0041] Figure 6 A structural block diagram of a device for predicting high-temperature anomaly areas based on gravity and magnetic data is provided for the embodiments of the present application.

[0042] Figure 7 A structural block diagram of an electronic device is provided for the embodiments of the present application. DETAILED DESCRIPTION

[0043] In order to make the objectives, technical solutions and advantages of the present application clearer, further detailed description will be given to the present application in combination with the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are only used to explain the present application and are not used to limit the present application.

[0044] Therefore, the detailed description of the embodiments of the present application provided in the accompanying drawings is not intended to limit the scope of the claimed present application, but only represents selected embodiments of the present application. All other embodiments obtained by those of ordinary skill in the art based on the embodiments in the present application without making creative efforts are within the scope of protection of the present application.

[0045] Wherein, in the present application, relational terms such as first and second and the like can merely be used to distinguish one entity or action from another, without necessarily requiring or implying that these entities or actions are in any way mutually exclusive, directly or indirectly, unless the context clearly suggests otherwise. Moreover, the terms "comprises", "comprising", or any other variation thereof, are intended to cover a non-exclusive inclusion, such that a process, method, article, or apparatus that comprises a list of elements does not include only those elements but can include other elements not expressly listed or inherent to such process, method, article, or apparatus. An element proceeded by "comprises a... " does not, without more constraints, exclude the presence of additional identical elements in the process, method, article, or apparatus that comprises the element.

[0046] Some embodiments of the present application will be described in detail with reference to the drawings, in which some embodiments of the present application are shown. The following embodiments and features of the embodiments can be combined with each other, without conflict.

[0047] At present, the research on the characteristics of the geothermal field of a sedimentary basin mainly adopts two methods of borehole temperature measurement and geothermal flow calculation. Among them, the borehole temperature measurement method is to measure the temperature of the drilling fluid to display the original temperature of the underground rock formation. The calculation of geothermal flow needs to obtain the geothermal gradient value and the rock thermal conductivity in the corresponding depth section. However, both methods need to be supported by drilling data and cannot be applied to sedimentary basins with low exploration degree and lack of drilling data.

[0048] Therefore, the embodiment of the present application provides a method for predicting high geothermal anomaly areas by gravity and magnetic data to solve at least one of the problems existing in the prior art.

[0049] Embodiment one

[0050] Figure 1 A flowchart of a method for predicting high geothermal anomaly areas by gravity and magnetic data provided by the embodiment of the present application is shown in Figure 1 The method for predicting high geothermal anomaly areas by gravity and magnetic data comprises the following steps:

[0051] S1, geologically interpret the gravity data, and interpret the geological significance of the gravity data anomaly; calculate the gravity average data according to the gravity data, determine the gravity gradient zone by using the gravity average data, and analyze the tectonic fracture zone according to the gravity gradient zone and regional tectonic evolution events.

[0052] Further, as shown in Figure 2 The gravity data can be original free air gravity data.

[0053] The original free-air gravity data is obtained by downloading the EGM2008 original free-air gravity data of the Powder River Basin with a grid spacing of 0.01' x 0.01' from the International Center for Global Earth Models (ICGEM) website.

[0054] It should be noted that the geological interpretation of gravity data is a common geophysical exploration method. For example, the geological interpretation of gravity data can help understand the density of underground rocks, geological structure and geological characteristics. Among them, the gravity data anomaly can reflect the density change of underground rocks, the geological structure of underground rocks and the different anomalies of the geological characteristics of underground rocks in the gravity field.

[0055] In one example, the density of underground rocks is different, which will cause the change of gravity field. For example, high-density rocks (such as granite) will produce positive gravity data anomaly, while low-density rocks (such as volcanic rocks) will cause negative gravity data anomaly. By analyzing the gravity data anomaly, the density change of underground rocks can be inferred, and then the geological structure and rock type can be inferred.

[0056] In another example, the geological structure of underground rocks will cause the change of gravity field. For example, the geological structure such as fault zone, fold, magmatic intrusion, etc. shows gravity data anomaly in gravity data. Through the analysis of gravity anomaly data, the spatial distribution and shape of gravity anomaly can be obtained, so as to identify the location, shape and nature of the geological structure of underground rocks.

[0057] In the embodiment of the present application, taking the Powder River Basin as an example, the original free-air gravity data of the Powder River Basin is geologically interpreted, and the middle part of the Powder River Basin is a low gravity anomaly, indicating that the basin basement is deep, the southern part of the Powder River Basin shows a weak low gravity anomaly, indicating that the basin basement is relatively shallow compared with the western part of the basin, and the eastern part of the Powder River Basin is a high gravity anomaly, indicating that the basin basement is shallow.

[0058] In one example, the gravity average data is calculated according to the gravity data, and the gravity gradient zone is determined by using the gravity average data. The gravity average data is obtained by filtering or smoothing the gravity data, which can reduce the noise and sudden change in the gravity data. Therefore, the gravity average data is smoother than the gravity data, and it is easier to determine the gravity gradient zone, thereby helping to reduce random errors in the data.

[0059] In addition, the tectonic fracture zone can be analyzed by analyzing the spatial distribution and shape of the gravity gradient zone and combining the geological tectonic evolution history of the Powder River Basin to infer the location, nature and evolution process of the tectonic fracture zone.

[0060] Tectonic fracture zone is usually a manifestation of tectonic activity zone, which can be a fault zone, shear zone and other geological structure features, and has important indicative significance for the geological evolution and structural features of the pocono basin. Through comprehensive analysis of gravity data and tectonic evolution events, the formation mechanism and geological significance of the tectonic fracture zone can be better understood.

[0061] S2, geologically interpreting the magnetic force data and explaining the geological significance of the magnetic force data anomaly; calculating the magnetic force residual data according to the magnetic force data, and determining the magnetic anomaly area by using the magnetic force residual data and the tectonic fracture zone.

[0062] For example, in the embodiment of the present application, according to the magnetic force residual data and the tectonic fracture zone, it can be determined that there is a high magnetic anomaly area extending from the north of the pocono basin to the northwest, and there is a high magnetic anomaly area extending from the middle of the pocono basin to the northeast.

[0063] It should be noted that the magnetic force data is an important non-destructive detection means in geological interpretation. After the geological interpretation of the magnetic force data, information about the properties of underground rocks, geological structures, geological features, etc. can be obtained by analyzing the magnetic force data. For example, in an example, the magnetic anomaly data in the magnetic force data can reflect the magnetic difference of the underground rocks, such as rock type, mineral content, etc. By analyzing the size, shape, direction and other characteristics of the anomaly, it can be predicted that there may be a magnetic anomaly area, such as a fault zone, a fold, a magmatic intrusion, etc.

[0064] In the embodiment of the present application, the magmatic intrusion is taken as an example for illustration. When the magma intrudes underground, it will change the magnetic characteristics of the surrounding rocks. This magnetic change will cause disturbance of the underground magnetic field, thereby forming a magnetic anomaly area.

[0065] S3, determining the distribution range of the magmatic intrusion area by using the magnetic anomaly area and the regional tectonic evolution event, and predicting the high geothermal anomaly area according to the distribution range of the magmatic intrusion and the geothermal limit data.

[0066] The present application proposes a method for predicting a high geothermal anomaly area based on the analysis of regional tectonic evolution events and basin structure, etc. geological background, using gravity field characteristics (gravity data, gravity average data, gravity gradient zone) to analyze the tectonic fracture zone of the basin structure, using magnetic field characteristics (magnetic force data, magnetic force residual data) to analyze the magnetic anomaly area and the distribution range of the magmatic intrusion, and then using limited geothermal data as a constraint to predict the high geothermal anomaly area of the basin, and realize the analysis of the geothermal field characteristics of the whole basin. The present application can not only break through the constraint of lacking drilling data in low exploration basin, but also has important significance for evaluating the geothermal resources and oil and gas potential of the system.

[0067] In some embodiments, step S1 further comprises:

[0068] The regional tectonic evolution information is collected to extract the regional tectonic evolution events. The regional tectonic evolution events include regional tectonic compression events and / or regional tectonic stretching events.

[0069] In one example, the regional tectonic compression events include plate collision, orogenic movement, rock deformation, tectonic folding, and seismic activity, etc. These regional tectonic compression events are common geological processes on the Earth's surface and crust, reflecting the complexity and diversity of the Earth's tectonic evolution. Studying these regional tectonic compression events helps to understand the internal structure of the region, tectonic movement, and the distribution and formation of geological resources.

[0070] In another example, the regional tectonic stretching events include plate expansion, volcanic eruption, fault zone formation, mantle column uplift, etc. These regional tectonic stretching events are common geological processes on the Earth's surface and crust, reflecting the complexity and diversity of the Earth's tectonic evolution. Studying these regional tectonic stretching events helps to understand the internal structure of the region, tectonic movement, and the distribution and formation of geological resources.

[0071] In the embodiments of the present application, taking the Po River Basin as an example, the tectonic evolution information of the Po River Basin is collected, and two sets of basement fault evolution events in the southeast and north-west directions of the Po River Basin can be obtained, which can provide a channel for the upwelling of mantle material and heat flow.

[0072] In some embodiments, the tectonic fracture zone is analyzed according to the gravity gradient zone and the regional tectonic evolution event, comprising:

[0073] The gravity gradient characteristic parameters of the gravity gradient zone are obtained. The gravity gradient characteristic parameters include at least one of the gravity gradient anomaly value, the gravity gradient vector, the gravity gradient modulus value, and the gravity gradient vector rotation angle.

[0074] The tectonic fracture zone is analyzed from the regional tectonic evolution event according to the gravity gradient characteristic parameters.

[0075] In the embodiments of the present application, taking the Po River Basin as an example, two gravity gradient zones in the northeast and northwest directions of the basin can be identified according to the gravity average data, wherein the gravity gradient zone in the northwest direction is faulted and broken by the gravity gradient zone in the northeast direction. Combined with the geological background of the Po River Basin, it can be determined that there are two tectonic fracture zones in the northeast and northwest directions of the Po River Basin.

[0076] In some embodiments, the gravity gradient zone includes a vertical gradient zone and a horizontal gradient zone. For example, in one example, the vertical gradient zone of the gravity field in three-dimensional space can be calculated by the following formula:

[0077]

[0078] wherein g(z) is a function of gravity average data at height z, h is a height variation.

[0079] In some embodiments, the step of calculating gravity average data from the gravity data comprises:

[0080] calculating the gravity average data from the gravity data based on a sliding average method.

[0081] For example, in one example, the step of calculating the gravity average data from the gravity data based on a sliding average method comprises:

[0082] Step 1, determining the size of the sliding average window;

[0083] Step 2, selecting a starting point as the starting position of the sliding average, and starting from the starting point, calculating the gravity average data of the center point in the window according to the size of the selected window and the gravity data in the window.

[0084] For example, in one example, for original gravity data of m columns and n rows, line distance Δx, and point distance Δy, the size of the sliding average window is set to (m*Δx)*(n*Δy), wherein m is the number of lines in the x direction, and n is the number of points in the y direction. The calculation formula of the gravity average data (i.e., the average value of the gravity data) in the window is:

[0085]

[0086] g ij is the gravity data of the ij point in the (m*Δx)*(n*Δy) window, is the gravity average data of the center point in the (m*Δx)*(n*Δy) window, m>0, n>0.

[0087] In some embodiments, the step of calculating magnetic force residual data from the magnetic force data comprises:

[0088] calculating the magnetic force average data from the magnetic force data based on a sliding average method; and determining the magnetic force residual data according to the difference between the magnetic force data and the magnetic force average data.

[0089] Further, as shown in Figure 4 the magnetic force data can be original magnetic force data.

[0090] The magnetic force data is obtained by downloading EGM2008 original magnetic force data of the Powder Basin with a grid interval of 0.01'x0.01' from an International Center for Global Earth Models (ICGEM) website.

[0091] For example, in one example, the magnetic force average data is calculated according to the magnetic force data based on a sliding average method, including:

[0092] Step 1, determining the size of a sliding average window;

[0093] Step 2, selecting a starting point as a starting position of the sliding average, and starting from the starting point, sequentially calculating the magnetic force average data of a center point in a window according to the size of the selected window and the magnetic force data in the window.

[0094] For example, in one example, for original magnetic force data with m columns and n rows, a line interval of Δx and a point interval of Δy, the window size of the sliding average is set to (m*Δx)*(n*Δy), m is the number of lines in the x direction, and n is the number of points in the y direction. The calculation formula of the magnetic force average data in the window is:

[0095]

[0096] T ij is the magnetic force data of the ij point in the (m*Δx)*(n*Δy) window, is the magnetic force average data of the center point in the (m*Δx)*(n*Δy) window, m>0, n>0.

[0097] In another example, the magnetic force residual data is obtained according to the magnetic force data and the magnetic force average data, including the step of: subtracting the magnetic force average data from the magnetic force data to obtain the magnetic force residual data.

[0098] For example, further, in some embodiments, the magnetic force data, the magnetic force average data and the magnetic force residual data are interpreted in terms of geology.

[0099] wherein the magnetic force residual data is denoted as ΔT is denoted as a magnetic force residual field, ΔT kl is denoted as a magnetic force data field, is denoted as a magnetic force average field. The formula for obtaining the magnetic force average field is:

[0100]

[0101] In some embodiments, the magnetic force anomaly area is determined by using the magnetic force residual data and the tectonic fracture zone, including: establishing the correlation between the magnetic force residual data and the tectonic fracture zone according to the geological interpretation of the magnetic force residual data; and determining the magnetic force anomaly area according to the correlation between the magnetic force residual data and the tectonic fracture zone.

[0102] In some embodiments, the distribution range of the magma intrusion area is determined by using the magnetic anomaly area and regional tectonic evolution event, and the high geothermal anomaly area is predicted according to the distribution range of the magma intrusion body and the geothermal limit data, comprising:

[0103] The magnetic anomaly characteristic parameters of the magnetic anomaly area are acquired; wherein the magnetic anomaly characteristic parameters include at least one parameter in magnetic field strength anomaly value, magnetic deflection angle, magnetic gradient vector rotation angle and magnetic polarity anomaly value.

[0104] The magma intrusion area is determined from the regional tectonic evolution event according to the magnetic anomaly characteristic parameters, and the high geothermal anomaly area is predicted by using the distribution range of the magma intrusion body.

[0105] The present application can not only break through the constraint of lacking drilling data in the low exploration area, but also has important significance for the system to evaluate the geothermal resources and oil and gas potential in the low exploration area.

[0106] Embodiment two

[0107] Based on the foregoing description, the present application will be further described in combination with specific scenarios:

[0108] Figure 2 The EGM2008 original free air gravity map of the Powder River Basin provided by the embodiments of the present application, Figure 3 The 49*49 window sliding average field map of the EMG2008 free air gravity data of the Powder River Basin provided by the embodiments of the present application, Figure 4 The EGM2008 original magnetic map of the Powder River Basin provided by the embodiments of the present application, Figure 5 The 39*39 window sliding average field map of the EMG2008 free air gravity data of the Powder River Basin provided by the embodiments of the present application, as Figure 1 The method for predicting a high geothermal anomaly area by using gravity and magnetic data provided by the embodiments of the present application comprises the following steps:

[0109] S1, geological interpretation is performed on gravity data, and the geological significance of gravity data anomaly is interpreted; gravity average data is calculated according to the gravity data, gravity gradient zones are determined by using the gravity average data, and tectonic fracture zones are analyzed according to the gravity gradient zones and regional tectonic evolution events.

[0110] Further, as Figure 2 The gravity data can be original free air gravity data.

[0111] The acquisition method of the original free air gravity data comprises: downloading the EGM2008 original free air gravity data of the Powder River Basin with a grid interval of 0.01'x0.01' from the international global earth model center (ICGEM) website.

[0112] It should be noted that the geological interpretation of gravity data is a common geophysical exploration method. For example, the geological interpretation of gravity data can help understand the density of underground rocks, geological structure and geological characteristics. Among them, the gravity data anomaly can reflect the density change of underground rocks, the geological structure of underground rocks and the different anomalies of underground rocks in the gravity field.

[0113] In one example, the density difference of underground rocks causes the change of gravity field. For example, high-density rocks (such as granite) will produce positive gravity data anomalies, while low-density rocks (such as volcanic rocks) will cause negative gravity data anomalies. By analyzing the gravity data anomaly, the density change of underground rocks can be inferred, and then the geological structure and rock type can be inferred.

[0114] In another example, the geological structure of underground rocks will cause the change of gravity field. For example, the geological structure such as fault zone, fold, magmatic intrusion, etc. shows gravity data anomaly in gravity data. Through the analysis of gravity anomaly data, the spatial distribution and shape of gravity anomaly can be obtained, so as to identify the location, shape and nature of the geological structure of underground rocks.

[0115] In the embodiment of the present application, taking the Powder River Basin as an example, the original free air gravity data of the Powder River Basin is interpreted geologically, and the middle part of the Powder River Basin is a low gravity anomaly, indicating that the basin basement is deep, the southern part of the Powder River Basin shows a weak low gravity anomaly, indicating that the basin basement is relatively shallow compared with the western part of the basin, and the eastern part of the Powder River Basin is a high gravity anomaly, indicating that the basin basement is shallow.

[0116] In one example, the gravity average data is calculated according to the gravity data, and the gravity gradient zone is determined by using the gravity average data. The gravity average data is obtained by filtering or smoothing the gravity data, which can reduce the noise and sudden change in the gravity data. Therefore, the gravity average data is smoother than the gravity data, and it is easier to determine the gravity gradient zone, thereby helping to reduce random errors in the data.

[0117] In addition, the tectonic fracture zone can be analyzed by analyzing the spatial distribution and shape of the gravity gradient zone and the regional tectonic evolution event. Combined with the geological tectonic evolution history of the Powder River Basin, the location, nature and evolution process of the tectonic fracture zone can be inferred.

[0118] The tectonic fracture zone is usually a manifestation of the geological tectonic activity zone, which may be a fault zone, a shear zone and other geological structure characteristics, and has important indicative significance for the geological evolution and tectonic characteristics of the Powder River Basin. By comprehensively analyzing the gravity data and tectonic evolution events, the formation mechanism and geological significance of the tectonic fracture zone can be better understood.

[0119] S2, geologically interpret the magnetic force data, and interpret the geological significance of the magnetic force data anomaly; calculate magnetic force residual data according to the magnetic force data, and determine the magnetic force anomaly area by using the magnetic force residual data and the tectonic fracture zone.

[0120] For example, in the embodiment of the present application, according to the magnetic force residual data and the tectonic fracture zone, it can be determined that there is a high magnetic force anomaly area extending from the north of the basin to the northwest direction in the Powell basin, and there is a high magnetic force anomaly area extending from the middle of the basin to the northeast direction in the Powell basin.

[0121] It should be noted that the magnetic force data is an important non-destructive detection means in geological interpretation. After the geological interpretation of the magnetic force data, information about the properties of underground rocks, geological structure, geological characteristics, etc. can be obtained by analyzing the magnetic force data. For example, in an example, the magnetic force anomaly data in the magnetic force data can reflect the magnetic difference of the underground rocks, such as rock type, mineral content, etc. By analyzing the characteristics such as size, shape, direction, etc. of the anomaly, it can be predicted that there may be a magnetic force anomaly area, such as a fracture zone, a fold, a magmatic intrusion, etc.

[0122] In the embodiment of the present application, the magmatic intrusion is taken as an example for illustration. When the magma intrudes underground, it will change the magnetic characteristics of the surrounding rocks. This magnetic change will cause disturbance of the underground magnetic field, thereby forming a magnetic force anomaly area.

[0123] S3, determine the distribution range of the magmatic intrusion area by using the magnetic force anomaly area and the regional tectonic evolution event, and predict the high geothermal anomaly area according to the distribution range of the magmatic intrusion and the geothermal limit data.

[0124] The present application proposes a method for predicting a high geothermal anomaly area from gravity and magnetic data. Based on the analysis of regional tectonic evolution events and basin structure, etc. geological background, the method uses gravity field characteristics (gravity data, gravity average data, gravity gradient zone) to analyze the tectonic fracture zone of the basin structure, uses magnetic field characteristics (magnetic force data, magnetic force residual data) to analyze the magnetic force anomaly area and the distribution range of the magmatic intrusion, and then uses limited geothermal data as a constraint to predict the high geothermal anomaly area of the basin, and realizes the analysis of the geothermal field characteristics of the whole basin. The present application can not only break through the constraint of the lack of drilling data in the low exploration basin, but also has important significance for the evaluation of the geothermal resources and oil and gas potential of the system.

[0125] In some embodiments, step S1 further includes:

[0126] Collect regional tectonic evolution information, and extract the regional tectonic evolution event. The regional tectonic evolution event includes a regional tectonic compression event and / or a regional tectonic stretching event.

[0127] In one example, the regional tectonic compression events include plate collision, orogeny, rock deformation, tectonic folding, and seismic activity, etc. These regional tectonic compression events are common geological processes in the earth's surface and crust, reflecting the complexity and diversity of the evolution of the earth's tectonics. Studying these regional tectonic compression events helps to understand the internal structure of the region, tectonic movement, and the distribution and formation of geological resources.

[0128] In another example, the regional tectonic extension events include plate extension, volcanic eruption, fault zone formation, mantle column uplift, etc. These regional tectonic extension events are common geological processes in the earth's surface and crust, reflecting the complexity and diversity of the evolution of the earth's tectonics. Studying these regional tectonic extension events helps to understand the internal structure of the region, tectonic movement, and the distribution and formation of geological resources.

[0129] In the embodiments of the present application, taking the Po River Basin as an example, collecting the tectonic evolution information of the Po River Basin can obtain two groups of basement fault evolution events in the southeast and north-west directions of the Po River Basin, which can provide a channel for the upwelling of mantle material and heat flow.

[0130] In some embodiments, the tectonic fracture zone is analyzed according to the gravity gradient zone and the regional tectonic evolution event, comprising:

[0131] Obtaining the gravity gradient characteristic parameters of the gravity gradient zone. The gravity gradient characteristic parameters include at least one of the gravity gradient anomaly value, the gravity gradient vector, the gravity gradient modulus value, and the gravity gradient vector rotation angle.

[0132] According to the gravity gradient characteristic parameters, the tectonic fracture zone is analyzed from the regional tectonic evolution event.

[0133] In the embodiments of the present application, taking the Po River Basin as an example, according to the gravity average data, it can be identified that there are two gravity gradient zones in the northeast and northwest directions of the basin, and the gravity gradient zone in the northwest direction is faulted and broken by the gravity gradient zone in the northeast direction. Combined with the geological background of the Po River Basin, it can be determined that there are two tectonic fracture zones in the northeast and northwest directions of the Po River Basin.

[0134] The gravity gradient zone includes a vertical gradient zone and a horizontal gradient zone. For example, in one example, the vertical gradient zone of the gravity field in three-dimensional space can be calculated by the following formula:

[0135]

[0136] Where g(z) is a function of the gravity average data at height z, and h is the height variation.

[0137] In some embodiments, the gravity average data is calculated according to the gravity data, comprising:

[0138] The gravity average data is calculated according to the gravity data based on a sliding average method.

[0139] For example, in one example, the step of calculating the gravity average data according to the gravity data based on a sliding average method comprises:

[0140] Step 1, determining the size of the sliding average window for calculation;

[0141] Step 2, selecting a starting point as the starting position of the sliding average, and starting from the starting point, calculating the gravity average data of the center point in the window according to the size of the selected window and the gravity data in the window.

[0142] For example, in one example, for the original gravity data of M columns and N rows, with a line distance of Δx and a point distance of Δy, the size of the sliding average window is set to (m*Δx)*(n*Δy), where m is the number of lines in the x direction and n is the number of points in the y direction. The calculation formula of the gravity average data (i.e., the average value of the gravity data) in the window is:

[0143]

[0144] g ij is the gravity data of the ij point in the (m*Δx)*(n*Δy) window, is the gravity average value of the center point in the (m*Δx)*(n*Δy) window, m>0, n>0.

[0145] In some embodiments, the magnetic force residual data is calculated according to the magnetic force data, comprising:

[0146] The magnetic force average data is calculated according to the magnetic force data based on a sliding average method, and the magnetic force residual data is determined according to the difference between the magnetic force data and the magnetic force average data.

[0147] Further, as shown in Figure 4 , the magnetic force data can be original magnetic force data.

[0148] The magnetic force data can be obtained by downloading the EGM2008 original magnetic force data of the Powder River Basin with a grid spacing of 0.01'x0.01' from the International Global Geomagnetic Model Center (ICGEM) website.

[0149] For example, in one example, the step of calculating the magnetic force average data according to the magnetic force data based on a sliding average method comprises:

[0150] Step 1, determining the size of the sliding average window for calculation;

[0151] Step 2, selecting a starting point as the starting position of the sliding average, and starting from the starting point, sequentially calculating the magnetic force average data of the center point in the window according to the size of the selected window and the magnetic force data in the window.

[0152] For example, in one example, for the original magnetic force data of M columns and N rows, line spacing Δx, and point spacing Δy, the window size of the sliding average is set to (m*Δx)*(n*Δy), m is the number of lines in the x direction, and n is the number of points in the y direction. The calculation formula of the magnetic force average data in the window is:

[0153]

[0154] T ij is the magnetic force data of the ij point in the (m*Δx)*(n*Δy) window, is the magnetic force average data of the center point in the (m*Δx)*(n*Δy) window, m>0, n>0.

[0155] In another example, the magnetic force residual data is obtained according to the magnetic force data and the magnetic force average data, including the step of: subtracting the magnetic force average data from the magnetic force data to obtain the magnetic force residual data.

[0156] For example, further, in some embodiments, the magnetic force data, the magnetic force average data, and the magnetic force residual data are interpreted for geological significance.

[0157] wherein δΔT kl is the magnetic force residual field, ΔT kl is the magnetic force data field, is the magnetic force average field. The formula for obtaining the magnetic force average field is:

[0158]

[0159] In some embodiments, the magnetic force anomaly area is determined using the magnetic force residual data and the tectonic fracture zone, including: establishing the correlation between the magnetic force residual data and the tectonic fracture zone according to the geological interpretation of the magnetic force residual data; and determining the magnetic force anomaly area according to the correlation between the magnetic force residual data and the tectonic fracture zone.

[0160] In some embodiments, the distribution range of the magmatic intrusion area is determined using the magnetic force anomaly area and the regional tectonic evolution event, and the high geothermal anomaly area is predicted according to the distribution range of the magmatic intrusion body and the geothermal limit data, including:

[0161] The magnetic force anomaly characteristic parameters of the magnetic force anomaly area are obtained; wherein the magnetic force anomaly characteristic parameters include at least one of the magnetic field strength anomaly value, the magnetic deflection angle, the magnetic gradient vector rotation angle, and the magnetic polarity anomaly value.

[0162] According to the magnetic anomaly characteristic parameters, the magmatic intrusion area is determined from the regional tectonic evolution events, and the high ground temperature anomaly area is predicted by using the distribution range of the magmatic intrusion body.

[0163] The present application provides a method for predicting a high ground temperature anomaly area based on gravity and magnetic data, which is based on the analysis of regional tectonic evolution events and basin structure, etc. geological background, using gravity field characteristics (gravity data, gravity average data, gravity gradient zone) to analyze the tectonic fracture zone of the basin structure, using magnetic field characteristics (magnetic force data, magnetic residual data) to analyze the magnetic anomaly area and the distribution range of the magmatic intrusion body, and then using the limited geothermal data as a constraint to predict the high geothermal anomaly area of the basin, and realize the analysis of the geothermal field characteristics of the whole basin. The present application can not only break through the constraint of the lack of drilling data in the low exploration basin, but also has important significance for the evaluation of the geothermal resources and oil and gas potential of the system.

[0164] Example three

[0165] Figure 6 The structural block diagram of a device for predicting a high ground temperature anomaly area based on gravity and magnetic data provided by the present application is shown in Figure 6 The device for predicting a high ground temperature anomaly area based on gravity and magnetic data provided by the present application includes a tectonic fracture zone analysis module, a magnetic anomaly area determination module, and a high ground temperature anomaly area prediction module.

[0166] The tectonic fracture zone analysis module is configured to perform geological interpretation on gravity data and interpret the geological significance of gravity data anomalies; calculate gravity average data from the gravity data, determine gravity gradient zones using the gravity average data, and analyze tectonic fracture zones according to the gravity gradient zones and regional tectonic evolution information.

[0167] Further, the gravity data can be original free air gravity data.

[0168] The acquisition method of the original free air gravity data includes downloading the EGM2008 original free air gravity data of the Wuchuan Basin with a grid spacing of 0.01' x 0.01' from the International Global Earth Model Center (ICGEM) website.

[0169] It should be noted that the geological interpretation of gravity data is a common geophysical exploration method. For example, the geological interpretation of gravity data can help understand the density of underground rocks, geological structure and geological characteristics. Among them, gravity data anomalies usually reflect the density changes of underground rocks, the geological structure of underground rocks and the different anomalies of the geological characteristics of underground rocks in the gravity field.

[0170] In one example, the density of underground rocks varies, causing changes in the gravity field. For example, high-density rocks (such as granite) produce positive gravity data anomalies, while low-density rocks (such as volcanic rocks) cause negative gravity data anomalies. By analyzing the gravity data anomalies, the density variations of underground rocks can be inferred, and thus the geological structure and rock type can be inferred.

[0171] In another example, the geological structure of underground rocks causes changes in the gravity field. For example, geological structures such as fault zones, folds, and magma intrusions appear as gravity data anomalies in gravity data. By analyzing the gravity anomaly data, the spatial distribution and morphology of the gravity anomaly can be obtained, thereby identifying the location, morphology, and nature of the geological structure of the underground rocks.

[0172] In an embodiment of the present application, taking the Powder River Basin as an example, the original free air gravity data of the Powder River Basin is geologically interpreted. The middle of the Powder River Basin is a low gravity anomaly, indicating that the basin basement is relatively deep. The southern part of the Powder River Basin shows a weak low gravity anomaly, indicating that the basin basement is relatively shallow compared to the western part of the basin. The eastern part of the Powder River Basin is a high gravity anomaly, indicating that the basin basement is relatively shallow.

[0173] In one example, gravity average data is calculated from gravity data, and gravity gradient zones are determined using gravity average data. The gravity average data is obtained by filtering or smoothing the gravity data, which can reduce noise and sudden changes in the gravity data. Therefore, the gravity average data is smoother than the gravity data, and it is easier to determine the gravity gradient zone, thereby helping to reduce random errors in the data.

[0174] In addition, the tectonic fracture zone can be analyzed by analyzing the spatial distribution and morphology of the gravity gradient zone and the regional tectonic evolution event. By combining the geological evolution history of the Powder River Basin, the location, nature, and evolution process of the tectonic fracture zone can be inferred.

[0175] The tectonic fracture zone is usually a manifestation of the geological structure activity zone, which can be a fault zone, shear zone, or other geological structure features, and has important indicative significance for the geological evolution and structural features of the Powder River Basin. By comprehensively analyzing the gravity data and tectonic evolution events, the formation mechanism and geological significance of the tectonic fracture zone can be better understood.

[0176] In an embodiment of the present application, the magnetic anomaly area determination module is configured to geologically interpret the magnetic data and interpret the geological significance of the magnetic data anomalies; calculate the magnetic residual data from the magnetic data, and determine the magnetic anomaly area using the magnetic residual data and the tectonic fracture zone.

[0177] For example, in an embodiment of the present application, according to the magnetic residual data and the construction of the fracture zone, it can be determined that there is a high magnetic anomaly area extending from the north of the basin to the northwest in the Powder River Basin, and there is a high magnetic anomaly area extending from the middle of the basin to the northeast in the middle of the Powder River Basin.

[0178] It should be noted that magnetic data is an important non-destructive detection means in geological interpretation. After the geological interpretation of the magnetic data, information about the properties of the underground rocks, geological structure, geological characteristics, etc. can be obtained by analyzing the magnetic data. For example, in an example, the magnetic anomaly data in the magnetic data can reflect the magnetic difference of the underground rocks, such as rock type, mineral content, etc. By analyzing the size, shape, direction, etc. of the anomaly, it can be predicted that there may be a magnetic anomaly area, such as a fracture zone, a fold, a magmatic intrusion, etc.

[0179] In an embodiment of the present application, a magmatic intrusion is taken as an example for illustration. When the magma intrudes underground, it will change the magnetic characteristics of the surrounding rocks. This magnetic change will cause disturbance of the underground magnetic field, thereby forming a magnetic anomaly area.

[0180] In an embodiment of the present application, the high geothermal anomaly area prediction module is configured to determine a magmatic intrusion area according to the magnetic anomaly area and the regional tectonic evolution information, and to predict a high geothermal anomaly area using the distribution range of the magmatic intrusion.

[0181] The present application proposes a method for predicting a high geothermal anomaly area from gravity and magnetic data. Based on the analysis of regional tectonic evolution events and basin structure, etc. geological background, the method uses gravity field characteristics (gravity data, gravity average data, gravity gradient zone) to analyze the tectonic fracture zone of the basin structure, uses magnetic field characteristics (magnetic data, magnetic residual data) to analyze the magnetic anomaly area and the distribution range of the magmatic intrusion, and then uses the limited geothermal data as a constraint to predict the high geothermal anomaly area of the basin, and realizes the analysis of the geothermal field characteristics of the whole basin. The present application can not only break through the constraint of the lack of drilling data in the low exploration basin, but also has important significance for the evaluation of the geothermal resources and oil and gas potential of the system.

[0182] Embodiment four

[0183] Figure 7 A structural block diagram of an electronic device provided in an embodiment of the present application is shown in FIG. 1. Figure 7 As shown in FIG. 1, the electronic device 100 includes a memory 101, at least one processor 102, a computer program 103 stored in the memory 101 and executable on the at least one processor 102, and at least one communication bus 104.

[0184] The memory 101 can be configured to store a computer program 103, and the processor 102 can implement the method for predicting high-temperature anomaly areas by gravity and magnetic data in Embodiment I by running or executing the computer program stored in the memory 101 and calling data stored in the memory 101.

[0185] The memory 101 can mainly include a program storage area and a data storage area. The program storage area can store an operating system, at least one application program (such as a sound playing function, an image playing function, etc.) required by a function, etc. The data storage area can store data (such as audio data) created according to the use of the electronic device 100, etc. In addition, the memory 101 can include a non-volatile memory, such as a hard disk, a memory, a plug-in hard disk, a smart media card (SMC), a secure digital (SD) card, a flash card, at least one disk storage device, a flash memory device, or other non-volatile solid-state storage device.

[0186] The at least one processor 102 can be a central processing unit (CPU), and can also be other general-purpose processors, digital signal processors (DSPs), application specific integrated circuits (ASICs), field-programmable gate arrays (FPGAs) or other programmable logic devices, discrete gates or transistor logic devices, discrete hardware components, etc. The processor 102 can be a microprocessor or the processor 102 can also be any conventional processor, etc. The processor 102 is a control center of the electronic device 100, and connects all parts of the electronic device 100 through various interfaces and lines.

[0187] In one example, the memory 101 in the electronic device 100 stores a plurality of instructions to implement a seismic wave matching tracking method, and the processor 102 can execute the plurality of instructions to implement:

[0188] S1, geologically interpreting gravity data and interpreting the geological significance of gravity data anomalies; calculating gravity average data according to the gravity data, determining a gravity gradient zone using the gravity average data, and analyzing a tectonic fracture zone according to the gravity gradient zone and regional tectonic evolution events.

[0189] Further, as shown in Figure 2 The gravity data can be original free-air gravity data.

[0190] The original free-air gravity data is obtained by downloading the EGM2008 original free-air gravity data of the Powder River Basin with a grid spacing of 0.01' x 0.01' from the International Center for Global Earth Models (ICGEM) website.

[0191] It should be noted that the geological interpretation of gravity data is a common geophysical exploration method. For example, the geological interpretation of gravity data can help understand the density of underground rocks, geological structure and geological characteristics. Among them, the gravity data anomaly usually reflects the density change of the underground rock, the geological structure of the underground rock and the different anomalies of the geological characteristics of the underground rock in the gravity field.

[0192] In one example, the density of underground rocks is different, which will cause the change of the gravity field. For example, high-density rocks (such as granite) will produce positive gravity data anomalies, while low-density rocks (such as volcanic rocks) will cause negative gravity data anomalies. By analyzing the gravity data anomaly, the density change of the underground rock can be inferred, and then the geological structure and rock type can be inferred.

[0193] In another example, the geological structure of the underground rock will cause the change of the gravity field. For example, the geological structure such as fault zone, fold, magmatic intrusion, etc. shows gravity data anomaly in the gravity data. Through the analysis of the gravity anomaly data, the spatial distribution and shape of the gravity anomaly can be obtained, so as to identify the location, shape and nature of the geological structure of the underground rock.

[0194] In the embodiment of the application, the original free-air gravity data of the Powder River Basin is taken as an example for geological interpretation. The middle part of the Powder River Basin is a low gravity anomaly, indicating that the basin basement is deep. The southern part of the Powder River Basin shows a weak low gravity anomaly, indicating that the basin basement is relatively shallow compared with the western part of the basin. The eastern part of the Powder River Basin is a high gravity anomaly, indicating that the basin basement is shallow.

[0195] In one example, the gravity average data is calculated according to the gravity data, and the gravity gradient zone is determined by using the gravity average data. The gravity average data is obtained by filtering or smoothing the gravity data, which can reduce the noise and sudden change in the gravity data. Therefore, the gravity average data is smoother than the gravity data, and it is easier to determine the gravity gradient zone, thereby helping to reduce random errors in the data.

[0196] In addition, the tectonic fracture zone is analyzed by combining the gravity gradient zone and the regional tectonic evolution event. By analyzing the spatial distribution and shape of the gravity gradient zone, and combining the geological tectonic evolution history of the Powder River Basin, the location, nature and evolution process of the tectonic fracture zone can be inferred.

[0197] Tectonic fracture zone is usually a manifestation of tectonic activity zone, which can be a fault zone, shear zone and other geological structure features, and has important indicative significance for the geological evolution and structural characteristics of the Powell basin. By comprehensively analyzing gravity data and tectonic evolution events, the formation mechanism and geological significance of the tectonic fracture zone can be better understood.

[0198] S2, geologically interpreting the magnetic data and explaining the geological significance of the magnetic data anomaly; calculating the magnetic residual data according to the magnetic data, and determining the magnetic anomaly area by using the magnetic residual data and the tectonic fracture zone.

[0199] It should be noted that magnetic data is an important non-destructive detection means in geological interpretation. After the geological interpretation of the magnetic data, information about the properties of underground rocks, geological structures, geological features, etc. can be obtained by analyzing the magnetic data. For example, in an example, the magnetic anomaly data in the magnetic data can reflect the magnetic differences of the underground rocks, such as rock types, mineral content, etc. By analyzing the size, shape, direction and other characteristics of the anomaly, it can be predicted that there may be a magnetic anomaly area, such as a fault zone, a fold, a magmatic intrusion, etc.

[0200] In the embodiments of the present application, the magmatic intrusion is taken as an example for illustration. When the magma intrudes underground, it will change the magnetic characteristics of the surrounding rocks. This magnetic change will cause disturbance of the underground magnetic field, thereby forming a magnetic anomaly area.

[0201] S3, determining the distribution range of the magmatic intrusion area by using the magnetic anomaly area and the regional tectonic evolution event, and predicting the high geothermal anomaly area according to the distribution range of the magmatic intrusion and the geothermal limit data.

[0202] The present application proposes a method for predicting a high geothermal anomaly area based on the analysis of regional tectonic evolution events and basin structure, etc. geological background, using gravity field characteristics (gravity data, gravity average data, gravity gradient zone) to analyze the tectonic fracture zone of the basin structure, using magnetic field characteristics (magnetic data, magnetic residual data) to analyze the magnetic anomaly area and the distribution range of the magmatic intrusion, and then using limited geothermal data as a constraint to predict the high geothermal anomaly area of the basin, and to realize the analysis of the geothermal field characteristics of the whole basin. The present application can not only break through the constraint of the lack of drilling data in the low exploration basin, but also has important significance for the evaluation of the geothermal resources and oil and gas potential of the system.

[0203] Example five

[0204] The modules / units integrated in the electronic device 100, if implemented in the form of software function units and sold or used as independent products, can be stored in a computer-readable storage medium. Based on such understanding, all or part of the processes in the above-mentioned embodiment methods can also be completed by a computer program instructing related hardware, and the computer program can be stored in a computer-readable storage medium. When the computer program is executed by a processor, the steps of the above-mentioned method embodiment can be implemented.

[0205] The computer program includes computer program code, which can be in the form of source code, object code, executable files, or some intermediate forms. The computer-readable medium can include any entity or device capable of carrying computer program code, recording medium, U disk, mobile hard disk, magnetic disk, optical disk, computer memory, and read-only memory (ROM).

[0206] Those skilled in the art should understand that the embodiments of the present application can be provided as a method, a system, or a computer program product. Therefore, the present application can take the form of a complete hardware embodiment, a complete software embodiment, or an embodiment combining software and hardware aspects. Moreover, the present application can take the form of a computer program product implemented on one or more computer-readable storage media (including but not limited to disk storage, CD-ROM, optical storage, etc.) containing computer-usable program code.

[0207] The present application is described with reference to flowcharts and / or block diagrams of the method, device (system), and computer program product according to the embodiments of the present application. It should be understood that each flow and / or block in the flowcharts and / or block diagrams, and the combination of flows and / or blocks in the flowcharts and / or block diagrams can be implemented by computer program instructions. These computer program instructions can be provided to a general-purpose computer, a special-purpose computer, an embedded processor, or other programmable data processing devices to produce a machine, so that the instructions executed by the computer or other programmable data processing devices produce a device that implements the functions specified in the flowcharts and / or block diagrams. Figure 1 The functions specified in one or more flows and / or blocks Figure 1 The device that implements the functions specified in one or more blocks.

[0208] These computer program instructions can also be stored in a computer-readable storage medium that can guide the computer or other programmable data processing devices to work in a specific way, so that the instructions stored in the computer-readable storage medium produce a product including instruction devices that implement the functions specified in the flowcharts and / or block diagrams. Figure 1 The functions specified in one or more flows and / or blocks Figure 1 The device that implements the functions specified in one or more blocks.

[0209] These computer program instructions can also be loaded into a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart block or blocks. Figure 1 one or more flowcharts and / or blocks Figure 1 one or more flowcharts and / or blocks

[0210] In the description of the present specification, the description of the terms "one embodiment", "an example", "a specific example" and the like means that the specific features, structures, materials or characteristics described in connection with the embodiment or example are included in at least one embodiment or example of the present application. In the present specification, the illustrative expressions of the above terms do not necessarily refer to the same embodiment or example. Moreover, the specific features, structures, materials or characteristics described can be combined in any one or more embodiments or examples in a suitable manner.

[0211] Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of the present application but not to limit it, and although the present application has been described in detail with reference to the above embodiments, those skilled in the art should understand that the specific embodiments of the present application can be modified or replaced by the equivalent, without departing from the spirit and scope of the present application, and any modification or equivalent replacement should be covered in the protection scope of the claims of the present application.

Claims

1. A method for predicting high geothermal anomaly areas using gravity and magnetic data, characterized in that, The method includes the following steps: S1. Perform geological interpretation on gravity data and explain the geological significance of the gravity data anomalies; calculate gravity average data based on the gravity data, determine gravity gradient zones using the gravity average data, and analyze tectonic fracture zones based on the gravity gradient zones and regional tectonic evolution events. S2. Perform geological interpretation on the magnetic data and explain the geological significance of the magnetic data anomalies; calculate the residual magnetic data based on the magnetic data, and use the residual magnetic data and the tectonic fracture zone to determine the magnetic anomaly area; S3. Determine the distribution range of the magma intrusion zone using the magnetic anomaly zone and the regional tectonic evolution events, and predict the high geothermal anomaly zone based on the distribution range of the magma intrusion body and the geothermal limit data.

2. The method according to claim 1, characterized in that, Step S1 is preceded by: Collect regional tectonic evolution information and extract the regional tectonic evolution events; The regional tectonic evolution events include regional tectonic compression events and / or regional tectonic stretching events.

3. The method according to claim 2, characterized in that, Based on the gravity gradient zone and the analytical construction of regional tectonic evolution events, a fracture zone is constructed, including: Obtain the gravity gradient characteristic parameters of the gravity gradient zone; wherein, the gravity gradient characteristic parameters include at least one of gravity gradient outliers, gravity gradient vector, gravity gradient magnitude, and gravity gradient vector rotation angle; The tectonic fracture zone is analyzed from the regional tectonic evolution events based on the gravity gradient characteristic parameters.

4. The method according to claim 1, characterized in that, The gravity average data is calculated based on the gravity data, including: The average gravity data is calculated based on the moving average method and the gravity data.

5. The method according to claim 1, characterized in that, The remaining magnetic force data is calculated based on the magnetic force data, including: The average magnetic force data is calculated based on the moving average method and the magnetic force data. The remaining magnetic force data is determined based on the difference between the magnetic force data and the average magnetic force data.

6. The method according to claim 5, characterized in that, Determining magnetic anomaly zones using the residual magnetic data and the structural fracture zone includes: Based on the geological interpretation of the magnetic residual data, establish the correlation between the magnetic residual data and the tectonic fracture zone; The magnetic anomaly zone is determined based on the correlation between the remaining magnetic force data and the structural fracture zone.

7. The method according to claim 6, characterized in that, The distribution range of magmatic intrusion zones is determined using the magnetic anomaly zones and the regional tectonic evolution events. High geothermal anomaly zones are predicted based on the distribution range of the magmatic intrusions and geothermal constraint data, including: Obtain the magnetic anomaly characteristic parameters of the magnetic anomaly region; wherein, the magnetic anomaly characteristic parameters include at least one of the following parameters: magnetic field strength anomaly value, magnetic deflection angle, magnetic gradient vector rotation angle, and magnetic polarity anomaly value; The magma intrusion zone and its distribution range are determined from the regional tectonic evolution events based on the magnetic anomaly characteristic parameters. High geothermal anomaly zones are predicted based on the distribution range of the magmatic intrusion and the geothermal constraints.

8. A device for predicting high geothermal anomaly zones based on gravity and magnetic data, characterized in that, include: A fracture zone analysis module is constructed and configured to perform geological interpretation of gravity data and interpret the geological significance of the gravity data anomalies; The gravity average data is calculated based on the gravity data, the gravity gradient zone is determined using the gravity average data, and the fracture zone is analyzed and constructed based on the gravity gradient zone and regional tectonic evolution information. The magnetic anomaly zone identification module is configured to perform geological interpretation of magnetic data and interpret the geological significance of the magnetic data anomalies; Calculate the remaining magnetic data based on the magnetic data, and use the remaining magnetic data and the structural fracture zone to determine the magnetic anomaly zone; The high geothermal anomaly prediction module is configured to determine the distribution range of magma intrusion zones using the magnetic anomaly zones and the regional tectonic evolution events, and to predict high geothermal anomaly zones based on the distribution range of the magma intrusion bodies and geothermal constraint data.

9. An electronic device comprising a memory and a processor, wherein the memory stores a computer program, Its features are, When the processor executes the computer program, it implements the steps of the method as described in any one of claims 1 to 7.

10. A computer-readable storage medium, characterized in that, The computer-readable storage medium stores at least one instruction that, when executed by a processor, implements the steps of the method as claimed in any one of claims 1 to 7.

Images