An airborne magnetotelluric data processing method, device, equipment and storage medium

Abstract

The application provides an airborne magnetotelluric data processing method, device and equipment and a storage medium, relates to the technical field of geophysical exploration, and comprises the following steps: arranging at least two ground base stations in a target exploration area and determining the plane coordinates of the ground base stations; acquiring horizontal electric field time series collected by each base station; synchronously collecting magnetic field time series and space coordinates of an air measurement point through an air mobile platform; projecting the space coordinates to the ground to obtain a projection point, and determining a dynamic weight coefficient according to the horizontal distance from the projection point to each base station; generating a virtual ground electric field by weighting and fusing the horizontal electric field and the corresponding weight coefficient; pairing the virtual ground electric field and the magnetic field time series, converting to the frequency domain, solving the tensor impedance and determining the electrical parameter of the target exploration area. The application is used for solving the problem of spatial position mismatch between the air mobile measurement point and the ground fixed base station under the existing single base station architecture.

Description

Technical Field

[0001] This invention relates to the field of geophysical exploration technology, and more specifically, to an airborne magnetotelluric data processing method, apparatus, equipment, and storage medium. Background Technology

[0002] Airborne magnetotellurics is a highly efficient non-contact exploration method in the field of geophysical exploration for deep resource exploration and geological structure interpretation. Its core principle is to use the natural magnetotelluric field as a field source, measure the electric and magnetic field components of the electromagnetic field, solve the tensor impedance of the underground medium according to Maxwell's equations, and then invert the underground electrical structure. It is widely used in mineral resource exploration, oil and gas reservoir exploration, geological hazard investigation, and deep geological research.

[0003] Currently, various technical solutions have been developed in the industry for the research and engineering application of airborne magnetotelluric surveying. The core of these solutions is based on a ground-based base station collecting electric field components and an airborne platform collecting magnetic field components, with subsequent technical optimizations and improvements. While these solutions have improved and upgraded airborne magnetotelluric surveying technology from different dimensions, achieving enhancements in data acquisition signal-to-noise ratio, inversion accuracy, and power spectrum estimation reliability, none have broken through the traditional single-ground-base station acquisition architecture and have failed to solve the technical problem of spatial mismatch between airborne mobile measuring points and ground-based fixed base stations. Summary of the Invention

[0004] The purpose of this invention is to provide an airborne magnetotelluric data processing method, apparatus, device, and storage medium to correct for positional errors in aerial mobile measurements, adapt to the non-uniformity of electromagnetic fields, and improve operational efficiency and adaptability to complex terrain. To achieve the above objectives, the technical solution adopted by this invention is as follows:

[0005] In a first aspect, this application provides a method for processing airborne magnetotelluric data, including:

[0006] Deploy at least two ground base stations within the target exploration area and determine the planar coordinates of each ground base station;

[0007] The horizontal electric field time series of the magnetotelluric field collected by each ground base station is obtained, and the horizontal electric field time series is collected according to the preset sampling time.

[0008] The magnetic field time series of the aerial measuring points is collected synchronously through an aerial mobile platform, and the spatial coordinates of the aerial measuring points at each sampling time are collected.

[0009] The spatial coordinates of each sampling time are projected onto the ground to obtain projection points. Based on the horizontal distance from the projection points to each ground base station, the dynamic weight coefficients of each ground base station are determined. The dynamic weight coefficients are inversely correlated with the horizontal distance.

[0010] The horizontal electric field in the horizontal electric field time series is weighted and fused with the dynamic weight coefficients at the corresponding sampling time to generate the virtual ground electric field at the projection point;

[0011] The virtual ground electric field at each sampling moment is paired with the magnetic field time series of the aerial measuring point and then converted to the frequency domain. The tensor impedance is solved in the frequency domain to determine the electrical parameters of the target exploration area.

[0012] Secondly, this application also provides an airborne magnetotelluric data processing device, comprising:

[0013] The first acquisition module is used to acquire the planar coordinates of at least two ground base stations deployed in the target exploration area, and the horizontal electric field time series of the magnetotelluric field collected by each ground base station, wherein the horizontal electric field time series is acquired according to a preset sampling time.

[0014] The second acquisition module is used to acquire the magnetic field time series of the aerial measuring points synchronously collected by the aerial mobile platform, as well as the spatial coordinates of the aerial measuring points at each sampling time.

[0015] The weight determination module is used to project the spatial coordinates of each sampling time onto the ground to obtain a projection point, and determine the dynamic weight coefficient of each ground base station based on the horizontal distance from the projection point to each ground base station, wherein the dynamic weight coefficient is inversely correlated with the horizontal distance;

[0016] The virtual electric field synthesis module is used to weight and fuse the horizontal electric field in the horizontal electric field time series with the dynamic weight coefficients at the corresponding sampling time to generate the virtual ground electric field at the projection point.

[0017] The impedance calculation module is used to pair the virtual ground electric field at each sampling moment with the magnetic field time series of the aerial measuring point and convert it to the frequency domain. In the frequency domain, the tensor impedance is solved to determine the electrical parameters of the target exploration area.

[0018] Thirdly, this application also provides an airborne magnetotelluric data processing device, comprising:

[0019] Memory, used to store computer programs;

[0020] A processor is used to implement the steps of the airborne magnetotelluric data processing method when executing the computer program.

[0021] Fourthly, this application also provides a readable storage medium storing a computer program, which, when executed by a processor, implements the steps of the above-described airborne magnetotelluric data processing method.

[0022] The beneficial effects of this invention are as follows:

[0023] 1. This invention allows an aerial platform to fly continuously over a wide area by deploying multiple ground base stations. It can operate as long as there is signal coverage from any two base stations, without having to hover around a single base station. This significantly increases the data collection area per flight and reduces the frequency of base station relocation and field operation costs.

[0024] 2. This invention effectively corrects the position mismatch error caused by the spatial separation of aerial measuring points and base stations by deploying multiple base stations within the target exploration area and assigning dynamic weight coefficients to the multiple base stations based on the real-time spatial location of the aerial measuring points, thereby synthesizing a virtual ground electric field.

[0025] Other features and advantages of the invention will be set forth in the following description, and will be apparent in part from the description, or may be learned by practicing embodiments of the invention. The objects and other advantages of the invention may be realized and obtained by means of the structures particularly pointed out in the written description, claims, and drawings. Attached Figure Description

[0026] To more clearly illustrate the technical solutions of the embodiments of the present invention, the accompanying drawings used in the embodiments will be briefly introduced below. It should be understood that the following drawings only show some embodiments of the present invention and should not be regarded as a limitation on the scope. For those skilled in the art, other related drawings can be obtained based on these drawings without creative effort.

[0027] Figure 1 This is a schematic diagram of the airborne magnetotelluric data processing method described in an embodiment of the present invention;

[0028] Figure 2 This is a schematic diagram illustrating the positional relationship between the ground base station and the airborne mobile platform in an embodiment of the present invention;

[0029] Figure 3 This is a schematic diagram of the airborne magnetotelluric data processing device described in an embodiment of the present invention.

[0030] Figure 4 This is a schematic diagram of the airborne magnetotelluric data processing equipment described in an embodiment of the present invention.

[0031] Marked in the image:

[0032] 800. Airborne magnetotelluric data processing equipment; 801. Processor; 802. Memory; 803. Multimedia components; 804. I / O interface; 805. Communication components. Detailed Implementation

[0033] To make the objectives, technical solutions, and advantages of the embodiments of the present invention clearer, the technical solutions of the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some, not all, of the embodiments of the present invention. The components of the embodiments of the present invention described and shown in the accompanying drawings can generally be arranged and designed in various different configurations. Therefore, the following detailed description of the embodiments of the present invention provided in the accompanying drawings is not intended to limit the scope of the claimed invention, but merely to illustrate selected embodiments of the invention. All other embodiments obtained by those skilled in the art based on the embodiments of the present invention without inventive effort are within the scope of protection of the present invention.

[0034] It should be noted that similar reference numerals and letters in the following figures indicate similar items; therefore, once an item is defined in one figure, it does not need to be further defined and explained in subsequent figures. Furthermore, in the description of this invention, terms such as "first," "second," etc., are used only to distinguish descriptions and should not be construed as indicating or implying relative importance.

[0035] Example 1:

[0036] Therefore, this embodiment provides an airborne magnetotelluric data processing method that is suitable for exploration tasks under complex terrain conditions.

[0037] See Figure 1 , Figure 2 The figure shows that this method includes:

[0038] S1. Deploy at least two ground base stations within the target exploration area and determine the planar coordinates of each ground base station;

[0039] Specifically, step S1 includes:

[0040] S11. Based on the exploration depth, geological complexity, and survey area of ​​the target exploration area, determine the deployment locations of at least two ground base stations, and ensure that the connection between at least two ground base stations covers the target exploration area;

[0041] In this embodiment, the deployment of two ground base stations (base station A and base station B) is used as an example for illustration. Specifically, the deployment distance between the two base stations is set to several hundred meters to several kilometers according to actual needs, and the connection between the two base stations is guaranteed to cover the main survey line area within the target exploration area.

[0042] S12. Determine the planar coordinates of each ground base station, where the planar coordinates of base station A are: The planar coordinates of base station B are

[0043] S13. Configure each ground base station with an electric field sensor, a data acquisition instrument, a clock synchronization module, and a power supply module; wherein, the electric field sensor is used to sense the horizontal electric field component of the earth's electromagnetic field, the data acquisition instrument is used to convert the analog signals output by the sensor into digital signals and record the time series, the clock synchronization module is used to achieve time synchronization, and the power supply module provides the base station with a continuous and stable operating power supply.

[0044] Based on the above embodiments, this method further includes:

[0045] S2. Obtain the horizontal electric field time series of the magnetotelluric field collected by each ground base station, wherein the horizontal electric field time series is collected according to a preset sampling time;

[0046] Specifically, step S2 includes:

[0047] S21. Obtain the first and second directional components of the horizontal electric field collected by the electric field sensor of each ground base station;

[0048] For base station A, the first direction component of the recorded horizontal electric field is acquired. Second direction component For base station B, the first direction component of the recorded horizontal electric field is acquired. Second direction component ,in, Indicates the sampling time.

[0049] S22. Arrange the first directional component and the second directional component into a time series according to a preset sampling time to obtain the horizontal electric field time series; specifically, the horizontal electric field time series of base station A is expressed as follows: and , The time series of the horizontal electric field of base station B is represented as { and .

[0050] Based on the above embodiments, this method further includes:

[0051] S3. Synchronously collect the magnetic field time series of the aerial measuring points through an aerial mobile platform, and collect the spatial coordinates of the aerial measuring points at each sampling time;

[0052] In this embodiment, the aerial mobile platform can be an unmanned aerial vehicle (UAV) or a manned aircraft. Specifically, the aerial mobile platform integrates a three-axis fluxgate magnetometer, a GNSS positioning device, a data acquisition unit, and a time synchronization module. The three-axis magnetometer is used to sense the three components of the magnetotelluric field at the aerial measurement point, and the GNSS positioning device is used to measure the spatial position of the aerial mobile platform in real time.

[0053] At each sampling time The three-axis fluxgate magnetometer collects the three components of the magnetic field at the aerial measuring point, which are respectively ,in, Represents the magnetic field components in the east-west direction. This represents the magnetic field component in the north-south direction. This represents the vertical magnetic field component. Simultaneously, the positioning device records the three-dimensional spatial coordinates of the aerial mobile platform. ,in and These are the planar coordinate components, This refers to altitude.

[0054] Based on the above embodiments, this method further includes:

[0055] S4. Project the spatial coordinates of each sampling time onto the ground to obtain projection points. Determine the dynamic weight coefficients of each ground base station based on the horizontal distance from the projection points to each ground base station. The dynamic weight coefficients are inversely correlated with the horizontal distance.

[0056] Specifically, step S4 includes:

[0057] S41. Project the spatial coordinates at each sampling time point vertically onto the ground to obtain a projection point, the coordinates of which are: This indicates the position where the aerial measuring point intersects the ground directly below it.

[0058] S42. Calculate the horizontal distance between the projection point and each ground base station respectively;

[0059] In this embodiment, the projection points are calculated respectively. Horizontal distances to ground base station A and base station B:

[0060] ;

[0061] In the formula, Represents the projection point Distance to ground base station A Represents the projection point Distance to ground base station B.

[0062] S43. Correct each horizontal distance using a preset smoothing factor to obtain the corrected distance for each ground base station;

[0063] To prevent computational overflow caused by the distance between the projected point and the base station being zero when they coincide, this embodiment introduces a positive decimal smoothing factor. The horizontal distance is corrected by the positive decimal smoothing factor. The distance is determined based on the measurement accuracy and the distance between base stations, and is usually 0.01~0.05km.

[0064] The corrected distances for base stations A and B are respectively )and ).

[0065] Specifically, step S4 also includes:

[0066] S44. For the corrected distance of each ground base station, determine the weight base value of each ground base station according to the inverse relationship between the smaller the corrected distance and the weight base value;

[0067] Based on the principle of distance-inverse weighting, the smaller the correction distance of a base station, the stronger the correlation between its electric field data and the electric field at the airborne measurement point; therefore, it should be assigned a larger weight base value. Thus, in this embodiment, the reciprocal of the correction distance of each ground base station is used as the weight base value for the base station. For base station A and base station B, their weight base values ​​are respectively... and .

[0068] S45. Sum the weight base values ​​of all ground base stations to obtain the total weight base value;

[0069] S46. Calculate the ratio of the weight base value of each ground base station to the sum of the weight base values ​​to obtain the dynamic weight coefficient of the corresponding ground base station:

[0070] ;

[0071] In the formula, For base station A, the dynamic weighting coefficient is... This represents the dynamic weighting coefficient of base station B.

[0072] Furthermore, the dynamic weighting coefficients of all ground base stations satisfy the normalization condition:

[0073] ;

[0074] In this embodiment, each sampling time point corresponds to a dynamic weighting coefficient of a ground base station.

[0075] Based on the above embodiments, this method further includes:

[0076] S5. The horizontal electric field in the horizontal electric field time series is weighted and fused with the dynamic weight coefficients at the corresponding sampling time to generate the virtual ground electric field at the projection point;

[0077] Specifically, step S5 includes:

[0078] S51. Weight the first directional component of each ground base station with the corresponding dynamic weight coefficient to obtain the first directional contribution value;

[0079] Specifically, for each sampling time The first-direction contribution value of base station A is The first-direction contribution value of base station B is ;

[0080] S52. Accumulate the first-direction contribution values ​​of all ground base stations to obtain the first-direction component of the virtual ground electric field. :

[0081] ;

[0082] S53. Weight the second direction component of each ground base station with the corresponding dynamic weight coefficient to obtain the second direction contribution value;

[0083] Specifically, for each sampling time The second-direction contribution value of base station A is The second-direction contribution value of base station B is ;

[0084] S54. Accumulate the second-direction contribution values ​​of all ground base stations to obtain the second-direction component of the virtual ground electric field. :

[0085] ;

[0086] In this embodiment, at each sampling time All were synthesized to obtain projection points The virtual ground electric field at that location has two orthogonal components, namely: and The virtual ground electric field replaces the traditional method of directly using the electric field data of a single base station. It reconstructs the real ground electric field response directly below the aerial measuring point from a physical perspective, correcting the position mismatch error caused by the aerial measuring point deviating from the single base station.

[0087] Based on the above embodiments, this method further includes:

[0088] S6. Pair the virtual ground electric field at each sampling time with the magnetic field time series of the aerial measuring point and convert it to the frequency domain. Solve the tensor impedance in the frequency domain to determine the electrical parameters of the target exploration area.

[0089] Specifically, step S6 includes:

[0090] S61. Pair the virtual ground electric field at each sampling moment with the magnetic field time series of the air measurement point according to the sampling point to obtain time-domain paired data;

[0091] Specifically, the virtual ground electric field components and With the three components of the magnetic field Pair the data one-to-one according to the same time label to obtain time-domain paired data. .

[0092] S62. Perform a Fourier transform on the time-domain paired data to obtain the electric field component and magnetic field component at each frequency point;

[0093] Each component in the time-domain paired dataset is transformed from the time domain to the frequency domain using Fourier transform. Specifically, the frequency domain component of the virtual ground electric field is... and The frequency domain component of the magnetic field at the air measurement point is .

[0094] S63. At each frequency point, establish a set of equations about the tensor impedance components based on the linear relationship between the electric and magnetic fields in the magnetotelluric field;

[0095] ;

[0096] In the formula, , , and The four components of the tensor impedance to be solved;

[0097] In this embodiment, a system of linear equations is established for the four unknown complex impedance components at each frequency point.

[0098] S64. Solve the system of equations by least squares method or robust estimation method to obtain each component of tensor impedance, and calculate the electrical parameters of the target exploration area based on the components of tensor impedance.

[0099] In this embodiment, the least squares method or robust estimation method is used to solve the system of equations to obtain the tensor impedance at each frequency point. Furthermore, the apparent resistivity and phase of the target exploration area are calculated based on the tensor impedance components:

[0100] ;

[0101] In the formula, Indicates apparent resistivity. Indicates phase, Indicates the imaginary part. Indicates the real part, It represents the vacuum permeability.

[0102] The calculated apparent resistivity and phase at each frequency point are used to construct apparent resistivity-frequency curves and phase-frequency curves, which reflect the electrical stratification characteristics of the underground medium from shallow to deep.

[0103] Example 2:

[0104] like Figure 3 As shown, this embodiment provides an airborne magnetotelluric data processing device, the device comprising:

[0105] The ground base station data acquisition module is used to deploy at least two ground base stations within the target exploration area and determine the planar coordinates of each ground base station.

[0106] The horizontal electric field acquisition module is used to acquire the horizontal electric field time series of the magnetotelluric field collected by each ground base station. The horizontal electric field time series is acquired according to a preset sampling time.

[0107] The aerial data acquisition module is used to synchronously acquire the magnetic field time series of aerial measuring points through an aerial mobile platform, and to acquire the spatial coordinates of the aerial measuring points at each sampling time.

[0108] The weight determination module is used to project the spatial coordinates of each sampling time onto the ground to obtain a projection point, and determine the dynamic weight coefficient of each ground base station based on the horizontal distance from the projection point to each ground base station, wherein the dynamic weight coefficient is inversely correlated with the horizontal distance;

[0109] The virtual electric field synthesis module is used to weight and fuse the horizontal electric field in the horizontal electric field time series with the dynamic weight coefficients at the corresponding sampling time to generate the virtual ground electric field at the projection point.

[0110] The impedance calculation module is used to pair the virtual ground electric field at each sampling moment with the magnetic field time series of the aerial measuring point and convert it to the frequency domain. In the frequency domain, the tensor impedance is solved to determine the electrical parameters of the target exploration area.

[0111] Based on the above embodiments, the ground base station data acquisition module includes:

[0112] The deployment location determination unit is used to determine the deployment locations of at least two ground base stations based on the exploration depth, geological complexity, and survey area of ​​the target exploration area, and to ensure that the line connecting at least two ground base stations covers the target exploration area.

[0113] Planar coordinate determination unit, used to determine the planar coordinates of each ground base station;

[0114] The base station configuration unit is used to configure an electric field sensor, a data acquisition instrument, a clock synchronization module, and a power supply module for each ground base station.

[0115] Based on the above embodiments, the horizontal electric field acquisition module includes:

[0116] The component acquisition unit is used to acquire the first directional component and the second directional component of the horizontal electric field collected by the electric field sensor of each ground base station.

[0117] The time series arrangement unit is used to arrange the first directional component and the second directional component into a time series according to a preset sampling time to obtain the horizontal electric field time series.

[0118] Based on the above embodiments, the weight determination module includes:

[0119] A vertical projection unit is used to vertically project the spatial coordinates at each sampling time onto the ground to obtain the projection point;

[0120] A horizontal distance calculation unit is used to calculate the horizontal distance between the projection point and each ground base station respectively;

[0121] The distance correction unit is used to correct each horizontal distance using a preset smoothing factor to obtain the corrected distance for each ground base station.

[0122] Based on the above embodiments, the weight determination module further includes:

[0123] The weight base value determination unit is used to determine the weight base value of each ground base station based on the principle that the smaller the correction distance, the more inversely proportional it is to the weight base value.

[0124] The weighted base value summation unit is used to sum the weighted base values ​​of all ground base stations to obtain the total weighted base value;

[0125] The dynamic weight coefficient calculation unit is used to calculate the ratio of the weight base value of each ground base station to the sum of the weight base values, so as to obtain the dynamic weight coefficient of the corresponding ground base station.

[0126] Based on the above embodiments, the virtual electric field synthesis module includes:

[0127] The first directional weighting unit is used to weight the first directional component of each ground base station with the corresponding dynamic weight coefficient to obtain the first directional contribution value.

[0128] The first directional accumulation unit is used to accumulate the first directional contribution values ​​of all ground base stations as the first directional component of the virtual ground electric field.

[0129] The second directional weighting unit is used to weight the second directional component of each ground base station with the corresponding dynamic weight coefficient to obtain the second directional contribution value.

[0130] The second direction accumulation unit is used to accumulate the second direction contribution values ​​of all ground base stations as the second direction component of the virtual ground electric field.

[0131] Based on the above embodiments, the impedance calculation module includes:

[0132] The time-domain pairing unit is used to pair the virtual ground electric field at each sampling moment with the magnetic field time series of the air measurement point according to the sampling point to obtain time-domain pairing data;

[0133] The Fourier transform unit is used to perform Fourier transform on the time-domain paired data to obtain the electric field component and magnetic field component at each frequency point.

[0134] The equation-setting unit is used to establish a set of equations about the tensor impedance components at each frequency point, based on the linear relationship between the electric and magnetic fields in the geomagnetic field.

[0135] The impedance solving unit is used to solve the system of equations by least squares method or robust estimation method to obtain each component of tensor impedance, and to calculate the electrical parameters of the target exploration area based on the components of tensor impedance.

[0136] It should be noted that the specific manner in which each module performs its operation in the apparatus described in the above embodiments has been described in detail in the embodiments of the method, and will not be elaborated here.

[0137] Example 3:

[0138] Corresponding to the above method embodiments, this embodiment also provides an airborne magnetotelluric data processing device. The airborne magnetotelluric data processing device described below and the airborne magnetotelluric data processing method described above can be referred to in correspondence.

[0139] Figure 4 This is a block diagram illustrating an airborne magnetotelluric data processing device 800 according to an exemplary embodiment. Figure 4 As shown, the airborne magnetotelluric data processing device 800 may include a processor 801 and a memory 802. The airborne magnetotelluric data processing device 800 may also include one or more of a multimedia component 803, an I / O interface 804, and a communication component 805.

[0140] The processor 801 controls the overall operation of the airborne magnetotelluric data processing device 800 to complete all or part of the steps in the aforementioned airborne magnetotelluric data processing method. The memory 802 stores various types of data to support the operation of the airborne magnetotelluric data processing device 800. This data may include, for example, instructions for any application or method operating on the airborne magnetotelluric data processing device 800, as well as application-related data such as contact data, sent and received messages, images, audio, video, etc. The memory 802 can be implemented by any type of volatile or non-volatile storage device or a combination thereof, such as Static Random Access Memory (SRAM), Electrically Erasable Programmable Read-Only Memory (EEPROM), Erasable Programmable Read-Only Memory (EPROM), Programmable Read-Only Memory (PROM), Read-Only Memory (ROM), magnetic storage, flash memory, magnetic disk, or optical disk. The multimedia component 803 may include a screen and an audio component. The screen may be, for example, a touchscreen, and the audio component is used to output and / or input audio signals. For example, the audio component may include a microphone for receiving external audio signals. The received audio signals may be further stored in the memory 802 or transmitted via the communication component 805. The audio component also includes at least one speaker for outputting audio signals. I / O interface 804 provides an interface between processor 801 and other interface modules, such as keyboards, mice, and buttons. These buttons can be virtual or physical. Communication component 805 is used for wired or wireless communication between the airborne magnetotelluric data processing device 800 and other devices. Wireless communication includes, for example, Wi-Fi, Bluetooth, Near Field Communication (NFC), 2G, 3G, or 4G, or a combination thereof. Therefore, the corresponding communication component 805 may include a Wi-Fi module, a Bluetooth module, and an NFC module.

[0141] In an exemplary embodiment, the airborne magnetotelluric data processing device 800 may be implemented by one or more application-specific integrated circuits (ASICs), digital signal processors (DSPs), digital signal processing devices (DSPDs), programmable logic devices (PLDs), field-programmable gate arrays (FPGAs), controllers, microcontrollers, microprocessors, or other electronic components to perform the airborne magnetotelluric data processing method described above.

[0142] In another exemplary embodiment, a computer-readable storage medium including program instructions is also provided, which, when executed by a processor, implement the steps of the airborne magnetotelluric data processing method described above. For example, the computer-readable storage medium may be the memory 802 including the program instructions described above, which may be executed by the processor 801 of the airborne magnetotelluric data processing device 800 to complete the airborne magnetotelluric data processing method described above.

[0143] Example 4:

[0144] Corresponding to the above method embodiments, this embodiment also provides a readable storage medium. The readable storage medium described below can be referred to in conjunction with the airborne magnetotelluric data processing method described above.

[0145] A readable storage medium storing a computer program, which, when executed by a processor, implements the steps of the airborne magnetotelluric data processing method described in the above method embodiments.

[0146] Specifically, the readable storage medium can be a USB flash drive, a portable hard drive, a read-only memory (ROM), a random access memory (RAM), a magnetic disk, or an optical disk, or any other readable storage medium capable of storing program code.

[0147] The above description is merely a preferred embodiment of the present invention and is not intended to limit the invention. Various modifications and variations can be made to the present invention by those skilled in the art. Any modifications, equivalent substitutions, improvements, etc., made within the spirit and principles of the present invention should be included within the scope of protection of the present invention.

[0148] The above description is merely a specific embodiment of the present invention, but the scope of protection of the present invention is not limited thereto. Any variations or substitutions that can be easily conceived by those skilled in the art within the technical scope disclosed in the present invention should be included within the scope of protection of the present invention. Therefore, the scope of protection of the present invention should be determined by the scope of the claims.

Claims

1. A method for processing airborne magnetotelluric data, characterized in that, include: Deploy at least two ground base stations within the target exploration area and determine the planar coordinates of each ground base station; The horizontal electric field time series of the magnetotelluric field collected by each ground base station is obtained, and the horizontal electric field time series is collected according to the preset sampling time. The magnetic field time series of the aerial measuring points is collected synchronously through an aerial mobile platform, and the spatial coordinates of the aerial measuring points at each sampling time are collected. The spatial coordinates of each sampling time are projected onto the ground to obtain projection points. Based on the horizontal distance from the projection points to each ground base station, the dynamic weight coefficients of each ground base station are determined. The dynamic weight coefficients are inversely correlated with the horizontal distance. The horizontal electric field in the horizontal electric field time series is weighted and fused with the dynamic weight coefficients at the corresponding sampling time to generate the virtual ground electric field at the projection point; The virtual ground electric field at each sampling moment is paired with the magnetic field time series of the aerial measuring point and then converted to the frequency domain. The tensor impedance is solved in the frequency domain to determine the electrical parameters of the target exploration area.

2. The airborne magnetotelluric data processing method according to claim 1, characterized in that, Deploy at least two ground base stations within the target exploration area and determine the planar coordinates of each ground base station, including: Based on the exploration depth, geological complexity, and survey area of ​​the target exploration area, determine the deployment locations of at least two ground base stations, and ensure that the line connecting at least two ground base stations covers the target exploration area. Determine the planar coordinates of each ground base station; Each ground base station is equipped with an electric field sensor, a data acquisition instrument, a clock synchronization module, and a power supply module.

3. The airborne magnetotelluric data processing method according to claim 2, characterized in that, Obtain the horizontal electric field time series of the magnetotelluric field collected by each ground base station, including: Acquire the first and second directional components of the horizontal electric field collected by the electric field sensor of each ground base station; The first directional component and the second directional component are arranged into a time series according to the preset sampling time to obtain the horizontal electric field time series.

4. The airborne magnetotelluric data processing method according to claim 1, characterized in that, The spatial coordinates of each sampling time are projected onto the ground to obtain projection points. Based on the horizontal distance from the projection points to each ground base station, the following is included: The spatial coordinates at each sampling time are vertically projected onto the ground to obtain the projection point; Calculate the horizontal distance between the projection point and each ground base station; The corrected distance for each ground base station is obtained by correcting each horizontal distance using a preset smoothing factor.

5. The airborne magnetotelluric data processing method according to claim 4, characterized in that, Determine the dynamic weighting coefficients for each ground base station, including: For each ground base station, the weight base value is determined according to the inverse relationship between the smaller the correction distance and the weight base value. The total weighted base values ​​are obtained by summing the weighted base values ​​of all ground base stations. The dynamic weight coefficient of the corresponding ground base station is obtained by calculating the ratio of the weight base value of each ground base station to the sum of the weight base values.

6. The airborne magnetotelluric data processing method according to claim 1, characterized in that, The horizontal electric field in the time series of the horizontal electric field is weighted and fused with the dynamic weight coefficients at the corresponding sampling time to generate the virtual ground electric field at the projection point, including: The first directional component of each ground base station is weighted with the corresponding dynamic weight coefficient to obtain the first directional contribution value. The first-direction contribution values ​​of all ground base stations are summed up to form the first-direction component of the virtual ground electric field. The second-direction component of each ground base station is weighted with the corresponding dynamic weight coefficient to obtain the second-direction contribution value. The second-direction contribution values ​​of all ground base stations are summed to form the second-direction component of the virtual ground electric field.

7. The airborne magnetotelluric data processing method according to claim 1, characterized in that, The virtual ground electric field at each sampling time is paired with the magnetic field time series of the aerial measuring points and then converted to the frequency domain. The tensor impedance is then solved in the frequency domain to determine the electrical parameters of the target exploration area, including: The virtual ground electric field at each sampling moment is paired with the magnetic field time series of the air measurement point according to the one-to-one correspondence of the sampling points to obtain time-domain paired data; Perform a Fourier transform on the time-domain paired data to obtain the electric field component and magnetic field component at each frequency point; At each frequency point, a set of equations about the tensor impedance components is established based on the linear relationship between the electric and magnetic fields in the magnetotelluric field. The equations are solved by least squares method or robust estimation method to obtain the components of tensor impedance, and the electrical parameters of the target exploration area are calculated based on the components of tensor impedance.

8. An airborne magnetotelluric data processing device, characterized in that, include: The ground base station data acquisition module is used to deploy at least two ground base stations within the target exploration area and determine the planar coordinates of each ground base station. The horizontal electric field acquisition module is used to acquire the horizontal electric field time series of the magnetotelluric field collected by each ground base station. The horizontal electric field time series is acquired according to the preset sampling time. The aerial data acquisition module is used to synchronously acquire the magnetic field time series of aerial measuring points through an aerial mobile platform, and to acquire the spatial coordinates of the aerial measuring points at each sampling time. The weight determination module is used to project the spatial coordinates of each sampling time onto the ground to obtain a projection point, and determine the dynamic weight coefficient of each ground base station based on the horizontal distance from the projection point to each ground base station, wherein the dynamic weight coefficient is inversely correlated with the horizontal distance; The virtual electric field synthesis module is used to weight and fuse the horizontal electric field in the horizontal electric field time series with the dynamic weight coefficients at the corresponding sampling time to generate the virtual ground electric field at the projection point. The impedance calculation module is used to pair the virtual ground electric field at each sampling time with the magnetic field time series of the aerial measuring point and convert it to the frequency domain. In the frequency domain, the tensor impedance is solved to determine the electrical parameters of the target exploration area.

9. An airborne magnetotelluric data processing device, characterized in that, include: Memory, used to store computer programs; A processor, configured to execute the computer program to implement the steps of the airborne magnetotelluric data processing method as described in any one of claims 1 to 7.

10. A readable storage medium, characterized in that: The readable storage medium stores a computer program that, when executed by a processor, implements the steps of the airborne magnetotelluric data processing method as described in any one of claims 1 to 7.

Images