A method, system and device for calculating the three-dimensional shape of an iron ore body

Abstract

The present invention discloses a method, system and equipment for calculating the three-dimensional morphology of an iron ore body. The method comprises: obtaining magnetic anomaly data and a magnetic anomaly area of ​​a known hidden magnetite body in a target area; obtaining a first top surface burial depth of the magnetite body based on the magnetic anomaly area or through drilling data and obtaining a three-dimensional distribution of magnetization intensity of the magnetic anomaly area; delineating a range of magnetic bodies having a magnetization intensity less than a preset magnetization intensity according to the first top surface burial depth of the magnetite body and the three-dimensional distribution of magnetization intensity; calculating a minimum value of the difference between the second top surface depth and the first top surface burial depth according to the second top surface depth of the iron ore body calculated according to the preset magnetization intensity; extracting a magnetization intensity distribution range inverted greater than the minimum value of the depth difference from a magnetization intensity inversion result, and calculating the morphological value of the iron ore body according to the magnetization intensity distribution range.

Description

Technical Field

[0001] The present invention relates to the field of geophysical exploration, and more specifically, to a method, system and device for calculating the three-dimensional morphology of an iron ore body. Background Art

[0002] Domestic iron ore, especially rich iron ore, falls far short of demand. Magnetite is the primary occurrence of iron ore. As exploration progresses, shallow deposits are well understood, and deep iron ore exploration within covered areas and known deposits holds enormous potential. The exploration and mining of concealed magnetite deposits are becoming increasingly important, requiring an understanding of their three-dimensional morphology. Concealed magnetite often exhibits magnetic characteristics distinct from other geological bodies, making it possible to calculate its three-dimensional morphology using aeromagnetic methods. Conventional inversion methods for determining magnetite morphology are closely tied to the interpreter's experience, resulting in low efficiency, strong ambiguity, and difficulty in determining the 3D morphology of the magnetite ore bodies.

[0003] In view of this, overcoming the technical defects of the above-mentioned prior art is an urgent problem to be solved in this technical field. Summary of the Invention

[0004] In response to the above-mentioned defects or improvement needs of the prior art, the present invention provides a method, system and equipment for calculating the three-dimensional morphology of iron ore bodies. Its purpose is to solve the practical difficulty of calculating the three-dimensional morphology of iron ore bodies by aeromagnetic method, avoid the uncertainty of delineating iron ore bodies by three-dimensional inversion method, and provide important information required for iron ore exploration.

[0005] To achieve the above object, according to one aspect of the present invention, a method for calculating the three-dimensional morphology of an iron ore body is provided, the method comprising:

[0006] Obtain magnetic anomaly data and magnetic anomaly areas of known concealed magnetite bodies in the target area;

[0007] Based on the magnetic anomaly area, the first top surface burial depth of the magnetite body is calculated or obtained through drilling data, and the three-dimensional distribution of magnetization intensity in the magnetic anomaly area is obtained;

[0008] Determining a range of magnetic bodies having a magnetization intensity less than a preset magnetization intensity according to the buried depth of the first top surface of the magnetite body and the three-dimensional distribution of the magnetization intensity;

[0009] a second top surface depth of the iron ore body calculated according to the preset magnetization intensity, and calculating a minimum difference between the second top surface depth and the first top surface burial depth;

[0010] A magnetization intensity distribution range greater than the minimum inversion value of the depth difference is extracted from the magnetization intensity inversion result, and the morphological value of the iron ore body is calculated according to the magnetization intensity distribution range.

[0011] Preferably, the acquisition of magnetic anomaly data of a known hidden magnetite body in the target area includes: the amplitude and shape of the magnetic anomaly caused by the iron ore body.

[0012] Preferably, the method for obtaining a magnetic anomaly area within the range includes:

[0013] Obtain magnetic anomaly areas caused by concealed magnetite bodies in magnetization pole data.

[0014] Preferably, the method for calculating the first top surface burial depth of the magnetite body based on the magnetic anomaly area includes:

[0015] A magnetic anomaly curve is prepared according to the magnetic anomaly data, and a first top surface burial depth is calculated using tangents of characteristic points on the magnetic anomaly curve and the horizontal coordinates of intersections between the tangents;

[0016] Alternatively, the first top surface burial depth is obtained according to known drilling data of the iron ore body.

[0017] Preferably, the magnetic anomaly curve is drawn with 5 tangents and 4 intersection points according to the tangent method.

[0018] Preferably, the method for obtaining the three-dimensional distribution of magnetization intensity in the magnetic anomaly region includes:

[0019] Establishing a depth-weighted focusing inversion objective function based on the magnetic anomaly region;

[0020] The three-dimensional distribution of magnetization intensity in the magnetic anomaly area is obtained based on the depth-weighted focusing inversion objective function.

[0021] Preferably, the method for calculating the first top surface burial depth using tangents of characteristic points on the magnetic anomaly curve and the horizontal coordinates of intersections between the tangents includes:

[0022] h=K(x1-x2+x3-x4) / 2=K(b1+b2) / 2

[0023] Wherein: h is the buried depth of the first top surface, K is the correction coefficient, x is the abscissa of the intersection between the tangent lines, the order of the abscissas from right to left is x1, x2, x3, x4, b1 is the difference between x1 and x2, and b2 is the difference between x3 and x4;

[0024] Alternatively, the actual buried depth of the iron ore body is obtained through drilling data as the first top surface buried depth.

[0025] Preferably, the method for defining the range of magnetic bodies having a magnetization intensity less than a preset magnetization intensity based on the buried depth of the first top surface of the magnetite body and the three-dimensional distribution of magnetization intensity comprises:

[0026] Assume that there are n first top surface buried depths of the magnetite body calculated based on the magnetic anomaly area, where the first top surface buried depth of point i is h i The vertical projection depth of the top surface of the magnetite body at the preset magnetization intensity z value is z. i The vertical projection depth of the top surface of the magnetite body at point i in other magnetization intensity g values ​​is g i , the preset magnetization intensity z value should satisfy:

[0027] As a further improvement and supplement to the above solution, the present invention also includes the following additional technical features.

[0028] According to another aspect of the present invention, a system for calculating the three-dimensional morphology of an iron ore body is provided, the system comprising:

[0029] A data module is used to obtain magnetic anomaly data and magnetic anomaly areas of known hidden magnetite bodies in the target area;

[0030] a calculation module, configured to calculate a first top surface burial depth of the magnetite body based on the magnetic anomaly region and obtain a three-dimensional distribution of magnetization intensity in the magnetic anomaly region;

[0031] a delineation module, configured to delineate a range of magnetic bodies having a magnetization intensity less than a preset magnetization intensity according to the buried depth of the first top surface of the magnetite body and the three-dimensional distribution of the magnetization intensity;

[0032] a statistical module, configured to calculate a minimum value of a difference between a second top surface depth of the iron ore body obtained by calculating the second top surface depth and the first top surface burial depth according to the preset magnetization intensity;

[0033] The extraction module is used to extract the magnetization intensity distribution range inverted from the magnetization intensity inversion result, which is greater than the minimum depth difference, and calculate the morphological value of the iron ore body according to the magnetization intensity distribution range.

[0034] According to another aspect of the present invention, there is provided an apparatus for calculating the three-dimensional shape of an iron ore body, the apparatus comprising:

[0035] one or more processors;

[0036] A storage device is used to store one or more programs. When the one or more programs are executed by the one or more processors, the one or more processors implement the method for calculating the three-dimensional morphology of the iron ore body as described in the first aspect.

[0037] In general, the above technical solutions conceived by the present invention have the following beneficial effects compared with the prior art:

[0038] Conventional 3D inversion methods for delineating magnetic bodies often rely on subjective judgment, requiring multiple attempts, resulting in considerable uncertainty and inefficiency. This invention, based on computer-generated 3D inversion technology, magnetic body tangent calculation techniques, and logical operations, allows for automated computer calculation of iron ore morphology, resulting in greater efficiency and certainty compared to conventional methods.

[0039] The present invention manually adds actual ore body buried depth points as ore body shape constraints, accurately depicts the shape of the ore body outside the borehole under the borehole constraint, introduces more evidence than the existing technology, and increases credibility. BRIEF DESCRIPTION OF THE DRAWINGS

[0040] To more clearly illustrate the technical solutions of the embodiments of the present invention, the following briefly introduces the drawings required for use in the embodiments of the present invention. Obviously, the drawings described below are only some embodiments of the present invention. Those skilled in the art can also derive other drawings based on these drawings without inventive effort.

[0041] Figure 1 This is a schematic flow chart of a method for calculating the three-dimensional morphology of an iron ore body provided in the first embodiment;

[0042] Figure 2 is a schematic diagram of the tangent method for the magnetic anomaly curve in the first embodiment;

[0043] Figure 3 Schematic diagram of the top surface of the magnetite body defined by three different preset magnetization intensities in the first embodiment;

[0044] Figure 4 Schematic diagram of a cross-sectional curve for calculating the top buried depth of a magnetite body and the isosurface of a magnetite body using the tangent method in the first embodiment;

[0045] Figure 5 This is a schematic diagram of the system for calculating the three-dimensional morphology of the iron ore body in the second embodiment;

[0046] Figure 6 This is a schematic diagram of the equipment for calculating the three-dimensional shape of the iron ore body in the third embodiment. DETAILED DESCRIPTION

[0047] In order to make the objectives, technical solutions and advantages of the present invention more clearly understood, the present invention is further described in detail below with reference to the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are merely for the purpose of explaining the present invention and are not intended to limit the present invention. In addition, the technical features involved in the various embodiments of the present invention described below may be combined with each other as long as they do not conflict with each other.

[0048] Example 1

[0049] This embodiment provides a method for calculating the three-dimensional shape of an iron ore body, the method comprising the following steps: Figure 1 As shown:

[0050] S101: Acquire magnetic anomaly data and magnetic anomaly areas of known hidden magnetite bodies in the target area.

[0051] The magnetic anomaly data includes the amplitude and shape of the magnetic anomaly caused by the iron ore body. The range of the hidden magnetite body in the target area is determined based on the magnetic anomaly data, and the magnetic anomaly area caused by the hidden magnetite body is obtained based on the aeromagnetic pole data.

[0052] S102: Calculating the first top surface burial depth of the magnetite body based on the magnetic anomaly region or obtaining drilling data and obtaining a three-dimensional distribution of magnetization intensity in the magnetic anomaly region.

[0053] Establishing a depth-weighted focusing inversion objective function based on the magnetic anomaly region;

[0054] The objective function of the magnetic three-dimensional inversion is as follows:

[0055] ||Wd(G*m–d)||2+μ*||Wm(m–mref)||=min

[0056] Where ||Wd(G*m–d)||2 is the data fitting function, ||Wm(m–mref)|| is the model objective function, Wd is the data space weighting matrix, Wm is the model space weighting matrix, μ is the Lagrangian operator or regularization factor, G is the magnetic forward and inversion calculation kernel function, and m is the physical property model vector.

[0057] The three-dimensional distribution of magnetization intensity in the region is obtained based on the depth-weighted focusing inversion objective function. The spatial morphology of the ferromagnetic body (actual iron ore body) can be delineated using the magnetization intensity z value and the three-dimensional distribution of magnetization intensity.

[0058] S103: Determine a range of magnetic bodies having a magnetization intensity less than a preset magnetization intensity according to the buried depth of the first top surface of the magnetite body and the three-dimensional distribution of the magnetization intensity.

[0059] The method for calculating the first top surface burial depth using the tangent lines of the characteristic points on the magnetic anomaly curve and the horizontal coordinates of the intersection points between the tangent lines includes:

[0060] h=K(x1-x2+x3-x4) / 2=K(b1+b2) / 2

[0061] Wherein: h is the buried depth of the first top surface, K is the correction coefficient, x is the abscissa of the intersection between the tangent lines, the order of the abscissas from right to left is x1, x2, x3, x4, b1 is the difference between x1 and x2, and b2 is the difference between x3 and x4;

[0062] Alternatively, the actual buried depth of the iron ore body is obtained through drilling data as the first top surface buried depth.

[0063] The tangent method is a method of calculating the minimum burial depth (top burial depth) or center burial depth, width and other basic parameters of the magnetic body by using the relationship between the coordinates of the intersection points of the tangents passing through some characteristic points (such as extreme points, inflection points, etc.) on the magnetic anomaly curve, thereby determining the occurrence elements of the magnetic body. This method is applicable to any Za (vertical magnetic anomaly) and ΔT (total magnetic field anomaly) curves. When calculating the burial depth of a magnetic body, first make 5 tangents to the Za or ΔT curves. Three of the horizontal tangents pass through the maximum and minimum points respectively, and the other two tangents pass through the two inflection points of the curve. The 5 tangents intersect at 4 points, and their horizontal coordinates are x1, x2, x3, and x4 respectively. The magnetic anomaly tangent method is shown as follows. Figure 2 shown.

[0064] S104: Calculating a minimum value of a difference between a second top surface depth of the iron ore body obtained by calculating the second top surface depth and the first top surface burial depth according to the preset magnetization intensity.

[0065] The method for defining the range of magnetic bodies having a magnetization intensity less than a preset magnetization intensity according to the buried depth of the first top surface of the magnetite body and the three-dimensional distribution of the magnetization intensity includes:

[0066] Assume that there are n first top surface buried depths of the magnetite body calculated based on the magnetic anomaly area, where the first top surface buried depth of point i is h i The vertical projection depth of the top surface of the magnetite body at the preset magnetization intensity z value is z. i The vertical projection depth of the top surface of the magnetite body at point i in other magnetization intensity g values ​​is g i , the preset magnetization intensity z value should satisfy:

[0067] S105: extracting a magnetization intensity distribution range greater than the minimum inverted depth difference from the magnetization intensity inversion result, and calculating the morphological value of the iron ore body according to the magnetization intensity distribution range.

[0068] In combination with this embodiment, there is also a preferred implementation scheme. Specifically, obtaining magnetic anomaly data of a known hidden magnetite body in the target area includes: the amplitude and shape of the magnetic anomaly caused by the iron ore body.

[0069] In conjunction with this embodiment, there is also a preferred implementation scheme. Specifically, the method for obtaining the magnetic anomaly area within the range includes:

[0070] Obtain magnetic anomaly areas caused by concealed magnetite bodies in magnetization pole data.

[0071] In conjunction with this embodiment, there is also a preferred implementation scheme. Specifically, the method for calculating the buried depth of the first top surface of the magnetite body based on the magnetic anomaly area includes:

[0072] A magnetic anomaly curve is prepared according to the magnetic anomaly data, and a first top surface burial depth is calculated using tangents of characteristic points on the magnetic anomaly curve and the horizontal coordinates of intersections between the tangents;

[0073] Alternatively, the first top surface burial depth is obtained according to known drilling data of the iron ore body.

[0074] In combination with this embodiment, there is also a preferred implementation scheme. Specifically, the magnetic anomaly curve is drawn with 5 tangents and 4 intersection points according to the tangent method.

[0075] In conjunction with this embodiment, there is also a preferred implementation scheme. Specifically, the method for obtaining the three-dimensional distribution of magnetization intensity in the magnetic anomaly area includes:

[0076] Establishing a depth-weighted focusing inversion objective function based on the magnetic anomaly region;

[0077] The three-dimensional distribution of magnetization intensity in the magnetic anomaly area is obtained based on the depth-weighted focusing inversion objective function.

[0078] In conjunction with this embodiment, there is also a preferred implementation scheme. Specifically, the method of calculating the first top surface burial depth using the tangent lines of the characteristic points on the magnetic anomaly curve and the horizontal coordinates of the intersection points between the tangent lines includes:

[0079] h=K(x1-x2+x3-x4) / 2=K(b1+b2) / 2

[0080] Wherein: h is the buried depth of the first top surface, K is the correction coefficient, x is the abscissa of the intersection between the tangent lines, the order of the abscissas from right to left is x1, x2, x3, x4, b1 is the difference between x1 and x2, and b2 is the difference between x3 and x4;

[0081] Alternatively, the actual buried depth of the iron ore body is obtained through drilling data as the first top surface buried depth.

[0082] In conjunction with this embodiment, there is also a preferred implementation scheme. Specifically, the method of delineating the range of magnetic bodies having a magnetization intensity less than a preset magnetization intensity based on the buried depth of the first top surface of the magnetite body and the three-dimensional distribution of the magnetization intensity includes:

[0083] Assume that there are n first top surface buried depths of the magnetite body calculated based on the magnetic anomaly area, where the first top surface buried depth of point i is h i The vertical projection depth of the top surface of the magnetite body at the preset magnetization intensity z value is z. i The vertical projection depth of the top surface of the magnetite body at point i in other magnetization intensity g values ​​is g i , the preset magnetization intensity z value should satisfy:

[0084] This formula can determine whether the depth of the top surface of the magnetite body delineated by using the preset magnetization intensity z value is closer to the depth point calculated by the tangent method than the depth of the top surface of the magnetite body delineated by other magnetization intensity values.

[0085] like Figure 3 and Figure 4 As shown, M1-M3 are the top surfaces of three magnetite bodies delineated by different parameters. The blue dots represent the buried depth of the top surface of the magnetite body calculated by the tangent method. The curves represent the cross-sectional curves of the isosurfaces of the magnetite body delineated by different magnetic susceptibilities. Figure 4 The vertical sections of the top surfaces of the three magnetite bodies with different colors are represented as Figure 3 The three curves in Figure 4 The blue points are the depth values ​​at different positions calculated by the tangent method on the cross section. Figure 3 It can be seen that the M2 curve calculated using a certain parameter is closer to the top surface burial depth value calculated by the tangent method than the M1 and M3 curves, and then the magnetization intensity z value used by the M2 curve is used to calculate the morphological value of the iron ore body.

[0086] This embodiment utilizes computer-generated 3D inversion technology, magnetic tangent calculation techniques, and logical operations to automatically calculate the morphology of iron ore bodies, resulting in greater efficiency and certainty than conventional methods. Furthermore, this embodiment manually incorporates actual ore body depths as morphological constraints, accurately depicting the morphology of the ore body outside the borehole within the constraints of the borehole. This incorporates more evidence and enhances credibility compared to existing techniques.

[0087] Example 2:

[0088] This embodiment 2 provides a system for calculating the three-dimensional shape of an iron ore body, such as Figure 5 As shown, the system includes:

[0089] A data module is used to obtain magnetic anomaly data and magnetic anomaly areas of known hidden magnetite bodies in the target area;

[0090] a calculation module, configured to calculate a first top surface burial depth of the magnetite body based on the magnetic anomaly region and obtain a three-dimensional distribution of magnetization intensity in the magnetic anomaly region;

[0091] a delineation module, configured to delineate a range of magnetic bodies having a magnetization intensity less than a preset magnetization intensity according to the buried depth of the first top surface of the magnetite body and the three-dimensional distribution of the magnetization intensity;

[0092] a statistical module, configured to calculate a minimum value of a difference between a second top surface depth of the iron ore body obtained by calculating the second top surface depth and the first top surface burial depth according to the preset magnetization intensity;

[0093] The extraction module is used to extract the magnetization intensity distribution range inverted from the magnetization intensity inversion result, which is greater than the minimum depth difference, and calculate the morphological value of the iron ore body according to the magnetization intensity distribution range.

[0094] Example 3:

[0095] A device for calculating the three-dimensional shape of an iron ore body, such as Figure 6 As shown, the equipment includes:

[0096] one or more processors;

[0097] A storage device is used to store one or more programs. When the one or more programs are executed by one or more processors, the one or more processors implement the method for calculating the three-dimensional morphology of an iron ore body as described in any one of the first embodiments.

[0098] Figure 6 This is a schematic diagram of the structure of the equipment for calculating the three-dimensional shape of the iron ore body provided in the third embodiment. Figure 6 A block diagram of an exemplary device for calculating the three-dimensional shape of an iron ore body suitable for implementing an embodiment of the present invention is shown. Figure 6 The device for calculating the three-dimensional shape of an iron ore body shown is merely an example and should not limit the functions and scope of use of the embodiments of the present invention.

[0099] like Figure 6 As shown, the device for calculating the three-dimensional shape of an iron ore body is presented as a general device. Components of the device for calculating the three-dimensional shape of an iron ore body may include, but are not limited to, one or more processors or processing units, memory, and a bus connecting different system components (including the memory and processing units).

[0100] The term "bus" refers to one or more of several types of bus structures, including a memory bus or memory controller, a peripheral bus, an accelerated graphics port, a processor, or a local bus using any of a variety of bus architectures. Examples of these architectures include, but are not limited to, the Industry Standard Architecture (ISA) bus, the Micro Channel Architecture (MAC) bus, the Enhanced ISA bus, the Video Electronics Standards Association (VESA) local bus, and the Peripheral Component Interconnect (PCI) bus.

[0101] The device for calculating the three-dimensional shape of an iron ore body typically includes a variety of computer system readable media. These media can be any available media that can be accessed by the device for correcting the intelligent well logging interpretation model, including volatile and non-volatile media, removable and non-removable media.

[0102] The memory may include computer system readable media in the form of volatile memory, such as random access memory (RAM) 30 and / or cache memory. The apparatus for calculating the three-dimensional morphology of the iron ore body may further include other removable / non-removable, volatile / non-volatile computer system storage media. By way of example only, the storage system may be used to read and write non-removable, non-volatile magnetic media ( Figure 6 Not shown, often called a "hard drive"). Although Figure 6 Not shown, a disk drive for reading and writing to a removable non-volatile disk (e.g., a "floppy disk"), and an optical drive for reading and writing to a removable non-volatile optical disk (e.g., a CD-ROM, DVD-ROM, or other optical medium) may be provided. In these cases, each drive may be connected to the bus via one or more data medium interfaces. The memory may include at least one program product having a set (e.g., at least one) of program modules configured to perform the functions of various embodiments of the present invention.

[0103] A program / utility having a set (at least one) of program modules, which may be stored, for example, in a memory, includes, but is not limited to, an operating system, one or more application programs, other program modules, and program data, each of which, or some combination thereof, may include an implementation of a network environment. The program modules generally implement the functions and / or methods of the embodiments described herein.

[0104] The device for calculating the three-dimensional shape of the iron ore body can also communicate with one or more external devices (such as keyboards, pointing devices, displays, etc.), and can also communicate with one or more devices that enable users to interact with the device for calculating the three-dimensional shape of the iron ore body, and / or communicate with any device that enables the device for calculating the three-dimensional shape of the iron ore body to communicate with one or more other devices (such as network cards, modems, etc.). Such communication can be carried out through an input / output (I / O) interface. In addition, the device for correcting the intelligent logging interpretation model can also communicate with one or more networks (such as local area networks (LANs), wide area networks (WANs) and / or public networks, such as the Internet) through a network adapter. Figure 6As shown, the network adapter communicates with other modules of the apparatus for calculating the three-dimensional morphology of an iron ore body via a bus. It should be understood that, although not shown in the figures, other hardware and / or software modules may be used in conjunction with the apparatus for calculating the three-dimensional morphology of an iron ore body, including but not limited to: microcode, device drivers, redundant processing units, external disk drive arrays, RAID systems, tape drives, and data backup storage systems.

[0105] The processing unit executes various functional applications and data processing by running programs stored in the memory, such as implementing the method for calculating the three-dimensional morphology of an iron ore body provided in any embodiment of the present invention. Specifically, the method comprises: obtaining magnetic anomaly data and a magnetic anomaly region of a known hidden magnetite body in a target area; calculating a first top surface burial depth of the magnetite body based on the magnetic anomaly region or using drilling data, and obtaining a three-dimensional distribution of magnetization intensity in the magnetic anomaly region; delineating a range of magnetic bodies having a magnetization intensity less than a preset magnetization intensity based on the first top surface burial depth of the magnetite body and the three-dimensional distribution of magnetization intensity; calculating a minimum difference between the second top surface depth and the first top surface burial depth obtained from the second top surface depth of the iron ore body according to the preset magnetization intensity; extracting a magnetization intensity distribution range inverted greater than the minimum difference in depth from the magnetization intensity inversion result, and calculating the morphological value of the iron ore body according to the magnetization intensity distribution range.

[0106] The above description is only a preferred embodiment of the present invention and is not intended to limit the present invention. Any modifications, equivalent substitutions and improvements made within the spirit and principles of the present invention should be included in the scope of protection of the present invention.

Claims

1. A method for calculating the three-dimensional morphology of an iron ore body, characterized in that the method include: Obtain magnetic anomaly data and magnetic anomaly areas of known concealed magnetite bodies in the target area; Based on the magnetic anomaly area, the first top surface burial depth of the magnetite body is calculated or obtained through drilling data, and the three-dimensional distribution of magnetization intensity in the magnetic anomaly area is obtained; Determining a range of magnetic bodies having a magnetization intensity less than a preset magnetization intensity according to the buried depth of the first top surface of the magnetite body and the three-dimensional distribution of the magnetization intensity; a second top surface depth of the iron ore body calculated according to the preset magnetization intensity, and calculating a minimum difference between the second top surface depth and the first top surface burial depth; A magnetization intensity distribution range greater than the minimum inversion value of the depth difference is extracted from the magnetization intensity inversion result, and the morphological value of the iron ore body is calculated according to the magnetization intensity distribution range.

2. The method for calculating the three-dimensional morphology of an iron ore body according to claim 1, wherein: The acquisition of magnetic anomaly data of a known hidden magnetite body in the target area includes: the amplitude and shape of the magnetic anomaly caused by the iron ore body.

3. The method for calculating the three-dimensional morphology of an iron ore body according to claim 2, wherein: Methods for obtaining magnetic anomaly areas include: A magnetic anomaly region caused by the concealed magnetite body is obtained in the magnetization pole data.

4. The method for calculating the three-dimensional morphology of an iron ore body according to claim 3, wherein: The method for calculating the first top surface burial depth of the magnetite body based on the magnetic anomaly area includes: A magnetic anomaly curve is prepared according to the magnetic anomaly data, and a first top surface burial depth is calculated using tangents of characteristic points on the magnetic anomaly curve and the horizontal coordinates of intersections between the tangents; Alternatively, the first top surface burial depth is obtained according to known drilling data of the iron ore body.

5. The method for calculating the three-dimensional morphology of an iron ore body according to claim 4, wherein: The magnetic anomaly curve is drawn with 5 tangent lines and 4 intersection points according to the tangent method.

6. The method for calculating the three-dimensional morphology of an iron ore body according to claim 5, wherein: The method for obtaining the three-dimensional distribution of magnetization intensity in the magnetic anomaly area includes: Establishing a depth-weighted focusing inversion objective function based on the magnetic anomaly region; The three-dimensional distribution of magnetization intensity in the magnetic anomaly area is obtained based on the depth-weighted focusing inversion objective function.

7. The method for calculating the three-dimensional morphology of an iron ore body according to claim 6, wherein: The method for calculating the first top surface burial depth by using the tangent lines of the characteristic points on the magnetic anomaly curve and the horizontal coordinates of the intersection points between the tangent lines includes: h=K(x1-x2+x3-x4) / 2=K(b1+b2) / 2 Wherein: h is the buried depth of the first top surface, K is the correction coefficient, x is the abscissa of the intersection between the tangent lines, the order of the abscissas from right to left is x1, x2, x3, x4, b1 is the difference between x1 and x2, and b2 is the difference between x3 and x4; Alternatively, the actual buried depth of the iron ore body is obtained through drilling data as the first top surface buried depth.

8. The method for calculating the three-dimensional morphology of an iron ore body according to claim 7, wherein: The method for defining a range of magnetic bodies having a magnetization intensity less than a preset magnetization intensity based on the buried depth of the first top surface of the magnetite body and the three-dimensional distribution of magnetization intensity includes: Assume that there are n first top surface buried depths of the magnetite body calculated based on the magnetic anomaly area, where the first top surface buried depth of point i is h i The vertical projection depth of the top surface of the magnetite body at the preset magnetization intensity z value is z. i The vertical projection depth of the top surface of the magnetite body at point i in other magnetization intensity g values ​​is g i , the preset magnetization intensity z value should satisfy:

9. A system for calculating the three-dimensional morphology of an iron ore body, characterized in that: The system includes: A data module is used to obtain magnetic anomaly data and magnetic anomaly areas of known hidden magnetite bodies in the target area; a calculation module, configured to calculate a first top surface burial depth of the magnetite body based on the magnetic anomaly region and obtain a three-dimensional distribution of magnetization intensity in the magnetic anomaly region; a delineation module, configured to delineate a range of magnetic bodies having a magnetization intensity less than a preset magnetization intensity according to the buried depth of the first top surface of the magnetite body and the three-dimensional distribution of the magnetization intensity; a statistical module, configured to calculate a minimum value of a difference between a second top surface depth of the iron ore body obtained by calculating the second top surface depth and the first top surface burial depth according to the preset magnetization intensity; The extraction module is used to extract the magnetization intensity distribution range inverted from the magnetization intensity inversion result, which is greater than the minimum depth difference, and calculate the morphological value of the iron ore body according to the magnetization intensity distribution range.

10. A device for calculating the three-dimensional shape of an iron ore body, characterized in that the device include: one or more processors; A storage device for storing one or more programs, when the one or more programs are executed by the one or more processors, enables the one or more processors to implement the method for calculating the three-dimensional morphology of an iron ore body as described in any one of claims 1-8.

Images