High-density gravity-magnetic force composite acquisition processing method and device

Abstract

The application discloses a high-density gravity-magnetic force composite acquisition processing method, which comprises the following steps: establishing a first measuring network and a second measuring network in a measuring area, wherein the first measuring network and the second measuring network are gravity-magnetic force general point measuring networks, the measuring network density of the first measuring network is different from that of the second measuring network, and the second measuring network comprises two measuring line orthogonal measuring networks; collecting gravity-magnetic force data in the first measuring network and the second measuring network respectively, and performing preliminary processing on the gravity-magnetic force data, wherein the preliminary processing comprises preliminary merging and interference value correction; and performing secondary merging on the data of the first measuring network and the second measuring network after the preliminary processing to generate a composite data body. The method realizes high-density gravity-magnetic force data acquisition processing, and has important practical and popularization and application values in complex oil and gas target gravity-magnetic force exploration. The application further discloses a high-density gravity-magnetic force composite acquisition processing device.

Description

Technical Field

[0001] This invention relates to the field of geophysical exploration technology, and in particular to a high-density gravity and magnetic composite acquisition and processing method and apparatus. Background Technology

[0002] Gravity and magnetic exploration are two important methods in oil and gas exploration, playing a crucial role in the study of regional tectonics, basement lithology, rift basins, local structural zones, buried hills, rift valleys, and piedmont zones. With the continuous development of oil and gas exploration, deep and complex targets, basement lithology, and volcanic reservoir targets have become important areas of focus in oil and gas exploration in my country and worldwide. Developing and innovating new high-density gravity and magnetic data acquisition, processing, and interpretation technologies for complex targets is of significant practical importance. Currently, gravity and magnetic data acquisition in oil and gas exploration conventionally uses 1:50,000 high-precision gravity and magnetic methods, typically with a 500m×500m area and completed in a single acquisition. However, this acquisition method has relatively low accuracy and suffers from problems such as large random errors, low verification rates, and significant interference effects.

[0003] Therefore, there is still a need to improve existing methods and devices for combined gravity and magnetic force acquisition and processing. Summary of the Invention

[0004] In view of this, the purpose of this invention is to propose a high-density gravity and magnetic composite acquisition and processing method and apparatus. The method and apparatus of this invention eliminates interference through composite acquisition and achieves high-density gravity and magnetic data acquisition and processing by combining multiple corrections, which is suitable for the exploration of complex oil and gas targets.

[0005] To achieve the above objectives, one aspect of the present invention provides a high-density gravity and magnetic force composite acquisition and processing method, comprising the following steps:

[0006] A first and a second measurement network are established in the measurement area. Both the first and the second measurement networks are gravity and magnetic field measurement networks. The measurement network densities of the first and the second measurement networks are different. The second measurement network includes measurement networks with two orthogonal directions.

[0007] Gravity and magnetic force data were collected in the first and second measurement networks respectively. The gravity and magnetic force data collected in the first and second measurement networks were preliminarily processed, including preliminary merging and interference value correction. The data from the first and second measurement networks after the preliminary processing were merged again to generate a composite data volume.

[0008] In some implementations, the method further includes:

[0009] Data gridding is performed on the composite data volume to generate a gridded data volume;

[0010] The gridded data volume is subjected to displacement direction filtering correction to generate the first correction data.

[0011] In some implementations, the method further includes:

[0012] Magnetic anomaly displacement correction is performed on the first correction data to generate a high-density gravity and magnetic force composite data volume.

[0013] In some implementations, displacement direction filtering correction includes:

[0014] The gridded data volume is analyzed and verified using gravity and magnetic anomaly diagrams to obtain error data. The error data is then subjected to directional filtering along the survey line to obtain directional filtering data. The directional filtering data is multiplied by the directional filtering potential change coefficient to generate the filtering correction data value. The filtering correction data value is then removed from the gridded data volume to generate the first correction data.

[0015] In some implementations, the survey line is used as the displacement region, the directional filtering displacement coefficient within the displacement region is 1, the directional filtering displacement coefficient outside the displacement region is 0, and the directional filtering displacement coefficient between the displacement region and the displacement region transitions between 1 and 0.

[0016] In some implementations, magnetic anomaly displacement correction includes:

[0017] Elevation data and magnetic anomaly data are obtained based on gravity and magnetic data. Residual magnetic anomaly data are further calculated based on the magnetic anomaly data. The magnetic topographic influence area is delineated by superimposing elevation data and magnetic anomaly data and the magnetic displacement coefficient is determined. The magnetic displacement coefficient is multiplied by the residual magnetic anomaly to generate a magnetic topographic related displacement correction value. The magnetic topographic related displacement correction value is subtracted from the first correction data to generate a high-density gravity and magnetic composite data volume.

[0018] In some implementations, the magnetic displacement coefficient of the area affected by the magnetic topography is taken as 1, while the magnetic displacement coefficient of other areas is taken as 0.

[0019] In some implementations, the azimuth angles of the first and second surveying networks are different.

[0020] In some implementations, the mesh density of the second mesh is higher than that of the first mesh.

[0021] In some implementations, the density of the first survey network is 500m × 500m, and the density of the second survey network is 200m × 500m.

[0022] In some implementations, the initial merging includes:

[0023] At the overlapping location of the two orthogonal measurement networks in the second measurement network, the data of the two orthogonal measurement networks are summed and averaged, and recorded as the data at that overlapping location;

[0024] At the overlapping position of the first and second measurement networks, the data from the first and second measurement networks are summed and averaged, and this average is recorded as the data at that overlapping position.

[0025] In some implementations, interference value correction includes:

[0026] The gravity and magnetic data of the first and second measurement networks are respectively processed into grids to generate first gridded data and second gridded data. The first gridded data and the second gridded data use the same grid spacing to facilitate interpolation. Interference points in the second gridded data are identified in conjunction with interference source identification, and the data corresponding to the interference points are replaced with the corresponding data in the first gridded data.

[0027] In some implementations, the interference source is a human-caused interference source.

[0028] Another aspect of this invention provides a high-density gravity and magnetic force composite acquisition and processing device, comprising:

[0029] The measurement network construction module is configured to establish a first measurement network and a second measurement network in the measurement area. Both the first and second measurement networks are gravity and magnetic field measurement networks. The measurement network densities of the first and second measurement networks are different. The second measurement network includes two measurement networks with orthogonal measurement lines.

[0030] The data acquisition module collects gravity and magnetic data in the first and second measurement networks respectively. The gravity and magnetic data includes the vertical coordinate, horizontal coordinate, elevation, gravity field value, and magnetic field value of the measurement point.

[0031] The first processing module is configured to perform preliminary processing on the gravity and magnetic data collected from the first and second measurement networks. The preliminary processing includes preliminary merging and interference value correction. The data from the first and second measurement networks after the preliminary processing are then merged again to generate a composite data volume.

[0032] In some embodiments, the apparatus further includes:

[0033] The second processing module is configured to perform data gridding on the composite data volume to generate a gridded data volume.

[0034] The first correction module is configured to perform displacement direction filtering correction on the gridded data volume to generate first correction data.

[0035] In some embodiments, the apparatus further includes:

[0036] The second correction module is configured to perform magnetic anomaly displacement correction on the first correction data to generate a high-density gravity and magnetic force composite data volume.

[0037] In some implementations, the first correction module obtains error data by performing gravity and magnetic anomaly analysis on the gridded data volume, performs directional filtering on the error data along the survey line direction to obtain directional filtered data, multiplies the directional filtered data with the directional filtered potential change coefficient to generate filtered correction data value, and removes the filtered correction data value from the gridded data volume to generate the first correction data.

[0038] In some implementations, the survey line is used as the displacement region, the directional filtering displacement coefficient within the displacement region is 1, the directional filtering displacement coefficient outside the displacement region is 0, and the directional filtering displacement coefficient between the displacement region and the displacement region transitions between 1 and 0.

[0039] In some implementations, the second correction module acquires elevation data and magnetic anomaly data based on gravity and magnetic data, further calculates residual magnetic anomaly data based on the magnetic anomaly data, delineates the magnetic topographic influence area by superimposing elevation data and magnetic anomaly data and determines the magnetic displacement coefficient, multiplies the magnetic displacement coefficient with the residual magnetic anomaly to generate a magnetic topographic related displacement correction value, and subtracts the magnetic topographic related displacement correction value from the first correction data to generate a high-density gravity and magnetic composite data volume.

[0040] In some implementations, the magnetic displacement coefficient of the area affected by the magnetic topography is taken as 1, while the magnetic displacement coefficient of other areas is taken as 0.

[0041] The present invention has at least the following beneficial technical effects:

[0042] This invention establishes a first and a second survey network with different densities and survey line directions, and further establishes two orthogonal survey networks within the second survey network. This enables composite acquisition of gravity and magnetic data in the same measurement area, achieving high-density gravity and magnetic data acquisition, and effectively avoiding errors caused by near-surface human interference.

[0043] The gravity and magnetic composite acquisition and processing method provided by this invention is applicable to the exploration of complex oil and gas targets. It has important practical and promotional application value in the gravity and magnetic exploration of complex oil and gas targets. It is a further extension of high-precision gravity and magnetic exploration technology and an important innovation of existing gravity and magnetic acquisition and processing methods. It provides an important new method for gravity and magnetic data acquisition and processing.

[0044] The apparatus provided by this invention has the same beneficial effects as the method. Attached Figure Description

[0045] To more clearly illustrate the technical solutions in the embodiments of the present invention or the prior art, the drawings used in the description of the embodiments or the prior art will be briefly introduced below. Obviously, the drawings described below are only some embodiments of the present invention. For those skilled in the art, other embodiments can be obtained based on these drawings without creative effort.

[0046] Figure 1 This is a schematic diagram of an embodiment of the high-density gravity and magnetic force composite acquisition and processing method provided by the present invention;

[0047] Figure 2 This is a schematic diagram of an embodiment of the high-density gravity and magnetic force composite acquisition and processing device provided by the present invention. Detailed Implementation

[0048] To make the objectives, technical solutions, and advantages of the present invention clearer, the embodiments of the present invention will be further described in detail below with reference to specific examples and the accompanying drawings.

[0049] It should be noted that all uses of "first" and "second" in the embodiments of the present invention are for the purpose of distinguishing two entities or parameters with the same name but different names. It is clear that "first" and "second" are only for the convenience of expression and should not be construed as limiting the embodiments of the present invention. Subsequent embodiments will not explain this in detail.

[0050] Based on the above objectives, the first aspect of the present invention proposes a high-density gravity and magnetic force composite acquisition and processing method. Figure 1 The diagram shown is a schematic representation of an embodiment of the high-density gravity and magnetic force composite acquisition and processing method provided by the present invention. Figure 1 As shown, the high-density gravity and magnetic force composite acquisition and processing method of this invention includes the following steps:

[0051] S1. Establish a first measurement network and a second measurement network in the measurement area. Both the first and second measurement networks are gravity and magnetic force general point measurement networks. The measurement network densities of the first and second measurement networks are different. The second measurement network includes measurement networks with two orthogonal direction measurement lines.

[0052] S2. Collect gravity and magnetic force data in the first and second measurement networks respectively, and perform preliminary processing on the gravity and magnetic force data collected in the first and second measurement networks. The preliminary processing includes preliminary merging and interference value correction. Then, merge the data from the first and second measurement networks after the preliminary processing to generate a composite data volume.

[0053] S1 establishes a first and a second survey network within the measurement area. Both networks are gravity and magnetic spectral survey networks. The scale and network density of the first and second networks differ. Network density represents the number of survey points per unit area, typically expressed as the product of the line spacing and the point spacing; for example, 200m × 500m represents a line spacing of 200m and a point spacing of 500m. Preferably, the scale and network density of the first network are smaller than those of the second network. For instance, the scale of the first network is 1:50,000, and its network density is 500m × 500m, while the scale of the second network is approximately 1:20,000, and its network density is 200m × 500m. The first network is a unidirectional network, for example, with an azimuth of 0°. The second network comprises two orthogonal survey networks; preferably, the azimuths of the two networks are 45° and 315°, respectively. Secondly, those skilled in the art should understand that the scale, density, and azimuth of the survey network given in this embodiment are merely illustrative and can be adjusted according to the actual exploration area, all of which fall within the protection scope of this invention.

[0054] After the measurement network is established, proceed to S2 and collect gravity and magnetic data in the first and second measurement networks respectively. The first and second measurement networks can be used to collect data in different years, with the first measurement network being implemented in an earlier year and the second measurement network in a later year. The gravity and magnetic data include the vertical coordinates, horizontal coordinates, elevations, gravity field values, and magnetic field values ​​of the measurement points.

[0055] Furthermore, the collected gravity and magnetic data undergo preliminary processing, including data merging and interference correction to generate a composite data volume. In geophysical exploration, gravity and magnetic exploration data generally require only one data point per coordinate location. However, due to the layout of the first and second survey networks in this embodiment, some survey point locations overlap. Furthermore, the orthogonal intersection of the two survey networks in the second survey network also results in some overlapping survey point locations. These duplicate data points refer to data points with two or more data points at a single coordinate location. For subsequent data processing, which requires only one data point per coordinate location, it is necessary to select duplicate data points. Therefore, data merging is performed. Merging includes summing and averaging the data from the first and second survey networks at their overlapping locations, recording this as the field value of the overlapping survey point; and summing and averaging the data from the two orthogonal survey networks at their overlapping locations in the second survey network, recording this as the field value of the intersection point.

[0056] Furthermore, regarding interference value correction, that is, interference correction is performed on the distortion of the field value of the measuring point affected by the newly added interference source in the collected data. In this embodiment, since the second measuring network is a large scale and the first measuring network is a small scale, the interference source mainly affects the field value of the measuring point of the second measuring network. The interference source is a human interference source, such as electromagnetic interference.

[0057] In this embodiment, the gravity and magnetic data collected in the field are irregular, for example, when the terrain is very complex or when it is impossible to measure at certain points, resulting in an irregular distribution of actual measuring points. However, further processing of the gravity and magnetic data requires the data to be distributed according to a grid. Therefore, data gridding is used to interpolate the measured data, and the values ​​at the grid nodes are calculated using the values ​​at the irregularly distributed interpolation nodes to meet the needs of data processing. Thus, by independently gridding the gravity and magnetic data of the first and second measurement networks to generate first and second gridded data with the same grid spacing, the irregularly spaced data can be converted into regularly spaced data using an interpolation algorithm. The two gridded data are compared, and interference points in the second gridded data are identified in conjunction with the interference source. The data corresponding to the interference points is replaced with the corresponding data in the first gridded data to complete the interference value correction and generate a composite data volume.

[0058] In some preferred embodiments of the present invention, the method of the present invention may further include:

[0059] S3. Perform data gridding on the composite data volume to generate a gridded data volume;

[0060] S4. Perform displacement direction filtering correction on the gridded data volume to generate the first correction data;

[0061] Furthermore, S3 performs data gridding on the composite data volume to generate a gridded data volume, which is a data file after interference correction of high-density composite acquisition data, and serves as the gridded data volume for subsequent correction processing.

[0062] Further, S4 performs displacement direction filtering correction on the gridded data volume to generate first corrected data. This includes analyzing and verifying the magnetic anomaly map of the gridded data volume to obtain error data; performing displacement direction filtering on the magnetic anomaly data of the survey line segment containing error data distributed along the survey line direction; performing directional filtering on the grid data with the survey line direction as the directional filtering orientation (northeast, northwest); taking the survey line segment as the displacement region; taking the directional filtering displacement coefficient within the displacement region as 1.0 and the directional filtering displacement coefficient outside the region as 0; and transitioning the directional filtering displacement coefficient between the region and the outside with a coefficient of 1.0 to 0.0; multiplying the directional filtering data by the directional filtering displacement coefficient to obtain the displacement direction filtering corrected data value; and removing the filtered corrected data value from the gridded data volume to generate the first corrected data. This first corrected data is the (gravity)magnetic data after displacement direction filtering and denoising processing.

[0063] In some other preferred embodiments of the present invention, the method of the present invention may further include:

[0064] S5. Perform magnetic anomaly displacement correction on the first correction data to generate a high-density gravity and magnetic composite data volume. Further, S5 performs magnetic anomaly displacement correction on the first correction data to generate a high-density gravity and magnetic composite data volume. Magnetic anomaly displacement correction includes: gridding the elevation data in the gravity and magnetic data to obtain an elevation data grid; processing the magnetic field values ​​in the gravity and magnetic data based on the International Geomagnetic Reference Field to obtain magnetic anomaly data; gridding the magnetic anomaly data to generate a magnetic anomaly grid; superimposing the contour lines in the elevation data grid with the contour lines in the magnetic anomaly grid; delineating the existing magnetic topographic influence area and determining the magnetic displacement coefficient; the magnetic displacement coefficient in the magnetic topographic influence area is set to 1, and the magnetic displacement coefficient in other areas is set to 0; using a field value extending 1000m above as a background field value to obtain a residual magnetic anomaly; multiplying the magnetic displacement coefficient by the residual magnetic anomaly to generate a magnetic topographic related displacement correction value; and subtracting the topographic magnetic related displacement correction value from the first correction data to generate a high-density gravity and magnetic composite data volume.

[0065] The gravity and magnetic composite acquisition and processing method provided by this invention eliminates interference through composite acquisition and achieves high-density gravity and magnetic data acquisition and processing by combining multiple corrections. It is suitable for exploration of complex oil and gas targets and has important practical and promotional application value in gravity and magnetic exploration of complex oil and gas targets. It is a further extension of high-precision gravity and magnetic exploration technology and an important innovation of existing gravity and magnetic acquisition and processing methods, providing an important new method for gravity and magnetic data acquisition and processing.

[0066] The method of this invention performs multiple corrections on gravity and magnetic data, including displacement direction filtering correction and magnetic anomaly correction. It provides targeted processing for error data and sets displacement coefficient correction values ​​according to regional characteristics, effectively avoiding erroneous corrections caused by using conventional correction methods that apply a one-size-fits-all approach to all data.

[0067] The specific embodiments of the present invention are further described below with reference to specific examples.

[0068] It should be noted that the steps in each embodiment of the above-mentioned high-density gravity and magnetic force composite acquisition and processing method can be interleaved, substituted, added, or deleted. Therefore, these reasonable permutations and combinations of the high-density gravity and magnetic force composite acquisition and processing method should also fall within the protection scope of this invention, and the protection scope of this invention should not be limited to the embodiments.

[0069] Based on the above objectives, a second aspect of the present invention provides a high-density gravity and magnetic force composite acquisition and processing device. Figure 2 The diagram shown is a schematic representation of an embodiment of the high-density gravity and magnetic force composite acquisition and processing device provided by the present invention. Figure 2 As shown, the high-density gravity and magnetic force composite acquisition and processing device of this invention includes the following modules:

[0070] The measurement network construction module 011 is configured to establish a first measurement network and a second measurement network in the measurement area. Both the first and second measurement networks are gravity and magnetic field measurement networks. The scale and measurement network density of the first and second measurement networks are different. The second measurement network includes two measurement networks with orthogonal measurement lines.

[0071] Data acquisition module 012 is configured to acquire gravity and magnetic data in the first and second measurement networks respectively. The gravity and magnetic data include the vertical coordinate, horizontal coordinate, elevation, gravity field value, and magnetic field value of the measurement point.

[0072] The first processing module 013 is configured to perform preliminary processing on gravity and magnetic data. The preliminary processing includes preliminary merging and interference value correction. The data from the first and second measurement networks after the preliminary processing are merged again to generate a composite data volume.

[0073] In a preferred embodiment of the present invention, the device further includes the following modules:

[0074] The second processing module is configured to perform data gridding on the composite data volume to generate a gridded data volume.

[0075] The first correction module is configured to perform directional filtering correction on the gridded data volume to generate first correction data. The first correction module obtains error data by analyzing and verifying the gravity and magnetic anomaly map of the gridded data volume, performs directional filtering on the error data along the survey line direction to obtain directional filtering data, multiplies the directional filtering data by the directional filtering displacement coefficient to generate filtered correction data values, and removes the filtered correction data values ​​from the gridded data volume to generate the first correction data. Specifically, the survey line is considered the displacement region; the directional filtering displacement coefficient within the displacement region is 1, and the directional filtering displacement coefficient outside the displacement region is 0, with the directional filtering displacement coefficient transitioning between 1 and 0 between the two regions.

[0076] The second correction module is configured to acquire elevation and magnetic anomaly data based on gravity and magnetic data. It further calculates residual magnetic anomaly data based on the magnetic anomaly data, delineates the magnetic topographic influence area by overlaying the elevation and magnetic anomaly data, and determines the magnetic displacement coefficient. The magnetic displacement coefficient is multiplied by the residual magnetic anomaly to generate a magnetic topographic-related displacement correction value. This value is then subtracted from the first correction data to generate a high-density gravity and magnetic composite data volume. The magnetic displacement coefficient is set to 1 in the magnetic topographic influence area and 0 in other areas.

[0077] Finally, it should be noted that those skilled in the art will understand that all or part of the processes in the above embodiments can be implemented by a computer program instructing related hardware. The program for the high-density gravity and magnetic force composite acquisition and processing method can be stored in a computer-readable storage medium. When executed, the program can include the processes of the embodiments of the above methods. The storage medium for the program can be a magnetic disk, optical disk, read-only memory (ROM), or random access memory (RAM), etc. The above computer program embodiments can achieve the same or similar effects as any of the corresponding foregoing method embodiments.

[0078] Furthermore, the method disclosed in the embodiments of the present invention can also be implemented as a computer program executed by a processor, which may be stored in a computer-readable storage medium. When the computer program is executed by the processor, it performs the functions defined in the method disclosed in the embodiments of the present invention.

[0079] Furthermore, the above-described method steps and system units can also be implemented using a controller and a computer-readable storage medium for storing a computer program that enables the controller to perform the functions of the above-described steps or units.

[0080] Those skilled in the art will also understand that the various exemplary logic blocks, modules, circuits, and algorithm steps described in conjunction with the disclosure herein can be implemented as electronic hardware, computer software, or a combination of both. To clearly illustrate this interchangeability between hardware and software, the functionality of various illustrative components, blocks, modules, circuits, and steps has been generally described. Whether this functionality is implemented as software or as hardware depends on the specific application and the design constraints imposed on the system as a whole. Those skilled in the art can implement the functionality in various ways for each specific application, but such implementation decisions should not be construed as departing from the scope of the embodiments disclosed herein.

[0081] In one or more exemplary designs, functionality may be implemented in hardware, software, firmware, or any combination thereof. If implemented in software, the functionality may be stored as one or more instructions or code on or transmitted via a computer-readable medium. Computer-readable media include computer storage media and communication media, including any medium that facilitates the transfer of a computer program from one location to another. Storage media may be any available medium accessible to a general-purpose or special-purpose computer. By way of example, and not limitation, computer-readable media may include RAM, ROM, EEPROM, CD-ROM or other optical disc storage devices, disk storage devices or other magnetic storage devices, or any other medium that may be used to carry or store the required program code in the form of instructions or data structures and is accessible to a general-purpose or special-purpose computer or a general-purpose or special-purpose processor. Furthermore, any connection may be appropriately referred to as computer-readable media. For example, if software is transmitted from a website, server, or other remote source using coaxial cable, fiber optic cable, twisted pair, digital subscriber line (DOL), or wireless technologies such as infrared, radio, and microwave, then the aforementioned coaxial cable, fiber optic cable, twisted pair, DOL, or wireless technologies such as infrared, radio, and microwave are all included in the definition of media. As used herein, disks and optical discs include compact discs (CDs), laser discs, optical discs, digital multifunction discs (DVDs), floppy disks, and Blu-ray discs, wherein disks typically reproduce data magnetically, while optical discs reproduce data optically using lasers. Combinations of the above should also be included within the scope of computer-readable media.

[0082] The above are exemplary embodiments disclosed in this invention. However, it should be noted that various changes and modifications can be made without departing from the scope of the embodiments of this invention as defined by the claims. The functions, steps, and / or actions of the methods according to the disclosed embodiments described herein do not need to be performed in any particular order. Furthermore, although the elements disclosed in the embodiments of this invention may be described or claimed individually, they may be understood as multiple unless explicitly limited to a singular number.

[0083] It should be understood that, as used herein, the singular form “a” is intended to include the plural form as well, unless the context clearly supports an exception. It should also be understood that, as used herein, “and / or” refers to any and all possible combinations of one or more of the associated listed items.

[0084] The embodiment numbers disclosed in the above embodiments of the present invention are for descriptive purposes only and do not represent the superiority or inferiority of the embodiments.

[0085] Those skilled in the art will understand that all or part of the steps of the above embodiments can be implemented by hardware or by a program instructing related hardware. The program can be stored in a computer-readable storage medium, such as a read-only memory, a disk, or an optical disk.

[0086] 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 the invention (including the claims) is limited to these examples. Within the framework of the invention, technical features of the above embodiments or different embodiments can be combined, and many other variations of different aspects of the invention exist, which are not provided in the details for the sake of brevity. Therefore, any omissions, modifications, equivalent substitutions, improvements, etc., made within the spirit and principles of the invention should be included within the protection scope of the invention.

Claims

1. A high-density gravity and magnetic force composite acquisition and processing method, characterized in that, Includes the following steps: A first measurement network and a second measurement network are established in the measurement area. Both the first measurement network and the second measurement network are gravity and magnetic field measurement networks. The measurement network densities of the first measurement network and the second measurement network are different. The second measurement network includes measurement networks with two orthogonal direction measurement lines. Gravity and magnetic force data are collected in the first and second measurement networks respectively, and the gravity and magnetic force data collected in the first and second measurement networks are preliminarily processed. The preliminary processing includes preliminary merging and interference value correction. The data from the first and second measurement networks after the preliminary processing are merged a second time to generate a composite data volume. The composite data volume is then meshed to generate a meshed data volume; The gridded data volume is subjected to displacement direction filtering correction to generate first correction data; The first correction data is subjected to magnetic anomaly displacement correction to generate a high-density gravity and magnetic force composite data volume; The displacement direction filtering correction includes: The gridded data volume is analyzed and verified using gravity and magnetic anomaly diagrams to obtain error data. The error data is then subjected to directional filtering along the survey line direction to obtain directional filtering data. The directional filtering data is multiplied by the directional filtering potential change coefficient to generate a filtering correction data value. The filtering correction data value is then removed from the gridded data volume to generate the first correction data. Taking the survey line as the displacement region, the directional filter displacement coefficient within the displacement region is 1, the directional filter displacement coefficient outside the displacement region is 0, and the directional filter displacement coefficient between the displacement region and the displacement region transitions between 1 and 0. The magnetic anomaly displacement correction includes: Elevation data and magnetic anomaly data are obtained based on gravity and magnetic data. Residual magnetic anomaly data are further calculated based on the magnetic anomaly data. The magnetic topographic influence area is delineated and the magnetic displacement coefficient is determined by overlaying the elevation data and the magnetic anomaly data. The magnetic displacement coefficient is multiplied by the residual magnetic anomaly to generate a magnetic topographic related displacement correction value. The magnetic topographic related displacement correction value is subtracted from the first correction data to generate a high-density gravity and magnetic composite data volume.

2. The high-density gravity and magnetic force composite acquisition and processing method according to claim 1, characterized in that, The magnetic displacement coefficient in the area affected by the magnetic topography is 1, while the magnetic displacement coefficient in other areas is 0.

3. The high-density gravity and magnetic force composite acquisition and processing method according to claim 1, characterized in that, The first and second survey networks have different azimuth angles.

4. The high-density gravity and magnetic force composite acquisition and processing method according to claim 1, characterized in that, The density of the second measurement network is higher than that of the first measurement network.

5. The high-density gravity and magnetic force composite acquisition and processing method according to claim 4, characterized in that, The first survey network has a survey network density of 500m × 500m, and the second survey network has a survey network density of 200m × 500m.

6. The high-density gravity and magnetic force composite acquisition and processing method according to claim 1, characterized in that, The preliminary merger includes: At the overlapping position of two orthogonal measurement networks in the second measurement network, the data of the two orthogonal measurement networks are summed and averaged, and recorded as the data at the overlapping position; At the overlapping position of the first and second measurement networks, the data from the first and second measurement networks are summed and averaged, and this average is recorded as the data at that overlapping position.

7. The high-density gravity and magnetic force composite acquisition and processing method according to claim 1, characterized in that, The interference value correction includes: The gravity and magnetic data of the first and second measurement networks are respectively processed into grids to generate first gridded data and second gridded data. The first gridded data and the second gridded data use the same grid spacing to facilitate interpolation. Interference points in the second gridded data are identified in conjunction with interference sources, and the data corresponding to the interference points are replaced with the corresponding data in the first gridded data.

8. The high-density gravity and magnetic force composite acquisition and processing method according to claim 4, characterized in that, The source of interference is human-made interference.

9. A high-density gravity and magnetic force composite acquisition and processing device, characterized in that, include: A measurement network construction module is configured to establish a first measurement network and a second measurement network in the measurement area. Both the first measurement network and the second measurement network are gravity and magnetic force general point measurement networks. The measurement network densities of the first measurement network and the second measurement network are different. The second measurement network includes measurement networks with two orthogonal direction measurement lines. The data acquisition module collects gravity and magnetic force data in the first and second measurement networks, respectively. The gravity and magnetic force data includes the vertical coordinate, horizontal coordinate, elevation, gravity field value, and magnetic field value of the measurement point. The first processing module is configured to perform preliminary processing on the gravity and magnetic data collected in the first and second measurement networks. The preliminary processing includes preliminary merging and interference value correction. The data from the first and second measurement networks after the preliminary processing are merged a second time to generate a composite data volume. Also includes: The second processing module is configured to perform data gridding on the composite data volume to generate a gridded data volume. A first correction module is configured to perform displacement direction filtering correction on the gridded data volume to generate first correction data. The second correction module is configured to perform magnetic anomaly displacement correction on the first correction data to generate a high-density gravity and magnetic force composite data volume. The first correction module obtains error data by performing gravity and magnetic anomaly diagram analysis on the gridded data volume, performs directional filtering on the error data along the survey line direction to obtain directional filtering data, multiplies the directional filtering data with the directional filtering potential change coefficient to generate a filtered correction data value, and removes the filtered correction data value from the gridded data volume to generate the first correction data. Taking the survey line as the displacement region, the directional filter displacement coefficient within the displacement region is 1, the directional filter displacement coefficient outside the displacement region is 0, and the directional filter displacement coefficient between the displacement region and the displacement region transitions between 1 and 0. The second correction module acquires elevation data and magnetic anomaly data based on gravity and magnetic data. It further calculates residual magnetic anomaly data based on the magnetic anomaly data. By superimposing the elevation data and the magnetic anomaly data, it delineates the magnetic topographic influence area and determines the magnetic displacement coefficient. The magnetic displacement coefficient is multiplied by the residual magnetic anomaly to generate a magnetic topographic-related displacement correction value. The magnetic topographic-related displacement correction value is subtracted from the first correction data to generate a high-density gravity and magnetic composite data volume.

10. The high-density gravity and magnetic force composite acquisition and processing device according to claim 9, characterized in that, The magnetic displacement coefficient in the area affected by the magnetic topography is 1, while the magnetic displacement coefficient in other areas is 0.

Images