Methods, devices, electronic equipment and storage media for deep geothermal reservoir well site target area exploration

Abstract

This application relates to the field of geological exploration, particularly to methods, apparatus, electronic equipment, and storage media for well site target area surveying in deep geothermal reservoirs. The method includes: acquiring gravity data of the area to be surveyed; determining buried hills and / or uplift zones based on the gravity data, and designating these zones as preferred target areas; determining the magnetotelluric signals of the preferred target areas using magnetotelluric sounding; and determining the electrical data of the preferred target areas based on the magnetotelluric signals. This application improves the efficiency of well site selection.

Description

Technical Field

[0001] This application relates to the field of geological exploration, and in particular to methods, apparatus, electronic equipment and storage media for the exploration of target areas at deep geothermal reservoir well sites. Background Technology

[0002] Deep carbonate reservoirs are widely distributed in the North China Plain, forming an important geothermal reservoir layer. Their burial depths are extremely uneven, ranging from hundreds to thousands of meters. Currently, due to limitations in extraction costs and technical capabilities, the main depth for geothermal resource extraction is shallower than 4000 meters. Selecting carbonate reservoirs with suitable burial depths is the primary objective for geothermal well selection in the North China Plain. Traditional well location selection methods largely rely on the experience of technical personnel. However, these traditional methods are relatively inefficient. Summary of the Invention

[0003] To improve the efficiency of well location selection, this application provides a method, apparatus, electronic equipment, and storage medium for deep geothermal reservoir well location target area exploration.

[0004] Firstly, this application provides a method for exploring target areas at deep thermal reservoir well locations, employing the following technical solution:

[0005] A method for well location target area exploration in deep geothermal reservoirs includes:

[0006] Obtain gravity data for the area to be tested;

[0007] Based on the gravity data, the buried hill zone and / or the raised area are determined, and the buried hill zone and / or the raised area are selected as the preferred target area;

[0008] The magnetotelluric signal of the preferred target area is determined based on magnetotelluric sounding.

[0009] The electrical data of the preferred target area are determined based on the magnetotelluric signal.

[0010] By adopting the above technical solution, gravity data of the area to be tested is obtained, and gravity anomalies are determined based on the gravity data. High gravity anomalies are buried hills and / or uplift areas. The entire area to be tested is initially screened using gravity data to identify potential target areas for energy storage. Then, the electrical data within the target areas are analyzed to further analyze the underground structure of the target areas, thereby determining the well depth and location. This preliminary filtering method improves the efficiency of well location selection.

[0011] In another possible implementation, after determining the magnetotelluric signal of the preferred target area based on magnetotelluric sounding, the method further includes:

[0012] The magnetotelluric signals of the preferred target area are preprocessed, including noise suppression, data filtering, static correction, terrain correction, and resolution matching.

[0013] In another possible implementation, determining the electrical data of the preferred target area based on the magnetotelluric signal includes: performing data inversion on the magnetotelluric signal, wherein the magnetotelluric signal used for data inversion is a TE mode magnetotelluric signal.

[0014] In another possible implementation, the data inversion of the magnetotelluric signal includes:

[0015] The magnetotelluric signal is used to generate an initial model through one-dimensional inversion based on the Occam inversion method;

[0016] Based on the initial model, two-dimensional inversion is performed to generate electrical data for the preferred target area.

[0017] In another possible implementation, the gravity data of the region to be measured is a Bouguer gravity anomaly map of the region to be measured.

[0018] In another possible implementation, the surveying method further includes:

[0019] The gravity data of the preferred target area is extended upward, and the first-order horizontal derivatives in different directions are determined to generate the first-order derivative maps of Bouguer gravity anomalies in different directions and the upward extension map of Bouguer gravity anomalies.

[0020] Based on the derivative diagram, the Bouguer gravity anomaly upward extension diagram, and the generated gravity-inferred fracture structure diagram.

[0021] By adopting the above technical solution, the distribution of fracture structures in the preferred target area can be reflected by generating a gravity-inferred fracture structure map. In addition to electrical data, the preferred target area can be interpreted in another way, which facilitates the selection of well locations by the staff.

[0022] Secondly, this application provides a deep geothermal reservoir well location target area exploration device, which adopts the following technical solution:

[0023] A deep geothermal reservoir well location target area exploration device, comprising:

[0024] The acquisition module is used to acquire gravity data for the area to be tested.

[0025] The target area determination module is used to determine the buried hill and / or protrusion area based on the gravity data, and to select the buried hill and / or protrusion area as the preferred target area.

[0026] The re-measurement module is used to determine the magnetotelluric signal of the preferred target area based on magnetotelluric sounding.

[0027] The evaluation module is used to determine the electrical data of the preferred target area based on the magnetotelluric signal.

[0028] In one possible implementation, the surveying device further includes a preprocessing module, which is used to: after determining the magnetotelluric signal of the preferred target area based on magnetotelluric sounding, preprocess the magnetotelluric signal of the preferred target area, the preprocessing including: noise suppression, data filtering, static correction, terrain correction and resolution matching.

[0029] In one possible implementation, when the evaluation module determines the electrical data of the preferred target area based on the magnetotelluric signal, it is specifically used to: perform data inversion on the magnetotelluric signal, wherein the magnetotelluric signal used for data inversion is a TE mode magnetotelluric signal.

[0030] In one possible implementation, the evaluation module, when performing data inversion, is specifically used to: perform one-dimensional inversion of the magnetotelluric signal based on the Occam inversion method to generate an initial model;

[0031] Based on the initial model, two-dimensional inversion is performed to generate electrical data for the preferred target area.

[0032] In one possible implementation, the gravity data of the area to be measured is a Bouguer gravity anomaly map of the area to be measured.

[0033] In one possible implementation, the surveying device further includes a generation module, which is used to: extend the gravity data of the preferred target area upward and determine the first-order horizontal derivatives in different directions to generate a first-order derivative map of Bouguer gravity anomaly in different directions and an upward extension map of Bouguer gravity anomaly.

[0034] Based on the derivative diagram, the Bouguer gravity anomaly upward extension diagram, and the generated gravity-inferred fracture structure diagram.

[0035] By adopting the above technical solution, the distribution of fracture structures in the preferred target area can be reflected by generating a gravity-inferred fracture structure map. In addition to electrical data, the preferred target area can be interpreted in another way, which facilitates the selection of well locations by the staff.

[0036] Thirdly, this application provides an electronic device that adopts the following technical solution:

[0037] An electronic device comprising:

[0038] At least one processor;

[0039] Memory;

[0040] At least one application, wherein the at least one application is stored in memory and configured to be executed by at least one processor, the at least one application being configured to: execute the aforementioned method for deep thermal reservoir well location target area exploration.

[0041] Fourthly, this application provides a computer-readable storage medium, which adopts the following technical solution:

[0042] A computer-readable storage medium includes: a computer program stored thereon that can be loaded by a processor and executed as described above for a deep thermal reservoir well location target area exploration method.

[0043] In summary, this application includes the following beneficial technical effects:

[0044] By acquiring gravity data of the area to be tested, gravity anomalies are identified based on the gravity data. High gravity anomalies are identified as buried hills and / or uplift areas. The entire area to be tested is initially screened using gravity data to identify potential target areas for energy storage. Then, the electrical data within the target areas are analyzed to further analyze the underground structure of the target areas, thereby determining the well depth and location. This preliminary filtering method improves the efficiency of well location selection. Attached Figure Description

[0045] Figure 1 This is a flowchart illustrating the deep geothermal reservoir well location target area exploration method according to an embodiment of this application;

[0046] Figure 2 This is a block diagram of the deep thermal reservoir well location target area exploration device according to an embodiment of this application;

[0047] Figure 3 This is a schematic diagram of an electronic device according to an embodiment of this application. Detailed Implementation

[0048] The following is in conjunction with the appendix Figure 1-3 This application will be described in further detail.

[0049] To make the objectives, technical solutions, and advantages of the embodiments of this application clearer, the technical solutions of the embodiments of this application will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some embodiments of this application, not all embodiments. Based on the embodiments of this application, all other embodiments obtained by those skilled in the art without creative effort are within the scope of protection of this application.

[0050] Furthermore, the term "and / or" in this article is merely a description of the relationship between related objects, indicating that three relationships can exist. For example, A and / or B can represent: A existing alone, A and B existing simultaneously, or B existing alone. Additionally, the character " / " in this article, unless otherwise specified, generally indicates that the preceding and following related objects have an "or" relationship.

[0051] To facilitate understanding of the technical solutions proposed in this application, several elements that will be introduced in the description of this application are first presented here. It should be understood that the following description is only for the purpose of understanding these elements and the content of the embodiments of this application, and does not necessarily cover all possible situations.

[0052] (1) Gravity anomaly: Earth’s space gravity anomaly refers to the gravity difference obtained by subtracting the normal gravity field value from the measured gravity value on the natural surface of the Earth after height correction. Space gravity anomaly is caused by the uneven distribution of the density of the crust and the upper mantle, reflecting the difference between the actual shape and mass distribution of the Earth and the theoretical reference ellipsoid.

[0053] (2) Buried hills: Buried hills are an important deep geological unit, referring to ancient landforms or tectonic uplifts that were buried underground and formed before the basin received sediments. They were subsequently covered and buried by new strata during the basin's evolution. Buried hills are cut by faults, resulting in a large number of fractures. These fractures extend to deeper layers under the action of weathering and erosion, forming expanded reservoir spaces, which are important reservoirs. Multiple phases and multiple directions of fractures organically connect various reservoir spaces and can serve as effective conductors. Therefore, buried hill fractures are rich, productive, and complex reservoir systems.

[0054] (3) Uplifted areas: Uplifted areas and slope areas are the dominant migration directions of oil and gas and are key areas for oil and gas exploration.

[0055] This application provides a method for deep geothermal reservoir well location target area exploration, referring to... Figure 1 Performed by an electronic device, the method includes:

[0056] Step S101: Obtain gravity data for the area to be tested.

[0057] The gravity data consists of regional Bouguer gravity anomaly maps of the area to be measured, ranging from 1:50,000 to 1:200,000. Studying density differences between different geological bodies is a geophysical prerequisite for gravity exploration. The division of density interfaces is closely related to the division of tectonic layers. Measuring the density of strata (rocks) and analyzing the results are important aspects of gravity exploration.

[0058] Based on measured density data of exposed rocks in the area under test and adjacent areas, this study comprehensively analyzes the density characteristics of stratigraphic boundaries and inter-strata, elucidating the density characteristics of strata and intrusive rocks and their relationship with gravity anomalies. This provides a reliable geophysical basis for the interpretation of gravity data and reservoir exploration in the area under test. The delineation of density boundaries is closely related to regional tectonics and reservoir structures, and calculating the undulations or depth variations of density boundaries is of great significance in regional tectonic studies and reservoir exploration.

[0059] Taking the North China Plain as an example, based on the density of the strata and the characteristics of rock formation, the strata in the region are divided into four density layers from oldest to newest.

[0060] The first density layer comprises Neoarchean and Paleoproterozoic metamorphic rocks, with an average density of 2.71 × 10³ kg / m³. 3 It is a high-density layer. It is a set of high-altitude crystalline metamorphic rocks that form a crystalline basement. Its protolith is mainly intermediate to basic volcanic rocks, and the constituent rocks are mostly gneiss, granulite, and various granulites and schists containing amphibole, pyroxene, etc.

[0061] Second density layer: average density is 2.66 × 10³ kg / m³ 3 This is a medium-density layer, comprising various formations from the Middle to Neoproterozoic and Paleozoic strata. It is a set of marine and transitional marine sedimentary rocks, including various limestones, dolomite, sandstone, shale, and mudstones. It is predominantly composed of carbonate rocks such as limestone and dolomite.

[0062] Third density layer: average density 2.47 × 10³ kg / m³ 3 This is a low-density layer. It includes various Mesozoic, Neogene, and Paleogene formations. It is a set of terrestrial sedimentary rocks, volcanic clastic rocks, and locally distributed volcanic rocks. The diagenesis is relatively recent, with high porosity and low density.

[0063] Fourth density layer: average density only 2.05 × 10³ kg / m³ 3 This is the lowest density layer, mainly composed of loose Quaternary sediments.

[0064] In addition, the average density of the deep silica-alumina layer in the Earth's crust is 2.78 × 10³ kg / m³. 3 The average density of the silicon-magnesium layer is 3.18 × 10³ kg / m³. 3 The average density of the upper mantle is 3.48 × 10³ kg / m³. 3 The layers are separated by the Conrad and Mohorovičić discontinuities, which are the most prominent density interfaces in the Earth's crust. Their variations are the main cause of regional variations in the gravitational field.

[0065] Carbonate rocks, as the second density layer of the North China Plain, are mainly composed of Paleozoic Ordovician and Cambrian limestones, as well as Middle and Upper Proterozoic dolomite from the Wumishan and Gaoyuzhuang Formations of the Jixian System.

[0066] Step S102: Determine the buried hill and / or raised area based on gravity data, and use the buried hill and / or raised area as the preferred target area.

[0067] Among them, the high gravity area corresponds to the uplift area or the buried hill zone of the structure, and this area is preferred as the drilling target area. The high gravity area is the area where the gravity is greater than the preset value.

[0068] Step S103: Determine the magnetotelluric signal of the preferred target area based on magnetotelluric sounding.

[0069] Among them, magnetotelluric sounding (MT) is an important method in the field of electromagnetic exploration. This method uses the Earth's natural alternating electromagnetic field to explore the underground electrical structure. When electromagnetic waves propagate in underground rock strata, electromagnetic waves of different frequencies have different penetration depths. By observing these electromagnetic fields on the ground, we can understand the distribution of electrical properties in the underground medium, thereby realizing the detection of the electrical structure inside the Earth.

[0070] In magnetotelluric sounding theory, the electromagnetic waves generated by the field source are considered as plane waves incident vertically downwards. The propagation and penetration depth of electromagnetic waves in the subsurface medium mainly depend on the frequency (o) of the electromagnetic wave and the conductivity (σ) of the subsurface medium. The lower the frequency of the electromagnetic field and the higher the conductivity of the medium, the less energy the electromagnetic wave loses during propagation in the subsurface medium, and the deeper its penetration depth. This constitutes the physical basis of magnetotelluric sounding.

[0071] During the exploration, the magnetotelluric sounding data acquisition device is placed on the ground to observe the changing electromagnetic field. Specifically, multiple test points are set up in the selected target area at equal intervals of 500 meters and a depth of 5 km to conduct multi-point exploration tests.

[0072] Then, the signal is converted from the time domain to the frequency domain. The impedance tensor is determined based on the frequency domain signal, and the resistivity and phase are determined based on the impedance tensor. In the two-dimensional geological model, the conductivity varies only along a single cross-section, while the conductivity along the structural strike direction is defined as constant. Under this assumption, the electric and magnetic fields are orthogonal, and the equivalent impedance tensor can be decomposed into two independent modes: the TE mode, which describes the response of current flowing along the structural direction, and the TM mode, which describes the response of current perpendicular to the structural direction.

[0073] Step S104: Determine the electrical data of the preferred target area based on the magnetotelluric signal.

[0074] Fault structures often serve as conduits for the migration of deep geothermal energy to the surface. Because fractured structures are prone to water accumulation or filling with clay, the resistivity at fault zones is generally low. Therefore, analyzing the distribution of underground resistivity can reveal the distribution of fault structures. Faults provide favorable migration channels for geothermal energy, making the areas surrounding fault zones promising areas for geothermal development. Studies have found numerous fault structures distributed along the sides of Mesozoic basins, making these areas important for geothermal exploration and development. Therefore, magnetotelluric sounding can be used to determine the distribution of stratigraphic resistivity. After determining the regional electrical trends of tectonic units and different frequency bands, electrical data for the optimal target area are obtained through inversion using segmented electrical trend data along the profile.

[0075] The electrical data of the selected target area characterizes the deep electrical structure in detail, revealing the deep extension of the main fault zones along the profile, the characteristics of the deep electrical structure, and the development of low-resistivity layers within the shell. This enables the inference of thermal reservoir analysis data based on resistivity values, thereby facilitating the rapid location of wells by the staff.

[0076] Based on the deep geothermal reservoir well location target area exploration method provided in this application, gravity data of the area to be tested is obtained, and gravity anomalies are determined based on the gravity data. Among them, high gravity anomalies are buried hills and / or uplift areas. The entire area to be tested is initially screened using gravity data to determine the preferred target areas that may contain energy storage. Then, the electrical data in the preferred target areas are analyzed to further analyze the underground structure of the preferred target areas, thereby determining the well depth and well location. The efficiency of well location selection is improved through the preliminary filtering method.

[0077] In a specific embodiment, after step S103, step S105 is further included: preprocessing the magnetotelluric signal of the preferred target area, the preprocessing including: noise suppression, data filtering, static correction, terrain correction and resolution matching.

[0078] Specifically, noise suppression refers to removing noise from magnetotelluric signals; data filtering refers to selecting data from magnetotelluric signals that meet the requirements according to set criteria; static correction refers to translating and realigning magnetotelluric signal curves affected by static effects to remove the influence of shallow surface inhomogeneities; topographic correction refers to correcting receiving points at different heights on the survey line to the reference plane to eliminate the gravity effects caused by topographic undulations; and resolution matching refers to improving the resolution of magnetotelluric signal curves to the same resolution as another high-resolution curve through filtering and linear regression to increase the density of longitudinal data.

[0079] In one specific embodiment, in step S104, determining the electrical data of the preferred target area based on the magnetotelluric signal includes: performing data inversion on the magnetotelluric signal, wherein the magnetotelluric signal used for data inversion is a TE mode magnetotelluric signal. Here, data inversion refers to inferring the formation resistivity based on the obtained magnetotelluric signal.

[0080] In a specific embodiment, step S104 involves data inversion of the magnetotelluric signal, including: performing one-dimensional inversion of the magnetotelluric signal based on the Occam inversion method to generate an initial model; and performing two-dimensional inversion based on the initial model to generate electrical data of the preferred target area.

[0081] In a specific embodiment, the survey method further includes step S105 (not shown in the figure), which is set after step S102, and specifically includes: extending the gravity data of the preferred target area upward and determining the first-order derivatives in different directions to generate the first-order derivative map of Bouguer gravity anomaly in different directions and the upward extension map of Bouguer gravity anomaly; and generating a gravity comprehensive inference fracture structure map based on the derivative map, the upward extension map of Bouguer gravity anomaly, and the gravity comprehensive inference fracture structure map.

[0082] Specifically, this includes generating a gravity-inferred fracture structure map by utilizing the linear gradient zone of Bouguer gravity anomaly, the boundary line of different gravity field characteristics, the linear extreme zone of the horizontal derivative, and the linear bending, truncation, local expansion, or convergence fracture structures of gravity anomalies.

[0083] Reference Figure 2 The above embodiments describe a method for deep geothermal reservoir well location target area exploration from the perspective of process flow. The following embodiments describe a deep geothermal reservoir well location target area exploration device 100 from the perspective of virtual module or virtual unit. For details, please refer to the following embodiments.

[0084] A deep geothermal reservoir well location target area exploration device 100, comprising:

[0085] Module 1001 is used to acquire gravity data of the area to be tested.

[0086] The target area determination module 1002 is used to determine the buried hill and / or the raised area based on gravity data, and to select the buried hill and / or the raised area as the preferred target area.

[0087] The re-measurement module 1003 is used to determine the magnetotelluric signal of the preferred target area based on magnetotelluric sounding.

[0088] Evaluation module 1004 is used to determine the electrical data of the preferred target area based on magnetotelluric signals.

[0089] In one possible implementation, the surveying device further includes a preprocessing module, which is used to: after determining the magnetotelluric signal of the preferred target area based on magnetotelluric sounding, preprocess the magnetotelluric signal of the preferred target area, including noise suppression, data filtering, static correction, terrain correction, and resolution matching.

[0090] In one possible implementation, when the evaluation module 1004 determines the electrical data of the preferred target area based on the magnetotelluric signal, it is specifically used to: perform data inversion on the magnetotelluric signal, wherein the magnetotelluric signal used for data inversion is a TE mode magnetotelluric signal.

[0091] In one possible implementation, the evaluation module 1004, when performing data inversion, is specifically used to: perform one-dimensional inversion of the magnetotelluric signal based on the Occam inversion method to generate an initial model;

[0092] Based on the initial model, two-dimensional inversion is performed to generate electrical data for the preferred target area.

[0093] In one possible implementation, the gravity data for the region to be measured is a Bouguer gravity anomaly map of the region to be measured.

[0094] In one possible implementation, the surveying device further includes a generation module, which is used to: extend the gravity data of the preferred target area upward and determine the first-order horizontal derivatives in different directions to generate a first-order derivative map of Bouguer gravity anomaly in different directions and an upward extension map of Bouguer gravity anomaly.

[0095] Based on the derivative diagram, the Bouguer gravity anomaly upward extension diagram, and the generated gravity comprehensive inference fracture structure diagram.

[0096] This application also describes an electronic device from the perspective of a physical device, such as... Figure 3 As shown, Figure 3 The illustrated electronic device 1100 includes a processor 1101 and a memory 1103. The processor 1101 and the memory 1103 are connected, for example, via a bus 1102. Optionally, the electronic device 1100 may also include a transceiver 1104. It should be noted that in practical applications, the transceiver 1104 is not limited to one unit, and the structure of this electronic device 1100 does not constitute a limitation on the embodiments of this application.

[0097] Processor 1101 may be a CPU (Central Processing Unit), a general-purpose processor, a DSP (Digital Signal Processor), an ASIC (Application Specific Integrated Circuit), an FPGA (Field Programmable Gate Array), or other programmable logic devices, transistor logic devices, hardware components, or any combination thereof. It can implement or execute the various exemplary logic blocks, modules, and circuits described in conjunction with the disclosure of this application. Processor 1101 may also be a combination that implements computational functions, such as including one or more microprocessor combinations, a combination of a DSP and a microprocessor, etc.

[0098] Bus 1102 may include a pathway for transmitting information between the aforementioned components. Bus 1102 may be a PCI (Peripheral Component Interconnect) bus or an EISA (Extended Industry Standard Architecture) bus, etc. Bus 1102 can be divided into address bus, data bus, control bus, etc. For ease of representation, Figure 3 The bus is represented by a single thick line, but this does not mean that there is only one bus or one type of bus.

[0099] The memory 1103 may be a ROM (Read Only Memory) or other type of static storage device capable of storing static information and instructions, RAM (Random Access Memory) or other type of dynamic storage device capable of storing information and instructions, or an EEPROM (Electrically Erasable Programmable Read Only Memory), CD-ROM (Compact Disc Read Only Memory) or other optical disc storage, optical disc storage (including compressed optical discs, laser discs, optical discs, digital universal optical discs, Blu-ray discs, etc.), magnetic disk storage media or other magnetic storage devices, or any other medium capable of carrying or storing desired program code in the form of instructions or data structures and accessible by a computer, but not limited thereto.

[0100] The memory 1103 is used to store application code that executes the solution of this application, and its execution is controlled by the processor 1101. The processor 1101 is used to execute the application code stored in the memory 1103 to implement the content shown in the foregoing method embodiments.

[0101] Electronic devices include, but are not limited to: mobile terminals such as mobile phones, laptops, digital radio receivers, PDAs (personal digital assistants), PADs (tablet computers), PMPs (portable multimedia players), and in-vehicle terminals (such as in-vehicle navigation terminals), as well as fixed terminals such as digital TVs and desktop computers. Servers can also be included. Figure 3 The electronic device shown is merely an example and should not impose any limitation on the functionality and scope of use of the embodiments of this application.

[0102] It should be understood that although the steps in the flowcharts of the accompanying figures are shown sequentially as indicated by the arrows, these steps are not necessarily executed in the order indicated by the arrows. Unless explicitly stated herein, there is no strict order restriction on the execution of these steps, and they can be executed in other orders. Moreover, at least some steps in the flowcharts of the accompanying figures may include multiple sub-steps or multiple stages. These sub-steps or stages are not necessarily completed at the same time, but can be executed at different times, and their execution order is not necessarily sequential, but can be performed alternately or in turn with other steps or at least some of the sub-steps or stages of other steps.

[0103] The above are only some embodiments of this application. It should be noted that for those skilled in the art, several improvements and modifications can be made without departing from the principle of this application, and these improvements and modifications should also be considered within the scope of protection of this application.

Claims

1. A method for exploring target areas of deep thermal reservoir well locations, characterized in that, include: Obtain regional Bouguer gravity anomaly maps of the area to be tested at scales of 1:50,000 to 1:200,000. Based on the Bouguer gravity anomaly map, gravity anomaly points are determined, buried hills and / or uplifted areas are identified, and the buried hills and / or uplifted areas are selected as preferred target areas. The magnetotelluric signal of the preferred target area is determined based on magnetotelluric sounding. The magnetotelluric signal is converted from a time domain signal to a frequency domain signal. The impedance tensor is determined based on the frequency domain signal, and the TE mode magnetotelluric signal is decomposed from the impedance tensor. The magnetotelluric signal of the TE mode is used to generate an initial model by performing a one-dimensional inversion based on the Occam inversion method; Based on the initial model, two-dimensional inversion is performed to generate electrical data of the preferred target area, so as to analyze the underground structure of the preferred target area and determine the well depth and well location; The gravity data of the preferred target area is extended upward, and the first-order horizontal derivatives in different directions are determined to generate the first-order derivative maps of Bouguer gravity anomalies in different directions and the upward extension map of Bouguer gravity anomalies. Gravity-inferred fracture structure maps are generated by utilizing the linear gradient zones of Bouguer gravity anomalies, the boundaries of different gravity field characteristics, the linear extreme zones of horizontal derivatives, and the linear bending, truncation, local expansion, or convergence fracture structures of gravity anomalies in the same direction.

2. The method according to claim 1, characterized in that, After determining the magnetotelluric signal of the preferred target area based on magnetotelluric sounding, the process further includes: The magnetotelluric signals of the preferred target area are preprocessed, including noise suppression, data filtering, static correction, terrain correction, and resolution matching.

3. A deep geothermal reservoir well location target area exploration device, characterized in that, include: The acquisition module is used to acquire regional Bouguer gravity anomaly maps of the area to be tested, ranging from 1:50,000 to 1:200,000. The target area determination module is used to determine gravity anomaly points based on the Bouguer gravity anomaly map, determine buried hills and / or protrusions, and select the buried hills and / or protrusions as preferred target areas. The re-measurement module is used to determine the magnetotelluric signal of the preferred target area based on magnetotelluric sounding; convert the magnetotelluric signal from the time domain signal to the frequency domain signal, determine the impedance tensor based on the frequency domain signal, and decompose the TE mode magnetotelluric signal from the impedance tensor. The evaluation module is used to perform one-dimensional inversion on the magnetotelluric signal of the TE mode based on the Occam inversion method to generate an initial model; and to perform two-dimensional inversion on the initial model to generate electrical data of the preferred target area, so as to analyze the subsurface structure of the preferred target area and determine the well depth and well location. The generation module is used to extend the gravity data of the preferred target area upward and determine the first-order horizontal derivative in different directions to generate the first-order derivative map of Bouguer gravity anomaly in different directions and the upward extension map of Bouguer gravity anomaly. Gravity-inferred fracture structure maps are generated by utilizing the linear gradient zones of Bouguer gravity anomalies, the boundaries of different gravity field characteristics, the linear extreme zones of horizontal derivatives, and the linear bending, truncation, local expansion, or convergence fracture structures of gravity anomalies in the same direction.

4. An electronic device, characterized in that, The electronic device includes: At least one processor; Memory; At least one application, wherein the at least one application is stored in memory and configured to be executed by at least one processor, the at least one application being configured to: perform the deep thermal reservoir well location target area exploration method as described in claim 1 or 2.

5. A computer-readable storage medium having a computer program stored thereon, characterized in that, When the computer program is executed in the computer, it causes the computer to perform the deep thermal reservoir well location target area exploration method as described in claim 1 or 2.

Images