Low-latitude magnetic anomaly reversal method and device, electronic equipment and storage medium

Abstract

The application provides a low-latitude magnetic anomaly reduction-to-pole method and device, electronic equipment and storage medium, comprising: obtaining frequency domain magnetic anomaly data through Fourier transform on acquired space domain magnetic anomaly data; separating initial amplitude spectrum and initial phase spectrum by decomposing frequency domain initial pole factor; modifying the initial amplitude spectrum through a suppression modification algorithm, and setting an adaptive modification domain; recombining the modified amplitude spectrum and the separated initial phase spectrum to construct a recombined pole factor, and performing pole reduction operation on the frequency domain magnetic anomaly data and the recombined pole factor to obtain pole-reduced frequency domain magnetic anomaly data; and obtaining a pole reduction result by inverse Fourier transform of the pole-reduced frequency domain magnetic anomaly data to the space domain. The application has accurate and reliable reduction-to-pole results in low-latitude regions, higher precision, and effectively eliminates the stretching phenomenon of magnetic anomalies along the direction perpendicular to the magnetic declination angle.

Description

Technical Field

[0001] This application belongs to the field of geophysical exploration technology, and in particular relates to a low-latitude magnetic anomaly polarization method, device, electronic equipment and storage medium. Background Technology

[0002] When the background magnetic field forms an angle with the horizontal direction, the magnetic anomaly will undergo tilt magnetization, complicating the distribution characteristics of the magnetic anomaly and making the relationship between the magnetic anomaly field and its source less intuitive. This greatly increases the difficulty of accurately interpreting the magnetic anomaly field, necessitating polarization treatment.

[0003] In low-latitude regions, conventional pole-setting of magnetic anomalies introduces significant errors into the results, leading to stretching of contour lines perpendicular to magnetic declination and banded interference, thus reducing the accuracy of the pole-setting results. Obtaining high-precision magnetic anomaly pole-setting results efficiently is a key challenge in improving the accuracy and utilization of magnetic data interpretation in low-latitude regions. Summary of the Invention

[0004] In view of this, this application aims to propose a method, device, electronic device and storage medium for polarization of magnetic anomalies in low latitudes, so as to solve the problems of amplification and stretching of high-frequency interference along the direction perpendicular to the magnetic declination in the polarization process of magnetic anomaly data in low latitude regions.

[0005] To achieve the above objectives, the technical solution of this application is implemented as follows:

[0006] In a first aspect, this application provides a low-latitude magnetic anomaly polarization method, including:

[0007] The acquired spatial domain magnetic anomaly data is transformed into frequency domain magnetic anomaly data using Fourier transform.

[0008] The initial amplitude spectrum and the initial phase spectrum are obtained by decomposing the frequency domain initialization polar factor, and the initial amplitude spectrum is modified by the suppression modification algorithm, and an adaptive modification domain is set.

[0009] The modified amplitude spectrum is recombined with the separated initial phase spectrum to construct a recombined polarization factor. The frequency domain magnetic anomaly data and the recombined polarization factor are then subjected to polarization processing to obtain polarized frequency domain magnetic anomaly data.

[0010] The frequency domain magnetic anomaly data after polarization is transformed to the spatial domain through inverse Fourier transform to obtain the polarization result.

[0011] Secondly, based on the same inventive concept, this application also provides a low-latitude magnetic anomaly polarization device, comprising:

[0012] The Fourier transform module is configured to obtain frequency domain magnetic anomaly data by Fourier transforming the acquired spatial domain magnetic anomaly data.

[0013] The amplitude spectrum modification module is configured to separate the initial amplitude spectrum and the initial phase spectrum by decomposing the frequency domain initial polar factor, and to modify the initial amplitude spectrum by a suppression modification algorithm, and to set an adaptive modification domain.

[0014] The polarization factor recombination module is configured to recombine the modified amplitude spectrum with the separated initial phase spectrum to construct a recombined polarization factor, and to perform polarization processing on the frequency domain magnetic anomaly data and the recombined polarization factor to obtain polarized frequency domain magnetic anomaly data.

[0015] The inverse Fourier transform module is configured to transform the frequency domain magnetic anomaly data after polarization to the spatial domain through inverse Fourier transform to obtain the polarization result.

[0016] Thirdly, based on the same inventive concept, this application also provides an electronic device, including a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor executes the program to implement the method described in the first aspect.

[0017] Fourthly, based on the same inventive concept, this application also provides a non-transitory computer-readable storage medium, wherein the non-transitory computer-readable storage medium stores computer instructions for causing the computer to perform the method as described in the first aspect.

[0018] Compared with the prior art, the low-latitude magnetic anomaly polarization method, apparatus, electronic device, and storage medium described in this application have the following advantages:

[0019] The low-latitude magnetic anomaly polarization method described in this application is based on the characteristics of the polarization factor, specifically modifies the polarization factor, and introduces an adaptive modification domain to suppress the amplification effect of conventional polarization methods, thereby achieving the purpose of improving the low-latitude polarization effect and increasing the accuracy of the results. Attached Figure Description

[0020] The accompanying drawings, which form part of this application, are used to provide a further understanding of this application. The illustrative embodiments and descriptions of this application are used to explain this application and do not constitute an undue limitation of this application. In the drawings:

[0021] Figure 1 This is a flowchart of a low-latitude magnetic anomaly polarization method according to an embodiment of this application;

[0022] Figure 2 This is a schematic diagram of the conventional amplitude spectrum at different magnetic tilt angles as described in the embodiments of this application;

[0023] Figure 3 This is a schematic diagram of the new polarization factor amplitude spectrum and adaptive modification domain under different magnetic tilt angles as described in the embodiments of this application;

[0024] Figure 4 This is a schematic diagram comparing the single-model forward magnetic field characteristics and polarization results described in the embodiments of this application;

[0025] Figure 5 This is a schematic diagram comparing the multi-model forward magnetic field characteristics and polarization results described in the embodiments of this application;

[0026] Figure 6 This is a schematic diagram of the measured magnetic anomaly data described in the embodiments of this application;

[0027] Figure 7 This is a schematic diagram of the measured magnetic anomaly data polarization results described in the embodiments of this application;

[0028] Figure 8 This is a schematic diagram of a low-latitude magnetic anomaly polarization device according to an embodiment of this application;

[0029] Figure 9 This is a schematic diagram of the hardware structure of the electronic device described in an embodiment of this application. Detailed Implementation

[0030] To make the objectives, technical solutions, and advantages of this application clearer, the following detailed description is provided in conjunction with specific embodiments and the accompanying drawings.

[0031] It should be noted that, unless otherwise defined, the technical or scientific terms used in the embodiments of this application should have the ordinary meaning understood by one of ordinary skill in the art to which this application pertains. The terms "first," "second," and similar terms used in the embodiments of this application do not indicate any order, quantity, or importance, but are merely used to distinguish different components. Terms such as "comprising" or "including" mean that the element or object preceding the word encompasses the elements or objects listed after the word and their equivalents, without excluding other elements or objects. Terms such as "connected" or "linked" are not limited to physical or mechanical connections, but can include electrical connections, whether direct or indirect. Terms such as "upper," "lower," "left," and "right" are only used to indicate relative positional relationships; when the absolute position of the described object changes, the relative positional relationship may also change accordingly.

[0032] As described in the background art, this application aims to solve the problems of amplification and stretching of high-frequency interference along the direction perpendicular to magnetic declination in the polarization processing of magnetic anomaly data in low-latitude regions (generally referring to regions or sea areas with latitude below 20°), thereby improving the accuracy of magnetic data interpretation. This application proposes a magnetic anomaly polarization method applicable to low-latitude regions. Through verification with single-model and multi-model forward modeling data and analysis of measured data, this method demonstrates accurate and reliable polarization results in low-latitude regions, with higher precision, and effectively eliminates the stretching phenomenon of magnetic anomalies along the direction perpendicular to magnetic declination.

[0033] The embodiments of this application are described in detail below with reference to the accompanying drawings.

[0034] Please see Figure 1 As shown, this embodiment provides a low-latitude magnetic anomaly polarization method, which specifically includes the following steps:

[0035] Step S101: Obtain frequency domain magnetic anomaly data by Fourier transforming the acquired spatial domain magnetic anomaly data.

[0036] Specifically, in this embodiment, two magnetic anomaly forward modeling models are constructed, including a single model and a multi-model model. The single model uses a cuboid anomaly with a length, width, and height of 20 m, 20 m, and 2 m respectively, and a center burial depth of 2 m. It calculates a 65×65 magnetic anomaly field, with a grid spacing of 1 m in both the X and Y directions. The magnetization model forward models the background magnetic field magnetic tilt angle. Magnetic declination The multi-model consists of multiple individual unit models, including three cuboid models and two spherical models. The dimensions of the three cuboids are (30 m, 20 m, 2 m), (20 m, 30 m, 2 m), and (20 m, 20 m, 2 m), respectively, with a center depth of 2 m. The two spherical models each have a radius of 10 m and a center depth of 15 m. The model forward models the background magnetic field tilt angle. Magnetic declination The grid size is 100×100.

[0037] Earth's magnetic field induction calculations were performed on the single and multiple models constructed above to obtain the model's magnetic induction field, which was used as the initial magnetic anomaly data. Gaussian white noise with a mean of 0 and a variance of 1 was added to the initial magnetic anomaly data of the single and multiple models, and regularization filtering was performed on it to simulate the magnetic anomaly data in actual work as closely as possible.

[0038] By applying the cosine attenuation method to expand the edges of the forward model data, it is easier to perform a fast Fourier transform to obtain the frequency domain magnetic anomaly data. The specific formula is as follows:

[0039]

[0040] In the formula: This is frequency domain magnetic anomaly data. This is a Fourier transform operation. This is spatial domain magnetic anomaly data.

[0041] Step S102: Decompose the frequency domain initial polar factor to separate the initial amplitude spectrum and the initial phase spectrum, and modify the initial amplitude spectrum by the suppression modification algorithm, and set the adaptive modification domain.

[0042] Specifically, in this embodiment, an amplitude spectrum suppression and modification algorithm is adopted to stabilize the amplitude spectrum. Based on the different magnetic inclination angles, an adaptive amplitude spectrum modification domain is introduced to suppress the amplitude spectrum to the optimal degree.

[0043] Furthermore, phase spectrum separation and amplitude spectrum separation are performed on the polarization factor, and the separation formulas are as follows:

[0044]

[0045] in, Represents the amplitude spectrum of the polarization factor (see Figure 2 ), , Represents the phase spectrum of the polarization factor. , and These are the initial magnetic anomaly tilt angle and the initial magnetic anomaly declination angle, respectively. The value represents the angle, ranging from 0 to 360°, and i represents the imaginary part.

[0046] The separated initial amplitude spectrum is specifically modified to suppress its unique amplification effect, and an adaptive modification domain is introduced to optimize the adverse effects of the modified polarization factor on the polarization results. The modification formula is as follows:

[0047]

[0048] In the formula, For adaptive transformation of the domain angle (see Figure 3 ), As a modifying factor, The amplitude spectrum after modification. This represents an angle, with a value ranging from 0 to 360°.

[0049] It should be noted that, from the perspective of the transformation domain The value is determined by the initial magnetic anomaly tilt angle. and initial magnetic anomaly declination The magnitude is determined jointly; at the same time, the amplitude spectrum modification effect From the initial magnetic anomaly magnetic tilt angle and initial magnetic anomaly declination , transformation domain perspective To be determined jointly.

[0050] In this step, by modifying the polarization factor, the stretching phenomenon of magnetic anomalies in low-latitude regions along the direction perpendicular to magnetic declination was effectively suppressed, improving the accuracy of the polarization results. Simultaneously, the adaptive selection method automatically adjusts the modification effect and modification domain based on the initial magnetic inclination and declination, reducing the adverse interference of human experience on the polarization results and improving the stability of the method.

[0051] Therefore, the polarization method proposed in this embodiment can achieve the effect of adaptively selecting the modification parameters based on the differences in initial data, greatly reducing the adverse effects of polarization processing differences caused by human experience and improving processing efficiency.

[0052] Step S103: Reconstruct the modified amplitude spectrum and the separated initial phase spectrum to construct the reconstructed polarization factor. Perform polarization processing on the frequency domain magnetic anomaly data and the reconstructed polarization factor to obtain the polarized frequency domain magnetic anomaly data.

[0053] Specifically, in this embodiment, the modified amplitude spectrum in step S102 is recombined with the separated initial phase spectrum to construct a new polarization factor (i.e., the recombined polarization factor). The specific formula is as follows:

[0054]

[0055] In the formula, Indicates the initial magnetic declination of the magnetic anomaly. This represents the amplitude spectrum after modification. Represents the original amplitude spectrum;

[0056] The frequency-domain magnetic anomaly data obtained in step S101 is multiplied with the recombined polarization factor to obtain the frequency-domain magnetic anomaly data after polarization (i.e., the magnetic anomaly spectrum). The specific formula is as follows:

[0057]

[0058] In the formula, As a recombination polar factor, This is frequency domain magnetic anomaly data. This is the frequency domain magnetic anomaly data after polarization.

[0059] Step S104: The frequency domain magnetic anomaly data after polarization is transformed into the spatial domain through inverse Fourier transform to obtain the polarization result.

[0060] Specifically, in this embodiment, the frequency domain magnetic anomaly data after pole-shifting is inversely Fourier transformed to the spatial domain to obtain the spatial domain magnetic anomaly data after pole-shifting, which is the pole-shifting processing result.

[0061]

[0062] In the formula, This refers to the spatial domain magnetic anomaly data after pole conversion. This is the inverse Fourier transform operation. This is the frequency domain magnetic anomaly data after polarization.

[0063] Careful analysis of the single-model and multi-model polarization results (see...) Figure 4 and Figure 5 As can be seen, the magnetic gradient is significant after polarization, the magnetic contour lines are only weakly stretched along the direction perpendicular to the magnetic declination, and the convergence phenomenon is obvious. Two roughly symmetrical negative anomaly regions are displayed outside the positive anomaly region, with obvious anomaly distribution characteristics that are close to the theoretical values ​​under perpendicular magnetization. The anomaly geometry information corresponds one-to-one with the simulation model, thus achieving the goal of low-latitude polarization.

[0064] The low-latitude magnetic anomaly polarization method described in this embodiment adaptively modifies the conventional frequency domain polarization factor to suppress the amplitude spectrum amplification effect of the conventional polarization factor and eliminate the stretching phenomenon of magnetic contour lines, thereby obtaining stable and reliable magnetic anomaly data under vertical magnetization. The processing is simple and convenient, with low computational load, and the results are reliable and stable. It can well meet the requirements of high-precision polarization processing in low-latitude regions and has good reliability, practicality and versatility.

[0065] Based on the above method, this embodiment applies the new method to the polarization processing of magnetic anomaly data from a low-latitude Indian Ocean area. A comparison before and after polarization is shown in [link to example]. Figure 6 and Figure 7 It can be seen that after processing, both positive and negative anomalies correspond and their center positions are not significantly different. There is no obvious stretching of magnetic isopleths along the direction perpendicular to the magnetic declination. The magnetic anomaly distribution characteristics are obvious, and the geometric information of the underground magnetic anomaly body is effectively enhanced, which demonstrates the reliability and stability of the polarization method proposed in this embodiment.

[0066] It should be noted that the above description describes some embodiments of this application. Other embodiments are within the scope of the appended claims. In some cases, the actions or steps recorded in the claims can be performed in a different order than that shown in the above embodiments and still achieve the desired result. Furthermore, the processes depicted in the drawings do not necessarily require a specific or sequential order to achieve the desired result. In some embodiments, multitasking and parallel processing are also possible or may be advantageous.

[0067] Based on the same inventive concept, and corresponding to any of the above embodiments, the embodiments of this application also provide a low-latitude magnetic anomaly polarization device.

[0068] like Figure 8 As shown, the low-latitude magnetic anomaly polarization device includes:

[0069] Fourier transform module 11 is configured to obtain frequency domain magnetic anomaly data by Fourier transforming the acquired spatial domain magnetic anomaly data.

[0070] The amplitude spectrum modification module 12 is configured to separate the initial amplitude spectrum and the initial phase spectrum by decomposing the frequency domain initial polar factor, and to modify the initial amplitude spectrum by a suppression modification algorithm, and to set an adaptive modification domain.

[0071] The polarization factor recombination module 13 is configured to recombine the modified amplitude spectrum with the separated initial phase spectrum to construct a recombined polarization factor, and to perform polarization processing on the frequency domain magnetic anomaly data and the recombined polarization factor to obtain polarized frequency domain magnetic anomaly data.

[0072] The inverse Fourier transform module 14 is configured to transform the frequency domain magnetic anomaly data after polarization to the spatial domain through inverse Fourier transform to obtain the polarization result.

[0073] For ease of description, the above apparatus is described in terms of its functions, divided into various modules. Of course, in implementing the embodiments of this application, the functions of each module can be implemented in one or more software and / or hardware.

[0074] The apparatus of the above embodiments is used to implement the corresponding method in any of the foregoing embodiments and has the beneficial effects of the corresponding method embodiments, which will not be repeated here.

[0075] Based on the same inventive concept, corresponding to the methods of any of the above embodiments, embodiments of this application also provide an electronic device, including a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor executes the program to implement the methods described in any of the above embodiments.

[0076] Figure 9 This embodiment illustrates a more specific hardware structure of an electronic device, which may include a processor 1010, a memory 1020, an input / output interface 1030, a communication interface 1040, and a bus 1050. The processor 1010, memory 1020, input / output interface 1030, and communication interface 1040 are interconnected internally via the bus 1050.

[0077] The processor 1010 can be implemented using a general-purpose CPU (Central Processing Unit), microprocessor, application-specific integrated circuit (ASIC), or one or more integrated circuits, and is used to execute relevant programs to implement the technical solutions provided in the embodiments of this specification.

[0078] The memory 1020 can be implemented in the form of ROM (Read Only Memory), RAM (Random Access Memory), static storage device, dynamic storage device, etc. The memory 1020 can store the operating system and other applications. When the technical solutions provided in the embodiments of this specification are implemented by software or firmware, the relevant program code is stored in the memory 1020 and is called and executed by the processor 1010.

[0079] The input / output interface 1030 is used to connect input / output modules to realize information input and output. The input / output modules can be configured as components in the device (not shown in the figure) or externally connected to the device to provide corresponding functions. Input devices may include keyboards, mice, touch screens, microphones, various sensors, etc., and output devices may include displays, speakers, vibrators, indicator lights, etc.

[0080] The communication interface 1040 is used to connect a communication module (not shown in the figure) to enable communication between this device and other devices. The communication module can communicate via wired means (such as USB, Ethernet cable, etc.) or wireless means (such as mobile network, WIFI, Bluetooth, etc.).

[0081] Bus 1050 includes a pathway for transmitting information between various components of the device, such as processor 1010, memory 1020, input / output interface 1030, and communication interface 1040.

[0082] It should be noted that although the above-described device only shows the processor 1010, memory 1020, input / output interface 1030, communication interface 1040, and bus 1050, in specific implementations, the device may also include other components necessary for normal operation. Furthermore, those skilled in the art will understand that the above-described device may only include the components necessary for implementing the embodiments of this specification, and not necessarily all the components shown in the figures.

[0083] The electronic devices described above are used to implement the corresponding methods in any of the foregoing embodiments and have the beneficial effects of the corresponding method embodiments, which will not be repeated here.

[0084] Based on the same inventive concept, corresponding to the methods of any of the above embodiments, this application also provides a non-transitory computer-readable storage medium that stores computer instructions for causing the computer to perform the methods described in any of the above embodiments.

[0085] The computer-readable medium of this embodiment includes permanent and non-permanent, removable and non-removable media, and information storage can be implemented by any method or technology. Information can be computer-readable instructions, data structures, program modules, or other data. Examples of computer storage media include, but are not limited to, phase-change memory (PRAM), static random access memory (SRAM), dynamic random access memory (DRAM), other types of random access memory (RAM), read-only memory (ROM), electrically erasable programmable read-only memory (EEPROM), flash memory or other memory technologies, CD-ROM, digital versatile optical disc (DVD) or other optical storage, magnetic tape, magnetic disk storage or other magnetic storage devices, or any other non-transfer medium that can be used to store information accessible by a computing device.

[0086] The computer instructions stored in the storage medium of the above embodiments are used to cause the computer to perform the methods described in any of the above embodiments, and have the beneficial effects of the corresponding method embodiments, which will not be repeated here.

[0087] Those skilled in the art should understand that the discussion of any of the above embodiments is merely exemplary and is not intended to imply that the scope of this application (including the claims) is limited to these examples; within the framework of this application, the technical features of the above embodiments or different embodiments can also be combined, the steps can be implemented in any order, and there are many other variations of different aspects of the embodiments of this application as described above, which are not provided in the details for the sake of brevity.

[0088] Although this application has been described in conjunction with specific embodiments thereof, many substitutions, modifications, and variations of these embodiments will be apparent to those skilled in the art from the foregoing description. For example, other memory architectures (e.g., dynamic RAM (DRAM)) may be used with the embodiments discussed.

[0089] The embodiments of this application are intended to cover all such substitutions, modifications, and variations that fall within the broad scope of the appended claims. Therefore, any omissions, modifications, equivalent substitutions, improvements, etc., made within the spirit and principles of the embodiments of this application should be included within the protection scope of this application.

Claims

1. A low latitude magnetic anomaly reduction to the pole method, characterized by, The method comprises the following steps: The acquired spatial domain magnetic anomaly data is subjected to Fourier transform to obtain frequency domain magnetic anomaly data; The initial amplitude spectrum and the initial phase spectrum are separated by decomposing the frequency domain initial polar factor, the initial amplitude spectrum is modified by a suppression modification algorithm, and an adaptive modification domain is set; The modified amplitude spectrum and the separated initial phase spectrum are recombined to construct a recombined polar factor, and the frequency domain magnetic anomaly data is subjected to polar reduction processing with the recombined polar factor to obtain polar-reduced frequency domain magnetic anomaly data; The polar-reduced frequency domain magnetic anomaly data is subjected to inverse Fourier transform to the spatial domain to obtain a polar reduction result; The modification formula is as follows: ； wherein denotes the angle of the adaptive reshaping domain, denotes the reshaping factor, denotes the reshaped amplitude spectrum; The recombination factor is represented as follows: ； In the formula, denotes the initial magnetic anomaly magnetic declination, denotes the amplitude spectrum after the modification, denotes the original amplitude spectrum.

2. The method of claim 1, wherein, The method for acquiring the spatial domain magnetic anomaly data comprises the following steps: Two magnetic anomaly forward models, including a single model and a multi-model, are constructed, wherein the multi-model is a plurality of single models combined; The model magnetic induction field output by the earth magnetic field induction processing of the single model and the multi-model is taken as initial spatial domain magnetic anomaly data, and the initial spatial domain magnetic anomaly data is subjected to simulation processing to obtain spatial domain magnetic anomaly data.

3. The method of claim 1, wherein, The separation formula is as follows: ； wherein, denotes the amplitude spectrum of the polarizing factor, , denotes the phase spectrum of the polarizing factor, , and denote the initial magnetic anomaly magnetic inclination and the initial magnetic anomaly magnetic declination, respectively, denotes the angle, which takes values from 0 to 360°, denotes the imaginary part.

4. The method of claim 1, wherein, The polar reduction operation formula is as follows: ； In the formula, denotes the frequency domain magnetic anomaly data, denotes the recombinationization factor, denotes the frequency domain magnetic anomaly data after recombinationization.

5. The method of claim 4, wherein, The inverse Fourier transform formula is as follows: ； In the formula, represents the spatial domain magnetic anomaly data after reduction to the pole, is an inverse Fourier transform operation, is the frequency domain magnetic anomaly data after reduction to the pole.

6. A low latitude magnetic anomaly reversal device, characterized by, The method comprises the following steps: The Fourier transform module is configured to subject the acquired spatial domain magnetic anomaly data to Fourier transform to obtain frequency domain magnetic anomaly data; The amplitude spectrum modification module is configured to separate the initial amplitude spectrum and the initial phase spectrum by decomposing the frequency domain initial polar factor, modify the initial amplitude spectrum by a suppression modification algorithm, and set an adaptive modification domain; The polar factor recombination module is configured to recombine the modified amplitude spectrum and the separated initial phase spectrum to construct a recombined polar factor, and subject the frequency domain magnetic anomaly data to polar reduction processing with the recombined polar factor to obtain polar-reduced frequency domain magnetic anomaly data; The inverse Fourier transform module is configured to subject the polar-reduced frequency domain magnetic anomaly data to inverse Fourier transform to the spatial domain to obtain a polar reduction result; The modification formula is as follows: ； wherein denotes the angle of the adaptive reshaping domain, denotes the reshaping factor, denotes the reshaped amplitude spectrum; The recombination factor is represented as follows: ； In the formula, denotes the initial magnetic anomaly magnetic declination, denotes the amplitude spectrum after the modification, denotes the original amplitude spectrum. 7.An electronic device comprising a memory, a processor, and a computer program stored on the memory and executable on the processor, wherein the processor implements the method of any one of claims 1-5 when executing the program.

8. A non-transitory computer-readable storage medium, comprising: The non-transitory computer readable storage medium stores computer instructions for causing a computer to execute the method of any one of claims 1-5. The non-transitory computer readable storage medium stores computer instructions for causing a computer to execute the method of any one of claims 1-5.

Images