Method, system and equipment for calculating three-dimensional form of iron ore body

Abstract

The invention discloses a method, a system and equipment for calculating a three-dimensional form of an iron ore body. The method comprises the following steps: acquiring magnetic anomaly data and a magnetic anomaly region of a known hidden magnetite body in a target region; the first top surface burial depth of the magnetite body is calculated or obtained through drilling data based on the magnetic anomaly area, and magnetization intensity three-dimensional distribution of the magnetic anomaly area is obtained; delineating a magnetic body range smaller than preset magnetization intensity according to the first top surface burial depth of the magnetite body and the magnetization intensity three-dimensional distribution; calculating the second top surface depth of the iron ore body according to the preset magnetization intensity, and counting the minimum value of the difference between the second top surface depth and the first top surface burial depth; and extracting a magnetization intensity distribution range which is greater than the inversion of the minimum value of the depth difference in the magnetization intensity inversion result, and calculating the form numerical 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, a system and a device for calculating the three-dimensional shape of an iron ore body. Background Art

[0002] Domestic iron ore, especially high-grade iron ore, far from meets the demand. Iron ore mainly occurs in the form of magnetite. With the development of prospecting work, the shallow ore is well understood, and there is great potential for prospecting deep iron ore in the covered area and known ore deposits. Prospecting and mining of hidden magnetite are becoming increasingly important, all of which require understanding the three-dimensional shape of the magnetite body. Hidden magnetite often has magnetic characteristics different from other geological bodies, making it possible to calculate the three-dimensional shape of hidden iron ore by airborne magnetic survey. Conventional inversion methods for determining the shape of magnetite are closely related to the interpretation experience of inversion personnel, with low efficiency, strong non-uniqueness, and it is difficult to obtain the three-dimensional shape of the magnetite ore body.

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

[0004] In view of the above-mentioned defects or improvement requirements of the existing technology, the present invention provides a method, a system and a device for calculating the three-dimensional shape of an iron ore body, aiming to solve the practical difficulty that it is difficult to calculate the three-dimensional shape of an iron ore body by airborne magnetic survey, avoid the uncertainty of delineating an iron ore body by three-dimensional inversion methods, and provide important information required for iron ore exploration.

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

[0006] Obtaining magnetic anomaly data and magnetic anomaly regions of known hidden magnetite ore bodies in a target area;

[0007] Calculating the buried depth of the first top surface of the magnetite ore body based on the magnetic anomaly region or obtaining the three-dimensional magnetization intensity distribution of the magnetic anomaly region through drilling data;

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

[0009] Calculating the depth of the second top surface of the iron ore body according to the preset magnetization intensity, and statistically obtaining the minimum value of the difference between the depth of the second top surface and the buried depth of the first top surface;

[0010] Extracting the magnetization intensity distribution range inverted greater than the minimum value of the difference in depth, and calculating the morphological value of the iron ore body according to the magnetization intensity distribution range.

[0011] Preferably, the obtaining of the magnetic anomaly data of the known concealed magnetite ore bodies 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 the magnetic anomaly area within the range includes:

[0013] Obtaining the magnetic anomaly area caused by the concealed magnetite ore body from the magnetization pole data.

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

[0015] Making a magnetic anomaly curve according to the magnetic anomaly data, and calculating the buried depth of the first top surface by using the tangent of the characteristic points on the magnetic anomaly curve and the abscissa of the intersection points between the tangents;

[0016] Or, obtaining the buried depth of the first top surface according to the known drilling data of the iron ore body.

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

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

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

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

[0021] Preferably, the method for calculating the buried depth of the first top surface by using the tangent of the characteristic points on the magnetic anomaly curve and the abscissa of the intersection points between the tangents includes:

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

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

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

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

[0026] Suppose there are n values for the first top surface buried depth of the magnetite ore body calculated based on the magnetic anomaly region, and the first top surface buried depth at the i-th point is h i , the vertical projection depth of the i-th point on the top surface of the magnetite ore body outlined by the preset magnetization intensity z value is z i , the vertical projection depth of the i-th point on the top surface of the magnetite ore body outlined by 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 further includes the following additional technical features.

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

[0029] A data module for obtaining magnetic anomaly data and magnetic anomaly regions of known buried magnetite ore bodies in a target area;

[0030] A calculation module for calculating the first top surface buried depth of the magnetite ore body based on the magnetic anomaly region and obtaining the three-dimensional magnetization intensity distribution of the magnetic anomaly region;

[0031] A delineation module for delineating the range of magnetic bodies with a magnetization intensity less than the preset magnetization intensity according to the first top surface buried depth of the magnetite ore body and the three-dimensional magnetization intensity distribution;

[0032] A statistics module for calculating the second top surface depth of the iron ore body according to the preset magnetization intensity and statistically obtaining the minimum value of the difference between the second top surface depth and the first top surface buried depth;

[0033] An extraction module for extracting the magnetization intensity distribution range inverted by the magnetization intensity inversion result greater than the minimum value of the depth difference inversion, and calculating the morphological values of the iron ore body according to the magnetization intensity distribution range.

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

[0035] One or more processors;

[0036] A storage device for storing one or more programs, which when executed by the one or more processors cause the one or more processors to implement the method for calculating the three-dimensional shape of an iron ore body as described in the first aspect.

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

[0038] Conventional three-dimensional inversion methods for delineating magnetic bodies often rely mainly on subjective judgment, which requires multiple attempts and has considerable uncertainty and inefficiency. Based on computer three-dimensional automatic inversion technology, the tangent method calculation technology for magnetic bodies, and logical operations, the calculation of the shape of iron ore bodies can be realized by computer automatic operation, which has high efficiency and certainty compared with conventional methods.

[0039] In the present invention, actual ore body burial depth points are manually added as constraints on the ore body shape, and the shape of the ore body outside the drill holes is accurately depicted under drill hole constraints, introducing more evidence than the prior art and increasing the credibility. BRIEF DESCRIPTION OF THE DRAWINGS

[0040] In order to more clearly illustrate the technical solutions of the embodiments of the present invention, the drawings required for use in the embodiments of the present invention will be briefly introduced below. Obviously, the drawings described below are only some embodiments of the present invention, and those of ordinary skill in the art can obtain other drawings based on these drawings without creative efforts.

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

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

[0043] Figure 3 It is a schematic diagram of the top surface of a magnetite ore body delineated by three different preset magnetization intensities in the first embodiment;

[0044] Figure 4 It is a schematic cross-sectional curve diagram of the burial depth of the top surface of a magnetite ore body calculated by the tangent method and the equipotential surface of the magnetite ore body in the first embodiment;

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

[0046] Figure 6 It is a schematic diagram of a device for calculating the three-dimensional shape of an iron ore body in the third embodiment. DETAILED DESCRIPTION OF THE EMBODIMENTS

[0047] In order to make the objectives, technical solutions and advantages of the present invention clearer, the present invention will be further described in detail below with reference to the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are only used to explain the present invention and are not used to limit the present invention. In addition, the technical features involved in the various embodiments of the present invention described below can be combined with each other as long as they do not conflict with each other.

[0048] Embodiment 1

[0049] Embodiment 1 provides a method for calculating the three-dimensional shape of an iron ore body. The method includes the following steps, as Figure 1 shown:

[0050] S101: Obtain the magnetic anomaly data and the magnetic anomaly region of the known concealed magnetite ore body in the target area.

[0051] The magnetic anomaly data includes: the amplitude and shape of the magnetic anomaly caused by the iron ore body. Determine the range of the concealed magnetite ore body in the target area according to the magnetic anomaly data, and then obtain the magnetic anomaly region caused by the concealed magnetite ore body from the aeromagnetic reduction to the pole data.

[0052] S102: Calculate the first top surface burial depth of the magnetite ore body based on the magnetic anomaly region or obtain it through drilling data, and obtain the three-dimensional magnetization intensity distribution of the magnetic anomaly region.

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

[0054] The magnetic three-dimensional inversion objective function 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 Lagrange operator or regularization factor, G is the magnetic forward and inverse calculation kernel function, and m is the physical property model vector.

[0057] Obtain the three-dimensional magnetization intensity distribution of the region based on the depth-weighted focusing inversion objective function. Using the z value of the magnetization intensity and the three-dimensional magnetization intensity distribution, the spatial shape of the strongly magnetic body (the actual iron ore body) can be delineated.

[0058] S103: Delineate the range of the magnetic body with a magnetization intensity less than the preset magnetization intensity according to the first top surface burial depth of the magnetite ore body and the three-dimensional magnetization intensity distribution.

[0059] The method for calculating the first top surface burial depth by using the tangent of the characteristic points on the magnetic anomaly curve and the abscissa of the intersection point between the tangents includes:

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

[0061] Where: h is the buried depth of the first top surface, K is the correction coefficient, x is the abscissa of the intersection point between the tangents, and the abscissas from right to left are x1, x2, x3, x4. b1 is the difference between x1 and x2, and b2 is the difference between x3 and x4;

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

[0063] The tangent method uses 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 to calculate the basic parameters such as the minimum buried depth (top layer buried depth), central buried depth, and width of the magnetic body, so as to determine 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 buried depth of the magnetic body, first make 5 tangents to the curves such as Za or ΔT. Among them, 3 horizontal tangents pass through the maximum point and the minimum point respectively, and the other 2 tangents pass through the 2 inflection points of the curve. The 5 tangents intersect at 4 points, and their abscissas are x1, x2, x3, x4 respectively. The schematic diagram of the magnetic anomaly tangent method is as Figure 2 shown.

[0064] S104: According to the second top surface depth of the iron ore body calculated according to the preset magnetization intensity, count the minimum value of the difference between the second top surface depth and the buried depth of the first top surface.

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

[0066] Assume that there are n calculated buried depths of the first top surface of the magnetite ore body based on the magnetic anomaly region, where the buried depth of the first top surface at the i-th point is h i , the vertical projection depth of the i-th point on the top surface of the magnetite ore body delineated by the preset magnetization intensity z value is z i , the vertical projection depth of the i-th point on the top surface of the magnetite ore body delineated by other magnetization intensity g values is g i , the preset magnetization intensity z value should satisfy:

[0067] S105: Extract the magnetization intensity distribution range inverted by the magnetization intensity greater than the minimum value of the depth difference inversion, and calculate the morphological values of the iron ore body according to the magnetization intensity distribution range.

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

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

[0070] Obtain the magnetic anomaly area caused by the concealed magnetite ore body from the magnetization pole data.

[0071] Combined 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 ore body based on the magnetic anomaly area includes:

[0072] Make a magnetic anomaly curve according to the magnetic anomaly data, and calculate the buried depth of the first top surface by using the tangent of the characteristic points on the magnetic anomaly curve and the abscissa of the intersection points between the tangents;

[0073] Or, obtain the buried depth of the first top surface according to the known drilling data of the iron ore body.

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

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

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

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

[0078] Combined with this embodiment, there is also a preferred implementation scheme. Specifically, the method for calculating the buried depth of the first top surface by using the tangent of the characteristic points on the magnetic anomaly curve and the abscissa of the intersection points between the tangents includes:

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

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

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

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

[0083] Suppose there are n calculated first top surface burial depths of magnetite ore bodies based on the magnetic anomaly region, and the first top surface burial depth at the i-th point is h i , the vertical projection depth of the i-th point on the top surface of the magnetite ore body delineated by the preset magnetization intensity z value is z i , the vertical projection depth of the i-th point on the top surface of the magnetite ore body delineated by 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 ore body delineated by the preset magnetization intensity z value is overall closer to the depth point calculated by the tangent method than the depth of the top surface of the magnetite ore body delineated by other magnetization intensity values.

[0085] As Figure 3 and Figure 4 shown, M1 - M3 are the top surfaces of three magnetite ore bodies delineated with different parameters, the blue dots represent the calculated top surface burial depth of the magnetite ore body by the tangent method, and the curves represent the cross-sectional curves of the isosurfaces of the magnetite ore body delineated with different magnetic susceptibilities. Figure 4 The vertical cross-sections of the three top surfaces of the magnetite ore bodies with different colors in Figure 3 correspond to the three curves shown in Figure 4 , and the blue dots in Figure 3 are the depth values at different positions calculated by the tangent method on the cross-section. It can be seen from

[0086] that the M2 curve calculated with a certain parameter is closer to the calculated top surface burial depth value by the tangent method compared to the M1 and M3 curves. Then, the magnetization intensity z value used by the M2 curve is used to calculate the morphological values of the iron ore body.

[0087] Example Two:

[0088] This Example Two provides a system for calculating the three-dimensional morphology of an iron ore body. As Figure 5 shown, the system includes:

[0089] A data module for obtaining magnetic anomaly data and magnetic anomaly regions of known buried magnetite ore bodies in the target area;

[0090] A calculation module for calculating the first top surface burial depth of the magnetite ore body based on the magnetic anomaly region and obtaining the three-dimensional magnetization intensity distribution of the magnetic anomaly region;

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

[0092] A statistics module, configured to calculate the minimum value of the difference between the second top surface depth of the iron ore body calculated according to the preset magnetization intensity and the depth of the first top surface;

[0093] An extraction module, configured to extract the range of the magnetization intensity distribution inverted by the magnetization intensity inversion result that is greater than the minimum value of the difference in depth, and calculate the morphological value of the iron ore body according to the range of the magnetization intensity distribution.

[0094] Embodiment 3:

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

[0096] One or more processors;

[0097] A storage device, configured to store one or more programs, which when executed by one or more processors, cause the one or more processors to implement the method for calculating the three-dimensional shape of an iron ore body according to any one of Embodiment 1.

[0098] Figure 6 This is a schematic structural diagram of the device for calculating the three-dimensional shape of an iron ore body provided in Embodiment 3. Figure 6 It shows a block diagram of an exemplary device for calculating the three-dimensional shape of an iron ore body suitable for implementing the embodiments of the present invention. Figure 6 The shown device for calculating the three-dimensional shape of an iron ore body is only an example and should not impose any limitation on the functions and the scope of use of the embodiments of the present invention.

[0099] As Figure 6 shown, the device for calculating the three-dimensional shape of an iron ore body is presented in the form of a general-purpose device. The 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, a memory, and a bus connecting different system components (including the memory and the processing unit).

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

[0101] Devices for calculating the three-dimensional shape of iron ore bodies typically include a variety of computer system-readable media. These media can be any available media accessible by the devices that can be modified by intelligent logging interpretation models, including volatile and non-volatile media, removable and non-removable media.

[0102] The memory can include computer system-readable media in the form of volatile memory, such as random access memory (RAM) 30 and / or cache memory. The device for calculating the three-dimensional shape of iron ore bodies can further include other removable / non-removable, volatile / non-volatile computer system storage media. By way of example only, a storage system can be used for reading and writing on non-removable, non-volatile magnetic media ( Figure 6 not shown, commonly referred to as a "hard disk drive"). Although Figure 6 not shown in, a disk drive for reading and writing on removable non-volatile disks (such as "floppy disks") can be provided, as well as an optical disk drive for reading and writing on removable non-volatile optical disks (such as CD-ROM, DVD-ROM or other optical media). In these cases, each drive can be connected to the bus through one or more data media interfaces. The memory can include at least one program product having a set (such as at least one) of program modules configured to perform the functions of the embodiments of the present invention.

[0103] A program / utility with a set (at least one) of program modules can be stored, for example, in the memory. Such program modules include, but are not limited to, an operating system, one or more application programs, other program modules, and program data. Each or some combination of these examples may include the implementation of a network environment. The program modules typically perform the functions and / or methods in the embodiments described in the present invention.

[0104] The device for calculating the three-dimensional shape of iron ore bodies can also communicate with one or more external devices (such as a keyboard, a pointing device, a display, etc.), and can also communicate with one or more devices that enable a user to interact with the device for calculating the three-dimensional shape of iron ore bodies, and / or communicate with any device that enables the device for calculating the three-dimensional shape of iron ore bodies to communicate with one or more other devices (such as a network card, a modem, etc.). Such communication can be carried out through an input / output (I / O) interface. And, the device for modifying the intelligent logging interpretation model can also communicate with one or more networks (such as a local area network (LAN), a wide area network (WAN) and / or a public network, such as the Internet) through a network adapter. As Figure 6As shown, the network adapter communicates with other modules of the device for calculating the three-dimensional shape of the iron ore body through a bus. It should be understood that although not shown in the figure, other hardware and / or software modules can be used in combination with the device for calculating the three-dimensional shape of the 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, etc.

[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 shape of the iron ore body provided in any embodiment of the present invention. That is: obtaining the magnetic anomaly data and magnetic anomaly area of the known buried magnetite ore body in the target area; calculating the first top surface burial depth of the magnetite ore body based on the magnetic anomaly area or obtaining the three-dimensional distribution of the magnetization intensity of the magnetic anomaly area through drilling data; delineating the range of magnetic bodies with a magnetization intensity less than a preset magnetization intensity according to the first top surface burial depth of the magnetite ore body and the three-dimensional distribution of the magnetization intensity; calculating the second top surface depth of the iron ore body according to the preset magnetization intensity, and statistically obtaining the minimum value of the difference between the second top surface depth and the first top surface burial depth; extracting the range of the magnetization intensity distribution inverted by the magnetization intensity inversion result that is greater than the depth difference minimum value, and calculating the morphological value of the iron ore body according to the range of the magnetization intensity distribution.

[0106] The above are only the preferred embodiments of the present invention, and are not intended to limit the present invention. Any modifications, equivalent replacements, and improvements made within the spirit and principle of the present invention shall be included in the protection scope of the present invention.

Claims

1. A method for calculating the three-dimensional shape of an iron ore body, characterized in that the method Including: Obtaining magnetic anomaly data and magnetic anomaly regions of known concealed magnetite ore bodies in the target area; Calculating the first top surface burial depth of the magnetite ore body based on the magnetic anomaly region or obtaining the three-dimensional magnetization intensity distribution of the magnetic anomaly region through drilling data; Defining the range of magnetic bodies with magnetization intensity less than the preset magnetization intensity according to the first top surface burial depth of the magnetite ore body and the three-dimensional magnetization intensity distribution; Calculating the second top surface depth of the iron ore body according to the preset magnetization intensity, and statistically obtaining the minimum value of the difference between the second top surface depth and the first top surface burial depth; Extracting the magnetization intensity distribution range inverted by the magnetization intensity inversion result greater than the minimum value of the depth difference, and calculating the morphological values of the iron ore body according to the magnetization intensity distribution range.

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

3. The method for calculating the three-dimensional shape of an iron ore body according to claim 2, wherein The method for obtaining the magnetic anomaly region within the range includes: Obtaining the magnetic anomaly region caused by the concealed magnetite ore body from the magnetized pole data.

4. The method for calculating the three-dimensional shape of an iron ore body according to claim 3, characterized in that The method for calculating the first top surface burial depth of the magnetite ore body based on the magnetic anomaly region includes: Making a magnetic anomaly curve according to the magnetic anomaly data, and calculating the first top surface burial depth by using the tangents of characteristic points on the magnetic anomaly curve and the abscissas of the intersection points between the tangents; Or, obtaining the first top surface burial depth according to the known drilling data of the iron ore body.

5. The method for calculating the three-dimensional shape of an iron ore body according to claim 4, characterized in that Five tangents and four intersection points are drawn on the magnetic anomaly curve according to the tangent method.

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

7. The method for calculating the three-dimensional shape of an iron ore body according to claim 6, wherein The method for calculating the first top surface burial depth by using the tangents of characteristic points on the magnetic anomaly curve and the abscissas of the intersection points between the tangents includes: h = K(x1 - x2 + x3 - x4) / 2 = K(b1 + b2) / 2 Where: h is the first top surface burial depth, K is a correction coefficient, x is the abscissa of the intersection point between the tangents, and the abscissas are x1, x2, x3, x4 in the order from right to left, b1 is the difference between x1 and x2, and b2 is the difference between x3 and x4; Or, obtaining the actual burial depth of the iron ore body through drilling data as the first top surface burial depth.

8. The method for calculating the three-dimensional shape of an iron ore body according to claim 7, characterized in that, The method for defining the range of magnetic bodies with magnetization intensity less than the preset magnetization intensity according to the first top surface burial depth of the magnetite ore body and the three-dimensional magnetization intensity distribution includes: Suppose there are n calculated first top surface burial depths of the magnetite ore body based on the magnetic anomaly region, and the first top surface burial depth of the i-th point is h i , the vertical projection depth of the i-th point on the top surface of the magnetite ore body delineated by the preset magnetization intensity z value is z i , the vertical projection depth of the i-th point on the top surface of the magnetite ore body delineated by other magnetization intensity g values is g i , the preset magnetization intensity z value should satisfy:

9. A system for calculating the three-dimensional shape of an iron ore body, characterized in that, The system includes: A data module for obtaining magnetic anomaly data and magnetic anomaly regions of known concealed magnetite ore bodies in the target area; A calculation module for calculating the first top surface burial depth of the magnetite ore body based on the magnetic anomaly region and obtaining the three-dimensional magnetization intensity distribution of the magnetic anomaly region; A defining module for defining the range of magnetic bodies with magnetization intensity less than the preset magnetization intensity according to the first top surface burial depth of the magnetite ore body and the three-dimensional magnetization intensity distribution; A statistical module for calculating the second top surface depth of the iron ore body according to the preset magnetization intensity and statistically obtaining the minimum value of the difference between the second top surface depth and the first top surface burial depth; An extraction module, configured to extract the magnetization intensity distribution range in the magnetization intensity inversion result that is greater than the minimum value of the depth difference, and calculate the morphological value of the iron ore body according to the magnetization intensity distribution range.

10. An apparatus for calculating the three-dimensional shape of an iron ore body, characterized in that the apparatus Comprising: One or more processors; A storage device, configured to store one or more programs, which, when executed by the one or more processors, cause the one or more processors to implement the method for calculating the three-dimensional morphology of the iron ore body as described in any one of claims 1-8.

Images