A method and system for estimating lateral resolution of gravity and magnetic anomalies based on high-frequency information

Abstract

The application discloses a kind of estimation method and system based on high frequency information gravity-magnetic anomaly lateral resolution, comprising the following steps: obtaining gravity-magnetic anomaly data, griding gravity-magnetic anomaly data, obtaining griding gravity-magnetic anomaly data;Low-pass filtering is carried out to the griding gravity-magnetic anomaly data, and regional gravity-magnetic anomaly data are obtained;Based on the griding gravity-magnetic anomaly data and regional gravity-magnetic anomaly data, subtract to obtain residual gravity-magnetic anomaly data;Change geometric scale factor, obtain residual gravity-magnetic anomaly data of different scales;The residual gravity-magnetic anomaly data of different scales are compared and analyzed, and the preset scale that presents residual gravity-magnetic anomaly distribution characteristics is obtained, and the estimation of gravity-magnetic anomaly lateral resolution is realized.The lateral resolution of gravity anomaly or magnetic anomaly estimated by the method can evaluate the related geological problems that can be solved by existing gravity-magnetic data, and can improve the efficiency of gravity-magnetic data processing and interpretation.

Description

Technical Field

[0001] This application belongs to the field of geophysical exploration technology, specifically relating to an estimation method and system for the lateral resolution of gravity and magnetic anomalies based on high-frequency information. Background Technology

[0002] Gravity and magnetic anomaly interpretation relies on the characteristics identified by these anomalies for geological interpretation. Gravity anomalies reflect the density characteristics of geological bodies, while magnetic anomalies reflect their magnetic properties. Gravity and magnetic anomalies have high planar resolution and are well-suited for identifying the planar distribution of subsurface lithology and structures. The larger the scale of gravity and magnetic field measurements, the higher the planar resolution of gravity and magnetic anomalies, i.e., the higher their lateral resolution. However, gravity and magnetic measurement data are classified, making it difficult for researchers to obtain field data, especially large-scale field measurements. In gravity and magnetic data processing and interpretation research, researchers often obtain gridded data from unknown field measurement scales, resulting in unknown lateral resolution and affecting the accuracy of gravity and magnetic data processing and interpretation. Therefore, estimating the lateral resolution of gravity and magnetic data under unknown field measurement scales is of great significance. Summary of the Invention

[0003] This application proposes a method and system for estimating the lateral resolution of gravity and magnetic anomalies based on high-frequency information, solving the aforementioned problems in the prior art. Its purpose is to estimate the lateral resolution of gravity and magnetic anomalies based on high-frequency information acquired from gravity and magnetic data, thereby improving the efficiency of gravity and magnetic data processing and interpretation, and allowing direct application to gravity and magnetic detection research.

[0004] To achieve the above objectives, this application provides the following solution:

[0005] A method for estimating the lateral resolution of gravity and magnetic anomalies based on high-frequency information includes the following steps:

[0006] S1: Acquire gravity and magnetic anomaly data, and grid the gravity and magnetic anomaly data to obtain gridded gravity and magnetic anomaly data;

[0007] S2: Perform low-pass filtering on the gridded gravity and magnetic anomaly data to obtain regional gravity and magnetic anomaly data;

[0008] S3: Subtract the gridded gravity and magnetic anomaly data from the regional gravity and magnetic anomaly data to obtain the remaining gravity and magnetic anomaly data;

[0009] S4: Change the geometric scale factor of the gridded gravity and magnetic anomaly data, and repeat S2 and S3 to obtain residual gravity and magnetic anomaly data at different scales;

[0010] S5: Compare and analyze the residual gravity and magnetic anomaly data at different scales to obtain a preset scale that presents the distribution characteristics of the residual gravity and magnetic anomaly. Based on the preset scale, estimate the lateral resolution of the gravity and magnetic anomaly.

[0011] Preferably, the gravity and magnetic anomaly data includes gravity anomaly data and magnetic anomaly data;

[0012] The gridded gravity and magnetic anomaly data includes gridded gravity anomaly data and gridded magnetic anomaly data;

[0013] The regional gravity and magnetic anomaly data includes gravity anomaly data and regional magnetic anomaly data;

[0014] The residual gravity and magnetic anomaly data includes residual gravity anomaly data and residual magnetic anomaly data.

[0015] Preferably, the low-pass filter employs a regularized stable filter factor.

[0016] Preferably, the method for obtaining the regularized stable filter factor is as follows:

[0017] The gravity field observation surface is set as the horizontal plane with h = 0, and the residual gravity and magnetic anomaly is Δg(x,y);

[0018] Based on the horizontal plane and the residual gravity and magnetic anomaly Δg(x,y), the gravity and magnetic anomaly spectrum is obtained;

[0019] The regularized stable filter factor is obtained by multiplying the gravity and magnetic anomaly spectrum with the stability factor.

[0020] Preferably, the formula for the regularized stable filter factor is:

[0021]

[0022] Where, λ x It is a parameter representing the fundamental wavelength, equal to the line spacing of the survey area multiplied by the number of lines; β≥2, f0=1 / λ0, λ0 is the wavelength of the smallest wavenumber component constituting the local interference anomaly, wavenumber u=m / λ x v = n / λ y f = (u 2 +v 2 ) 1 / 2 m and n are the wavenumbers on the u-axis and v-axis of the two-dimensional coordinate system in the frequency domain, respectively, and λ x , λ y This is the fundamental wavelength.

[0023] This application also discloses an estimation system for the lateral resolution of gravity and magnetic anomalies based on high-frequency information, including a gridding module, a regional anomaly data acquisition module, a residual anomaly data acquisition module, a parameter changing module, and a comparison module;

[0024] The gridding module is used to acquire gravity and magnetic anomaly data, and to grid the gravity and magnetic anomaly data to obtain gridded gravity and magnetic anomaly data.

[0025] The regional anomaly data acquisition module is used to perform low-pass filtering on the gridded gravity and magnetic anomaly data to obtain regional gravity and magnetic anomaly data.

[0026] The remaining abnormal data acquisition module is used to subtract the gridded gravity and magnetic anomaly data from the regional gravity and magnetic anomaly data to obtain the remaining gravity and magnetic anomaly data.

[0027] The parameter changing module is used to change the geometric scale factor of the gridded gravity and magnetic anomaly data, repeating the operations of the regional anomaly data acquisition module and the residual anomaly data acquisition module to obtain residual gravity and magnetic anomaly data at different scales.

[0028] The comparison module is used to compare and analyze the residual gravity and magnetic anomaly data at different scales to obtain a preset scale that presents the distribution characteristics of the residual gravity and magnetic anomalies. Based on the preset scale, the lateral resolution of the gravity and magnetic anomalies is estimated.

[0029] Preferably, the gravity and magnetic anomaly data includes gravity anomaly data and magnetic anomaly data;

[0030] The gridded gravity and magnetic anomaly data includes gridded gravity anomaly data and gridded magnetic anomaly data;

[0031] The regional gravity and magnetic anomaly data includes gravity anomaly data and regional magnetic anomaly data;

[0032] The residual gravity and magnetic anomaly data includes residual gravity anomaly data and residual magnetic anomaly data.

[0033] Preferably, the low-pass filter employs a regularized stable filter factor.

[0034] Preferably, the regional anomaly data acquisition module includes a setting unit, an anomaly spectrum acquisition unit, and a regularization unit;

[0035] The setting unit is used to set the gravity field observation surface to a horizontal plane with h = 0 and the residual gravity and magnetic anomaly to Δg(x,y);

[0036] The abnormal spectrum acquisition unit is used to obtain the gravity and magnetic anomaly spectrum based on the horizontal plane and the residual gravity and magnetic anomaly Δg(x,y);

[0037] The regularization unit is used to obtain a regularized stable filter factor based on the product of the gravity and magnetic anomaly spectrum and the stability factor.

[0038] Preferably, the formula for the regularized stable filter factor is:

[0039]

[0040] Where, λ x It is a parameter representing the fundamental wavelength, equal to the line spacing of the survey area multiplied by the number of lines; β≥2, f0=1 / λ0, λ0 is the wavelength of the smallest wavenumber component constituting the local interference anomaly, wavenumber u=m / λ x v = n / λ y f = (u 2 +v 2 ) 1 / 2 m and n are the wavenumbers on the u-axis and v-axis of the two-dimensional coordinate system in the frequency domain, respectively, and λ x , λ y This is the fundamental wavelength.

[0041] The beneficial effects of this application are as follows: The technical solution of this application has the advantage that it does not require knowledge of the scale of gravity and magnetic field measurements, and the lateral resolution of gravity or magnetic anomalies can be estimated. The lateral resolution of gravity or magnetic anomalies estimated using this method can be used to assess relevant geological problems that existing gravity and magnetic data can solve, thereby improving the efficiency of gravity and magnetic data processing and interpretation. This application has broad application potential and practical value. Attached Figure Description

[0042] To more clearly illustrate the technical solutions of this application, the drawings used in the embodiments are briefly introduced below. Obviously, the drawings described below are only some embodiments of this application. For those skilled in the art, other drawings can be obtained based on these drawings without creative effort.

[0043] Figure 1 This is a flowchart of the estimation method for the lateral resolution of gravity and magnetic anomalies based on high-frequency information, as described in Embodiment 1 of this application.

[0044] Figure 2 This is a Bouguer gravity anomaly map of block A calculated based on gravity field measurement data in Embodiment 1 of this application;

[0045] Figure 3 This is a residual gravity anomaly map of block 1A in Embodiment 1 of this application with a geometric scale factor of 4km;

[0046] Figure 4 This is a residual gravity anomaly map of block 1A in Embodiment 1 of this application with a geometric scale factor of 6 km;

[0047] Figure 5 This is a residual gravity anomaly map of block 1A in Embodiment 1 of this application with a geometric scale factor of 8 km;

[0048] Figure 6 This is a residual gravity anomaly map of block 1A in Embodiment 1 of this application with a geometric scale factor of 10km;

[0049] Figure 7 This is a residual gravity anomaly map of block 1A in Embodiment 1 of this application with a geometric scale factor of 12km;

[0050] Figure 8 This is a residual gravity anomaly map of block 1A in Embodiment 1 of this application with a geometric scale factor of 14km;

[0051] Figure 9 This is a residual gravity anomaly map of block 1A in Embodiment 1 of this application, with a geometric scale factor of 16km. Detailed Implementation

[0052] 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, and 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.

[0053] To make the above-mentioned objectives, features and advantages of this application more apparent and understandable, the application will be further described in detail below with reference to the accompanying drawings and specific embodiments.

[0054] Example 1: An estimation method for lateral resolution of gravity and magnetic anomalies based on high-frequency information

[0055] like Figure 1 As shown, a method for estimating the lateral resolution of gravity and magnetic anomalies based on high-frequency information includes the following steps:

[0056] S1: Acquire gravity and magnetic anomaly data, and grid the gravity and magnetic anomaly data to obtain gridded gravity and magnetic anomaly data;

[0057] S2: Low-pass filtering is applied to the gridded gravity and magnetic anomaly data to obtain regional gravity and magnetic anomaly data;

[0058] S3: Subtract the gridded gravity and magnetic anomaly data from the regional gravity and magnetic anomaly data to obtain the remaining gravity and magnetic anomaly data;

[0059] S4: Change the geometric scale factor of the gridded gravity and magnetic anomaly data, and repeat S2 and S3 to obtain residual gravity and magnetic anomaly data at different scales.

[0060] S5: Compare and analyze residual gravity and magnetic anomaly data at different scales to obtain a preset scale that presents the distribution characteristics of residual gravity and magnetic anomalies. Based on the obtained preset scale, estimate the lateral resolution of gravity and magnetic anomalies.

[0061] Specifically, gravity and magnetic anomaly data include gravity anomaly data and magnetic anomaly data;

[0062] The gridded gravity and magnetic anomaly data includes gridded gravity anomaly data Δg and magnetic anomaly data ΔT;

[0063] Regional gravity and magnetic anomaly data include gravity anomaly data Δg 区域 and regional magnetic anomaly data ΔT 区域 ;

[0064] Residual gravity and magnetic anomaly data include residual gravity anomaly data Δg 剩余 and residual magnetic anomaly data ΔT 剩余 .

[0065] The low-pass filter uses a regularized stable filter factor.

[0066] The filtering parameters λ0 and f0 of the regularized stable filtering factor are scale-matched with the local anomaly field to be eliminated, and are directly measured from the original anomaly profile map and the plane contour map.

[0067] The method for obtaining the regularized stable filter factor is as follows:

[0068] The gravity field observation surface is set as the horizontal plane with h = 0, and the residual gravity and magnetic anomaly is Δg(x,y);

[0069] Based on the gravity field observation surface and the residual gravity and magnetic anomalies, the gravity and magnetic anomaly spectrum Δg(u,v) is obtained;

[0070] The regularized stable filter factor is obtained by multiplying the gravity and magnetic anomaly spectrum Δg(u,v) with the stability factor.

[0071] The formula for the regularized stable filter factor is:

[0072]

[0073] Where, λ x It is a parameter representing the fundamental wavelength, equal to the line spacing of the survey area multiplied by the number of lines; β≥2, f0=1 / λ0, λ0 is the wavelength of the smallest wavenumber component constituting the local interference anomaly, wavenumber u=m / λ x v = n / λ y f = (u 2 +v 2 ) 1 / 2 m and n are the wavenumbers on the u-axis and v-axis of the two-dimensional coordinate system in the frequency domain, respectively, and λ x, λ y λ0 is the fundamental wavelength. Specifically, λ0 is the maximum length, including both positive and negative anomaly segments, measured along the vertical anomaly direction from among numerous local interference anomalies synthesized from unstable high wavenumber anomalies within the survey area. It is equivalent to the wavelength of the smallest wavenumber component that makes up the local interference anomaly.

[0074] Specifically, in the wavenumber domain, in order to suppress local anomaly interference and achieve stable operation during anomaly continuation, differentiation, and interface burial depth inversion, it is necessary to multiply the gravity and magnetic anomaly spectrum by a stabilization factor to obtain a regularized stabilization filter factor.

[0075] Specifically, in S4, the geometric scale factor of the gridded gravity and magnetic anomaly data is changed, and S2 and S3 are repeated to obtain residual gravity and magnetic anomaly data at different scales: Figure 2 This is a Bouguer gravity anomaly map of Block A calculated based on gravity field measurement data (line spacing 4km, point spacing 1km, contour unit 10). -5 m / s 2 Based on the centerline and point spacing distribution information in the Technical Specification for Onshore Gravity Exploration (SY / T5819-2010), the scale of this field survey is 1:200,000.

[0076] Figure 3 Residual gravity anomaly map for block 1A (geometric scale factor 4km, contour spacing 0.5×10⁻⁶). -5 m / s 2 The residual gravity anomaly in this figure is randomly distributed and cannot be used as the minimum scale for estimating the lateral resolution of the anomaly.

[0077] Figure 4 Residual gravity anomaly map for block 1A (geometric scale factor 6 km, contour spacing 0.5 × 10⁻⁶). -5 m / s 2 The residual gravity anomaly in this figure is randomly distributed and cannot be used as the minimum scale for estimating the lateral resolution of the anomaly.

[0078] Figure 5 Residual gravity anomaly map for block 1A (geometric scale factor 8 km, contour spacing 0.5 × 10⁻⁶). -5 m / s 2 The distribution characteristics of the residual gravity anomaly in this figure are not very obvious and should not be used as the minimum scale for estimating the lateral resolution of the anomaly.

[0079] Figure 6 Residual gravity anomaly map for block 1A (geometric scale factor 10km, contour spacing 0.5×10⁻⁶). -5 m / s 2The figure shows a clear distribution of residual gravity anomalies and can be used as the smallest scale for estimating the lateral resolution of the anomalies, which is consistent with the lateral resolution of gravity data at a scale of 1:200,000.

[0080] Figure 7 Residual gravity anomaly map for block 1A (geometric scale factor 12km, contour spacing 0.5×10⁻⁶). -5 m / s 2 This figure shows a clear distribution of residual gravity anomalies, but it is not used as the minimum scale for estimating the lateral resolution of the anomalies.

[0081] Figure 8 Residual gravity anomaly map for block 1A (geometric scale factor 14km, contour spacing 0.5×10⁻⁶). -5 m / s 2 This figure shows a clear distribution of residual gravity anomalies, but it is not used as the minimum scale for estimating the lateral resolution of the anomalies.

[0082] Figure 9 Residual gravity anomaly map for block 1A (geometric scale factor 16km, contour spacing 0.5×10⁻⁶). -5 m / s 2 The figure shows a clear distribution of residual gravity anomalies, but it is not used as the minimum scale for estimating the lateral resolution of the anomalies.

[0083] Example 2: An estimation system for lateral resolution of gravity and magnetic anomalies based on high-frequency information

[0084] An estimation system based on the lateral resolution of gravity and magnetic anomalies using high-frequency information includes a gridding module, a regional anomaly data acquisition module, a residual anomaly data acquisition module, a parameter changing module, and a comparison module.

[0085] A gridding module is used to acquire gravity and magnetic anomaly data, and to grid the gravity and magnetic anomaly data to obtain the gridded gravity and magnetic anomaly data.

[0086] The regional anomaly data acquisition module is used to perform low-pass filtering on the gridded gravity and magnetic anomaly data to obtain regional gravity and magnetic anomaly data.

[0087] The remaining anomaly data acquisition module is used to obtain the remaining gravity and magnetic anomaly data by subtracting the gridded gravity and magnetic anomaly data from the regional gravity and magnetic anomaly data.

[0088] The parameter changing module is used to change the geometric scale factor of the gridded gravity and magnetic anomaly data. It repeats the operations of the regional anomaly data acquisition module and the residual anomaly data acquisition module to obtain residual gravity and magnetic anomaly data at different scales.

[0089] The comparison module is used to compare and analyze residual gravity and magnetic anomaly data at different scales, obtain a preset scale that presents the distribution characteristics of residual gravity and magnetic anomalies, and estimate the lateral resolution of gravity and magnetic anomalies based on the preset scale.

[0090] Specifically, gravity and magnetic anomaly data include gravity anomaly data and magnetic anomaly data;

[0091] The gridded gravity and magnetic anomaly data includes gridded gravity anomaly data Δg and magnetic anomaly data ΔT;

[0092] Regional gravity and magnetic anomaly data include gravity anomaly data Δg 区域 and regional magnetic anomaly data ΔT 区域 ;

[0093] Residual gravity and magnetic anomaly data include residual gravity anomaly data Δg 剩余 and residual magnetic anomaly data ΔT 剩余 .

[0094] The low-pass filter uses a regularized stable filter factor.

[0095] The filtering parameters λ0 and f0 of the regularized stable filtering factor are scale-matched with the local anomaly field to be eliminated, and are directly measured from the original anomaly profile map and the plane contour map.

[0096] The regional anomaly data acquisition module includes a setting unit, an anomaly spectrum acquisition unit, and a regularization unit;

[0097] The setting unit is used to set the gravity field observation surface to the horizontal plane with h = 0 and the residual gravity and magnetic anomaly to Δg(x,y);

[0098] Anomaly spectrum acquisition unit is used to obtain the gravity and magnetic anomaly spectrum based on the horizontal plane and the residual gravity and magnetic anomaly Δg(x,y);

[0099] The regularization unit is used to obtain the regularized stable filter factor based on the product of the gravity and magnetic anomaly spectrum and the stability factor.

[0100] The formula for the regularized stable filter factor is:

[0101]

[0102] Where, λ x It is a parameter representing the fundamental wavelength, equal to the line spacing of the survey area multiplied by the number of lines; β≥2, f0=1 / λ0, λ0 is the wavelength of the smallest wavenumber component constituting the local interference anomaly, wavenumber u=m / λ x v = n / λ y f = (u 2 +v 2 ) 1 / 2m and n are the wavenumbers on the u-axis and v-axis of the two-dimensional coordinate system in the frequency domain, respectively, and λ x , λ y This is the fundamental wavelength.

[0103] Specifically, in the wavenumber domain, in order to suppress local anomaly interference and achieve stable operation during anomaly continuation, differentiation, and interface burial depth inversion, it is necessary to multiply the gravity and magnetic anomaly spectrum by a stabilization factor to obtain a regularized stabilization filter factor.

[0104] Specifically, by changing the geometric scale factor of the gridded gravity and magnetic anomaly data in the parameter module, and repeating the operations of the regional anomaly data acquisition module and the residual anomaly data acquisition module, residual gravity and magnetic anomaly data at different scales can be obtained.

[0105] The embodiments described above are merely preferred embodiments of this application and are not intended to limit the scope of this application. Any modifications and improvements made to the technical solutions of this application by those skilled in the art without departing from the spirit of this application shall fall within the protection scope defined by the claims of this application.

Claims

1. A method for estimating the lateral resolution of gravity and magnetic anomalies based on high-frequency information, characterized in that, Includes the following steps: S1: Acquire gravity and magnetic anomaly data, and grid the gravity and magnetic anomaly data to obtain gridded gravity and magnetic anomaly data; S2: Perform low-pass filtering on the gridded gravity and magnetic anomaly data to obtain regional gravity and magnetic anomaly data; S3: Subtract the gridded gravity and magnetic anomaly data from the regional gravity and magnetic anomaly data to obtain the remaining gravity and magnetic anomaly data; S4: Change the geometric scale factor of the gridded gravity and magnetic anomaly data, and repeat S2 and S3 to obtain residual gravity and magnetic anomaly data at different scales; S5: Compare and analyze the residual gravity and magnetic anomaly data at different scales to obtain a preset scale that presents the distribution characteristics of the residual gravity and magnetic anomaly. Based on the preset scale, estimate the lateral resolution of the gravity and magnetic anomaly. The low-pass filter employs a regularized stable filter factor; The formula for the regularized stable filter factor is: Where, λ x It is a parameter representing the fundamental wavelength, equal to the line spacing of the survey area multiplied by the number of lines; β≥2, f0=1 / λ0, λ0 is the wavelength of the smallest wavenumber component constituting the local interference anomaly, and the wavenumber u=m / λ x v=n / λ y f=(u 2 +v 2 ) 1 / 2 m and n are the wavenumbers on the u-axis and v-axis of the two-dimensional coordinate system in the frequency domain, respectively, and λ x , λ y This is the fundamental wavelength.

2. The estimation method for lateral resolution of gravity and magnetic anomalies based on high-frequency information according to claim 1, characterized in that, The gravity and magnetic anomaly data includes gravity anomaly data and magnetic anomaly data; The gridded gravity and magnetic anomaly data includes gridded gravity anomaly data and gridded magnetic anomaly data; The regional gravity and magnetic anomaly data includes regional gravity anomaly data and regional magnetic anomaly data; The residual gravity and magnetic anomaly data includes residual gravity anomaly data and residual magnetic anomaly data.

3. The estimation method for lateral resolution of gravity and magnetic anomalies based on high-frequency information according to claim 1, characterized in that, The method for obtaining the regularized stable filter factor is as follows: The gravity field observation surface is set as the horizontal plane with h=0, and the residual gravity and magnetic anomaly is Δg(x,y); Based on the horizontal plane and the residual gravity and magnetic anomaly Δg(x,y), the gravity and magnetic anomaly spectrum is obtained; The regularized stable filter factor is obtained by multiplying the gravity and magnetic anomaly spectrum with the stability factor.

4. An estimation system for the lateral resolution of gravity and magnetic anomalies based on high-frequency information, characterized in that, It includes a gridding module, a regional anomaly data acquisition module, a remaining anomaly data acquisition module, a parameter modification module, and a comparison module; The gridding module is used to acquire gravity and magnetic anomaly data, and to grid the gravity and magnetic anomaly data to obtain gridded gravity and magnetic anomaly data. The regional anomaly data acquisition module is used to perform low-pass filtering on the gridded gravity and magnetic anomaly data to obtain regional gravity and magnetic anomaly data. The remaining abnormal data acquisition module is used to subtract the gridded gravity and magnetic anomaly data from the regional gravity and magnetic anomaly data to obtain the remaining gravity and magnetic anomaly data. The parameter changing module is used to change the geometric scale factor of the gridded gravity and magnetic anomaly data, repeating the operations of the regional anomaly data acquisition module and the residual anomaly data acquisition module to obtain residual gravity and magnetic anomaly data at different scales. The comparison module is used to compare and analyze the residual gravity and magnetic anomaly data at different scales to obtain a preset scale that presents the distribution characteristics of the residual gravity and magnetic anomaly, and based on the preset scale, to estimate the lateral resolution of the gravity and magnetic anomaly. The low-pass filter employs a regularized stable filter factor; The formula for the regularized stable filter factor is: Where, λ x It is a parameter representing the fundamental wavelength, equal to the line spacing of the survey area multiplied by the number of lines; β≥2, f0=1 / λ0, λ0 is the wavelength of the smallest wavenumber component constituting the local interference anomaly, and the wavenumber u=m / λ x v=n / λ y f=(u 2 +v 2 ) 1 / 2 m and n are the wavenumbers on the u-axis and v-axis of the two-dimensional coordinate system in the frequency domain, respectively, and λ x , λ y This is the fundamental wavelength.

5. The estimation system for lateral resolution of gravity and magnetic anomalies based on high-frequency information as described in claim 4, characterized in that, The gravity and magnetic anomaly data includes gravity anomaly data and magnetic anomaly data; The gridded gravity and magnetic anomaly data includes gridded gravity anomaly data and gridded magnetic anomaly data; The regional gravity and magnetic anomaly data includes regional gravity anomaly data and regional magnetic anomaly data; The residual gravity and magnetic anomaly data includes residual gravity anomaly data and residual magnetic anomaly data.

6. The estimation system for lateral resolution of gravity and magnetic anomalies based on high-frequency information as described in claim 4, characterized in that, The regional anomaly data acquisition module includes a setting unit, an anomaly spectrum acquisition unit, and a regularization unit; The setting unit is used to set the gravity field observation surface to a horizontal plane with h=0 and the residual gravity and magnetic anomaly to Δg(x,y); The abnormal spectrum acquisition unit is used to obtain the gravity and magnetic anomaly spectrum based on the horizontal plane and the residual gravity and magnetic anomaly Δg(x,y); The regularization unit is used to obtain a regularized stable filter factor based on the product of the gravity and magnetic anomaly spectrum and the stability factor.

Images