A method, device and equipment for exploring hot dry rock target area based on aeromagnetic data

Abstract

The present invention relates to a method, a device and equipment for exploring a hot dry rock target area based on aeromagnetic data. The method comprises: acquiring aeromagnetic data of an exploration area; determining at least rock mass distribution information, Curie point depth fluctuation information, crust temperature distribution information and crust fracture distribution information based on the aeromagnetic data; determining the hot dry rock target area information of the exploration area by using the determined rock mass distribution information, Curie point depth fluctuation information, crust temperature distribution information and crust fracture distribution information; wherein the hot dry rock target area information comprises: possible distribution areas of the hot dry rock target area and the probability corresponding to each distribution area.

Description

Technical Field

[0001] The present invention relates to the field of geological exploration technology, and in particular to a method, a device, equipment and a storage medium for exploring a hot dry rock target area based on aeromagnetic data. Background Art

[0002] Hot dry rock refers to rock (body) that is high in temperature underground but lacks fluid due to low porosity and permeability. Hot dry rock is a clean energy source, and its development and utilization is of great significance for alleviating energy shortages and reducing environmental pollution.

[0003] Hot dry rock target areas refer to specific areas with the potential for developing hot dry rock resources. These areas usually have specific geological structures and geothermal geological conditions, and hot dry rocks may exist in these areas. Therefore, determining hot dry rock target areas is one of the key steps to achieve sustainable development and utilization of hot dry rock resources. Summary of the invention

[0004] In order to solve the above-mentioned problems in the prior art, the present invention provides a method, device, equipment and storage medium for surveying hot dry rock target areas based on aeromagnetic data, so as to efficiently and conveniently determine the hot dry rock target areas.

[0005] In a first aspect, an embodiment of the present application provides a method for surveying a hot dry rock target area based on aeromagnetic data, comprising: acquiring aeromagnetic data of a survey area;

[0006] Based on the above aeromagnetic data, at least rock mass distribution information, Curie point depth fluctuation information, crustal temperature distribution information and crustal fracture distribution information are determined;

[0007] Determine the hot dry rock target area information of the exploration area by using the determined rock mass distribution information, Curie point depth fluctuation information, crust temperature distribution information and crust fracture distribution information;

[0008] The above hot dry rock target area information includes:

[0009] Possible distribution areas of hot dry rock targets and the corresponding probabilities of each distribution area.

[0010] In an optional implementation of the first aspect, the determination of the hot dry rock target area information of the exploration area by using the determined rock mass distribution information, the Curie point depth fluctuation information, the crust temperature distribution information and the crust fracture distribution information includes:

[0011] Determine a first area corresponding to which the rock mass distribution information meets a first preset condition;

[0012] Determine a second area corresponding to which the Curie surface depth fluctuation information meets a second preset condition;

[0013] Determine a third region corresponding to which the crust temperature distribution information meets a third preset condition;

[0014] Determine a fourth region corresponding to which the crustal fracture distribution information meets a fourth preset condition;

[0015] Based on the first area, the second area, the third area and the fourth area, possible distribution areas of the hot dry rock target area are determined.

[0016] In an optional implementation of the first aspect, the first area corresponds to a first weight value, the second area corresponds to a second weight value, and the third area and the fourth area correspond to a third weight value;

[0017] The first weight value is greater than the second weight value, and the second weight value is greater than the third weight value;

[0018] And, based on the first area, the second area, the third area and the fourth area, the possible distribution area of ​​the hot dry rock target area is determined, including:

[0019] In response to determining that the weight value of a certain area is greater than a preset weight value, the area is determined to be a distribution area, and the probability corresponding to the distribution area is determined according to the size of the weight value.

[0020] In an optional implementation of the first aspect, the rock mass distribution information is determined in the following manner:

[0021] Based on the above aeromagnetic data, determine the three-dimensional distribution information of crustal magnetization intensity;

[0022] The above rock mass distribution information is determined based on the three-dimensional distribution information of the crustal magnetization intensity.

[0023] In an optional implementation manner of the first aspect, the Curie surface depth fluctuation information is determined in the following manner:

[0024] Based on the above aeromagnetic data, the temperature interface that corresponds to the demagnetization temperature of ilmenite is determined;

[0025] Based on the determined temperature interface, the Curie surface depth fluctuation information is determined.

[0026] In an optional implementation of the first aspect, the crust temperature distribution information is determined in the following manner:

[0027] Establish a one-dimensional steady-state heat conduction equation;

[0028] Based on the above aerial survey data and the above one-dimensional steady-state heat conduction equation, the crust temperature distribution information is determined.

[0029] In an optional implementation of the first aspect, the crust fracture distribution information is determined in the following manner:

[0030] The above aeromagnetic data are processed to enhance the THVH value;

[0031] The above crustal fracture distribution information is determined based on the aeromagnetic data after enhancement processing.

[0032] In a second aspect, an embodiment of the present application provides a device for surveying a hot dry rock target area based on aeromagnetic data, comprising: an acquisition unit, for acquiring aeromagnetic data of a survey area;

[0033] A first determination unit is used to determine at least rock mass distribution information, Curie point depth fluctuation information, crust temperature distribution information and crust fracture distribution information based on the aeromagnetic data;

[0034] A second determination unit is used to determine the hot dry rock target area information of the exploration area by using the determined rock mass distribution information, the Curie point depth fluctuation information, the crust temperature distribution information and the crust fracture distribution information;

[0035] The above hot dry rock target area information includes:

[0036] Possible distribution areas of hot dry rock targets and the corresponding probabilities of each distribution area.

[0037] In a third aspect, an embodiment of the present application provides an electronic device, comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor implements the method for surveying hot dry rock target areas based on aeromagnetic data as described in any aspect of the first invention when executing the computer program.

[0038] In a fourth aspect, an embodiment of the present application provides a computer-readable medium having a computer program stored thereon, which, when executed by a processor, implements the method for exploring a hot dry rock target area based on aeromagnetic data as described in any aspect of the first invention.

[0039] The beneficial effects of the present invention include: by acquiring aeromagnetic data of the survey area, the rock mass distribution information, Curie point depth fluctuation information, crust temperature distribution information and crust fracture distribution information of the survey area can be determined by using the aeromagnetic data, and the hot dry rock target area information of the survey area can be determined by using the rock mass distribution information, Curie point depth fluctuation information, crust temperature distribution information and crust fracture distribution information at the same time. In this way, the hot dry rock target area information of the survey area can be determined more efficiently and accurately, and a large amount of manpower and material costs can be saved.

[0040] Moreover, the determined hot dry rock target area information also includes the possible distribution areas of the hot dry rock target areas and the corresponding probabilities of each distribution area. This can also facilitate users to conduct more targeted on-site verification of each distribution area in the future. BRIEF DESCRIPTION OF THE DRAWINGS

[0041] Figure 1 A flowchart of the steps of a method for surveying hot dry rock target areas based on aeromagnetic data provided by an embodiment of the present invention;

[0042] Figure 2 A schematic diagram of the distribution of a hot dry rock target area provided by an embodiment of the present invention;

[0043] Figure 3 A module block diagram of a device for surveying hot dry rock target areas based on aeromagnetic data provided by an embodiment of the present invention;

[0044] Figure 4 A module block diagram of an electronic device provided by an embodiment of the present invention. DETAILED DESCRIPTION

[0045] In the following description, specific details such as specific system structures, technologies, etc. are provided for the purpose of illustration rather than limitation, so as to provide a thorough understanding of the embodiments of the present application. However, it should be clear to those skilled in the art that the present application may also be implemented in other embodiments without these specific details. In other cases, detailed descriptions of well-known systems, devices, circuits, and methods are omitted to prevent unnecessary details from obstructing the description of the present application.

[0046] It should be understood that when used in the present specification and the appended claims, the term "comprising" indicates the presence of described features, integers, steps, operations, elements and / or components, but does not exclude the presence or addition of one or more other features, integers, steps, operations, elements, components and / or combinations thereof.

[0047] It should also be understood that the term “and / or” used in the specification and appended claims refers to any and all possible combinations of one or more of the associated listed items, and includes these combinations.

[0048] As used in this specification and the appended claims, the term "if" may be interpreted as "when" or "upon" or "in response to determining" or "in response to detecting" depending on the context.

[0049] In addition, in the description of the present application specification and the appended claims, the terms "first", "second", "third", etc. are only used to distinguish the descriptions and cannot be understood as indicating or implying relative importance.

[0050] References to "one embodiment" or "some embodiments" etc. described in the specification of this application mean that one or more embodiments of the present application include specific features, structures or characteristics described in conjunction with the embodiment. Therefore, the statements "in one embodiment", "in some embodiments", "in some other embodiments", "in some other embodiments", etc. that appear in different places in this specification do not necessarily refer to the same embodiment, but mean "one or more but not all embodiments", unless otherwise specifically emphasized in other ways. The terms "including", "comprising", "having" and their variations all mean "including but not limited to", unless otherwise specifically emphasized in other ways.

[0051] At present, the steps for determining the hot dry rock target area mainly include the following aspects:

[0052] First, data collection. For example, you can collect existing geological, geothermal hot spring, coal field, mining exploration data in the work area, as well as geophysical exploration data such as earthquake, gravity, aeromagnetic, and magnetotelluric sounding. Perform a comprehensive analysis of these data, and determine the areas where hot dry rocks may exist based on the large geological background.

[0053] Second, geological survey and geophysical exploration. For example, field geological surveys can be conducted in areas where hot dry rocks may exist, including detailed investigations of geological features, structures, lithology, geophysics, and geochemistry.

[0054] Third, data processing and interpretation. For example, the collected geological and geophysical data are processed and interpreted to extract information related to hot dry rocks, such as the rock's magnetic, electrical, density and other physical properties.

[0055] Fourth, anomaly identification and delineation. For example, based on the results of data processing and interpretation, abnormal areas where hot dry rocks may exist can be identified. These abnormal areas usually manifest as geomagnetic anomalies, geothermal anomalies, gravity anomalies, etc. These abnormal areas are delineated to preliminarily determine the location and scope of the hot dry rock target area.

[0056] Fifth, drilling verification. For example, drilling verification can be carried out in the preliminarily determined hot dry rock target area. Through drilling sampling and temperature measurement, the temperature, lithology, heat reservoir type and other information of underground rocks can be directly obtained to further confirm the existence and distribution of hot dry rocks.

[0057] Sixth, comprehensive assessment and development decision-making. For example, comprehensive assessment of hot dry rock target areas is conducted based on drilling verification results and other geological, geophysical and other data.

[0058] However, the method of determining the hot dry rock target area in the related technology requires a lot of time due to the need for data collection in the early stage, and also requires field geological surveys, which makes the geological exploration efficiency low. In other words, the method of determining the hot dry rock target area in the related technology is inefficient, and in areas with complex terrain, it may be impossible to survey.

[0059] To be more specific, the currently commonly used exploration methods for hot dry rock target areas mainly include: magnetotelluric method, high-density electrical method, gravity method and direct drilling method. In the shallow crust, magnetotelluric method and high-density electrical method are mainly sensitive to fluids but not to temperature. Moreover, these two methods have a long construction period, sparse data sampling and high cost, and are not suitable for large-scale and large-scale hot dry rock surveys in the early stage. The gravity method is also insensitive to temperature, and its resolution is low, which is less effective for hot dry rock exploration. The drilling method can directly measure the material properties and temperature of underground rocks with high precision, but the drilling cost is very high and is not suitable for the early hot dry rock survey. In other words, the applicability of the relevant detection methods is not high.

[0060] The method for exploring hot dry rock target areas based on aeromagnetic data provided by the present disclosure can directly determine the aeromagnetic data, and then use the aeromagnetic data to determine the rock mass distribution information, Curie point depth fluctuation information, crust temperature distribution information and crust fracture distribution information, and use the determined rock mass distribution information, Curie point depth fluctuation information, crust temperature distribution information and crust fracture distribution information to determine the hot dry rock target area information of the exploration area. In this way, the hot dry rock target area information of the exploration area can be determined efficiently and at a low cost.

[0061] See also Figure 1 The embodiment of the present application provides a method for surveying a hot dry rock target area based on aeromagnetic data, including: steps 101 to 104.

[0062] Step 101: Acquire aeromagnetic data of the survey area.

[0063] As an example, aeromagnetic data is the earth's magnetic field data obtained using an aeromagnetic measurement system. For example, data can be obtained by measuring magnetometer equipment carried on an aircraft or an unmanned aerial vehicle, thereby obtaining aeromagnetic data of the survey area.

[0064] Step 102: Based on the aeromagnetic data, at least rock mass distribution information, Curie point depth fluctuation information, crust temperature distribution information and crust fracture distribution information are determined.

[0065] As an example, rock mass distribution information can indicate the distribution of different rock types, helping to assess their potential as mineral resources or exploitable geological materials.

[0066] As an example, the Curie depth fluctuation information may indicate different attributes of the stratum interfaces in the stratum, for example, geological interfaces and geophysical interfaces. These interfaces are formed due to factors such as the physical properties, chemical composition, structure and construction of underground rocks.

[0067] As an example, crust temperature distribution information refers to the temperature distribution inside the earth's crust.

[0068] As an example, crustal fracture distribution information refers to the distribution of different fractures in the crust.

[0069] It can be seen that in the present disclosure, at least the rock mass distribution information, Curie point depth fluctuation information, crustal temperature distribution information and crustal fracture distribution information can be determined through aeromagnetic data, and then the determined rock mass distribution information, Curie point depth fluctuation information, crustal temperature distribution information and crustal fracture distribution information can be further analyzed to determine the hot dry rock target area information.

[0070] This method does not require users to read a large amount of data in advance, and does not require users to conduct field surveys. It can improve the efficiency of determining rock distribution information, Curie point depth fluctuation information, crustal temperature distribution information and crustal fracture distribution information. In this way, the hot dry rock target area information can be determined more efficiently.

[0071] Step 103: Determine the hot dry rock target area information of the exploration area using the determined rock mass distribution information, Curie point depth fluctuation information, crust temperature distribution information and crust fracture distribution information.

[0072] Here, the above hot dry rock target area information includes:

[0073] Possible distribution areas of hot dry rock targets and the corresponding probabilities of each distribution area.

[0074] Generally speaking, the distribution area of ​​medium-strong magnetic bodies in the middle and upper crust, the distribution area of ​​shallow Curie point, the distribution area of ​​high temperature in the middle and upper crust, and the distribution area of ​​deep and large faults are favorable hot dry rock targets (it can also be understood that the area where hot dry rock targets are usually located is generally in the distribution area of ​​medium-strong magnetic bodies in the middle and upper crust, the distribution area of ​​shallow Curie point, the distribution area of ​​high temperature in the middle and upper crust, and the distribution area of ​​deep and large faults). Therefore, based on the rock mass distribution information, the Curie point depth fluctuation information, the crust temperature distribution information, and the crust fault distribution information, the possible distribution area of ​​the hot dry rock target area can be accurately determined.

[0075] Since the rock mass distribution information, Curie point depth fluctuation information, crust temperature distribution information and crust fracture distribution information have different influencing factors (or can be understood as different control factors) on the hot dry rock target area, the probability corresponding to each distribution area can also be determined based on the rock mass distribution information, Curie point depth fluctuation information, crust temperature distribution information and crust fracture distribution information. Through this numerical method, users can more intuitively understand the situation of the hot dry rock target area in the exploration area.

[0076] It can be seen that in the present disclosure, by acquiring aeromagnetic data of the survey area, the aeromagnetic data can be used to determine the rock mass distribution information, Curie point depth fluctuation information, crust temperature distribution information and crust fracture distribution information of the survey area, and the rock mass distribution information, Curie point depth fluctuation information, crust temperature distribution information and crust fracture distribution information can be used simultaneously to determine the hot dry rock target area information of the survey area. This method can more efficiently and accurately determine the hot dry rock target area information of the survey area, and also save a lot of manpower and material costs.

[0077] Moreover, the determined hot dry rock target area information also includes the possible distribution areas of the hot dry rock target areas and the corresponding probabilities of each distribution area. This can also facilitate users to conduct more targeted on-site verification of each distribution area in the future.

[0078] In some embodiments, in the above step 103, "using the determined rock mass distribution information, Curie point depth fluctuation information, crust temperature distribution information and crust fracture distribution information to determine the hot dry rock target area information of the exploration area" may specifically include:

[0079] Determine a first area corresponding to which the rock mass distribution information meets a first preset condition;

[0080] Determine a second area corresponding to which the Curie surface depth fluctuation information meets a second preset condition;

[0081] Determine a third region corresponding to which the crust temperature distribution information meets a third preset condition;

[0082] Determine a fourth region corresponding to which the crustal fracture distribution information meets a fourth preset condition;

[0083] Based on the above-mentioned first area, second area, third area and fourth area, the possible distribution areas of the above-mentioned hot dry rock target area are determined.

[0084] In some possible implementations, an overlapping area among the first area, the second area, the third area, and the fourth area may be used as a distribution area.

[0085] As an example, the first preset condition can be understood as the magnetic strength distribution information. It should be noted that the rock mass distribution information is related to the magnetic body distribution, and the first area can be understood as the medium-strong magnetic body distribution area in the middle and upper crust. That is, it can be understood that based on the first preset condition and the rock mass distribution information, the medium-strong magnetic body distribution area (first area) in the middle and upper crust can be determined.

[0086] The second preset condition can be understood as depth information. Of course, the specific value of the depth information can be limited according to the actual situation, and the second area can be understood as a shallow Curie surface distribution area, that is, based on the Curie surface depth fluctuation information and the second preset condition, the shallow Curie surface distribution area (second area) can be determined.

[0087] The third preset condition can be understood as temperature information. Of course, the specific value of the temperature information can also be limited according to actual conditions. The third region can be understood as the high temperature distribution region in the middle and upper crust, that is, based on the crust temperature distribution information and the third preset condition, the high temperature distribution region in the middle and upper crust (the third region) can be determined.

[0088] The fourth preset condition can be understood as the depth and fracture degree information. Of course, the specific values ​​of the depth and fracture degree can be limited according to the actual situation. The fourth area can be understood as the deep and large fracture distribution area, that is, based on the crust fracture distribution information and the fourth preset condition, the deep and large fracture distribution area (fourth area) can be determined.

[0089] Generally speaking, the distribution areas of medium-strong magnetic bodies in the middle and upper crust, the shallow Curie distribution areas, the high temperature distribution areas in the middle and upper crust, and the deep and large fault distribution areas are favorable hot dry rock targets. Therefore, by determining the first area, the second area, the third area and the fourth area, the possible distribution areas of the hot dry rock targets can be determined more efficiently and conveniently.

[0090] In some embodiments, the first area corresponds to a first weight value, the second area corresponds to a second weight value, and the third area and the fourth area correspond to a third weight value;

[0091] The first weight value is greater than the second weight value, and the second weight value is greater than the third weight value;

[0092] Furthermore, the determination of the possible distribution area of ​​the hot dry rock target area based on the first area, the second area, the third area and the fourth area may specifically include:

[0093] In response to determining that the weight value of a certain area is greater than a preset weight value, the area is determined to be a distribution area, and the probability corresponding to the distribution area is determined according to the size of the weight value.

[0094] As an example, specific values ​​of the first weight value, the second weight value, and the third weight value may be set according to actual needs.

[0095] Here, by setting the first weight value, the second weight value and the third weight value, the weights corresponding to different factors may be different, so that the distribution area may be determined more accurately.

[0096] For example, the high temperature distribution area in the middle and upper crust is the primary control factor, the distribution of medium-strong magnetic bodies is the secondary control factor, and the Curie point depth and deep large fault distribution area are the tertiary control factors. That is, the high temperature distribution area in the middle and upper crust is most likely to form a hot dry rock target area, while the distribution of medium-strong magnetic bodies is also more likely to form a hot dry rock target area, and the Curie point depth and deep large fault distribution area are also likely to form a hot dry rock target area.

[0097] It can also be understood as: the area where the hot dry rock target area is distributed is most likely the high-temperature distribution area in the middle and upper crust, and there is a certain probability that this area is the medium-strong magnetic body distribution area, and there is also a certain probability that it is the Curie point depth and deep and large fault distribution area.

[0098] That is, by setting different weight values ​​for different areas, the distribution area can be determined more accurately and the probability corresponding to the distribution area can be determined efficiently.

[0099] It should be understood that the specific numerical value of the preset weight value can be limited according to actual conditions, and the specific numerical value of the preset weight value is not limited here.

[0100] In order to facilitate understanding of the concepts in this disclosure, Figure 2 To explain, Figure 2 It can be understood as a probability diagram corresponding to the distribution area, combined with Figure 2 It can be seen that users can clearly know which areas are likely to be hot dry rock target areas, as well as the corresponding probability distribution, so that they can conduct targeted exploration of hot dry rock target areas.

[0101] In some embodiments, the rock mass distribution information may be determined by:

[0102] Determine the three-dimensional distribution information of crustal magnetization intensity based on aeromagnetic data;

[0103] The above rock mass distribution information is determined based on the three-dimensional distribution information of the crustal magnetization intensity.

[0104] Specifically, by performing aeromagnetic oblique magnetization three-dimensional magnetization intensity inversion, the location of the magnetic source body can be obtained, and thus the rock mass distribution information can be obtained.

[0105] Generally speaking, if the latitude span of the exploration area is large, its magnetic dip and declination may vary greatly, and processing methods such as polarization may introduce some human errors. At this time, in order to more accurately determine the rock distribution information, the magnetic dip and declination can be introduced into the inversion of the three-dimensional magnetic anomaly, and the deep magnetic source body can be three-dimensionally characterized. Then, the rock distribution information can be determined using the magnetic source body completed with the three-dimensional characterization.

[0106] As an example, the objective function of the three-dimensional focusing inversion of magnetic anomaly is as follows (Formula 1):

[0107] P α (m)=Φ(m)+αS(m)

[0108] =||W d (Am-d)|| 2 +α||W m W e (mm apr )|| 2

[0109] In formula 1, m is the model space; Φ(m) is the data fitting functional; S(m) is the model stabilization functional; α is the regularization factor, which can be used to balance the data fitting functional and the model stabilization functional; W d and W m are the weight matrices of data and model respectively; d is the observed data; W e To focus the functional matrix, the minimum support functional matrix is ​​used here; A is the three-dimensional magnetic anomaly forward operator of the oblique magnetization cube, and its expression can be as follows (Formula 2):

[0110]

[0111] In Formula 2, J is the modulus of magnetization, μ0 is the vacuum permeability, π is the circumference of a circle, (x, y, z) is the coordinate of the observation point, (ξ, η, ζ) is the coordinate of the cubic model grid, a, b and c are the lengths of the grid unit in the ξ, η and ζ directions respectively, and the grid center coordinates are (x0, y0, z0). The expressions of k1, k2, k3, k4, k5 and k6 in Formula 2 can be as follows (Formula 3):

[0112]

[0113] Where I is the magnetic inclination and D is the magnetic declination. In the case of unknown remanent magnetism, the direction of the geomagnetic field is set to be consistent with the direction of the magnetization intensity. By using Formula 2 and Formula 3, the magnetic field value of any cube grid model at any observation point can be forward modeled, and the forward modeling process will be called as a subroutine in each iterative inversion process.

[0114] It should be understood that in the actual inversion process, since it is difficult to evaluate the error of the observed data, the data weighting function W d Set to a constant diagonal matrix. W m As a model weighting matrix, its purpose is to overcome the skin effect of the inversion process, that is, to prevent the inversion results from being close to the surface. Since the comprehensive sensitivity matrix has a good weighting effect on the model of the entire space, this matrix can be selected as the model weighting method, and its expression can be as follows (Formula 4):

[0115]

[0116] Here, A is the cubic magnetic field forward operator under uniform magnetization conditions, that is, Formula 2.

[0117] W e As a focusing functional matrix, its purpose is to make the non-zero space of the inversion result smaller, so as to obtain a magnetic source distribution with clear boundaries. Here, the minimum support functional can be used to achieve this goal. The expression of this functional can be (Formula 5):

[0118]

[0119] Among them, β is a small value, and its main functions include: preventing the denominator from being zero, the division formula from being meaningless, and focusing the magnetization intensity model.

[0120] It should be understood that by combining Formulas 2-5 and using a preconditioned conjugate gradient algorithm, the objective function (Formula 1) can be solved and a high-precision three-dimensional distribution of crustal magnetization intensity can be obtained.

[0121] Of course, the specific solution process can be carried out by computer. For example, a C language algorithm solution program can be written according to the actual operation situation to solve the objective function (Formula 1), thereby obtaining the three-dimensional distribution of crustal magnetization intensity and determining the rock distribution information in combination with the three-dimensional distribution of crustal magnetization intensity.

[0122] In some embodiments, the Curie surface depth fluctuation information may be determined as follows:

[0123] Based on aeromagnetic data, the temperature interface corresponding to the demagnetization temperature of ilmenite was determined;

[0124] Based on the determined temperature interface, the Curie surface depth fluctuation information is determined.

[0125] As an example, the Curie point is the temperature interface at which ferromagnetic minerals in the lithosphere demagnetize, and below this interface, as the temperature exceeds the Curie point of the ferromagnetic minerals in the rock, these ferromagnetic minerals transform into paramagnetic minerals.

[0126] As an example, the demagnetization temperature of ilmenite is generally considered to be about 580°C. This temperature interface can be obtained by the central point method of aeromagnetics. For a magnetic layer that extends infinitely in the horizontal direction, the radial average amplitude spectrum of the total magnetic field generated by it can be expressed as (Formula 6):

[0127]

[0128] In Formula 6, A ΔT (|k|) and A M (|k|) are the radial average amplitude spectra of the total magnetic anomaly and the magnetization intensity of the magnetic source body; |k| is the circular wave number; Z t is the upper interface of the magnetic layer; Z b is the lower interface of the magnetic layer, that is, the Curie point to be calculated; B is a constant, which has no effect on the calculation result of the Curie point.

[0129] The crustal magnetization intensity generally presents fractal characteristics, which can be expressed by formula 7, which is as follows:

[0130] A M (|k|)∝|k| (β-1) / 2

[0131] Here, β is the three-dimensional power spectrum index.

[0132] Combining Formula 6 and Formula 7, by taking the natural logarithm of both sides of Formula 6 and simplifying the process in different wave number domains, the following fitting formulas (Formula 8 and Formula 9) can be obtained:

[0133] Formula 8: ln[|k| (β-1) / 2 ×A ΔT (|k|)]≈lnD-|k|·Z t

[0134] Formula 9: ln[|k| (β-1) / 2 ×A ΔT (|k|) / |k|]≈lnE-|k|·Z0

[0135] Here, D and E are constants, and their values ​​do not affect the Curie point results. Z0 is the depth of the center of the magnetic layer. By using formula 8 and formula 9, the least squares principle is used to fit the straight line in the medium and high wave number domain and the low wave number domain respectively, and the corresponding Z t and Z0.

[0136] Finally, the bottom depth of the magnetic source body, that is, the depth Z of the inner surface b It can be obtained using the following simple formula (Formula 10):

[0137] Formula 10: Z b =2·Z0-Z t .

[0138] That is, the Curie surface depth fluctuation information can be obtained through formula 6-10.

[0139] In some embodiments, the crust temperature distribution information can be determined as follows:

[0140] Establish a one-dimensional steady-state heat conduction equation;

[0141] Based on the above aerial survey data and the above one-dimensional steady-state heat conduction equation, the crust temperature distribution information is determined.

[0142] As an example, by solving the one-dimensional steady-state heat conduction equation, the crustal temperature structure can be obtained. For a temperature-dependent thermal conductivity model and a model in which the heat production rate decays exponentially with increasing depth, the one-dimensional steady-state heat conduction equation can be expressed as follows (Equation 11):

[0143] Formula 11:

[0144] Where T is temperature; z is depth, positive downward; H0 is surface heat production, h r is the attenuation factor of heat production with depth. k is thermal conductivity, which is a function of temperature T(z). Its specific expression can be shown as follows (Formula 12):

[0145] Formula 12:

[0146] As an example, solving the one-dimensional steady-state heat conduction equation (Formula 11) requires two boundary conditions, namely the surface temperature and the Curie temperature. By substituting these two boundary conditions into the equation (Formula 11) and solving the equation (Formula 11) using the finite difference Runge-Kutta method, the temperature distribution inside the crust can be obtained, thereby obtaining the crust temperature distribution information.

[0147] In some embodiments, the crust fracture distribution information can be determined as follows:

[0148] The above aeromagnetic data are processed to enhance the THVH value;

[0149] The above crustal fracture distribution information is determined based on the aeromagnetic data after enhancement processing.

[0150] As an example, by performing THVH boundary signal enhancement processing on aeromagnetic data, structural boundaries and fracture divisions can be effectively characterized.

[0151] As an example, the THVH method can first calculate the horizontal derivative of the gravity and magnetic anomaly field to ensure the accuracy of the identification result in the horizontal position; then calculate the first-order vertical derivative of the horizontal derivative to enhance the vertical resolution, and finally calculate its tilt gradient to balance the boundary information presented by the strong and weak amplitudes. The THVH method expression can be as follows (Formula 13):

[0152] Formula 13:

[0153] in, THDR is the horizontal derivative of aeromagnetic anomalies. Formula 13 reduces the impact of high-frequency noise to a certain extent, significantly improves the vertical resolution, and also improves the ability to identify detailed information. The new boundary recognition technology based on the THVH method can not only balance the expression of strong and weak information, but also depict complex superimposed and hidden structures, which is more intuitive and visible. In other words, in this way, the fault distribution information can be obtained more accurately.

[0154] See also Figure 3 Based on the same inventive concept, the embodiment of the present application provides a device 300 for surveying a hot dry rock target area based on aeromagnetic data, comprising:

[0155] An acquisition unit 301 is used to acquire aeromagnetic data of the survey area;

[0156] A first determining unit 302 is used to determine at least rock mass distribution information, Curie point depth fluctuation information, crust temperature distribution information and crust fracture distribution information based on the above-mentioned aeromagnetic data;

[0157] The second determination unit 303 is used to determine the hot dry rock target area information of the exploration area by using the determined rock mass distribution information, the Curie point depth fluctuation information, the crust temperature distribution information and the crust fracture distribution information;

[0158] The above hot dry rock target area information includes:

[0159] Possible distribution areas of hot dry rock targets and the corresponding probabilities of each distribution area.

[0160] Optionally, in one embodiment, the second determining unit 303 is further configured to: determine a first area corresponding to the rock mass distribution information meeting a first preset condition;

[0161] Determine a second area corresponding to which the Curie surface depth fluctuation information meets a second preset condition;

[0162] Determine a third region corresponding to which the crust temperature distribution information meets a third preset condition;

[0163] Determine a fourth region corresponding to which the crustal fracture distribution information meets a fourth preset condition;

[0164] Based on the first area, the second area, the third area and the fourth area, possible distribution areas of the hot dry rock target area are determined.

[0165] Optionally, in one embodiment, the first area corresponds to a first weight value, the second area corresponds to a second weight value, and the third area and the fourth area correspond to a third weight value;

[0166] The above-mentioned first weight value is greater than the above-mentioned second weight value, and the above-mentioned second weight value is greater than the above-mentioned third weight value; and the second determination unit 303 is specifically used to: in response to determining that the weight value of a certain area is greater than a preset weight value, determine that the area is a distribution area, and determine the probability corresponding to the distribution area according to the size of the weight value.

[0167] Optionally, in one embodiment, the device 300 for surveying a hot dry rock target area using aeromagnetic data is further used to determine rock mass distribution information in the following manner:

[0168] Based on the above aeromagnetic data, determine the three-dimensional distribution information of crustal magnetization intensity;

[0169] The above rock mass distribution information is determined based on the three-dimensional distribution information of the crustal magnetization intensity.

[0170] Optionally, in one embodiment, the device 300 for surveying a hot dry rock target area using aeromagnetic data is further used to determine the Curie surface depth fluctuation information in the following manner:

[0171] Based on the above aeromagnetic data, the temperature interface that corresponds to the demagnetization temperature of ilmenite is determined;

[0172] Based on the determined temperature interface, the Curie surface depth fluctuation information is determined.

[0173] Optionally, in one embodiment, the device 300 for surveying a hot dry rock target area using aeromagnetic data is further used to determine the crust temperature distribution information in the following manner:

[0174] Establish a one-dimensional steady-state heat conduction equation;

[0175] Based on the above aerial survey data and the above one-dimensional steady-state heat conduction equation, the crust temperature distribution information is determined.

[0176] Optionally, in one embodiment, the device 300 for surveying a hot dry rock target area using aeromagnetic data is further used to determine the crustal fracture distribution information in the following manner:

[0177] The above aeromagnetic data are processed to enhance the THVH value;

[0178] The above crustal fracture distribution information is determined based on the aeromagnetic data after enhancement processing.

[0179] See also Figure 4 Based on the same concept, the embodiment of the present application also provides an electronic device 400 that applies the above method for surveying hot dry rock target areas based on aeromagnetic data. The electronic device includes: a memory 402, a processor 401, and a computer program 403 stored in the memory and executable on the processor.

[0180] The electronic device 400 may be a personal computer, a notebook computer, or the like.

[0181] Those skilled in the art will understand that Figure 4 The electronic device 400 is merely an example and does not constitute a limitation on the electronic device 400 , and may include more or less components than those shown in the figure, or may combine certain components, or may include different components.

[0182] The processor 401 may be a central processing unit (CPU), or other general-purpose processors, digital signal processors (DSP), application-specific integrated circuits (ASIC), field-programmable gate arrays (FPGA) or other programmable logic devices, discrete gate or transistor logic devices, discrete hardware components, etc. A general-purpose processor may be a microprocessor or any conventional processor, etc.

[0183] In some embodiments, the memory 402 may be an internal storage unit of the electronic device 400, such as a hard disk or memory of the electronic device 400. In other embodiments, the memory 402 may also be an external storage device of the electronic device 400, such as a plug-in hard disk, a smart memory card (Smart Media Card, SMC), a secure digital (Secure Digital, SD) card, a flash card (Flash Card), etc. equipped on the electronic device 400. Further, the memory 402 may also include both an internal storage unit of the electronic device 400 and an external storage device.

[0184] It should be noted that the information interaction, execution process, etc. between the above-mentioned devices / units are based on the same concept as the method embodiment of the present application. Their specific functions and technical effects can be found in the method embodiment part and will not be repeated here.

[0185] The technicians in the relevant field can clearly understand that for the convenience and simplicity of description, only the division of the above-mentioned functional units and modules is used as an example for illustration. In practical applications, the above-mentioned function allocation can be completed by different functional units and modules as needed, that is, the internal structure of the above-mentioned device can be divided into different functional units or modules to complete all or part of the functions described above. The functional units and modules in the embodiment can be integrated in a processing unit, or each unit can exist physically separately, or two or more units can be integrated in one unit. The above-mentioned integrated unit can be implemented in the form of hardware or in the form of software functional units. In addition, the specific names of the functional units and modules are only for the convenience of distinguishing each other, and are not used to limit the scope of protection of this application. The specific working process of the units and modules in the above-mentioned system can refer to the corresponding process in the aforementioned method embodiment, which will not be repeated here.

[0186] An embodiment of the present application further provides a computer-readable storage medium, wherein the computer-readable storage medium stores a computer program, and when the computer program is executed by a processor, the steps in the above-mentioned method embodiments can be implemented.

[0187] An embodiment of the present application provides a computer program product. When the computer program product runs on a mobile terminal, the mobile terminal can implement the steps in the above-mentioned method embodiments when executing the computer program product.

[0188] If the above-mentioned integrated unit is implemented in the form of a software functional unit and sold or used as an independent product, it can be stored in a computer-readable storage medium. Based on this understanding, the present application implements all or part of the processes in the above-mentioned embodiment method, which can be completed by instructing the relevant hardware through a computer program. The above-mentioned computer program can be stored in a computer-readable storage medium, and the computer program can implement the steps of the above-mentioned various method embodiments when executed by the processor. Among them, the above-mentioned computer program includes computer program code, and the above-mentioned computer program code can be in source code form, object code form, executable file or some intermediate form. The above-mentioned computer-readable medium may at least include: any entity or device that can carry the computer program code to the camera / electronic device, recording medium, computer memory, read-only memory (ROM, Read-Only Memory), random access memory (RAM, RandomAccess Memory), electric carrier signal, telecommunication signal and software distribution medium. For example, a USB flash drive, a mobile hard disk, a magnetic disk or an optical disk. In some jurisdictions, according to legislation and patent practice, computer-readable media cannot be electric carrier signals and telecommunication signals.

[0189] In the above embodiments, the description of each embodiment has its own emphasis. For parts that are not described or recorded in detail in a certain embodiment, reference can be made to the relevant descriptions of other embodiments.

[0190] Those of ordinary skill in the art will appreciate that the units and algorithm steps of each example described in conjunction with the embodiments disclosed herein can be implemented in electronic hardware, or a combination of computer software and electronic hardware. Whether these functions are performed in hardware or software depends on the specific application and design constraints of the technical solution. Professional and technical personnel can use different methods to implement the described functions for each specific application, but such implementation should not be considered to be beyond the scope of this application.

[0191] In the embodiments provided in the present application, it should be understood that the disclosed devices / network equipment and methods can be implemented in other ways. For example, the device / network equipment embodiments described above are merely schematic. For example, the division of the above modules or units is only a logical function division. There may be other division methods in actual implementation, such as multiple units or components can be combined or integrated into another system, or some features can be ignored or not executed. Another point is that the mutual coupling or direct coupling or communication connection shown or discussed can be through some interfaces, indirect coupling or communication connection of devices or units, which can be electrical, mechanical or other forms.

[0192] The units described above as separate components may or may not be physically separated, and the components shown as units may or may not be physical units, that is, they may be located in one place or distributed on multiple network units. Some or all of the units may be selected according to actual needs to achieve the purpose of the solution of this embodiment.

[0193] The above embodiments are only used to illustrate the technical solutions of the present application, rather than to limit them. Although the present application has been described in detail with reference to the above embodiments, those skilled in the art should understand that they can still modify the technical solutions described in the above embodiments, or replace some of the technical features therein by equivalents. These modifications or replacements do not deviate the essence of the corresponding technical solutions from the spirit and scope of the technical solutions of the embodiments of the present application, and should all be included in the protection scope of the present application.

Claims

1. A method for exploring hot dry rock target areas based on aeromagnetic data, characterized in that: include: Obtain aeromagnetic data of the survey area; Based on the aeromagnetic data, at least rock mass distribution information, Curie point depth fluctuation information, crust temperature distribution information and crust fracture distribution information are determined; Determine the hot dry rock target area information of the exploration area by using the determined rock mass distribution information, the Curie point depth fluctuation information, the crust temperature distribution information and the crust fracture distribution information; The hot dry rock target area information includes: Possible distribution areas of hot dry rock targets and the corresponding probabilities of each distribution area; Among them, the rock mass distribution information is determined by the following methods: Based on the aeromagnetic data, determine the three-dimensional distribution information of crustal magnetization intensity; Determining the rock mass distribution information according to the three-dimensional distribution information of the crustal magnetization intensity; Wherein, the Curie surface depth fluctuation information is determined by the following method: Based on the aeromagnetic data, determining a temperature interface that meets the demagnetization temperature of ilmenite; Determining Curie surface depth fluctuation information based on the determined temperature interface; The crust temperature distribution information is determined in the following manner: Establish a one-dimensional steady-state heat conduction equation; Determining crustal temperature distribution information based on the aeromagnetic data and the one-dimensional steady-state heat conduction equation; The crustal fracture distribution information is determined in the following manner: Performing THVH value enhancement processing on the aeromagnetic data; Determining the crustal fracture distribution information according to the enhanced aeromagnetic data; The step of performing THVH value enhancement processing on the aeromagnetic data includes: Based on the horizontal derivative of the gravity and magnetic anomaly field, the first-order vertical derivative of the horizontal derivative and the tilt gradient value, the aeromagnetic data is enhanced; Among them, the THVH value is determined based on the following formula: ， in, , Used for the horizontal derivative of gravity and magnetic anomaly fields.

2. The method according to claim 1, characterized in that The method of determining the hot dry rock target area information of the exploration area by using the determined rock mass distribution information, the Curie point depth fluctuation information, the crust temperature distribution information and the crust fracture distribution information includes: Determine a first area corresponding to which the rock mass distribution information meets a first preset condition; Determine a second area corresponding to which the Curie surface depth fluctuation information meets a second preset condition; Determine a third region corresponding to which the crust temperature distribution information meets a third preset condition; Determining a fourth region corresponding to which the crustal fracture distribution information meets a fourth preset condition; Based on the first area, the second area, the third area and the fourth area, a possible distribution area of ​​the hot dry rock target area is determined.

3. The method according to claim 2, characterized in that The first area corresponds to a first weight value, the second area corresponds to a second weight value, and the third area and the fourth area correspond to a third weight value; The first weight value is greater than the second weight value, and the second weight value is greater than the third weight value; And, determining the possible distribution area of ​​the hot dry rock target area based on the first area, the second area, the third area and the fourth area includes: In response to determining that the weight value of a certain area is greater than a preset weight value, the area is determined to be a distribution area, and the probability corresponding to the distribution area is determined according to the size of the weight value.

4. A device for exploring hot dry rock target areas based on aeromagnetic data, characterized in that: include: An acquisition unit, used to acquire aeromagnetic data of the survey area; A first determination unit is used to determine at least rock mass distribution information, Curie point depth fluctuation information, crust temperature distribution information and crust fracture distribution information based on the aeromagnetic data; A second determination unit is used to determine the hot dry rock target area information of the exploration area by using the determined rock mass distribution information, the Curie point depth fluctuation information, the crust temperature distribution information and the crust fracture distribution information; The hot dry rock target area information includes: Possible distribution areas of hot dry rock targets and the corresponding probabilities of each distribution area; The device is also used to determine rock mass distribution information in the following manner: Based on the aeromagnetic data, determine the three-dimensional distribution information of the crustal magnetization intensity; according to the three-dimensional distribution information of the crustal magnetization intensity, determine the rock mass distribution information; Wherein, the device is further used to determine the Curie surface depth fluctuation information in the following manner: Based on the aeromagnetic data, determining a temperature interface that meets the demagnetization temperature of ilmenite; Determining Curie surface depth fluctuation information based on the determined temperature interface; Wherein, the device is also used to determine the crust temperature distribution information in the following manner: Establish a one-dimensional steady-state heat conduction equation; Determining crustal temperature distribution information based on the aeromagnetic data and the one-dimensional steady-state heat conduction equation; Wherein, the device is also used to determine the crust fracture distribution information in the following manner: Performing THVH value enhancement processing on the aeromagnetic data; Determining the crustal fracture distribution information according to the enhanced aeromagnetic data; The step of performing THVH value enhancement processing on the aeromagnetic data includes: Based on the horizontal derivative of the gravity and magnetic anomaly field, the first-order vertical derivative of the horizontal derivative and the tilt gradient value, the aeromagnetic data is enhanced; Among them, the THVH value is determined based on the following formula: ， in, , Used for the horizontal derivative of gravity and magnetic anomaly fields.

5. An electronic device comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, characterized in that: When the processor executes the computer program, the method according to any one of claims 1 to 3 is implemented.

6. A computer readable medium having a computer program stored thereon, characterized in that: When the program is executed by a processor, the method according to any one of claims 1 to 3 is implemented.

Images